We study the influence of the turbulent transport on ice chemistry in protoplanetary disks, focusing on carbon- and nitrogen-bearing molecules. Chemical rate equations are solved with the diffusion ...term, mimicking the turbulent mixing in the vertical direction. Turbulence can bring ice-coated dust grains from the midplane to the warm irradiated disk surface, and the ice mantles are reprocessed by photoreactions, thermal desorption, and surface reactions. The upward transport decreases the abundance of methanol and ammonia ices at r lap 30 AU because warm dust temperature prohibits their reformation on grain surfaces. This reprocessing could explain the smaller abundances of carbon and nitrogen bearing molecules in cometary coma than those in low-mass protostellar envelopes. We also show the effect of mixing on the synthesis of complex organic molecules (COMs) in two ways: (1) transport of ices from the midplane to the disk surface and (2) transport of atomic hydrogen from the surface to the midplane. The former enhances the COMs formation in the disk surface, while the latter suppresses it in the midplane. Then, when mixing is strong, COMs are predominantly formed in the disk surface, while their parent molecules are (re)formed in the midplane. This cycle expands the COMs distribution both vertically and radially outward compared with that in the non-turbulent model. We derive the timescale of the sink mechanism by which CO and N sub(2) are converted to less volatile molecules to be depleted from the gas phase and find that the vertical mixing suppresses this mechanism in the inner disks.
We investigate nitrogen-isotope fractionation in forming and evolving molecular clouds using gas-ice astrochemical simulations. We find that the bulk gas can become depleted in heavy nitrogen (15N) ...due to the formation of 15N-enriched ices. Around the chemical transition from atomic nitrogen to N2, N15N is selectively photodissociated, which results in the enrichment of 15N in atomic nitrogen. As 15N-enriched atomic nitrogen is converted to ammonia ice via grain surface reactions, the bulk gas is depleted in 15N. The level of 15N depletion in the bulk gas can be up to a factor of two compared to the elemental nitrogen-isotope ratio, depending on the photodesorption yield of ammonia ice. Once the nitrogen isotopes are differentially partitioned between gas and solids in a molecular cloud, the condition should remain in the later stages of star formation (e.g., prestellar core) as long as the sublimation of ammonia ice is inefficient. Our model suggests that all of the N-bearing molecules in the cold gas of star-forming regions can be depleted in 15N, which is at least qualitatively consistent with the observations toward prestellar core L1544. In our models, icy species show both 15N and deuterium fractionation. The fractionation pattern within ice mantles is different between 15N and deuterium, reflecting their fractionation mechanisms; while the concentration of deuterium almost monotonically increases from the lower layers of the ice mantles to the upper layers, the concentration of 15N reaches the maximum at a certain depth and declines toward the surface.
Hot corino chemistry and warm carbon chain chemistry (WCCC) are driven by gas-grain interactions in star-forming cores: radical-radical recombination reactions to form complex organic molecules ...(COMs) in the ice mantle, sublimation of CH4 and COMs, and their subsequent gas-phase reactions. These chemical features are expected to depend on the composition of the ice mantle, which is set in the prestellar phase. We calculated the gas-grain chemical reaction network considering a layered ice mantle structure in star-forming cores to investigate how the hot corino chemistry and WCCC depend on the physical condition of the static phase before the onset of gravitational collapse. We found that WCCC becomes more active if the temperature is lower, or the visual extinction is lower in the static phase, or the static phase is longer. The dependence of hot corino chemistry on the static-phase condition is more complex. While CH3OH is less abundant in the models with a warmer static phase, some COMs are formed efficiently in those warm models because there are various formation paths of COMs. If the visual extinction is lower, photolysis makes COMs less abundant in the static phase. Once the collapse starts and visual extinction increases, however, COMs can be formed efficiently. The duration of the static phase does not largely affect COM abundances. The chemical diversity between prototypical hot corinos and hybrid sources, in which both COMs and carbon chains are reasonably abundant, can be explained by the variation of prestellar conditions. Deficiency of gaseous COMs in prototypical WCCC sources is, however, hard to reproduce within our models.
Abstract Understanding in which chemical forms phosphorus exists in star- and planet-forming regions and how phosphorus is delivered to planets are of great interest from the viewpoint of the origin ...of life on Earth. Phosphine (PH 3 ) is thought to be a key species to understanding phosphorus chemistry, but never has been detected in star- and planet-forming regions. We performed sensitive observations of the ortho-PH 3 1 0 − 0 0 transition (266.944 GHz) toward the low-mass prestellar core L1544 with the Atacama Compact Array stand-alone mode of the Atacama Large Millimeter/submillimeter Array. The line was not detected down to 3 σ levels in 0.07 km s −1 channels of 18 mK. The nondetection provides the upper limit to the gas-phase PH 3 abundance of 5 × 10 −12 with respect to H 2 in the central part of the core. Based on the gas-ice astrochemical modeling, we find the scaling relationship between the gas-phase PH 3 abundance and the volatile (gas and ice with larger volatility than water) P elemental abundance for given physical conditions. This characteristic and well-constrained physical properties of L1544 allow us to constrain the upper limit to the volatile P elemental abundance of 5 × 10 −9 , which is a factor of 60 lower than the overall P abundance in the interstellar medium. Then the majority of P should exist in refractory forms. The volatile P elemental abundance of L1544 is smaller than that in the coma of comet 67P/C-G, implying that the conversion of refractory phosphorus to volatile phosphorus could have occurred along the trail from the presolar core to the protosolar disk through, e.g., sputtering by accretion/outflow shocks.
Abstract
Interstellar chemistry in low-metallicity environments is crucial to understand chemical processes in the past metal-poor universe. Recent studies of interstellar molecules in nearby ...low-metallicity galaxies have suggested that metallicity has a significant effect on the chemistry of star-forming cores. Here we report the first detection of a hot molecular core in the extreme outer Galaxy, which is an excellent laboratory to study star formation and the interstellar medium in a Galactic low-metallicity environment. The target star-forming region, WB 89–789, is located at a galactocentric distance of 19 kpc. Our Atacama Large Millimeter/submillimeter Array observations in 241–246, 256–261, 337–341, and 349–353 GHz have detected a variety of carbon-, oxygen-, nitrogen-, sulfur-, and silicon-bearing species, including complex organic molecules (COMs) containing up to nine atoms, toward a warm (>100 K) and compact (<0.03 pc) region associated with a protostar (∼8 × 10
3
L
☉
). Deuterated species such as HDO, HDCO, D
2
CO, and CH
2
DOH are also detected. A comparison of fractional abundances of COMs relative to CH
3
OH between the outer Galactic hot core and an inner Galactic counterpart shows a remarkable similarity. On the other hand, the molecular abundances in the present source do not resemble those of low-metallicity hot cores in the Large Magellanic Cloud. The results suggest that great molecular complexity exists even in the primordial environment of the extreme outer Galaxy. The detection of another embedded protostar associated with high-velocity SiO outflows is also reported.
We investigate deuterium chemistry coupled with the nuclear spin-state chemistry of H2 and in protoplanetary disks. Multiple paths of deuterium fractionation are found; exchange reactions with D ...atoms, such as HCO+ + D, are effective in addition to those with HD. In a disk model with grain sizes appropriate for dark clouds, the freeze-out of molecules is severe in the outer midplane, while the disk surface is shielded from UV radiation. Gaseous molecules, including DCO+, thus become abundant at the disk surface, which tends to make their column density distribution relatively flat. If the dust grains have grown to millimeter size, the freeze-out rate of neutral species is reduced and the abundances of gaseous molecules, including DCO+ and N2D+, are enhanced in the cold midplane. Turbulent diffusion transports D atoms and radicals at the disk surface to the midplane, and stable ice species in the midplane to the disk surface. The effects of turbulence on chemistry are thus multifold; while DCO+ and N2D+ abundances increase or decrease depending on the regions, HCN and DCN in the gas and ice are greatly reduced at the innermost radii, compared to the model without turbulence. When cosmic rays penetrate the disk, the ortho-to-para ratio (OPR) of H2 is found to be thermal in the disk, except in the cold ( 10 K) midplane. We also analyze the OPR of and H2D+, as well as the main reactions of H2D+, DCO+, and N2D+, in order to analytically derive their abundances in the cold midplane.
Abstract Carbon isotope fractionation of CO has been reported in the disk around TW Hya, where elemental carbon is more abundant than elemental oxygen (C/O elem > 1). We investigated the effects of ...the C/O elem ratio on carbon fractionation using astrochemical models that incorporate isotope-selective photodissociation and isotope exchange reactions. The 12 CO/ 13 CO ratio could be lower than the elemental carbon isotope ratio due to isotope exchange reactions when the C/O elem ratio exceeds unity. The observed 12 CO/ 13 CO and H 12 CN/H 13 CN ratios around TW Hya could be reproduced when the C/O elem ratio is 2–5. In the vicinity of the lower boundary of the warm molecular layer, the formation of ices leads to the gas-phase C/O elem ratio approaching unity, irrespective of the total (gas + ice) C/O elem ratio. This phenomenon reduces the variation in the 12 CO/ 13 CO ratio across different C/O elem ratios.
Abstract
Observations have revealed that the elemental abundances of carbon and oxygen in the warm molecular layers of some protoplanetary disks are depleted compared to those in the interstellar ...medium by a factor of ∼10–100. Meanwhile, little is known about nitrogen. To investigate the time evolution of nitrogen, carbon, and oxygen elemental abundances in disks, we develop a one-dimensional plane-parallel model that incorporates dust settling, turbulent diffusion of dust and ices, as well as gas-ice chemistry including the chemistry driven by stellar UV/X-rays and galactic cosmic rays. We find that gaseous CO in the warm molecular layer is converted to CO
2
ice and locked up near the midplane via the combination of turbulent mixing (i.e., the vertical cold finger effect) and ice chemistry driven by stellar UV photons. On the other hand, gaseous N
2
, the main nitrogen reservoir in the warm molecular layer, is less processed by ice chemistry and exists as it is. Then, nitrogen depletion occurs solely through the vertical cold finger effect of N
2
. As the binding energy of N
2
is lower than that of CO and CO
2
, the degree of nitrogen depletion is smaller than that of carbon and oxygen depletion, leading to higher elemental abundance of nitrogen than that of carbon and oxygen. This evolution occurs within 1 Myr and proceeds further, when the
α
parameter for the diffusion coefficient is ≳10
−3
. Consequently, the N
2
H
+
/CO column density ratio increases with time. How the vertical transport affects the midplane ice composition is briefly discussed.
ABSTRACT We investigate the chemistry of ion molecules in protoplanetary disks, motivated by the detection of the N2H+ ring around TW Hya. While the ring inner radius coincides with the CO snow line, ...it is not apparent why N2H+ is abundant outside the CO snow line in spite of the similar sublimation temperatures of CO and N2. Using the full gas-grain network model, we reproduced the N2H+ ring in a disk model with millimeter grains. The chemical conversion of CO and N2 to less volatile species (sink effect hereinafter) is found to affect the N2H+ distribution. Since the efficiency of the sink depends on various parameters such as activation barriers of grain-surface reactions, which are not well constrained, we also constructed the no-sink model; the total (gas and ice) CO and N2 abundances are set constant, and their gaseous abundances are given by the balance between adsorption and desorption. Abundances of molecular ions in the no-sink model are calculated by analytical formulae, which are derived by analyzing the full-network model. The N2H+ ring is reproduced by the no-sink model, as well. The 2D (R-Z) distribution of N2H+, however, is different among the full-network model and no-sink model. The column density of N2H+ in the no-sink model depends sensitively on the desorption rate of CO and N2 and the cosmic-ray flux. We also found that N2H+ abundance can peak at the temperature slightly below the CO sublimation, even if the desorption energies of CO and N2 are the same.
ABSTRACT We investigate the chemistry in a radiation-hydrodynamics model of a star-forming core that evolves from a cold (∼10 K) prestellar core to the main accretion phase in ∼105 years. A ...rotationally supported gravitationally unstable disk is formed around a protostar. We extract the temporal variation of physical parameters in ∼1.5 × 103 SPH particles that end up in the disk, and perform post-processing calculations of the gas-grain chemistry adopting a three-phase model. Inside the disk, the SPH particles migrate both inward and outward. Since a significant fraction of volatiles such as CO can be trapped in the water-dominant ice in the three-phase model, the ice mantle composition depends not only on the current position in the disk, but also on whether the dust grain has ever experienced higher temperatures than the water sublimation temperature. Stable molecules such as H2O, CH4, NH3, and CH3OH are already abundant at the onset of gravitational collapse and are simply sublimated as the fluid parcels migrate inside the water snow line. On the other hand, various molecules such as carbon chains and complex organic molecules (COMs) are formed in the disk. The COMs abundance sensitively depends on the outcomes of photodissociation and diffusion rates of photofragments in bulk ice mantle. As for S-bearing species, H2S ice is abundant in the collapse phase. In the warm regions in the disk, H2S is sublimated to be destroyed, while SO, H2CS, OCS, and SO2 become abundant.