Episodic accretion: the interplay of infall and disc instabilities

Using zoom-simulations carried out with the adaptive mesh-refinement code RAMSES with a dynamic range of up to $2^{27} \approx 1.34 \times 10^8$ we investigate the accretion profiles around six stars embedded in different environments inside a (40 pc)$^3$ giant molecular cloud, the role of mass infall and disc instabilities on the accretion profile, and thus on the luminosity of the forming protostar. Our results show that the environment in which the protostar is embedded determines the overall accretion profile of the protostar. Infall on to the circumstellar disc may trigger gravitational disc instabilities in the disc at distances of around ~10 to ~50 au leading to rapid transport of angular momentum and strong accretion bursts. These bursts typically last for about ~10 to a ~100 yr, consistent with typical orbital times at the location of the instability, and enhance the luminosity of the protostar. Calculations with the stellar evolution code mesa show that the accretion bursts induce significant changes in the protostellar proper- ties, such as the stellar temperature and radius. We apply the obtained protostellar properties to produce synthetic observables with RADMC3D and predict that accretion bursts lead to ob- servable enhancements around 20 to 200 $\mu$m in the spectral energy distribution of Class 0 type young stellar objects.Comment: 17 pages, 14 figures, accepted by MNRA

TIPSY: Trajectory of Infalling Particles in Streamers around Young stars. Dynamical analysis of the streamers around S CrA and HL Tau

Context. Elongated trails of infalling gas, often referred to as "streamers," have recently been observed around young stellar objects (YSOs) at different evolutionary stages. This asymmetric infall of material can significantly alter star and planet formation processes, especially in the more evolved YSOs. Aims. In order to ascertain the infalling nature of observed streamer-like structures and then systematically characterize their dynamics, we developed the code TIPSY (Trajectory of Infalling Particles in Streamers around Young stars). Methods. Using TIPSY, the streamer molecular line emission is first isolated from the disk emission. Then the streamer emission, which is effectively a point cloud in three-dimensional (3D) position-position-velocity space, is simplified to a curve-like representation. The observed streamer curve is then compared to the theoretical trajectories of infalling material. The best-fit trajectories are used to constrain streamer features, such as the specific energy, the specific angular momenta, the infall timescale, and the 3D morphology. Results. We used TIPSY to fit molecular-line ALMA observations of streamers around a Class II binary system, S CrA, and a Class I/II protostar, HL Tau. Our results indicate that both of the streamers are consistent with infalling motion. TIPSY results and mass estimates suggest that S CrA and HL Tau are accreting material at a rate of $\gtrsim27$ M$_{jupiter}$ Myr$^{-1}$ and $\gtrsim5$ M$_{jupiter}$ Myr$^{-1}$, respectively, which can significantly increase the mass budget available to form planets. Conclusions. TIPSY can be used to assess whether the morphology and kinematics of observed streamers are consistent with infalling motion and to characterize their dynamics, which is crucial for quantifying their impact on the protostellar systems.Comment: Accepted in Astronomy & Astrophysic

From bubbles and filaments to cores and disks: gas gathering and growth of structure leading to the formation of stellar systems

35 pages, 20 figures. To appear in Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Small typographic changes, and color in referencesThe study of the development of structures on multiple scales in the cold interstellar medium has experienced rapid expansion in the past decade, on both the observational and the theoretical front. Spectral line studies at (sub-)millimeter wavelengths over a wide range of physical scales have provided unique probes of the kinematics of dense gas in star-forming regions, and have been complemented by extensive, high dynamic range dust continuum surveys of the column density structure of molecular cloud complexes, while dust polarization maps have highlighted the role of magnetic fields. This has been accompanied by increasingly sophisticated numerical simulations including new physics (e.g., supernova driving, cosmic rays, non-ideal magneto-hydrodynamics, radiation pressure) and new techniques such as zoom-in simulations allowing multi-scale studies. Taken together, these new data have emphasized the anisotropic growth of dense structures on all scales, from giant ISM bubbles driven by stellar feedback on $\sim$50--100 pc scales through parsec-scale molecular filaments down to $< 0.1$ pc dense cores and $< 1000$ au protostellar disks. Combining observations and theory, we present a coherent picture for the formation and evolution of these structures and synthesize a comprehensive physical scenario for the initial conditions and early stages of star and disk formation

From bubbles and filaments to cores and disks: gas gathering and growth of structure leading to the formation of stellar systems

35 pages, 20 figures. To appear in Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Small typographic changes, and color in referencesThe study of the development of structures on multiple scales in the cold interstellar medium has experienced rapid expansion in the past decade, on both the observational and the theoretical front. Spectral line studies at (sub-)millimeter wavelengths over a wide range of physical scales have provided unique probes of the kinematics of dense gas in star-forming regions, and have been complemented by extensive, high dynamic range dust continuum surveys of the column density structure of molecular cloud complexes, while dust polarization maps have highlighted the role of magnetic fields. This has been accompanied by increasingly sophisticated numerical simulations including new physics (e.g., supernova driving, cosmic rays, non-ideal magneto-hydrodynamics, radiation pressure) and new techniques such as zoom-in simulations allowing multi-scale studies. Taken together, these new data have emphasized the anisotropic growth of dense structures on all scales, from giant ISM bubbles driven by stellar feedback on $\sim$50--100 pc scales through parsec-scale molecular filaments down to $< 0.1$ pc dense cores and $< 1000$ au protostellar disks. Combining observations and theory, we present a coherent picture for the formation and evolution of these structures and synthesize a comprehensive physical scenario for the initial conditions and early stages of star and disk formation

Probing the physics of star formation (ProPStar)

Context. The detections of narrow channels of accretion toward protostellar disks, known as streamers, have increased in number in the last few years. However, it is unclear whether streamers are a common feature around protostars that were previously missed, or if they are a rare phenomenon. Aims. Our goals are to obtain the incidence of streamers toward a region of clustered star formation and to trace the origins of their gas to determine whether they originate within the filamentary structure of molecular clouds or from beyond. Methods. We used combined observations of the nearby NGC 1333 star-forming region, carried out with the NOEMA interferometer and the IRAM 30m single dish. Our observations cover the area between the systems IRAS 4 and SVS 13. We traced the chemically fresh gas within NGC 1333 with HC3N molecular gas emission and the structure of the fibers in this region with N2H+ emission. We fit multiple velocity components in both maps and used clustering algorithms to recover velocity-coherent structures. Results. We find streamer candidates toward 7 out of 16 young stellar objects within our field of view. This represents an incidence of approximately 40% of young stellar objects with streamer candidates in a clustered star-forming region. The incidence increases to about 60% when we only considered embedded protostars. All streamers are found in HC3N emission. Conclusions. Given the different velocities between HC3N and N2H+ emission, and because by construction, N2H+ traces the fiber structure, we suggest that the gas that forms the streamers comes from outside the fibers. This implies that streamers can connect cloud material that falls onto the filaments with protostellar disk scales

Probing the physics of star formation (ProPStar)

Context. Turbulence is a key component of molecular cloud structure. It is usually described by a cascade of energy down to the dissipation scale. The power spectrum for subsonic incompressible turbulence is ∝k−5/3, while for supersonic turbulence it is ∝k−2. Aims. We determine the power spectrum in an actively star-forming molecular cloud, from parsec scales down to the expected magnetohydrodynamic (MHD) wave cutoff (dissipation scale). Methods. We analyzed observations of the nearby NGC 1333 star-forming region in three different tracers to cover the different scales from ∼10 pc down to 20 mpc. The largest scales are covered with the low-density gas tracer 13CO (1–0) obtained with a single dish, the intermediate scales are covered with single-dish observations of the C18O (3–2) line, while the smallest scales are covered in H13CO+ (1–0) and HNC (1–0) with a combination of NOEMA interferometer and IRAM 30m single-dish observations. The complementarity of these observations enables us to generate a combined power spectrum covering more than two orders of magnitude in spatial scale. Results. We derive the power spectrum in an active star-forming region spanning more than 2 decades of spatial scales. The power spectrum of the intensity maps shows a single power-law behavior, with an exponent of 2.9 ± 0.1 and no evidence of dissipation. Moreover, there is evidence that the power spectrum of the ions to have more power at smaller scales than the neutrals, which is opposite to the theoretical expectations. Conclusions. We show new possibilities for studying the dissipation of energy at small scales in star-forming regions provided by interferometric observations

Probing the physics of star formation (ProPStar)

Context. Electron fraction and cosmic-ray ionization rates in star-forming regions are important quantities in astrochemical modeling and are critical to the degree of coupling between neutrals, ions, and electrons, which regulates the dynamics of the magnetic field. However, these are difficult quantities to estimate. Aims. We aim to derive the electron fraction and cosmic-ray ionization rate maps of an active star-forming region. Methods. We combined observations of the nearby NGC 1333 star-forming region carried out with the NOEMA interferometer and IRAM 30 m single dish to generate high spatial dynamic range maps of different molecular transitions. We used the DCO+ and H13CO+ ratio (in addition to complementary data) to estimate the electron fraction and produce cosmic-ray ionization rate maps. Results. We derived the first large-area electron fraction and cosmic-ray ionization rate resolved maps in a star-forming region, with typical values of 10−65 and 10−16.5 s−1, respectively. The maps present clear evidence of enhanced values around embedded young stellar objects (YSOs). This provides strong evidence for locally accelerated cosmic rays. We also found a strong enhancement toward the northwest region in the map that might be related either to an interaction with a bubble or to locally generated cosmic rays by YSOs. We used the typical electron fraction and derived a magnetohydrodynamic (MHD) turbulence dissipation scale of 0.054 pc, which could be tested with future observations. Conclusions. We found a higher cosmic-ray ionization rate compared to the canonical value for N(H2) = 1021−1023 cm−2 of 10−17 s−1 in the region, and it is likely generated by the accreting YSOs. The high value of the electron fraction suggests that new disks will form from gas in the ideal-MHD limit. This indicates that local enhancements of ζ(H2), due to YSOs, should be taken into account in the analysis of clustered star formation

Isotopic enrichment of forming planetary systems from supernova pollution

Heating by short-lived radioisotopes (SLRs) such as aluminum-26 and iron-60 fundamentally shaped the thermal history and interior structure of Solar System planetesimals during the early stages of planetary formation. The subsequent thermo-mechanical evolution, such as internal differentiation or rapid volatile degassing, yields important implications for the final structure, composition and evolution of terrestrial planets. SLR-driven heating in the Solar System is sensitive to the absolute abundance and homogeneity of SLRs within the protoplanetary disk present during the condensation of the first solids. In order to explain the diverse compositions found for extrasolar planets, it is important to understand the distribution of SLRs in active planet formation regions (star clusters) during their first few Myr of evolution. By constraining the range of possible effects, we show how the imprint of SLRs can be extrapolated to exoplanetary systems and derive statistical predictions for the distribution of aluminum-26 and iron-60 based on N-body simulations of typical to large clusters (1000-10000 stars) with a range of initial conditions. We quantify the pollution of protoplanetary disks by supernova ejecta and show that the likelihood of enrichment levels similar to or higher than the Solar System can vary considerably, depending on the cluster morphology. Furthermore, many enriched systems show an excess in radiogenic heating compared to Solar System levels, which implies that the formation and evolution of planetesimals could vary significantly depending on the birth environment of their host stars.Comment: 15 pages, 8 figures, 4 tables; accepted for publication in MNRAS; associated video files can be found at http://timlichtenberg.net/2016_enrichment.htm