Francesca
Lo Presti
,
Anna L.
Pellegrino
and
Graziella
Malandrino
*
Dipartimento di Scienze Chimiche, Università degli Studi di Catania, INSTM UdR Catania, Viale Andrea Doria 6, I-95125 Catania, Italy. E-mail: graziella.malandrino@unict.it
First published on 22nd April 2022
Atmospheric pressure (AP) vapor phase processes such as spatial atomic layer deposition (S-ALD) and AP-metalorganic chemical vapor deposition (AP-MOCVD) are becoming increasingly appealing for their use in a variety of academic and industrial applications. Evaluation of precursor vapour pressures is crucial for their application in AP processes and to this aim the Langmuir equation has been applied as a simple and straightforward method for estimating the vapor pressure and vaporization enthalpy of various metalorganic precursors. Using benzoic acid as a calibration reference, the vapour pressure–temperature curves for several alkaline-earth β-diketonate fluorinated compounds, with molecular formula “M(hfa)2·L” (with M = Mg, Ca, Sr, Ba; Hhfa = 1,1,1,5,5,5-hexafluoroacetylacetone and L = diglyme, triglyme, and tetraglyme) are derived from their termogravimetric curves. Thus, the enthalpy of vaporization of all complexes has been estimated using the Clausius–Clapeyron equation. As a proof of concept, preliminary results on the use of [Mg(hfa)2·2H2O]·2diglyme and [Ca(hfa)2·diglyme·H2O] or [Ca(hfa)2·triglyme] as precursors for AP-MOCVD deposition of MgF2 and CaF2 in the form of thin films are presented. This approach may be used to easily determine vapor pressures of complexes and thus evaluate “a priori” the suitability of a compound as precursor for AP-MOCVD and/or spatial ALD.
Great efforts have been devoted to overcome those issues, such as low volatility and poor thermal stability due to the lack of the saturation of the metal centre coordination sphere resulting in the formation of non-volatile oligomeric structures. These drawbacks have stimulated the synthesis of new metalorganic β-diketonate compounds, named “second-generation” precursors, in which the addition of neutral ligands such as polyethers prevents the oligomerization and water-coordination, thus improving thermal stability and volatility.10,11
In addition, along with the development of new and promising precursors, it is also critical to investigate their physico-chemical properties, particularly their thermal behaviour and vapor phase stability, which are essential characteristics for their application in metal–organic chemical vapor deposition (MOCVD) procedures. The above-mentioned “second-generation” precursors have been widely applied in MOCVD approaches under reduced pressure conditions,12–15 and their thermal properties make them also suitable for ALD processes.
Recently, another form of ALD processes is emerging, i.e. the spatial ALD (SALD),16,17 which has the great advantage of a much higher growth rate per cycle, but this variant is regularly carried out under atmospheric pressure.
These conditions make the thermal requirements of precursors even more stringent since they should possess a high volatility at atmospheric pressure. Under the light of these premises, optimizing a simple approach to evaluate the vapor pressures of complexes is fundamental to have a reference point on the suitability of a given precursor for application in AP-vapor phase processes. This information would avoid a full trial-and-error approach in the optimization of film deposition through SALD.
Herein, a simple and functional approach has been implemented to evaluate the vapor pressures of the fluorinated β-diketonate alkaline-earth adducts of the type “M(hfa)2·glyme” [with Hhfa = 1,1,1,5,5,5-hexafluoroacetylacetone, M = Mg, Ca, Sr and Ba; glyme = diglyme bis(2-methoxyethyl)ether, triglyme (2,5,8,11-tetraoxadodecane) and tetraglyme (2,5,8,11,14-pentaoxapentadecane)] for application in atmospheric pressure vapor phase processes. Thermogravimetric (TGA) analyses have been used to firstly confirm thermal stability and volatility of these complexes. Starting from the experimental data of the thermogravimetric analyses, we deeply discuss the possibility of using the Langmuir equation for the direct estimation of vapor pressure and enthalpy of vaporization for the principal alkaline-earth “second-generation” β-diketonate fluorinated precursors.
Furthermore, as a proof of concept, the [Mg(hfa)2·2H2O]·2diglyme and the calcium complexes, [Ca(hfa)2·diglyme·H2O] and [Ca(hfa)2·triglyme], have been tested for the fabrication of nanostructured MgF2 and CaF2, respectively, thin films on silicon substrates through atmospheric pressure processes in a customized horizontal hot-wall reactor.
Alkaline-earth metal | Ligands | ||
---|---|---|---|
Diglyme | Triglyme | Tetraglyme | |
[Mg(hfa)2·2H2O]·2diglyme19 | |||
[Ca(hfa)2·diglyme·H2O]24 | [Ca(hfa)2·triglyme]21,27 | [Ca(hfa)2·tetraglyme]18,25 | |
[Sr(hfa)2·diglyme·H2O]24 | [Sr(hfa)2·triglyme·H2O]228 | [Sr(hfa)2·tetraglyme]18 | |
[Ba(hfa)2·(diglyme)2]24 | [Ba(hfa)2·triglyme·H2O]29 | [Ba(hfa)2·tetraglyme]18,26 |
All the alkaline-earth metals forming the above-mentioned precursors, have a +2 oxidation state, but their coordination number tends to increase (Mg < Ca < Sr < Ba) with the increasing of the ionic radius. Nevertheless, there is essentially no preferential coordination geometry, at difference with transition metals, with the exception of a preferred spherical distribution of the donor atoms.30
The nature of the adducts has been assessed through comparison of their melting points and FT-IR spectra with literature data. All the FT-IR spectra show a peak at 1660 cm−1 associated with the CO stretching, whose shift, compared to the free ligand, confirms coordination of the hfa, hence formation of the complexes. Around 3500 cm−1 a broad band is observed in those adducts which have coordinated water molecules, while in the range 1100–750 cm−1 the peaks associated with the glyme coordination are observed. As an exemplificative case, the FT-IR spectra of the Ba adducts are reported in Fig. 1, while the other spectra are reported in the ESI in Fig. S1–S3.† In the case of the Ba adduct spectra (Fig. 1), the saturated signal at 2900 cm−1 is actually due to the C–H stretching of the nujol overlapping the C–H stretching signal of the glyme.
(1) |
P = υ·k | (2) |
(3) |
For a defined calibration standard of mass M, by plotting the logarithm of P (Pa) calculated from eqn (3)vs. υ derived directly from the TGA/DTG experiment, the value of k can be estimated.34
It is worth noting that the Antoine equation is closely linked to a specific temperature range, thus, the results will be more precise the smaller the temperature range employed to calculate them. Therefore, for each compound the vapor pressure can be calculated using eqn (2) if k is known.
The Clausius–Clapeyron equation35 (eqn (4)) can additionally be used to determine the enthalpy of vaporization of the substances under investigation:
(4) |
In the eqn (4), P is the vapor pressure (Pa), A is a constant, ΔH is the enthalpy of vaporization (kJ mol−1) and R is the gas constant (J K−1 mol−1). Using the vapor pressure from eqn (3), by plotting the ln of P (Pa) versus 1/T (K−1) the ΔH may be extrapolated from the slope of the straight line. It is worth noting that the Clausius–Clapeyron equation only applies to substances that evaporate and sublimate in a zero-order process.
Using eqn (2), the k value was calculated from three distinct experiments and found to be 839325 Pa kg1/2 m−2 s−1 K−1/2 on average (the relative plot is reported in Fig. S4†). This value is similar to previous literature estimates; the small variations are due to the fact that k is highly dependent on the instrumentation and operating conditions.39,41,42
Regardless, a back-test has been performed to verify the accuracy of the result. Therefore, the vapor pressure of another chemical, salicylic acid, whose vapor pressure values are well-known in the literature36,43 have been estimated. More specifically, salicylic acid vapor pressure has been assessed in the temperature range 416–477 K, using eqn (3) and the Antoine constant from ref. 3, and has been compared to the values obtained using eqn (2) with the k value derived from the benzoic acid analysis. Fig. S5† depicts the results, which indicate a good match with the reference data. The salicylic acid enthalpy of vaporization calculated using vapor pressure data has been compared to other references. The enthalpy of vaporization obtained is 76.03 ± 0.55 kJ, and compares very well with the value of 79 kJ obtained in ref. 36. TGA/DTG and ln of P vs. 1/T of salicylic acid are shown in Fig. S6.†
The TGA plot displays two weight loss regions: the first, in the temperature range 308–383 K, is associated with dehydration, and the second, in the temperature range from 402 to 461 K, is related to the quantitative evaporation of the anhydrous complex, [Mg(hfa)2·2diglyme], leaving a very low residue.
Fig. 4a shows the graph of the vapour pressure vs. temperature, while the graph in Fig. 4b, reporting the lnP in function of 1/T, allowed the estimation of the enthalpy of vaporization calculated using the Clausius–Claypeyron method. The estimated ΔH for a temperature range of 402–461 K is 55.1 ± 0.3 kJ mol−1, whereas the extrapolated vapor pressure at 423 K is 365.1 Pa.
Fig. 5 TGA curves (a) and DTG plot (b) of “Ca(hfa)2·L” compounds (with L = diglyme, triglyme, tetraglyme). Scheme of the [Ca(hfa)2·triglyme] is reportd as inset in (b). |
Regarding the [Ca(hfa)2·tetraglyme] structure, the Ca2+ is octacoordinated by the four oxygens from the hfa units and four of the five oxygen atoms from the tetraglyme ligand.25
The TGA/DTG plots for the Ca adducts with different glyme ligands are reported in Fig. 5a and b.
The presence of a minimum from 333 to 358 K for the [Ca(hfa)2·diglyme·H2O] (orange curve) is observed on the TGA plot, corresponding to the elimination of the water molecule, while the main weight loss at higher temperatures refers to the vaporization of the anhydrous complex [Ca(hfa)2·diglyme], leaving a residue of less than 1%. The thermal behaviour of [Ca(hfa)2·triglyme] (brown line) and [Ca(hfa)2·tetraglyme] (black line) is remarkably comparable in the temperature ranges of 423–529 K and 433–523 K, respectively. These curves show only one-step major loss and no residues, confirming quantitative vaporization processes for both complexes. Fig. 6 shows the vapor pressure trends as a function of temperature.
Fig. 6 Vapor pressure/temperature (a) and lnP vs. 1/T (b) plots of Ca(hfa)2·L compounds (with L = diglyme, triglyme, tetraglyme). |
According to the lnP vs. 1/T plot, the ΔH values for [Ca(hfa)2·diglyme·H2O], [Ca(hfa)2·triglyme], and [Ca(hfa)2·tetraglyme] are 64.4 ± 0.2 kJ mol−1, 68.6 ± 0.2 kJ mol−1 and 61.9 ± 0.2 kJ mol−1, respectively. The results show that the “Ca(hfa)2glyme” adducts have similar vaporization enthalpies and thus comparable volatilities under atmospheric conditions.
Fig. 7 TGA curves (a) and DTG plots (b) of “Sr(hfa)2·L” adducts (with L = triglyme, tetraglyme). Scheme structure of the [Sr(hfa)2·tetraglyme] is reported as inset in (b). |
The loss of coordinated water molecules from the [Sr(hfa)2·triglyme·H2O]2 complex (light green line) results in a slight weight loss at 345–373 K. Then, in the range of 466–528 K, the dehydrated species [Sr(hfa)2·triglyme] evaporates with one-step weight loss, leaving no residue. The [Sr(hfa)2·tetraglyme] adduct, on the other hand, shows a one-step weight loss of 98% in the temperature range of 454–545 K. Hence, the enthalpy of vaporization estimated from the vapor pressure vs. temperature plots (Fig. S7†), is equal to 60.4 ± 0.4 kJ mol−1 and 74.3 ± 0.1 kJ mol−1 for the [Sr(hfa)2·triglyme·H2O]2 and the [Sr(hfa)2·tetraglyme], respectively. In comparison to the compound with a longer chain ligand, it appears that the Sr adduct with triglyme has a lower enthalpy of vaporization, and hence a better thermal behaviour, with a trend similar to calcium ones. As a consequence, the predicted vapor pressure for both compounds at 423 K is 73.4 Pa for the [Sr(hfa)2·triglyme·H2O]2 and 34.8 Pa for the [Sr(hfa)2·tetraglyme] species.
Consequently, only the [Ba(hfa)2·triglyme·H2O] and [Ba(hfa)2·tetraglyme] are the Ba adducts under investigation in this work. The Ba(II) ion in the [Ba(hfa)2·triglyme·H2O] structure is nine coordinated by eight oxygens from the hfa and triglyme and an oxygen from the water molecule completing the coordination sphere.29
When tetraglyme is used, the 9-coordination of Ba(II) is reached with the four oxygens of the hfa and all the five oxygens from the tetraglyme ligand (inset in Fig. 8b), yielding the anhydrous volatile [Ba(hfa)2·tetraglyme].18,26 The thermal behaviour of these two adducts is investigated by their TGA analysis (Fig. 8a) and its derivative curve (Fig. 8b).
Fig. 8 TGA curves (a) and DTG plots (b) of Ba(hfa)2·L adducts (with L = triglyme, tetraglyme). Scheme structure of the [Ba(hfa)2·tetraglyme] is reportedas inset in (b). |
For the [Ba(hfa)2·triglyme·H2O] adduct (light blue line), the initial weight loss in the range between 358–373 K is due to the dehydration process, leaving the [Ba(hfa)2·triglyme] anhydrous species. The major vaporization process starts at a temperature of 438 K and ends at 542 K with a quite high residue of about 20%. It is worth to note that even though a certain residue is left behind, the ΔH and vapor pressure may still be evaluated considering that the product of dehydratation of the [Ba(hfa)2·triglyme·H2O] has been reported to evaporate congruently.27 Instead, the [Ba(hfa)2·tetraglyme] adduct (dark blue line) has a one-step 88% weight loss in the temperature range between 432–553 K. Fig. S8† depicts the vapor pressure/temperature plots.
The enthalpy of vaporization calculated from the lnP vs. 1/T plot for the [Ba(hfa)2·triglyme·H2O] and [Ba(hfa)2·tetraglyme] is 53.6 ± 0.3 kJ mol−1 and 82.1 ± 0.1 kJ mol−1, respectively. In the case of the Ba-tetraglyme adduct, the longer glyme yields a slightly less volatile precursor.
In literature there are few reports on the synthesis of MgF2 and CaF2 thin films through reduced pressure chemical vapor deposition from the “second-generation” [M(hfa)2·polyether] adducts.9,10,19,24,42,44–46 Recently, a solution process through sol–gel/spin coating approach has been tested as well for the deposition of CaF2.9,45 To the best of our knowledge, the present study represents the first report on the synthesis of MgF2 and CaF2 thin film through atmospheric pressure process starting from the [Mg(hfa)2·2H2O]·2diglyme and the [Ca(hfa)2·diglyme·H2O] or [Ca(hfa)2·triglyme] complexes, respectively.
Structural characterization, performed through XRD analysis, is reported in Fig. 9 for the film deposited at 450 °C on Si (100) substrate. The peaks at 2θ = 35.28°, 40.56°, 43.96°, 53.66°, and 56.14° of the pattern in Fig. 9a correspond to the 101, 111, 210, 211 and 220 reflections of the MgF2 phase (ICDD No. 070-2268). The CaF2 formation (Fig. 9b) is assessed by the presence of peaks at 2θ = 28.28°, 47.04° and 55.78° associated with the 111, 220 and 311 reflections, respectively (ICDD No. 35-0816). It is worth mentioning that a pure polycrystalline CaF2 phase is produced regardless of the applied precursor, either [Ca(hfa)2·diglyme·H2O] or [Ca(hfa)2·triglyme]. The only difference between the two patterns is the relative intensities of the peaks which do not exactly match the intensities reported in the database. Surface morphology investigations have been carried out through field–emission Scanning Electron Microscopy (FE-SEM) analysis in both plan view and cross-sectional configurations (Fig. 10). The FE-SEM images display the formation of uniform and compact layers over the whole surface with very small grains barely visible for both the MgF2 film and CaF2 one grown from the [Ca(hfa)2·diglyme·H2O] precursors. The homogeneity of the fluoride films is also confirmed by cross-sectional images, in which the estimated thickness is in the order of 750–800 nm in all cases. The surface of the film grown from the [Ca(hfa)2·triglyme] adduct, shown in Fig. 10c, is more structured and consists of smaller grains, but the film thickness is equivalent to that deposited with the [Ca(hfa)2·diglyme·H2O] source. Interestingly, the morphology presently observed is quite similar to the one recently found for lanthanide doped CaF2 thin films obtained at 450 °C on silicon using a reduced pressure MOCVD approach.45 This finding points to a direct correlation of the CaF2 film morphology with the kinetic growth process on silicon substrate for a given temperature range.
Finally, compositional characterization, through energy dispersive X-ray analysis (EDX), has been carried out on the MgF2 and CaF2 films in order to assess the compositional purity of the layers. In Fig. S9,† the EDX spectrum of the CaF2 film displays the presence of Ca Kα peak at 3.69 keV and F Kα peak at 0.68 keV. The peaks of Si Kα at 1.74 keV and O at 0.53 keV arise from the silicon substrate. Notably, the absence of the peak at 0.28 keV, due to the C Kα, excludes any carbon contamination. Analogously, the EDX spectrum of the MgF2 layer (Fig. S10†) shows the presence of the Mg Kα peak at 1.25 keV and F Kα peak. Also in this case, no C contamination is observed.
This requirement is crucial for the development of this method and may represent a limitation for the growth affirmation of the spatial ALD in specific applications to alkaline-earth containing films. The glyme adducts of the alkaline-earth metal hexafluoroacetylacetonate complexes are the most volatile, and thus appealing, precursors for chemical vapor phase processes under reduced pressure. Even though the characterizations of their thermal properties have been reported for some of these compounds, up to date only few studies have been reported on the evaluation of their vapor pressures, which are essential for applications in atmospheric pressure processes.
In particular, Tsymbarenko et al. devoted their attention to study the structures, thermal stability and volatility of Ca, Sr, and Ba hexafluoroacetylacetonates mixed ligand complexes with various polyethers.27 For these metal–organic adducts, the authors have investigated their volatility and thermal stability in the gas phase through mass spectrometry (MS) and density functional theory (DFT) calculations.27
Starting from present thermogravimetric analyses, and as confirmed by the studies of Tsymbarenko et al.27 and Luo et al.,48 all the investigated compounds evaporate congruently. Based on these data, we discuss the possibility of using the Langmuir equation for the direct estimation of the vapor pressure and the enthalpy of vaporization for the principal alkaline-earth “second-generation” fluorinated β-diketonate precursors. This approach has been used for the first time by the group of Kunte et al.36 for the estimation of the vapor pressure of two distinct titanium metalorganic precursors using benzoic acid as a standard.
The presently reported Langmuir procedure has allowed to easily estimate the vapor pressure and vaporization enthalpy of alkaline-earth metalorganic compounds. The slope of lnP versus 1/T, with a regression coefficient close to unity, has been used to estimate the enthalpy of vaporization. The vapor pressure at 423 K has then been extrapolated. The trend as a function of the type of ligands coordinating one single metal centre for each set of alkaline-earth precursors is reported in Table 2 and schematized in Fig. 11.
Compound | Temperature range (K) | ΔH (kJ mol−1) | Estimated V. pressure (Pa) at 423 K | Effective ionic radius for 6-fold coordination (Å) | Observed coordination number of the M2+ ion |
---|---|---|---|---|---|
[Mg(hfa)2·2H2O]·2diglyme | 402–461 | 55.1 ± 0.3 | 365.1 | 0.57 | 6 |
[Ca(hfa)2·diglyme·H2O] | 423–503 | 64.4 ± 0.2 | 89.2 | 1.00 | 8 |
[Ca(hfa)2·triglyme] | 423–529 | 68.6 ± 0.2 | 92.4 | ||
[Ca(hfa)2·tetraglyme] | 433–523 | 61.9 ± 0.2 | 78.9 | ||
[Sr(hfa)2·triglyme·H2O]2 | 466–528 | 60.4 ± 0.4 | 73.4 | 1.18 | 9 |
[Sr(hfa)2·tetraglyme] | 454–545 | 74.3 ± 0.1 | 34.8 | ||
[Ba(hfa)2·triglyme·H2O] | 438–542 | 53.6 ± 0.3 | 61.9 | 1.35 | 9 |
[Ba(hfa)2·tetraglyme] | 432–553 | 82.1 ± 0.1 | 23.7 |
It is interesting to comment on the trend observed for the vapor pressure and vaporization enthalpy of the different complexes. (a) A decrease in the vapor pressure of the different complexes has been observed on going from Mg2+ to Ba2+ ion. This trend parallels the larger ion dimension and consequently the higher coordination number on going from Mg2+ (optimal coordination is 6) to Ba2+, which needs at least a 9-coordination. To compare the various alkaline-earth ion dimensions, the Shannon 6-coordinated ionic radii,49 and related effective coordination in the most stable complexes, are reported in Table 2. (b) Related to the vapor pressure values, it is the variation of the vaporization enthalpy observed on going from the lightest to the heaviest among the considered AE metal complexes. These values are slightly smaller than those previously reported by Tsymbarenko et al.:27e.g. 61.9 ± 0.2 vs. 131.7 ± 1.8 (kJ mol−1), 74.3 ± 0.1 vs. 138.5 ± 8.9 (kJ mol−1), 82.1 ± 0.1 vs. 140.9 ± 4.5 (kJ mol−1), for the [Ca(hfa)2·tetraglyme] [Sr(hfa)2·tetraglyme] and [Ba(hfa)2·tetraglyme], respectively. Various factors may be responsible for these findings: (i) significantly different temperature ranges are used in the present and previous evaluation of ΔH values; (ii) completely different approaches have been used to determine these quantities. Significant differences in the ΔHs have also been observed for other adducts, for example for the Ba(β-diket)2·18-crown-6 compounds due to the different Knudsen and static method applied.50 (c) Polyether plays a key role on the vapor pressure of the complexes since the vapor pressure value decreases on increasing the polyether length. Thus, for a given AE ion, the tetraglyme adduct is usually the least volatile with respect to the triglyme or diglyme ones. In this context, a comment deserves the trend observed for the “Ca(hfa)2glyme” vaporization enthalpies and the relationship between vaporization enthalpies and vapor pressures. This trend may be explained considering that in the Clausius–Clapeyron equation (eqn (4)), in addition to the slope, which determines the ΔH, also the intercept should be considered. In addition, vapor pressures of compounds, being derived from eqn (1), are strictly dependent not only on their vaporization rate dm/dt, but also on their molecular mass.
Finally, the potentiality of the most volatile compounds in ALD and SALD processes may be envisaged considering the promising results obtained for the fluorinated β-diketonate Mg and Ca alkaline-earth adducts as single-source precursors in AP-MOCVD processes of binary fluoride thin films. In addition, this type of precursors would represent a great advantage for the ALD of AE metal fluorides since most of the actual reports use an independent source for fluorine, such as TiF4 and TaF5.51
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2dt00479h |
This journal is © The Royal Society of Chemistry 2022 |