Source author record

Patrick Hennebelle

Patrick Hennebelle appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.SR astro-ph.EP astro-ph.IM astro-ph.CO astro-ph astro-ph.HE physics.comp-ph

Catalog footprint

What is connected

65works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Dissecting the dust distribution and polarization around two B213 young stellar objects with ALMA

The earliest stages of disk formation and dust evolution during the protostellar phase remain poorly constrained. Millimeter dust emission and its polarization provide key insights into the physical processes and material distribution at the envelope-disk interface. We present ALMA polarimetric observations at 1.4 mm and 3 mm of two young stellar objects in Taurus, IRAS 04166+2706 (K04166) and IRAS 04169+2702 (K04169), probing scales from 25 au to 3000 au. We model the Stokes I emission to separate disk and envelope contributions and analyze the polarization properties to identify the dominant polarization mechanisms. K04166 shows extended Stokes I and polarized emission tracing a tentative hourglass magnetic field morphology in its envelope. In the inner envelope and disk (< 100 au), the properties of the polarized emission change, suggesting either a toroidal magnetic field or the presence of large grains. In contrast, K04169 exhibits compact Stokes I and polarized emission consistent with self-scattering from the disk. Both disks are extremely compact, but only K04166 retains a substantial envelope. Our multiscale ALMA polarimetric observations reveal a transition from magnetically aligned grains in envelopes to self-scattering in disks within the transition region of 20-50 au. These results provide important clues on dust grain growth and magnetic field morphology at the disk-envelope scales. Despite being embedded in the same filament, the two sources display striking differences, indicating that K04166 is a young embedded object with a substantial envelope threaded by relatively organized magnetic fields. Meanwhile, K04169 is more evolved, likely to be a young T-Tauri star. However, in both disks, the presence of large grains already suggests a scenario of early dust evolution in disks of the Class 0 stage.

preprint2024arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Influence of protostellar jets and HII regions on the formation and evolution of stellar clusters

Context. Understanding the conditions in which stars and stellar clusters form is of great importance. In particular the role that stellar feedback may have is still hampered by large uncertainties. Aims. We investigate the role played by ionising radiation and protostellar outflows during the formation and evolution of a stellar cluster. To self-consistently take into account gas accretion, we start with clumps of tens of parsecs in size. Methods. Using an adaptive mesh refinement code, we run magneto-hydrodynamical numerical simulations aiming at describing the collapse of massive clumps with either no stellar feedback or taking into account ionising radiation and/or protostellar jets. Results. Stellar feedback substantially modifies the protostellar cluster properties, in several ways. We confirm that protostellar outflows reduce the star formation rate by a factor of a few, although the outflows do not stop accretion and likely enough do not modify the final cluster mass. On the other hand, ionising radiation, once sufficiently massive stars have formed, efficiently expels the remaining gas and reduces the final cluster mass by a factor of several. We found that while HII radiation and jets barely change the distribution of high density gas, the latter increases, at a few places, the dense gas velocity dispersion again by a factor of several. As we are starting from a relatively large scale, we found that the clusters whose mass and size are respectively on the order of a few 1000 M and a fraction of parsec, present a significant level of rotation. Moreover we found that the sink particles which mimic the stars themselves, tend to have rotation axis aligned with the cluster large scale rotation. Finally, computing the classical Q parameter used to quantify stellar cluster structure, we infer that when jets are included in the calculation, [...]

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2021arXiv

An observational correlation between magnetic field, angular momentum and fragmentation in the envelopes of Class 0 protostars?

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.

preprint2021arXiv

Discs and outflows in the early phases of massive star formation: influence of magnetic fields and ambipolar diffusion

We study mass accretion and ejection in the vicinity of massive star forming cores using high-resolution (5 au) 3D AMR numerical simulations. We investigate the mechanisms at the origin of outflows and characterise the properties of the disc forming around massive protostars. We include both protostellar radiative feedback via PMS evolutionary tracks and magnetic ambipolar diffusion. We studied 3 different cases: purely hydrodynamical, ideal MHD, and ambipolar diffusion. In the resistive models, we investigate the effects the initial amplitude of both magnetic field and rotation have on the properties of the massive protostellar system. We use simple criteria to identify the outflow and disc material and follow their evolution as the central star accretes mass up to 20 solar mass. The outflow is completely different when magnetic fields are introduced, so that magnetic processes are the main driver of the outflow up to stellar masses of ~20 solar mass. The disc properties depend on the physics included. The disc formed in the ideal and resistive runs show opposite properties in terms of plasma beta and of magnetic fields topology. While the disc in the ideal case is dominated by the magnetic pressure and the toroidal magnetic fields, the one formed in the resistive runs is dominated by the thermal pressure and has essentially vertical magnetic fields in the inner regions (R<200 au). We find that magnetic processes dominate the early evolution of massive protostellar systems (<20 solar mass) and shapes the accretion/ejection as well as the disc formation. Ambipolar diffusion is mainly at work at disc scales and regulates its properties. Our finding for the outflow and disc properties are reminiscent of low-mass star formation, suggesting that accretion and ejection in young massive and low-mass protostars are regulated by the same physical processes at the early stages.

preprint2020arXiv

Formation and Evolution of Disks around Young Stellar Objects

Recent observations have suggested that circumstellar disks may commonly form around young stellar objects. Although the formation of circumstellar disks can be a natural result of the conservation of angular momentum in the parent cloud, theoretical studies instead show disk formation to be difficult from dense molecular cores magnetized to a realistic level, owing to efficient magnetic braking that transports a large fraction of the angular momentum away from the circumstellar region. We review recent progress in the formation and early evolution of disks around young stellar objects of both low-mass and high-mass, with an emphasis on mechanisms that may bridge the gap between observation and theory, including non-ideal MHD effects and asymmetric perturbations in the collapsing core (e.g., magnetic field misalignment and turbulence). We also address the associated processes of outflow launching and the formation of multiple systems, and discuss possible implications in properties of protoplanetary disks.

preprint2020arXiv

Large-scale turbulent driving regulates star formation in high-redshift gas-rich galaxies

The question of what regulates star formation is a long standing issue. To investigate this issue, we run simulations of a kiloparsec cube section of a galaxy with three kinds of stellar feedback: the formation of HII regions, the explosion of supernovae, and the UV heating. We show that stellar feedback is sufficient to reduce the averaged star formation rate (SFR) to the level of the Schmidt- Kennicutt law in Milky-Way like galaxies but not in high-redshift gas rich galaxies suggesting that another type of support should be added. We investigate whether an external driving of the turbulence such as the one created by the large galactic scales could diminish the SFR at the observed level. Assuming that the Toomre parameter is close to 1 as suggested by the observations, we infer a typical turbulent forcing that we argue should be applied parallel to the plane of the galactic disc. When this forcing is applied in our simulations, the SFR within our simulations closely follows the Schmidt- Kennicutt relation. We found that the velocity dispersion is strongly anisotropic with the velocity dispersion alongside the galactic plane being up to 10 times larger than the perpendicular velocity.

preprint2020arXiv

Physical Processes in Star Formation

Star formation is a complex multi-scale phenomenon that is of significant importance for astrophysics in general. Stars and star formation are key pillars in observational astronomy from local star forming regions in the Milky Way up to high-redshift galaxies. From a theoretical perspective, star formation and feedback processes (radiation, winds, and supernovae) play a pivotal role in advancing our understanding of the physical processes at work, both individually and of their interactions. In this review we will give an overview of the main processes that are important for the understanding of star formation. We start with an observationally motivated view on star formation from a global perspective and outline the general paradigm of the life-cycle of molecular clouds, in which star formation is the key process to close the cycle. After that we focus on the thermal and chemical aspects in star forming regions, discuss turbulence and magnetic fields as well as gravitational forces. Finally, we review the most important stellar feedback mechanisms.

preprint2020arXiv

Protostellar disk formation by a non-rotating, non-axisymmetric collapsing cloud: model and comparison with observations

Planet-forming disks are fundamental objects thought to be inherited from large scale rotation, through the conservation of angular momentum during the collapse of a prestellar dense core. We investigate the possibility for a protostellar disk to be formed from a motionless dense core which contains non-axisymmetric density fluctuations. The rotation is thus generated locally by the asymmetry of the collapse. We study the evolution of the angular momentum in a non-axisymmetric collapse of a dense core from an analytical point of view. To test the theory, we perform three-dimensional simulations of a collapsing prestellar dense core using adaptative mesh refinement. We start from a non-axisymmetrical situation, considering a dense core with random density perturbations that follow a turbulence spectrum. We analyse the emerging disk comparing the angular momentum it contains with the one expected from our analytic development. We study the velocity gradients at different scales in the simulation as it is done with observations. We show that the angular momentum in the frame of a stellar object which is not located at the center of mass of the core is not conserved, due to inertial forces. Our simulations of such non-axisymmetrical collapse quickly produce accretion disks at the small scales in the core. The analysis of the kinematics at different scales in the simulated core reveals projected velocity gradients of amplitudes similar to the ones observed in protostellar cores, and which directions vary, sometimes even reversing when small and large scales are compared. These complex kinematics patterns appear in recent observations, and could be a discriminating feature with models where rotation is inherited from large scales. Our results from simulations without initial rotation are more consistent with these recent observations than when solid-body rotation is initially [abridged]

preprint2020arXiv

The molecular cloud lifecycle

Giant molecular clouds (GMCs) and their stellar offspring are the building blocks of galaxies. The physical characteristics of GMCs and their evolution are tightly connected to galaxy evolution. The macroscopic properties of the interstellar medium propagate into the properties of GMCs condensing out of it, with correlations between e.g. the galactic and GMC scale gas pressures, surface densities and volume densities. That way, the galactic environment sets the initial conditions for star formation within GMCs. After the onset of massive star formation, stellar feedback from e.g. photoionisation, stellar winds, and supernovae eventually contributes to dispersing the parent cloud, depositing energy, momentum and metals into the surrounding medium, thereby changing the properties of galaxies. This cycling of matter between gas and stars, governed by star formation and feedback, is therefore a major driver of galaxy evolution. Much of the recent debate has focused on the durations of the various evolutionary phases that constitute this cycle in galaxies, and what these can teach us about the physical mechanisms driving the cycle. We review results from observational, theoretical, and numerical work to build a dynamical picture of the evolutionary lifecycle of GMC evolution, star formation, and feedback in galaxies.

preprint2020arXiv

The Origin of the Stellar Mass Distribution and Multiplicity

In this chapter, we review some historical understanding and recent advances on the Initial Mass Function (IMF) and the Core Mass Function (CMF), both in terms of observations and theories. We focus mostly on star formation in clustered environment since this is suggested by observations to be the dominant mode of star formation. The statistical properties and the fragmentation behaviour of turbulent gas is discussed, and we also discuss the formation of binaries and small multiple systems.

preprint2020arXiv

What determines the formation and characteristics of protoplanetary discs?

Planets form in protoplanetary discs. Their masses, distribution, and orbits sensitively depend on the structure of the protoplanetary discs. However, what sets the initial structure of the discs in terms of mass, radius and accretion rate is still unknown. We perform non-ideal MHD numerical simulations using the adaptive mesh refinement code Ramses, of a collapsing, one solar mass, molecular core to study the disc formation and early, up to 100 kyr, evolution, paying great attention to the impact of numerical resolution and accretion scheme. We found that while the mass of the central object is almost independent of the numerical parameters such as the resolution and the accretion scheme onto the sink particle, the disc mass, and to a lower extent its size, heavily depend on the accretion scheme, which we found, is itself resolution dependent. This implies that the accretion onto the star and through the disc are largely decoupled. For a relatively large domain of initial conditions (except at low magnetisation), we found that the properties of the disc do not change too significantly. In particular both the level of initial rotation and turbulence do not influence the disc properties provide the core is sufficiently magnetized. After a short relaxation phase, the disc settles in a stationary state. It then slowly grows in size but not in mass. The disc itself is weakly magnetized but its immediate surrounding is on the contrary highly magnetized. Our results show that the disc properties directly depend on the inner boundary condition, i.e. the accretion scheme onto the central object, suggesting that the disc mass is eventually controlled by the small scale accretion process, possibly the star-disc interaction. Because of ambipolar diffusion and its significant resistivity, the disc diversity remains limited and except for low magnetisation, their properties are (abridged).

preprint2016arXiv

A chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core collapse calculations

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.

preprint2016arXiv

Feedback in Clouds II: UV Photoionisation and the first supernova in a massive cloud

Molecular cloud structure is regulated by stellar feedback in various forms. Two of the most important feedback processes are UV photoionisation and supernovae from massive stars. However, the precise response of the cloud to these processes, and the interaction between them, remains an open question. In particular, we wish to know under which conditions the cloud can be dispersed by feedback, which in turn can give us hints as to how feedback regulates the star formation inside the cloud. We perform a suite of radiative magnetohydrodynamic simulations of a 10^5 solar mass cloud with embedded sources of ionising radiation and supernovae, including multiple supernovae and a hypernova model. A UV source corresponding to 10% of the mass of the cloud is required to disperse the cloud, suggesting that the star formation efficiency should be on the order of 10%. A single supernova is unable to significantly affect the evolution of the cloud. However, energetic hypernovae and multiple supernovae are able to add significant quantities of momentum to the cloud, approximately 10^{43} g cm/s of momentum per 10^{51} ergs of supernova energy. This is on the lower range of estimates in other works, since dense gas clumps that remain embedded inside the HII region cause rapid cooling in the supernova blast. We argue that supernovae alone are unable to regulate star formation in molecular clouds, and that strong pre-supernova feedback is required to allow supernova blastwaves to propagate efficiently into the interstellar medium

preprint2016arXiv

Formation of proto-cluster: a virialized structure from gravo-turbulent collapse I. Simulation of cluster formation in collapsing molecular cloud

Stars are often observed to form in clusters. It is therefore important to understand how such a region of concentrated mass is assembled out of the diffuse medium and its properties eventually prescribe the important physical mechanisms and determine the characteristics of the stellar cluster. We study the formation of a gaseous proto-cluster inside a molecular cloud by performing high resolution MHD simulations and associate its internal properties to those of the parent cloud by varying the level of the initial turbulence of the cloud, with a view to better characterize the subsequent stellar cluster formation. The gaseous proto-cluster is formed out of global collapse of a molecular cloud, and has non-negligible rotation due to angular momentum conservation during the collapse of the object. Most of the star formation occurs in this region which occupies only a small volume fraction of the whole cloud. We identify such regions in simulations and compare the gas and sink particles to observations. The gaseous proto-cluster inferred from simulation results present a mass-size relation that is compatible with observations. We stress that the stellar cluster radius, although clearly correlated with the gas cluster radius, depends sensitively on its definition. Energy analysis is performed to confirm that the gaseous proto-cluster is a product of gravo-turbulent reprocessing and that the support of turbulent and rotational energy against self-gravity yields a state of global virial equilibrium although collapse is occurring at smaller scale and the cluster is forming stars actively. This object then serves as the antecedent of the stellar cluster, to which the energy properties are passed on.

preprint2016arXiv

Formation of proto-cluster: a virialized structure from gravo-turbulent collapse II. A two-dimensional analytical model for rotating and accreting system

Most stars are born in the gaseous proto-cluster environment. The knowledge of this intermediate stage gives more accurate constraints on star formation characteristics. We demonstrate that a virialized globally supported structure, in which star formation happens, is formed out of a collapsing molecular cloud, and derive a mapping from the parent cloud parameters to the proto-cluster to predict its properties, with a view to confront analytical calculations with observations and simulations. The virial theorem is decomposed into two dimensions to account for the rotation and the flattened geometry. Equilibrium is found by balancing rotation, turbulence and self-gravity, while turbulence is maintained by accretion driving and dissipates in one crossing time. The angular momentum and the accretion rate of the proto-cluster are estimated from the parent cloud properties. The two-dimensional virial model predicts the size and velocity dispersion given the mass of the proto-cluster and that of the parent cloud. The gaseous proto-clusters lie on a sequence of equilibrium with the trend $R \sim M^{0.5}$, with limited variations depending on the evolutionary stage, the parent cloud, and the parameters not well known like turbulence driving efficiency by accretion and the turbulence anisotropy. The model reproduces successfully observations and simulation results. The properties of proto-clusters follow universal relations and they can be derived from that of the parent cloud. Using simple estimates to infer the peak position of the core mass function (CMF) we find a weak dependence on the cluster mass suggesting that the physical conditions inside proto-clusters may contribute to set a CMF, and by extension an IMF, that looks independent of the environment.

preprint2016arXiv

Photoionisation Feedback in a Self-Gravitating, Magnetised, Turbulent Cloud

We present a new set of analytic models for the expansion of HII regions powered by UV photoionisation from massive stars and compare them to a new suite of radiative magnetohydrodynamic simulations of turbulent, self-gravitating molecular clouds. To perform these simulations we use the Eulerian adaptive mesh magnetohydrodynamics code RAMSES-RT, including radiative transfer of UV photons. Our analytic models successfully predict the global behaviour of the HII region provided the density and velocity structure of the cloud is known. We give estimates for the HII region behaviour based on a power law fit to the density field assuming that the system is virialised. We give a radius at which the ionisation front should stop expanding ("stall"). If this radius is smaller than the distance to the edge of the cloud, the HII region will be trapped by the cloud. This effect is more severe in collapsing clouds than in virialised clouds, since the density in the former increases dramatically over time, with much larger photon emission rates needed for the HII region to escape a collapsing cloud. We also measure the response of Jeans unstable gas to the HII regions to predict the impact of UV radiation on star formation in the cloud. We find that the mass in unstable gas can be explained by a model in which the clouds are evaporated by UV photons, suggesting that the net feedback on star formation should be negative

preprint2016arXiv

Spiral-driven accretion in protoplanetary discs - II Self-similar solutions

Accretion discs are ubiquitous in the universe and it is a crucial issue to understand how angular momentum and mass are being radially transported in these objects. Here, we study the role played by non-linear spiral patterns within hydrodynamical and non self-gravitating accretion disc assuming that external disturbances such as infall onto the disc may trigger them. To do so, we computed self-similar solutions that describe discs in which a spiral wave propagates. Such solutions present both shocks and critical sonic points that we carefully analyze. For all allowed temperatures and for several spiral shocks, we calculated the wave structure. In particular we inferred the angle of the spiral patern, the stress it exerts on the disc as well as the associated flux of mass and angular momentum as a function of temperature. We quantified the rate of angular momentum transport by means of the dimensionless $α$ parameter. For the thickest disc we considered (corresponding to $h/r$ values of about 1/3), we found values of $α$ as high as $0.1$, and scaling with the temperature $T$ such that $α\propto T^{3/2} \propto (h/r)^3$. The spiral angle scales with the temperature as $\arctan(r/h)$. The existence of these solutions suggests that perturbations occurring at disc outer boundaries, such as for example perturbations due to infall motions, can propagate deep inside the disc and therefore should not be ignored, even when considering small radii.

preprint2016arXiv

The effect of ambipolar diffusion on low-density molecular ISM filaments

The filamentary structure of the molecular interstellar medium and the potential link of this morphology to star formation have been brought into focus recently by high resolution observational surveys. An especially puzzling matter is that local interstellar filaments appear to have the same thickness, independent of their column density. This requires a theoretical understanding of their formation process and the physics that governs their evolution. In this work we explore a scenario in which filaments are dissipative structures of the large-scale interstellar turbulence cascade and ion-neutral friction (also called ambipolar diffusion) is affecting their sizes by preventing small-scale compressions. We employ high-resolution, 3D MHD simulations, performed with the grid code RAMSES, to investigate non-ideal MHD turbulence as a filament formation mechanism. We focus the analysis on the mass and thickness distributions of the resulting filamentary structures. Simulations of both driven and decaying MHD turbulence show that the morphologies of the density and the magnetic field are different when ambipolar diffusion is included in the models. In particular, the densest structures are broader and more massive as an effect of ion-neutral friction and the power spectra of both the velocity and the density steepen at a smaller wavenumber. The comparison between ideal and non-ideal MHD simulations shows that ambipolar diffusion causes a shift of the filament thickness distribution towards higher values. However, none of the distributions exhibit the pronounced peak found in the observed local filaments. Limitations in dynamical range and the absence of self-gravity in these numerical experiments do not allow us to conclude at this time whether this is due to the different filament selection or due to the physics inherent of the filament formation.

preprint2015arXiv

Ambipolar diffusion in low-mass star formation. I. General comparison with the ideal MHD case

In this paper, we provide a more accurate description of the evolution of the magnetic flux redistribution during prestellar core collapse by including resistive terms in the magnetohydrodynamics (MHD) equations. We focus more particularly on the impact of ambipolar diffusion. We use the adaptive mesh refinement code RAMSES to carry out such calculations. The resistivities required to calculate the ambipolar diffusion terms were computed using a reduced chemical network of charged, neutral and grain species. The inclusion of ambipolar diffusion leads to the formation of a magnetic diffusion barrier in the vicinity of the core, preventing accumulation of magnetic flux in and around the core and amplification of the field above 0.1G. The mass and radius of the first Larson core remain similar between ideal and non-ideal MHD models. This diffusion plateau has crucial consequences on magnetic braking processes, allowing the formation of disk structures. Magnetically supported outflows launched in ideal MHD models are weakened when using non-ideal MHD. Contrary to ideal MHD misalignment between the initial rotation axis and the magnetic field direction does not significantly affect the results for a given mu, showing that the physical dissipation truly dominate over numerical diffusion. We demonstrate severe limits of the ideal MHD formalism, which yield unphysical behaviours in the long-term evolution of the system. This includes counter rotation inside the outflow, interchange instabilities, and flux redistribution triggered by numerical diffusion, none observed in non-ideal MHD. Disks with Keplerian velocity profiles form in all our non-ideal MHD simulations, with final mass and size which depend on the initial magnetisation. This ranges from a few 0.01 solar masses and 20-30 au for the most magnetised case (mu=2) to 0.2 solar masses and 40-80 au for a lower magnetisation (mu=5).

preprint2015arXiv

Cosmic-ray acceleration in young protostars

The main signature of the interaction between cosmic rays and molecular clouds is the high ionisation degree. This decreases towards the densest parts of a cloud, where star formation is expected, because of energy losses and magnetic effects. However recent observations hint to high levels of ionisation in protostellar systems, therefore leading to an apparent contradiction that could be explained by the presence of energetic particles accelerated within young protostars. Our modelling consists of a set of conditions that has to be satisfied in order to have an efficient particle acceleration through the diffusive shock acceleration mechanism. We find that jet shocks can be strong accelerators of protons which can be boosted up to relativistic energies. Another possibly efficient acceleration site is located at protostellar surfaces, where shocks caused by impacting material during the collapse phase are strong enough to accelerate protons. Our results demonstrate the possibility of accelerating particles during the early phase of a proto-Solar-like system and can be used as an argument to support available observations. The existence of an internal source of energetic particles can have a strong and unforeseen impact on the star and planet formation process as well as on the formation of pre-biotic molecules.

preprint2015arXiv

Cosmic-ray propagation at small scale: a support for protostellar disc formation

As long as magnetic fields remain frozen into the gas, the magnetic braking prevents the formation of protostellar discs. This condition is subordinate to the ionisation fraction characterising the inmost parts of a collapsing cloud. The ionisation level is established by the number and the energy of the cosmic rays able to reach these regions. Adopting the method developed in our previous studies, we computed how cosmic rays are attenuated as a function of column density and magnetic field strength. We applied our formalism to low- and high-mass star formation models obtained by numerical simulations of gravitational collapse that include rotation and turbulence. In general, we found that the decoupling between gas and magnetic fields, condition allowing the collapse to go ahead, occurs only when the cosmic-ray attenuation is taken into account with respect to a calculation in which the cosmic-ray ionisation rate is kept constant. We also found that the extent of the decoupling zone also depends on the dust grain size distribution and is larger if large grains (of radius about 0.1 microns) are formed by compression and coagulation during cloud collapse. The decoupling region disappears for the high-mass case due to magnetic field diffusion that is caused by turbulence and that is not included in the low-mass models. We infer that a simultaneous study of the cosmic-ray propagation during the cloud's collapse may lead to values of the gas resistivity in the innermost few hundred AU around a forming protostar that is higher than generally assumed.

preprint2015arXiv

H2 distribution during formation of multiphase molecular clouds

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.

preprint2015arXiv

Spiral-driven accretion in protoplanetary discs - I. 2D models

We numerically investigate the dynamics of a 2D non-magnetised protoplanetary disc surrounded by an inflow coming from an external envelope. We find that the accretion shock between the disc and the inflow is unstable, leading to the generation of large-amplitude spiral density waves. These spiral waves propagate over long distances, down to radii at least ten times smaller than the accretion shock radius. We measure spiral-driven outward angular momentum transport with 1e-4 < alpha < 1e-2 for an inflow accretion rate Mout>1e-8 Msun/yr. We conclude that the interaction of the disc with its envelope leads to long-lived spiral density waves and radial angular momentum transport with rates that cannot be neglected in young non-magnetised protostellar discs.

preprint2015arXiv

The 21-SPONGE HI Absorption Survey I: Techniques and Initial Results

We present methods and results from "21-cm Spectral Line Observations of Neutral Gas with the EVLA" (21-SPONGE), a large survey for Galactic neutral hydrogen (HI) absorption with the Karl G. Jansky Very Large Array (VLA). With the upgraded capabilities of the VLA, we reach median root-mean-square (RMS) noise in optical depth of $σ_τ=9\times 10^{-4}$ per $0.42\rm\,km\,s^{-1}$ channel for the 31 sources presented here. Upon completion, 21-SPONGE will be the largest HI absorption survey with this high sensitivity. We discuss the observations and data reduction strategies, as well as line fitting techniques. We prove that the VLA bandpass is stable enough to detect broad, shallow lines associated with warm HI, and show that bandpass observations can be combined in time to reduce spectral noise. In combination with matching HI emission profiles from the Arecibo Observatory ($\sim3.5'$ angular resolution), we estimate excitation (or spin) temperatures ($\rm T_s$) and column densities for Gaussian components fitted to sightlines along which we detect HI absorption (30/31). We measure temperatures up to $\rm T_s\sim1500\rm\,K$ for individual lines, showing that we can probe the thermally unstable interstellar medium (ISM) directly. However, we detect fewer of these thermally unstable components than expected from previous observational studies. We probe a wide range in column density between $\sim10^{16}$ and $>10^{21}\rm\,cm^{-2}$ for individual HI clouds. In addition, we reproduce the trend between cold gas fraction and average $\rm T_s$ found by synthetic observations of a hydrodynamic ISM simulation by Kim et al. (2014). Finally, we investigate methods for estimating HI $\rm T_s$ and discuss their biases.

preprint2014arXiv

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

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.

preprint2014arXiv

Autonomous Gaussian Decomposition

We present a new algorithm, named Autonomous Gaussian Decomposition (AGD), for automatically decomposing spectra into Gaussian components. AGD uses derivative spectroscopy and machine learning to provide optimized guesses for the number of Gaussian components in the data, and also their locations, widths, and amplitudes. We test AGD and find that it produces results comparable to human-derived solutions on 21cm absorption spectra from the 21cm SPectral line Observations of Neutral Gas with the EVLA (21-SPONGE) survey. We use AGD with Monte Carlo methods to derive the HI line completeness as a function of peak optical depth and velocity width for the 21-SPONGE data, and also show that the results of AGD are stable against varying observational noise intensity. The autonomy and computational efficiency of the method over traditional manual Gaussian fits allow for truly unbiased comparisons between observations and simulations, and for the ability to scale up and interpret the very large data volumes from the upcoming Square Kilometer Array and pathfinder telescopes.

preprint2014arXiv

Cosmic rays as regulators of molecular cloud properties

Cosmic rays are the main agents in controlling the chemical evolution and setting the ambipolar diffusion time of a molecular cloud. We summarise the processes causing the energy degradation of cosmic rays due to their interaction with molecular hydrogen, focusing on the magnetic effects that influence their propagation. Making use of magnetic field configurations generated by numerical simulations, we show that the increase of the field line density in the collapse region results in a reduction of the cosmic-ray ionisation rate. As a consequence the ionisation fraction decreases, facilitating the decoupling between the gas and the magnetic field.

preprint2014arXiv

Excitation temperature of the warm neutral medium as a new probe of the Lyman-α radiation field

We use the Karl G. Jansky Very Large Array (VLA) to conduct a high-sensitivity survey of neutral hydrogen (HI) absorption in the Milky Way. In combination with corresponding HI emission spectra obtained mostly with the Arecibo Observatory, we detect a widespread warm neutral medium (WNM) component with excitation temperature <Ts>= 7200 (+1800,-1200) K (68% confidence). This temperature lies above theoretical predictions based on collisional excitation alone, implying that Ly-α scattering, the most probable additional source of excitation, is more important in the interstellar medium (ISM) than previously assumed. Our results demonstrate that HI absorption can be used to constrain the Ly-α radiation field, a critical quantity for studying the energy balance in the ISM and intergalactic medium yet notoriously difficult to model because of its complicated radiative transfer, in and around galaxies nearby and at high redshift.

preprint2014arXiv

Mutual influence of supernovae and molecular clouds

Context. Molecular clouds are known to be turbulent and strongly affected by stellar feedback. Moreover, stellar feedback is believed to drive turbulence at large scales in galaxies. Aims. We study the role played by supernovae in molecular clouds and the influence of the magnetic field on this process. Methods. We performed three-dimensional numerical simulations of supernova explosions, in and near turbulent self-gravitating molecular clouds. In order to study the influence of the magnetic field, we performed both hydrodynamical and magnetohydrodynamical simulations. We also ran a series of simple uniform density medium simulations and developed a simple analytical model. Results. We find that the total amount of momentum that is delivered during supernova explosions typically varies by a factor of about 2, even when the gas density changes by 3 orders of magnitude. However, the amount of momentum delivered to the dense gas varies by almost a factor of 10 if the supernova explodes within or outside the molecular cloud. The magnetic field has little influence on the total amount of momentum injected by the supernova explosions but increases the momentum injected into the dense gas. Conclusions. Supernovae that explode inside molecular clouds remove a significant fraction of the cloud mass. Supernovae that explode outside have a limited influence on the cloud. It is thus essential to know sufficiently well the correlation between supernovae and the surrounding dense material in order to know whether supernovae can regulate star formation effectively.

preprint2014arXiv

Shell instability of a collapsing dense core

Understanding the formation of binary and multiple stellar systems largely comes down to studying the circumstances for the fragmentation of a condensing core during the first stages of the collapse. However, the probability of fragmentation and the number of fragments seem to be determined to a large degree by the initial conditions. In this work we study the fate of the linear perturbations of a homogeneous gas sphere both analytically and numerically. In particular, we investigate the stability of the well-known homologous solution that describes the collapse of a uniform spherical cloud. The difficulty of the mathematical singularity in the perturbation equations is surpassed here by explicitly introducing a weak shock next to the sonic point. In parallel, we perform adaptive mesh refinement (AMR) numerical simulations of the linear stages of the collapse and compared the growth rates obtained by each method. With this combination of analytical and numerical tools, we explore the behavior of both spherically symmetric and non-axisymmetric perturbations. The numerical experiments provide the linear growth rates as a function of the core's initial virial parameter and as a function of the azimuthal wave number of the perturbation. The overlapping regime of the numerical experiments and the analytical predictions is the situation of a cold and large cloud, and in this regime the analytically calculated growth rates agree very well with the ones obtained from the simulations. The use of a weak shock as part of the perturbation allows us to find a physically acceptable solution to the equations for a continuous range of growth rates. The numerical simulations agree very well with the analytical prediction for the most unstable cores, while they impose a limit of a virial parameter of 0.1 for core fragmentation in the absence of rotation.

preprint2014arXiv

Simulations of magnetized multiphase galactic disk regulated by supernovae explosions

What exactly controls star formation in the Galaxy remains controversial. In particular, the role of feedback and magnetic field are still partially understood. We investigate the role played by supernovae feedback and magnetic field onto the star formation and the structure of the Galactic disk. We perform numerical simulations of the turbulent, magnetized, self-gravitating, multi-phase, supernovae regulated ISM within a 1 kpc stratified box. We implemented various schemes for the supernovae. This goes from a random distribution at a fixed rate to distributions for which the supernovae are spatially and temporally correlated to the formation of stars. To study the influence of magnetic field on star formation, we perform both hydrodynamical and magneto-hydrodynamical simulations. We find that supernovae feedback has a drastic influence on the galactic evolution. The supernovae distribution is playing a very significant role. When the supernovae are not correlated with star formation events, they do not modify significantly the very high star formation rate obtained without feedback. When the supernovae follow the accretion, the star formation rate can be reduced by a factor up to 30. Magnetic field is also playing a significant role. It reduces the star formation rate by a factor up to 2-3 and reduces the number of collapse sites by a factor of about 2. The exact correlation between the supernovae and the dense gas appears to have significant consequences on the galactic disk evolution and the star formation. This implies that small scale studies are necessary to understand and quantify the feedback efficiency. Magnetic field does influence the star formation at galactic scales by reducing the star formation rate and the number of star formation sites.

preprint2014arXiv

The Origin and Universality of the Stellar Initial Mass Function

We review current theories for the origin of the Stellar Initial Mass Function (IMF) with particular focus on the extent to which the IMF can be considered universal across various environments. To place the issue in an observational context, we summarize the techniques used to determine the IMF for different stellar populations, the uncertainties affecting the results, and the evidence for systematic departures from universality under extreme circumstances. We next consider theories for the formation of prestellar cores by turbulent fragmentation and the possible impact of various thermal, hydrodynamic and magneto-hydrodynamic instabilities. We address the conversion of prestellar cores into stars and evaluate the roles played by different processes: competitive accretion, dynamical fragmentation, ejection and starvation, filament fragmentation and filamentary accretion flows, disk formation and fragmentation, critical scales imposed by thermodynamics, and magnetic braking. We present explanations for the characteristic shapes of the Present-Day Prestellar Core Mass Function and the IMF and consider what significance can be attached to their apparent similarity. Substantial computational advances have occurred in recent years, and we review the numerical simulations that have been performed to predict the IMF directly and discuss the influence of dynamics, time-dependent phenomena, and initial conditions.

preprint2014arXiv

The role of cosmic rays on magnetic field diffusion and the formation of protostellar discs

The formation of protostellar discs is severely hampered by magnetic braking, as long as magnetic fields remain frozen in the gas. The latter condition depends on the levels of ionisation that characterise the innermost regions of a collapsing cloud. The chemistry of dense cloud cores and, in particular, the ionisation fraction is largely controlled by cosmic rays. The aim of this paper is to evaluate whether the attenuation of the flux of cosmic rays expected in the regions around a forming protostar is sufficient to decouple the field from the gas, thereby influencing the formation of centrifugally supported disc. We adopted the method developed in a former study to compute the attenuation of the cosmic-ray flux as a function of the column density and the field strength in clouds threaded by poloidal and toroidal magnetic fields. We applied this formalism to models of low- and high-mass star formation extracted from numerical simulations of gravitational collapse that include rotation and turbulence. For each model we determine the size of the magnetic decoupling zone, where collapse or rotation motion becomes unaffected by the local magnetic field. In general, we find that decoupling only occurs when the attenuation of cosmic rays is taken into account with respect to a calculation in which the cosmic-ray ionisation rate is kept constant. The extent of the decoupling zone also depends on the dust grain size distribution and is larger if large grains (of radius $\sim 10^{-5}$ cm) are formed by compression and coagulation during cloud collapse. We conclude that a realistic treatment of cosmic-ray propagation and attenuation during cloud collapse may lead to a value of the resistivity of the gas in the innermost few hundred AU around a forming protostar that is higher than generally assumed.

preprint2014arXiv

The W43-MM1 mini-starburst ridge, a test for star formation efficiency models

Context: Star formation efficiency (SFE) theories are currently based on statistical distributions of turbulent cloud structures and a simple model of star formation from cores. They remain poorly tested, especially at the highest densities. Aims: We investigate the effects of gas density on the SFE through measurements of the core formation efficiency (CFE). With a total mass of $\sim2\times10^4$ M$_\odot$, the W43-MM1 ridge is one of the most convincing candidate precursor of starburst clusters and thus one of the best place to investigate star formation. Methods: We used high-angular resolution maps obtained at 3 mm and 1 mm within W43-MM1 with the IRAM Plateau de Bure Interferometer to reveal a cluster of 11 massive dense cores (MDCs), and, one of the most massive protostellar cores known. An Herschel column density image provided the mass distribution of the cloud gas. We then measured the 'instantaneous' CFE and estimated the SFE and the star formation rate (SFR) within subregions of the W43-MM1 ridge. Results: The high SFE found in the ridge ($\sim$6% enclosed in $\sim$8 pc$^3$) confirms its ability to form a starburst cluster. There is however a clear lack of dense cores in the northern part of the ridge, which may be currently assembling. The CFE and the SFE are observed to increase with volume gas density while the SFR steeply decreases with the virial parameter, $α_{vir}$. Statistical models of the SFR may well describe the outskirts of the W43-MM1 ridge but struggle to reproduce its inner part, which corresponds to measurements at low $α_{vir}$. It may be that ridges do not follow the log-normal density distribution, Larson relations, and stationary conditions forced in the statistical SFR models.

preprint2013arXiv

Analytical theory for the initial mass function: III time dependence and star formation rate

The present paper extends our previous theory of the stellar initial mass function (IMF) by including the time-dependence, and by including the impact of magnetic field. The predicted mass spectra are similar to the time independent ones with slightly shallower slopes at large masses and peak locations shifted toward smaller masses by a factor of a few. Assuming that star-forming clumps follow Larson type relations, we obtain core mass functions in good agreement with the observationally derived IMF, in particular when taking into account the thermodynamics of the gas. The time-dependent theory directly yields an analytical expression for the star formation rate (SFR) at cloud scales. The SFR values agree well with the observational determinations of various Galactic molecular clouds. Furthermore, we show that the SFR does not simply depend linearly on density, as sometimes claimed in the literature, but depends also strongly on the clump mass/size, which yields the observed scatter. We stress, however, that {\it any} SFR theory depends, explicitly or implicitly, on very uncertain assumptions like clump boundaries or the mass of the most massive stars that can form in a given clump, making the final determinations uncertain by a factor of a few. Finally, we derive a fully time-dependent model for the IMF by considering a clump, or a distribution of clumps accreting at a constant rate and thus whose physical properties evolve with time. In spite of its simplicity, this model reproduces reasonably well various features observed in numerical simulations of converging flows. Based on this general theory, we present a paradigm for star formation and the IMF.

preprint2013arXiv

Cosmic-ray ionisation in collapsing clouds

Cosmic rays (CR) play an important role in dense molecular cores, affecting their thermal and dynamical evolution and initiating the chemistry. Several studies have shown that the formation of protostellar discs in collapsing clouds is severely hampered by the braking torque exerted by the entrained magnetic field on the infalling gas, as long as the field remains frozen to the gas. We examine the possibility that the concentration and twisting of the field lines in the inner region of collapse can produce a significant reduction of the ionisation fraction. To check whether the CR ionisation rate (CRir) can fall below the critical value required to maintain good coupling, we first study the propagation of CRs in a model of a static magnetised cloud varying the relative strength of the toroidal/poloidal components and the mass-to-flux ratio. We then follow the path of CRs using realistic magnetic field configurations generated by numerical simulations of a rotating collapsing core. We find that an increment of the toroidal component of the magnetic field, or, in general, a more twisted configuration of the field lines, results in a decrease in the CR flux. This is mainly due to the magnetic mirroring effect that is stronger where larger variations in the field direction are present. In particular, we find a decrease of the CRir below 10^-18 s-1 in the central 300-400 AU, where density is higher than about 10^9 cm-3. This very low value of the CRir is attained in the cases of intermediate and low magnetisation (mass-to-flux ratio lambda=5 and 17, respectively) and for toroidal fields larger than about 40% of the total field. Magnetic field effects can significantly reduce the ionisation fraction in collapsing clouds. We provide a handy fitting formula to compute approximately the attenuation of the CRir in a molecular cloud as a function of the density and the magnetic configuration.

preprint2013arXiv

Ion-neutral friction and accretion-driven turbulence in self-gravitating filaments

Recent Herschel observations have confirmed that filaments are ubiquitous in molecular clouds and suggest that irrespectively of the column density, there is a characteristic width of about 0.1 pc whose physical origin remains unclear. We develop an analytical model that can be applied to self-gravitating accreting filaments. It is based on one hand on the virial equilibrium of the central part of the filament and on the other hand on energy balance between the turbulence driven by accretion onto the filament and dissipation. We consider two dissipation mechanisms the turbulent cascade and the ion-neutral friction. Our model predicts that the width of the filament inner part is almost independent of the column density and leads to values comparable to what is inferred observationally if dissipation is due to ion-neutral friction. On the contrary turbulent dissipation leads to a structure that is bigger and depends significantly on the column density. Our model provides a reasonable physical explanation which could explain the observed filament width when they are self-gravitating. It predicts the correct order or magnitude though hampered by some uncertainties.

preprint2013arXiv

On the origin of non self-gravitating filaments in the ISM

{Filaments are ubiquitous in the interstellar medium as recently emphasized by Herschel, yet their physical origin remains elusive} {It is therefore important to understand the physics of molecular clouds to investigate how filaments form and what is the role played by various processes such as turbulence and magnetic field.} {We use ideal MHD simulations to study the formation of clumps in various conditions including different magnetization and Mach numbers as well as two completely different setup. We then perform several analysis to compute the shape of the clumps and their link to velocities and forces using various approaches.} {We find that on average, clumps in MHD simulations are more filamentary that clumps in hydrodynamical simulations. Detailed analyses reveal that the filaments are in general preferentially aligned with the strain which means that these structures simply result from the strech induced by turbulence. Moreover filaments tend to be confined by the Lorentz force which therefore lead them to survive longer in magnetized flows. We show that they have in all simulations a typical thickness equal to a few grid cells suggesting that they are primarily associated to the energy dissipation within the flow. We estimate the order of magnitude of the dissipation length associated to the ion-neutral friction and conclude that in well UV shielded regions it is of the order of 0.1 pc and therefore could possibly set the typical size of non self-gravitating filaments.} {Filaments are ubiquitous because they are the results of the very generic turbulent strain and because magnetic field help to keep them coherent. We suggest that energy dissipation is playing a determinant role in their formation.}

preprint2013arXiv

The influence of turbulence during magnetized core collapse and its consequences on low-mass star formation

[Abridged] Theoretical and numerical studies of star formation have shown that magnetic field (B) has a strong influence on both disk formation and fragmentation; even a relatively low B can prevent these processes. However, very few studies investigated the combined effects of B and turbulence. We study the effects of turbulence in magnetized core collapse, focusing on the magnetic diffusion, the orientation of the angular momentum (J) of the protostellar core, and on its consequences on disk formation, fragmentation and outflows. We perform 3D, AMR, MHD simulations of magnetically supercritical collapsing dense cores of 5 Msun using the MHD code RAMSES. A turbulent velocity field is imposed as initial conditions, characterised by a Kolmogorov power spectrum. Different levels of turbulence and magnetization are investigated, as well as 3 realisations for the turbulent velocity field. Magnetic diffusion, orientation of the rotation axis with respect to B, transport of J, disk formation, fragmentation and outflows formation are studied. The turbulent velocity field imposed as initial conditions contains a non-zero J, responsible for a misalignment of the rotation axis. Turbulence is also responsible for an effective turbulent diffusivity in the vicinity of the core. Both effects are responsible for a significant decrease of the magnetic braking, and facilitate the formation of early massive disks for not too high magnetization. Fragmentation can occur even with mu ~ 5 at late time in contrast with 1 Msun cores for which fragmentation is prevented for such values of mu. Slow asymmetric outflows are launched. They carry a mass which is comparable to the mass within the core. Because of misalignment and turbulent diffusion, massive disk formation is possible though their mass and size are still reduced compared to the hydrodynamical case. We find that for mu >= 5, fragmentation can happen.

preprint2012arXiv

Early stages of cluster formation: fragmentation of massive dense cores down to ~1000 AU

In order to study the fragmentation of massive dense cores, which constitute the cluster cradles, we observed with the PdBI in the most extended configuration the continuum at 1.3 mm and the CO(2-1) emission of four massive cores. We detect dust condensations down to ~0.3 Msun and separate millimeter sources down to 0.4" or ~1000 AU, comparable to the sensitivities and separations reached in optical/infrared studies of clusters. The CO(2-1) high angular resolution images reveal high-velocity knots usually aligned with previously known outflow directions. This, in combination with additional cores from the literature observed at similar mass sensitivity and spatial resolution, allowed us to build a sample of 18 protoclusters with luminosities spanning 3 orders of magnitude. Among the 18 regions, ~30% show no signs of fragmentation, while 50% split up into ~4 millimeter sources. We compiled a list of properties for the 18 massive dense cores, such as bolometric luminosity, total mass, and mean density, and found no correlation of any of these parameters with the fragmentation level. In order to investigate the combined effects of magnetic field, radiative feedback and turbulence in the fragmentation process, we compared our observations to radiation magneto-hydrodynamic simulations, and found that the low-fragmented regions are well reproduced in the magnetized core case, while the highly-fragmented regions are consistent with cores where turbulence dominates over the magnetic field. Overall, our study suggests that the fragmentation in massive dense cores could be determined by the initial magnetic field/turbulence balance in each particular core.

preprint2012arXiv

Formation of low-mass stars and brown dwarfs

These lectures attempt to expose the most important ideas, which have been proposed to explain the formation of stars with particular emphasis on the formation of brown dwarfs and low-mass stars. We first describe the important physical processes which trigger the collapse of a self-gravitating piece of fluid and regulate the star formation rate in molecular clouds. Then we review the various theories which have been proposed along the years to explain the origin of the stellar initial mass function paying particular attention to four models, namely the competitive accretion and the theories based respectively on stopped accretion, MHD shocks and turbulent dispersion. As it is yet unsettled whether the brown dwarfs form as low-mass stars, we present the theory of brown dwarfs based on disk fragmentation stressing all the uncertainties due to the radiative feedback and magnetic field. Finally, we describe the results of large scale simulations performed to explain the collapse and fragmentation of molecular clouds.

preprint2012arXiv

Formation of proto-clusters and star formation within clusters: apparent universality of the initial mass function ?

It is believed that the majority of stars form in clusters. Therefore it is likely that the gas physical conditions that prevail in forming clusters, largely determine the properties of stars that form and in particular the initial mass function. We develop an analytical model to account for the formation of low mass clusters and the formation of stars within clusters. The formation of clusters is determined by an accretion rate, the virial equilibrium and energy as well as thermal balance. For this latter both molecular and dust cooling are considered using published rates. The star distribution is computed within the cluster using the physical conditions inferred from this model and the Hennebelle & Chabrier theory. Our model reproduces well the mass-size relation of low mass clusters (up to few $\simeq 10^3$ M$_\odot$ of stars corresponding to about 5 times more gas) and an initial mass function which is $i)$ very close to the Chabrier's IMF, $ii)$ weakly dependent on the mass of the clusters, $iii)$ relatively robust to (i.e. not too steeply dependent on) variations of physical quantities as accretion rate, radiation and cosmic rays abundances. The weak dependence of the mass distribution of stars with the cluster mass results from the compensation between varying clusters densities, velocity dispersions and temperatures all inferred from first physical principles. This constitutes a possible explanation for the apparent universality of the IMF within the Galaxy though variations with the local conditions could certainly be observed.

preprint2012arXiv

Protostellar disk formation and transport of angular momentum during magnetized core collapse

Theoretical studies of collapsing clouds have found that even a relatively weak magnetic field (B) may prevent the formation of disks and their fragmentation. However, most previous studies have been limited to cases where B and the rotation axis of the cloud are aligned. We study the transport of angular momentum, and its effects on disk formation, for non-aligned initial configurations and a range magnetic intensities. We perform 3D AMR MHD simulations of magnetically supercritical collapsing dense cores using the code Ramses. We compute the contributions of the processes transporting angular momentum (J), in the envelope and the region of the disk. We clearly define what could be defined as centrifugally supported disks and study their properties. At variance with earlier analyses, we show that the transport of J acts less efficiently in collapsing cores with non-aligned rotation axis and B. Analytically, this result can be understood by taking into account the bending of field lines occurring during the gravitational collapse. For the transport of J, we conclude that magnetic braking in the mean direction of B tends to dominate over both the gravitational and outflow transport of J. We find that massive disks, containing at least 10% of the initial core mass, can form during the earliest stages of star formation even for mass-to-flux ratios as small as 3 to 5 times the critical value. At higher field intensities, the early formation of massive disks is prevented. Given the ubiquity of Class I disks, and because the early formation of massive disks can take place at moderate magnetic intensities, we speculate that for stronger fields, disks will form later, when most of the envelope will have been accreted. In addition, we speculate that some observed early massive disks may actually be outflow cavities, mistaken for disks by projection effects. (Abridged version of the abstract.)

preprint2012arXiv

Theories of the massive star formation: a (short) review

We briefly review the recent numerical works that have been performed to understand the formation of massive stars. After a brief description of the classical works, we review more specifically $i)$ the problem of building stars more massive than 20 $M_\odot$ and $ii)$ how to prevent the massive cores to fragment in many objects. Multi-D simulations succeed in circumventing the radiative pressure leading to the formation of massive stars although some questions are still debated regarding how is accretion exactly proceeding. While the core fragmentation is slightly reduced by the radiative feedback and the magnetic field when they are treated separately, it is almost entirely suppressed when both of them are included. This is because, magnetic field by removing angular momentum focusses the flow in a compact region. This makes the radiative feedback very efficient leading to a significant increase of the temperature.

preprint2012arXiv

Turbulent molecular clouds

Stars form within molecular clouds but our understanding of this fundamental process remains hampered by the complexity of the physics that drives their evolution. We review our observational and theoretical knowledge of molecular clouds trying to confront the two approaches wherever possible. After a broad presentation of the cold interstellar medium and molecular clouds, we emphasize the dynamical processes with special focus to turbulence and its impact on cloud evolution. We then review our knowledge of the velocity, density and magnetic fields. We end by openings towards new chemistry models and the links between molecular cloud structure and star--formation rates.

preprint2012arXiv

UV-driven chemistry in simulations of the interstellar medium. I. Post-processed chemistry with the Meudon PDR code

Our main purpose is to estimate the effect of assuming uniform density on the line-of-sight in PDR chemistry models, compared to a more realistic distribution for which total gas densities may well vary by several orders of magnitude. A secondary goal of this paper is to estimate the amount of molecular hydrogen which is not properly traced by the CO (J = 1 -> 0) line, the so-called "dark molecular gas". We use results from a magnetohydrodynamical (MHD) simulation as a model for the density structures found in a turbulent diffuse ISM with no star-formation activity. The Meudon PDR code is then applied to a number of lines of sight through this model, to derive their chemical structures. It is found that, compared to the uniform density assumption, maximal chemical abundances for H2, CO, CH and CN are increased by a factor 2 to 4 when taking into account density fluctuations on the line of sight. The correlations between column densities of CO, CH and CN with respect to those of H2 are also found to be in better overall agreement with observations. For instance, at N(H2) > 2.10^{20} cm-2, while observations suggest that d[log N(CO)]=d[log N(H2)] = 3.07 +/- 0.73, we find d[log N(CO)]=d[log N(H2)] =14 when assuming uniform density, and d[log N(CO)]=d[log N(H2)] = 5.2 when including density fluctuations.

preprint2011arXiv

Analytical star formation rate from gravoturbulent fragmentation

We present an analytical determination of the star formation rate (SFR) in molecular clouds, based on a time-dependent extension of our analytical theory of the stellar initial mass function (IMF). The theory yields SFR's in good agreement with observations, suggesting that turbulence {\it is} the dominant, initial process responsible for star formation. In contrast to previous SFR theories, the present one does not invoke an ad-hoc density threshold for star formation; instead, the SFR {\it continuously} increases with gas density, naturally yielding two different characteristic regimes, thus two different slopes in the SFR vs gas density relationship, in agreement with observational determinations. Besides the complete SFR derivation, we also provide a simplified expression, which reproduces reasonably well the complete calculations and can easily be used for quick determinations of SFR's in cloud environments. A key property at the heart of both our complete and simplified theory is that the SFR involves a {\it density-dependent dynamical time}, characteristic of each collapsing (prestellar) overdense region in the cloud, instead of one single mean or critical freefall timescale. Unfortunately, the SFR also depends on some ill determined parameters, such as the core-to-star mass conversion efficiency and the crossing timescale. Although we provide estimates for these parameters, their uncertainty hampers a precise quantitative determination of the SFR, within less than a factor of a few.

preprint2011arXiv

Collapse of massive magnetized dense cores using radiation-magneto-hydrodynamics: early fragmentation inhibition

We report the results of radiation-magneto-hydrodynamics calculations in the context of high mass star formation, using for the first time a self-consistent model for photon emission (i.e. via thermal emission and in radiative shocks) and with the high resolution necessary to resolve properly magnetic braking effects and radiative shocks on scales <100 AU. We investigate the combined effects of magnetic field, turbulence, and radiative transfer on the early phases of the collapse and the fragmentation of massive dense cores. We identify a new mechanism that inhibits initial fragmentation of massive dense cores, where magnetic field and radiative transfer interplay. We show that this interplay becomes stronger as the magnetic field strength increases. Magnetic braking is transporting angular momentum outwards and is lowering the rotational support and is thus increasing the infall velocity. This enhances the radiative feedback owing to the accretion shock on the first core. We speculate that highly magnetized massive dense cores are good candidates for isolated massive star formation, while moderately magnetized massive dense cores are more appropriate to form OB associations or small star clusters.

preprint2011arXiv

Dimensional argument for the impact of turbulent support on the stellar initial mass function

We present a simple dimensional argument to illustrate the impact of nonthermal support from turbulent velocity dispersion on the shape of the prestellar core mass function (CMF), precursor of the stellar initial mass function (IMF). The argument demonstrates the need to invoke such support to recover the Salpeter slope in the high-mass part of the CMF/IMF, whereas pure thermal support leads to a much steeper slope. This simple dimensional argument clearly highlights the results obtained in the complete Hennebelle-Chabrier theory of the IMF.

preprint2011arXiv

Molecular Cloud Evolution IV: Magnetic Fields, Ambipolar Diffusion, and the Star Formation Efficiency

We investigate the formation and evolution of giant molecular clouds (GMCs) by the collision of convergent warm neutral medium (WNM) streams in the interstellar medium, in the presence of magnetic fields and ambipolar diffusion (AD), focusing on the evolution of the star formation rate (SFR) and efficiency (SFE), as well as of the mass-to-magnetic-flux ratio (M2FR) in the forming clouds. We find that: 1) Clouds formed by supercritical inflow streams proceed directly to collapse, while clouds formed by subcritical streams first contract and then re-expand, oscillating on the scale of tens of Myr. 2) Our suite of simulations with initial magnetic field strength of 2, 3, and 4 $μ\G$ show that only supercritical or marginal critical streams lead to reasonable star forming rates. 3) The GMC's M2FR is a generally increasing function of time, whose growth rate depends on the details of how mass is added to the GMC from the WNM. 4) The M2FR is a highly fluctuating function of position in the clouds. 5) In our simulations, the SFE approaches stationarity, because mass is added to the GMC at a similar rate at which it converts mass to stars. In such an approximately stationary regime, the SFE provides a proxy of the supercritical mass fraction in the cloud. 6) We observe the occurrence of buoyancy of the low-M2FR regions within the gravitationally-contracting GMCs, so that the latter naturally segregate into a high-density, high-M2FR "core" and a low-density, low-M2FR "envelope", without the intervention of AD. (Abridged)

preprint2011arXiv

Radiation hydrodynamics with Adaptive Mesh Refinement and application to prestellar core collapse. I Methods

Radiative transfer has a strong impact on the collapse and the fragmentation of prestellar dense cores. We present the radiation-hydrodynamics solver we designed for the RAMSES code. The method is designed for astrophysical purposes, and in particular for protostellar collapse. We present the solver, using the co-moving frame to evaluate the radiative quantities. We use the popular flux limited diffusion approximation, under the grey approximation (one group of photon). The solver is based on the second-order Godunov scheme of RAMSES for its hyperbolic part, and on an implicit scheme for the radiation diffusion and the coupling between radiation and matter. We report in details our methodology to integrate the RHD solver into RAMSES. We test successfully the method against several conventional tests. For validation in 3D, we perform calculations of the collapse of an isolated 1 M_sun prestellar dense core, without rotation. We compare successfully the results with previous studies using different models for radiation and hydrodynamics. We have developed a full radiation hydrodynamics solver in the RAMSES code, that handles adaptive mesh refinement grids. The method is a combination of an explicit scheme and an implicit scheme, accurate to the second-order in space. Our method is well suited for star formation purposes. Results of multidimensional dense core collapse calculations with rotation are presented in a companion paper.

preprint2011arXiv

Theories of the initial mass function

We review the various theories which have been proposed along the years to explain the origin of the stellar initial mass function. We pay particular attention to four models, namely the competitive accretion and the theories based respectively on stopped accretion, MHD shocks and turbulent dispersion. In each case, we derive the main assumptions and calculations that support each theory and stress their respective successes and failures or difficulties.

preprint2010arXiv

Accretion-Driven Turbulence as Universal Process: Galaxies, Molecular Clouds, and Protostellar Disks

Complex turbulent motions are ubiquitously observed in many astrophysical systems. Their origin, however, is still poorly understood. When cosmic structures form, they grow in mass via accretion from the surrounding environment. We propose that this accretion is able to drive internal turbulent motions in a wide range of astrophysical objects and study this process in the case of galaxies, molecular clouds and protoplanetary disks. We use a combination of numerical simulations and analytical arguments to predict the level of turbulence as a function of the accretion rate, the dissipation scale, and the density contrast, and compare with observational data. We find that in Milky Way type galaxies the observed level of turbulence in the interstellar medium can be explained by accretion, provided that the galaxies gain mass at a rate comparable to the rate at which they form stars. This process is particularly relevant in the extended outer disks beyond the star-forming radius. We also calculate the rate at which molecular clouds grow in mass when they build up from the atomic component of the galactic gas and find that their internal turbulence is likely to be driven by accretion as well. It is the very process of cloud formation that excites turbulent motions on small scales by establishing the turbulent cascade. In the case of T Tauri disks, we show that accretion can drive subsonic turbulence at the observed level if the rate at which gas falls onto the disk is comparable to the rate at which disk material accretes onto the central star. This also explains the observed relation of accretion rate and stellar mass, dM/dt ~ M^1.8. The efficiency required to convert infall motion into turbulence is of the order of a few percent in all three cases. We conclude that accretion-driven turbulence is a universal concept with far-reaching implications for a wide range of astrophysical objects.

preprint2010arXiv

Interstellar OH+, H2O+ and H3O+ along the sight-line to G10.6-0.4

We report the detection of absorption lines by the reactive ions OH+, H2O+ and H3O+ along the line of sight to the submillimeter continuum source G10.6$-$0.4 (W31C). We used the Herschel HIFI instrument in dual beam switch mode to observe the ground state rotational transitions of OH+ at 971 GHz, H2O+ at 1115 and 607 GHz, and H3O+ at 984 GHz. The resultant spectra show deep absorption over a broad velocity range that originates in the interstellar matter along the line of sight to G10.6$-$0.4 as well as in the molecular gas directly associated with that source. The OH+ spectrum reaches saturation over most velocities corresponding to the foreground gas, while the opacity of the H2O+ lines remains lower than 1 in the same velocity range, and the H3O+ line shows only weak absorption. For LSR velocities between 7 and 50 kms$^{-1}$ we estimate total column densities of $N$(OH+) $> 2.5 \times 10^{14}$ cm$^{-2}$, $N$(H2O+) $\sim 6 \times 10^{13}$ cm$^{-2}$ and $N$(H3O+) $\sim 4.0 \times 10^{13}$ cm$^{-2}$. These detections confirm the role of O$^+$ and OH$^+$ in initiating the oxygen chemistry in diffuse molecular gas and strengthen our understanding of the gas phase production of water. The high ratio of the OH+ by the H2O+ column density implies that these species predominantly trace low-density gas with a small fraction of hydrogen in molecular form.

preprint2010arXiv

Outflows and mass accretion in collapsing dense cores with misaligned rotation axis and magnetic field

Outflows and jets are intimately related to the formation of stars, and play an important role in redistributing mass, energy and angular momentum within the dense core and parent cloud. The interplay between magnetic field and rotation is responsible for launching these outflows, whose formation has been generally carried out for idealized systems where the angle $α$ between the rotation axis and large-scale magnetic field is zero. Here we explore, through three-dimensional ideal magneto-hydrodynamic simulations, the effects of a non-zero $α$ on the formation of outflows during the collapse of dense pre-stellar cores. We find that mass ejection is less efficient for increasing angle $α$, and that outflows are essentially suppressed for $α\sim90^{\circ}$. An important consequence is a corresponding increase of the mass accreted onto the adiabatic (first) core. In addition, mean flow velocities tend to increase with $α$, and misaligned configurations produce clumpy, heterogeneous outflows that undergo precession, and are more prone to instabilities.

preprint2010arXiv

Radiative, magnetic and numerical feedbacks on small-scale fragmentation

Radiative feedback and magnetic field are understood to have a strong impact on the protostellar collapse. We present high resolution numerical calculations of the collapse of a 1 solar mass dense core in solid body rotation, including both radiative transfer and magnetic field. Using typical parameters for low-mass cores, we study thoroughly the effect of radiative transfer and magnetic field on the first core formation and fragmentation. We show that including the two aforementioned physical processes does not correspond to the simple picture of adding them separately. The interplay between the two is extremely strong, via the magnetic braking and the radiation from the accretion shock.

preprint2010arXiv

Star formation: statistical measure of the correlation between the prestellar core mass function and the stellar initial mass function

We present a simple statistical analysis of recent numerical simulations exploring the correlation between the core mass function obtained from the fragmentation of a molecular cloud and the stellar mass function which forms from these collapsing cores. Our analysis shows that the distributions of bound cores and sink particles obtained in the simulations are consistent with the sinks being formed predominantly from their parent core mass reservoir, with a statistical dispersion of the order of one third of the core mass. Such a characteristic dispersion suggests that the stellar initial mass function is relatively tightly correlated to the parent core mass function, leading to two similar distributions, as observed. This in turn argues in favor of the IMF being essentially determined at the early stages of core formation and being only weakly affected by the various environmental factors beyond the initial core mass reservoir, at least in the mass range explored in the present study. Accordingly, the final IMF of a star forming region should be determined reasonably accurately, statistically speaking, from the initial core mass function, provided some uniform efficiency factor. The calculations also show that these statistical fluctuations, due e.g. to variations among the core properties, broaden the low-mass tail of the IMF compared with the parent CMF, providing an explanation for the fact that this latter appears to underestimate the number of "pre brown dwarf" cores compared with the observationally-derived brown dwarf IMF.

preprint2010arXiv

The Angular Momentum of Magnetized Molecular Cloud Cores: A 2D-3D Comparison

We study the rotational properties of magnetized and self-gravitating molecular cloud cores formed in 2 very high resolution 3D molecular cloud simulations.The simulations have been performed using the code RAMSES at an effective resolution of 4096^3.One simulation represents a mildly magnetically-supercritical cloud and the other a strongly magnetically-supercritical cloud.A noticeable difference between the 2 simulations is the core formation efficiency (CFE) of the high density cores.In the strongly supercritical simulations the CFE is ~33 % per free-fall time of the cloud tff,cl, whereas in the mildly supercritical simulations this value goes down to ~6%/tff,cl. A comparison of the intrinsic specific angular momentum j3D distributions of the cores with the distribitions of j2D derived using synthetic 2D velocity maps of the cores,shows that the synthetic observations tend to overestimate the true value of j by a factor of ~10.The origin of this discrepancy lies in the fact that contrary to the intrinsic determination which sums up the individual gas parcels contributions to j, the determination of j using the observational procedure which is based on a measurement on the global velocity gradient under the hypothesis of uniform rotation smoothes out the complex fluctuations present in the 3D velocity field. Our results provide a natural explanation for the discrepancy by a factor ~10 observed between the intrinsic 3D distributions of j and the corresponding distributions derived in real observations.We suggest that measurements of j which are based on the measurement of the observed global velocity gradients may need to be reduced by a factor of ~10 in order to derive a more accurate estimate of j in the cores.

preprint2010arXiv

The thermally-unstable warm neutral medium: key for modeling the interstellar medium

We present 21-cm absorption measurements towards 12 radio continuum sources with previously identified thermally-unstable warm neutral medium (WNM). These observations were obtained with the Expanded Very Large Array (EVLA) and were complemented with the HI emission spectra obtained with the Arecibo Observatory. Out of 12 sources, HI absorption was detected along 5 lines of sight (seven new absorption features in total), resulting in a detection rate of ~42%. While our observations are sensitive to the WNM with a spin temperature T_s<3000 K, we detected only two wide absorption lines with T_s=400-900 K. These temperatures lie above the range allowed for the cold neutral medium (CNM) by the thermal equilbrium models and signify the thermally unstable WNM. Several absorption features have an optical depth of only a few x10^{-3}. While this is close or lower than what is theoretically expected for the CNM, we show that these weak lines are important for constraining the fraction of the thermally unstable WNM. Our observations demonstrate that, for the first time, high bandpass stability can be achieved with the VLA, allowing detection of absorption lines with a peak optical depth of ~10^{-3}.

preprint2009arXiv

Analytical theory for the initial mass function: II. Properties of the flow

Recently, Hennebelle and Chabrier (2008) derived an analytical theory for the mass spectrum of non self-gravitating clumps associated with overdensities in molecular clouds and for the initial mass function of gravitationally bound prestellar cores, as produced by the turbulent collapse of the cloud. In this companion paper, we examine the effects of the non-isothermality of the flow, of the turbulence forcing and of local fluctuation of the velocity dispersion, on the mass function. In particular, we investigate the influence of a polytropic equation of state and of the effective adiabatic exponent $γ$ and find that it has a drastic influence on the low mass part of the IMF. We also consider a barotropic equation of state (i.e. a piecewise polytropic eos) that mimics the thermal behaviour of the molecular gas and compare the prediction of our theory with the results of numerical simulations and with the observationally-derived IMF, for cloud parameters which satisfy Larson's type relations. We find that for clouds whose density is, at all scales, almost an order of magnitude larger than the density inferred for the CO clumps in the Galaxy, a good agreement is obtained between the theory and the observed IMF, suggesting that star formation preferentially occurs in high density environments. We derive an analytical expression for the IMF which generalizes the expression previously obtained for the isothermal case. This easy-to-implement analytical IMF should serve as a template to compare observational or numerical results with the theory.

preprint2009arXiv

Clump morphology and evolution in MHD simulations of molecular cloud formation

Abridged: We study the properties of clumps formed in three-dimensional weakly magnetized magneto-hydrodynamic simulations of converging flows in the thermally bistable, warm neutral medium (WNM). We find that: (1) Similarly to the situation in the classical two-phase medium, cold, dense clumps form through dynamically-triggered thermal instability in the compressed layer between the convergent flows, and are often characterised by a sharp density jump at their boundaries though not always. (2) However, the clumps are bounded by phase-transition fronts rather than by contact discontinuities, and thus they grow in size and mass mainly by accretion of WNM material through their boundaries. (3) The clump boundaries generally consist of thin layers of thermally unstable gas, but these layers are often widened by the turbulence, and penetrate deep into the clumps. (4) The clumps are approximately in both ram and thermal pressure balance with their surroundings, a condition which causes their internal Mach numbers to be comparable to the bulk Mach number of the colliding WNM flows. (5) The clumps typically have mean temperatures 20 < T < 50 K, corresponding to the wide range of densities they contain (20 < n < 5000 pcc) under a nearly-isothermal equation of state. (6) The turbulent ram pressure fluctuations of the WNM induce density fluctuations that then serve as seeds for local gravitational collapse within the clumps. (7) The velocity and magnetic fields tend to be aligned with each other within the clumps, although both are significantly fluctuating, suggesting that the velocity tends to stretch and align the magnetic field with it. (8) The typical mean field strength in the clumps is a few times larger than that in the WNM. (9) The magnetic field strength has a mean value of B ~ 6 mu G ...

Patrick Hennebelle

What is connected

Connect this record

See the researcher in context

Building this map preview

65 published item(s)

Dissecting the dust distribution and polarization around two B213 young stellar objects with ALMA

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

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

Influence of protostellar jets and HII regions on the formation and evolution of stellar clusters

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

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

An observational correlation between magnetic field, angular momentum and fragmentation in the envelopes of Class 0 protostars?

Discs and outflows in the early phases of massive star formation: influence of magnetic fields and ambipolar diffusion

Formation and Evolution of Disks around Young Stellar Objects

Large-scale turbulent driving regulates star formation in high-redshift gas-rich galaxies

Physical Processes in Star Formation

Protostellar disk formation by a non-rotating, non-axisymmetric collapsing cloud: model and comparison with observations

The molecular cloud lifecycle

The Origin of the Stellar Mass Distribution and Multiplicity

What determines the formation and characteristics of protoplanetary discs?

A chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core collapse calculations

Feedback in Clouds II: UV Photoionisation and the first supernova in a massive cloud

Formation of proto-cluster: a virialized structure from gravo-turbulent collapse I. Simulation of cluster formation in collapsing molecular cloud

Formation of proto-cluster: a virialized structure from gravo-turbulent collapse II. A two-dimensional analytical model for rotating and accreting system

Photoionisation Feedback in a Self-Gravitating, Magnetised, Turbulent Cloud

Spiral-driven accretion in protoplanetary discs - II Self-similar solutions

The effect of ambipolar diffusion on low-density molecular ISM filaments

Ambipolar diffusion in low-mass star formation. I. General comparison with the ideal MHD case

Cosmic-ray acceleration in young protostars

Cosmic-ray propagation at small scale: a support for protostellar disc formation

H2 distribution during formation of multiphase molecular clouds

Spiral-driven accretion in protoplanetary discs - I. 2D models

The 21-SPONGE HI Absorption Survey I: Techniques and Initial Results

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

Autonomous Gaussian Decomposition

Cosmic rays as regulators of molecular cloud properties

Excitation temperature of the warm neutral medium as a new probe of the Lyman-α radiation field

Mutual influence of supernovae and molecular clouds

Shell instability of a collapsing dense core

Simulations of magnetized multiphase galactic disk regulated by supernovae explosions

The Origin and Universality of the Stellar Initial Mass Function

The role of cosmic rays on magnetic field diffusion and the formation of protostellar discs

The W43-MM1 mini-starburst ridge, a test for star formation efficiency models

Analytical theory for the initial mass function: III time dependence and star formation rate

Cosmic-ray ionisation in collapsing clouds

Ion-neutral friction and accretion-driven turbulence in self-gravitating filaments

On the origin of non self-gravitating filaments in the ISM

The influence of turbulence during magnetized core collapse and its consequences on low-mass star formation

Early stages of cluster formation: fragmentation of massive dense cores down to ~1000 AU

Formation of low-mass stars and brown dwarfs

Formation of proto-clusters and star formation within clusters: apparent universality of the initial mass function ?

Protostellar disk formation and transport of angular momentum during magnetized core collapse

Theories of the massive star formation: a (short) review

Turbulent molecular clouds

UV-driven chemistry in simulations of the interstellar medium. I. Post-processed chemistry with the Meudon PDR code

Analytical star formation rate from gravoturbulent fragmentation

Collapse of massive magnetized dense cores using radiation-magneto-hydrodynamics: early fragmentation inhibition

Dimensional argument for the impact of turbulent support on the stellar initial mass function

Molecular Cloud Evolution IV: Magnetic Fields, Ambipolar Diffusion, and the Star Formation Efficiency

Radiation hydrodynamics with Adaptive Mesh Refinement and application to prestellar core collapse. I Methods

Theories of the initial mass function

Accretion-Driven Turbulence as Universal Process: Galaxies, Molecular Clouds, and Protostellar Disks

Interstellar OH+, H2O+ and H3O+ along the sight-line to G10.6-0.4

Outflows and mass accretion in collapsing dense cores with misaligned rotation axis and magnetic field

Radiative, magnetic and numerical feedbacks on small-scale fragmentation

Star formation: statistical measure of the correlation between the prestellar core mass function and the stellar initial mass function

The Angular Momentum of Magnetized Molecular Cloud Cores: A 2D-3D Comparison

The thermally-unstable warm neutral medium: key for modeling the interstellar medium

Analytical theory for the initial mass function: II. Properties of the flow

Clump morphology and evolution in MHD simulations of molecular cloud formation