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 parsecs 1,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 protostars 3. 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 components 3. 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\ua0reveal copious emission from hot molecules, explaining the origin of the ‘green fuzzies’ 4–7 discovered nearly two decades ago by the Spitzer Space Telescope 8. 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

JOYS+: mid-infrared detection of gas-phase SO2_2 emission in a low-mass protostar. The case of NGC 1333 IRAS2A: hot core or accretion shock?

JWST/MIRI has sharpened our infrared eyes toward the star formation process. This paper presents the first mid-infrared detection of gaseous SO$_2$ emission in an embedded low-mass protostellar system. MIRI-MRS observations of the low-mass protostellar binary NGC 1333 IRAS2A are presented from the JWST Observations of Young protoStars (JOYS+) program, revealing emission from the SO$_2~\nu_3$ asymmetric stretching mode at 7.35 micron. The results are compared to those derived from high-angular resolution SO$_2$ data obtained with ALMA. The SO$_2$ emission from the $\nu_3$ band is predominantly located on $\sim50-100$ au scales around the main component of the binary, IRAS2A1. A rotational temperature of $92\pm8$ K is derived from the $\nu_3$ lines. This is in good agreement with the rotational temperature derived from pure rotational lines in the vibrational ground state (i.e., $\nu=0$) with ALMA ($104\pm5$ K). However, the emission of the $\nu_3$ lines is not in LTE given that the total number of molecules predicted by a LTE model is found to be a factor $2\times10^4$ higher than what is derived for the $\nu=0$ state. This difference can be explained by a vibrational temperature that is $\sim100$ K higher than the derived rotational temperature of the $\nu=0$ state. The brightness temperature derived from the continuum around the $\nu_3$ band of SO$_2$ is $\sim180$ K, which confirms that the $\nu_3=1$ level is not collisionally populated but rather infrared pumped by scattered radiation. This is also consistent with the non-detection of the $\nu_2$ bending mode at 18-20 micron. Given the rotational temperature, the extent of the emission ($\sim100$ au in radius), and the narrow line widths in the ALMA data (3.5 km/s), the SO$_2$ in IRAS2A likely originates from ice sublimation in the central hot core around the protostar rather than from an accretion shock at the disk-envelope boundary.Comment: 19 pages, 17 figures, accepted for publication in A&A, abstract abbreviate

JOYS: JWST Observations of Young protoStars: Outflows and accretion in the high-mass star-forming region IRAS23385+605

Aims: The JWST program JOYS (JWST Observations of Young protoStars) aims at characterizing the physical and chemical properties of young high- and low-mass star-forming regions, in particular the unique mid-infrared diagnostics of the warmer gas and solid-state components. We present early results from the high-mass star formation region IRAS23385+6053. Methods: The JOYS program uses the MIRI MRS with its IFU to investigate a sample of high- and low-mass star-forming protostellar systems. Results: The 5 to 28mum MIRI spectrum of IRAS23385+6053 shows a plethora of features. While the general spectrum is typical for an embedded protostar, we see many atomic and molecular gas lines boosted by the higher spectral resolution and sensitivity compared to previous space missions. Furthermore, ice and dust absorption features are also present. Here, we focus on the continuum emission, outflow tracers like the H2, [FeII] and [NeII] lines as well as the potential accretion tracer Humphreys alpha HI(7--6). The short-wavelength MIRI data resolve two continuum sources A and B, where mid-infrared source A is associated with the main mm continuum peak. The combination of mid-infrared and mm data reveals a young cluster in its making. Combining the mid-infrared outflow tracer H2, [FeII] and [NeII] with mm SiO data shows a complex interplay of at least three molecular outflows driven by protostars in the forming cluster. Furthermore, the Humphreys alpha line is detected at a 3-4sigma level towards the mid-infrared sources A and B. Following Rigliaco et al. (2015), one can roughly estimate accretion luminosities and corresponding accretion rates between ~2.6x10^-6 and ~0.9x10^-4 M_sun/yr. This is discussed in the context of the observed outflow rates. Conclusions: The analysis of the MIRI MRS observations for this young high-mass star-forming region reveals connected outflow and accretion signatures.Comment: 12 pages, 9 figures, accepted for Astronomy & Astrophysics, the paper is also available at https://www2.mpia-hd.mpg.de/homes/beuther/papers.htm

Fragmentation and disk formation during high-mass star formation: The IRAM NOEMA (Northern Extended Millimeter Array) large program CORE

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation

Multi-scale view of star formation in IRAS 21078+5211: From clump fragmentation to disk wind

In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we derive a mass infall rate of about 10^{-4}~M_{sun}~yr^{-1}. The most massive cores (labeled 1, 2, and 3) are found at the center of the cluster, and these are the only ones that present a signature of protostellar activity in terms of emission from high-excitation molecular lines or a molecular outflow. We reveal an extended (size about 0.1~pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by a (previously discovered) compact (size <= 1000~au) radio jet, ejected by a YSO embedded in core 1 (named YSO-1), since the molecular outflow and the radio jet are almost parallel and have a comparable momentum rate. By means of high-excitation lines, we find a large (14~km/s over 500~au) Vlsr gradient at the position of YSO-1, oriented approximately perpendicular to the radio jet. Assuming this is an edge-on, rotating disk and fitting a Keplerian rotation pattern, we determine the YSO-1 mass to be 5.6+/-2.0~M_{sun}. The water masers (previously observed with VLBI) emerge within 100-300~au from YSO-1 and are unique tracers of the jet kinematics. Their three-dimensional (3D) velocity pattern reveals that the gas flows along, and rotates about, the jet axis. We show that the 3D maser velocities are fully consistent with the magneto-centrifugal disk-wind models predicting a cylindrical rotating jet. Under this hypothesis, we determine the jet radius to be about 16~au and the corresponding launching radius and terminal velocity to be about 2.2~au and 200~km/s, respectively

Disk fragmentation in high-mass star formation. High-resolution observations towards AFGL 2591-VLA 3

Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk around the high-mass star forming region AFGL 2591-VLA 3 which was hypothesized to be fragmenting based on the observations that show multiple outflow directions. We use a new set of high-resolution (0.19 arcsec) IRAM/NOEMA observations at 843 micron towards VLA 3 which allow us to resolve its disk, characterize the fragmentation, and study its kinematics. In addition to the 843 micron continuum emission, our spectral setup targets warm dense gas and outflow tracers such as HCN, HC$_3$N and SO$_2$, as well as vibrationally excited HCN lines. The high resolution continuum and line emission maps reveal multiple fragments with subsolar masses within the inner 1000 AU of VLA 3. Furthermore, the velocity field of the inner disk observed at 843 micron shows a similar behavior to that of the larger scale velocity field studied in the CORE project at 1.37 mm. We present the first observational evidence for disk fragmentation towards AFGL 2591-VLA 3, a source that was thought to be a single high-mass core. While the fragments themselves are low-mass, the rotation of the disk is dominated by the protostar with a mass of 10.3$\pm 1.8~M_{\odot}$. These data also show that NOEMA Band 4 can obtain the highest currently achievable spatial resolution at (sub-)mm wavelengths in observations of strong northern sources

Chemical complexity in high-mass star formation: An observational and modeling case study of the AFGL 2591 VLA 3 hot core

Aims. In order to understand the observed molecular diversity in high-mass star-forming regions, we have to determine the underlying physical and chemical structure of those regions at high angular resolution and over a range of evolutionary stages. We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth

Transcriptional Regulation of Rod Photoreceptor Homeostasis Revealed by In Vivo NRL Targetome Analysis

A stringent control of homeostasis is critical for functional maintenance and survival of neurons. In the mammalian retina, the basic motif leucine zipper transcription factor NRL determines rod versus cone photoreceptor cell fate and activates the expression of many rod-specific genes. Here, we report an integrated analysis of NRL-centered gene regulatory network by coupling chromatin immunoprecipitation followed by high-throughput sequencing (ChIP–Seq) data from Illumina and ABI platforms with global expression profiling and in vivo knockdown studies. We identified approximately 300 direct NRL target genes. Of these, 22 NRL targets are associated with human retinal dystrophies, whereas 95 mapped to regions of as yet uncloned retinal disease loci. In silico analysis of NRL ChIP–Seq peak sequences revealed an enrichment of distinct sets of transcription factor binding sites. Specifically, we discovered that genes involved in photoreceptor function include binding sites for both NRL and homeodomain protein CRX. Evaluation of 26 ChIP–Seq regions validated their enhancer functions in reporter assays. In vivo knockdown of 16 NRL target genes resulted in death or abnormal morphology of rod photoreceptors, suggesting their importance in maintaining retinal function. We also identified histone demethylase Kdm5b as a novel secondary node in NRL transcriptional hierarchy. Exon array analysis of flow-sorted photoreceptors in which Kdm5b was knocked down by shRNA indicated its role in regulating rod-expressed genes. Our studies identify candidate genes for retinal dystrophies, define cis-regulatory module(s) for photoreceptor-expressed genes and provide a framework for decoding transcriptional regulatory networks that dictate rod homeostasis

Fragmentation and disk formation during high-mass star formation. IRAM NOEMA (Northern Extended Millimeter Array) large program CORE

Context. High-mass stars form in clusters, but neither the early fragmentation processes nor the detailed physical processes leading to the most massive stars are well understood. Aims. We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during highmass star formation on spatial scales of individual cores. Methods. Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30 m telescope, we have observed in the IRAM large program CORE the 1.37 mm continuum and spectral line emission at high angular resolution (?0.4??) for a sample of 20 well-known high-mass star-forming regions with distances below 5.5 kpc and luminosities larger than 104 L?. Results. We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000 AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions. A large sample of high-mass star-forming regions at high spatial resolution allows us to study the fragmentation properties of young cluster-forming regions. The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation

Kinematics and stability of high-mass protostellar disk candidates at sub-arcsecond resolution

Context. The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars (M & 8 M) are still not well understood. A growing number of case studies have found massive young stellar objects (MYSOs) to harbour disk-like structures, painting a picture that the formation of high-mass stars may proceed through disk accretion, similar to that of lower mass stars. However, the properties of such structures have yet to be uniformly and systematically characterised. Massive disks are prone to fragmentation via gravitational instabilities due to high gas densities and accretion rates. Therefore, it is important to study the stability of such disks in order to put into context the role of disk fragmentation in setting the final stellar mass distribution in high-mass star forming regions. Aims. The aim of this work is to uniformly study the kinematic properties of a large sample of MYSOs and characterise the stability of possible circumstellar disks against gravitational fragmentation. Methods. We have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous (L > 104 L) protostellar objects in the 1.37 mm wavelength regime in both continuum and spectral line emission, reaching 0.400 resolution (800 au at 2 kpc). Results. We present the gas kinematics of the full sample and detect dense gas emission surrounding 15 regions within the CORE sample. Using the dense gas tracer CH3CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. The extent of the CH3CN emission tracing the disk candidates varies from 1800 − 8500 au. Analysing the free-fall to rotational timescales, we find that the sources are rotationally supported. The rotation profiles of some disk candidates are well described by differential rotation while for others the profiles are poorly resolved. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of ∼ 10 − 25 M. Modelling the level population of CH3CN (12K − 11K) K = 0 − 6 lines we present temperature maps and find median temperature in the range 70–210 K with a diversity in distributions. Radial profiles of the specific angular momentum (j) for the best disk candidates span a range of 1–2 orders of magnitude, on average ∼ 10−3 km s−1 pc, and follow j ∝ r 1.7, consistent with a poorly resolved rotating and infalling envelope/disk model. Studying the Toomre stability of the disk candidates, we find almost all (11 out of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations, as a result of their high disk to stellar mass ratio. In particular, disks with masses greater than ∼ 10 − 20% of the mass of their host (proto)stars are Toomre unstable, and more luminous YSOs tend to have disks that are more massive compared to their host star and hence more prone to fragmentation. Conclusions. In this work, we show that most disk structures around high-mass YSOs are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars