No Arabic abstract
(Abridged) The terrestrial planets, comets, and meteorites are significantly enriched in 15N compared to the Sun and Jupiter. While the solar and jovian nitrogen isotope ratio is believed to represent the composition of the protosolar nebula, a still unidentified process has caused 15N-enrichment in the solids. Several mechanisms have been proposed to explain the variations, including chemical fractionation. However, observational results that constrain the fractionation models are scarce. While there is evidence of 15N-enrichment in prestellar cores, it is unclear how the signature evolves into the protostellar phases. Our aim is to measure the 14N/15N ratio around three nearby, embedded low-to-intermediate-mass protostars. Isotopologues of HCN and HNC were used to probe the 14N/15N ratio. A selection of H13CN, HC15N, HN13C, and H15NC transitions was observed with the APEX telescope. The 14N/15N ratios were derived from the integrated intensities assuming a standard 12C/13C ratio. The assumption of optically thin emission was verified using radiative transfer modeling and hyperfine structure fitting. Two sources, IRAS 16293A and R CrA IRS7B, show 15N-enrichment by a factor of around 1.5-2.5 in both HCN and HNC with respect to the solar composition. Solar composition cannot be excluded for the third source, OMC-3 MMS6. Furthermore, there are indications of a trend toward increasing 14N/15N ratios with increasing outer envelope temperature. The enhanced 15N abundances in HCN and HNC found in two Class~0 sources (14N/15N of 160-290) and the tentative trend toward a temperature-dependent 14N/15N ratio are consistent with the chemical fractionation scenario, but 14N/15N ratios from additional tracers are indispensable for testing the models. Spatially resolved observations are needed to distinguish between chemical fractionation and isotope-selective photochemistry.
Aims: The two stable isotopes of nitrogen, 14N and 15N, exhibit a range of abundance ratios both inside and outside the solar system. The elemental ratio in the solar neighborhood is 440. Recent ALMA observations showed HCN/HC15N ratios from 83 to 156 in six T Tauri and Herbig disks and a CN/C15 N ratio of 323 +/- 30 in one T Tauri star. We aim to determine the dominant mechanism responsible for these enhancements of 15N: low-temperature exchange reactions or isotope-selective photodissociation of N2. Methods: Using the thermochemical code DALI, we model the nitrogen isotope chemistry in circumstellar disks with a 2D axisymmetric geometry. Our chemical network is the first to include both fractionation mechanisms for nitrogen. The model produces abundance profiles and isotope ratios for several key N-bearing species. We study how these isotope ratios depend on various disk parameters. Results: The formation of CN and HCN is closely coupled to the vibrational excitation of H2 in the UV-irradiated surface layers of the disk. Isotope fractionation is completely dominated by isotope-selective photodissociation of N2. The column density ratio of HCN over HC15N in the disks inner 100 au does not depend strongly on the disk mass, the flaring angle or the stellar spectrum, but it is sensitive to the grain size distribution. For larger grains, self-shielding of N2 becomes more important relative to dust extinction, leading to stronger isotope fractionation. Between disk radii of ~50 and 200 au, the models predict HCN/HC15N and CN/C15N abundance ratios consistent with observations of disks and comets. The HCN/HC15N and CN/C15N column density ratios in the models are a factor of 2-3 higher than those inferred from the ALMA observations.
We perform a comparative numerical hydrodynamics study of embedded protostellar disks formed as a result of the gravitational collapse of cloud cores of distinct mass (M_cl=0.2--1.7 M_sun) and ratio of rotational to gravitational energy (beta=0.0028--0.023). An increase in M_cl and/or beta leads to the formation of protostellar disks that are more susceptible to gravitational instability. Disk fragmentation occurs in most models but its effect is often limited to the very early stage, with the fragments being either dispersed or driven onto the forming star during tens of orbital periods. Only cloud cores with high enough M_cl or beta may eventually form wide-separation binary/multiple systems with low mass ratios and brown dwarf or sub-solar mass companions. It is feasible that such systems may eventually break up, giving birth to rogue brown dwarfs. Protostellar disks of {it equal} age formed from cloud cores of greater mass (but equal beta) are generally denser, hotter, larger, and more massive. On the other hand, protostellar disks formed from cloud cores of higher beta (but equal M_cl) are generally thinner and colder but larger and more massive. In all models, the difference between the irradiation temperature and midplane temperature triangle T is small, except for the innermost regions of young disks, dense fragments, and disks outer edge where triangle T is negative and may reach a factor of two or even more. Gravitationally unstable, embedded disks show radial pulsations, the amplitude of which increases along the line of increasing M_cl and beta but tends to diminish as the envelope clears. We find that single stars with a disk-to-star mass ratio of order unity can be formed only from high-beta cloud cores, but such massive disks are unstable and quickly fragment into binary/multiple systems.
We study the relationship between the mid-infrared and sub-mm variability of deeply embedded protostars using the multi-epoch data from the Wide Infrared Survey Explorer ($WISE$/NEOWISE) and the ongoing James Clerk Maxwell Telescope (JCMT) transient survey. Our search for signs of stochastic (random) and/or secular (roughly monotonic in time) variability in a sample of 59 young stellar objects (YSOs) revealed that 35 are variable in at least one of the two surveys. This variability is dominated by secular changes. Of those objects with secular variability, 14 objects ($22%$ of the sample) show correlated secular variability over mid-IR and sub-mm wavelengths. Variable accretion is the likely mechanism responsible for this type of variability. Fluxes of YSOs that vary in both wavelengths follow a relation of $log_{10} F_{4.6}(t)=eta log_{10} F_{850}(t)$ between the mid-IR and sub-mm, with $eta=5.53pm0.29$. This relationship arises from the fact that sub-mm fluxes respond to the dust temperature in the larger envelope whereas the mid-IR emissivity is more directly proportional to the accretion luminosity. The exact scaling relation, however, depends on the structure of the envelope, the importance of viscous heating in the disc, and dust opacity laws.
The physical evolution of Young Stellar Objects (YSOs) is accompanied by an enrichment of the molecular complexity, mainly triggered by the heating and energetic processing of the astrophysical ices. In this paper, a study of how the ice column density varies across the protostellar evolution has been performed. Tabulated data of H$_2$O, CO$_2$, CH$_3$OH, HCOOH observed by ground- and space-based telescopes toward 27 early-stage YSOs were taken from the literature. The observational data shows that ice column density and spectral index ($alpha$), used to classify the evolutionary stage, are well correlated. A 2D continuum radiative transfer simulation containing bare and grains covered by ices at different levels of cosmic-ray processing were used to calculate the Spectral Energy Distributions (SEDs) in different angle inclinations between face-on and edge-on configuration. The H$_2$O:CO$_2$ ice mixture was used to address the H$_2$O and CO$_2$ column density variation whereas the CH$_3$OH and HCOOH are a byproduct of the virgin ice after the energetic processing. The simulated spectra were used to calculate the ice column densities of YSOs in an evolutionary sequence. As a result, the models show that the ice column density variation of HCOOH with $alpha$ can be justified by the envelope dissipation and ice energetic processing. On the other hand, the ice column densities are mostly overestimated in the cases of H$_2$O, CO$_2$ and CH$_3$OH, even though the physical and cosmic-ray processing effects are taken into account.
Protostellar flares are rapid magnetic energy release events associated with formation of hot plasma in protostars. In the previous models of protostellar flares, the interaction between a protostellar magnetosphere with the surrounding disk plays crucial roles in building-up and releasing the magnetic energy. However, it remains unclear if protostars indeed have magnetospheres because vigorous disk accretion and strong disk magnetic fields in the protostellar phase may destroy the magnetosphere. Considering this possibility, we investigate the energy accumulation and release processes in the absence of a magnetosphere using a three-dimensional magnetohydrodynamic simulation. Our simulation reveals that protostellar flares are repeatedly produced even in such a case. Unlike in the magnetospheric models, the protostar accumulates magnetic energy by acquiring large-scale magnetic fields from the disk by accretion. Protostellar flares occur when a portion of the large-scale magnetic fields are removed from the protostar as a result of magnetic reconnection. Protostellar flares in the simulation are consistent with observations; the released magnetic energy (up to $sim 3times 10^{38}$ erg) is large enough to drive observed flares, and the flares produce hot ejecta. The expelled magnetic fields enhance accretion, and the energy build-up and release processes are repeated as a result. The magnetic flux removal via reconnection leads to redistribution of magnetic fields in the inner disk. We therefore consider that protostellar flares will play an important role in the evolution of the disk magnetic fields in the vicinity of protostars.