-
Clearing the Air: Updated CS Fluorescence Efficiencies Don't Solve the NUV/Sub-mm Cometary CS Abundance Discrepancy
Authors:
Steven Bromley,
John Noonan,
Barbora Stachova,
Juraj Orszagh,
Dennis Bodewits
Abstract:
Observations of carbon monosulfide (CS), a possible photodissociation fragment of CS$_2$, has a long history serving as a remote proxy for atomic sulfur, and more broadly, one of the sulfur reservoirs in cometary bodies. Interpretations of CS fluorescence emissions in ultraviolet observations of comets have, to-date, relied on a murky and ill-documented lineage of calculations for CS fluorescence…
▽ More
Observations of carbon monosulfide (CS), a possible photodissociation fragment of CS$_2$, has a long history serving as a remote proxy for atomic sulfur, and more broadly, one of the sulfur reservoirs in cometary bodies. Interpretations of CS fluorescence emissions in ultraviolet observations of comets have, to-date, relied on a murky and ill-documented lineage of calculations for CS fluorescence efficiencies that can be difficult to assess. We report new fluorescence efficiencies of the CS radical, utilizing a rovibrational structure with vibrational states up to $v = 8$ and rotational states up to $N = 100$. A new set of band transition rates is produced through benchmarking comparisons to laboratory electron impact experiments, and subsequently utilized in the fluorescence calculations. We confirm historical reports of the ultraviolet $(0,0)$ band fluorescence efficiency, and conclude that blends and an expanded rovibrational structure have little impact on the fluorescence properties. Combined, the present results affirm the accuracy of the historical record of sulfur abundances derived via ultraviolet CS emissions in cometary observations with e.g., IUE and HST. Benchmark comparisons to IUE observations of comet Bradfield show favorable agreement with the theoretical models. We expand the model to the time-domain, and find that thermalization of the ground state and/or `hot' initial distributions of CS lead to noticable changes in the shape of the $(0,0)$ band. The impact on existing sulfur abundances derived via ultraviolet CS observations, and new connections to sulfur reservoirs in protoplanetary disks are discussed. The model code and computed fluorescence efficiencies are made publicly available.
△ Less
Submitted 27 June, 2024;
originally announced June 2024.
-
Bottom-up dust nucleation theory in oxygen-rich evolved stars II. Magnesium and calcium aluminate clusters
Authors:
David Gobrecht,
S. Rasoul Hashemi,
John M. C. Plane,
Stefan T. Bromley,
Gunnar Nyman,
Leen Decin
Abstract:
Spinel (MgAl$_{2}$O$_{4}$) and krotite (CaAl$_{2}$O$_{4}$) are alternative candidates to alumina (Al$_2$O$_3$) as primary dust condensates in the atmospheres of oxygen-rich evolved stars. Moreover, spinel was proposed as a potential carrier of the circumstellar 13 $μ$m feature. However, the formation of nucleating spinel clusters is challenging; in particular, the inclusion of Mg constitutes a kin…
▽ More
Spinel (MgAl$_{2}$O$_{4}$) and krotite (CaAl$_{2}$O$_{4}$) are alternative candidates to alumina (Al$_2$O$_3$) as primary dust condensates in the atmospheres of oxygen-rich evolved stars. Moreover, spinel was proposed as a potential carrier of the circumstellar 13 $μ$m feature. However, the formation of nucleating spinel clusters is challenging; in particular, the inclusion of Mg constitutes a kinetic bottleneck. We aim to understand the initial steps of cosmic dust formation (i.e. nucleation) in oxygen-rich environments using a quantum-chemical bottom-up approach. Starting with an elemental gas-phase composition, we constructed a detailed chemical-kinetic network that describes the formation and destruction of magnesium-, calcium-, and aluminium-bearing molecules as well as the smallest dust-forming (MgAl$_{2}$O$_{4}$)$_1$ and (CaAl$_{2}$O$_{4}$)$_1$ monomer clusters. Different formation scenarios with exothermic pathways were explored, including the alumina (Al$_2$O$_3$) cluster chemistry studied in Paper I of this series. The resulting extensive network was applied to two model stars, a semi-regular variable and a Mira-type star, and to different circumstellar gas trajectories, including a non-pulsating outflow and a pulsating model. We employed global optimisation techniques to find the most favourable (MgAl$_2$O$_4$)$_n$, (CaAl$_2$O$_4$)$_n$, and mixed (Mg$_x$Ca$_{(1-x)}$Al$_2$O$_4$)$_n$ isomers, with $n$=1$-$7 and x$\in$[0..1], and we used high level quantum-chemical methods to determine their potential energies. The growth of larger clusters with $n$=2$-$7 is described by the temperature-dependent Gibbs free energies.In the considered stellar outflow models, spinel clusters do not form in significant amounts. However, we find that in the Mira-type non-pulsating model CaAl$_2$O$_3$(OH)$_2$, a hydroxylated form of the calcium aluminate krotite monomer forms ...
△ Less
Submitted 12 October, 2023;
originally announced October 2023.
-
A machine learning potential for simulating infrared spectra of nanosilicate clusters
Authors:
Zeyuan Tang,
Stefan T. Bromley,
Bjørk Hammer
Abstract:
The use of machine learning (ML) in chemical physics has enabled the construction of interatomic potentials having the accuracy of ab initio methods and a computational cost comparable to that of classical force fields. Training an ML model requires an efficient method for the generation of training data. Here we apply an accurate and efficient protocol to collect training data for constructing a…
▽ More
The use of machine learning (ML) in chemical physics has enabled the construction of interatomic potentials having the accuracy of ab initio methods and a computational cost comparable to that of classical force fields. Training an ML model requires an efficient method for the generation of training data. Here we apply an accurate and efficient protocol to collect training data for constructing a neural network based ML interatomic potential for nanosilicate clusters. Initial training data are taken from normal modes and farthest point sampling. Later on, the set of training data is extended via an active learning strategy in which new data are identified by the disagreement between an ensemble of ML models. The whole process is further accelerated by parallel sampling over structures. We use the ML model to run molecular dynamics (MD) simulations of nanosilicate clusters with various sizes, from which infrared spectra with anharmonicity included can be extracted. Such spectroscopic data are needed for understanding the properties of silicate dust grains in the interstellar medium (ISM) and in circumstellar environments.
△ Less
Submitted 25 May, 2023;
originally announced May 2023.
-
Relativistic Atomic Structure of Au IV and the Os Isoelectronic Sequence: opacity data for kilonova ejecta
Authors:
Z. S. Taghadomi,
Y. Wan,
A. Flowers,
P. C. Stancil,
B. M. McLaughlin,
S. Bromley,
J. P. Marler,
C. E. Sosolik,
S. Loch
Abstract:
Direct detection of gravitational waves (GW) on Aug. 17, 2017, propagating from a binary neutron star merger, opened the era of multimessenger astronomy. The ejected material from neutron star mergers, or kilonova, is a good candidate for optical and near infrared follow-up observations after the detection of GWs. The kilonova from the ejecta of GW1780817 provided the first evidence for the astrop…
▽ More
Direct detection of gravitational waves (GW) on Aug. 17, 2017, propagating from a binary neutron star merger, opened the era of multimessenger astronomy. The ejected material from neutron star mergers, or kilonova, is a good candidate for optical and near infrared follow-up observations after the detection of GWs. The kilonova from the ejecta of GW1780817 provided the first evidence for the astrophysical site of the synthesis of heavy nuclei through the rapid neutron capture process or rprocess. Since properties of the emission are largely affected by opacities of the ejected material, enhancements in the available, but limited rprocess atomic data have been motivated recently. However, given the complexity of the electronic structure of these heavy elements, considerable efforts are still needed to converge to a reliable set of atomic structure data. The aim of this work is to alleviate this situation for low charge state elements in the Os-like isoelectronic sequence. In this regard, the general purpose relativistic atomic structure packages (GRASP0 and GRASP2K) were used to obtain energy levels and transition probabilities (E1 and M1). We provide line lists and expansion opacities for a range of r-process elements. We focus here on the Os isoelectronic sequence (Os I, Ir II, Pt III, Au IV, Hg V). The results are benchmarked against existing experimental data and prior calculations, and predictions of emission spectra relevant to kilonovae are provided. Finestructure (M1) lines in the infrared potentially observable by the James Webb Space Telescope are highlighted.
△ Less
Submitted 5 July, 2022;
originally announced July 2022.
-
Bottom-up dust nucleation theory in oxygen-rich evolved stars I. Aluminium oxide clusters
Authors:
David Gobrecht,
John M. C. Plane,
Stefan T. Bromley,
Leen Decin,
Sergio Cristallo,
Sanjay Sekaran
Abstract:
Aluminum oxide (alumina, Al$_{2}$O$_{3}$) is a promising candidate as a primary dust condensate in the atmospheres of oxygen-rich evolved stars. Therefore, alumina \textit{seed} particles might trigger the onset of stellar dust formation and of stellar mass loss in the wind. However, the formation of alumina dust grains is not well understood.} {To shed light on the initial steps of cosmic dust fo…
▽ More
Aluminum oxide (alumina, Al$_{2}$O$_{3}$) is a promising candidate as a primary dust condensate in the atmospheres of oxygen-rich evolved stars. Therefore, alumina \textit{seed} particles might trigger the onset of stellar dust formation and of stellar mass loss in the wind. However, the formation of alumina dust grains is not well understood.} {To shed light on the initial steps of cosmic dust formation (i.e. nucleation) in oxygen-rich environments by a quantum-chemical bottom-up approach.} {Starting with an elemental gas-phase composition, we construct a detailed chemical-kinetic network describing the formation and destruction of aluminium-bearing molecules and dust-forming (Al$_{2}$O$_{3}$)$_n$ clusters up to the size of dimers ($n$=2) coagulating to tetramers ($n=$4). Intermediary species include the prevalent gas-phase molecules AlO and AlOH, and Al$_x$O$_y$ clusters with $x=$1$-$5, $y=$1$-$6. The resulting extensive network is applied to two model stars, representing a semi-regular variable and a Mira-type, and to different circumstellar gas trajectories including a non-pulsating outflow and a pulsating model. The growth of larger-sized (Al$_{2}$O$_{3}$)$_n$ clusters with $n=$4$-$10 is described by the temperature-dependent Gibbs free energies of the most favourable structures (i.e. the global minima clusters) as derived from global optimisation techniques and calculated by density functional theory. We provide energies, bond characteristics, electrostatic properties and vibrational spectra of the clusters as a function of size $n$ and compare these to corundum corresponding to the crystalline bulk limit ($n\rightarrow \infty$). {The circumstellar aluminium gas-phase chemistry in oxygen-rich giants is primarily controlled by AlOH and AlO, which are tightly coupled by the reactions AlO+H$_2$, AlO+H$_2$O, and their reverse. Models of ...}
△ Less
Submitted 21 October, 2021;
originally announced October 2021.
-
Atomic iron and nickel in the coma of C/1996 B2 (Hyakutake): production rates, emission mechanisms, and possible parents
Authors:
Steven Bromley,
Brynna Neff,
Stuart Loch,
Joan Marler,
Juraj Országh,
Kumar Venkataramani,
Dennis Bodewits
Abstract:
Two papers recently reported the detection of gaseous nickel and iron in the comae of over 20 comets from observations collected over two decades, including interstellar comet 2I/Borisov. To evaluate the state of the laboratory data in support of these identifications, we re-analyzed archived spectra of comet C/1996 B2 (Hyakutake), one of the nearest and brightest comets of the last century, using…
▽ More
Two papers recently reported the detection of gaseous nickel and iron in the comae of over 20 comets from observations collected over two decades, including interstellar comet 2I/Borisov. To evaluate the state of the laboratory data in support of these identifications, we re-analyzed archived spectra of comet C/1996 B2 (Hyakutake), one of the nearest and brightest comets of the last century, using a combined experimental and computational approach. We developed a new, many-level fluorescence model that indicates that the fluorescence emission of Fe I and Ni I vary greatly with heliocentric velocity. Combining this model with laboratory spectra of an Fe-Ni plasma, we identified 22 lines of Fe I and 14 lines of Ni I in the spectrum of Hyakutake. Using Haser models, we estimate the nickel and iron production rates as Q(Ni) = 2.6 - 4.1 x 10^22 s^-1 and Q(Fe) = 0.4 - 2.8 x 10^23 s^-1. From derived column densities, the Ni/Fe abundance ratio log10[Ni/Fe] = -0.15 +/- 0.07 deviates significantly from solar abundance ratios, and it is consistent with the ratios observed in solar system comets. Possible production and emission mechanisms are analyzed in context of existing laboratory measurements. Based on the observed spatial distributions, excellent fluorescence model agreement, and Ni/Fe ratio, our findings support an origin consisting of a short-lived unknown parent followed by fluorescence emission. Our models suggest that the strong heliocentric velocity dependence of the fluorescence efficiencies can provide a meaningful test of the physical process responsible for the Fe I and Ni I emission.
△ Less
Submitted 12 October, 2021; v1 submitted 8 June, 2021;
originally announced June 2021.
-
FUV Observations of the Inner Coma of 46P/Wirtanen
Authors:
John W. Noonan,
Walter M. Harris,
Steven Bromley,
Davide Farnocchia,
Jian-Yang Li,
Kathleen E. Mandt,
Joel Wm. Parker,
Kumar Venkataramani,
Dennis Bodewits
Abstract:
Far ultraviolet observations of comets yield information about the energetic processes that dissociate the sublimated gases from their primitive surfaces. Understanding which emission processes are dominant, their effects on the observed cometary spectrum, and how to properly invert the spectrum back to composition of the presumably pristine surface ices of a comet nuclei are all critical componen…
▽ More
Far ultraviolet observations of comets yield information about the energetic processes that dissociate the sublimated gases from their primitive surfaces. Understanding which emission processes are dominant, their effects on the observed cometary spectrum, and how to properly invert the spectrum back to composition of the presumably pristine surface ices of a comet nuclei are all critical components for proper interpretation and analysis of comets. The close approach of comet 46P/Wirtanen in 2018-2019 provided a unique opportunity to study the inner most parts of a cometary coma with the Hubble Space Telescope Cosmic Origins Spectrograph, rarely accessible with remote observations, at length scales (100's of km) and wavelengths (900-1430 Angstroms) previously probed only by the European Space Agency's Rosetta spacecraft. Our observations show a complex picture for the inner coma; atomic production rates for H and O that show water is the dominant source of both, an abundance of atomic sulfur that is difficult to explain with the lifetimes of common sulfur parent molecules, and a density distribution that is poorly fit with both Haser and vectorial models.
△ Less
Submitted 8 December, 2020;
originally announced December 2020.
-
Nucleation of small silicon carbide dust clusters in AGB stars
Authors:
David Gobrecht,
Sergio Cristallo,
Luciano Piersanti,
Stefan T. Bromley
Abstract:
Silicon carbide (SiC) grains are a major dust component in carbon-rich AGB stars. The formation pathways of these grains are, however, not fully understood.\ We calculate ground states and energetically low-lying structures of (SiC)$_n$, $n=1,16$ clusters by means of simulated annealing (SA) and Monte Carlo simulations of seed structures and subsequent quantum-mechanical calculations on the densit…
▽ More
Silicon carbide (SiC) grains are a major dust component in carbon-rich AGB stars. The formation pathways of these grains are, however, not fully understood.\ We calculate ground states and energetically low-lying structures of (SiC)$_n$, $n=1,16$ clusters by means of simulated annealing (SA) and Monte Carlo simulations of seed structures and subsequent quantum-mechanical calculations on the density functional level of theory. We derive the infrared (IR) spectra of these clusters and compare the IR signatures to observational and laboratory data.\ According to energetic considerations, we evaluate the viability of SiC cluster growth at several densities and temperatures, characterising various locations and evolutionary states in circumstellar envelopes.\ We discover new, energetically low-lying structures for Si$_{4}$C$_{4}$, Si$_{5}$C$_{5}$, Si$_{15}$C$_{15}$ and Si$_{16}$C$_{16}$, and new ground states for Si$_{10}$C$_{10}$ and Si$_{15}$C$_{15}$. The clusters with carbon-segregated substructures tend to be more stable by 4-9 eV than their bulk-like isomers with alternating Si-C bonds. However, we find ground states with cage ("bucky"-like) geometries for Si$_{12}$C$_{12}$ and Si$_{16}$C$_{16}$ and low-lying, stable cage structures for n $\ge$ 12. The latter findings indicate thus a regime of clusters sizes that differs from small clusters as well as from large-scale crystals. Thus, and owing to their stability and geometry, the latter clusters may mark a transition from a quantum-confined cluster regime to crystalline, solid bulk-material.
The calculated vibrational IR spectra of the ground-state SiC clusters shows significant emission. They include the 10-13 $μ$m wavelength range and the 11.3 $μ$m feature inferred from laboratory measurements and observations, respectively, though the overall intensities are rather low.
△ Less
Submitted 19 April, 2017;
originally announced April 2017.
-
Dust formation in the oxygen-rich AGB star IK Tau
Authors:
D. Gobrecht,
I. Cherchneff,
A. Sarangi,
J. M. C. Plane,
S. T. Bromley
Abstract:
We model the synthesis of molecules and dust in the inner wind of the oxygen-rich Mira-type star IK Tau, by considering the effects of periodic shocks induced by the stellar pulsation on the gas, and by following the non-equilibrium chemistry in the shocked gas layers between 1 and 10 Rstar. We consider a complete set of molecules and dust clusters, and combine the nucleation phase of dust formati…
▽ More
We model the synthesis of molecules and dust in the inner wind of the oxygen-rich Mira-type star IK Tau, by considering the effects of periodic shocks induced by the stellar pulsation on the gas, and by following the non-equilibrium chemistry in the shocked gas layers between 1 and 10 Rstar. We consider a complete set of molecules and dust clusters, and combine the nucleation phase of dust formation with the condensation of these clusters into dust grains. Our derived molecular abundances and dust properties are compared to the most recent observational data. The chemistry is described by using a chemical kinetic network of reactions and the condensation mechanism is described by a Brownian formalism. The shocks drive an active non-equilibrium chemistry in the dust formation zone of IK Tau where the collision destruction of CO in the post-shock gas triggers the formation of C-bearing species such as HCN and CS. Most of the modelled molecular abundances agree well with the latest values derived from Herschel data. Clusters of alumina are produced within 2 Rstar and lead to a population of alumina grains close to the stellar surface. Clusters of silicates form at larger radii (r > 3 Rstar), where their nucleation is triggered by the formation of HSiO and H2SiO. They efficiently condense and reach their final grain size distribution between ~ 6 and 8 Rstar, with a major population of medium size grains peaking at~ 0.02 microns. This two dust-shell configuration agrees with recent interferometric observations. The derived dust-to-gas mass ratio for IK Tau is in the range 1-6x10^-3 and agrees with values derived from observations of O-rich Mira-type stars. Our results confirm the importance of periodic shocks in chemically shaping the inner wind of AGB stars and providing gas conditions conducive to the efficient synthesis of molecules and dust by non-equilibrium processes.
△ Less
Submitted 16 October, 2015; v1 submitted 25 September, 2015;
originally announced September 2015.
-
Dust in brown dwarfs and extra-solar planets IV. Assessing TiO2 and SiO nucleation for cloud formation modeling
Authors:
E. K. H. Lee,
Ch. Helling,
H. Giles,
S. T. Bromley
Abstract:
Clouds form in atmospheres of brown dwarfs and planets. The cloud particle formation processes are similar to the dust formation process studied in circumstellar shells of AGB stars and in Supernovae. Cloud formation modelling in substellar objects requires gravitational settling and element replenishment in addition to element depletion. All processes depend on the local conditions, and a simulta…
▽ More
Clouds form in atmospheres of brown dwarfs and planets. The cloud particle formation processes are similar to the dust formation process studied in circumstellar shells of AGB stars and in Supernovae. Cloud formation modelling in substellar objects requires gravitational settling and element replenishment in addition to element depletion. All processes depend on the local conditions, and a simultaneous treatment is required. We apply new material data in order to assess our cloud formation model results regarding the treatment of the formation of condensation seeds. We re-address the question of the primary nucleation species in view of new (TiO2)_N-cluster data and new SiO vapour pressure data. We apply the density functional theory using the computational chemistry package Gaussian 09 to derive updated thermodynamical data for (TiO2)_N-clusters as input for our TiO2 seed formation model. We test different nucleation treatments and their effect on the overall cloud structure by solving a system of dust moment equations and element conservation or a pre-scribed Drift-Phoenix atmosphere structure. Updated Gibbs free energies for the (TiO2)_N-clusters are presented, and a slightly temperature dependent surface tension for T=500 ... 2000K with an average value of sigma_infty = 480.6 erg 1/cm2. The TiO2-seed formation rate changes only slightly with the updated cluster data. A considerably larger effect on the rate of seed formation, and hence on grain size and dust number density, results from a switch to SiO-nucleation. Despite the higher abundance of SiO over TiO2 in the gas phase, TiO2 remains considerably more efficient in forming condensation seeds by homogeneous nucleation followed by heterogeneous grain growth. The paper discussed the effect on the cloud structure in more detail.
△ Less
Submitted 24 October, 2014;
originally announced October 2014.