-
PRODIGE -- Planet-forming disks in Taurus with NOEMA
Authors:
R. Franceschi,
Th. Henning,
G. V. Smirnov-Pinchukov,
D. A. Semenov,
K. Schwarz,
A. Dutrey,
E. Chapillon,
U. Gorti,
S. Guilloteau,
V. Piétu,
S. van Terwisga,
L. Bouscasse,
P. Caselli,
G. Gieser,
T. -H. Hsieh,
A. Lopez-Sepulcre,
D. M. Segura-Cox,
J. E. Pineda,
M. J. Maureira,
M. T. Valdivia-Mena
Abstract:
We aim to constrain the gas density and temperature distributions as well as gas masses in several T Tauri protoplanetary disks located in Taurus. We use the 12CO, 13CO, and C18O (2-1) isotopologue emission observed at 0.9 with the IRAM NOrthern Extended Millimeter Array (NOEMA) as part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution PIs: P. Caselli & Th. Hennin…
▽ More
We aim to constrain the gas density and temperature distributions as well as gas masses in several T Tauri protoplanetary disks located in Taurus. We use the 12CO, 13CO, and C18O (2-1) isotopologue emission observed at 0.9 with the IRAM NOrthern Extended Millimeter Array (NOEMA) as part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution PIs: P. Caselli & Th. Henning). Our sample consists of Class II disks with no evidence of strong radial substructures. We use thesedata to constrain the thermal and chemical structure of these disks through theoretical models for gas emission. To fit the combined optically thick and thin CO line data in Fourier space, we developed the DiskCheF code, which includes the parameterized disk physical structure, machine-learning (ML) accelerated chemistry, and the RADMC-3D line radiative transfer module. A key novelty of DiskCheF is the fast and feasible ML-based chemistry trained on the extended grid of the disk physical-chemical models precomputed with the ANDES2 code. This ML approach allows complex chemical kinetics models to be included in a time-consuming disk fitting without the need to run a chemical code. We present a novel approach to incorporate chemistry into disk modeling without the need to explicitly calculate a chemical network every time. Using this new disk modeling tool, we successfully fit the 12CO, 13CO, and C18O (2-1) data from the CI, CY, DL, DM, DN, and IQ Tau disks. The combination of optically thin and optically thick CO lines allows us to simultaneously constrain the disk temperature and mass distribution, and derive the CO-based gas masses. These values are in reasonable agreement with the disk dust masses rescaled by a factor of 100 as well as with other indirect gas measurements.
△ Less
Submitted 24 June, 2024;
originally announced June 2024.
-
Probing the Physics of Star-Formation (ProPStar): II. The first systematic search for streamers toward protostars
Authors:
María Teresa Valdivia-Mena,
Jaime E. Pineda,
Paola Caselli,
Dominique M. Segura-Cox,
Anika Schmiedeke,
Silvia Spezzano,
Stella Offner,
Alexei V. Ivlev,
Michael Küffmeier,
Nichol Cunningham,
Roberto Neri,
María José Maureira
Abstract:
The detection of narrow channels of accretion toward protostellar disks, known as streamers, have increased in number in the last few years. However, it is unclear if streamers are a common feature around protostars that were previously missed, or if they are a rare phenomenon. Our goals are to obtain the incidence of streamers toward a region of clustered star formation and to trace the origins o…
▽ More
The detection of narrow channels of accretion toward protostellar disks, known as streamers, have increased in number in the last few years. However, it is unclear if streamers are a common feature around protostars that were previously missed, or if they are a rare phenomenon. 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 if they originate within the filamentary structure of molecular clouds or from beyond. 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 IRAS 4 and SVS 13 systems. 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. 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 when looking at a clustered star forming region. The incidence increases to about 60\% when we considered only embedded protostars. All streamers are found in HC3N emission. Given the different velocities between HC3N and N2H+ emission, and the fact that, 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 to the filaments with protostellar disk scales.
△ Less
Submitted 2 April, 2024;
originally announced April 2024.
-
PRODIGE -- Envelope to Disk with NOEMA III. The origin of complex organic molecule emission in SVS13A
Authors:
T. -H. Hsieh,
J. E. Pineda,
D. M. Segura-Cox,
P. Caselli,
M. T. Valdivia-Mena,
C. Gieser,
M. J. Maureira,
A. Lopez-Sepulcre,
L. Bouscasse,
R. Neri,
Th. Möller,
A. Dutrey,
A. Fuente,
D. Semenov,
E. Chapillon,
N. Cunningham,
Th. Henning,
V. Pietu,
I. Jimenez-Serra,
S. Marino,
C. Ceccarelli
Abstract:
Complex Organic Molecules (COMs) have been found toward low-mass protostars but the origins of the COM emission are still unclear. It can be associated with, for example, hot corinos, outflows, and/or accretion shock/disk atmosphere. We have conducted NOEMA observations toward SVS13A from the PROtostars & DIsks: Global Evolution (PRODIGE) program. Our previous \ce{DCN} observations reveal a possib…
▽ More
Complex Organic Molecules (COMs) have been found toward low-mass protostars but the origins of the COM emission are still unclear. It can be associated with, for example, hot corinos, outflows, and/or accretion shock/disk atmosphere. We have conducted NOEMA observations toward SVS13A from the PROtostars & DIsks: Global Evolution (PRODIGE) program. Our previous \ce{DCN} observations reveal a possible infalling streamer, which may affect the chemistry of the central protobinary by inducing accretion outbursts and/or shocked gas. Here, we further analyze six O-bearing COMs: CH3OH, aGg'-(CH2OH)2, C2H5OH, CH2(OH)CHO, CH3CHO, and CH3OCHO. Although the COM emission is not spatially resolved, we constrain the source sizes to $\lesssim0.3-0.4$ arcsec (90$-$120 au) by conducting uv-domain Gaussian fitting. Interestingly, the high-spectral resolution data reveal complex line profiles with multiple peaks showing differences between these six O-bearing COMs. The LTE fitting unveils differences in excitation temperatures and emitting areas among these COMs. We further conduct multiple-velocity-component LTE fitting to decompose the line emission into different kinematic components. Up to 6 velocity components are found from the LTE modeling. The temperature, column density, and source size of these components from each COM are obtained. We find a variety in excitation temperatures ($100-500$ K) and source sizes (D$\sim10-70$ au) from these kinematic components from different COMs. The emission of each COM can trace several components and different COMs most likely trace different regions. Given this complex structure, we suggest that the central region is inhomogeneous and unlikely to be heated by only protostellar radiation. We conclude that accretion shocks induced by the large-scale infalling streamer likely exist and contribute to the complexity of the COM emission.
△ Less
Submitted 25 March, 2024;
originally announced March 2024.
-
Probing the physics of star formation (ProPStar): I. First resolved maps of the electron fraction and cosmic-ray ionization rate in NGC 1333
Authors:
Jaime E. Pineda,
Olli Sipilä,
Dominique M. Segura-Cox,
Maria Teresa Valdivia-Mena,
Roberto Neri,
Michael Kuffmeier,
Alexei V. Ivlev,
Stella S. R. Offner,
Maria Jose Maureira,
Paola Caselli,
Silvia Spezzano,
Nichol Cunningham,
Anika Schmiedeke,
Mike Chen
Abstract:
Electron fraction and cosmic-ray ionization rates (CRIR) 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. We aim to derive the electron fraction and CRIR maps of an active star-forming re…
▽ More
Electron fraction and cosmic-ray ionization rates (CRIR) 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. We aim to derive the electron fraction and CRIR maps of an active star-forming region. 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 H$^{13}$CO$^+$ ratio (in addition to complementary data) to estimate the electron fraction and produce cosmic-ray ionization rate maps. We derived the first large-area electron fraction and CRIR resolved maps in a star-forming region, with typical values of $10^{-6.5}$ 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 MHD turbulence dissipation scale of 0.054 pc, which could be tested with future observations. We found a higher cosmic-ray ionization rate compared to the canonical value for $N({\rm H_2})=10^{21}-10^{23}$ 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 $ζ({\rm H_2})$, due to YSOs, should be taken into account in the analysis of clustered star formation.
△ Less
Submitted 25 February, 2024;
originally announced February 2024.
-
PRODIGE -- Planet-forming disks in Taurus with NOEMA. I. Overview and first results for 12CO, 13CO, and C18O
Authors:
D. Semenov,
Th. Henning,
S. Guilloteau,
G. Smirnov-Pinchukov,
A. Dutrey,
E. Chapillon,
V. Pietu,
R. Franceschi,
K. Schwarz,
S. van Terwisga,
L. Bouscasse,
P. Caselli,
C. Ceccarelli,
N. Cunningham,
A. Fuente,
C. Gieser,
T. -H. Hsieh,
A. Lopez-Sepulcre,
D. M. Segura-Cox,
J. E. Pineda,
M. J. Maureira,
Th. Moeller,
M. Tafalla,
M. T. Valdivia-Mena
Abstract:
We are performing a line survey of 8 planet-forming Class II disks in Taurus with the IRAM NOrthern Extended Millimeter Array (NOEMA), as a part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution; PIs: P. Caselli and Th. Henning). Compact and extended disks around T Tauri stars CI, CY, DG, DL, DM, DN, IQ Tau, and UZ Tau E are observed in ~80 lines from >20 C-, O,-…
▽ More
We are performing a line survey of 8 planet-forming Class II disks in Taurus with the IRAM NOrthern Extended Millimeter Array (NOEMA), as a part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution; PIs: P. Caselli and Th. Henning). Compact and extended disks around T Tauri stars CI, CY, DG, DL, DM, DN, IQ Tau, and UZ Tau E are observed in ~80 lines from >20 C-, O,- N-, and S-bearing species. The observations in four spectral settings at 210-280 GHz with $1σ$ rms sensitivity of ~ 8-12 mJy/beam at 0.9" and 0.3 km/s resolution will be completed in 2024. The uv-visibilities are fitted with the DiskFit model to obtain key stellar and disk properties. In this paper, the combined $^{12}$CO, $^{13}$CO and C$^{18}$O $J = 2-1$ data are presented. We find that the CO fluxes and disk masses inferred from dust continuum tentatively correlate with the CO emission sizes. We constrain dynamical stellar masses, geometries, temperatures, the CO column densities and gas masses for each disk. The best-fit temperatures at 100 au are ~ 17-37 K, and decrease radially with the power-law exponent q ~ 0.05-0.76. The inferred CO column densities decrease radially with the power-law exponent p ~ 0.2-3.1. The gas masses estimated from $^{13}$CO (2-1) are ~ $0.001-0.2 M_\textrm{Sun}$. The best-fit CO column densities point to severe CO freeze-out in the disks. The DL Tau disk is an outlier, and has either stronger CO depletion or lower gas mass than the rest of the sample. The CO isotopologue ratios are roughly consistent with the observed values in disks and the low-mass star-forming regions.
△ Less
Submitted 27 February, 2024; v1 submitted 22 February, 2024;
originally announced February 2024.
-
Infall of material onto the filaments in Barnard 5
Authors:
Spandan Choudhury,
Jaime E. Pineda,
Paola Caselli,
Michael Chun-Yuan Chen,
Stella S. R. Offner,
Maria Teresa Valdivia-Mena
Abstract:
Aims. We aim to study the structure and kinematics of the two filaments inside the subsonic core Barnard 5 in Perseus using high-resolution ($\approx$ 2400 au) NH3 data and a multi-component fit analysis.
Methods. We used observations of NH3 (1,1) and (2,2) inversion transitions using the Very Large Array (VLA) and the Green Bank Telescope (GBT). We smoothed the data to a beam of 8'' to reliably…
▽ More
Aims. We aim to study the structure and kinematics of the two filaments inside the subsonic core Barnard 5 in Perseus using high-resolution ($\approx$ 2400 au) NH3 data and a multi-component fit analysis.
Methods. We used observations of NH3 (1,1) and (2,2) inversion transitions using the Very Large Array (VLA) and the Green Bank Telescope (GBT). We smoothed the data to a beam of 8'' to reliably fit multiple velocity components towards the two filamentary structures identified in B5.
Results. Along with the core and cloud components, which dominate the flux in the line of sight, we detected two components towards the two filaments showing signs of infall. We also detected two additional components that can possibly trace new material falling into the subsonic core of B5.
Conclusions. Following comparison with previous simulations of filament formation scenarios in planar geometry, we conclude that either the formation of the B5 filaments is likely to be rather cylindrically symmetrical or the filaments are magnetically supported. We also estimate infall rates of $1.6\times10^{-4}\,M_\odot\,yr^{-1}$ and $1.8\times10^{-4}\,M_\odot\,yr^{-1}$ (upper limits) for the material being accreted onto the two filaments. At these rates, the filament masses can change significantly during the core lifetime. We also estimate an upper limit of $3.5\times10^{-5}\,M_\odot\,yr^{-1}$ for the rate of possible infall onto the core itself. Accretion of new material onto cores indicates the need for a significant update to current core evolution models, where cores are assumed to evolve in isolation.
△ Less
Submitted 28 December, 2023;
originally announced December 2023.
-
Flow of gas detected from beyond the filaments to protostellar scales in Barnard 5
Authors:
M. T. Valdivia-Mena,
J. E. Pineda,
D. M. Segura-Cox,
P. Caselli,
A. Schmiedeke,
S. Choudhury,
S. S. R. Offner,
R. Neri,
A. Goodman,
G. A. Fuller
Abstract:
The infall of gas from outside natal cores has proven to feed protostars after the main accretion phase (Class 0). This changes our view of star formation to a picture that includes asymmetric accretion (streamers), and a larger role of the environment. However, the connection between streamers and the filaments that prevail in star-forming regions is unknown. We investigate the flow of material t…
▽ More
The infall of gas from outside natal cores has proven to feed protostars after the main accretion phase (Class 0). This changes our view of star formation to a picture that includes asymmetric accretion (streamers), and a larger role of the environment. However, the connection between streamers and the filaments that prevail in star-forming regions is unknown. We investigate the flow of material toward the filaments within Barnard 5 (B5) and the infall from the envelope to the protostellar disk of the embedded protostar B5-IRS1. Our goal is to follow the flow of material from the larger, dense core scale, to the protostellar disk scale. We present new HC$_3$N line data from the NOEMA and 30m telescopes covering the coherence zone of B5, together with ALMA H$_2$CO and C$^{18}$O maps toward the protostellar envelope. We fit multiple Gaussian components to the lines so as to decompose their individual physical components. We investigate the HC$_3$N velocity gradients to determine the direction of chemically-fresh gas flow. At envelope scales, we use a clustering algorithm to disentangle the different kinematic components within H$_2$CO emission. At dense core scales, HC$_3$N traces the infall from the B5 region toward the filaments. HC$_3$N velocity gradients are consistent with accretion toward the filament spines plus flow along them. We found a $\sim2800$ au streamer in H$_2$CO emission which is blueshifted with respect to the protostar and deposits gas at outer disk scales. The strongest velocity gradients at large scales curve toward the position of the streamer at small scales, suggesting a connection between both flows. Our analysis suggests that the gas can flow from the dense core to the protostar. This implies that the mass available for a protostar is not limited to its envelope, and can receiving chemically-unprocessed gas after the main accretion phase.
△ Less
Submitted 1 August, 2023; v1 submitted 26 July, 2023;
originally announced July 2023.
-
PRODIGE -- Envelope to Disk with NOEMA II. Small-scale temperature structure and a streamer feeding the SVS13A protobinary using CH3CN and DCN
Authors:
T. -H. Hsieh,
D. M. Segura-Cox,
J. E. Pineda,
P. Caselli,
L. Bouscasse,
R. Neri,
A. Lopez-Sepulcre,
M. T. Valdivia-Mena,
M. J. Maureira,
Th. Henning,
G. V. Smirnov-Pinchukov,
D. Semenov,
Th. Möller,
N. Cunningham,
A. Fuente,
S. Marino,
A. Dutrey,
M. Tafalla,
E. Chapillon,
C. Ceccarelli,
B. Zhao
Abstract:
Aims. We present high sensitivity and high-spectral resolution NOEMA observations of the Class 0/I binary system SVS13A, composed of the low-mass protostars VLA4A and VLA4B with a separation of ~90 au. VLA4A is undergoing an accretion burst that enriches the chemistry of the surrounding gas. This gives us an excellent opportunity to probe the chemical and physical conditions as well as the accreti…
▽ More
Aims. We present high sensitivity and high-spectral resolution NOEMA observations of the Class 0/I binary system SVS13A, composed of the low-mass protostars VLA4A and VLA4B with a separation of ~90 au. VLA4A is undergoing an accretion burst that enriches the chemistry of the surrounding gas. This gives us an excellent opportunity to probe the chemical and physical conditions as well as the accretion process. Methods. We observe the (12K-11K) lines of CH3CN and CH313CN, the DCN (3-2) line, and the C18O (2-1) line toward SVS13A using NOEMA. Results. We find complex line profiles at disk scales which cannot be explained by a single component or pure Keplerian motion. By adopting two velocity components to model the complex line profiles, we find that the temperatures and densities are significantly different between these two components. This suggests that the physical conditions of the emitting gas traced via CH3CN can change dramatically within the circumbinary disk. In addition, combining our observations of DCN (3-2) with previous ALMA high-angular-resolution observations, we find that the binary system (or VLA4A) might be fed by an infalling streamer from envelope scales (~700 au). If this is the case, this streamer contributes to the accretion of material onto the system with a rate of at least 1.4x10-6 Msun yr-1. Conclusions. We conclude that the CH3CN emission in SVS13A traces hot gas from a complex structure. This complexity might be affected by a streamer that is possibly infalling and funneling material into the central region.
△ Less
Submitted 25 December, 2022; v1 submitted 9 November, 2022;
originally announced November 2022.
-
PRODIGE -- Envelope to disk with NOEMA I. A 3000 au streamer feeding a Class I protostar
Authors:
M. T. Valdivia-Mena,
J. E. Pineda,
D. M. Segura-Cox,
P. Caselli,
R. Neri,
A. López-Sepulcre,
N. Cunningham,
L. Bouscasse,
D. Semenov,
Th. Henning,
V. Piétu,
E. Chapillon,
A. Dutrey,
A. Fuente,
S. Guilloteau,
T. H. Hsieh,
I. Jiménez-Serra,
S. Marino,
M. J. Maureira,
G. V. Smirnov-Pinchukov,
M. Tafalla,
B. Zhao
Abstract:
Context. In the past few years, there has been a rise in the detection of streamers, asymmetric flows of material directed toward the protostellar disk with material from outside the star's natal core. It is unclear how they affect the process of mass accretion, in particular beyond the Class 0 phase. Aims. We investigate the gas kinematics around Per-emb-50, a Class I source in the crowded star-f…
▽ More
Context. In the past few years, there has been a rise in the detection of streamers, asymmetric flows of material directed toward the protostellar disk with material from outside the star's natal core. It is unclear how they affect the process of mass accretion, in particular beyond the Class 0 phase. Aims. We investigate the gas kinematics around Per-emb-50, a Class I source in the crowded star-forming region NGC 1333. Our goal is to study how the mass infall proceeds from envelope to disk scales in this source. Results. We discover a streamer delivering material toward Per-emb-50 in H$_2$CO and C$^{18}$O emission. The streamer's emission can be well described by the analytic solutions for an infalling parcel of gas along a streamline with conserved angular momentum, both in the image plane and along the line of sight velocities. The streamer has a mean infall rate of $1.3 \times 10^{ -6}$ M$_{ \odot}$ yr$^{ -1}$, $5 -10$ times higher than the current accretion rate of the protostar. SO and SO$_2$ emission reveal asymmetric infall motions in the inner envelope, additional to the streamer around Per-emb-50. Furthermore, the presence of SO$_2$ could mark the impact zone of the infalling material. Conclusions. The streamer delivers sufficient mass to sustain the protostellar accretion rate and might produce an accretion burst, which would explain the protostar's high luminosity with respect to other Class I sources. Our results highlight the importance of late infall for protostellar evolution: streamers might provide a significant amount of mass for stellar accretion after the Class 0 phase.
△ Less
Submitted 1 August, 2022;
originally announced August 2022.
-
ALMA resolves molecular clouds in the metal poor Magellanic Bridge A
Authors:
M. T. Valdivia-Mena,
M. Rubio,
A. D. Bolatto,
H. P. Saldaño,
C. Verdugo
Abstract:
(Abridged)We characterize gas and dust emission in Magellanic Bridge A, which has the highest 870$μ$m excess of emission found in single dish surveys. Using the ALMA telescope, we mapped the Magellanic Bridge A molecular cloud with sub-parsec resolution, in 1.3 mm continuum and CO(2-1) line emission. We also map the cloud in 870$μ$m continuum and CO(2-1) line emission at ~6 pc resolution with APEX…
▽ More
(Abridged)We characterize gas and dust emission in Magellanic Bridge A, which has the highest 870$μ$m excess of emission found in single dish surveys. Using the ALMA telescope, we mapped the Magellanic Bridge A molecular cloud with sub-parsec resolution, in 1.3 mm continuum and CO(2-1) line emission. We also map the cloud in 870$μ$m continuum and CO(2-1) line emission at ~6 pc resolution with APEX. We combine the ALMA and APEX CO(2-1) line cubes to study the molecular gas emission. Magellanic Bridge A breaks up into two distinct molecular clouds in dust and CO(2-1) emission, which we call North and South. Dust emission in the North source, according to our best parameters from fitting the far-infrarred fluxes, is ~3 K colder than in the South source in correspondence to its less developed star formation. Both dust sources present large submillimeter excesses in LABOCA data: according to our best fits the excess over the modified blackbody (MBB) fit to the Spitzer/Herschel continuum are ~7 and ~3 for the North and South sources respectively. Nonetheless, we do not detect the corresponding 1.3 mm continuum with ALMA. Our limits are compatible with the extrapolation of the MBB fits and therefore we cannot independently confirm the excess at this longer wavelength. The CO(2-1) emission is in two parsec-sized clouds with virial masses around 400 and 700 Mo each. Their volume densities are ~700-2600 cm$^{-3}$, larger than typical bulk densities of Galactic molecular clouds. The CO-to-H2 conversion factor is 6.5 and 15 M$_{\odot}$ (K km s$^{-1}$ pc$^2$)$^{-1}$ for the North and South clouds, calculated using their respective virial masses and CO(2-1) luminosities. Gas mass estimates from our MBB fits to dust emission yields masses $M\sim1.3\times10^3$ M$_{\odot}$ and $2.9\times10^3$ M$_{\odot}$ for North and South respectively, a factor ~4 larger than the virial masses we infer from CO.
△ Less
Submitted 2 July, 2020;
originally announced July 2020.