No Arabic abstract
Snowlines of major volatiles regulate the gas and solid C/N/O ratios in the planet-forming midplanes of protoplanetary disks. Snow surfaces are the 2D extensions of snowlines in the outer disk regions, where radiative heating results in a decreasing temperature with disk height. CO and N$_2$ are two of the most abundant carriers of C, N and O. N$_2$H$^+$ can be used to probe the snow surfaces of both molecules, because it is destroyed by CO and formed from N$_2$. Here we present Atacama Large Millimeter/submillimeter Array (ALMA) observations of N$_2$H$^+$ at 0.2$$-0.4$$ resolution in the disks around LkCa 15, GM Aur, DM Tau, V4046 Sgr, AS 209, and IM Lup. We find two distinctive emission morphologies: N$_2$H$^+$ is either present in a bright, narrow ring surrounded by extended tenuous emission, or in a broad ring. These emission patterns can be explained by two different kinds of vertical temperature structures. Bright, narrow N$_2$H$^+$ rings are expected in disks with a thick Vertically Isothermal Region above the Midplane (VIRaM) layer (LkCa 15, GM Aur, DM Tau) where the N$_2$H$^+$ emission peaks between the CO and N$_2$ snowlines. Broad N$_2$H$^+$ rings come from disks with a thin VIRaM layer (V4046 Sgr, AS 209, IM Lup). We use a simple model to extract the first sets of CO and N$_2$ snowline pairs and corresponding freeze-out temperatures towards the disks with a thick VIRaM layer. The results reveal a range of N$_2$ and CO snowline radii towards stars of similar spectral type, demonstrating the need for empirically determined snowlines in disks.
{CircumStellar Envelopes (CSEs) of stars are complex chemical objects for which theoretical models encounter difficulties in elaborating a comprehensive overview of the occurring chemical processes. Along with photodissociation, ion-neutral reactions and dissociative recombination might play an important role in controlling molecular growth in outer CSEs. The aim of this work is to provide experimental insights into pathways of photochemistry-driven molecular growth within outer CSEs to draw a more complete picture of the chemical processes occurring within these molecule-rich environments. A simplified CSE environment was therefore reproduced in the laboratory through gas-phase experiments exposing relevant gas mixtures to an Extreme UltraViolet (EUV) photon source. This photochemical reactor should ultimately allow us to investigate chemical processes and their resulting products occurring under conditions akin to outer CSEs. We used a recently developed EUV lamp coupled to the APSIS photochemical cell to irradiate CSE relevant gas mixtures of H$_2$, CO and N$_2$, at one wavelength, 73.6 nm. The detection and identification of chemical species in the photochemical reactor was achieved through in-situ mass spectrometry analysis of neutral and cationic molecules. We find that exposing CO-N$_2$-H$_2$ gas mixtures to EUV photons at 73.6 nm induces photochemical reactions that yield the formation of complex, neutral and ionic species. Our work shows that N$_2$H$^+$ can be formed through photochemistry along with highly oxygenated ion molecules like HCO$^+$ in CSE environments. We also observe neutral N-rich organic species including triazole and aromatic molecules. These results confirm the suitability of our experimental setting to investigate photochemical reactions and provide fundamental insights into the mechanisms of molecular growth in the outer CSEs.
The earliest atmospheres of rocky planets originate from extensive volatile release during magma ocean epochs that occur during assembly of the planet. These establish the initial distribution of the major volatile elements between different chemical reservoirs that subsequently evolve via geological cycles. Current theoretical techniques are limited in exploring the anticipated range of compositional and thermal scenarios of early planetary evolution, even though these are of prime importance to aid astronomical inferences on the environmental context and geological history of extrasolar planets. Here, we present a coupled numerical framework that links an evolutionary, vertically-resolved model of the planetary silicate mantle with a radiative-convective model of the atmosphere. Using this method we investigate the early evolution of idealized Earth-sized rocky planets with end-member, clear-sky atmospheres dominated by either H$_2$, H$_2$O, CO$_2$, CH$_4$, CO, O$_2$, or N$_2$. We find central metrics of early planetary evolution, such as energy gradient, sequence of mantle solidification, surface pressure, or vertical stratification of the atmosphere, to be intimately controlled by the dominant volatile and outgassing history of the planet. Thermal sequences fall into three general classes with increasing cooling timescale: CO, N$_2$, and O$_2$ with minimal effect, H$_2$O, CO$_2$, and CH$_4$ with intermediate influence, and H$_2$ with several orders of magnitude increase in solidification time and atmosphere vertical stratification. Our numerical experiments exemplify the capabilities of the presented modeling framework and link the interior and atmospheric evolution of rocky exoplanets with multi-wavelength astronomical observations.
Context. Studying gas chemistry in protoplanetary disks is key to understanding the process of planet formation. Sulfur chemistry in particular is poorly understood in interstellar environments, and the location of the main reservoirs remains unknown. Protoplanetary disks in Taurus are ideal targets for studying the evolution of the composition of planet forming systems. Aims. We aim to elucidate the chemical origin of sulfur-bearing molecular emission in protoplanetary disks, with a special focus on H$_2$S emission, and to identify candidate species that could become the main molecular sulfur reservoirs in protoplanetary systems. Methods. We used IRAM 30m observations of nine gas-rich young stellar objects (YSOs) in Taurus to perform a survey of sulfur-bearing and oxygen-bearing molecular species. In this paper we present our results for the CS 3-2 ($ u_0$ = 146.969 GHz), H$_2$CO 2$_{11}$-1$_{10}$ ($ u_0$ = 150.498 GHz), and H$_2$S 1$_{10}$-1$_{01}$ ($ u_0$ = 168,763 GHz) emission lines. Results. We detected H$_2$S emission in four sources out of the nine observed, significantly increasing the number of detections toward YSOs. We also detected H$_2$CO and CS in six out of the nine. We identify a tentative correlation between H$_2$S 1$_{10}$-1$_{01}$ and H$_2$CO 2$_{11}$-1$_{10}$ as well as a tentative correlation between H$_2$S 1$_{10}$-1$_{01}$ and H$_2$O 8$_{18}$-7$_{07}$. By assuming local thermodynamical equilibrium, we computed column densities for the sources in the sample, with N(o-H$_2$S) values ranging between $2.6times10^{12}$ cm$^{-2}$ and $1.5times10^{13}$ cm$^{-2}$.
The modelling of emission spectra of molecules seen in interstellar clouds requires the knowledge of collisional rate coefficients. Among the commonly observed species, N$_2$H$^+$ is of particular interest since it was shown to be a good probe of the physical conditions of cold molecular clouds. Thus, we have calculated hyperfine-structure resolved excitation rate coefficients of N$_2$H$^+$(X$^1Sigma^+$) by H$_2(j=0)$, the most abundant collisional partner in the cold interstellar medium. The calculations are based on a new potential energy surface, obtained from highly correlated {it ab initio} calculations. State-to-state rate coefficients between the first hyperfine levels were calculated, for temperatures ranging from 5 K to 70 K. By comparison with previously published N$_2$H$^+$-He rate coefficients, we found significant differences which cannot be reproduced by a simple scaling relationship. As a first application, we also performed radiative transfer calculations, for physical conditions typical of cold molecular clouds. We found that the simulated line intensities significantly increase when using the new H$_2$ rate coefficients, by comparison with the predictions based on the He rate coefficients. In particular, we revisited the modelling of the N$_2$H$^+$ emission in the LDN 183 core, using the new collisional data, and found that all three of the density, gas kinetic temperature and N$_2$H$^+$ abundance had to be revised.
ALMA observations of protoplanetary disks confirm earlier indications that there is a clear difference between the dust and gas radial extents. The origin of this difference is still debated, with both radial drift of the dust and optical depth effects suggested in the literature. In this work, the feedback of realistic dust particle distributions onto the gas chemistry and molecular emissivity is investigated, with a particular focus on CO isotopologues. The radial dust grain size distribution is determined using dust evolution models that include growth, fragmentation and radial drift. A new version of the code DALI is used to take into account how dust surface area and density influence the disk thermal structure, molecular abundances and excitation. The difference of dust and gas radial sizes is largely due to differences in the optical depth of CO lines and millimeter continuum, without the need to invoke radial drift. The effect of radial drift is primarily visible in the sharp outer edge of the continuum intensity profile. The gas outer radius probed by $^{12}$CO emission can easily differ by a factor of $sim 2$ between the models for a turbulent $alpha$ ranging between typical values. Grain growth and settling concur in thermally decoupling the gas and dust components, due to the low collision rate with large grains. As a result, the gas can be much colder than the dust at intermediate heights, reducing the CO excitation and emission, especially for low turbulence values. Also, due to disk mid-plane shadowing, a second CO thermal desorption (rather than photodesorption) front can occur in the warmer outer mid-plane disk. The models are compared to ALMA observations of HD 163296 as a test case. In order to reproduce the observed CO snowline of the system, a binding energy for CO typical of ice mixtures needs to be used rather than the lower pure CO value.