Abstract
Introduction:
Traditionally, the condensation sequence of circumstellar dust is predicted based on the thermodynamic stabilities of specific condensates in the macroscopic bulk phase. However, at the (sub-) nanometer scale clusters with non-crystalline structures and significantly different properties are energetically favoured.
Methods:
For this reason, we study the thermodynamic stabilities of metal oxide clusters with generic stoichiometries of M2O3 and M3O4, where M represents a metal atom. With an upper size limit of 50 atoms, we consider clusters with sizes n = 1–10 for (M2O3)n, and n = 1–7 for (M3O4)n. The M2O3 clusters comprise alumina (Al2O3), Mg-rich pyroxene (MgSiO3) and a size-limited sample of titanates (CaTiO3), whereas the M3O4 clusters include spinel (MgAl2O4), Mg-rich olivine (Mg2SiO4) and calcium aluminates (CaAl2O4).
Results:
We find that, apart from the alumina monomer, the aluminum-bearing clusters (Al2O3)n, n = 1–10, and (MgAl2O4)n, n = 1–7, are favoured over their silicate counterparts (MgSiO3)n, n = 1–10 and (Mg2SiO4)n, n = 1–7. Also, we find that calcium aluminate clusters, CaAl2O4, are energetically more favourable than magnesium aluminate clusters, MgAl2O4. Furthermore, for a limited data set of (CaTiO3)n, n = 1–2, clusters we find significantly larger stabilities than for the other considered (M2O3)n clusters, namely Al2O3 and MgSiO3.
Discussion:
Future investigations, in particular on titanates and on Ca-rich silicates, are required to draw a more thorough and complete picture of the condensation sequence at the (sub-)nanoscale.
1 Introduction
Asymptotic Giant Branch (AGB) stars with their highly dynamical atmospheres and circumstellar envelopes are unique astrochemical laboratories (). Shell burning and mixing processes inside of AGB stars lead to a changing elemental composition at the atmosphere over time (). In AGB circumstellar envelopes the presence of large-amplitude stellar pulsations, active dust formation from the gas-phase, and maser emission indicate departures from thermodynamic equilibrium conditions, which require non-equilibrium modelling ().
In equilibrium, the chemistry is controlled by the C/O ratio, which is a consequence of the triple bonded CO molecule corresponding to the diatomic species with the highest bond energy (11.3 eV) known. As a result, an oxygen-rich chemistry with C/O < 1 and carbon-rich environments with C/O > 1 are expected, as the CO molecule locks the lesser abundant element (C or O). The observationally confirmed presence of HCN and CS in oxygen-dominated atmospheres (; ), as well as the detection of H2O vapour and SiO in Carbon-rich envelopes challenge this traditional dichotomy (; ). With regard to the dust populations the situation is different. To the best of our knowledge there is no confirmed active carbonaceous dust formation in oxygen-dominated AGB stellar atmospheres and no oxide dust synthesis in carbon-rich AGB envelopes.
Circumstellar dust species in oxygen-rich conditions comprise minerals and amorphous solids that are made of metal oxides with the stoichiometric formula M2O3 including corundum (Al2O3), enstatite (MgSiO3), perovskite (CaTiO3). This is also true for metal oxides with the generic formula M3O4 that includes MgAl2O4 (spinel), CaAl2O4 (krotite), Mg2SiO4 (forsterite). Evidence for the presence of alumina and spinel is found in spectral features observed at around 13 microns (; ), whereas the silicates, here enstatite and forsterite, show Si-O stretching modes at 10 micron and O-Si-O bending modes at around 18 microns ().
Iron-bearing solid oxides like hematite (Fe2O3), ferrosilit (FeSiO3) or fayalite (Fe2SiO4) fall also in the M2O3 and M3O4 stoichiometric families, but are not expected to act as seed particles for dust formation, owing to the too large opacities of the Fe-bearing seeds ().
Traditionally, the thermodynamic stabilities of condensates are assessed based on a top-down approach starting with the macroscopic crystalline structure of a specific mineral. In classical nucleation theories (CNTs) the cohesive energy of a microscopic particle like a cluster is determined by an attractive volume term and a repulsive surface term resulting in a spherical particle as lowest-energy geometry. Depending on the material in question, top-down approaches like CNTs can predict reliable particle energies down to typically tens of nanometers, in rare cases even down to a few nanometers.
However, at the nanoscale and below, the lowest energy structures are significantly different from top-down derived geometries. Owing to finite size and quantum effects the properties of the small clusters including atomic coordination, bond lengths, formal charges, dipole moments are notably different from the macroscopic bulk phase. Notably, the most favourable cluster structures are non-crystalline. Therefore, the drawbacks of CNTs at the (sub-)nanoscale comprise the spherical cluster structures, fully coordinated atoms in the homogeneous interior of the cluster, growth by monomeric additions, bimolecular association reactions, and most prominently, unrealistic potential energies.
In this study we will focus on the thermochemistry of bottom-up generated, (sub-)nanometer sized clusters as precursors of dust grains in oxygen-rich circumstellar atmospheres and envelopes. Apart from Al2O3 the clusters presented in this study are ternary oxides comprising two different metal elements. In contrast, binary oxides contain just one type of metal and are therefore chemically and structurally less complex. Still, it is a challenging task to find the lowest-energy isomers referred to as global minima (GM) structures in the case of binary oxide clusters.
Smaller gas phase dust precursors in oxygen-rich environments include simple diatomic and triatomic molecules with well-defined spectroscopic constants, rotational and vibrational transitions. In the case of silicates, a very likely molecular precursor is SiO (), whereas the aluminates including alumina are preceded by AlO and AlOH (; ; ). Titanantes originate from TiO and TiO2 molecules showing relatively strong Ti-O bonds (; ). Moreover, the hydroxyl radical, OH, and water vapour, H2O, are considered to act as the most efficient oxidizer of the metal oxide molecules and clusters ().
In this study we will focus on the thermochemistry of (sub-)nanometersized clusters with M2O3 and M3O4 stoichiometries as precursors of dust grains in oxygen-rich circumstellar atmospheres and envelopes. Moreover, this study aims at overcoming the drawbacks of classical nucleation theories by employing a bottom-up approach using lowest-energy cluster configuration with the most accurate yet affordable quantum calculations, instead of using top-down generated spherical particle structures.
This paper is organized as follows. In Section 2, we describe the methods used to compute the most favourable (M2O3)n and (M3O4)n cluster structures. The results are presented in Section 3. Section 4 discusses our results in the light of previous studies and summarizes our findings.
2 Methods
The energetically most favourable cluster isomers, or GM candidates, of the considered dust precursors of alumina, spinel, Ca aluminates, Mg-rich silicates and calcium titanates were recently investigated. For our purposes we consult the GM candidates for (Al2O3)n, n = 1–10 (), MgAl2O4, CaAl2O4 (), (MgSiO3)n, n = 1–10, (Mg2SiO4)n, n = 1–7 (), and (CaTiO3)n, n = 1,2 ().
A comparison among the trioxides Al2O3, MgSiO3, and CaTiO3 is straightforward as these cluster contain the same number of atoms, metal atoms and oxygen atoms, respectively, per formula unit or cluster size n. The same reasoning is true for a comparison among clusters of MgAl2O4, CaAl2O4 and Mg2SiO4, showing seven atoms per formula unit of which three are metals and four are oxygens. For consistency we apply the B3LYP density functional in combination with a cc-pVTZ basis set including a vibrational frequency analysis for all calculations presented in this study (). This density functional basis set combination was chosen as a sensible compromise between computational cost and desired accuracy. Moreover, this combination has shown a reasonable agreement with experimental results for transition metal oxides such as titania and vanadia (; ).
We use the RRHO (Rigid Rotor Harmonic Oscillator) approximation to compute the partition functions of the cluster leading the thermodynamic potential of interest, namely the enthalpy of formation, entropy and Gibbs Free energy of formation. For consistency, all presented quantum-chemical density functional calculations were computed with the Gaussian16 programme suite ().
3 Results
In Figure 1 the thermodynamic stabilities represented by the normalized Gibbs free energy Δf G(T) of formation of the (M2O3)n clusters are shown as a function of the cluster size n. The energies of CaTiO3 clusters are colour-coded in orange, those of Al2O3 clusters in purple and those of MgSiO3 clusters in green, respectively. Solid lines correspond to a temperature of T = 0 K, dash-dotted lines to a temperature of T = 1000 K, and dashed lines to a temperature of T = 2000 K.
FIGURE 1
Clearly, the (CaTiO3)n family is most favourable with regard to the studied M2O3 clusters for all considered temperatures. Alumina clusters have stabilities that are more than 200 kJ/mol less per formula unit than the calcium titanates, but are more favourable than the Mg-rich pyroxene clusters by about 100–150 kJ/mol, except for the monomer. The monomers Al2O3 and MgSiO3 show almost identical Gibbs free energies at a temperature of 2000 K. Moreover, we note a comparatively lower stability for (Al2O3)n clusters at sizes n = 8 and 10 and for (MgSiO3)n at size n = 5.
In Figure 2 the normalised Gibbs free energies Δf G(T) of formation of the (M3O4)n clusters are shown as a function of the cluster size n. The energies of CaAl2O4 clusters are colour-coded in red, those of MgAl2O4 clusters in purple and those of Mg2SiO4 clusters in green, respectively. As in Figure 1, solid lines correspond to a temperature of T = 0 K, dash-dotted lines to a temperature of T = 1000 K, and dashed lines to a temperature of T = 2000 K. The calcium aluminate clusters constitute the most favourable clusters among the M3O4 set for all selected temperatures. The energy difference to their Mg-rich counterpart clusters, (MgAl2O4)n, is in the range of ∼150–200 kJ/mol. This energy difference is larger than the difference of ∼100–150 kJ/mol between the spinel clusters, (MgAl2O4)n, to the less favourable Mg-rich olivine silicate clusters. For cluster size n = 5 we find a comparatively low stability of the aluminate clusters (MgAl2O4)5 and (CaAl2O4)5.
FIGURE 2
Note that the Gibbs free energies of formation in Figures 1, 2 are not scaled to the atomic heats of formation as in the JANAF tables (). In this way, the enthalpy of formation Δf H (0 K) and Gibbs free energy of formation Δf H (0 K) correspond to the unscaled total binding energy of the cluster at 0 K. The corresponding thermochemical tables in the standard JANAF format can be found in the Supplementary Appendix A1. Note that the thermochemical tables of the cluster families (MgAl2O4)n and (CaAl2O4)n can be found in the supporting information of .
An illustration of the applicability of the calculated Gibbs free energies of formation is shown in the form of chemical equilibrium abundances in the Supplementary Appendix A3. The corresponding computations are performed using the software GGchem () at a pressure of 1 bar using a solar elemental composition (). The clusters of the M2O3 and M3O4 stoichiometries were treated separately in order to avoid interference between these families. In the M2O3 family, clusters larger than the monomer (n > 1) are present in the form of (CaTiO3)2 for temperatures of T ≤ 1500 K, (Al2O3)10 for T ≤ 1300 K, and (MgSiO3)10 for T ≤ 1200 K. The M3O4 family shows no significant amount of (MgAl2O4)n clusters, which can be explained by the presence of the thermodynamically favoured, Al-bearing (CaAl2O4)n clusters at T ≤ 1400 K and the Mg-bearing silicate clusters (Mg2SiO4)n at T ≤ 1200 K. These equilibrium abundances reflect the results presented in Figures 1 and 2 weighted by elemental abundances and indicate trends in a potential sequence of condensation.
However, we emphasize that the equilibrium abundances shown here are not representative for the actual existence or presence of the corresponding clusters in stellar atmospheres, since chemical equilibrium cannot trace the complex chemical kinetics of cluster formation. For example, formation routes that are prohibited by energy barriers or spin cannot be identified or predicted in chemical equilibrium considerations.
4 Discussion and summary
Alumina and ternary aluminium-bearing oxides are thermodynamically favoured over their silicate counterparts. Furthermore, calculations constrained to small cluster sizes indicate that titanium-bearing ternary oxides are even more favourable compared with aluminates and silicates.
The thermodynamic stabilities of the cluster families seem to follow an opposite trend with respect to the elemental abundances of their constituent metals, where Ti is about an order less abundant than Al and about two orders of magnitude less abundant than Si, according to solar abundances () and AGB stellar evolution models (). This trend is also seen, when comparing Mg-containing clusters with the corresponding Ca-bearing counterparts, which are more favourable and Ca is about an order of magnitude less abundant than Mg.
Apart from alumina (Al2O3) the investigated M2O3 and M3O4 clusters are ternary metal oxides comprising two different metals. A comparison with binary oxide clusters such as titania (TiO2), silica (SiO2), and magnesia (MgO) is not included in our study, since their metal-to-oxygen ratio does not agree with a M2O3 or a M3O4 stoichiometry, respectively, making such comparisons not straight-forward. Moreover, the listed binary oxide clusters can be regarded as an integral part of the presented M2O3 and M3O4 clusters. TiO2 is contained in CaTiO3, whereas SiO2 and MgO represent formal constituents of MgSiO3 and Mg2SiO4. So, albeit clusters of TiO2, SiO and MgO are tightly linked to the ternary oxide clusters of the presented study, they are not explicitly part of this study. However, we briefly discus their energetic and kinetic viabilities as well as limitations in AGB circumstellar envelopes with respect to homogeneous nucleation.
SiO is an abundant molecule in all chemical types of AGB stars (see e.g. ). However, the SiO nucleation via (SiO)n clusters was found to be inefficient and negligible under circumstellar conditions (). Moreover, with increasing size the most favourable (SiO)n cluster exhibit segregated Si islands, further impeding the cluster growth by monomeric SiO additions. It should also be noted here that macroscopic solid SiO quickly segregates into islands of amorphous silica (SiO2) and silicon (). Therefore, SiO cannot be regarded as a stable condensate. However, at the (sub-)nanoscale (SiO)n clusters can exist. SiO is a key molecule for the formation of silicates. explored a mechanism, where the SiO molecule is dimerized to Si2O2, which subsequently oxidized and enriched with Mg atoms in an alternating manner to form Mg4Si2O8, the dimer of Mg-rich olivine. The dimerisation of SiO represents the energetic bottleneck in their scheme at 1000 K, whereas the formation of the trimer, Si3O3, could proceed via SiO2 and Si2O3 in terms of energy. However, the oxidation of SiO to SiO2 is kinetically hampered and also the dimerization of SiO2 to Si2O4 represents a bottleneck reaction in circumstellar envelopes as has been confirmed by .
To our knowledge MgO does not exist in AGB circumstellar atmospheres, and Mg is predominantly expected to be in atomic form. Even if the MgO monomers were present in circumstellar envelopes, its homogeneous nucleation is hampered by particularly favourable “magic” cluster sizes that act as bottlenecks for cluster growth (; ). For Mg-bearing titanates and aluminates the inclusion of magnesium constitutes a kinetic bottleneck (; ). In these studies, it was also found that reactions with calcium atoms are significantly more efficient than those with Mg atoms and can, under certain cirumstances, lead to the formation of condensation seeds.
TiO2 exists as a strongly bond gas phase molecule, i.e. as a monomer, but also in a macroscopic crystalline phase as anatase or rutile making it a prime candidate for nucleation. Its nucleation path via stoichiometric (TiO2)n clusters follows indeed an energetically downhill process (; ). However, a chemical-kinetic assessment of the dimerization of TiO2 to (TiO2)2, representing the first reaction in this scenario, shows a slow and ineffective reaction rate with a negative temperature-dependence (). It implies that the rate is very slow at high temperatures and requires a third body to proceed. Moreover, homogeneous TiO2 nucleation is limited by the availability of Ti having a low abundance.
Formation routes towards the alumina monomer, Al2O3, are thermodynamically and kinetically hampered. The synthesis of the triplet Al2O3 monomer synthesis is prohibited by an unfavourable, i.e. an endothermic oxidation of Al2O2 (). Instead, a viable pathway leading to the formation of the alumina dimer, (Al2O3)2, was identified. Alternative solutions to form (Al2O3)n clusters might involve Al(OH)3 as suggested by .
In summary, SiO nucleation is hampered by increasing atomic segregation and slow growth rates, MgO nucleation is hindered by the absence of the MgO monomer and energy bottlenecks at magic cluster sizes, and homogeneous TiO2 nucleation is viable, but constrained by the overall low Ti abundance. Given the chemical wealth of oxygen-rich AGB circumstellar envelopes, it is likely that ternary oxide clusters form. In particular, the above mentionned drawbacks of homogeneous nucleation of SiO, MgO and TiO2, respectively, indicate that the formation of a ternary oxide involving an additional metal is required for subsequent cluster growth.
For a comparison between the binary oxide clusters of SiO, TiO2, MgO and Al2O3 we refer to the study of . Although a (chemical-)kinetic treatment of the studied clusters would be desirable, an investigation of the presented clusters is computationally very expensive or even prohibitive for cluster sizes n > 3. The reaction rates involving larger sized clusters are very expensive to compute at high level rate theories like statistical Ramsperger-Kassel-Marcus (RRKM) and transition state theories. For this reason, reaction systems comprising more than ∼12 atoms the rates are typically assessed with theories of approximative nature including classical and kinetic nucleation, collisional rates, capture rates and detailed balance.
Spectral identification of the presented clusters is challenging. Owing to their large number of degrees of freedom the presented clusters exhibit many, partly overlapping and blurred spectral features in contrast to small, di- and tri-atomic molecules with sharp, discrete and well-defined transitions. In addition, the nucleating clusters might not have a longevity to be observed. However, recent studies have attempted to identify the spectral signatures of different cluster families in theoretical frameworks (; ; ; ; ) as well as observationally (; ).
Calcium silicate clusters, in particular the pyroxenes and olivines (CaSiO3)n and (Ca2SiO4)n, are entirely missing for both stoichiometries, M2O3 and M3O4, which constitutes one of the major limitations of the present study. Moreover, with the exception of the monomer and the dimer of perovskite (CaTiO3), the calcium titanates clusters are not part of this study. We aim to investigate the cluster families of calcium silicates and titanates in a forthcoming study.
Therefore, no final conclusion can be drawn. However, our results indicate the following trends:
(i) alumina and aluminate clusters are energetically more favourable than their silicate counterparts
(ii) Owing to their thermodynamic stability calcium titanates are promising candidates as primary condensation seeds
(iii) Clusters comprising calcium are significantly more stable than their magnesium-containing counterparts
Other stoichiometric clusters of minerals like geikielite (MgTiO3) and quandilite (Mg2TiO4) were not considered in this study. These condensates are rare and were, to the best of our knowledge, not found in meteorites. For consistency and completeness, however, these clusters as well as the aforementioned calcium silicates CaSiO3 and Ca2SiO4 should be investigated in detail to draw firmer conclusions on the cluster condensation sequence among titanates, aluminates and silicates.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
DG: Writing – review and editing, Writing – original draft, Formal Analysis, Conceptualization, Validation, Data curation, Investigation.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. The computations involved the Swedish National Infrastructure for Computing (SNIC) at Chalmers Centre for Computational Science and Engineering (C3SE) partially funded by the Swedish Research Council through grant no. 2018-05973. The computations were partially enabled by resources provided by the National Academic Infrastructure for Supercomputing in Sweden (NAISS) partially funded by the Swedish Research Council through grant agreement no. 2022-06725.
Acknowledgments
I am deeply indebted to Tom Millar and Mauro Pirarba for enabling a publication in frontiers.
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The handling editor declared a past co-authorship with one of the authors DG.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fspas.2025.1632593/full#supplementary-material
References
1
AnderssonS.GobrechtD.ValeroR. (2023). Mechanisms of SiO oxidation: implications for dust formation. Front. Astron. Space Sci.10, 1135156. 10.3389/fspas.2023.1135156
2
AsplundM.GrevesseN.SauvalA. J.ScottP. (2009). The chemical composition of the sun. ARA&A47, 481–522. 10.1146/annurev.astro.46.060407.145222
3
BaeyensR.BaratS.DecinL.DesertJ. M.GobrechtD.HellingCh.et al (2024). JWST proposal. Cycle 3, ID. #6045. Available online at: https://ui.adsabs.harvard.edu/abs/2024jwst.prop.6045B/abstract.
4
BaudryA.WongK. T.EtokaS.RichardsA. M. S.MüllerH. S. P.HerpinF.et al (2023). ATOMIUM: probing the inner wind of evolved O-rich stars with new, highly excited H2O and OH lines. A&A674, A125. 10.1051/0004-6361/202245193
5
BeckeA. D. (1993). J. Chem. Phys.98, 1372. 10.1063/1.464913
6
BhattJ. S.FordJ. (2007). Investigation of MgO as a candidate for the primary nucleating dust species around M stars. MNRAS382 (1), 291–298. 10.1111/j.1365-2966.2007.12358.x
7
BoulangierJ.GobrechtD.DecinL.de KoterA.YatesJ. (2019). Developing a self-consistent AGB wind model – II. Non-classical, non-equilibrium polymer nucleation in a chemical mixture. MNRAS489, 4890–4911. 10.1093/mnras/stz2358
8
BromleyS. T.Gómez MartínJ. C.PlaneJ. M. C. (2016). Phys. Chem. Chem. Phys.18, 26913–26922. 10.1039/C6CP03629E
9
ChaseM. (1998). NIST-JANAF thermochemical tables. 4th Edition. American Institute of Physics. Available online at: https://www.nist.gov/publications/nist-janaf-thermochemical-tables-4th-edition.
10
CristalloS.StranieroO.PiersantiL.GobrechtD. (2015). Evolution, nucleosynthesis, and yields of agb stars at different metallicities. iii. intermediate-mass models, revised low-mass models, and the ph-fruity interface. ApJS219, 40. 10.1088/0067-0049/219/2/40
11
DanilovichT.GottliebC. A.DecinL.RichardsA. M. S.LeeK. L. K.KamińskiT.et al (2020). Rotational spectra of vibrationally excited AlO and TiO in oxygen-rich stars. ApJ904, 110. 10.3847/1538-4357/abc079
12
DecinL.De BeckE.BrünkenS.MüllerH. S. P.MentenK. M.KimH.et al (2010). A&A516, A69. 10.1051/0004-6361/201014136
13
DecinL.RichardsA. M. S.WatersL. B. F. M.DanilovichT.GobrechtD.KhouriT.et al (2017). Study of the aluminium content in AGB winds using ALMA: indications for the presence of gas-phase (Al2O3)nclusters. A&A608, A55. 10.1051/0004-6361/201730782
14
FirthR. A.BellK. M.FortenberryR. C. (2025). Formation of AlO, AlOH, and Al(OH)3 in the interstellar medium and circumstellar envelopes of AGB stars. ACS Earth Space Chem.8 (5), 974–982. 10.1021/acsearthspacechem.3c00335
15
FonfriaJ. P.Fernandez-LopezM.AgundezM.Sánchez-ContrerasC.CurielS.CernicharoJ. (2014). The complex dust formation zone of the AGB star IRC+10216 probed with CARMA 0.25 arcsec angular resolution molecular observations. Mon. Not. R. Astron. Soc.445, 3289–3308. 10.1093/mnras/stu1968
16
FriedeB.JansenM. (1996). Some comments on so-called ‘silicon monoxide’. J. Non-Crystalline Solids204 (2), 202–203. 10.1016/S0022-3093(96)00555-8
17
FrischM. J.TrucksG. W.SchlegelH. B.ScuseriaG. E.RobbM. A.CheesemanJ. R.et al (2016). Gaussian’16 revision C.01. Wallingford, CT: Gaussian Inc.Available online at: https://gaussian.com/citation/.
18
GobrechtD.CherchneffI.SarangiA.PlaneJ. M. C.BromleyS. T. (2016). Dust formation in the oxygen-rich AGB star IK tauri. A&A585, A6. 10.1051/0004-6361/201425363
19
GobrechtD.PlaneJ. M. C.BromleyS. T.DecinL.CristalloS.SekaranS. (2022). Bottom-up dust nucleation theory in oxygen-rich evolved stars: I. Aluminium oxide clusters. A&A658, A167. 10.1051/0004-6361/202141976
20
GobrechtD.HashemiS. R.PlaneJ. M. C.BromleyS. T.NymanG.DecinL. (2023). Bottom-up dust nucleation theory in oxygen-rich evolved stars: II. Magnesium and calcium aluminate clusters. A&A680, A18. 10.1051/0004-6361/202347546
21
GoumansT. P. M.BromleyS. T. (2012). MNRAS420, 3344. 10.1111/j.1365-2966.2011.20255.x
22
HerwigF. (2005). Evolution of asymptotic giant branch stars. ARA&A43, 435–479. 10.1146/annurev.astro.43.072103.150600
23
HöfnerS.OlofssonH. (2018). A&ARv26, 1. 10.1007/s00159-017-0106-5
24
KamińskiT.WongK. T.SchmidtM. R.MüllerH. S. P.GottliebC. A.CherchneffI.et al (2016). An observational study of dust nucleation in mira (O Ceti): I. Vaciable features of AlO and other Al-bearing species⋆. A&A592, A42. 10.1051/0004-6361/201628664
25
KamińskiT.MüllerH. S. P.SchmidtM. R. (2017). A&A599, A59. 10.1051/0004-6361/201629838
26
KöhlerT. M.GailH. P.SedlmayrE. (1997). A&A320, 553. Available online at: https://adsabs.harvard.edu/full/1997A%26A.320.553K.
27
Lamiel-GarciaO.KoK. C.LeeJ. Y.BromleyS. T.IllasF. (2017). J. Chem. Theory Comput.13, 4. 10.1021/acs.jctc.7b00085
28
Lecoq-MolinosH.GobrechtD.SindelJ. P.DecinL. (2024). Vanadium oxide clusters in substellar atmospheres: a quantum chemical study. A&A690, A34. 10.1051/0004-6361/202347693
29
Macià-EscatllarA.LazaukasT.WoodleyS. M.BromleyS. T. (2019). Structure and properties of nanosilicates with olivine (Mg2SiO4)N and pyroxene (MgSiO3)N compositions. ACS Earth Space Chem.3, 2390–2403. 10.1021/acsearthspacechem.9b00139
30
Mariñoso GuiuJ.Macià-EscatllarA.BromleyS. T. (2021). ACS Earth and. Space Chem.5, 4. 10.1021/acsearthspacechem.0c00341
31
NeufeldD. A.TollsV.AgúndezM.González-AlfonsoE.DecinL.DanielF.et al (2013). ApJL767, 1–5. 10.1088/2041-8205/767/1/L3
32
NuthJ. A.DonnB. (1981). Vibrational disequilibrium in low pressure clouds. ApJ247, 925. 10.1086/159101
33
PlaneJ. M. C. (2013). Philos. Trans. Roy. Soc. A: math. Phys. Eng. Sci.371, 20120335. 10.1098/rsta.2012.0335
34
PlaneJ. M. C.RobertsonS. (2022). Master equation modelling of non-equilibrium chemistry in stellar outflows. Faraday Discuss.238, 461–474. 10.1039/D2FD00025C
35
PoschT.KerschbaumF.MutschkeH.FabianD.DorschnerJ.HronJ. (1999). A&A352, 609. Available online at: https://ui.adsabs.harvard.edu/abs/1999A%26A.352.609P/abstract.
36
RamstedtS.SchöierF. L.OlofssonH. (2009). A&A499, 2–515. 10.1051/0004-6361/200911730
37
SchöierF. L.RamstedtS.OlofssonH.LindqvistM.BiegingJ. H.MarvelK. B. (2013). The abundance of HCN in circumstellar envelopes of AGB stars of different chemical type. A&A550, A78. 10.1051/0004-6361/201220400
38
SindelJ. P.GobrechtD.HellingC.DecinL. (2022). Revisiting fundamental properties of TiO2 nanoclusters as condensation seeds in astrophysical environments. A&A668, A35. 10.1051/0004-6361/202243306
39
SindelJ. P.HellingC.GobrechtD.ChubbK. L.DecinL. (2023). Infrared spectra of TiO2 clusters for hot Jupiter atmospheres. A&A680, A65. 10.1051/0004-6361/202346347
40
SloanG. C.KraemerK. E.GoebelJ. H.PriceS. D. (2003). Guilt by association: the 13 micron dust emission feature and its correlation to other gas and dust features. ApJ594, 483–495. 10.1086/376857
41
WoitkeP. (2006). Too little radiation pressure on dust in the winds of oxygen-rich AGB stars. A&A460, L9–L12. 10.1051/0004-6361:20066322
42
WoitkeP.HellingC.HunterG. H.MillardJ. D.TurnerG. E.WortersM.et al (2018). Equilibrium chemistry Down to 100 K: impact of silicates and phyllosilicates on the carbon to oxygen ratio. A&A614, A1–A28. 10.1051/0004-6361/201732193
43
WoolfN. J.NeyE. P. (1969). Circumstellar infrared emission from cool stars. ApJ155, L181. 10.1086/180331
Summary
Keywords
nucleation, clusters, dust, circumstellar, metal oxides, silicates, alumina
Citation
Gobrecht D (2025) Condensation sequence of circumstellar cluster seeds (CSCCS). Front. Astron. Space Sci. 12:1632593. doi: 10.3389/fspas.2025.1632593
Received
21 May 2025
Accepted
15 September 2025
Published
10 October 2025
Volume
12 - 2025
Edited by
Xiaohu Li, Chinese Academy of Sciences, China
Updates
Copyright
© 2025 Gobrecht.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: David Gobrecht, david.gobrecht@ipi.ch
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.