Showing 1–50 of 229 results for author: Hennebelle, P

The evolution of the $M_{\mathrm{d}}-M_{\star}$ and $\dot M-M_{\star}$ correlations traces protoplanetary disc dispersal

Authors: Alice Somigliana, Leonardo Testi, Giovanni Rosotti, Claudia Toci, Giuseppe Lodato, Rossella Anania, Benoît Tabone, Marco Tazzari, Ralf Klessen, Ugo Lebreuilly, Patrick Hennebelle, Sergio Molinari

Abstract: (Abridged) Observational surveys of entire star-forming regions have provided evidence of power-law correlations between the disc properties and the stellar mass, especially the disc mass (${M_d \propto M_*}^{λらむだ_m}$) and the accretion rate ($\dot M \propto {M_*}^{λらむだ_{acc}}$). Whether the secular disc evolution affects said correlations is still debated: while the purely viscous scenario has been pro… ▽ More (Abridged) Observational surveys of entire star-forming regions have provided evidence of power-law correlations between the disc properties and the stellar mass, especially the disc mass (${M_d \propto M_*}^{λらむだ_m}$) and the accretion rate ($\dot M \propto {M_*}^{λらむだ_{acc}}$). Whether the secular disc evolution affects said correlations is still debated: while the purely viscous scenario has been probed, other mechanisms could impact differently. We study the evolution of the slopes $λらむだ_m$ and $λらむだ_{acc}$ in the wind-driven and hybrid case and compare it to the viscous prediction, using a combination of analytical calculations and numerical simulations (performed with the 1D population synthesis code Diskpop, that we also present and release). Assuming $M_d(0) \propto {M_*}^{λらむだ_{m, 0}}$ and $\dot M(0) \propto {M_*}^{λらむだ_{acc, 0}}$ as initial conditions, we find that viscous and hybrid accretion preserve the shape of the correlations and evolve their slope; on the other hand, MHD winds change the shape of the correlations, bending them according to the scaling of the accretion timescale with the stellar mass. We also show how a spread in the initial conditions conceals this behaviour. We then analyse the impact of disc dispersal, and find that the currently available sample sizes ($\sim 30$ discs at 5 Myr) introduce stochastic oscillations in the slopes evolution, which dominate over the physical signatures. Increasing the sample size could mitigate this issue: $\sim 140$ discs at 5 Myr, corresponding to the complete Upper Sco sample, would give small enough error bars to use the evolution of the slopes as a proxy for the driving mechanism of disc evolution. Finally, we discuss how the observational claim of steepening slopes necessarily leads to an initially steeper $M_d - M_*$ correlation with respect to $\dot M - M_*$. △ Less

Submitted 30 July, 2024; originally announced July 2024.

Comments: 14 pages, 5 figures and 3 tables (plus Appendix). Accepted for publication in A&A. Includes Diskpop and popcorn software release

ALMA-IMF XV: The core mass function in the high-mass star-formation regime

Authors: F. Louvet, P. Sanhueza, A. Stutz, A. Men'shchikov, F. Motte, R. Galván-Madrid, S. Bontemps, Y. Pouteau, A. Ginsburg, T. Csengeri, J. Di Francesco, P. Dell'Ova, M. González, P. Didelon, J. Braine, N. Cunningham, B. Thomasson, P. Lesaffre, P. Hennebelle, M. Bonfand, A. Gusdorf, R. H. Álverez-Gutiérrez, T. Nony, G. Busquet, F. Olguin , et al. (16 additional authors not shown)

Abstract: The stellar initial mass function (IMF) is critical to our understanding of star formation and the effects of young stars on their environment. On large scales, it enables us to use tracers such as UV or Halpha emission to estimate the star formation rate of a system and interpret unresolved star clusters across the universe. So far, there is little firm evidence of large-scale variations of the I… ▽ More The stellar initial mass function (IMF) is critical to our understanding of star formation and the effects of young stars on their environment. On large scales, it enables us to use tracers such as UV or Halpha emission to estimate the star formation rate of a system and interpret unresolved star clusters across the universe. So far, there is little firm evidence of large-scale variations of the IMF, which is thus generally considered universal. Stars form from cores and it is now possible to estimate core masses and compare the core mass function (CMF) with the IMF, which it presumably produces. The goal of the ALMA-IMF large program is to measure the core mass function at high linear resolution (2700 au) in 15 typical Milky Way protoclusters spanning a mass range of 2500 to 32700 Msun. In this work, we used two different core extraction algorithms to extract about 680 gravitationally bound cores from these 15 protoclusters. We adopt per core temperature using the temperature estimate from the PPMAP Bayesian method. A power-law fit to the CMF of the sub-sample of cores above the 1.64 Msun completeness limit, 330 cores, through the maximum likelihood estimate technique yields a slope of 1.97 +/- 0.06, significantly flatter than the 2.35 Salpeter slope. Assuming a self-similar mapping between the CMF and the IMF, this result implies that these 15 high-mass protoclusters will generate atypical IMFs. This sample is the largest to date produced and analysed self-consistently, derived at matched physical resolution, with per-core temperature estimates and cores as massive as 150 Msun. We provide the raw source extraction catalogues and the source derived size, temperature, mass, and spectral indices in the 15 protoclusters. △ Less

Submitted 26 July, 2024; originally announced July 2024.

Comments: Accepted by A&A

On the relation between magnetic field strength and gas density in the interstellar medium: A multiscale analysis

Authors: David J. Whitworth, Sundar Srinivasan, Ralph E. Pudritz, Mordecai M. Mac Low, Rowan J. Smith, Aina Palau, Kate Pattle, Gwendoline Eadie, Hector Robinson, Rachel Pillsworth, James Wadsley, Noe Brucy, Ugo Lebreuilly, Patrick Hennebelle, Philipp Girichidis, Fred A. Gent, Jessy Marin, Lylon Sánchez Valido, Vianey Camacho, Ralf S. Klessen, Enrique Vázquez-Semadeni

Abstract: The relation between magnetic field strength B and gas density n in the interstellar medium is of fundamental importance to many areas of astrophysics, from protostellar disks to galaxy evolution. We present and compare Bayesian analyses of the B - n relation for a comprehensive observational data set, as well as a large body of numerical MHD simulations. We extend the original Zeeman relation of… ▽ More The relation between magnetic field strength B and gas density n in the interstellar medium is of fundamental importance to many areas of astrophysics, from protostellar disks to galaxy evolution. We present and compare Bayesian analyses of the B - n relation for a comprehensive observational data set, as well as a large body of numerical MHD simulations. We extend the original Zeeman relation of Crutcher et al. (2010) with a large body of magnetic data that includes 700 observations with the Davis-Chandrasekhar-Fermi method. By using a new multiparameter Bayesian analysis we present a new, more general, time-averaged observational relation: B \propto n^{0.27 \pm 0.017} for n \leq n_0 and B \propto n^{0.54 \pm 0.18} for n \geq n_0, with n_0 = 924^(+145-144) cm^-3. We perform a separate analysis on 19 numerical magnetohydrodynamics simulations that cover a wide range of scales, resolutions, initial conditions, and completed with a variety of codes: arepo, flash, pencil, and ramses. The power law exponents derived from the simulations depend on several physical factors including: dynamo effects, time scales, turbulence, and the initial seed field strength. In particular, early-time simulations where the density, velocity and magnetic fields are unevolved do not match the observational scalings. The simulations that trace the observed density range best, the evolved dwarf galaxy and Milky Way like galaxy simulations, settle into a near consistent exponent of = 0.5 in the dense gas, with variability in the diffuse gas exponent. △ Less

Submitted 25 July, 2024; originally announced July 2024.

Comments: 23 figures, 28 pages, submitted to MNRAS, Comments welcome

Synthetic Modelling of Polarized Dust Emission in Intermediate-Mass YSOs: I: Constraining the Role of Iron Inclusions and Inelastic Relaxation on Grain Alignment with ALMA Polarization

Authors: Nguyen Chau Giang, V. J. M. Le Gouellec, Thiem Hoang, A. J. Maury, P. Hennebelle

Abstract: Iron inclusions embedded inside dust grains play a crucial role in both internal alignment via Barnett relaxation and external alignment via the MAgnetically Enhanced RAdiative Torque (MRAT) mechanism. Moreover, inelastic relaxation is predicted to dominate over Barnett relaxation in driving the internal alignment of micron-sized and very large grains above $10μみゅーm$ (VLGs). Yet, a detailed modeling… ▽ More Iron inclusions embedded inside dust grains play a crucial role in both internal alignment via Barnett relaxation and external alignment via the MAgnetically Enhanced RAdiative Torque (MRAT) mechanism. Moreover, inelastic relaxation is predicted to dominate over Barnett relaxation in driving the internal alignment of micron-sized and very large grains above $10μみゅーm$ (VLGs). Yet, a detailed modeling of polarized thermal dust emission from Class 0/I Young Stellar Objects (YSOs) taking into account these effects and their observational constraints are still lacking. In this paper, we perform synthetic dust polarization modeling for MHD simulations of an intermediate-mass YSOs using our updated POLARIS code developed by Giang et al. (2022). Synthetic polarization results are then post-processed with CASA to confront ALMA observations. We found that to reproduce the high polarization degree of $p \sim 5-30\%$ observed in protostellar envelopes by ALMA, micron-sized and VLGs must contain iron inclusions with $N_{\rm cl} \sim 5 - 10^{4}$ iron atoms per cluster, assuming $30\%$ of iron abundance locked inside dust grains under the cluster form. Inside the inner $\sim 500$ au region, it requires inelastic relaxation, larger iron inclusions of $N_{\rm cl} \sim 10^{2}-10^{4}$, and grain growth beyond $\geq 10μみゅーm$ to reproduce $\sim 3-10\%$ of polarization observed by ALMA. But given such a combination, the internal alignment and MRAT efficiency acting on VLGs still decrease toward the center, which produces the decrease of $p(\%)$ with increasing gas density, reaching $p \sim 1\%$ inside the disk. Additionally, our results show that the high polarization by ALMA toward YSOs is not totally an artifact caused by the interferometric filtering, but it could be explained by the efficient magnetic alignment of SPM grains by MRAT mechanism. △ Less

Submitted 14 July, 2024; originally announced July 2024.

Comments: 26 pages, submitted to ApJ

Description of turbulent dynamics in the interstellar medium: Multifractal microcanonical analysis: II. Sparse filtering of Herschel observation maps and visualization of filamentary structures at different length scales

Authors: A. Rashidi, H. Yahia, S. Bontemps, N. Schneider, L. Bonne, P. Hennebelle, J. Scholtys, G. Attuel, A. Turiel, R. Simon, A. Cailly, A. Zebadua, A. Cherif, C. Lacroix, M. Martin, A. El Aouni, C. Sakka, S. K. Maji

Abstract: We present significant improvements to our previous work on noise reduction in {\sl Herschel} observation maps by defining sparse filtering tools capable of handling, in a unified formalism, a significantly improved noise reduction as well as a deconvolution in order to reduce effects introduced by the limited instrumental response (beam). We implement greater flexibility by allowing a wider choic… ▽ More We present significant improvements to our previous work on noise reduction in {\sl Herschel} observation maps by defining sparse filtering tools capable of handling, in a unified formalism, a significantly improved noise reduction as well as a deconvolution in order to reduce effects introduced by the limited instrumental response (beam). We implement greater flexibility by allowing a wider choice of parsimonious priors in the noise-reduction process. More precisely, we introduce a sparse filtering and deconvolution approach approach of type $l^2$-$l^p$, with $p > 0$ variable and apply it to a larger set of molecular clouds using {\sl Herschel} 250 $μみゅー$m data in order to demonstrate their wide range of application. In the {\sl Herschel} data, we are able to use this approach to highlight extremely fine filamentary structures and obtain singularity spectra that tend to show a significantly less $\log$-normal behavior and a filamentary nature in the less dense regions. We also use high-resolution adaptive magneto-hydrodynamic simulation data to assess the quality of deconvolution in such a simulated beaming framework. △ Less

Submitted 5 June, 2024; v1 submitted 2 June, 2024; originally announced June 2024.

Comments: 24 pages, 32 figures, official acceptance date on 16/05/2024 (Astronomy & Astrophysics)

MSC Class: 85-02

Journal ref: Astronomy & Astrophysics 2024

Dust evolution during the protostellar collapse: influence on the coupling between the neutral gas and the magnetic field

Authors: Valentin Vallucci-Goy, Ugo Lebreuilly, Patrick Hennebelle

Abstract: The coupling between the magnetic field and the gas during the collapsing phase of star-forming cores is strongly affected by the dust size distribution, which is expected to evolve. We aim to investigate the influence of key parameters on the evolution of the dust distribution as well as on the magnetic resistivities during the protostellar collapse. We perform collapsing single zone simulations… ▽ More The coupling between the magnetic field and the gas during the collapsing phase of star-forming cores is strongly affected by the dust size distribution, which is expected to evolve. We aim to investigate the influence of key parameters on the evolution of the dust distribution as well as on the magnetic resistivities during the protostellar collapse. We perform collapsing single zone simulations with shark. The code computes the evolution of the dust distribution, accounting for different grain growth and destruction processes. It also computes the magnetic resistivities. We find that the dust distribution significantly evolves during the protostellar collapse, shaping the magnetic resistivities. The peak size of the distribution, the population of small grains and consequently the magnetic resistivities are controlled by both coagulation and fragmentation rates. Under standard assumptions, the small grains coagulate very early as they collide by ambipolar drift, yielding magnetic resistivities orders of magnitude away from the non-evolving dust case. In particular, the ambipolar resistivity ηいーた_AD is very high prior to nH=10^10 cm^-3, and as a consequence magnetic braking should be ineffective. In this case, large size protoplanetary disks should result, which is inconsistent with recent observations. To alleviate this tension, we identify mechanisms to reduce the ambipolar resistivity during the protostellar collapse. The most promising are namely: electrostatic repulsion and grain-grain erosion. The evolution of the magnetic resistivities during the protostellar collapse and consequently the shape of the magnetic field in the early life of the protoplanetary disk strongly depends on the possibility to repopulate the small grains or to prevent their early coagulation. Therefore, it is crucial to better constrain the collision outcomes and the dust grain elastic properties. △ Less

Submitted 11 June, 2024; v1 submitted 31 May, 2024; originally announced June 2024.

Comments: 20 pages, 7 figures, accepted in A&A

Inefficient star formation in high Mach number environments II. Numerical simulations and comparison with analytical models

Authors: Noé Brucy, Patrick Hennebelle, Tine Colman, Ralf S. Klessen, Corentin Le Yhuelic

Abstract: Predicting the star formation rate (SFR) in galaxies is crucial to understand their evolution and morphology. To do so requires a fine understanding of how dense structures of gas are created and collapse. In that, turbulence and gravity play a major role. Within the gravo-turbulent framework, we assume that turbulence shapes the ISM, creating density fluctuations that, if gravitationally unstable… ▽ More Predicting the star formation rate (SFR) in galaxies is crucial to understand their evolution and morphology. To do so requires a fine understanding of how dense structures of gas are created and collapse. In that, turbulence and gravity play a major role. Within the gravo-turbulent framework, we assume that turbulence shapes the ISM, creating density fluctuations that, if gravitationally unstable, will collapse and form stars. The goal of this work is to quantify how different regimes of turbulence, characterized by the strength and compressibility of the driving, shape the density field. We are interested in the outcome in terms of SFR and how it compares with existing analytical models for the SFR. We run a series of hydrodynamical simulations of turbulent gas. The simulations are first conducted without gravity, so that the density and velocity are shaped by the turbulence driving. Gravity is then switched on, and the SFR is measured and compared with analytical models. The physics included in these simulations is very close to the one assumed in the classical gravo-turbulent SFR analytical models, which makes the comparison straightforward. We found that the existing analytical models convincingly agree with simulations at low Mach number, but we measure a much lower SFR in the simulation with a high Mach number. We develop, in a companion paper, an updated physically-motivated SFR model that reproduces well the inefficient high Mach regime of the simulations. Our work demonstrates that accurate estimations of the turbulent-driven replenishment time of dense structures and the dense gas spatial distribution are necessary to correctly predict the SFR in the high Mach regime. The inefficient high-Mach regime is a possible explanation for the low SFR found in dense and turbulent environments such as the centers of our Milky Way and other galaxies. △ Less

Submitted 25 July, 2024; v1 submitted 26 April, 2024; originally announced April 2024.

Comments: 13 pages, 14 figures, accepted for publication in A&A

Inefficient star formation in high Mach number environments I. The turbulent support analytical model

Authors: Patrick Hennebelle, Noé Brucy, Tine Colman

Abstract: The star formation rate (SFR), the number of stars formed per unit of time, is a fundamental quantity in the evolution of the Universe. While turbulence is believed to play a crucial role in setting the SFR, the exact mechanism remains unclear. Turbulence promotes star formation by compressing the gas, but also slows it down by stabilizing the gas against gravity. Most widely-used analytical model… ▽ More The star formation rate (SFR), the number of stars formed per unit of time, is a fundamental quantity in the evolution of the Universe. While turbulence is believed to play a crucial role in setting the SFR, the exact mechanism remains unclear. Turbulence promotes star formation by compressing the gas, but also slows it down by stabilizing the gas against gravity. Most widely-used analytical models rely on questionable assumptions, including: $i)$ integrating over the density PDF, a one-point statistical description that ignores spatial correlation, $ii)$ selecting self-gravitating gas based on a density threshold that often ignores turbulent dispersion, $iii)$ assuming the freefall time as the timescale for estimating SFR without considering the need to rejuvenate the density PDF, $iv)$ assuming the density PDF to be lognormal. Improving upon the only existing model that incorporates the spatial correlation of the density field, we present a new analytical model. We calculate the time needed to rejuvenate density fluctuations of a given density and spatial scale, revealing that it is generally much longer than the freefall time, rendering the latter inappropriate for use. We make specific predictions regarding the role of the Mach number, $ M $, and the driving scale of turbulence divided by the mean Jeans length. At low to moderate Mach numbers, turbulence does not reduce and may even slightly promote star formation by broadening the PDF. However, at higher Mach numbers, most density fluctuations are stabilized by turbulent dispersion, leading to a steep drop in the SFR as the Mach number increases. △ Less

Submitted 25 July, 2024; v1 submitted 26 April, 2024; originally announced April 2024.

Comments: 16 pages, 7 figures, Accepted for publication in A&A

Formation of low mass protostars and their circumstellar disks

Authors: Adnan Ali Ahmad, Matthias González, Patrick Hennebelle, Benoît Commerçon

Abstract: The birth process of circumstellar disks remains poorly constrained due to observational and numerical challenges. Recent numerical works have shown that the small-scale physics, often wrapped into a sub-grid model, play a crucial role in disk formation and evolution. This calls for a combined approach in which both the protostar and circumstellar disk are studied in concert. We aim to elucidate t… ▽ More The birth process of circumstellar disks remains poorly constrained due to observational and numerical challenges. Recent numerical works have shown that the small-scale physics, often wrapped into a sub-grid model, play a crucial role in disk formation and evolution. This calls for a combined approach in which both the protostar and circumstellar disk are studied in concert. We aim to elucidate the small scale physics and constrain sub-grid parameters commonly chosen in the literature by resolving the star-disk interaction. We carry out a set of very high resolution 3D radiative-hydrodynamics simulations that self-consistently describe the collapse of a turbulent dense molecular cloud core to stellar densities. We study the birth of the protostar, the circumstellar disk, and its early evolution (< 6 yr after protostellar formation). Following the second gravitational collapse, the nascent protostar quickly reaches breakup velocity and sheds its surface material, thus forming a hot ($\sim 10^{3}$ K), dense, and highly flared circumstellar disk. The protostar is embedded within the disk, such that material can flow without crossing any shock fronts. The circumstellar disk mass quickly exceeds that of the protostar, and its kinematics are dominated by self-gravity. Accretion onto the disk is highly anisotropic, and accretion onto the protostar mainly occurs through material that slides on the disk surface. The polar mass flux is negligible in comparison. The radiative behavior also displays a strong anisotropy, as the polar accretion shock is shown to be supercritical whereas its equatorial counterpart is subcritical. We also find a remarkable convergence of our results with respect to initial conditions. These results reveal the structure and kinematics in the smallest spatial scales relevant to protostellar and circumstellar disk evolution. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A, 17 pages, 15 figures

Journal ref: A&A 687, A90 (2024)

The Physical Origin of the Stellar Initial Mass Function

Authors: Patrick Hennebelle, Michael Y. Grudić

Abstract: Stars are amongst the most fundamental structures of our Universe. They comprise most of the baryonic and luminous mass of galaxies, synthethise heavy elements, and injec\ t mass, momentum, and energy into the interstellar medium. They are also home to the planets. Since stellar properties are primarily decided by their mass, the so-called \ stellar initial mass function (IMF) is critical to the s… ▽ More Stars are amongst the most fundamental structures of our Universe. They comprise most of the baryonic and luminous mass of galaxies, synthethise heavy elements, and injec\ t mass, momentum, and energy into the interstellar medium. They are also home to the planets. Since stellar properties are primarily decided by their mass, the so-called \ stellar initial mass function (IMF) is critical to the structuring of our Universe. We review the various physical processes, and theories which have been put forward as well as the numerical simulations which have been carried out to explain the origin of the stellar initial mass function. Key messages from this review are: (1) Gravity and turbulence most likely determine the power-law, high-mass part of the IMF. (2) Depending of the Mach number and the density distribution, several regimes are possible, including $Γがんま_{IMF} \simeq 0$, -0.8, -1 or -1.3 where $d N / d \log M \propto M^{Γがんま_{IMF}}$. These regimes are likely universal, however the transition between these regimes is not. (3) Protostellar jets can play a regulating influence on the IMF by injecting momentum into collapsing clumps and unbinding gas. (4) The peak of the IMF may be a consequence of dust opacity and molecular hydrogen physics at the origin of the first hydrostatic core. This depends weakly on large scale environmental conditions such as radiation, magnetic field, turbulence or metallicity. This likely constitutes one of the reason of the relative universality of the IMF. △ Less

Submitted 14 April, 2024; v1 submitted 10 April, 2024; originally announced April 2024.

Comments: To appear in Annual Reviews of Astronomy and Astrophysics

Signatures of magnetic braking in Class 0 protostars ? Exploring the gas kinematics in magnetized models of low-mass star formation

Authors: N. Añez-Lopez, U. Lebreuilly, A. Maury, P. Hennebelle

Abstract: Only indirect evidence of the role of magnetic braking in regulating gravitational collapse and the formation of circumstellar disks was found from observational work, such as compact disk sizes and the launching of high-velocity collimated protostellar jets. More direct tests of the magnetic braking shaping the angular momentum (AM) of the gas in Class 0 protostars are crucially needed. In the pr… ▽ More Only indirect evidence of the role of magnetic braking in regulating gravitational collapse and the formation of circumstellar disks was found from observational work, such as compact disk sizes and the launching of high-velocity collimated protostellar jets. More direct tests of the magnetic braking shaping the angular momentum (AM) of the gas in Class 0 protostars are crucially needed. In the present work we have used non-ideal MHD models of protostellar collapse and synthetic observations of molecular gas spectral emission that we analyze to test whether possible kinematic signatures of the magnetic braking in the gas velocity field can be captured from maps of the molecular gas emission in protostellar envelopes. By comparing the 3D Specific AM of models with varying turbulent energy and magnetization, we show that, in the numerical models of protostellar evolution explored, the increase in magnetization and its consequences on the spatial redistribution of SAM modifies the shapes of the radial profiles of SAM. We show that widely used observational methods fail to quantitatively capture the magnitude of SAM of the gas in protostellar envelopes, and that no method allows to measure the differences in radial evolution of SAM due to different magnetization at all envelope radii. This is especially true in the more magnetized cases. However, our analysis suggests that the detection of symmetric patterns and organized velocity fields, in the moment-1 maps of the molecular line emission as well as monotonous radial profiles of the SAM showing a power-law decline, should be suggestive of a less magnetized scenario. Protostellar cores where efficient magnetic braking is at work are more likely to present a highly asymmetric velocity field, and more prone to show complex radial profiles of their specific angular momentum measured in the equatorial plane. △ Less

Submitted 17 April, 2024; v1 submitted 26 March, 2024; originally announced March 2024.

Comments: Accepted A&A, 19 pages, 10 figures

Testing kinematic distances under a realistic Galactic potential

Authors: Glen H. Hunter, Mattia C. Sormani, Jan P. Beckmann, Eugene Vasiliev, Simon C. O. Glover, Ralf S. Klessen, Juan D. Soler, Noé Brucy, Philipp Girichidis, Junia Göller, Loke Ohlin, Robin Tress, Sergio Molinari, Ortwin Gerhard, Milena Benedettini, Rowan Smith, Patrick Hennebelle, Leonardo Testi

Abstract: Obtaining reliable distance estimates to gas clouds within the Milky Way is challenging in the absence of certain tracers. The kinematic distance approach has been used as an alternative, derived from the assumption of circular trajectories around the Galactic centre. Consequently, significant errors are expected in regions where gas flow deviates from purely circular motions. We aim to quantify t… ▽ More Obtaining reliable distance estimates to gas clouds within the Milky Way is challenging in the absence of certain tracers. The kinematic distance approach has been used as an alternative, derived from the assumption of circular trajectories around the Galactic centre. Consequently, significant errors are expected in regions where gas flow deviates from purely circular motions. We aim to quantify the systematic errors that arise from the kinematic distance method in the presence of a Galactic potential that is non-axisymmetric. We investigate how these errors differ in certain regions of the Galaxy and how they relate to the underlying dynamics. We perform 2D isothermal hydrodynamical simulation of the gas disk with the moving-mesh code Arepo, adding the capability of using an external potential provided by the Agama library for galactic dynamics. We introduce a new analytic potential of the Milky Way, taking elements from existing models and adjusting parameters to match recent observational constraints. We find significant errors in the kinematic distance estimate for gas close to the Sun, along sight lines towards the Galactic centre and anti-centre, and significant deviations associated with the Galactic bar. Kinematic distance errors are low within the spiral arms as gas resides close to local potential minima and the resulting line-of-sight velocity is close to what is expected for an axisymmetric potential. Interarm regions exhibit large deviations at any given Galactic radius. This is caused by the gas being sped up or slowed down as it travels into or out of the spiral arm. We are able to define 'zones of avoidance' in the lv-diagram, where the kinematic distance method is particularly unreliable and should only be used with caution. We report a power law relation between the kinematic distance error and the deviation of the project line-of-sight velocity from circular motion. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: 23 pages, 22 figures, 6 tables, submitted to A&A

Magnetic field morphology and evolution in the Central Molecular Zone and its effect on gas dynamics

Authors: R. G. Tress, M. C. Sormani, P. Girichidis, S. C. O. Glover, R. S. Klessen, R. J. Smith, E. Sobacchi, L. Armillotta, A. T. Barnes, C. Battersby, K. R. J. Bogue, N. Brucy, L. Colzi, C. Federrath, P. García, A. Ginsburg, J. Göller, H P. Hatchfield, C. Henkel, P. Hennebelle, J. D. Henshaw, M. Hirschmann, Y. Hu, J. Kauffmann, J. M. D. Kruijssen , et al. (12 additional authors not shown)

Abstract: The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We f… ▽ More The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We find that: (1) The magnetic field is conveniently decomposed into a regular time-averaged component and an irregular turbulent component. The regular component aligns well with the velocity vectors of the gas everywhere, including within the bar lanes. (2) The field geometry transitions from parallel to the Galactic plane near $z=0$ to poloidal away from the plane. (3) The magneto-rotational instability (MRI) causes an in-plane inflow of matter from the CMZ gas ring towards the central few parsecs of $0.01-0.1$ M$_\odot$ yr$^{-1}$ that is absent in the unmagnetised simulations. However, the magnetic fields have no significant effect on the larger-scale bar-driven inflow that brings the gas from the Galactic disc into the CMZ. (4) A combination of bar inflow and MRI-driven turbulence can sustain a turbulent vertical velocity dispersion of $σしぐま_z \simeq 5$ km s$^{-1}$ on scales of $20$ pc in the CMZ ring. The MRI alone sustains a velocity dispersion of $σしぐま_z \simeq 3$ km s$^{-1}$. Both these numbers are lower than the observed velocity dispersion of gas in the CMZ, suggesting that other processes such as stellar feedback are necessary to explain the observations. (5) Dynamo action driven by differential rotation and the MRI amplifies the magnetic fields in the CMZ ring until they saturate at a value that scales with the average local density as $B \simeq 102 (n/10^3 cm^{-3})^{0.33}$ $μみゅー$G. Finally, we discuss the implications of our results within the observational context in the CMZ. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: submitted to A&A, comments welcome

Grain growth and its chemical impact in the first hydrostatic core phase

Authors: D. Navarro-Almaida, U. Lebreuilly, P. Hennebelle, A. Fuente, B. Commerçon, R. Le Gal, V. Wakelam, M. Gerin, P. Riviére-Marichalar, L. Beitia-Antero, Y. Ascasibar

Abstract: The first hydrostatic core (FHSC) phase is a brief stage in the protostellar evolution that is difficult to detect. Our goal is to characterize the chemical evolution of gas and dust during the formation of the FHSC. Moreover, we are interested in analyzing, for the first time with 3D magnetohydrodynamic (MHD) simulations, the role of grain growth in its chemistry. We postprocessed… ▽ More The first hydrostatic core (FHSC) phase is a brief stage in the protostellar evolution that is difficult to detect. Our goal is to characterize the chemical evolution of gas and dust during the formation of the FHSC. Moreover, we are interested in analyzing, for the first time with 3D magnetohydrodynamic (MHD) simulations, the role of grain growth in its chemistry. We postprocessed $2\times10^{5}$ tracer particles from a $\texttt{RAMSES}$ non-ideal MHD simulation using the codes $\texttt{NAUTILUS}$ and $\texttt{SHARK}$ to follow the chemistry and grain growth throughout the simulation. A great chemical inheritance is seen, as gas-phase abundances of most of the C, O, N, and S reservoirs in the hot corino at the end of the simulation match the ice-phase abundances from the prestellar phase. Additionally, interstellar complex organic molecules (iCOMs) such as methyl formate, acetaldehyde, and formamide are formed during the warm-up process. The typical grain size in the hot corino $(n_{\rm H}>10^{11}\ {\rm cm^{-3}})$ increases forty-fold during the last 30 kyr, with negligible effects on its chemical composition. At moderate densities $(10^{10}<n_{\rm H}<10^{11}\ {\rm cm^{-3}})$ and cool temperatures $15<T<50$ K, increasing grain sizes delay molecular depletion. Finally, at low densities $(n_{\rm H}\sim10^{7}\ {\rm cm^{-3}})$, grains do not grow significantly. We also compared our results with a two-step model that reproduces well the abundances of C and O reservoirs, but not the N and S reservoirs. We conclude that the chemical composition of the FHSC is heavily determined by that of the parent prestellar core, chemo-MHD computations are needed for an accurate prediction of the abundances of the main N and S elemental reservoirs, and that the impact of grain growth in moderately dense areas delaying depletion permits the use of abundance ratios as grain growth proxies. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 24 pages, 26 figures, accepted for publication in Astronomy and Astrophysics

Cloud properties across spatial scales in simulations of the interstellar medium

Authors: Tine Colman, Noé Brucy, Philipp Girichidis, Simon C. O Glover, Milena Benedettini, Juan D. Soler, Robin G. Tress, Alessio Traficante, Patrick Hennebelle, Ralf S. Klessen, Sergio Molinari, Marc-Antoine Miville-Deschênes

Abstract: Molecular clouds (MC) are structures of dense gas in the interstellar medium (ISM), that extend from ten to a few hundred parsecs and form the main gas reservoir available for star formation. Hydrodynamical simulations of varying complexity are a promising way to investigate MC evolution and their properties. However, each simulation typically has a limited range in resolution and different cloud… ▽ More Molecular clouds (MC) are structures of dense gas in the interstellar medium (ISM), that extend from ten to a few hundred parsecs and form the main gas reservoir available for star formation. Hydrodynamical simulations of varying complexity are a promising way to investigate MC evolution and their properties. However, each simulation typically has a limited range in resolution and different cloud extraction algorithms are used, which complicates the comparison between simulations. In this work, we aim to extract clouds from different simulations covering a wide range of spatial scales. We compare their properties, such as size, shape, mass, internal velocity dispersion and virial state. We apply the Hop cloud detection algorithm on (M)HD numerical simulations of stratified ISM boxes and isolated galactic disk simulations that were produced using Flash Ramses and Arepo We find that the extracted clouds are complex in shape ranging from round objects to complex filamentary networks in all setups. Despite the wide range of scales, resolution, and sub-grid physics, we observe surprisingly robust trends in the investigated metrics. The mass spectrum matches in the overlap between simulations without rescaling and with a high-mass slope of $\mathrm{d} N/\mathrm{d}\ln M\propto-1$ in accordance with theoretical predictions. The internal velocity dispersion scales with the size of the cloud as $σしぐま\propto R^{0.75}$ for large clouds ($R\gtrsim3\,\mathrm{pc}$). For small clouds we find larger sigma compared to the power-law scaling, as seen in observations, which is due to supernova-driven turbulence. Almost all clouds are gravitationally unbound with the virial parameter scaling as $αあるふぁ_\mathrm{vir}\propto M^{-0.4}$, which is slightly flatter compared to observed scaling, but in agreement given the large scatter. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 22 pages, 16 figures, proposed for acceptance in A&A

Journal ref: A&A 686, A155 (2024)

Accuracy of ALMA estimates of young disk radii and masses. Predicted observations from numerical simulations

Authors: Ngo-Duy Tung, Leonardo Testi, Ugo Lebreuilly, Patrick Hennebelle, Anaëlle Maury, Ralf S. Klessen, Luca Cacciapuoti, Matthias González, Giovanni Rosotti, Sergio Molinari

Abstract: Protoplanetary disks, which are the natural consequence of the gravitational collapse of the dense molecular cloud cores, host the formation of the known planetary systems in our universe. Substantial efforts have been dedicated to investigating the properties of these disks in the more mature Class II stage, either via numerical simulations of disk evolution from a limited range of initial condit… ▽ More Protoplanetary disks, which are the natural consequence of the gravitational collapse of the dense molecular cloud cores, host the formation of the known planetary systems in our universe. Substantial efforts have been dedicated to investigating the properties of these disks in the more mature Class II stage, either via numerical simulations of disk evolution from a limited range of initial conditions or observations of their dust continuum and line emission from specific molecular tracers. The results coming from these two standpoints have been used to draw comparisons. However, few studies have investigated the main limitations at work when measuring the embedded Class 0/I disk properties from observations, especially in a statistical fashion. In this study, we provide a first attempt to compare the accuracy of some critical disk parameters in Class 0/I systems, as derived on real ALMA observational data, with the corresponding physical parameters that can be directly defined by theoreticians and modellers in numerical simulations. The approach we follow here is to provide full post-processing of the numerical simulations and apply it to the synthetic observations the same techniques used by observers to derive the physical parameters. We performed 3D Monte Carlo radiative transfer and mock interferometric observations of the disk populations formed in a magnetohydrodynamic (MHD) simulation model of disk formation through the collapse of massive clumps with the tools \textsc{Radmc-3d} and \textsc{Casa}, respectively, to obtain their synthetic observations. With these observations, we re-employed the techniques commonly used in disk modelling from their continuum emissions to infer the properties that would most likely be obtained with real interferometers. We then demonstrated how these properties may vary with respect to the gas kinematics analyses and dust continuum modelling. △ Less

Submitted 29 January, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: Accepted for publication in Astronomy & Astrophysics, 32 pages, 28 figures

Protoplanetary disk size under non-ideal magnetohydrodynamics: A general formalism with inclined magnetic field

Authors: Yueh-Ning Lee, Barshan Ray, Pierre Marchand, Patrick Hennebelle

Abstract: Many mechanisms have been proposed to alleviate the magnetic catastrophe, which prevents the Keplerian disk from forming inside a collapsing magnetized core. Such propositions include inclined field and non-ideal magnetohydrodynamics effects, and have been supported with numerical experiments. Models have been formulated for typical disk sizes when a field threads the rotating disk, parallel to th… ▽ More Many mechanisms have been proposed to alleviate the magnetic catastrophe, which prevents the Keplerian disk from forming inside a collapsing magnetized core. Such propositions include inclined field and non-ideal magnetohydrodynamics effects, and have been supported with numerical experiments. Models have been formulated for typical disk sizes when a field threads the rotating disk, parallel to the rotation axis, while observations at the core scales do not seem to show evident correlation between the directions of angular momentum and the magnetic field. In the present study, we propose a new model that considers both vertical and horizontal fields and discuss their effects on the protoplanetary disk size. △ Less

Submitted 5 January, 2024; originally announced January 2024.

Comments: Accepted for publication in ApJ Letters

Stellar Feedback in the Star Formation-Gas Density Relation: Comparison between Simulations and Observations

Authors: Paolo Suin, Annie Zavagno, Tine Colman, Patrick Hennebelle, Antoine Verliat, Delphine Russeil

Abstract: Context. The impact of stellar feedback on the Kennicutt-Schmidt law (KS law), which relates star formation rate (SFR) to surface gas density, is a topic of ongoing debate. The interpretation of individual cloud observations is challenging due to the various processes at play simultaneously and inherent biases. Therefore, a numerical investigation is necessary to understand the role of stellar fee… ▽ More Context. The impact of stellar feedback on the Kennicutt-Schmidt law (KS law), which relates star formation rate (SFR) to surface gas density, is a topic of ongoing debate. The interpretation of individual cloud observations is challenging due to the various processes at play simultaneously and inherent biases. Therefore, a numerical investigation is necessary to understand the role of stellar feedback and identify observable signatures. Aims. We investigate the role of stellar feedback on the KS law, aiming to identify distinct signatures that can be observed and analysed. Methods. We analyse MHD numerical simulations of a $10^4\,M_{\odot}$ cloud evolving under different feedback prescriptions. The set of simulations contains four types of feedback: with only protostellar jets, with ionising radiation from massive stars $(>8\,M_{\odot})$, with both of them and without any stellar feedback. To compare these simulations with the existing observational results, we analyse their evolution by adopting the same techniques applied in observational studies. Then, we simulate how the same analyses would change if the data were affected by typical observational biases. Conclusions. The presence of stellar feedback strongly influences the KS relation and the star formation efficiency per free-fall time ($εいぷしろん_\mathrm{ff}$). Its impact is primarily governed by its influence on the cloud's structure. Although the $εいぷしろん_\mathrm{ff}$ measured in our clouds results to be higher than what is usually observed in real clouds, upon applying prescriptions to mimic observational biases we recover good agreement with the expected values. Therefore, we can infer that observations tend to underestimate the total SFR. Moreover, this likely indicates that the physics included in our simulations is sufficient to reproduce the basic mechanisms contributing to set $εいぷしろん_\mathrm{ff}$. △ Less

Submitted 30 November, 2023; originally announced November 2023.

Comments: 16 pages, 13 figures, accepted for publication to A&A

Protostellar chimney flues: are jets and outflows lifting submillimetre dust grains from discs into envelopes?

Authors: L. Cacciapuoti, L. Testi, L. Podio, C. Codella, A. J. Maury, M. De Simone, P. Hennebelle, U. Lebreuilly, R. S. Klessen, S. Molinari

Abstract: Low dust opacity spectral indices ($βべーた< 1$) measured in the inner envelopes of class 0/I young stellar objects (age $\sim 10^{4-5}$ yr) have been interpreted as the presence of (sub-)millimetre dust grains in these environments. The density conditions and the lifetimes of collapsing envelopes have proven unfavorable for the growth of solids up to millimetre sizes. As an alternative, magneto-hydrody… ▽ More Low dust opacity spectral indices ($βべーた< 1$) measured in the inner envelopes of class 0/I young stellar objects (age $\sim 10^{4-5}$ yr) have been interpreted as the presence of (sub-)millimetre dust grains in these environments. The density conditions and the lifetimes of collapsing envelopes have proven unfavorable for the growth of solids up to millimetre sizes. As an alternative, magneto-hydrodynamical simulations suggest that protostellar jets and outflows might lift grains from circumstellar discs and diffuse them in the envelope. We reframe available data for the CALYPSO sample of Class 0/I sources and show tentative evidence for an anti-correlation between the value of $βべーた_{1-3mm}$ measured in the inner envelope and the mass loss rate of their jets and outflows, supporting a connection between the two. We discuss the implications that dust transport from the disc to the inner envelope might have for several aspects of planet formation. Finally, we urge for more accurate measurements of both correlated quantities and extension of this work to larger samples, necessary to further test the transport scenario. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: Accepted for publication on ApJ, 12 pages, 3 figures, 3 tables

A dusty streamer infalling onto the disk of a class I protostar. ALMA dual-band constraints on grain properties and mass infall rate

Authors: L. Cacciapuoti, E. Macias, A. Gupta, L. Testi, A. Miotello, C. Espaillat, M. Kuffmeier, S. van Terwisga, J. Tobin, S. Grant, C. F. Manara, D. Segura-Cox, J. Wendeborn, R. S. Klessen, A. J. Maury, U. Lebreuilly, P. Hennebelle, S. Molinari

Abstract: Observations of interstellar material infalling onto star- and planet-forming systems have become increasingly common thanks to recent advancements in radio interferometry. These structures replenish disks with fresh material, have the potential to significantly alter their dynamics, trigger the formation of substructures, induce shocks, and modify their physical and chemical properties. In this s… ▽ More Observations of interstellar material infalling onto star- and planet-forming systems have become increasingly common thanks to recent advancements in radio interferometry. These structures replenish disks with fresh material, have the potential to significantly alter their dynamics, trigger the formation of substructures, induce shocks, and modify their physical and chemical properties. In this study, we combine new ALMA band 3 and archival band 6 observations to characterize the dust content and mass infall rate of a 4,000 au arc-like structure infalling onto M512, a class I young stellar object located in the Lynds 1641 region of the Orion A molecular cloud. We measure for the first time spectral index maps and derive a dust opacity index profile along a streamer, constraining grain properties and its dust mass. We measure a spectral index $αあるふぁ\sim$ 3.2 across the entire structure, and a dust opacity index $βべーた\sim$ 1.6. Given grain properties consistent with the measured $βべーた$, the structure can host up to 245 M$_{\oplus}$ of dust, being comparable or even exceeding the mass of the inner, unresolved 600 au, which contains the protoplanetary disk of M512. Such a massive streamer can strongly affect the evolution of the star- and planet-forming inner system. Assuming typical ISM dust-to-gas ratio of 1%, free-fall timescales (50 kyr) imply total mass infall rates up to 1.5 $\cdot$ 10$^{-6}$ M$_{\odot}$/yr. M512 has been classified as an outbursting source with multi-epoch photometry, thus representing an interesting case study to explore the possible connection between infalling streamers and accretion outbursts. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: Accepted for publication in A&A, 13 pages, 5 figures, 3 tables

Environmental transition: overview of actions to reduce the environmental footprint of astronomy

Authors: Lucie Leboulleux, Faustine Cantalloube, Marie-Alice Foujols, Martin Giard, Jérôme Guilet, Jürgen Knödlseder, Alexandre Santerne, Lilia Todorov, Didier Barret, Olivier Berne, Aurélien Crida, Patrick Hennebelle, Quentin Kral, Eric Lagadec, Fabien Malbet, Julien Milli, Mamadou N'Diaye, Françoise Roques

Abstract: To keep current global warming below 1.5°C compared with the pre-industrial era, measures must be taken as quickly as possible in all spheres of society. Astronomy must also make its contribution. In this proceeding, and during the workshop to which it refers, different levers of actions are discussed through various examples: individual efforts, laboratory-level actions, impact evaluation and mit… ▽ More To keep current global warming below 1.5°C compared with the pre-industrial era, measures must be taken as quickly as possible in all spheres of society. Astronomy must also make its contribution. In this proceeding, and during the workshop to which it refers, different levers of actions are discussed through various examples: individual efforts, laboratory-level actions, impact evaluation and mitigation in major projects, institutional level, and involvement through collectives. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: 5 pages, 0 figures, SF2A conference 2023

Report number: 2023sf2a.conf..531L

Journal ref: SF2A proceedings 2023

Synthetic populations of protoplanetary disks. Impact of magnetic fields and radiative transfer

Authors: U. Lebreuilly, P. Hennebelle, T. Colman, A. Maury, N. -D. Tung, L. Testi, R. Klessen, S. Molinari, B. Commerçon, M. González, E. Pacetti, A. Somigliana, G. Rosotti

Abstract: Protostellar disks are the product of angular momentum conservation during the protostellar collapse. Understanding their formation is crucial because they are the birthplace of planets and because their formation is tightly related to star formation. Unfortunately, the initial properties of Class 0 disks and their evolution are still poorly constrained observationally and theoretically. We aim to… ▽ More Protostellar disks are the product of angular momentum conservation during the protostellar collapse. Understanding their formation is crucial because they are the birthplace of planets and because their formation is tightly related to star formation. Unfortunately, the initial properties of Class 0 disks and their evolution are still poorly constrained observationally and theoretically. We aim to better understand the mechanisms that set the statistics of disk properties as well as to study their formation in massive protostellar clumps. We also want to provide the community with synthetic disk populations to better interpret young disk observations. We use the ramses code to model star and disk formation in massive protostellar clumps with MHD including the effect of ambipolar diffusion and RT including the stellar radiative feedback. Those simulations, resolved up to the astronomical unit scale, allow to investigate the formation of disk populations. Magnetic fields play a crucial role in disk formation. A weaker initial field leads to larger and massive disks and weakens the stellar radiative feedback by increasing fragmentation. We find that ambipolar diffusion impacts disk and star formation and leads to very different disk magnetic properties. The stellar radiative feedback also have a strong influence, increasing the temperature and reducing fragmentation. Comparing our disk populations with observations reveals that our models with a mass-to-flux ratio of 10 seems to better reproduce observed disk sizes. This also sheds light on a tension between models and observations for the disk masses. The clump properties and physical modeling impact disk populations significantly. The tension between observations and models for disk mass estimates is critical to solve with synthetic observations in future years, in particular for our comprehension of planet formation. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 20 pages, 15 figures, accepted for publication in Astronomy & Astrophysics

The birth and early evolution of a low mass protostar

Authors: Adnan Ali Ahmad, Matthias González, Patrick Hennebelle, Benoît Commerçon

Abstract: Understanding the collapse of dense molecular cloud cores to stellar densities and the subsequent evolution of the protostar is of importance to model the feedback effects such an object has on its surrounding environment, as well as describing the conditions with which it enters the stellar evolutionary track. This process is fundamentally multi-scale and necessitates the use of robust numerical… ▽ More Understanding the collapse of dense molecular cloud cores to stellar densities and the subsequent evolution of the protostar is of importance to model the feedback effects such an object has on its surrounding environment, as well as describing the conditions with which it enters the stellar evolutionary track. This process is fundamentally multi-scale and necessitates the use of robust numerical simulations. We aim to model the birth and early evolution of a low-mass protostar. We also seek to describe the interior structure of the protostar and the radiative behavior of its accretion shock front. We carry out a high resolution numerical simulation of the collapse of a gravitationally unstable $1$ $\mathrm{M_{\odot}}$ dense molecular cloud core to stellar densities using three-dimensional radiation hydrodynamics under the gray flux-limited diffusion approximation. We follow the initial isothermal phase, the first adiabatic contraction, the second gravitational collapse triggered by the dissociation of $\mathrm{H}_{2}$ molecules, and $\approx 247$ days of the subsequent main accretion phase. We find that the sub-critical radiative behavior of the protostar's shock front causes it to swell as it accretes matter. We also find that the protostar is turbulent from the moment of its inception despite its radiative stability. This turbulence causes significant entropy mixing inside the protostar, which regulates the swelling. Furthermore, we find that the protostar is not fully ionized at birth, but the relative amount of ionized material within it increases as it accretes matter from its surroundings. Finally, we report in the appendix the results of the first 3D calculations involving a frequency-dependent treatment of radiative transfer, which has not produced any major differences with its gray counterpart. △ Less

Submitted 6 October, 2023; v1 submitted 2 October, 2023; originally announced October 2023.

Comments: Accepted for publication in A&A, 23 pages, 27 figures

Journal ref: A&A 680, A23 (2023)

The population of young low-mass stars in Trumpler 14

Authors: Dominika Itrich, Leonardo Testi, Giacomo Beccari, Carlo F. Manara, Megan Reiter, Thomas Preibisch, Anna F. McLeod, Giovanni Rosotti, Ralf Klessen, Sergio Molinari, Patrick Hennebelle

Abstract: Massive star-forming regions are thought to be the most common birth environments in the Galaxy and the only birth places of very massive stars. Their presence in the stellar cluster alters the conditions within the cluster impacting at the same time the evolution of other cluster members. In principle, copious amounts of ultraviolet radiation produced by massive stars can remove material from out… ▽ More Massive star-forming regions are thought to be the most common birth environments in the Galaxy and the only birth places of very massive stars. Their presence in the stellar cluster alters the conditions within the cluster impacting at the same time the evolution of other cluster members. In principle, copious amounts of ultraviolet radiation produced by massive stars can remove material from outer parts of the protoplanetary disks around low- and intermediate-mass stars in the process of external photoevaporation, effectively reducing the planet-formation capabilities of those disks. Here, we present deep VLT/MUSE observations of low-mass stars in Trumpler 14, one of the most massive, young, and compact clusters in the Carina Nebula Complex. We provide spectral and stellar properties of 717 sources and based on the distribution of stellar ages derive the cluster age of $\sim$1~Myr. The majority of the stars in our sample have masses $\leqslant$1~$M_\odot$, what makes our spectroscopic catalogue the most deep to date in term of masses, and proves that detailed investigations of low-mass stars are possible in the massive but distant regions. Spectroscopic studies of low-mass members of the whole Carina Nebula Complex are missing. Our work provides an important step forward towards filling this gap and set the stage for follow-up investigation of accretion properties in Trumpler 14. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Comments: Accepted for publication in A&A, 27 pages, 28 figures

Journal ref: A&A 685, A100 (2024)

Influence of protostellar outflows on star and protoplanetary disk formation in a massive star-forming clump

Authors: U. Lebreuilly, P. Hennebelle, A. Maury, M. González, A. Traficante, R. Klessen, L. Testi, S. Molinari

Abstract: Context. Due to the presence of magnetic fields, protostellar jets/outflows are a natural consequence of accretion onto protostars. They are expected to play an important role for star and protoplanetary disk formation. Aims. We aim to determine the influence of outflows on star and protoplanetary disk formation in star forming clumps. Methods. Using RAMSES, we perform the first magnetohydrodynami… ▽ More Context. Due to the presence of magnetic fields, protostellar jets/outflows are a natural consequence of accretion onto protostars. They are expected to play an important role for star and protoplanetary disk formation. Aims. We aim to determine the influence of outflows on star and protoplanetary disk formation in star forming clumps. Methods. Using RAMSES, we perform the first magnetohydrodynamics calculation of massive star-forming clumps with ambipolar diffusion, radiative transfer including the radiative feedback of protostars and protostellar outflows while systematically resolving the disk scales. We compare it to a model without outflows. Results. We find that protostellar outflows have a significant impact on both star and disk formation. They provide significant additional kinetic energy to the clump, with typical velocities of a few 10 km/s, impact the clump and disk temperatures, reduce the accretion rate onto the protostars and enhance fragmentation in the filaments. We find that they promote a more numerous stellar population. They do not impact much the low mass end of the IMF, which is probably controlled by the mass of the first Larson core, however, that they have an influence on its peak and high-mass end. Conclusions. Protostellar outflows appear to have a significant influence on both star and disk formation and should therefore be included in realistic simulations of star-forming environments. △ Less

Submitted 6 December, 2023; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: 16 pages, 14 figures. Accepted for publication at A&A

A deep-learning approach to the 3D reconstruction of dust density and temperature in star-forming regions

Authors: Victor F. Ksoll, Stefan Reissl, Ralf S. Klessen, Ian W. Stephens, Rowan J. Smith, Juan D. Soler, Alessio Traficante, Leonardo Testi, Patrick Hennebelle, Sergio Molinari

Abstract: Aims: We introduce a new deep-learning approach for the reconstruction of 3D dust density and temperature distributions from multi-wavelength dust emission observations on the scale of individual star-forming cloud cores (<0.2pc). Methods: We construct a training data set by processing cloud cores from the Cloud Factory simulations with the POLARIS radiative transfer code to produce synthetic du… ▽ More Aims: We introduce a new deep-learning approach for the reconstruction of 3D dust density and temperature distributions from multi-wavelength dust emission observations on the scale of individual star-forming cloud cores (<0.2pc). Methods: We construct a training data set by processing cloud cores from the Cloud Factory simulations with the POLARIS radiative transfer code to produce synthetic dust emission observations at 23 wavelengths between 12 and 1300 $μみゅー$m. We simplify the task by reconstructing the cloud structure along individual lines of sight and train a conditional invertible neural network (cINN) for this purpose. The cINN belongs to the group of normalising flow methods and is able to predict full posterior distributions for the target dust properties. We test different cINN setups, ranging from a scenario that includes all 23 wavelengths down to a more realistically limited case with observations at only seven wavelengths. We evaluate the predictive performance of these models on synthetic test data. Results: We report an excellent reconstruction performance for the 23-wavelengths cINN model, achieving median absolute relative errors of about 1.8% in $\log(n/m^{-3})$ and 1% in $\log(T_{dust}/K)$, respectively. We identify trends towards overestimation at the low end of the density range and towards underestimation at the high end of both the density and temperature values, which may be related to a bias in the training data. Limiting our coverage to a combination of only seven wavelengths, we still find a satisfactory performance with average absolute relative errors of about 3.3% and 2.5% in $\log(n/m^{-3})$ and $\log(T_{dust}/K)$. Conclusions: This proof-of-concept study shows that the cINN-based approach for 3D reconstruction of dust density and temperature is very promising and even compatible with a more realistically constrained wavelength coverage. △ Less

Submitted 9 February, 2024; v1 submitted 18 August, 2023; originally announced August 2023.

Comments: 38 pages, 25 figures, accepted for publication in Astronomy & Astrophysics on 26.01.2024

A comparison of the Milky Way's recent star formation revealed by dust thermal emission and high-mass stars

Authors: J. D. Soler, E. Zari, D. Elia, S. Molinari, C. Mininni, E. Schisano, A. Traficante, R. S. Klessen, S. C. O. Glover, P. Hennebelle, T. Colman, N. Frankel, T. Wenger

Abstract: We present a comparison of the Milky Way's star formation rate (SFR) surface density ($Σしぐま_{\rm SFR}$) obtained with two independent state-of-the-art observational methods. The first method infers $Σしぐま_{\rm SFR}$ from observations of the dust thermal emission from interstellar dust grains in far-infrared wavelengths registered in the Herschel infrared Galactic Plane Survey (Hi-GAL). The second method… ▽ More We present a comparison of the Milky Way's star formation rate (SFR) surface density ($Σしぐま_{\rm SFR}$) obtained with two independent state-of-the-art observational methods. The first method infers $Σしぐま_{\rm SFR}$ from observations of the dust thermal emission from interstellar dust grains in far-infrared wavelengths registered in the Herschel infrared Galactic Plane Survey (Hi-GAL). The second method determines $Σしぐま_{\rm SFR}$ by modeling the current population of O-, B-, and A-type stars in a 6 kpc $\times$ 6 kpc area around the Sun. We find an agreement between the two methods within a factor of two for the mean SFRs and the SFR surface density profiles. Given the broad differences between the observational techniques and the independent assumptions in the methods for computing the SFRs, this agreement constitutes a significant advance in our understanding of the star formation of our Galaxy and implies that the local SFR has been roughly constant over the past 10\,Myr. △ Less

Submitted 18 September, 2023; v1 submitted 2 August, 2023; originally announced August 2023.

Comments: Accepted for publication at Astronomy & Astrophysics (11SEP2023). 8 pages, 7 figures

Journal ref: A&A 678, A95 (2023)

FAUST IX. Multi-band, multi-scale dust study of L1527 IRS. Evidence for dust properties variations within the envelope of a Class 0/I YSO

Authors: L. Cacciapuoti, E. Macias, A. J. Maury, C. J. Chandler, N. Sakai, Ł. Tychoniec, S. Viti, A. Natta, M. De Simone, A. Miotello, C. Codella, C. Ceccarelli, L. Podio, D. Fedele, D. Johnstone, Y. Shirley, B. J. Liu, E. Bianchi, Z. E. Zhang, J. Pineda, L. Loinard, F. Ménard, U. Lebreuilly, R. S. Klessen, P. Hennebelle , et al. (3 additional authors not shown)

Abstract: Early dust grain growth in protostellar envelopes infalling on young discs has been suggested in recent studies, supporting the hypothesis that dust particles start to agglomerate already during the Class 0/I phase of young stellar objects (YSOs). If this early evolution were confirmed, it would impact the usually assumed initial conditions of planet formation, where only particles with sizes… ▽ More Early dust grain growth in protostellar envelopes infalling on young discs has been suggested in recent studies, supporting the hypothesis that dust particles start to agglomerate already during the Class 0/I phase of young stellar objects (YSOs). If this early evolution were confirmed, it would impact the usually assumed initial conditions of planet formation, where only particles with sizes $\lesssim 0.25 μみゅー$m are usually considered for protostellar envelopes. We aim to determine the maximum grain size of the dust population in the envelope of the Class 0/I protostar L1527 IRS, located in the Taurus star-forming region (140 pc). We use Atacama Large millimetre/sub-millimetre Array (ALMA) and Atacama Compact Array (ACA) archival data and present new observations, in an effort to both enhance the signal-to-noise ratio of the faint extended continuum emission and properly account for the compact emission from the inner disc. Using observations performed in four wavelength bands and extending the spatial range of previous studies, we aim to place tight constraints on the spectral ($αあるふぁ$) and dust emissivity ($βべーた$) indices in the envelope of L1527 IRS. We find a rather flat $αあるふぁ\sim$ 3.0 profile in the range 50-2000 au. Accounting for the envelope temperature profile, we derive values for the dust emissivity index, 0.9 < $βべーた$ < 1.6, and reveal a tentative, positive outward gradient. This could be interpreted as a distribution of mainly ISM-like grains at 2000 au, gradually progressing to (sub-)millimetre-sized dust grains in the inner envelope, where at R=300 au, $βべーた$ = 1.1 +/- 0.1. Our study supports a variation of the dust properties in the envelope of L1527 IRS. We discuss how this can be the result of in-situ grain growth, dust differential collapse from the parent core, or upward transport of disc large grains. △ Less

Submitted 21 November, 2023; v1 submitted 5 June, 2023; originally announced June 2023.

Comments: Accepted for publication in A&A. Contains 18 pages, 21 figures, 5 tables Replacement on Nov 22 to change title number of FAUST series from "X" to "IX."

Journal ref: A&A 676, A4 (2023)

Large-scale turbulent driving regulates star formation in high-redshift gas-rich galaxies II: Influence of the magnetic field and the turbulent compressive fraction

Authors: Noé Brucy, Patrick Hennebelle, Tine Colman, Simon Iteanu

Abstract: The observed star formation rate (SFR) in galaxies is well below what it should be if gravitational collapse alone were at play. It has recently been shown that one candidate that might regulate star formation, the feedback from massive stars, is suitable only if the mean column density at the kiloparsec scale is lower than $\approx 20 M_\odot\cdot\mathrm{pc}^{-2}$. On the other hand, intense larg… ▽ More The observed star formation rate (SFR) in galaxies is well below what it should be if gravitational collapse alone were at play. It has recently been shown that one candidate that might regulate star formation, the feedback from massive stars, is suitable only if the mean column density at the kiloparsec scale is lower than $\approx 20 M_\odot\cdot\mathrm{pc}^{-2}$. On the other hand, intense large-scale turbulent driving might slow down star formation in high-density environments to values that are compatible with observations. In this work, we explore the effect of the nature and strength of the turbulent driving, as well as the effect of the magnetic field. We performed a large series of feedback-regulated numerical simulations of the interstellar medium in which bidimensional large-scale turbulent driving was also applied. We determined the driving intensity needed to reproduce the Schmidt-Kennicutt relation for several gas column densities, magnetization, and driving compressibility. We confirm that in the absence of turbulent forcing and even with a substantial magnetic field, the SFR is too high, particularly at a high column density, compared to the Schmidt-Kennicutt relation. We find that the SFR outcome strongly depends on the initial magnetic field and on the compressibility of the turbulent driving. As a consequence, a higher magnetic field in high column density environment may lower the energy necessary to sustain a turbulence that is sufficiently intense to regulate star formation. Stellar feedback does not seem to be sufficient to regulate star formation in gas-rich galaxies where large-scale turbulent driving may be needed. The sources of this large-scale turbulence as well as its characteristics, such as its intensity, compressibility, and anisotropy, need to be understood and quantified. △ Less

Submitted 6 June, 2023; v1 submitted 29 May, 2023; originally announced May 2023.

Comments: 16 pages, 14 figures, accepted for publication in A&A

Journal ref: A&A 675, A144 (2023)

A panoptic view of the Taurus molecular cloud I. The cloud dynamics revealed by gas emission and 3D dust

Authors: J. D. Soler, C. Zucker, J. E. G. Peek, M. Heyer, P. F. Goldsmith, S. C. O. Glover, S. Molinari, R. S. Klessen, P. Hennebelle, L. Testi, T. Colman, M. Benedettini, D. Elia, C. Mininni, S. Pezzuto, E. Schisano, A. Traficante

Abstract: We present a study of the three-dimensional (3D) distribution of interstellar dust derived from stellar extinction observations toward the Taurus molecular cloud (MC) and its relation with the neutral atomic hydrogen (HI) emission at 21 cm wavelength and the carbon monoxide $^{12}$CO and $^{13}$CO emission in the $J=1\rightarrow0$ transition. We used the histogram of oriented gradients (HOG) metho… ▽ More We present a study of the three-dimensional (3D) distribution of interstellar dust derived from stellar extinction observations toward the Taurus molecular cloud (MC) and its relation with the neutral atomic hydrogen (HI) emission at 21 cm wavelength and the carbon monoxide $^{12}$CO and $^{13}$CO emission in the $J=1\rightarrow0$ transition. We used the histogram of oriented gradients (HOG) method to match the morphology in a 3D reconstruction of the dust density (3D dust) and the distribution of the gas tracers' emission. The result of the HOG analysis is a map of the relationship between the distances and radial velocities. The HOG comparison between the 3D dust and the HI emission indicates a morphological match at the distance of Taurus but an anti-correlation between the dust density and the HI emission, which uncovers a significant amount of cold HI within the Taurus MC. The HOG between the 3D dust and $^{12}$CO reveals a pattern in radial velocities and distances that is consistent with converging motions of the gas in the Taurus MC, with the near side of the cloud moving at higher velocities and the far side moving at lower velocities. This convergence of flows is likely triggered by the large-scale gas compression caused by the interaction of the Local Bubble and the Per-Tau shell, with Taurus lying at the intersection of the two bubble surfaces. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: 18 pages, 21 figures. Accepted for publication in Astronomy & Astrophysics (22MAY2023)

Journal ref: A&A 675, A206 (2023)

Spectral classification of young stars using conditional invertible neural networks I. Introducing and validating the method

Authors: Da Eun Kang, Victor F. Ksoll, Dominika Itrich, Leonardo Testi, Ralf S. Klessen, Patrick Hennebelle, Sergio Molinari

Abstract: Aims. We introduce a new deep learning tool that estimates stellar parameters (such as effective temperature, surface gravity, and extinction) of young low-mass stars by coupling the Phoenix stellar atmosphere model with a conditional invertible neural network (cINN). Our networks allow us to infer the posterior distribution of each stellar parameter from the optical spectrum. Methods. We discus… ▽ More Aims. We introduce a new deep learning tool that estimates stellar parameters (such as effective temperature, surface gravity, and extinction) of young low-mass stars by coupling the Phoenix stellar atmosphere model with a conditional invertible neural network (cINN). Our networks allow us to infer the posterior distribution of each stellar parameter from the optical spectrum. Methods. We discuss cINNs trained on three different Phoenix grids: Settl, NextGen, and Dusty. We evaluate the performance of these cINNs on unlearned Phoenix synthetic spectra and on the spectra of 36 Class III template stars with well-characterised stellar parameters. Results. We confirm that the cINNs estimate the considered stellar parameters almost perfectly when tested on unlearned Phoenix synthetic spectra. Applying our networks to Class III stars, we find good agreement with deviations of at most 5--10 per cent. The cINNs perform slightly better for earlier-type stars than for later-type stars like late M-type stars, but we conclude that estimations of effective temperature and surface gravity are reliable for all spectral types within the network's training range. Conclusions. Our networks are time-efficient tools applicable to large amounts of observations. Among the three networks, we recommend using the cINN trained on the Settl library (Settl-Net), as it provides the best performance across the largest range of temperature and gravity. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 29 pages, 19 figures, Accepted for publication by Astronomy & Astrophysics on 10. April

Journal ref: A&A 674, A175 (2023)

Physical conditions for dust grain alignment in Class 0 protostellar cores II. The role of the radiation field in models aligning/disrupting dust grains

Authors: Valentin J. M. Le Gouellec, Anaëlle J. Maury, Charles L. H. Hull, Antoine Verliat, Patrick Hennebelle, Valeska Valdivia

Abstract: The polarized dust emission observed in Class 0 protostellar cores at high angular resolution with ALMA has raised several concerns about the grain alignment conditions in these regions. We aim to study the role of the radiation field on the grain alignment mechanisms occurring in the interior (<1000 au) of Class 0 protostars. We produce synthetic observations of the polarized dust emission from a… ▽ More The polarized dust emission observed in Class 0 protostellar cores at high angular resolution with ALMA has raised several concerns about the grain alignment conditions in these regions. We aim to study the role of the radiation field on the grain alignment mechanisms occurring in the interior (<1000 au) of Class 0 protostars. We produce synthetic observations of the polarized dust emission from a MHD model of protostellar formation, using the POLARIS dust radiative transfer tool, which includes dust alignment with Radiative Torques Alignment (RATs). We test how the polarized dust emission from the model core depends on the irradiation conditions in the protostellar envelope, by varying the radiation due to accretion luminosity propagating from the central protostellar embryo throughout the envelope. The level of grain alignment efficiency obtained in the radiative transfer models is then compared to (sub-) millimeter ALMA dust polarization observations of Class 0 protostars. Our radiative transfer calculations have a central irradiation that reproduces the protostellar luminosities typically observed towards low- to intermediate-mass protostars, as well as super-paramagnetic grains, and grains >10 micron, which are required to bring the dust grain alignment efficiencies of the synthetic observations up to observed levels. Our radiative transfer calculations show that irradiation plays an important role in the mechanisms that dictate the size range of aligned grains in Class 0 protostars. Regions of the envelope that are preferentially irradiated harbor strong polarized dust emission but can be affected by the rotational disruption of dust grains. Episodes of high luminosity could affect grain alignment and trigger grain disruption mechanisms. [abridged] △ Less

Submitted 21 March, 2023; originally announced March 2023.

Journal ref: A&A 675, A133 (2023)

On the propagation of Alfvén waves in the dusty interstellar medium. Can we neglect dust inertia in molecular clouds?

Authors: P. Hennebelle, U. Lebreuilly

Abstract: Alfvén waves are fundamental magnetized modes which play an important role in the dynamics of magnetized flows such as the interstellar medium (ISM). In weakly ionised medium, their propagation critically depends on the ionisation rate but also on the charge carriers which, depending on gas density can be ions, electrons or dust grains. The latter in particular are well known to have a drastic inf… ▽ More Alfvén waves are fundamental magnetized modes which play an important role in the dynamics of magnetized flows such as the interstellar medium (ISM). In weakly ionised medium, their propagation critically depends on the ionisation rate but also on the charge carriers which, depending on gas density can be ions, electrons or dust grains. The latter in particular are well known to have a drastic influence on the magnetic resistivities in the dense ISM such as collapsing dense cores. Yet, in most calculations, for numerical reasons, the grain inertia is usually neglected. We investigate analytically the propagation of Alfvén waves both in a single-size and multi-size grain medium such as the ISM and we obtain exact expressions giving wavenumbers as a function of wave frequencies. These expressions are then solved analytically or numerically taking into account or neglecting grain inertia. Whereas at large wavelengths, neglecting grain inertia is a very good approximation, the situation is rather different for wavelengths shorter than a critical value, which broadly scales as $1/n$, $n$ being the gas density. More precisely, when inertia is neglected the waves do not propagate at short wavelengths or, due to the Hall effect, develop for one circular polarisation only, a whistler mode such that ${\cal R} _e (ωおめが) \propto k^2$, whereas the other polarisation presents a zero group velocity, i.e. ${\cal R} _e (ωおめが) \propto k^0$. When grain inertia is accounted for, the propagation of the two polarisations tend to be more symmetrical and the whistler mode is only present at density higher than $\simeq 10^8$ cm$^{-3}$. At lower density it is replaced by a mode having ${\cal R} _e (ωおめが) \propto k^{\simeq 1.2}$. Abridged △ Less

Submitted 17 March, 2023; originally announced March 2023.

Comments: submitted in A&A, comments welcome

Journal ref: A&A 674, A149 (2023)

Universal gravity-driven isothermal turbulence cascade in disk galaxies

Authors: Jérémy Fensch, Frédéric Bournaud, Noé Brucy, Yohan Dubois, Patrick Hennebelle, Joakim Rosdahl

Abstract: While interstellar gas is known to be supersonically turbulent, the injection processes of this turbulence are still unclear. Many studies suggest a dominant role of gravitational instabilities. However, their effect on galaxy morphology and large-scale dynamics vary across cosmic times, in particular due to the evolution of the gas fraction of galaxies. In this paper, we propose numerical simulat… ▽ More While interstellar gas is known to be supersonically turbulent, the injection processes of this turbulence are still unclear. Many studies suggest a dominant role of gravitational instabilities. However, their effect on galaxy morphology and large-scale dynamics vary across cosmic times, in particular due to the evolution of the gas fraction of galaxies. In this paper, we propose numerical simulations to follow the isothermal turbulent cascade of purely gravitationally-driven turbulence from its injection scale down to 0.095 pc for a gas-poor spiral disk and a gas-rich clumpy disk. To this purpose, and to lift the memory-footprint technical lock of sufficiently resolving the interstellar medium of a galaxy, we developed an encapsulated zoom method that allows us to probe self-consistently the self-generated turbulence cascade over three orders of magnitude on spatial scales. We follow this cascade for 10 Myrs. We find that the turbulent cascade follows the same scaling laws in both setups. Namely, in both cases the turbulence is close to equipartition between its compressive and solenoidal modes, the velocity power spectrum follows the Burgers' scaling and the density power spectrum is rather shallow, with a power-law slope of -0.7. Last, gravitationally-bound substructures follow a mass distribution with a -1.8 slope, similar to that of CO clumps. These simulations thus suggest a universality of gravity-driven isothermal turbulent cascade in disk galaxies across cosmic time. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 12 pages, 13 figures. Accepted for publication in A&A

Journal ref: A&A 672, A193 (2023)

On the distribution of the Cold Neutral Medium in galaxy discs

Authors: Rowan J. Smith, Robin Tress, Juan D. Soler, Ralf S. Klessen, Simon C. O. Glover, Patrick Hennebelle, Sergio Molinari, Mordecai-Mark Mac Low, David Whitworth

Abstract: The Cold Neutral Medium (CNM) is an important part of the galactic gas cycle and a precondition for the formation of molecular and star forming gas, yet its distribution is still not fully understood. In this work we present extremely high resolution simulations of spiral galaxies with time-dependent chemistry such that we can track the formation of the CNM, its distribution within the galaxy, and… ▽ More The Cold Neutral Medium (CNM) is an important part of the galactic gas cycle and a precondition for the formation of molecular and star forming gas, yet its distribution is still not fully understood. In this work we present extremely high resolution simulations of spiral galaxies with time-dependent chemistry such that we can track the formation of the CNM, its distribution within the galaxy, and its correlation with star formation. We find no strong radial dependence between the CNM fraction and total HI due to the decreasing interstellar radiation field counterbalancing the decreasing gas column density at larger galactic radii.However, the CNM fraction does increase in spiral arms where the CNM distribution is clumpy, rather than continuous, overlapping more closely with H2. The CNM doesn't extend out radially as far as HI, and the vertical scale height is smaller in the outer galaxy compared to HI with no flaring. The CNM column density scales with total midplane pressure and disappears from the gas phase below values of PT/kB =1000 K/cm3. We find that the star formation rate density follows a similar scaling law with CNM column density to the total gas Kennicutt-Schmidt law. In the outer galaxy we produce realistic vertical velocity dispersions in the HI purely from galactic dynamics but our models do not predict CNM at the extremely large radii observed in HI absorption studies of the Milky Way. We suggest that extended spiral arms might produce isolated clumps of CNM at these radii. △ Less

Submitted 19 July, 2023; v1 submitted 19 January, 2023; originally announced January 2023.

Comments: 13 pages, 20 figures, accepted by MNRAS

Dynamics of dust grains in turbulent molecular clouds. Conditions for decoupling and limits of different numerical implementations

Authors: Benoît Commerçon, Ugo Lebreuilly, Daniel J. Price, Francesco Lovascio, Guillaume Laibe, Patrick Hennebelle

Abstract: Dust grain dynamics in molecular clouds is regulated by its interplay with supersonic turbulent gas motions. The conditions under which dust grains decouple from the dynamics of gas remain poorly constrained. We first aim to investigate the critical dust grain size for dynamical decoupling, using both analytical predictions and numerical experiments. Second, we aim to set the range of validity of… ▽ More Dust grain dynamics in molecular clouds is regulated by its interplay with supersonic turbulent gas motions. The conditions under which dust grains decouple from the dynamics of gas remain poorly constrained. We first aim to investigate the critical dust grain size for dynamical decoupling, using both analytical predictions and numerical experiments. Second, we aim to set the range of validity of two fundamentally different numerical implementations for the evolution of dust and gas mixtures in turbulent molecular clouds. We carried out a suite of numerical experiments using two different schemes. First, we used a monofluid formalism in the terminal velocity approximation (TVA) on a Eulerian grid. Second, we used a two-fluid scheme, in which the dust dynamics is handled with Lagrangian super-particles, and the gas dynamics on a Eulerian grid. The monofluid results are in good agreement with the theoretical critical size for decoupling. We report dust dynamics decoupling for Stokes number St>0.1, that is, dust grains of $s>4~μみゅー$m in size. We find that the TVA is well suited for grain sizes of 10 $μみゅー$m in molecular clouds, in particular in the densest regions. However, the maximum dust enrichment measured in the low-density material where St>1 is questionable. In the Lagrangian dust experiments, we show that the results are affected by the numerics for all dust grain sizes. At St<<1, the dust dynamics is largely affected by artificial trapping in the high-density regions, leading to spurious variations of the dust concentration. At St>1, the maximum dust enrichment is regulated by the grid resolution used for the gas dynamics. The results of previous similar numerical work should therefore be revisited with respect to the limitations we highlight in this study. Dust enrichment of submicron dust grains is unlikely to occur in the densest parts of molecular clouds. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: 22 pages, 16 figures, accepted for publication in A&A

Journal ref: A&A 671, A128 (2023)

ALMA-IMF VI -- Investigating the origin of stellar masses: Core mass function evolution in the W43-MM2&MM3 mini-starburst

Authors: Y. Pouteau, F. Motte, T. Nony, M. Gonzalez, I. Joncour, J. -F. Robitaille, G. Busquet, R. Galvan-Madrid, A. Gusdorf, P. Hennebelle, A. Ginsburg, T. Csengeri, P. Sanhueza, P. Dell'Ova, A. M. Stutz, A. P. M. Towner, N. Cunningham, F. Louvet, A. Men'shchikov, M. Fernandez-Lopez, N. Schneider, M. Armante, J. Bally, T. Baug, M. Bonfand , et al. (13 additional authors not shown)

Abstract: Among the most central open questions regarding the initial mass function (IMF) of stars is the impact of environment on the shape of the core mass function (CMF) and thus potentially on the IMF. The ALMA-IMF Large Program aims to investigate the variations in the core distributions with cloud characteristics, as diagnostic observables of the formation process and evolution of clouds. The present… ▽ More Among the most central open questions regarding the initial mass function (IMF) of stars is the impact of environment on the shape of the core mass function (CMF) and thus potentially on the IMF. The ALMA-IMF Large Program aims to investigate the variations in the core distributions with cloud characteristics, as diagnostic observables of the formation process and evolution of clouds. The present study focuses on the W43-MM2&MM3 mini-starburst, whose CMF has recently been found to be top-heavy with respect to the Salpeter slope. W43-MM2&MM3 harbors a rich cluster that contains a statistically significant number of cores, which was previously characterized in Paper III. We applied a multi-scale decomposition technique to the ALMA 1.3 mm and 3 mm continuum images to define six subregions. For each subregion we characterized the high column density probability distribution function, n-PDF, and the shape of the cloud gas using the 1.3 mm image. Using the core catalog, we investigate correlations between the CMF and cloud and core properties. We classify the subregions into different stages of evolution, from quiescent to burst to post-burst, based on the surface number density of cores, number of outflows, and UCHii presence. The high-mass end of the subregion CMFs varies from being close to the Salpeter slope (quiescent) to top-heavy (burst and post-burst). Moreover, the second tail of the n-PDF varies from steep, to flat like observed for the high mass star-forming clouds. We found that subregions with flat second n-PDF tails display top-heavy CMFs. The CMF may evolve from Salpeter to top-heavy throughout the star formation process from the quiescent to the burst phase. This scenario raises the question of if the CMF might revert again to Salpeter as the cloud approaches the end of its star formation stage, a hypothesis that remains to be tested. △ Less

Submitted 20 February, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: Accepted for publication in A&A (december, 5th 2022) after language editing; 25 pages, 14 figures

Journal ref: A&A 674, A76 (2023)

Synthetic dust polarization emission maps at 353 GHz for an observer placed inside a Local Bubble-like cavity

Authors: E. Maconi, J. D. Soler, S. Reissl, P. Girichidis, R. S. Klessen, P. Hennebelle, S. Molinari, L. Testi, R. J. Smith, M. C. Sormani, J. W. Teh, A. Traficante, .

Abstract: We present a study of synthetic observations of polarized dust emission at 353 GHz as seen by an observer within a cavity in the interstellar medium (ISM). The cavity is selected from a magnetohydrodynamic simulation of the local ISM with time-dependent chemistry, star formation, and stellar feedback in form of supernova explosions with physical properties comparable to the Local Bubble ones. We f… ▽ More We present a study of synthetic observations of polarized dust emission at 353 GHz as seen by an observer within a cavity in the interstellar medium (ISM). The cavity is selected from a magnetohydrodynamic simulation of the local ISM with time-dependent chemistry, star formation, and stellar feedback in form of supernova explosions with physical properties comparable to the Local Bubble ones. We find that the local density enhancement together with the coherent magnetic field in the cavity walls makes the selected candidate a translucent polarization filter to the emission coming from beyond its domains. This underlines the importance of studying the Local Bubble in further detail. The magnetic field lines inferred from synthetic dust polarization data are qualitatively in agreement with the all-sky maps of polarized emission at 353 GHz from the Planck satellite in the latitudes interval 15deg <= |b| <= 65deg. As our numerical simulation allows us to track the Galactic midplane only out to distances of 250 pc, we exclude the region |b|<15deg from our analysis. At large Galactic latitudes, our model exhibits a high degree of small-scale structures. On the contrary, the observed polarization pattern around the Galactic poles is relatively coherent and regular, and we argue that the global toroidal magnetic field of the Milky Way is important for explaining the data at |b| > 65deg. We show that from our synthetic polarization maps, it is difficult to distinguish between an open and a closed Galactic cap using the inferred magnetic field morphology alone. △ Less

Submitted 19 June, 2023; v1 submitted 13 December, 2022; originally announced December 2022.

The Star Formation Rate of the Milky Way as seen by Herschel

Authors: D. Elia, S. Molinari, E. Schisano, J. D. Soler, M. Merello, D. Russeil, M. Veneziani, A. Zavagno, A. Noriega-Crespo, L. Olmi, M. Benedettini, P. Hennebelle, R. S. Klessen, S. Leurini, R. Paladini, S. Pezzuto, A. Traficante, D. J. Eden, P. G. Martin, M. Sormani, A. Coletta, T. Colman, R. Plume, Y. Maruccia, C. Mininni , et al. (1 additional authors not shown)

Abstract: We present a new derivation of the Milky Way's current star formation rate (SFR) based on the data of the Hi-GAL Galactic plane survey. We estimate the distribution of the SFR across the Galactic plane from the star-forming clumps identified in the Hi-GAL survey and calculate the total SFR from the sum of their contributions. The estimate of the global SFR amounts to $2.0 \pm 0.7$~M$_{\odot}$~yr… ▽ More We present a new derivation of the Milky Way's current star formation rate (SFR) based on the data of the Hi-GAL Galactic plane survey. We estimate the distribution of the SFR across the Galactic plane from the star-forming clumps identified in the Hi-GAL survey and calculate the total SFR from the sum of their contributions. The estimate of the global SFR amounts to $2.0 \pm 0.7$~M$_{\odot}$~yr$^{-1}$, of which $1.7 \pm 0.6$~M$_{\odot}$~yr$^{-1}$ coming from clumps with reliable heliocentric distance assignment. This value is in general agreement with estimates found in the literature of last decades. The profile of SFR density averaged in Galactocentric rings is found to be qualitatively similar to others previously computed, with a peak corresponding to the Central Molecular Zone and another one around Galactocentric radius $R_\mathrm{gal} \sim 5$~kpc, followed by an exponential decrease as $\log(Σしぐま_\mathrm{SFR}/[\mathrm{M}_\odot~\mathrm{yr}^{-1}~\mathrm{kpc}^{-2}])=a\,R_\mathrm{gal}/[\mathrm{kpc}]+b $, with $a=-0.28 \pm 0.01$. In this regard, the fraction of SFR produced within and outside the Solar circle is 84\% and 16\%, respectively; the fraction corresponding to the far outer Galaxy ($R_\mathrm{gal} > 13.5$~kpc) is only 1\%. We also find that, for $R_\mathrm{gal}>3$~kpc, our data follow a power law as a function of density, similarly to the Kennicutt-Schmidt relation. Finally, we compare the distribution of the SFR density across the face-on Galactic plane and those of median parameters, such as temperature, luminosity/mass ratio and bolometric temperature, describing the evolutionary stage of Hi-GAL clumps. We found no clear correlation between the SFR and the clump evolutionary stage. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: To be published in ApJ

How magnetic field and stellar radiative feedback influences the collapse and the stellar mass spectrum of a massive star forming clump

Authors: Patrick Hennebelle, Ugo Lebreuilly, Tine Colman, Davide Elia, Gary Fuller, Silvia Leurini, Thomas Nony, Eugenio Schisano, Juan D. Soler, Alessio Traficante, Ralf S. Klessen, Sergio Molinari, Leonardo Testi

Abstract: In spite of decades of theoretical efforts, the physical origin of the stellar initial mass function (IMF) is still debated. We aim at understanding the influence of various physical processes such as radiative stellar feedback, magnetic field and non-ideal magneto-hydrodynamics on the IMF. We present a series of numerical simulations of collapsing 1000 M$_\odot$ clumps taking into account radiati… ▽ More In spite of decades of theoretical efforts, the physical origin of the stellar initial mass function (IMF) is still debated. We aim at understanding the influence of various physical processes such as radiative stellar feedback, magnetic field and non-ideal magneto-hydrodynamics on the IMF. We present a series of numerical simulations of collapsing 1000 M$_\odot$ clumps taking into account radiative feedback and magnetic field with spatial resolution down to 1 AUえーゆー. Both ideal and non-ideal MHD runs are performed and various radiative feedback efficiencies are considered. We also develop analytical models that we confront to the numerical results. The sum of the luminosities produced by the stars in the calculations is computed and it compares well with the bolometric luminosities reported in observations of massive star forming clumps. The temperatures, velocities and densities are also found to be in good agreement with recent observations. The stellar mass spectrum inferred for the simulations is, generally speaking, not strictly universal and in particular varies with magnetic intensity. It is also influenced by the choice of the radiative feedback efficiency. In all simulations, a sharp drop in the stellar distribution is found at about $M_{min} \simeq$ 0.1 M$_\odot$, which is likely a consequence of the adiabatic behaviour induced by dust opacities at high densities. As a consequence, when the combination of magnetic and thermal support is not too large, the mass distribution presents a peak located at 0.3-0.5 M$_\odot$. When magnetic and thermal support are large, the mass distribution is better described by a plateau, i.e. $d N / d \log M \propto M^{-Γがんま}$, $Γがんま\simeq 0$. Abridged △ Less

Submitted 22 October, 2022; originally announced October 2022.

Comments: accepted for publication in A&A

Journal ref: A&A 668, A147 (2022)

3D chemical structure of the diffuse turbulent ISM II -- Origin of CH$^+$, new solution to an 80 years mystery

Authors: Benjamin Godard, Guillaume Pineau Des Forêts, Patrick Hennebelle, Elena Bellomi, Valeska Valdivia

Abstract: Aims: The large abundances of CH$^+$ in the diffuse interstellar medium (ISM) are a long standing issue of our understanding of the thermodynamical and chemical states of the gas. We investigate, here, the formation of CH+ in turbulent and multiphase environments, where the heating of the gas is almost solely driven by the photoelectric effect. Methods: The diffuse ISM is simulated using the magne… ▽ More Aims: The large abundances of CH$^+$ in the diffuse interstellar medium (ISM) are a long standing issue of our understanding of the thermodynamical and chemical states of the gas. We investigate, here, the formation of CH+ in turbulent and multiphase environments, where the heating of the gas is almost solely driven by the photoelectric effect. Methods: The diffuse ISM is simulated using the magnetohydrodynamic (MHD) code RAMSES which self-consistently computes the dynamical and thermal evolution of the gas along with the time-dependent evolutions of the abundances of H$^+$, H, and H$_2$. The rest of the chemistry, including the abundance of CH$^+$, is computed in post-processing, at equilibrium, under the constraint of out-ofequilibrium of H$^+$, H, and H$_2$. The comparison with the observations is performed taking into account an often neglected, yet paramount, piece of information, namely the length of the intercepted diffuse matter along the observed lines of sight. Results: The quasi totality of the mass of CH$^+$ originates from the unstable gas, in environments where the kinetic temperature is larger than 600 K, the density ranges between 0.6 and 10 cm$^{-3}$, the electronic fraction ranges between 3 x 10$^{-4}$ and 6 x 10$^{-3}$, and the molecular fraction is smaller than 0.4. Its formation is driven by warm and out-of-equilibrium H$_2$ initially formed in the cold neutral medium (CNM) and injected in more diffuse environments and even the warm neutral medium (WNM) through a combination of advection and thermal instability. The simulation which displays the tightest agreement with the HI-to-H$_2$ transition and the thermal pressure distribution observed in the Solar Neighborhood is found to naturally reproduce the observed abundances of CH$^+$, the dispersion of observations, the probability of occurrence of most of the lines of sight, the fraction of non-detections of CH$^+$, and the distribution of its line profiles. The amount of CH$^+$ and the statistical properties of the simulated lines of sight are set by the fraction of unstable gas rich in H$_2$ which is controlled, on Galactic scales, by the mean density of the diffuse ISM (or, equivalently, its total mass), the amplitude of the mean UV radiation field, and the strength of the turbulent forcing. Conclusions: This work offers a new and natural solution to an 80 years old chemical riddle. The almost ubiquitous presence of CH$^+$ in the diffuse ISM likely results from the exchanges of matter between the CNM and the WNM induced by the combination of turbulent advection and thermal instability, without the need to invoke ambipolar diffusion or regions of intermittent turbulent dissipation. Through two phase turbulent mixing, CH$^+$ might thus be a tracer of the H$_2$ mass loss rate of CNM clouds. △ Less

Submitted 21 September, 2022; originally announced September 2022.

Comments: Astronomy and Astrophysics, EDP Sciences, In press

Journal ref: A&A 669, A74 (2023)

Recent progress with observations and models to characterize the magnetic fields from star-forming cores to protostellar disks

Authors: Anaëlle Maury, Patrick Hennebelle, Josep Miquel Girart

Abstract: In this review paper, we aim at providing a global outlook on the progresses made in the recent years to characterize the role of magnetic fields during the embedded phases of the star formation process. Thanks to the development of observational capabilities and the parallel progress in numerical models capturing most of the important physics at work during star formation, it has recently become… ▽ More In this review paper, we aim at providing a global outlook on the progresses made in the recent years to characterize the role of magnetic fields during the embedded phases of the star formation process. Thanks to the development of observational capabilities and the parallel progress in numerical models capturing most of the important physics at work during star formation, it has recently become possible to confront detailed predictions of magnetized models to observational properties of the youngest protostars. We provide an overview of the most important consequences when adding magnetic fields to state-of-the-art models of protostellar formation, emphasizing their role to shape the resulting star(s) and their disk(s). We discuss the importance of magnetic field coupling to set the efficiency of magnetic processes, and provide a review of observational works putting constraints on the two main agents responsible for the coupling in star-forming cores: dust grains and ionized gas. We recall the physical processes and observational methods allowing to trace the magnetic field topology and its intensity in embedded protostars, and review the main steps, success and limitations in comparing real observations to synthetic observations from the non-ideal MHD models. Finally, we discuss the main threads of observational evidence that suggest a key role of magnetic fields for star and disk formation, and propose a scenario solving the angular momentum for star formation, also highlighting the remaining tensions that exist between models and observations. △ Less

Submitted 2 September, 2022; originally announced September 2022.

Comments: Accepted as review paper in Frontiers

Chemical Diversity in Protoplanetary Disks and Its Impact on the Formation History of Giant Planets

Authors: Elenia Pacetti, Diego Turrini, Eugenio Schisano, Sergio Molinari, Sergio Fonte, Romolo Politi, Patrick Hennebelle, Ralf Klessen, Leonardo Testi, Ugo Lebreuilly

Abstract: Giant planets can interact with multiple and chemically diverse environments in protoplanetary discs while they form and migrate to their final orbits. The way this interaction affects the accretion of gas and solids shapes the chemical composition of the planets and of their atmospheres. Here we investigate the effects of different chemical structures of the host protoplanetary disc on the planet… ▽ More Giant planets can interact with multiple and chemically diverse environments in protoplanetary discs while they form and migrate to their final orbits. The way this interaction affects the accretion of gas and solids shapes the chemical composition of the planets and of their atmospheres. Here we investigate the effects of different chemical structures of the host protoplanetary disc on the planetary composition. We consider both scenarios of molecular (inheritance from the pre-stellar cloud) and atomic (complete chemical reset) initial abundances in the disc. We focus on four elemental tracers of different volatility: C, O, N, and S. We explore the entire extension of possible formation regions suggested by observations by coupling the disc chemical scenarios with N-body simulations of forming and migrating giant planets. The planet formation process produces giant planets with chemical compositions significantly deviating from that of the host disc. We find that the C/N, N/O, and S/N ratios follow monotonic trends with the extent of migration. The C/O ratio shows a more complex behaviour, dependent on the planet accretion history and on the chemical structure of the formation environment. The comparison between S/N* and C/N* (where * indicates normalisation to the stellar value), constrains the relative contribution of gas and solids to the total metallicity. Giant planets whose metallicity is dominated by the contribution of the gas are characterised by N/O* > C/O* > C/N* and allow for constraining the disc chemical scenario. When the planetary metallicity is instead dominated by the contribution of the solids we find that C/N* > C/O* > N/O*. △ Less

Submitted 22 September, 2022; v1 submitted 29 June, 2022; originally announced June 2022.

Comments: 27 pages, 10 figures, 1 table. Published in The Astrophysical Journal

Journal ref: The Astrophysical Journal, 937-36, 2022

The signature of large scale turbulence driving on the structure of the interstellar medium

Authors: Tine Colman, Jean-François Robitaille, Patrick Hennebelle, Marc-Antoine Miville-Deschênes, Noé Brucy, Ralf S. Klessen, Simon C. O. Glover, Juan D. Soler, Davide Elia, Alessio Traficante, Sergio Molinari, Leonardo Testi

Abstract: The mechanisms that maintain turbulence in the interstellar medium (ISM) are still not identified. This work investigates how we can distinguish between two fundamental driving mechanisms: the accumulated effect of stellar feedback versus the energy injection from Galactic scales. We perform a series of numerical simulations describing a stratified star forming ISM subject to self-consistent stell… ▽ More The mechanisms that maintain turbulence in the interstellar medium (ISM) are still not identified. This work investigates how we can distinguish between two fundamental driving mechanisms: the accumulated effect of stellar feedback versus the energy injection from Galactic scales. We perform a series of numerical simulations describing a stratified star forming ISM subject to self-consistent stellar feedback. Large scale external turbulent driving of various intensities is added to mimic galactic driving mechanisms. We analyse the resulting column density maps with a technique called Multi-scale non-Gaussian segmentation that separates the coherent structures and the Gaussian background. This effectively discriminates between the various simulations and is a promising method to understand the ISM structure. In particular the power spectrum of the coherent structures flattens above 60 pc when turbulence is driven only by stellar feedback. When large-scale driving is applied, the turn-over shifts to larger scales. A systematic comparison with the Large Magellanic Cloud (LMC) is then performed. Only 1 out of 25 regions has a coherent power spectrum which is consistent with the feedback-only simulation. A detailed study of the turn-over scale leads us to conclude that regular stellar feedback is not enough to explain the observed ISM structure on scales larger than 60 pc. Extreme feedback in the form of supergiant shells likely plays an important role but cannot explain all the regions of the LMC. If we assume ISM structure is generated by turbulence, another large scale driving mechanism is needed to explain the entirety of the observations. △ Less

Submitted 1 June, 2022; originally announced June 2022.

Comments: 15 pages, 11 figures, to be published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 514, Issue 3, August 2022

The Galactic dynamics revealed by the filamentary structure in atomic hydrogen emission

Authors: Juan D. Soler, Marc-Antoine Miville-Deschênes, Sergio Molinari, Ralf S. Klessen, Patrick Hennebelle, Leonardo Testi, Naomi M. McClure-Griffiths, Henrik Beuther, Davide Elia, Eugenio Schisano, Alessio Traficante, Philipp Girichidis, Simon C. O. Glover, Rowan J. Smith, Mattia Sormani, Robin Treß

Abstract: We present a study of the filamentary structure in the atomic hydrogen (HI) emission at the 21 cm wavelength toward the Galactic plane using the observations in the HI4PI survey. Using the Hessian matrix method across radial velocity channels, we identified the filamentary structures and quantified their orientations using circular statistics. We found that the regions of the Milky Way's disk beyo… ▽ More We present a study of the filamentary structure in the atomic hydrogen (HI) emission at the 21 cm wavelength toward the Galactic plane using the observations in the HI4PI survey. Using the Hessian matrix method across radial velocity channels, we identified the filamentary structures and quantified their orientations using circular statistics. We found that the regions of the Milky Way's disk beyond 10 kpc and up to roughly 18 kpc from the Galactic center display HI filamentary structures predominantly parallel to the Galactic plane. For regions at lower Galactocentric radii, we found that the HI filaments are mostly perpendicular or do not have a preferred orientation with respect to the Galactic plane. We interpret these results as the imprint of supernova feedback in the inner Galaxy and Galactic rotation in the outer Milky Way. We found that the HI filamentary structures follow the Galactic warp and that they highlight some of the variations interpreted as the effect of the gravitational interaction with satellite galaxies. In addition, the mean scale height of the filamentary structures is lower than that sampled by the bulk of the HI emission, thus indicating that the cold and warm atomic hydrogen phases have different scale heights in the outer galaxy. Finally, we found that the fraction of the column density in HI filaments is almost constant up to approximately 18 kpc from the Galactic center. This is possibly a result of the roughly constant ratio between the cold and warm atomic hydrogen phases inferred from the HI absorption studies. Our results indicate that the HI filamentary structures provide insight into the dynamical processes shaping the Galactic disk. Their orientations record how and where the stellar energy input, the Galactic fountain process, the cosmic ray diffusion, and the gas accretion have molded the diffuse interstellar medium in the Galactic plane. △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: 33 pages. 39 figures. Accepted for publication in Astronomy & Astrophysics (06MAY2022)

Report number: AA/2022/43334

Journal ref: A&A 662, A96 (2022)

Cosmic-ray diffusion and the multi-phase interstellar medium in a dwarf galaxy. I. Large-scale properties and $γがんま$-ray luminosities

Authors: A. Nuñez-Castiñeyra, I. A. Grenier, F. Bournaud, Y. Dubois, F. R. Kamal Youssef, P. Hennebelle

Abstract: Dynamically, cosmic rays with energies above about one GeV/nucleon may be important agents of galaxy evolution. Their pressures compare with the thermal and magnetic ones impacting galactic gas accretion, fountains and galactic outflows, and alter the mass cycling between the gas phases, its efficiency depends on the properties of CR transport in the different media. We aim to study the dynamical… ▽ More Dynamically, cosmic rays with energies above about one GeV/nucleon may be important agents of galaxy evolution. Their pressures compare with the thermal and magnetic ones impacting galactic gas accretion, fountains and galactic outflows, and alter the mass cycling between the gas phases, its efficiency depends on the properties of CR transport in the different media. We aim to study the dynamical role of CRs in shaping the interstellar medium of a galaxy when changing their propagation mode. We perform MHD simulations with the AMR code RAMSES of the evolution of the same isolated galaxy (dwarf galaxy of $10^{11}$ M$_{\odot}$ down to 9-pc resolution) and compare the impact of the simplest cosmic-ray transport assumption of uniform diffusion. We have also updated the observational relation seen between the $γがんま$-ray luminosities and SFR of galaxies using the latest detection of Fermi LAT sources. We find that the radial and vertical distributions, and mass fractions of the gas in the different phases are marginally altered when changing CR transport. We observe positive feedback of CR on the amplification of the magnetic field in the inner half of the galaxy, except for fast isotropic diffusion. The increase in CR pressure for slow or anisotropic diffusion can suppress star formation by up to 50\%, but the dual effect of cosmic-ray pressure and magnetic amplification can reduce star formation by a factor 2.5. The $γがんま$-ray luminosities and SFR of the simulated galaxies are fully consistent with the trend seen in the observations in the case of anisotropic $10^{27.5-29}$ cm$^2$ s$^{-1}$ diffusion and for isotropic diffusion slower or equal to $3 \times 10^{28}$cm$^2$ s$^{-1}$. These results, therefore, do not confirm claims of very fast $10^{29-31}$ cm$^2$ s$^{-1}$ diffusion to match the Fermi LAT observations. △ Less

Submitted 17 May, 2022; originally announced May 2022.

Magnetically regulated collapse in the B335 protostar? II. Observational constraints on gas ionization and magnetic field coupling

Authors: Victoria Cabedo, Anaelle Maury, Josep Miquel Girart, Marco Padovani, Patrick Hennebelle, Martin Houde, Qizhou Zhang

Abstract: Non-ideal magnetohydrodynamic effects that rule the coupling of the magnetic field to the circumstellar gas during the low-mass star formation process depend heavily on the local physical conditions, such as the ionization fraction of the gas. The purpose of this work is to observationally characterize the level of ionization of the circumstellar gas at small envelope radii and investigate its rel… ▽ More Non-ideal magnetohydrodynamic effects that rule the coupling of the magnetic field to the circumstellar gas during the low-mass star formation process depend heavily on the local physical conditions, such as the ionization fraction of the gas. The purpose of this work is to observationally characterize the level of ionization of the circumstellar gas at small envelope radii and investigate its relation to the efficiency of the coupling between the star-forming gas and the magnetic field in the Class 0 protostar B335. We have obtained molecular line emission maps of B335 with ALMA, which we use to measure the deuteration fraction of the gas, its ionization fraction, and the cosmic-ray ionization rate, at envelope radii $\lesssim$1000 au. We find large fractions of ionized gas, $χかい_{e} \simeq 1-8 \times 10^{-6}$. Our observations also reveal an enhanced ionization that increases at small envelope radii, reaching values up to $ζぜーた_{CR} \simeq 10^{-14}$~s$^{-1}$ at a few hundred au from the central protostellar object. We show that this extreme ionization rate can be attributed to the presence of cosmic rays accelerated close to the protostar. We report the first resolved map of the cosmic-ray ionization rate at scales $\lesssim 1000$~au in a solar-type Class 0 protostar, finding remarkably high values. Our observations suggest that local acceleration of cosmic rays, and not the penetration of interstellar Galactic cosmic rays, may be responsible for the gas ionization in the inner envelope, potentially down to disk forming scales. If confirmed, our findings imply that protostellar disk properties may also be determined by local processes setting the coupling between the gas and the magnetic field, and not only by the amount of angular momentum available at large envelope scales and the magnetic field strength in protostellar cores. △ Less

Submitted 21 November, 2022; v1 submitted 21 April, 2022; originally announced April 2022.

Comments: A&A, Forthcoming article

Is the mm/submm dust polarization a robust tracer of the magnetic field topology in protostellar envelopes? A model exploration

Authors: Valeska Valdivia, Anaëlle Maury, Patrick Hennebelle

Abstract: High resolution (sub-)millimeter polarization observations have opened a new era in the understanding of how B-fields are organized in star forming regions, unveiling an intricate interplay between the B-fields and the gas in protostellar cores. However, to assess the role of the B-field in the process of solar-type star formation, it is key to be able to understand to what extent these polarized… ▽ More High resolution (sub-)millimeter polarization observations have opened a new era in the understanding of how B-fields are organized in star forming regions, unveiling an intricate interplay between the B-fields and the gas in protostellar cores. However, to assess the role of the B-field in the process of solar-type star formation, it is key to be able to understand to what extent these polarized dust emissions are good tracers of the B-field in the youngest protostellar objects. We present a thorough investigation of the fidelity and limitations of using dust polarized emission to map the B-field topologies in low-mass protostars. To assess the importance of these effects, we performed the analysis of B-field properties in 27 realizations of MHD models of star-forming cores. Assuming a uniform population of dust grains whose sizes follow the standard MRN, we analyze the synthetic polarized dust emission maps produced if these grains align with the local B-field thanks to B-RATs. We find that (sub-)millimeter polarized dust emission is a robust tracer of the B-field topologies in inner protostellar envelopes and is successful at capturing the details of the B-field spatial distribution down to radii ~100 au. Measurements of the los averaged B-field orientation using the polarized dust emission are precise to < 15° in about 75 - 95% of the independent lines of sight peering through protostellar envelopes. Large discrepancies between the integrated B-field mean orientation and the orientation reconstructed from the polarized dust emission are mostly observed in (i) lines of sight where the B-field is highly disorganized and (ii) lines of sight probing large column densities. Our analysis shows that high opacity of the thermal dust emission and low polarization fractions could be used to avoid utilizing the small fraction of measurements affected by large errors. △ Less

Submitted 9 September, 2022; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 668, A83 (2022)

ALMA-IMF III -- Investigating the origin of stellar masses: Top-heavy core mass function in the W43-MM2&MM3 mini-starburst

Authors: Y. Pouteau, F. Motte, T. Nony, R. Galván-Madrid, A. Men'shchikov, S. Bontemps, J. -F. Robitaille, F. Louvet, A. Ginsburg, F. Herpin, A. López-Sepulcre, P. Dell'Ova, A. Gusdorf, P. Sanhueza, A. M. Stutz, N. Brouillet, B. Thomasson, M. Armante, T. Baug, G. Busquet, T. Csengeri, N. Cunningham, M. Fernández-López, H. -L. Liu, F. Olguin , et al. (13 additional authors not shown)

Abstract: The ALMA-IMF Large Program observed the W43-MM2-MM3 ridge, whose 1.3mm and 3mm ALMA 12m array continuum images reach a 2500au spatial resolution. We used both the best-sensitivity and the line-free ALMA-IMF images, reduced the noise with the multi-resolution segmentation technique MnGSeg, and derived the most complete and most robust core catalog possible. Using two different extraction software p… ▽ More The ALMA-IMF Large Program observed the W43-MM2-MM3 ridge, whose 1.3mm and 3mm ALMA 12m array continuum images reach a 2500au spatial resolution. We used both the best-sensitivity and the line-free ALMA-IMF images, reduced the noise with the multi-resolution segmentation technique MnGSeg, and derived the most complete and most robust core catalog possible. Using two different extraction software packages, getsf and GExt2D, we identified 200 compact sources, whose 100 common sources have on average fluxes consistent to within 30%. We filtered sources with non-negligible free-free contamination and corrected fluxes from line contamination, resulting in a W43-MM2-MM3 catalog of 205 getsf cores. With a median deconvolved FWHM size of 3400au, core masses range from 0.1Msun to 70Msun and the getsf catalog is 90% complete down to 0.8Msun. The high-mass end of the core mass function (CMF) of W43-MM2-MM3 is top-heavy compared to the canonical IMF. Fitting the cumulative CMF with a single power law of the form N(>logM)\propto M^a, we measured a=-0.95\pm0.04, compared to the canonical a=-1.35 Salpeter IMF slope. The slope of the CMF is robust with respect to map processing, extraction software package, and reasonable variations in the assumptions taken to estimate core masses. We explore several assumptions on how cores transfer their mass to stars and sub-fragment to predict the IMF resulting from the W43-MM2-MM3 CMF. In stark contrast to the commonly accepted paradigm, our result argues against the universality of the CMF shape. More robust functions of the star-formation efficiency and core sub-fragmentation are required to better predict the resulting IMF, here suggested to remain top-heavy at the end of the star-formation phase. If confirmed, the IMFs emerging from starburst events could inherit their top-heavy shape from their parental CMFs, challenging the IMF universality. △ Less

Submitted 28 April, 2022; v1 submitted 7 March, 2022; originally announced March 2022.

Comments: 28 pages, 16 figures, 6 tables, A&A accepted on 02/08/22, revised following the language referee report

Journal ref: A&A 664, A26 (2022)

HyGAL: Characterizing the Galactic ISM with observations of hydrides and other small molecules -- I. Survey description and a first look toward W3(OH), W3 IRS5 and NGC 7538 IRS1

Authors: A. M. Jacob, D. A. Neufeld, P. Schilke, H. Wiesemeyer, W. Kim, S. Bialy, M. Busch, D. Elia, E. Falgarone, M. Gerin, B. Godard, R. Higgins, P. Hennebelle, N. Indriolo, D. C. Lis, K. M. Menten, A. Sanchez-Monge, V. Ossenkopf-Okada, M. R. Rugel, D. Seifried, P. Sonnentrucker, S. Walch, M. Wolfire, F. Wyrowski, V. Valdivia

Abstract: The HyGAL SOFIA legacy program surveys six hydride molecules -- ArH+, OH+, H2O+, SH, OH, and CH -- and two atomic constituents -- C+ and O -- within the diffuse interstellar medium (ISM) by means of absorption-line spectroscopy toward 25 bright Galactic background continuum sources. This detailed spectroscopic study is designed to exploit the unique value of specific hydrides as tracers and probes… ▽ More The HyGAL SOFIA legacy program surveys six hydride molecules -- ArH+, OH+, H2O+, SH, OH, and CH -- and two atomic constituents -- C+ and O -- within the diffuse interstellar medium (ISM) by means of absorption-line spectroscopy toward 25 bright Galactic background continuum sources. This detailed spectroscopic study is designed to exploit the unique value of specific hydrides as tracers and probes of different phases of the ISM, as demonstrated by recent studies with the Herschel Space Observatory. The observations performed under the HyGAL program will allow us to address several questions related to the lifecycle of molecular material in the ISM and the physical processes that impact its phase transition, such as: (1) What is the distribution function of the H2 fraction in the ISM? (2) How does the ionization rate due to low-energy cosmic-rays vary within the Galaxy? (3) What is the nature of interstellar turbulence, and what mechanisms lead to its dissipation? This overview discusses the observing strategy, synergies with ancillary and archival observations, the data reduction and analysis schemes adopted; and presents the first results obtained toward three of the survey targets, W3(OH), W3IRS5 and NGC7538IRS1. Robust measurements of the column densities of these hydrides -- obtained through widespread observations of absorption lines-- help address the questions raised, and there is a timely synergy between these observations and the development of theoretical models, particularly pertaining to the formation of H2 within the turbulent ISM. The provision of enhanced HyGAL data products will therefore serve as a legacy for future ISM studies. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Comments: 33 pages, 17 figures, 7 tables. Accepted for publication in the Astrophysical Journal