Paul-Alexis
Pavard
a,
Eva
Pugliese
ab,
Damien
Coutancier
a,
Valérie
Albin
c,
Nicolas
Casaretto
b,
Sophie
Bourcier
b,
Virginie
Lair
c,
Armelle
Ringuede
c,
Audrey
Auffrant
b and
Nathanaelle
Schneider
*a
aInstitut Photovoltaïque d’Île-de-France (IPVF), UMR 9006, CNRS, Ecole Polytechnique - IP Paris, Chimie Paristech - PSL, 18 Boulevard Thomas Gobert, Palaiseau, 91120, France. E-mail: nathanaelle.schneider@cnrs.fr
bLaboratoire de Chimie Moléculaire (LCM), CNRS, école Polytechnique, Institut Polytechnique de Paris, Route de Saclay, 91120 Palaiseau, France
cPSL University, Chimie ParisTech, CNRS, Institut de Recherche Chimie Paris (IRCP), 11, rue Pierre et Marie Curie, F-75231 Paris, France
First published on 26th May 2025
The ability to predict the thermal properties of molecular compounds is essential for their successful integration into vapor-phase processes such as atomic layer deposition (ALD) and chemical vapor deposition (CVD), as well as in catalysis and materials synthesis. In this work, we present a systematic study of 24 gallium amidinate complexes, designed to explore the relationship between molecular structure and thermal behavior. The series encompasses a range of structural variations: different ligand substituents, molecular symmetry, and co-ligands. Structural characterization, including in some cases single-crystal X-ray diffraction, was combined with detailed thermal analysis using thermogravimetric analysis (TGA) and differential scanning calorimetry (DSC) under both atmospheric and reduced pressure conditions. The results reveal clear correlations between thermal properties and ligand architecture, with features such as alkyl chain type, methyl group presence, and overall symmetry playing key roles in determining volatility and stability. Importantly, both symmetric and dissymmetric complexes were found to possess the desired thermal characteristics for vapor-phase deposition processes. Beyond offering valuable design principles for gallium precursors, the dataset generated herein provides a foundation for improving predictive models—empirical and AI-driven alike—towards the rational development of next-generation functional molecular compounds.
Based on this, predicting the thermal properties of molecular compounds is very important and remains a critical challenge. Empirical approaches, based on established structure–property relationships and chemical intuition, have historically guided the development of predictive rules.17–22 More recently, the emergence of artificial intelligence (AI) and machine learning has opened new avenues for forecasting thermal behaviour, offering data-driven models capable of identifying complex, non-linear patterns across diverse chemical families.23–26 Quantitative structure–property relationship (QSPR)27–30 or group-contribution (GC)31–34 approaches being empirical models, they inherently lack predictive accuracy for entirely novel materials that significantly differ from the compounds in the training dataset.35 In all cases, predictions must be anchored in experimental validation—either to confirm model accuracy or to generate robust datasets for further algorithm training and refinement.
To explore the relationship between molecular structure and thermal behaviour and evaluate the validity of empirical structure–property correlations across different structural motifs, we designed a systematic study involving a large set of molecular species. We selected gallium as metal centre and amidinate-based molecules as ligands,36–40 due to the specific interest in this class of compounds to generate gallium-based thin films.41,42 These include gallium complexes with one or two amidinate ligands, featuring various substituents in a symmetric or dissymmetric arrangement, and chloride or methyl as additional ligand. Whenever possible, advanced structural characterization was performed, including single-crystal X-ray diffraction, to provide unambiguous insight into molecular geometry and intermolecular interactions. Thermal properties were assessed both through direct thermal analysis—using thermogravimetric analysis (TGA) and differential scanning calorimetry (DSC) under atmospheric and reduced pressure—and indirectly through their behaviour during purification via distillation or sublimation. These combined data offer a robust basis to examine how specific structural features govern volatility and thermal stability.
C
N–R2) was reacted with one equivalent of the organolithium in diethyl ether at 0 °C and the resulting solid was washed with pentane and characterised by NMR spectroscopy and X-ray crystallography when possible. Nine lithium amidinate complexes were synthesised, with average to good yields (Y = 49 to 93%). To obtain the gallium complexes, the lithium amidinates were reacted with gallium trichloride in a transmetallation reaction. The mono ((Ln)GaCl2) and bis-amidinate ((Ln)2GaCl) complexes were both obtained by controlling the stoichiometry between the lithium amidinate and the gallium salt. The complexes were purified either by sublimation (for solids) or by vacuum distillation (for oils). For this step, the yields obtained varied noticeably (Y = 29 to 90%). Bis-amidinate gallium complexes could not be obtained from the L4Li, L5Li, L7Li and L9Li complexes, even upon heating as confirmed by the persistence of the lithium amidinate signals in the 1H-NMR spectra. This is most likely due to steric effects induced by the bulky ligands that restricts the coordination to one ligand.45–49 Lastly, the gallium methyl complexes were prepared by reaction of (Ln)xGaCl3−x with x equivalents of methyllithium (MeLi). The resulting (Ln)GaMe2 and (Ln)2GaMe complexes were either in oil forms and were purified by vacuum distillation or in solid forms and purified by sublimation ((L1)GaMe2, (L2)2GaMe and (L6)2GaMe). The obtained yields varied over a large range (Y = 13 to 80%), however most of the chlorinated gallium amidinate complexes could be transformed into their methyl counterparts. Note that some of these complexes were previously reported.41,43–45,50–52
The individual yields for all the gallium complexes are reported in Table 1, together with the purification temperatures by sublimation or distillation, which were recorded under a static vacuum of the order of 5 × 10−2 mbar. These temperatures do not appear to depend solely on the molar mass of the complex. The difference between mono- and bis-amidinate is pronounced; the (Ln)GaCl2 complexes have a vaporization temperature below 100 °C, unlike the (Ln)2GaCl complexes, which vaporise above 110 °C. The (L3)2GaCl complex, not being volatile enough, was purified by extraction. The steric hindrance of linear or branched alkyl chains is also an important parameter since (L1)GaCl2 (with CtBu), (L3)GaCl2 (with CnBu) and (L4)GaCl2 (with CsBu), have sublimation temperatures of 100 °C, 80 °C and 80 °C respectively. Fig. 1 shows an overlay of the different 1H-NMR spectra obtained for the complexes with the (L2) ligand. In this example, the ligand protons are characterised by three resonances a quadruplet, a doublet, and a singlet corresponding respectively to the CH and CH3 of the isopropyl groups and that of the endocyclic methyl. It highlights how from one lithium amidinate complex, four gallium amidinate complexes could be synthesised by tuning the number of equivalents of the gallium salt, and how these complexes are easily distinguishable by 1H-NMR analysis. In addition to the different chemical shifts in these five complexes, methyl groups bonded to the gallium led to a more shielded resonance (δ = −0.5–0.5 ppm). These chemical shifts are in agreement with literature precedents.43 More details on all of the reaction steps and characterizations can be found in the ESI.†
| L1Li | L2Li | L3Li | L4Li | L5Li | L6Li | L7Li | L8Li | L9Li | ||
|---|---|---|---|---|---|---|---|---|---|---|
| a For sublimation under a dynamic vacuum of 5 × 10−2 mbar. b For distillation under a dynamic vacuum of 5 × 10−2 mbar. | ||||||||||
|
(Ln)GaCl2 | 90% | 55% | 50% | 82% | 68% | 68% | 29% | 63% | |
| 100 °Ca | 50 °Ca | 80 °Ca | 80 °Ca | 100 °Ca | 85 °Ca | 50 °Cb | 70 °Ca | |||
|
(Ln)2GaCl | 37% | 62% | 62% | — | — | 32% | — | — | — |
| 120 °Ca | 115 °Ca | 150 °Ca | ||||||||
|
(Ln)GaMe2 | 34% | 55% | 29% | 58% | 41% | 49% | 30% | 27% | 50% |
| 50 °Ca | <25 °Cb | 50 °Cb | <25 °Cb | 75 °Cb | <25 °Cb | 40 °Cb | <25 °Cb | <25 °Cb | ||
|
(Ln)2GaMe | 54% | 81% | 13% | ||||||
| 150 °Ca | 90 °Ca | |||||||||
Gallium amidinate complexes were crystallised when possible by sublimation or by slow evaporation of the pure product in a highly volatile solvent such as pentane. The crystal structures with the (L2) amidinate are shown in Table 2, as a representative example. The other structures and more details can be found in the ESI (Tables S1–S6†).
Crystals suitable for single-crystal X-ray diffraction are generally more readily obtained from chlorinated gallium complexes than from their methylated counterparts. While the former typically form solids, the latter are often oily substances that can only be crystalized at low temperatures, as observed for the (L5)GaMe2 complex.
Gallium complexes all crystallise as monomers with amidinate ligands coordinated in κ2-coordination mode. The gallium metal centre is tetracoordinated in mono-amidinate complexes that exhibit a distorted tetrahedral geometry, or pentacoordinated in the case of bis-amidinate complexes. The parameter τ5 is useful to determine the geometry of distorted pentacoordinated complexes giving a value of 0 for a regular square-based pyramidal geometry to and 1 for a regular trigonal bi-pyramidal (see τ5 values in Table 2, and in ESI, Table S6†).53,54(Ln)2GaCl2 complexes present a distorted trigonal bi-pyramidal geometry according to τ5 value closer to 1 than to 0 in all the cases.
Based on its higher value for (L1)2GaCl and (L6)2GaCl complexes, their geometry is closer to a regular trigonal bipyramid that may be explained by the more important steric hindrance on the amidinate compared to the (L2) ligand. The NCN angle of the amidinate ligand also depends on the steric effects, in particular that of the endocyclic alkyl substituent. In general, the larger the substituent, the smaller the angle. This trend is obvious when comparing the (L5) and (L1) complexes with those of (L2) and (L6). Substituting isopropyl with tertbutyl groups at peripheral position seems to have a lower impact. Nevertheless, the NCN angle tends to decrease when increasing the bulkiness of the group (compare (L1)GaCl2 with (L5)GaCl2 and (L2)GaCl2 with (L6)GaCl2)55 (see Table S-6).
Comparing the mono and bis-amidinate complexes, an elongation of the Ga—N bonds (average of 2.05 Å for bis- and of 1.94 Å for mono-ligated complexes) together with a reduction of the bite angles (65.4° on average for (Ln)2GaClvs. 68.4° on average for ((Ln)GaCl2)) are observed. The N—C bond lengths of the amidinate ligand backbone are between 1.31 and 1.36 Å indicating a bond order of 1.5 for all the complexes and in agreement with previous reports.43
The TG/DSC curves under both atmospheric and reduced pressure conditions are provided in the ESI (Fig. S4-1, S4-2 and S4-3†). (Ln)2GaCl complexes were excluded from TG/DSC analyses due to their low volatility, observed during the purification step.
Mass loss profiles varied depending on the analysis conditions (atmospheric vs. reduced pressure) and the type of amidinate complex. The 20 gallium amidinate complexes analysed exhibited distinct thermal profiles, with %res values from −2% to 26%. As expected, TG/DSC analyses under reduced pressure generally resulted in lower residual percentages than those conducted at atmospheric pressure. Most complexes showed residual masses below 10%, except for (L6)2GaMe, (L2)2GaMe, and (L2)GaMe2.
DTG curves were particularly useful in identifying whether a complex vaporised as a single species or as multiple species. Reduced-pressure TG analyses generally led to a single-step vaporization for gallium amidinates, while atmospheric-pressure TG analyses of (Ln)GaCl2 complexes mostly resulted in thermal decomposition, despite the formation of volatile products, as indicated by their low residual masses. Black residue observed in crucibles after analyses further confirmed decomposition under atmospheric pressure conditions. (Ln)GaMe2 complexes demonstrated the most favourable behaviour across all conditions, with minimal decomposition. The only exception was (L2)GaMe2, which showed thermal degradation in the TG curve. The DSC curves of gallium amidinate complexes predominantly showed endothermic peaks: two for most solid complexes, a single peak for oil-like complexes. The endothermic peak attributed to the melting of solid complexes remained unaffected by pressure, as these compounds are low-compressibility materials, a behaviour also observed in cesium hafnium chloride.56 Based on the onset vaporization temperatures at both atmospheric and reduced pressure, methylated complexes exhibited higher volatility than their chlorinated counterparts. The vaporization energies, evaluated by integrating the peak area, do not represent the absolute values for gallium amidinate compounds, as the platinum rod calibration was not performed using a set of commercial compounds under identical experimental conditions (see ESI, Fig. S5†). However, the relative values remain valid, as the energy values were normalised per mole of compound. Thus, reduced-pressure TG/DSC required less energy for vaporization compared to atmospheric pressure. (Ln)GaCl2 complexes exhibited high vaporization energies under atmospheric pressure, likely contributing to multi-species volatilization and black residue formation due to decomposition. (Ln)GaMe2 complexes had the lowest vaporization energies among the tested compounds, indicating higher volatility.
The combination of TG/DSC analyses under different pressures enabled the identification of sublimable solid complexes by comparing their DSC curves. For example, in reduced-pressure TG/DSC of (L1)GaCl2, the presence of a single endothermic peak suggests that the complex undergoes sublimation. However, in atmospheric pressure TG/DSC, two endothermic peaks appear in the DSC signal, indicating a different thermal behaviour.
The purification vaporization temperatures (Tpurif), i.e. the temperature of appearance of a solid on the sublimation finger or the first liquid droplets during distillation, are slightly lower but align with the onset vaporization temperatures (Tonset) determined by TG/DSC. This confirms that volatility estimates from purification are consistent with TG/DSC results. This diversity of ligand and complexes further enables a discussion on the relationship between structure and thermal properties.
The first widely accepted correlation is between vaporization temperature and molecular weight. This trend is confirmed here: the lighter complexes ((Ln)GaMe2) are the most volatile, while the heavier ones ((Ln)2GaMe) exhibit the lowest volatility in our series. However, as shown in Fig. 2, volatility is not solely determined by molecular weight. For example, (L3)GaMe2, (L5)GaMe2 and (L8)GaMe2 complexes are among the heaviest but also show very good volatility.
![]() | ||
| Fig. 2 Influence of the molecular weight on the Tonset. Values, extracted from ATG/DSC analyses under [a] reduced pressure and at [b] atmospheric pressure. | ||
The symmetry of the molecular compound is also known to influence the thermal properties. To show it, we compare complexes with identical molecular weight and branching on the endocyclic carbon but different lateral substituents. Thus, (L3)-based complexes feature two iPr groups, while (L7)-based complexes have tBu and Et substituents, introducing dissymmetry into the complex. Although their onset vaporization temperatures are similar, (L7) exhibits lower residual mass after TG/DSC (%res < 10%), suggesting that the ligand dissymmetry may contribute to an enhanced volatility. Similar observations can be made on (L1)/(L8) and (L2)/(L9)-based complexes. Among all those, increased dissymmetry from lateral substituents appears to improve thermal properties, leading to higher volatility and lower residual mass in chlorinated complexes (see ESI, Fig. S6†). Branching also plays a significant role, as illustrated by the complexes bearing (L1), (L3), and (L4) ligands, which have the same molecular weight but differ in the substituent on the endocyclic carbon (tBu, nBu, and sBu). According to reduced-pressure TG/DSC, (L1)GaCl2 is the only complex displaying a sublimation endothermic peak (Ts = 137 °C), while the others exhibit a melting peak followed by a vaporization peak. Among their methylated counterparts, (L4)GaMe2 appears to be the most volatile, whereas (L1)GaMe2 shows the fastest vaporization rate. For this series of six complexes, the presence of a bulky tBu group on the endocyclic carbon leads to greater decomposition (see ESI, Fig. S7†). The comparison of (L8) and (L7)-based complexes further evaluates the combined effect of endocyclic substituent branching and symmetry. While (L8) and (L7)-based complexes show a large difference in onset vaporization temperature, their TG profiles remain similar (see ESI, Fig. S8†). These findings suggest that introducing bulky substituents on the endocyclic carbon in a dissymmetric complex enhances volatility, making these structures particularly relevant for applications requiring volatile precursors.
Finally, our modular synthesis approach enabled access to up to four gallium amidinate complexes for a given ligand. Due to the low volatility of bis-amidinate chlorinated complexes, only the three methylated derivatives were analysed by TG/DSC for the (L2), (L3) and (L6) amidinate ligands. Interestingly, (L2)2GaMe and (L2)GaCl2 exhibit similar vaporization temperatures despite a significant molecular weight difference (ΔM = 85 g mol−1). However, (L2)2GaMe appears to decompose, as indicated by changes in the TG slope between 100 and 125 °C under reduced-pressure TG/DSC. Similar observations can be made with (L3)-based complexes. Notably, bis-amidinate complexes seem to undergo partial modification under reduced pressure, whereas they remain stable at atmospheric pressure (see ESI, Fig. S9†). This highlights the influence of pressure conditions on their thermal behaviour, which should be systematically investigated for ALD/CVD precursor development.
Our findings highlight the significant influence of ligand architecture, including the nature of substituent chains, the presence of methyl groups, and overall molecular symmetry, on thermal behavior. Both symmetric and dissymmetric complexes yielded compounds with high volatility and stability, underscoring the versatility of the amidinate ligand framework for tailoring thermal properties.
This systematic study not only offers new insights into structure–thermal property relationships but also provides a valuable dataset for refining predictive models, including empirical and AI-based approaches. Such combined experimental and data-driven strategies will be key to accelerating the discovery and rational design of the next-generation molecular precursors.
Footnote |
| † Electronic supplementary information (ESI) available. CCDC 2444176–2444185. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d5dt00944h |
| This journal is © The Royal Society of Chemistry 2025 |