-
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.
-
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.
-
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.
-
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.
-
Calibrating the Davis-Chandrasekhar-Fermi method with numerical simulations: uncertainties in estimating the magnetic field strength from statistics of field orientations
Authors:
Junhao Liu,
Qizhou Zhang,
Benoit Commercon,
Valeska Valdivia,
Anaelle Maury,
Keping Qiu
Abstract:
The Davis-Chandrasekhar-Fermi (DCF) method is widely used to indirectly estimate the magnetic field strength from the plane-of-sky field orientation. In this work, we present a set of 3D MHD simulations and synthetic polarization images using radiative transfer of clustered massive star-forming regions. We apply the DCF method on the synthetic polarization maps to investigate its reliability in hi…
▽ More
The Davis-Chandrasekhar-Fermi (DCF) method is widely used to indirectly estimate the magnetic field strength from the plane-of-sky field orientation. In this work, we present a set of 3D MHD simulations and synthetic polarization images using radiative transfer of clustered massive star-forming regions. We apply the DCF method on the synthetic polarization maps to investigate its reliability in high-density molecular clumps and dense cores where self-gravity is significant. We investigate the validity of the assumptions of the DCF method step by step and compare the model and estimated field strength to derive the correction factors for the estimated uniform and total (rms) magnetic field strength at clump and core scales. The correction factors in different situations are catalogued. We find the DCF method works well in strong field cases. However, the magnetic field strength in weak field cases could be significantly overestimated by the DCF method when the turbulent magnetic energy is smaller than the turbulent kinetic energy. We investigate the accuracy of the angular dispersion function (ADF, a modified DCF method) method on the effects that may affect the measured angular dispersion and find that the ADF method correctly accounts for the ordered field structure, the beam-smoothing, and the interferometric filtering, but may not be applicable to account for the signal integration along the line of sight in most cases. Our results suggest that the DCF methods should be avoided to be applied below $\sim$0.1 pc scales if the effect of line-of-sight signal integration is not properly addressed.
△ Less
Submitted 12 July, 2021; v1 submitted 18 June, 2021;
originally announced June 2021.
-
An observational correlation between magnetic field, angular momentum and fragmentation in the envelopes of Class 0 protostars?
Authors:
Maud Galametz,
Anaelle Maury,
Josep M. Girart,
Ramprasad Rao,
Qizhou Zhang,
Mathilde Gaudel,
Valeska Valdivia,
Patrick Hennebelle,
Victoria Cabedo-Soto,
Eric Keto,
Shih-Ping Lai
Abstract:
To assess the potential role of magnetic fields in regulating the envelope rotation and the fragmentation of Class 0 protostars, we carried out observations of the dust polarized emission at 0.87 mm with the SMA, in the envelopes of a large sample of 20 Class 0 protostars. We estimate the mean magnetic field orientation over the central 1000 au envelope scales and compared it to that of the protos…
▽ More
To assess the potential role of magnetic fields in regulating the envelope rotation and the fragmentation of Class 0 protostars, we carried out observations of the dust polarized emission at 0.87 mm with the SMA, in the envelopes of a large sample of 20 Class 0 protostars. We estimate the mean magnetic field orientation over the central 1000 au envelope scales and compared it to that of the protostellar outflow in order to study the relation between their misalignment and the kinematics of the circumstellar gas. We discover a strong relationship between the misalignment of the magnetic field orientation with the outflow and the amount of angular momentum observed at similar scales in the protostellar envelope, revealing a potential link between the kinetic and the magnetic energy at envelope scales. The relation could be driven by favored B misalignments in more dynamical envelopes or a dependence of the envelope dynamics with the large-scale B initial configuration. Comparing the trend with the presence of fragmentation, we observe that single sources are mostly associated with conditions of low angular momentum in the inner envelope and good alignment of the magnetic field with protostellar outflows, at intermediate scales. Our results suggest that the properties of the magnetic field in protostellar envelopes bear a tight relationship with the rotating-infalling gas directly involved in the star and disk formation: we find that it may not only influence the fragmentation of protostellar cores into multiple stellar systems, but also set the conditions establishing the pristine properties of planet-forming disks.
△ Less
Submitted 27 January, 2021; v1 submitted 23 October, 2020;
originally announced October 2020.
-
A statistical analysis of dust polarization properties in ALMA observations of Class 0 protostellar cores
Authors:
V. J. M. Le Gouellec,
A. J. Maury,
V. Guillet,
C. L. H. Hull,
J. M. Girart,
A. Verliat,
R. Mignon-Risse,
V. Valdivia,
P. Hennebelle,
M. González,
F. Louvet
Abstract:
Recent observational progress has challenged the dust grain-alignment theories used to explain the polarized dust emission routinely observed in star-forming cores. In an effort to improve our understanding of the dust grain alignment mechanism(s), we have gathered a dozen ALMA maps of (sub)millimeter-wavelength polarized dust emission from Class 0 protostars, and carried out a comprehensive stati…
▽ More
Recent observational progress has challenged the dust grain-alignment theories used to explain the polarized dust emission routinely observed in star-forming cores. In an effort to improve our understanding of the dust grain alignment mechanism(s), we have gathered a dozen ALMA maps of (sub)millimeter-wavelength polarized dust emission from Class 0 protostars, and carried out a comprehensive statistical analysis of dust polarization quantities. We analyze the statistical properties of the polarization fraction P_frac and dispersion of polarization position angles S. More specifically, we investigate the relationship between S and P_frac as well as the evolution of the product S*P_frac as a function of the column density of the gas in the protostellar envelopes. We find a significant correlation in the polarized dust emission from protostellar envelopes seen with ALMA; the power-law index differs significantly from the one observed by Planck in star-forming clouds. The product S*P_frac, which is sensitive to the dust grain alignment efficiency, is approximately constant across three orders of magnitude in envelope column density. This suggests that the grain alignment mechanism producing the bulk of the polarized dust emission in star-forming cores may not depend systematically on the local conditions such as local gas density. Ultimately, our results suggest dust alignment mechanism(s) are efficient at producing dust polarized emission in the various local conditions typical of Class 0 protostars. The grain alignment efficiency found in these objects seems to be higher than the efficiency produced by the standard RAT alignment of paramagnetic grains. Further study will be needed to understand how more efficient grain alignment via, e.g., different irradiation conditions, dust grain characteristics, or additional grain alignment mechanisms can reproduce the observations.
△ Less
Submitted 15 September, 2020;
originally announced September 2020.
-
3D chemical structure of diffuse turbulent ISM. I. Statistics of the HI-to-H$_2$ transition
Authors:
Elena Bellomi,
Benjamin Godard,
Patrick Hennebelle,
Valeska Valdivia,
Guillaume Pineau des Forêts,
Pierre Lesaffre,
Michel Pérault
Abstract:
We studied the statistical properties of the HI-to-H$_2$ transition observed in absorption in the local diffuse and multiphase ISM to identify the physical processes controlling the probability of occurrence of any line of sight. The turbulent diffuse ISM is modeled using the RAMSES code, which includes detailed treatments of the magnetohydrodynamics, the thermal evolution of the gas, and the chem…
▽ More
We studied the statistical properties of the HI-to-H$_2$ transition observed in absorption in the local diffuse and multiphase ISM to identify the physical processes controlling the probability of occurrence of any line of sight. The turbulent diffuse ISM is modeled using the RAMSES code, which includes detailed treatments of the magnetohydrodynamics, the thermal evolution of the gas, and the chemistry of H$_2$. The impacts of the UV radiation field, the mean density, the turbulent forcing, the integral scale, the magnetic field, and the gravity on the molecular content of the gas are explored through a parametric study covering a wide range of physical conditions. The statistics of the HI-to-H$_2$ transition are interpreted through analytical prescriptions and compared with the observations using a modified and robust version of the Kolmogorov-Smirnov test. The results of one simulation, convolved with the distribution of distances of the observational sample, are able to explain most of the statistical properties of the HI-to-H$_2$ transition observed in the local ISM. The tightest agreement is obtained for a neutral diffuse gas modeled over ~200 pc, with a mean density of $1-2$ cm$^{-3}$, illuminated by the standard interstellar UV radiation field, and stirred up by a large-scale compressive turbulent forcing. Within this configuration, the 2D probability histogram of the column densities of H and H$_2$ is remarkably stable and is almost unaltered by gravity, the strength of the turbulent forcing, the resolution of the simulation, or the strength of the magnetic field $B_x$. The weak effect of the resolution and our analytical prescription suggest that the column densities of HI are likely built up in large-scale WNM and CNM structures correlated in density over ~20 pc and ~10 pc, respectively, while those of H$_2$ are built up in CNM structures between ~3 pc and ~10 pc.
△ Less
Submitted 25 October, 2020; v1 submitted 11 September, 2020;
originally announced September 2020.
-
Low dust emissivities and radial variations in the envelopes of Class 0 protostars: a signature of early grain growth?
Authors:
Maud Galametz,
Anaelle J. Maury,
Valeska Valdivia,
Leonardo Testi,
Arnaud Belloche,
Philippe Andre
Abstract:
Analyzing the properties of dust and its evolution in the early phases of star formation is crucial to put constraints on the collapse and accretion processes as well as on the pristine properties of planet-forming seeds. We analyze PdBI observations at 1.3 and 3.2 mm for 12 Class 0 protostars obtained as part of the CALYPSO survey and derive dust emissivity index (beta) profiles as a function of…
▽ More
Analyzing the properties of dust and its evolution in the early phases of star formation is crucial to put constraints on the collapse and accretion processes as well as on the pristine properties of planet-forming seeds. We analyze PdBI observations at 1.3 and 3.2 mm for 12 Class 0 protostars obtained as part of the CALYPSO survey and derive dust emissivity index (beta) profiles as a function of the envelope radius at 200-2000 au scales. Most of the protostellar envelopes show low dust emissivity indices decreasing toward the central regions. The decreasing trend remains after correction of the (potentially optically-thick) central region emission, with surprisingly low beta (lower than 1) values across most of the envelope radii of NGC1333-IRAS4A, NGC1333-IRAS4B, SVS13B, and Serpens-SMM4. We discuss the various processes that could explain such low and varying dust emissivity indices at envelope radii 200-2000 au. Our observations of extremely low dust emissivity indices could trace the presence of large (mm-size) grains in Class 0 envelopes and a radial increase of the dust grain size toward the inner envelope regions. While it is expected that large grains in young protostellar envelopes could be built via grain growth and coagulation, we stress that the typical timescales required to build mm grains in current coagulation models are at odds with the youth of our Class 0 protostars. Additional variations in the dust composition could also partly contribute to the low beta we observe. The steepness of the beta radial gradient depends strongly on the envelope mass, which might favor a scenario where large grains are built in high density protostellar disks and transported to the intermediate envelope radii, for example with the help of outflows and winds.
△ Less
Submitted 12 November, 2019; v1 submitted 10 October, 2019;
originally announced October 2019.
-
ESA Voyage 2050 white paper: A complete census of the gas phases in and around galaxies, far-UV spectropolarimetry as a prime tool for understanding galaxy evolution and star formation
Authors:
V. Lebouteiller,
C. Gry. H. Yan,
P. Richter,
B. Godard,
E. B. Jenkins,
D. Welty,
N. Lehner,
P. Guillard,
J. Roman-Duval,
E Roueff,
F. Leone,
D. Kunth,
J. C. Howk,
P. Boissé,
F. Boulanger,
E. Bron,
B. James,
J. Le Bourlot,
F. Le Petit,
M. Pieri,
V. Valdivia
Abstract:
(abridged) The far-UV wavelength range (912-2000A) provides access to atomic and molecular transitions of many species the interstellar medium (ISM), circumgalactic medium (CGM), and intergalactic medium, within phases spanning a wide range of ionization, density, temperature, and molecular gas fraction. Far-UV space telescopes have enabled detailed studies of the ISM in the Milky Way thanks to ab…
▽ More
(abridged) The far-UV wavelength range (912-2000A) provides access to atomic and molecular transitions of many species the interstellar medium (ISM), circumgalactic medium (CGM), and intergalactic medium, within phases spanning a wide range of ionization, density, temperature, and molecular gas fraction. Far-UV space telescopes have enabled detailed studies of the ISM in the Milky Way thanks to absorption features appearing in the UV spectra of hot stars and yielding fundamental insights into the composition and physical characteristics of all phases of the ISM along with the processes that influence them. However, we have yet to design a spectrometer able to observe the full UV domain at resolving power R>10^5 with a signal-to-noise ratio SNR>500. Such a resolution is necessary to resolve lines from both the cold molecular hydrogen and the warm metal ions with a turbulent velocity of about 1 km s-1, and to differentiate distinct velocity components. Future UV spectroscopic studies of the Milky Way ISM must revolutionize our understanding of the ISM as a dynamical, unstable, and magnetized medium, and rise to the challenge brought forward by current theories. Another interesting prospect is to transpose the same level of details that has been reached for the Milky Way to the ISM in external galaxies, in particular in metal-poor galaxies, where the ISM chemical composition, physical conditions, and topology change dramatically, with significant consequences on the star-formation properties. Finally, we need to be able to perform statistical analyses of background quasar lines of sight intersecting the CGM of galaxies at various redshifts and to comprehend the role of gas exchanges and flows for galaxy evolution.
△ Less
Submitted 6 September, 2019;
originally announced September 2019.
-
Indirect evidence of significant grain growth in young protostellar envelopes from polarized dust emission
Authors:
Valeska Valdivia,
Anaelle Maury,
Robert Brauer,
Patrick Hennebelle,
Maud Galametz,
Vincent Guillet,
Stefan Reissl
Abstract:
How and when in the star formation sequence do dust grains start to grow into pebbles is a cornerstone question to both star and planet formation. We compute the polarized radiative transfer from a model solar-type protostellar core, using the POLARIS code, aligning the dust grains with the local magnetic field, following the radiative torques (RATs) theory. We test the dependency of the resulting…
▽ More
How and when in the star formation sequence do dust grains start to grow into pebbles is a cornerstone question to both star and planet formation. We compute the polarized radiative transfer from a model solar-type protostellar core, using the POLARIS code, aligning the dust grains with the local magnetic field, following the radiative torques (RATs) theory. We test the dependency of the resulting dust polarized emission with the maximum grain size of the dust size distribution at the envelope scale, from amax = 1 micron to 50 micron. Our work shows that, in the framework of RAT alignment, large dust grains are required to produce polarized dust emission at levels similar to those currently observed in solar-type protostellar envelopes at millimeter wavelengths. Considering the current theoretical dificulties to align a large fraction of small ISM-like grains in the conditions typical of protostellar envelopes, our results suggest that grain growth (typically > 10 micron) might have already significantly progressed at scales 100-1000 au in the youngest objects, observed less than 10^5 years after the onset of collapse. Observations of dust polarized emission might open a new avenue to explore dust pristine properties and describe, for example, the initial conditions for the formation of planetesimals.
△ Less
Submitted 25 July, 2019;
originally announced July 2019.
-
SMA observations of the polarized dust emission in solar-type Class 0 protostars: the magnetic field properties at envelope scales
Authors:
Maud Galametz,
Anaelle Maury,
Josep M. Girart,
Ramprasad Rao,
Qizhou Zhang,
Mathilde Gaudel,
Valeska Valdivia,
Eric Keto,
Shih-Ping Lai
Abstract:
Although, from a theoretical point of view, magnetic fields are believed to have a significant role during the early stages of star formation, especially during the main accretion phase, the magnetic field strength and topology is poorly constrained in the youngest accreting Class 0 protostars that lead to the formation of solar-type stars. We carried out observations of the polarized dust continu…
▽ More
Although, from a theoretical point of view, magnetic fields are believed to have a significant role during the early stages of star formation, especially during the main accretion phase, the magnetic field strength and topology is poorly constrained in the youngest accreting Class 0 protostars that lead to the formation of solar-type stars. We carried out observations of the polarized dust continuum emission with the SMA interferometer at 0.87mm, in order to probe the structure of the magnetic field in a sample of 12 low-mass Class 0 envelopes, including both single protostars and multiple systems, in nearby clouds. Our SMA observations probe the envelope emission at scales of 600-5000 au with a spatial resolution ranging from 600 to 1500 au depending on the source distance. We report the detection of linearly polarized dust continuum emission in all of our targets, with average polarization fractions ranging from 2% to 10% in these protostellar envelopes. The polarization fraction decreases with the continuum flux density, which translates into a decrease with the H2 column density within an individual envelope. Our analysis show that the envelope-scale magnetic field is preferentially observed either aligned or perpendicular to the outflow direction. Interestingly, our results suggest for the first time a relation between the orientation of the magnetic field and the rotational energy of envelopes, with a larger occurrence of misalignment in sources where strong rotational motions are detected at hundreds to thousands of au scales. We also show that the best agreement between the magnetic field and outflow orientation is found in sources showing no small-scale multiplicity and no large disks at about 100 au scales.
△ Less
Submitted 16 April, 2018;
originally announced April 2018.
-
H$_2$ formation on interstellar dust grains: the viewpoints of theory, experiments, models and observations
Authors:
Valentine Wakelam,
Emeric Bron,
Stephanie Cazaux,
Francois Dulieu,
Cécile Gry,
Pierre Guillard,
Emilie Habart,
Liv Hornekær,
Sabine Morisset,
Gunnar Nyman,
Valerio Pirronello,
Stephen D. Price,
Valeska Valdivia,
Gianfranco Vidali,
Naoki Watanabe
Abstract:
Molecular hydrogen is the most abundant molecule in the universe. It is the first one to form and survive photo-dissociation in tenuous environments. Its formation involves catalytic reactions on the surface of interstellar grains. The micro-physics of the formation process has been investigated intensively in the last 20 years, in parallel of new astrophysical observational and modeling progresse…
▽ More
Molecular hydrogen is the most abundant molecule in the universe. It is the first one to form and survive photo-dissociation in tenuous environments. Its formation involves catalytic reactions on the surface of interstellar grains. The micro-physics of the formation process has been investigated intensively in the last 20 years, in parallel of new astrophysical observational and modeling progresses. In the perspectives of the probable revolution brought by the future satellite JWST, this article has been written to present what we think we know about the H$_2$ formation in a variety of interstellar environments.
△ Less
Submitted 28 November, 2017;
originally announced November 2017.
-
Chemical Evolution of Turbulent Multiphase Molecular Clouds
Authors:
Valeska Valdivia,
Patrick Hennebelle,
Benjamin Godard,
Maryvonne Gerin,
Pierre Lesaffre,
Jacques Le Bourlot
Abstract:
Molecular clouds are essentially made up of atomic and molecular hydrogen, which in spite of being the simplest molecule in the ISM plays a key role in the chemical evolution of molecular clouds. Since its formation time is very long, the H2 molecules can be transported by the turbulent motions within the cloud toward low density and warm regions, where its enhanced abundance can boost the abundan…
▽ More
Molecular clouds are essentially made up of atomic and molecular hydrogen, which in spite of being the simplest molecule in the ISM plays a key role in the chemical evolution of molecular clouds. Since its formation time is very long, the H2 molecules can be transported by the turbulent motions within the cloud toward low density and warm regions, where its enhanced abundance can boost the abundances of molecules with high endothermicities. We present high resolution simulations where we include the evolution of the molecular gas under the effect of the dynamics, and we analyze its impact on the abundance of CH+.
△ Less
Submitted 6 November, 2017;
originally announced November 2017.
-
Origin of CH+ in diffuse molecular clouds warm H2 and ion-neutral drift
Authors:
Valeska Valdivia,
Benjamin Godard,
Patrick Hennebelle,
Maryvonne Gerin,
Pierre Lesaffre,
Jacques Le Bourlot
Abstract:
This paper assesses the roles of the presence of warm H2, and the increased formation rate due to the ion-neutral drift. We performed ideal MHD simulations that include the heating and cooling of the multiphase ISM, and where we treat dynamically the formation of H2. In a post-processing step we compute the abundances of species at chemical equilibrium. We show that CH+ is efficiently formed at th…
▽ More
This paper assesses the roles of the presence of warm H2, and the increased formation rate due to the ion-neutral drift. We performed ideal MHD simulations that include the heating and cooling of the multiphase ISM, and where we treat dynamically the formation of H2. In a post-processing step we compute the abundances of species at chemical equilibrium. We show that CH+ is efficiently formed at the edge of clumps, in regions where the H2 fraction is low, but nevertheless higher than its equilibrium value, and where the gas temperature is high. We show that warm and out of equilibrium H2 increases the integrated column densities of CH+ by one order of magnitude, up to values still 3-10 times lower than those observed in the diffuse ISM. We balance the Lorentz force with the ion-neutral drag to estimate the ion-drift velocities (vd). We find that the vd distribution peaks around 0.04 km s-1, and that high vd are too rare to have a significant statistical impact on the abundances of CH+. Compared to previous works, our multiphase simulations reduce the spread in vd, and our self-consistent treatment of the ionisation leads to much reduced vd. Nevertheless, our resolution study shows that this velocity distribution is not converged: the ion-neutral drift has a higher impact on CH+ at higher resolution. On the other hand, our ideal MHD simulations do not include ambipolar diffusion, which would yield lower drift velocities. Within these limitations, we conclude that warm H2 is a key ingredient in the efficient formation of CH+ and that the ambipolar diffusion has very little influence on the abundance of CH+, mainly due to the small drift velocities obtained. However, we point out that small-scale processes and other non-thermal processes not included in our MHD simulation may be of crucial importance, and higher resolution studies with better controlled dissipation processes are needed.
△ Less
Submitted 16 December, 2016; v1 submitted 9 December, 2016;
originally announced December 2016.
-
H2 distribution during formation of multiphase molecular clouds
Authors:
Valeska Valdivia,
Patrick Hennebelle,
Maryvonne Gérin,
Pierre Lesaffre
Abstract:
H2 is the simplest and the most abundant molecule in the ISM, and its formation precedes the formation of other molecules. Understanding the dynamical influence of the environment and the interplay between the thermal processes related to the formation and destruction of H2 and the structure of the cloud is mandatory to understand correctly the observations of H2. We perform high resolution MHD co…
▽ More
H2 is the simplest and the most abundant molecule in the ISM, and its formation precedes the formation of other molecules. Understanding the dynamical influence of the environment and the interplay between the thermal processes related to the formation and destruction of H2 and the structure of the cloud is mandatory to understand correctly the observations of H2. We perform high resolution MHD colliding flow simulations with the AMR code RAMSES in which the physics of H2 has been included. We compare the simulation results with various observations including the column densities of excited rotational levels. Due to a combination of thermal pressure, ram pressure and gravity, the clouds produced at the converging point of HI streams are highly inhomogeneous. H2 molecules quickly form in relatively dense clumps and spread into the diffuse interclump gas. This in particular leads to the existence of significant abundances of H2 in the diffuse and warm gas that lies in between clumps. Simulations and observations show similar trends, specially for the HI-to-H2 transition. The abundances of excited rotational levels, calculated at equilibrium in the simulations are very similar to the observed abundances inferred from FUSE results. This is a direct consequence of the presence of the H2 enriched diffuse and warm gas. Our simulations show that H2 rapidly forms in the dense clumps and, due to the complex structure of molecular clouds, quickly spreads at lower densities. Consequently a significant fraction of warm H2 exists in the low density gas. This warm H2 leads to column densities of excited rotational levels close to the observed ones likely revealing the complex intermix between the warm and the cold gas in molecular clouds. This suggests that the 2-phase structure of molecular clouds is an essential ingredient to fully understand molecular hydrogen in these objects.
△ Less
Submitted 22 December, 2015; v1 submitted 17 December, 2015;
originally announced December 2015.
-
H2 distribution during 2-phase Molecular Cloud Formation
Authors:
V. Valdivia,
P. Hennebelle,
M. Gerin,
P. Lesaffre
Abstract:
We performed high-resolution, 3D MHD simulations and we compared to observations of translucent molecular clouds. We show that the observed populations of rotational levels of H2 can arise as a consequence of the multi-phase structure of the ISM.
We performed high-resolution, 3D MHD simulations and we compared to observations of translucent molecular clouds. We show that the observed populations of rotational levels of H2 can arise as a consequence of the multi-phase structure of the ISM.
△ Less
Submitted 12 December, 2015;
originally announced December 2015.
-
A fast tree-based method for estimating column densities in Adaptive Mesh Refinement codes Influence of UV radiation field on the structure of molecular clouds
Authors:
Valeska Valdivia,
Patrick Hennebelle
Abstract:
Context. Ultraviolet radiation plays a crucial role in molecular clouds. Radiation and matter are tightly coupled and their interplay influences the physical and chemical properties of gas. In particular, modeling the radiation propagation requires calculating column densities, which can be numerically expensive in high-resolution multidimensional simulations. Aims. Developing fast methods for est…
▽ More
Context. Ultraviolet radiation plays a crucial role in molecular clouds. Radiation and matter are tightly coupled and their interplay influences the physical and chemical properties of gas. In particular, modeling the radiation propagation requires calculating column densities, which can be numerically expensive in high-resolution multidimensional simulations. Aims. Developing fast methods for estimating column densities is mandatory if we are interested in the dynamical influence of the radiative transfer. In particular, we focus on the effect of the UV screening on the dynamics and on the statistical properties of molecular clouds. Methods. We have developed a tree-based method for a fast estimate of column densities, implemented in the adaptive mesh refinement code RAMSES. We performed numerical simulations using this method in order to analyze the influence of the screening on the clump formation. Results. We find that the accuracy for the extinction of the tree-based method is better than 10%, while the relative error for the column density can be much more. We describe the implementation of a method based on precalculating the geometrical terms that noticeably reduces the calculation time. To study the influence of the screening on the statistical properties of molecular clouds we present the probability distribution function (PDF) of gas and the associated temperature per density bin and the mass spectra for different density thresholds. Conclusions. The tree-based method is fast and accurate enough to be used during numerical simulations since no communication is needed between CPUs when using a fully threaded tree. It is then suitable to parallel computing. We show that the screening for far UV radiation mainly affects the dense gas, thereby favoring low temperatures and affecting the fragmentation.
△ Less
Submitted 2 October, 2014;
originally announced October 2014.