Pamburayi
Mpofu
a,
Houyem
Hafdi
a,
Pentti
Niiranen
a,
Jonas
Lauridsen
b,
Oscar
Alm
b,
Tommy
Larsson
b and
Henrik
Pedersen
*a
aDepartment of Physics, Chemistry and Biology, Linköping University, SE-581 83, Linköping, Sweden. E-mail: henrik.pedersen@liu.se
bSeco Tools AB, SE-737 82, Fagersta, Sweden
First published on 17th July 2024
Aluminum nitride (AlN) is a semiconductor with a very wide band gap and a potential dielectric material. Deposition of thin AlN films is routinely done by several techniques, including atomic layer deposition (ALD). In this study, we deposited AlN using ALD with trimethylaluminum (TMA) as the Al precursor and ammonia (NH3) with and without plasma activation as the N precursor in the temperature range from 100 to 400 °C while monitoring the surface reactions using mass spectrometry. Our results, combined with recent quantum chemical modelling, suggest that the surface chemistry of the deposition process is chemisorption of TMA followed by reductive elimination of the methyl groups to render mono methyl aluminum species. The NH3 chemisorption is done by ligand exchange to form CH4 and an –NH2 terminated surface.
ALD is based on sequential, self-limiting surface chemical reactions, enabled by carefully tuned deposition chemistry and an alternating precursor supply onto the substrate. This leads to the hallmark of ALD: conformal deposition on topologically complex surfaces. The alternating surface chemistry also facilitates studies of the surface chemistry and reaction mechanisms as the surface is only exposed to one type of precursor molecule at a time. Puurunen et al.8–10 have explored the surface chemistry of AlN using trimethylaluminum (TMA) and ammonia (NH3) by infrared (IR) spectroscopy and nuclear magnetic resonance (NMR) spectroscopy on porous silica and alumina. They found that the amount of hydrogen in amino (–NHx) groups increased with the number of reaction cycles on both alumina (Al2O3) and silica (SiO2) substrates, suggesting that more NHx groups were formed on the Al2O3 and SiO2 surfaces as the reactions progressed. The interaction of TMA with dehydroxylated SiO2 involved ligand exchange with surface hydroxyl (–OH) groups and dissociation in siloxane (–Si–O–Si–) bridges, leading to surface saturation with methyl (CH3) groups, followed by the replacement of these CH3 groups by NHx groups upon reaction with NH3, indicating a stepwise transformation of the surface chemistry. Additionally, upon reacting with silica, TMA generated –Si–CH3 groups on the surface, and the subsequent interaction with NH3, particularly at higher temperatures, led to the removal of over 90% of the initial CH3 groups, which aligned with the proposed reaction pathway for this process.
A recent modelling study of the surface reconstruction after the NH3 pulse in ALD of AlN, GaN and InN,11 suggests that, at typical ALD temperatures and pressures, the AlN surface is predominantly covered by –NH2 species at thermal equilibrium. The surface chemistry of TMA and NH3 in ALD of AlN was also recently mapped out using density functional theory (DFT),12 suggesting that the decomposition pathways of TMA are crucial for determining film growth rates. Three distinct decomposition pathways for TMA molecules on the surface were observed: (i) TMA undergoes reductive elimination to produce monomethyl aluminum on the surface and ethane gas; this pathway does not involve the removal of surface hydrogen atoms, (ii) methyl protonation followed by reductive elimination where TMA interacts with two surface hydrogen atoms, leading to the formation of monomethyl aluminum on the surface, along with methane and ethane gas, surface hydrogen atoms are transferred to the methyl groups during this decomposition pathway, and (iii) methyl protonation twice where TMA reacts with a surface hydrogen atom to generate an aluminum adatom, methane and ethane gas, this process involves the removal of hydrogen atoms from the surface during decomposition.
These pathways were found to be pivotal in determining surface coverage, hydrogen removal, and the overall AlN film-growth efficiency. The DFT study found that TMA would most likely adsorb as a Lewis adduct on the –NH2 covered surface and decompose via reductive elimination or methyl protonation pathways, with the proximity to previously adsorbed TMA molecules influencing the decomposition pathway. The ratio of the different decomposition pathways was suggested to affect the growth per cycle (GPC) by influencing the saturation coverage of the surface. The relative concentrations of TMA molecules decomposing via the three pathways determine the availability of surface sites for TMA adsorption. If the dominant pathway is methyl protonation, which removes surface hydrogen and allows for additional TMA adsorption sites, the GPC may increase due to higher surface saturation coverage. On the other hand, if the dominant pathway is reductive elimination, which does not involve hydrogen removal, the GPC may be lower as the surface may become saturated with hydrogen atoms, limiting additional TMA adsorption sites.
Herein, we use mass spectrometry to investigate the surface reactions in both thermal and plasma ALD of AlN, using TMA as the Al precursor and NH3, with and without plasma activation, as the N co-reactant. The use of in situ mass spectrometry for real-time analysis of gas-phase species during the reactions of this precursor combination has not been reported hence it is not fully understood. Our results add to the previous DFT study to form a surface chemical mechanism of ALD of AlN using TMA.
Film thicknesses were measured by spectroscopic ellipsometry (SE) that uses CompleteEASE software for data analysis and reporting of results. Ellipsometric spectra were recorded in the wavelength range of 210–1690 nm at four angles of incidence (45°, 55°, 65°and 75°) using a variable angle a J.A. Wollam RC2 Mueller matrix spectroscopic ellipsometer. The thicknesses of the films were then obtained by modeling the measured data using Cauchy dispersion function. Lower film thicknesses (e.g., below 100 nm) were measured by X-ray reflectivity (XRR) with the software PANalytical X’Pert reflectivity and a two-layer model to fit the data and calculate the film thickness. Crystallinity was measured in both θ–2θ and grazing incidence X-ray diffraction (GIXRD) geometries using a PANalytical EMPYREAN MRD diffractometer with a Cu X-ray tube, wavelength of 1.54 Å, operating at 45 kV and 40 mA.
X-ray photoelectron spectroscopy (XPS) was used to obtain elemental compositions and chemical bonding. The results were obtained using Kratos AXIS Ultra DLD, equipped with an Ar+ sputtering source. During acquisition, a monochromatic Al Kα X-ray radiation, wavelength 8.35 Å, with 150 W power (anode current = 10 mA and anode voltage = 15 kV) was used. The recorded photoelectron spectra were analyzed using the CasaXPS software package. By quantitative analysis in Casa XPS, the signals originating from the thin films could be de-convoluted and the chemical composition could be quantified. Gaussian–Lorentz functions and Shirley background were used to fit the experimental XPS data.
Optical emission spectroscopy (OES) utilized the Mechelle 900 Spectrometer with a cooled digital 12-bit SensiCam CCD camera system, covering a wavelength range from ultraviolet (UV) to near infrared (NIR), i.e., 200–1100 nm. Positioned perpendicular to the plasma source, the spectrometer captured light emitted from the NH3 plasma, directing it to a monochromator through an optical fiber cable. Measurements were conducted immediately after the ICP source in the gas flow direction, before reaching the substrate which is located approximately 65 cm further down which is the only available view port where we can observe the plasma in our reactor set up. Data analysis considered the spectrometer's consistent spectral resolution of 900, enabling the resolution of wavelengths separated by 1/900 of the wavelength.
Gaseous species in the exhaust from the ALD chamber were measured using a Hiden Analytical HPR-30 residual gas analysis (RGA) vacuum process sampling system that includes a HAL 201 RC mass spectrometer (MS) with a Faraday cup detector. The RGA system was connected approximately 80 cm from the substrate stage, with gas lines maintained at room temperature since no residue condensation along the lines was noted. The MS analysis involved collecting both scan spectra and multiple ion detection (MID) of specific ions with electron impact ionization set at 70 eV to ensure efficient reactant molecule fragmentation, enabling the detection and analysis of gas-phase reaction products. The system was differentially pumped and the operating pressure during typical measurements is kept in the order of 10−6 mbar to enable mean free paths that allow collisionless particle transport of the ions from the orifice to the detector. The obtained data was analysed using MASsoft 10 (Hidden Analytical) software.
The plasma process showed a temperature window where the GPC is not affected significantly by temperature between 250 and 325 °C (Fig. 2a). We therefore carried out saturation growth experiments for this process at 300 °C which is within this temperature window. The thermal process showed no clear temperature window where the GPC is not affected by the temperature (Fig. 2b) and very little films were deposited at temperatures below 300 °C. From this, we used 350 °C for studies of the thermal process. Nearly stoichiometric and polycrystalline films were deposited at 300 °C and 350 °C for the plasma and thermal processes respectively, see ESI,† Table S1 and Fig. S1, S2.
Fig. 2 Growth dependence on process temperature for the plasma (a) and thermal (b) processes using 0.2 s TMA and 12 s NH3 pulse times. Lines are a guide for the eyes. |
The optically active plasma species were studied by OES (Fig. 3), collected when the reactor walls were 300 °C, with a plasma power of 2800 W, and NH3 and Ar flow rates of 100 sccm each. Although the information obtained from OES is limited, relying on excitation processes and radiative transitions, it can give qualitative insight into the species produced in the plasma. The plasma emission lines from OES included NH at 336.1 nm, N2 at 357.1 nm, and N2+ at 391.4 nm and 427.6 nm.13,14 Additionally, the emission peaks at around 336–337 nm indicated the presence of excited NH and N2 radicals,14 corresponding to specific transitions, supporting the existence of NHx and N2 species.
A mass spectrometric survey scan was carried out to see the residual reaction species and byproducts leaving the reactor (Fig. 4). The main residual gases from the ALD reactor are H2, N, CH4, NH3, Ar, N2 and Al(CH3), while CH3 and C2H6 are detected as minor species. No other species were detected beyond m/z = 60. From these species, CH4 is neither a molecule delivered into the reactor nor a direct fragment of the dosed molecules (TMA and NH3) but is expected to be a product of surface chemical reactions. Hence, we studied and monitored its evolution with temperature over time as illustrated in Fig. 5. It can be noted that the signal from m/z = 16 increases rapidly with the NH3 pulse for both the thermal and the plasma process. As the temperature increases, the m/z = 16 signal is also detected with the TMA pulse. We note that the m/z = 16 signal could also emanate from NH2+, formed by decomposition of NH3 in the QMS. This is likely the explanation for the double peaks observed at the NH3 pulses, the first peak is the NH2+ from the NH3 pulse, the second is the CH4+ from the surface reactions.
Fig. 4 A survey scan showing average intensities of all detected masses during 5 ALD cycles of the TMA-NH3 process. No other species were detected beyond m/z 60. |
From Fig. 5, activation barriers were calculated using the Arrhenius equation, which relates the rate constant of a reaction to temperature and activation energy. The Arrhenius equation is given by:
k = A·e−(Ea/RT) |
From the slope of the tangent of the decay of the m/z = 16 signal, in Fig. 5, we extracted an approximation for the rate constants for the methane evolution from the surface reactions. By doing this for all three studied temperatures, we could approximate the activation energies (Ea) for the methane evolution, using Arrhenius plots shown in ESI,† Fig. S3. The extracted activation energies for methane evolution were 19–20 kJ mol−1 during the TMA pulse, 14 kJ mol−1 during the plasma activated NH3 pulse, and 17 kJ mol−1 during the thermally activated NH3 pulse. The shape of the Arrhenius plots suggests that the NH3 pulses are in a mass transport limited regime at T ≥ 300 °C, while the TMA pulses seem to be in a surface kinetics limited regime at all temperatures studied. However, considering the small differences in the mentioned activation energies, a comparison between the thermal and plasma processes could possibly show that the surface reactions may not be that different.
We thus suggest that the TMA molecules adsorb to the surface through their interaction with -NH2 surface groups, and this initial half-cycle can be defined by reaction (1):
‖– NH2 + Al(CH3)3 → ‖– NH2−xAl(CH3)3−x + xCH4(g) | (1) |
The saturation curve for the plasma NH3 pulse (Fig. 1a) shows differences to the thermal NH3 pulse (Fig. 1b), suggesting that the thermal reaction is not as strictly self-limited as the plasma reaction. However, such differences are not shown by the residual gas analysis, as both processes show the same pattern of evolution of methane with temperature (Fig. 5). The observation that no reaction occurs below 300 °C in the thermal process agrees with previous studies that showed that the NH3 surface reactions with TMA in thermal ALD only occur with reasonable rates at temperatures close to where decomposition of TMA takes place.16 The presence of hydrogen radicals in the NH3 plasma (Fig. 3), suggests that hydrogenation processes of the surface species could play a role in plasma ALD of AlN, which would limit TMA adsorption. This could explain why the plasma process has slightly lower GPCs than the thermal process. The difference in the saturation curves could also be a sign of difference in reaction activity. This is further supported by a higher nucleation delay in the thermal process compared to that in the plasma process (Fig. S4, ESI†).
As the NH3 plasma ALD process is expected to be further from thermal equilibrium, the surface termination may not be the same for the plasma process. The detection of NH radicals in the OES spectra (Fig. 3) indicates the presence of more reactive nitrogen species in the plasma.
If we assume that both the thermal and plasma process renders –NH2 surface species in the NH3 pulse, in accordance with recent modelling,11 the elimination of the –CH3 ligands is facilitated by protonation reactions, resulting in the generation of CH4 as follows:
‖– Al(CH3)x + x(NH3) (g) → ‖– Al–(NH2)x + xCH4(g) | (2) |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4tc01867b |
This journal is © The Royal Society of Chemistry 2024 |