-
Chemical inventory of the envelope of the Class I protostar L1551 IRS 5
Authors:
P. Marchand,
A. Coutens,
J. Scigliuto,
F. Cruz-Sáenz de Miera,
A. Andreu,
J. -C. Loison,
Á. Kóspál,
P. Ábrahám
Abstract:
Episodic accretion in protostars leads to luminosity outbursts that end up heating their surroundings. This rise in temperature pushes the snow lines back, enabling the desorption of chemical species from dust grain surfaces, which may significantly alter the chemical history of the accreting envelope. However, a limited number of extensive chemical surveys of eruptive young stars have been perfor…
▽ More
Episodic accretion in protostars leads to luminosity outbursts that end up heating their surroundings. This rise in temperature pushes the snow lines back, enabling the desorption of chemical species from dust grain surfaces, which may significantly alter the chemical history of the accreting envelope. However, a limited number of extensive chemical surveys of eruptive young stars have been performed thus far. In the present study, we carry out a large spectral survey of the binary Class I protostar L1551 IRS 5, known to be a FUor-like object, in the 3mm and 2mm bands with the IRAM-30m telescope. As a result, we detected more than 400 molecular lines. The source displays a great chemical richness with the detection of 75 species, including isotopologues. Among these species, there are 13 hydrocarbons, 25 N-bearing species, 30 O-bearing species, 15 S-bearing species, 12 deuterated molecules, and a total of 10 complex organic molecules (l-C4H2, CH3CCH, CH2DCCH, CH3CHO, CH3CN, CH3OCH3, CH3OCHO, CH3OH, CH2DOH, and HC5N). With the help of local thermodynamic equilibrium (LTE) and non-LTE models, we determined the column densities of most molecules as well as excitation and kinetic temperatures. While most of those molecules trace the cold envelope (< 20 K), the OCS and CH3OH emission arise from the warm (> 100 K) innermost (< 2'' ) regions. We compared the chemical inventory of L1551 IRS 5 and its column density ratios, including isotopic ratios, with other protostellar sources. A broad chemical diversity is seen among Class I objects. More observations with both single-dish telescopes and interferometers are needed to characterize the diversity in a larger sample of protostars, while more astrochemical models would help explain this diversity, in addition to the impact of luminosity outbursts on the chemistry of protostellar envelopes.
△ Less
Submitted 15 May, 2024; v1 submitted 14 May, 2024;
originally announced May 2024.
-
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.
-
Fast methods for tracking grain coagulation and ionization. III. Protostellar collapse with non-ideal MHD
Authors:
Pierre Marchand,
Ugo Lebreuilly,
Mordecai-Mark Mac Low,
Vincent Guillet
Abstract:
Dust grains influence many aspects of star formation, including planet formation, opacities for radiative transfer, chemistry, and the magnetic field via Ohmic, Hall, and ambipolar diffusion. The size distribution of the dust grains is the primary characteristic influencing all these aspects. Grain size increases by coagulation throughout the star formation process. We describe here numerical simu…
▽ More
Dust grains influence many aspects of star formation, including planet formation, opacities for radiative transfer, chemistry, and the magnetic field via Ohmic, Hall, and ambipolar diffusion. The size distribution of the dust grains is the primary characteristic influencing all these aspects. Grain size increases by coagulation throughout the star formation process. We describe here numerical simulations of protostellar collapse using methods described in earlier papers of this series. We compute the evolution of the grain size distribution from coagulation and the non-ideal magnetohydrodynamics effects self-consistently and at low numerical cost. We find that the coagulation efficiency is mostly affected by the time spent in high-density regions. Starting from sub-micron radii, grain sizes reach more than 100 μm in an inner protoplanetary disk that is only 1000 years old. We also show that the growth of grains significantly affects the resistivities, and indirectly the dynamics and angular momentum of the disk.
△ Less
Submitted 9 January, 2023; v1 submitted 4 January, 2023;
originally announced January 2023.
-
Protostellar collapse simulations in spherical geometry with dust coagulation and fragmentation
Authors:
Ugo Lebreuilly,
Valentin Vallucci-Goy,
Vincent Guillet,
Maxime Lombart,
Pierre Marchand
Abstract:
We model the coagulation and fragmentation of dust grains during the protostellar collapse with our newly developed shark code. It solves the gas-dust hydrodynamics in a spherical geometry and the coagulation/fragmentation equation. It also computes the ionization state of the cloud and the Ohmic, ambipolar and Hall resistivities. We find that the dust size distribution evolves significantly durin…
▽ More
We model the coagulation and fragmentation of dust grains during the protostellar collapse with our newly developed shark code. It solves the gas-dust hydrodynamics in a spherical geometry and the coagulation/fragmentation equation. It also computes the ionization state of the cloud and the Ohmic, ambipolar and Hall resistivities. We find that the dust size distribution evolves significantly during the collapse, large grain formation being controlled by the turbulent differential velocity. When turbulence is included, only ambipolar diffusion remains efficient at removing the small grains from the distribution, brownian motion is only efficient as a standalone process. The macroscopic gas-dust drift is negligible for grain growth and only dynamically significant near the first Larson core. At high density, we find that the coagulated distribution is unaffected by the initial choice of dust distribution. Strong magnetic fields are found to enhance the small grains depletion, causing an important increase of the ambipolar diffusion. This hints that the magnetic field strength could be regulated by the small grain population during the protostellar collapse. Fragmentation could be effective for bare silicates, but its modeling relies on the choice of ill-constrained parameters. It is also found to be negligible for icy grains. When fragmentation occurs, it strongly affects the magnetic resistivities profiles. Dust coagulation is a critical process that needs to be fully taken into account during the protostellar collapse. The onset and feedback of fragmentation remains uncertain and its modeling should be further investigated.
△ Less
Submitted 3 November, 2022;
originally announced November 2022.
-
Fast methods to track grain coagulation and ionization. II. Extension to thermal ionization
Authors:
Pierre Marchand,
Vincent Guillet,
Ugo Lebreuilly,
Mordecai-Mark Mac Low
Abstract:
Thermal ionization is a critical process at temperatures T > 10 3 K, particularly during star formation. An increase in ionization leads to a decrease in nonideal magnetohydrodynamics (MHD) resistivities, which has a significant impact on protoplanetary disks and protostar formation. We developed an extension of the fast computational ionization method presented in our recent paper to include ther…
▽ More
Thermal ionization is a critical process at temperatures T > 10 3 K, particularly during star formation. An increase in ionization leads to a decrease in nonideal magnetohydrodynamics (MHD) resistivities, which has a significant impact on protoplanetary disks and protostar formation. We developed an extension of the fast computational ionization method presented in our recent paper to include thermal ionization. The model can be used to inexpensively calculate the density of ions and electrons and the electric charge of each size of grains for an arbitrary size distribution. This tool should be particularly useful for the self-consistent calculation of nonideal MHD resistivities in multidimensional simulations, especially of protostellar collapse and protoplanetary disks.
△ Less
Submitted 30 June, 2022; v1 submitted 23 February, 2022;
originally announced February 2022.
-
Universal protoplanetary disk size under complete non-ideal magnetohydrodynamics: The interplay between ion-neutral friction, Hall effect, and the Ohmic dissipation
Authors:
Yueh-Ning Lee,
Pierre Marchand,
Yu-Hsuan Liu,
Patrick Hennebelle
Abstract:
The role of non-ideal magnetohydrodynamics has been proven critical during the formation of the protoplanetary disk, particularly in regulating its size. We provide a simple model to predict the disk size under the interplay among the ambipolar diffusion, the Hall effect, and the Ohmic dissipation. The model predicts a small disk size of around 20 AU, that depends only sub-linearly on disk paramet…
▽ More
The role of non-ideal magnetohydrodynamics has been proven critical during the formation of the protoplanetary disk, particularly in regulating its size. We provide a simple model to predict the disk size under the interplay among the ambipolar diffusion, the Hall effect, and the Ohmic dissipation. The model predicts a small disk size of around 20 AU, that depends only sub-linearly on disk parameters, for a wide range of initial conditions of sub-Solar mass and moderate magnetization. It is able to explain phenomena manifested in existing numerical simulations, including the bimodal disk behavior under parallel and anti-parallel alignment between the rotation and magnetic field. In the parallel configuration, the disk size decreases and eventually disappears. In the anti-parallel configuration, and the disk has an outer partition (or pseudo-disk) that is flat, shrinking , and short-lived, as well as a inner partition that grows slowly with mass and is long-lived. Even with significant initial magnetization, the vertical field in the disk can only dominate at the early stage when the mass is low, and the toroidal field eventually dominates in all disks.
△ Less
Submitted 13 September, 2021;
originally announced September 2021.
-
Fast methods for tracking grain coagulation and ionization. I. Analytic derivation
Authors:
Pierre Marchand,
Vincent Guillet,
Ugo Lebreuilly,
Mordecai-Mark Mac Low
Abstract:
Dust grains play a major role in many astrophysical contexts. They affect the chemical, magnetic, dynamical, and optical properties of their environment, from galaxies down to the interstellar medium, star-forming regions, and protoplanetary disks. Their coagulation leads to shifts in their size distribution and ultimately to the formation of planets. However, although the coagulation process is r…
▽ More
Dust grains play a major role in many astrophysical contexts. They affect the chemical, magnetic, dynamical, and optical properties of their environment, from galaxies down to the interstellar medium, star-forming regions, and protoplanetary disks. Their coagulation leads to shifts in their size distribution and ultimately to the formation of planets. However, although the coagulation process is reasonably uncomplicated to numerically implement by itself, it is difficult to couple it with multidimensional hydrodynamics numerical simulations because of its high computational cost. We propose here a simple method for tracking the coagulation of grains at far lower cost. Given an initial grain size distribution, the state of the distribution at time t is solely determined by the value of a single variable integrated along the trajectory, independently of the specific path taken by the grains. Although this method cannot account for processes other than coagulation, it is mathematically exact, fast, inexpensive, and can be used to evaluate the effect of grain coagulation in most astrophysical contexts. It is applicable to all coagulation kernels in which local physical conditions and grain properties can be separated. We also describe another method for calculating the average electric charge of grains and the density of ions and electrons in environments that are shielded from radiation fields, given the density and temperature of the gas, the cosmic-ray ionization rate, and the average mass of the ions. The equations we provide are fast to integrate numerically and can be used in multidimensional numerical simulations to self-consistently calculate on the fly the local resistivities that are required to model nonideal magnetohydrodynamics.
△ Less
Submitted 17 May, 2021; v1 submitted 25 February, 2021;
originally announced March 2021.
-
Protostellar collapse: regulation of the angular momentum and onset of an ionic precursor
Authors:
Pierre Marchand,
Kengo Tomida,
Kei Tanaka,
Benoît Commerçon,
Gilles Chabrier
Abstract:
Through the magnetic braking and the launching of protostellar outflows, magnetic fields play a major role in the regulation of angular momentum in star formation, which directly impacts the formation and evolution of protoplanetary disks and binary systems. The aim of this paper is to quantify those phenomena in the presence of non-ideal magnetohydrodynamics effects, namely the Ohmic and ambipola…
▽ More
Through the magnetic braking and the launching of protostellar outflows, magnetic fields play a major role in the regulation of angular momentum in star formation, which directly impacts the formation and evolution of protoplanetary disks and binary systems. The aim of this paper is to quantify those phenomena in the presence of non-ideal magnetohydrodynamics effects, namely the Ohmic and ambipola r diffusion. We perform three-dimensional simulations of protostellar collapses varying the mass of the prestellar dense core, the thermal support (the $α$ ratio) and the dust grain size-distribu tion. The mass mostly influences the magnetic braking in the pseudo-disk, while the thermal support impacts the accretion rate and hence the properties of the disk. Removing the grains smaller than 0. 1 $μ$m in the Mathis, Rumpl, Nordsieck (MRN) distribution enhances the ambipolar diffusion coefficient. Similarly to previous studies, we find that this change in the distribution reduces the magnet ic braking with an impact on the disk. The outflow is also significantly weakened. In either case, the magnetic braking largely dominates the outflow as a process to remove the angular momentum from t he disk. Finally, we report a large ionic precursor to the outflow with velocities of several km s$^{-1}$, which may be observable.
△ Less
Submitted 2 September, 2020;
originally announced September 2020.
-
Dust coagulation feedback on magnetohydrodynamic resistivities in protostellar collapse
Authors:
V. Guillet,
P. Hennebelle,
G. Pineau des Forêts,
A. Marcowith,
B. Commerçon,
P. Marchand
Abstract:
The degree of coupling between the gas and the magnetic field during the collapse of a core and the subsequent formation of a disk depends on the assumed dust size distribution. We study the impact of grain-grain coagulation on the evolution of magnetohydrodynamic (MHD) resistivities during the collapse of a prestellar core. We use a 1-D model to follow the evolution of the dust size distribution,…
▽ More
The degree of coupling between the gas and the magnetic field during the collapse of a core and the subsequent formation of a disk depends on the assumed dust size distribution. We study the impact of grain-grain coagulation on the evolution of magnetohydrodynamic (MHD) resistivities during the collapse of a prestellar core. We use a 1-D model to follow the evolution of the dust size distribution, out-of-equilibrium ionization state and gas chemistry during the collapse of a prestellar core. To compute the grain-grain collisional rate, we consider models for both random and systematic, size-dependent, velocities. We include grain growth through grain-grain coagulation and ice accretion, but ignore grain fragmentation. Starting with a MRN (Mathis et al. 1977) size distribution, we find that coagulation in grain-grain collisions generated by hydrodynamical turbulence is not efficient at removing the smallest grains, and as a consequence does not affect much the evolution of the Hall and ambipolar diffusion MHD resistivities which still severly drop during the collapse like in models without coagulation. The inclusion of systematic velocities, possibly induced by the presence of ambipolar diffusion, increases the coagulation rate between small and large grains, removing small grains earlier in the collapse and therefore limiting the drop in the Hall and ambipolar diffusion resistivities. At intermediate densities ($n_{\rm H} \sim 10^8\,$cm$^{-3}$), the Hall and ambipolar diffusion resistivities are found to be higher by 1 to 2 orders of magnitude in models with coagulation than in models where coagulation is ignored, and also higher than in a toy model without coagulation where all grains smaller than $0.1\,μ$m would have been removed in the parent cloud before the collapse. When grain drift velocities induced by ambipolar diffusion are included, dust coagulation ... (abridged)
△ Less
Submitted 8 July, 2020;
originally announced July 2020.
-
Impact of the Hall effect in star formation : improving the angular momentum conservation
Authors:
Pierre Marchand,
Kengo Tomida,
Benoît Commerçon,
Gilles Chabrier
Abstract:
We present here a minor modification of our numerical implementation of the Hall effect for the 2D Riemann solver used in Constrained Transport schemes, as described in Marchand et al. (2018). In the previous work, the tests showed that the angular momentum was not conserved during protostellar collapse simulations, with significant impact. By removing the whistler waves speed from the characteris…
▽ More
We present here a minor modification of our numerical implementation of the Hall effect for the 2D Riemann solver used in Constrained Transport schemes, as described in Marchand et al. (2018). In the previous work, the tests showed that the angular momentum was not conserved during protostellar collapse simulations, with significant impact. By removing the whistler waves speed from the characteristic speeds of non-magnetic variables in the 1D Riemann solver, we are able to improve the angular momentum conservation in our test-case by one order of magnitude, while keeping the second-order numerical convergence of the scheme. We also reproduce the simulations of Tsukamoto et al. (2015) with consistent resistivities, the three non-ideal MHD effects and initial rotation, and agree with their results. In this case, the violation of angular momentum conservation is negligible in regard to the total angular momentum and the angular momentum of the disk.
△ Less
Submitted 28 October, 2019; v1 submitted 19 September, 2019;
originally announced September 2019.
-
Impact of the Hall effect in star formation and the issue of angular momentum conservation
Authors:
Pierre Marchand,
Benoît Commerçon,
Gilles Chabrier
Abstract:
We present an implementation of the Hall term in the non-ideal magnetohydrodynamics equations into the adaptive-mesh-refinement code RAMSES to study its impact on star formation. Recent works show that the Hall effect heavily influences the regulation of the angular momentum in collapsing dense cores, strengthening or weakening the magnetic braking. Our method consists of a modification of the two…
▽ More
We present an implementation of the Hall term in the non-ideal magnetohydrodynamics equations into the adaptive-mesh-refinement code RAMSES to study its impact on star formation. Recent works show that the Hall effect heavily influences the regulation of the angular momentum in collapsing dense cores, strengthening or weakening the magnetic braking. Our method consists of a modification of the two-dimensional constrained transport scheme. Our scheme shows convergence of second-order in space and the frequency of the propagation of whistler waves is accurate. We confirm previous results, namely that during the collapse, the Hall effect generates a rotation of the fluid with a direction in the mid-plane that depends on the sign of the Hall resistivity, while counter-rotating envelopes develop on each side of the mid-plane. However, we find that the predictability of our numerical results is severely limited. The angular momentum is not conserved in any of our dense core-collapse simulations with the Hall effect: a large amount of angular momentum is generated within the first Larson core, a few hundred years after its formation, without compensation by the surrounding gas. This issue is not mentioned in previous studies and may be correlated to the formation of the accretion shock on the Larson core. We expect that this numerical effect could be a serious issue in star formation simulations.
△ Less
Submitted 27 August, 2018;
originally announced August 2018.
-
Magnetically self-regulated formation of early protoplanetary discs
Authors:
Patrick Hennebelle,
Benoit Commercon,
Gilles Chabrier,
Pierre Marchand
Abstract:
The formation of protoplanetary discs during the collapse of molecular dense cores is significantly influenced by angular momentum transport, notably by the magnetic torque. In turn, the evolution of the magnetic field is determined by dynamical processes and non-ideal MHD effects such as ambipolar diffusion. Considering simple relations between various timescales characteristic of the magnetized…
▽ More
The formation of protoplanetary discs during the collapse of molecular dense cores is significantly influenced by angular momentum transport, notably by the magnetic torque. In turn, the evolution of the magnetic field is determined by dynamical processes and non-ideal MHD effects such as ambipolar diffusion. Considering simple relations between various timescales characteristic of the magnetized collapse, we derive an expression for the early disc radius, $ r \simeq 18 \, {\rm AU} \, \left({η_{\rm AD} / 0.1 \, {\rm s}} \right)^{2/9} \left({B_z / 0.1\, {\rm G}} \right) ^{-4/9} \left({M / 0.1 \msol} \right) ^{1/3},$ where $M$ is the total disc plus protostar mass, $η_\mathrm{AD}$ is the ambipolar diffusion coefficient and $B_z$ is the magnetic field in the inner part of the core. This is about significantly smaller than the discs that would form if angular momentum was conserved.
The analytical predictions are confronted against a large sample of 3D, non-ideal MHD collapse calculations covering variations of a factor 100 in core mass, a factor 10 in the level of turbulence, a factor 5 in rotation, and magnetic mass-to-flux over critical mass-to-flux ratios 2 and 5. The disc radius estimates are found to agree with the numerical simulations within less than a factor 2.
A striking prediction of our analysis is the weak dependence of circumstellar disc radii upon the various relevant quantities, suggesting weak variations among class-0 disc sizes. In some cases, we note the onset of large spiral arms beyond this radius.
△ Less
Submitted 8 August, 2016;
originally announced August 2016.
-
A chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core collapse calculations
Authors:
Pierre Marchand,
Jacques Masson,
Gilles Chabrier,
Patrick Hennebelle,
Benoit Commerçon,
Neil Vaytet
Abstract:
We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities…
▽ More
We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of Potassium, Sodium and Hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to n_H = 10^12 cm^-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.
△ Less
Submitted 12 June, 2016; v1 submitted 19 April, 2016;
originally announced April 2016.