Abundant hydrocarbons in the disk around a very-low-mass star.

Very-low-mass stars (those less than 0.3 solar masses) host orbiting terrestrial planets more frequently than other types of stars. The compositions of those planets are largely unknown but are expected to relate to the protoplanetary disk in which they form. We used James Webb Space Telescope mid-infrared spectroscopy to investigate the chemical composition of the planet-forming disk around ISO-ChaI 147, a 0.11-solar-mass star. The inner disk has a carbon-rich chemistry; we identified emission from 13 carbon-bearing molecules, including ethane and benzene. The high column densities of hydrocarbons indicate that the observations probe deep into the disk. The high carbon-to-oxygen ratio indicates radial transport of material within the disk, which we predict would affect the bulk composition of any planets forming in the disk.

Emission lines due to ionizing radiation from a compact object in the remnant of Supernova 1987A.

The nearby Supernova 1987A was accompanied by a burst of neutrino emission, which indicates that a compact object (a neutron star or black hole) was formed in the explosion. There has been no direct observation of this compact object. In this work, we observe the supernova remnant with JWST spectroscopy, finding narrow infrared emission lines of argon and sulfur. The line emission is spatially unresolved and blueshifted in velocity relative to the supernova rest frame. We interpret the lines as gas illuminated by a source of ionizing photons located close to the center of the expanding ejecta. Photoionization models show that the line ratios are consistent with ionization by a cooling neutron star or a pulsar wind nebula. The velocity shift could be evidence for a neutron star natal kick.

Outflows from the youngest stars are mostly molecular.

The formation of stars and planets is accompanied not only by the build-up of matter, namely accretion, but also by its expulsion in the form of highly supersonic jets that can stretch for several parsecs1,2. As accretion and jet activity are correlated and because young stars acquire most of their mass rapidly early on, the most powerful jets are associated with the youngest protostars3. This period, however, coincides with the time when the protostar and its surroundings are hidden behind many magnitudes of visual extinction. Millimetre interferometers can probe this stage but only for the coolest components3. No information is provided on the hottest (greater than 1,000 K) constituents of the jet, that is, the atomic, ionized and high-temperature molecular gases that are thought to make up the jet's backbone. Detecting such a spine relies on observing in the infrared that can penetrate through the shroud of dust. Here we report near-infrared observations of Herbig-Haro 211 from the James Webb Space Telescope, an outflow from an analogue of our Sun when it was, at most, a few times 104 years old. These observations reveal copious emission from hot molecules, explaining the origin of the 'green fuzzies'4-7 discovered nearly two decades ago by the Spitzer Space Telescope8. This outflow is found to be propagating slowly in comparison to its more evolved counterparts and, surprisingly, almost no trace of atomic or ionized emission is seen, suggesting its spine is almost purely molecular.

Water in the terrestrial planet-forming zone of the PDS 70 disk.

Terrestrial and sub-Neptune planets are expected to form in the inner (less than 10 AU) regions of protoplanetary disks1. Water plays a key role in their formation2-4, although it is yet unclear whether water molecules are formed in situ or transported from the outer disk5,6. So far Spitzer Space Telescope observations have only provided water luminosity upper limits for dust-depleted inner disks7, similar to PDS 70, the first system with direct confirmation of protoplanet presence8,9. Here we report JWST observations of PDS 70, a benchmark target to search for water in a disk hosting a large (approximately 54 AU) planet-carved gap separating an inner and outer disk10,11. Our findings show water in the inner disk of PDS 70. This implies that potential terrestrial planets forming therein have access to a water reservoir. The column densities of water vapour suggest in-situ formation via a reaction sequence involving O, H2 and/or OH, and survival through water self-shielding5. This is also supported by the presence of CO2 emission, another molecule sensitive to ultraviolet photodissociation. Dust shielding, and replenishment of both gas and small dust from the outer disk, may also play a role in sustaining the water reservoir12. Our observations also reveal a strong variability of the mid-infrared spectral energy distribution, pointing to a change of inner disk geometry.

A cryogenic ice setup to simulate carbon atom reactions in interstellar ices.

The design, implementation, and performance of a customized carbon atom beam source for the purpose of investigating solid-state reaction routes in interstellar ices in molecular clouds are discussed. The source is integrated into an existing ultrahigh vacuum setup, SURFace REaction SImulation DEvice (SURFRESIDE2), which extends this double atom (H/D, O, and N) beamline apparatus with a third atom (C) beamline to a unique system that is fully suited to explore complex organic molecule solid-state formation under representative interstellar cloud conditions. The parameter space for this system is discussed, which includes the flux of the carbon atoms hitting the ice sample, their temperature, and the potential impact of temperature on ice reactions. Much effort has been put into constraining the beam size to within the limits of the sample size with the aim of reducing carbon pollution inside the setup. How the C-atom beam performs is quantitatively studied through the example experiment, C + 18O2, and supported by computationally derived activation barriers. The potential for this source to study the solid-state formation of interstellar complex organic molecules through C-atom reactions is discussed.

Abundant molecular oxygen in the coma of comet 67P/Churyumov-Gerasimenko.

The composition of the neutral gas comas of most comets is dominated by H2O, CO and CO2, typically comprising as much as 95 per cent of the total gas density. In addition, cometary comas have been found to contain a rich array of other molecules, including sulfuric compounds and complex hydrocarbons. Molecular oxygen (O2), however, despite its detection on other icy bodies such as the moons of Jupiter and Saturn, has remained undetected in cometary comas. Here we report in situ measurement of O2 in the coma of comet 67P/Churyumov-Gerasimenko, with local abundances ranging from one per cent to ten per cent relative to H2O and with a mean value of 3.80 ± 0.85 per cent. Our observations indicate that the O2/H2O ratio is isotropic in the coma and does not change systematically with heliocentric distance. This suggests that primordial O2 was incorporated into the nucleus during the comet's formation, which is unexpected given the low upper limits from remote sensing observations. Current Solar System formation models do not predict conditions that would allow this to occur.

Molecular dynamics simulations of CO2 formation in interstellar ices.

CO2 ice is one of the most abundant components in ice-coated interstellar ices besides H2O and CO, but the most favorable path to CO2 ice is still unclear. Molecular dynamics calculations on the ultraviolet photodissociation of different kinds of CO-H2O ice systems have been performed at 10 K in order to demonstrate that the reaction between CO and an OH molecule resulting from H2O photodissociation through the first excited state is a possible route to form CO2 ice. However, our calculations, which take into account different ice surface models, suggest that there is another product with a higher formation probability ((3.00 ± 0.07) × 10(-2)), which is the HOCO complex, whereas the formation of CO2 has a probability of only (3.6 ± 0.7) × 10(-4). The initial location of the CO is key to obtain reaction and form CO2: the CO needs to be located deep into the ice. The HOCO complex becomes trapped in the cold ice surface in the trans-HOCO minimum because it quickly loses its internal energy to the surrounding ice, preventing further reaction to H + CO2. Several laboratory experiments have been carried out recently, and they confirm that CO2 can also be formed through other, different routes. Here we compare our theoretical results with the data available from experiments studying the formation of CO2 through a similar pathway as ours, even though the initial conditions were not exactly the same. Our results also show that the HCO van der Waals complex can be formed through the interaction of CO with the H atom that is formed as a product of H2O photodissociation. Thus, the reaction of the H atom photofragment following H2O photodissociation with CO can be a possible route to form HCO ice.

Molecular dynamics simulations of D2O ice photodesorption.

Molecular dynamics (MD) calculations have been performed to study the ultraviolet (UV) photodissociation of D(2)O in an amorphous D(2)O ice surface at 10, 20, 60, and 90 K, in order to investigate the influence of isotope effects on the photodesorption processes. As for H(2)O, the main processes after UV photodissociation are trapping and desorption of either fragments or D(2)O molecules. Trapping mainly takes place in the deeper monolayers of the ice, whereas desorption occurs in the uppermost layers. There are three desorption processes: D atom, OD radical, and D(2)O molecule photodesorption. D(2)O desorption takes places either by direct desorption of a recombined D(2)O molecule, or when an energetic D atom produced by photodissociation kicks a surrounding D(2)O molecule out of the surface by transferring part of its momentum. Desorption probabilities are calculated for photoexcitation of D(2)O in the top four monolayers and are compared quantitatively with those for H(2)O obtained from previous MD simulations of UV photodissociation of amorphous water ice at different ice temperatures [Arasa et al., J. Chem. Phys. 132, 184510 (2010)]. The main conclusions are the same, but the average D atom photodesorption probability is smaller than that of the H atom (by about a factor of 0.9) because D has lower kinetic energy than H, whereas the average OD radical photodesorption probability is larger than that of OH (by about a factor of 2.5-2.9 depending on ice temperature) because OD has higher translational energy than OH for every ice temperature studied. The average D(2)O photodesorption probability is larger than that of H(2)O (by about a factor of 1.4-2.3 depending on ice temperature), and this is entirely due to a larger contribution of the D(2)O kick-out mechanism. This is an isotope effect: the kick-out mechanism is more efficient for D(2)O ice, because the D atom formed after D(2)O photodissociation has a larger momentum than photogenerated H atoms from H(2)O, and D transfers momentum more easily to D(2)O than H to H(2)O. The total (OD + D(2)O) yield has been compared with experiments and the total (OH + H(2)O) yield from previous simulations. We find better agreement when we compare experimental yields with calculated yields for D(2)O ice than when we compare with calculated yields for H(2)O ice.

Water formation by surface O3 hydrogenation.

Three solid state formation routes have been proposed in the past to explain the observed abundance of water in space: the hydrogenation reaction channels of atomic oxygen (O + H), molecular oxygen (O(2) + H), and ozone (O(3) + H). New data are presented here for the third scheme with a focus on the reactions O(3) + H, OH + H and OH + H(2), which were difficult to quantify in previous studies. A comprehensive set of H/D-atom addition experiments is presented for astronomically relevant temperatures. Starting from the hydrogenation/deuteration of solid O(3) ice, we find experimental evidence for H(2)O/D(2)O (and H(2)O(2)/D(2)O(2)) ice formation using reflection absorption infrared spectroscopy. The temperature and H/D-atom flux dependence are studied and this provides information on the mobility of ozone within the ice and possible isotope effects in the reaction scheme. The experiments show that the O(3) + H channel takes place through stages that interact with the O and O(2) hydrogenation reaction schemes. It is also found that the reaction OH + H(2) (OH + H), as an intermediate step, plays a prominent (less efficient) role. The main conclusion is that solid O(3) hydrogenation offers a potential reaction channel for the formation of water in space. Moreover, the nondetection of solid ozone in dense molecular clouds is consistent with the astrophysical picture in which O(3) + H is an efficient process under interstellar conditions.

Water formation at low temperatures by surface O2 hydrogenation I: Characterization of ice penetration.

Water is the main component of interstellar ice mantles, is abundant in the solar system and is a crucial ingredient for life. The formation of this molecule in the interstellar medium cannot be explained by gas-phase chemistry only and its surface hydrogenation formation routes at low temperatures (O, O(2), O(3) channels) are still unclear and most likely incomplete. In a previous paper we discussed an unexpected zeroth-order H(2)O production behavior in O(2) ice hydrogenation experiments compared to the first-order H(2)CO and CH(3)OH production behavior found in former studies on hydrogenation of CO ice. In this paper we experimentally investigate in detail how the structure of O(2) ice leads to this rare behavior in reaction order and production yield. In our experiments H atoms are added to a thick O(2) ice under fully controlled conditions, while the changes are followed by means of reflection absorption infrared spectroscopy (RAIRS). The H-atom penetration mechanism is systematically studied by varying the temperature, thickness and structure of the O(2) ice. We conclude that the competition between reaction and diffusion of the H atoms into the O(2) ice explains the unexpected H(2)O and H(2)O(2) formation behavior. In addition, we show that the proposed O(2) hydrogenation scheme is incomplete, suggesting that additional surface reactions should be considered. Indeed, the detection of newly formed O(3) in the ice upon H-atom exposure proves that the O(2) channel is not an isolated route. Furthermore, the addition of H(2) molecules is found not to have a measurable effect on the O(2) reaction channel.

Sticking of CO to crystalline and amorphous ice surfaces.

We present results of classical trajectory calculations on the sticking of hyperthermal CO to the basal plane (0001) face of crystalline ice Ih and to the surface of amorphous ice Ia. The calculations were performed for normal incidence at a surface temperature Ts = 90 K for ice Ia, and at Ts = 90 and 150 K for ice Ih. For both surfaces, the sticking probability can be fitted to a simple exponentially decaying function of the incidence energy, Ei: Ps = 1.0e(-Ei(kJ/mol)/90(kJ/mol)) at Ts = 90 K. The energy transfer from the impinging molecule to the crystalline and the amorphous surface is found to be quite efficient, in agreement with the results of molecular beam experiments on the scattering of the similar molecule, N2, from crystalline and amorphous ice. However, the energy transfer is less efficient for amorphous than for crystalline ice. Our calculations predict that the sticking probability decreases with Ts for CO scattering from crystalline ice, as the energy transfer from the impinging molecule to the warmer surfaces becomes less efficient. At high Ei (up to 193 kJ/mol), no surface penetration occurs in the case of crystalline ice. However, for CO colliding with the amorphous surface, a penetrating trajectory was observed to occur into a large water pore. The molecular dynamics calculations predict that the average potential energy of CO adsorbed to ice Ih is -10.1 +/- 0.2 and -8.4 +/- 0.2 kJ/mol for CO adsorbed to ice Ia. These values are in agreement with previous experimental and theoretical data. The distribution of the potential energy of CO adsorbed to ice Ia was found to be wider (with a standard deviation sigma of 2.4 kJ/mol) than that of CO interacting with ice Ih (sigma = 2.0 kJ/mol). In collisions with ice Ia, the CO molecules scatter at larger angles and over a wider distribution of angles than in collisions with ice Ih.

Substantial reservoirs of molecular hydrogen in the debris disks around young stars.

Circumstellar accretion disks transfer matter from molecular clouds to young stars and to the sites of planet formation. The disks observed around pre-main-sequence stars have properties consistent with those expected for the pre-solar nebula from which our own Solar System formed 4.5 Gyr ago. But the 'debris' disks that encircle more than 15% of nearby main-sequence stars appear to have very small amounts of gas, based on observations of the tracer molecule carbon monoxide: these observations have yielded gas/dust ratios much less than 0.1, whereas the interstellar value is about 100 (ref. 9). Here we report observations of the lowest rotational transitions of molecular hydrogen (H2) that reveal large quantities of gas in the debris disks around the stars beta Pictoris, 49 Ceti and HD135344. The gas masses calculated from the data are several hundreds to a thousand times greater than those estimated from the CO observations, and yield gas/dust ratios of the same order as the interstellar value.

Envelope structure of deeply embedded young stellar objects in the Serpens Molecular Cloud.

Aperture-synthesis and single-dish (sub-) millimeter molecular-line and continuum observations reveal in great detail the envelope structure of deeply embedded young stellar objects (SMM 1 = FIRS 1, SMM 2, SMM 3, SMM 4) in the densely star-forming Serpens Molecular Cloud. SMM 1, 3, and 4 show partially resolved (>2" = 800 AU) continuum emission in the beam of the Owens Valley Millimeter Array at lambda = 3.4-1.4 mm. The continuum visibilities accurately constrain the density structure in the envelopes, which can be described by a radial power law with slope -2.0 +/- 0.5 on scales of 300 to 8000 AU. Inferred envelope masses within a radius of 8000 AU are 8.7, 3.0, and 5.3 Msolar for SMM 1, 3, and 4, respectively. A point source with 20%-30% of the total flux at 1.1 mm is required to fit the observations on long baselines, corresponding to warm envelope material within approximately 100 AU or a circumstellar disk. No continuum emission is detected interferometrically toward SMM 2, corresponding to an upper limit of 0.2 Msolar assuming Td = 24 K. The lack of any compact dust emission suggests that the SMM 2 core does not contain a central protostar. Aperture-synthesis observations of the 13CO, C18O, HCO+, H13CO+, HCN, H13CN, N2H+ 1-0, SiO 2-1, and SO 2(2)-1(1) transitions reveal compact emission toward SMM 1, 3, and 4. SMM 2 shows only a number of clumps scattered throughout the primary field of view, supporting the conclusion that this core does not contain a central star. The compact molecular emission around SMM 1, 3, and 4 traces 5"-10" (2000-4000 AU) diameter cores that correspond to the densest regions of the envelopes, as well as material directly associated with the molecular outflow. Especially prominent are the optically thick HCN and HCO+ lines that show up brightly along the walls of the outflow cavities. SO and SiO trace shocked material, where their abundances may be enhanced by 1-2 orders of magnitude over dark-cloud values. A total of 31 molecular transitions have been observed with the James Clerk Maxwell and Caltech Submillimeter telescopes in the 230, 345, 490, and 690 GHz atmospheric windows toward all four sources, containing, among others, lines of CO, HCO+, HCN, H2CO, SiO, SO, and their isotopomers. These lines show 20-30 km s-1 wide line wings, deep and narrow (1-2 km s-1) self-absorption, and 2-3 km s-1 FWHM line cores. The presence of highly excited lines like 12CO 4-3 and 6-5, 13CO 6-5, and several H2CO transitions indicates the presence of material with temperatures > or approximately 100 K. Monte Carlo calculations of the molecular excitation and line transfer show that the envelope model derived from the dust emission can successfully reproduce the observed line intensities. The depletion of CO in the cold gas is modest compared to values inferred in objects like NGC 1333 IRAS 4, suggesting that the phase of large depletions through the entire envelope is short lived and may be influenced by the local star formation density. Emission in high-excitation lines of CO and H2CO requires the presence of a small amount of approximately 100 K material, comprising less than 1% of the total envelope mass and probably associated with the outflow or the innermost region of the envelope. The derived molecular abundances in the warm (Tkin > 20 K) envelope are similar to those found toward other class 0 YSOs like IRAS 16293-2422, though some species appear enhanced toward SMM 1. Taken together, the presented observations and analysis provide the first comprehensive view of the physical and chemical structure of the envelopes of deeply embedded young stellar objects in a clustered environment on scales between 1000 and 10,000 AU.

Envelope structure on 700 AU scales and the molecular outflows of low-mass young stellar objects.

Aperture synthesis observations of HCO+ J = 1-0, 13CO 1-0, and C18O 1-0 obtained with the Owens Valley Millimeter Array are used to probe the small-scale (5" approximately 700 AU) structure of the molecular envelopes of a well-defined sample of nine embedded low-mass young stellar objects in Taurus. The interferometer results can be understood in terms of: (1) a core of radius approximately or less than 1000 AU surrounding the central star, possibly flattened and rotating; (2) condensations scattered throughout the envelope that may be left over from the inhomogeneous structure of the original cloud core or that may have grown during collapse; and (3) material within the outflow or along the walls of the outflow cavity. Masses of the central cores are 0.001-0.1 M (solar), and agree well with dust continuum measurements. Averaged over the central 20" (3000 AU) region, an HCO+ abundance of 4 x 10(-8) is inferred, with a spread of a factor of 3 between the different sources. Reanalysis of previously presented single-dish data yields an HCO+ abundance of (5.0 +/- 1.7) x 10(-9), which may indicate an average increase by a factor of a few on the smaller scales sampled by the interferometer. Part of this apparent abundance variation could be explained by contributions from extended cloud emission to the single-dish C18O lines, and uncertainties in the assumed excitation temperatures and opacities. The properties of the molecular envelopes and outflows are further investigated through single-dish observations of 12CO J = 6-5, 4-3, and 3-2, 13CO 6-5 and 3-2, and C18O 3-2 and 2-1, obtained with the James Clerk Maxwell and IRAM 30 m telescopes, along with the Caltech Submillimeter Observatory. Ratios of the mid-J CO lines are used to estimate the excitation temperature, with values of 25-80 K derived for the gas near line centre. The outflow wings show a similar range, although Tex is enhanced by a factor of 2-3 in at least two sources. In contrast to the well-studied L1551 IRS 5 outflow, which extends over 10' (0.4 pc), seven of the remaining eight sources are found to drive 12CO 3-2 outflows over < or = 1' (0.04 pc); only L1527 IRS has a well-developed outflow of some 3'(0.12 pc). Estimates are obtained for the outflow kinetic luminosity, Lkin, and the flow momentum rate, FCO, applying corrections for line opacity and source inclination. The flow force FCO correlates with the envelope mass and with the 2.7 mm flux of the circumstellar disk. Only a weak correlation is seen with Lbol, while none is found with the relative age of the object as measured by integral Tmb(HCO+ 3-2)dV/Lbol. These trends support the hypothesis that outflows are driven by accretion through a disk, with a global mass infall rate determined by the mass and density of the envelope. The association of compact HCO+ emission with the walls of the outflow cavities indicates that outflows in turn influence the appearance of the envelopes. It is not yet clear, however, whether they are actively involved in sweeping up envelope material, or merely provide a low-opacity pathway for heating radiation to reach into the envelope.

Chemical evolution of star-forming regions.

Recent advances in the understanding of the chemical processes that occur during all stages of the formation of stars, from the collapse of molecular clouds to the assemblage of icy planetesimals in protoplanetary accretion disks, are reviewed. Observational studies of the circumstellar material within 100-10,000 AU of the young star with (sub)millimeter single-dish telescopes, millimeter interferometers, and ground-based as well as space-borne infrared observatories have only become possible within the past few years. Results are compared with detailed chemical models that emphasize the coupling of gas-phase and grain-surface chemistry. Molecules that are particularly sensitive to different routes of formation and that may be useful in distinguishing between a variety of environments and histories are outlined. In the cold, low-density prestellar cores, radicals and long unsaturated carbon chains are enhanced. During the cold collapse phase, most species freeze out onto the grains in the high-density inner region. Once young stars ignite, their surroundings are heated through radiation and/or shocks, whereupon new chemical characteristics appear. Evaporation of ices drives a ''hot core'' chemistry rich in organic molecules, whereas shocks propagating through the dense envelope release both refractory and volatile grain material, resulting in prominent SiO, OH, and H2O emission. The role of future instrumentation in further developing these chemical and temporal diagnostics is discussed.

A molecular line study of NGC 1333/IRAS 4.

Molecular line surveys and fully sampled spectral line maps at 1.3 and 0.87 mm are used to examine the physical and chemical characteristics of the extreme Class I sources IRAS 4A and 4B in the L1450/NGC 1333 molecular cloud complex. A very well collimated, jetlike molecular outflow emanates from IRAS 4A, with a dynamical age of a few thousand years. Symmetric, clumpy structure along the outflow lobes suggests that there is considerable variability in the mass-loss rate or wind velocity even at this young age. Molecular emission lines toward IRAS 4A and 4B are observed to be weak in the velocity range corresponding to quiescent material surrounding the young stellar objects (YSOs). Depletion factors of 10-20 are observed for all molecules, including CO, even for even for very conservative mass estimates from the measured millimeter and submillimeter dust continuum. However, abundances scaled with respect to CO are similar to other dark molecular cloud cores. Such depletions could be mimicked by high dust optical depths or increased grain emissivities at the observing frequencies of 230 and 345 GHz, but the millimeter and submillimeter spectral energy distributions suggest that this is unlikely over the single-dish size scales of 5000-10,000 AU. Dense, outflowing gas is found to be kinematically, but not spatially, distinct from the quiescent material on these size scales. If CO is used as a chemical standard for the high-velocity gas, we find substantial enhancements in the abundances of several molecules in outflowing material, most notably CS, SiO, and CH3OH. The SiO emission is kinematically well displaced from the bulk cloud velocity and likely arises from directly shocked material. As is the case for CO, however, the outflow features from more volatile species are centered near the cloud velocity and are often characterized by quite low rotational temperatures. We suggest that grain-grain collisions induced by velocity shear zones surrounding the outflow axes transiently desorb the grain mantles, resulting in large abundance enhancements of selected species. Similar results have recently been obtained in several other low-mass YSOs, where the outflowing gas is often both kinematically and spatially distinct, and are illustrative of the ability of accretion and outflow processes to simultaneously modify the composition of the gas and dust surrounding young stars.

Chemical evolution of circumstellar matter around young stellar objects.

Recent observational studies of the chemical composition of circumstellar matter around both high- and low-mass young stellar objects are reviewed. The molecular abundances are found to be a strong function of evolutionary state, but not of system mass or luminosity. The data are discussed with reference to recent theoretical models.