BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

P. Hennebelle

P. 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.GAastro-ph.SRastro-ph.IMastro-ph.EPastro-ph.HEastro-phastro-ph.CO

Catalog footprint

What is connected

47works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

47 published item(s)

preprint2025arXiv

ALMAGAL V. Relations between the core populations and the parent clump physical properties

Context. The fragmentation of massive molecular clumps into smaller, potentially star-forming cores plays a key role in the processes of high-mass star formation. The ALMAGAL project offers high-resolution data to investigate these processes across various evolutionary stages in the Galactic plane. Aims. This study aims at correlating the fragmentation properties of massive clumps, obtained from ALMA observations, with their global physical parameters (e.g., mass, surface density, and temperature) and evolutionary indicators (such as luminosity-to-mass ratio and bolometric temperature) obtained from Herschel observations. It seeks to assess whether the cores evolve in number and mass in tandem with their host clumps, and to determine the possible factors influencing the formation of massive cores (M > 24M_\odot). Methods. We analyzed the masses of 6348 fragments, estimated from 1.4 mm continuum data for 1007 ALMAGAL clumps. Leveraging this unprecedentedly large data set, we evaluated statistical relationships between clump parameters, estimated over about 0.1 pc scales, and fragment properties, corresponding to scales of a few 1000 au, while accounting for potential biases related to distance and observational resolution. Our results were further compared with predictions from numerical simulations. Results. The fragmentation level correlates preferentially with clump surface density, supporting a scenario of density-driven fragmentation, whereas it does not show any clear dependence on total clump mass. Both the mass of the most massive core and the core formation efficiency show a broad range and increase on average by an order of magnitude in the intervals spanned by evolutionary indicators such as clump dust temperature and the luminosity-to-mass ratio. This suggests that core growth continues throughout the clump evolution, favoring clump-fed over core-fed theoretical scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Anisotropy in the carbon monoxide (CO) line emission across the Milky Way's disk

We present a study of the $^{12}$CO(1-0) line emission anisotropy across the Milky Way's disk to examine the effect of stellar feedback and Galactic dynamics on the distribution of the dense interstellar medium. The Hessian matrix method is used to characterize the CO line emission distribution and identify the preferential orientation across line-of-sight velocity channels in the Dame et al. 2001 composite Galactic plane survey, which covers the Galactic latitude range $|b|<5^{\circ}$. The structures sampled with this tracer are predominantly parallel to the Galactic plane toward the inner Galaxy, in clear contrast with the predominantly perpendicular orientation of the structures traced by neutral atomic hydrogen (HI) emission toward the same regions. The analysis of the Galactic plane portions sampled at higher angular resolution with other surveys reveals that the alignment with the Galactic plane is also prevalent at smaller scales. We find no preferential orientation in the CO emission toward the outer Galaxy, in contrast with the preferential alignment with the Galactic plane displayed by HI in that portion of the Milky Way. We interpret these results as the combined effect of the decrease in mid-plane pressure with increasing Galactocentric radius and SN feedback lifting diffuse gas more efficiently than dense gas off the Galactic plane.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Evading the dust fragmentation barrier with the streaming instability in protoplanetary disks

Context: The streaming instability (SI) is a leading candidate for reaching solid densities sufficient to trigger the gravitational collapse needed for the formation of planetesimals. However, dust growth barriers appear to impede the ability to assemble sufficiently large dust particles to trigger strong clumping, providing a serious impediment to planetesimal formation. Aims: We aim to address the possibility to enhance dust clumping with dust growth in SI-produced structures, and to estimate the impact of the shift of the dust fragmentation threshold in regions where the SI has enhanced the dust density. Methods: We perform two-dimensional numerical simulations of the SI with a monodisperse description of dust growth, accounting for the impact of mass loading of the dust on the sound speed of the gas and dust mixture when computing dust collisional velocities. Results: Dust mass loading reduces collision velocities in high density regions, allowing dust particles to survive to larger sizes before shattering. In turn, dust clumping is boosted as particles grow in size, as long as they remain sufficiently coupled to the gas. Conclusions: This two-way synergy between dust growth and clumping, which depends on the initial dust-to-gas ratio and dust elastic properties, allows denser dust clumps to form and thus facilitates the onset of planetesimal formation.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Magnetic clumping of charged dust in the dense interstellar medium

Context: Dust grains undergo significant growth in star-forming environments, especially in dense regions prone to gravitational collapse. Although dust is generally assumed to represent $1 \%$ of the gas mass, dust density variations are expected on small scales due to differential dynamics with the gas, leading to enhanced coagulation rates in regions of dust enrichment. Aims: We aim to investigate the clumping of charged dust in the turbulent magnetized dense regions of the interstellar medium. Methods: We develop a dusty model that goes beyond the standard non-ideal MHD and use the code {\ttfamily shark} to perform multifluid 1D simulations of a single size charged dust species and neutral gas with large scale driven turbulence and including ion-neutral friction. Results: We identify a mechanism similar to the parametric instability that efficiently forms dust clumps even in presence of dissipative processes. Such strong clumping survives and is sustained when driving turbulence, and thus high levels of dust concentration are produced due to compressive magnetic effects in regions of shocks. Dust density enhancements are favored by a high transverse-to-longitudinal magnetic ratio which is controlled by: transverse Mach number and plasma parameter. We find that a substantial fraction of dust experiences a density increase of more than a factor of 10 under reasonable conditions, thus promoting dust growth. Conclusion: Our novel dusty non-ideal MHD model shows that dust grains (main charge carriers) are subject to small-scale compressive magnetic effects driven by a parametric instability - like mechanism in regions of shocks, and consequently experience high density enhancements in turbulent environments that go beyond those permitted by pure hydrodynamical processes, making in-situ formation of large grains (sub-mm) in protostellar envelopes a plausible scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

ALMA-IMF I -- Investigating the origin of stellar masses: Introduction to the Large Program and first results

The ALMA-IMF Large Program imaged a total noncontiguous area of 53pc2, covering 15 extreme, nearby protoclusters of the Milky Way. They were selected to span relevant early protocluster evolutionary stages. Our 1.3mm and 3mm observations provide continuum images that are homogeneously sensitive to point-like cores with masses of 0.2 and 0.6Msun, respectively, with a matched spatial resolution of 2000au. We also detect lines that probe the protocluster structure, kinematics, chemistry, and feedback over scales from clouds to filaments to cores. We classify ALMA-IMF protoclusters as Young, Intermediate, or Evolved based on the amount of dense gas in the cloud that has potentially been impacted by HII regions. The ALMA-IMF catalog contains 700 cores that span a mass range of 0.15-250Msun at a typical size of 2100au. We show that this core sample has no significant distance bias and can be used to build core mass functions at similar physical scales. Significant gas motions, which we highlight here in the G353.41 region, are traced down to core scales and can be used to look for inflowing gas streamers and to quantify the impact of the possible associated core mass growth on the shape of the CMF with time. Our first analysis does not reveal any significant evolution of the matter concentration from clouds to cores or from the youngest to more evolved protoclusters, indicating that cloud dynamical evolution and stellar feedback have for the moment only had a slight effect on the structure of high-density gas in our sample. Furthermore, the first-look analysis of the line richness toward bright cores indicates that the survey encompasses several tens of hot cores, of which we highlight the most massive in the G351.77 cloud. Their homogeneous characterization can be used to constrain the emerging molecular complexity in protostars of high to intermediate masses.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Angular momentum profiles of Class 0 protostellar envelopes

[abridged] Understanding how the infalling gas redistribute most of its initial angular momentum inherited from prestellar cores before reaching the stellar embryo is a key question. Disk formation has been naturally considered as a possible solution to this "angular momentum problem". However, how the initial angular momentum of protostellar cores is distributed and evolves during the main accretion phase and the beginning of disk formation has largely remained unconstrained up to now. In the framework of the IRAM CALYPSO survey, we used high dynamic range C$^{18}$O (2-1) and N$_2$H$^+$ (1-0) observations to quantify the distribution of specific angular momentum along the equatorial axis in a sample of 12 Class 0 protostellar envelopes from scales ~50 to 10000 au. The radial distributions of specific angular momentum in the CALYPSO sample suggest two distinct regimes within protostellar envelopes: the specific angular momentum decreases as $j \propto r^{1.6 \pm 0.2}$ down to ~1600 au and then tends to become relatively constant around 6 $\times$ 10$^{-4}$ km s$^{-1}$ pc down to ~50 au. The values of specific angular momentum measured in the inner Class 0 envelopes, namely that of the material directly involved in the star formation process ($<$1600 au), is on the same order of magnitude as what is inferred in small T-Tauri disks. Thus, disk formation appears to be a direct consequence of angular momentum conservation during the collapse. Our analysis reveals a dispersion of the directions of velocity gradients at envelope scales $>$1600 au, suggesting that they may not be related to rotational motions of the envelopes. We conclude that the specific angular momentum observed at these scales could find its origin in core-forming motions (infall, turbulence) or trace an imprint of the initial conditions for the formation of protostellar cores.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Searching for kinematic evidence of Keplerian disks around Class 0 protostars with CALYPSO

The formation of protoplanetary disks is not well understood. To understand how and when these disks are formed, it is crucial to characterize the kinematics of the youngest protostars at a high angular resolution. Here we study a sample of 16 Class 0 protostars to measure their rotation profile at scales from 50 to 500~au and search for Keplerian rotation. We used high-angular-resolution line observations obtained with the Plateau de Bure Interferometer as part of the CALYPSO large program. From $\mathrm{^{13}CO}$ ($J=2-1$), $\mathrm{C^{18}O}$ ($J=2-1$) and SO ($N_{j}=5_{6}-4_{5}$) moment maps, we find that seven sources show rotation about the jet axis at a few hundred au scales: SerpS-MM18, L1448-C, L1448-NB, L1527, NGC1333- IRAS2A, NGC1333-IRAS4B, and SVS13-B. We analyzed the kinematics of these sources in the $uv$ plane to derive the rotation profiles down to 50~au scales. We find evidence for Keplerian rotation in only two sources, L1527 and L1448-C. Overall, this suggests that Keplerian disks larger than 50 au are uncommon around Class 0 protostars. However, in some of the sources, the line emission could be optically thick and dominated by the envelope emission. Due to the optical thickness of these envelopes, some of the disks could have remained undetected in our observations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Probing the cold magnetized Universe with SPICA-POL (B-BOP)

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The dynamically young outflow of the Class 0 protostar Cha-MMS1

On the basis of its low luminosity, its chemical composition, and the absence of a large-scale outflow, the dense core Cha-MMS1 located in the Chamaeleon I molecular cloud was proposed as a first hydrostatic core (FHSC) candidate a decade ago. Our goal is to test this hypothesis by searching for a slow, compact outflow driven by Cha-MMS1 that would match the predictions of MHD simulations for this short phase of star formation. We use the Atacama Large Millimetre/submillimetre Array (ALMA) to map Cha-MMS1 at high angular resolution in CO 3-2 and 13CO 3-2 as well as in continuum emission. We report the detection of a bipolar outflow emanating from the central core, along a (projected) direction roughly parallel to the filament in which Cha-MMS1 is embedded and perpendicular to the large-scale magnetic field. The morphology of the outflow indicates that its axis lies close to the plane of the sky. We measure velocities corrected for inclination of more than 90km/s which is clearly incompatible with the expected properties of a FHSC outflow. Several properties of the outflow are determined and compared to previous studies of Class 0 and Class I protostars. The outflow of Cha-MMS1 has a much smaller momentum force than the outflows of other Class 0 protostars. In addition, we find a dynamical age of 200-3000yr indicating that Cha-MMS1 might be one of the youngest ever observed Class 0 protostars. While the existence of the outflow suggests the presence of a disk, no disk is detected in continuum emission and we derive an upper limit of 55au to its radius. We conclude that Cha-MMS1 has already gone through the FHSC phase and is a young Class 0 protostar, but it has not brought its outflow to full power yet.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

First image of the L1157 molecular jet by the CALYPSO IRAM-PdBI survey

Fast jets are thought to be a crucial ingredient of star formation because they might extract angular momentum from the disk and thus allow mass accretion onto the star. However, it is unclear whether jets are ubiquitous, and likewise, their contribution to mass and angular momentum extraction during protostar formation remains an open question. Our aim is to investigate the ejection process in the low-mass Class 0 protostar L1157. This source is associated with a spectacular bipolar outflow, and the recent detection of high-velocity SiO suggests the occurrence of a jet. Observations of CO 2-1 and SiO 5-4 at 0.8" resolution were obtained with the IRAM Plateau de Bure Interferometer as part of the CALYPSO large program. The jet and outflow structure were fit with a precession model. We derived the column density of CO and SiO, as well as the jet mass-loss rate and mechanical luminosity. High-velocity CO and SiO emission resolve for the first time the first 200 au of the outflow-driving molecular jet. The jet is strongly asymmetric, with the blue lobe 0.65 times slower than the red lobe. This suggests that the large-scale asymmetry of the outflow is directly linked to the jet velocity and that the asymmetry in the launching mechanism has been at work for the past 1800 yr. Velocity asymmetries are common in T Tauri stars, which suggests that the jet formation mechanism from Class 0 to Class II stages might be similar. Our model simultaneously fits the inner jet and the clumpy 0.2 pc scale outflow by assuming that the jet precesses counter-clockwise on a cone inclined by 73 degree to the line of sight with an opening angle of 8 degree on a period of 1640 yr. The estimated jet mass flux and mechanical luminosity are 7.7e-7 Msun/yr, and 0.9 Lsun, indicating that the jet could extract at least 25% of the gravitational energy released by the forming star.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Magnetic field morphology in nearby molecular clouds as revealed by starlight and submillimetre polarization

Within four nearby (d < 160 pc) molecular clouds, we statistically evaluate the structure of the interstellar magnetic field, projected on the plane of the sky and integrated along the line of sight, as inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz and from the optical and NIR polarization of background starlight. We compare the dispersion of the field orientation directly in vicinities with an area equivalent to that subtended by the Planck effective beam at 353 GHz (10') and using the second-order structure functions of the field orientation angles. We find that the average dispersion of the starlight-inferred field orientations within 10'-diameter vicinities is less than 20 deg, and that at these scales the mean field orientation is on average within 5 deg of that inferred from the submillimetre polarization observations in the considered regions. We also find that the dispersion of starlight polarization orientations and the polarization fractions within these vicinities are well reproduced by a Gaussian model of the turbulent structure of the magnetic field, in agreement with the findings reported by the Planck collaboration at scales greater than 10' and for comparable column densities. At scales greater than 10', we find differences of up to 14.7 deg between the second-order structure functions obtained from starlight and submillimetre polarization observations in the same positions in the plane of the sky, but comparison with a Gaussian model of the turbulent structure of the magnetic field indicates that these differences are small and are consistent with the difference in angular resolution between both techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Protostars: Forges of cosmic rays?

Galactic cosmic rays (CR) are particles presumably accelerated in supernova remnant shocks that propagate in the interstellar medium up to the densest parts of molecular clouds, losing energy and their ionisation efficiency because of the presence of magnetic fields and collisions with molecular hydrogen. Recent observations hint at high levels of ionisation and at the presence of synchrotron emission in protostellar systems, which leads to an apparent contradiction. We want to explain the origin of these CRs accelerated within young protostars as suggested by observations. Our modelling consists of a set of conditions that has to be satisfied in order to have an efficient CR acceleration through diffusive shock acceleration. We analyse three main acceleration sites, then we follow the propagation of these particles through the protostellar system up to the hot spot region. We find that jet shocks can be strong accelerators of CR protons, which can be boosted up to relativistic energies. Other promising acceleration sites are protostellar surfaces, where shocks caused by impacting material during the collapse phase are strong enough to accelerate CR protons. In contrast, accretion flow shocks are too weak to efficiently accelerate CRs. Though CR electrons are weakly accelerated, they can gain a strong boost to relativistic energies through re-acceleration in successive shocks. We suggest a mechanism able to accelerate both CR protons and electrons through the diffusive shock acceleration mechanism, which can be used to explain the high ionisation rate and the synchrotron emission observed towards protostellar sources. The existence of an internal source of energetic particles can have a strong and unforeseen impact on the ionisation of the protostellar disc, on the star and planet formation processes, and on the formation of pre-biotic molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

From forced collapse to H ii region expansion in Mon R2: Envelope density structure and age determination with Herschel

The surroundings of HII regions can have a profound influence on their development, morphology, and evolution. This paper explores the effect of the environment on H II regions in the MonR2 molecular cloud. We aim to investigate the density structure of envelopes surrounding HII regions and to determine their collapse and ionisation expansion ages. The Mon R2 molecular cloud is an ideal target since it hosts an H II region association. Column density and temperature images derived from Herschel data were used together to model the structure of HII bubbles and their surrounding envelopes. The resulting observational constraints were used to follow the development of the Mon R2 ionised regions with analytical calculations and numerical simulations. The four hot bubbles associated with H II regions are surrounded by dense, cold, and neutral gas envelopes. The radial density profiles are reminiscent of those of low-mass protostellar envelopes. The inner parts of envelopes of all four HII regions could be free-falling because they display shallow density profiles. As for their outer parts, the two compact HII regions show a density profile, which is typical of the equilibrium structure of an isothermal sphere. In contrast, the central UCHii region shows a steeper outer profile, that could be interpreted as material being forced to collapse. The size of the heated bubbles, the spectral type of the irradiating stars, and the mean initial neutral gas density are used to estimate the ionisation expansion time, texp, 0.1Myr,for the dense UCHII and compact HII regions and 0.35 Myr for the extended HII region. The envelope transition radii between the shallow and steeper density profiles are used to estimate the time elapsed since the formation of the first proto stellar embryo, Tinf : 1Myr, for the ultra-compact, 1.5 / 3Myr for the compact, and greater than 6Myr for the extended HII regions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Incorporating Ambipolar and Ohmic Diffusion in the AMR MHD code RAMSES

We have implemented non-ideal Magneto-Hydrodynamics (MHD) effects in the Adaptive Mesh Refinement (AMR) code RAMSES, namely ambipolar diffusion and Ohmic dissipation, as additional source terms in the ideal MHD equations. We describe in details how we have discretized these terms using the adaptive Cartesian mesh, and how the time step is diminished with respect to the ideal case, in order to perform a stable time integration. We have performed a large suite of test runs, featuring the Barenblatt diffusion test, the Ohmic diffusion test, the C-shock test and the Alfven wave test. For the latter, we have performed a careful truncation error analysis to estimate the magnitude of the numerical diffusion induced by our Godunov scheme, allowing us to estimate the spatial resolution that is required to address non-ideal MHD effects reliably. We show that our scheme is second-order accurate, and is therefore ideally suited to study non-ideal MHD effects in the context of star formation and molecular cloud dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Jet multiplicity in the proto-binary system NGC1333-IRAS4A. The detailed CALYPSO IRAM-PdBI view

Owing to the paucity of sub-arcsecond (sub)mm observations required to probe the innermost regions of newly forming protostars, several fundamental questions are still being debated, such as the existence and coevality of close multiple systems. We study the physical and chemical properties of the jets and protostellar sources in the NGC1333-IRAS4A proto-binary system using continuum emission and molecular tracers of shocked gas. We observed NGC1333-IRAS4A in the SiO(6-5), SO(6_5-5_4), and CO(2-1) lines and the continuum emission at 1.3, 1.4, and 3 mm using the IRAM Plateau de Bure Interferometer in the framework of the CALYPSO large program. We clearly disentangle for the first time the outflow emission from the two sources A1 and A2. The two protostellar jets have very different properties: the A1 jet is faster, has a short dynamical timescale (<10^3 yr), and is associated with H2 shocked emission, whereas the A2 jet, which dominates the large-scale emission, is associated with diffuse emission, bends, and emits at slower velocities. The observed bending of the A2 jet is consistent with the change of propagation direction observed at large scale and suggests jet precession on very short timescales (~200-600 yr). In addition, a chemically rich spectrum with emission from several COMs (e.g. HCOOH, CH3OCHO, CH3OCH3) is only detected towards A2. Finally, very high-velocity shocked emission (~50 km s^-1) is observed along the A1 jet. An LTE analysis shows that SiO, SO, and H2CO abundances in the gas phase are enhanced up to (3-4)x10^{-7}, (1.4-1.7)x10^{-6}, and (3-7.9)x10^{-7}, respectively. The intrinsic different properties of the jets and driving sources in NGC1333-IRAS4A suggest different evolutionary stages for the two protostars, with A1 being younger than A2, in a very early stage of star formation previous to the hot-corino phase.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Planck intermediate results. XXXV. Probing the role of the magnetic field in the formation of structure in molecular clouds

Within ten nearby (d < 450 pc) Gould Belt molecular clouds we evaluate statistically the relative orientation between the magnetic field projected on the plane of sky, inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz, and the gas column density structures, quantified by the gradient of the column density, $N_H$. The selected regions, covering several degrees in size, are analyzed at an effective angular resolution of 10' FWHM, thus sampling physical scales from 0.4 to 40 pc in the nearest cloud. The column densities in the selected regions range from $N_H \approx 10^{21}$ to $10^{23}$ cm$^{-2}$, and hence they correspond to the bulk of the molecular clouds. The relative orientation is evaluated pixel by pixel and analyzed in bins of column density using the novel statistical tool called "Histogram of Relative Orientations". Throughout this study, we assume that the polarized emission observed by Planck at 353 GHz is representative of the projected morphology of the magnetic field in each region, i.e., we assume a constant dust grain alignment efficiency, independent of the local environment. Within most clouds we find that the relative orientation changes progressively with increasing $N_H$, from preferentially parallel or having no preferred orientation to preferentially perpendicular. In simulations of magnetohydrodynamic turbulence in molecular clouds this trend in relative orientation is a signature of Alfvénic or sub-Alfvénic turbulence, implying that the magnetic field is significant for the gas dynamics at the scales probed by Planck. We compare the deduced magnetic field strength with estimates we obtain from other methods and discuss the implications of the Planck observations for the general picture of molecular cloud formation and evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fragmentation of massive dense cores down to ~1000 AU: Relation between fragmentation and density structure

In order to shed light on the main physical processes controlling fragmentation of massive dense cores, we present a uniform study of the density structure of 19 massive dense cores, selected to be at similar evolutionary stages, for which their relative fragmentation level was assessed in a previous work. We inferred the density structure of the 19 cores through a simultaneous fit of the radial intensity profiles at 450 and 850 micron (or 1.2 mm in two cases) and the Spectral Energy Distribution, assuming spherical symmetry and that the density and temperature of the cores decrease with radius following power-laws. We find a weak (inverse) trend of fragmentation level and density power-law index, with steeper density profiles tending to show lower fragmentation, and vice versa. In addition, we find a trend of fragmentation increasing with density within a given radius, which arises from a combination of flat density profile and high central density and is consistent with Jeans fragmentation. We considered the effects of rotational-to-gravitational energy ratio, non-thermal velocity dispersion, and turbulence mode on the density structure of the cores, and found that compressive turbulence seems to yield higher central densities. Finally, a possible explanation for the origin of cores with concentrated density profiles, which are the cores showing no fragmentation, could be related with a strong magnetic field, consistent with the outcome of radiation magnetohydrodynamic simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Variations of the stellar initial mass function in the progenitors of massive early-type galaxies and in extreme starburst environments

We examine variations of the stellar initial mass function (IMF) in extreme environments within the formalism derived by Hennebelle \& Chabrier. We focus on conditions encountered in progenitors of massive early type galaxies and starburst regions. We show that, when applying the concept of turbulent Jeans mass as the characteristic mass for fragmentation in a turbulent medium, instead of the standard thermal Jeans mass for purely gravitational fragmentation, the peak of the IMF in such environments is shifted towards smaller masses, leading to a bottom-heavy IMF, as suggested by various observations. In very dense and turbulent environments, we predict that the high-mass tail of the IMF can become even steeper than the standard Salpeter IMF, with a limit for the power law exponent $α\simeq -2.7$, in agreement with recent observational determinations. This steepening is a direct consequence of the high densities and Mach values in such regions but also of the time dependence of the fragmentation process, as incorporated in the Hennebelle-Chabrier theory. We provide analytical parametrizations of these IMFs in such environments, to be used in galaxy evolution calculations. We also calculate the star formation rates and the mass-to-light ratios expected under such extreme conditions and show that they agree well with the values inferred in starburst environments and massive high-redshift galaxies. This reinforces the paradigm of star formation as being a universal process, i.e. the direct outcome of gravitationally unstable fluctuations in a density field initially generated by large scale shock-dominated turbulence. This globally enables us to infer the variations of the stellar IMF and related properties for atypical galactic conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A sub-parsec resolution simulation of the Milky Way: Global structure of the ISM and properties of molecular clouds

We present a self-consistent hydrodynamical simulation of a Milky Way-like galaxy, at the resolution of 0.05 pc. The model includes star formation and a new implementation of stellar feedback through photo-ionization, radiative pressure and supernovae. The simulation resolves the structure of the interstellar medium at subparsec resolution for a few cloud lifetimes, and at 0.05 pc for about a cloud crossing time. Turbulence cascade and gravitation from the kpc scales are de facto included in smaller structures like molecular clouds. We show that the formation of a bar influences the dynamics of the central ~100 pc by creating resonances. At larger radii, the spiral arms host the formation of regularly spaced clouds: beads on a string and spurs. These instabilities pump turbulent energy into the gas, generally in the supersonic regime. Because of asymmetric drift, the supernovae explode outside of their gaseous nursery, which diminishes the effect of feedback on the structure of clouds. The evolution of clouds is thus mostly due to fragmentation and gas consumption, regulated mainly by supersonic turbulence. The transition from turbulence supported to self-gravitating gas is detected in the gas density probability distribution function at ~2000 cm^-3. The power spectrum density suggests that gravitation governs the hierarchical organisation of structures from the galactic scale down to a few parsecs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An Imprint of Molecular Cloud Magnetization in the Morphology of the Dust Polarized Emission

We describe a morphological imprint of magnetization found when considering the relative orientation of the magnetic field direction with respect to the density structures in simulated turbulent molecular clouds. This imprint was found using the Histogram of Relative Orientations (HRO): a new technique that utilizes the gradient to characterize the directionality of density and column density structures on multiple scales. We present results of the HRO analysis in three models of molecular clouds in which the initial magnetic field strength is varied, but an identical initial turbulent velocity field is introduced, which subsequently decays. The HRO analysis was applied to the simulated data cubes and mock-observations of the simulations produced by integrating the data cube along particular lines of sight. In the 3D analysis we describe the relative orientation of the magnetic field $\mathbf{B}$ with respect to the density structures, showing that: 1.The magnetic field shows a preferential orientation parallel to most of the density structures in the three simulated cubes. 2.The relative orientation changes from parallel to perpendicular in regions with density over a critical density $n_{T}$ in the highest magnetization case. 3.The change of relative orientation is largest for the highest magnetization and decreases in lower magnetization cases. This change in the relative orientation is also present in the projected maps. In conjunction with simulations HROs can be used to establish a link between the observed morphology in polarization maps and the physics included in simulations of molecular clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Global collapse of molecular clouds as a formation mechanism for the most massive stars

The relative importance of primordial molecular cloud fragmentation versus large-scale accretion still remains to be assessed in the context of massive core/star formation. Studying the kinematics of the dense gas surrounding massive-star progenitors can tell us the extent to which large-scale flow of material impacts the growth in mass of star-forming cores. Here we present a comprehensive dataset of the 5500(+/-800) Msun infrared dark cloud SDC335.579-0.272 (hereafter SDC335) which exhibits a network of cold, dense, parsec-long filaments. Atacama Large Millimeter Array (ALMA) Cycle 0 observations reveal two massive star-forming cores, MM1 and MM2, sitting at the centre of SDC335 where the filaments intersect. With a gas mass of 545(+770,-385) Msun contained within a source diameter of 0.05pc, MM1 is one of the most massive, compact protostellar cores ever observed in the Galaxy. As a whole, SDC335 could potentially form an OB cluster similar to the Trapezium cluster in Orion. ALMA and Mopra single-dish observations of the SDC335 dense gas furthermore reveal that the kinematics of this hub-filament system are consistent with a global collapse of the cloud. These molecular-line data point towards an infall velocity V_{inf} =0.7(+/-0.2) km/s, and a total mass infall rate \dot{M}_{inf} = 2.5(+/-1.0) x 10^{-3} Msun/yr towards the central pc-size region of SDC335. This infall rate brings 750(+/-300) Msun of gas to the centre of the cloud per free-fall time (t_{ff}=3x10^5 yr). This is enough to double the mass already present in the central pc-size region in 3.5(+2.2,-1.0) x t_{ff}. These values suggest that the global collapse of SDC335 over the past million year resulted in the formation of an early O-type star progenitor at the centre of the cloud's gravitational potential well.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

SPLASH: The Southern Parkes Large-Area Survey in Hydroxyl - First Science from the Pilot Region

SPLASH (the Southern Parkes Large-Area Survey in Hydroxyl) is a sensitive, unbiased and fully-sampled survey of the Southern Galactic Plane and Galactic Centre in all four ground-state transitions of the hydroxyl (OH) radical. The survey provides a deep census of 1612-, 1665-, 1667- and 1720-MHz OH absorption and emission from the Galactic ISM, and is also an unbiased search for maser sources in these transitions. We present here first results from the SPLASH pilot region, which covers Galactic longitudes 334 to 344 degrees and latitudes of -2 to +2 degrees. Diffuse OH is widely detected in all four transitions, with optical depths that are always small (averaged over the Parkes beam), and with departures from LTE common even in the 1665- and 1667-MHz main lines. To a 3$σ$ sensitivity of 30 mK, we find no evidence of OH envelopes extending beyond the CO-bright regions of molecular cloud complexes, and conclude that the similarity of the OH excitation temperature and the level of the continuum background is at least partly responsible for this. We detect masers and maser candidates in all four transitions, approximately 50 per cent of which are new detections. This implies that SPLASH will produce a substantial increase in the known population of ground-state OH masers in the Southern Galactic Plane.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The structure of the thermally bistable and turbulent atomic gas in the local interstellar medium

This paper is a numerical study of the condensation of the warm neutral medium (WNM) into cold neutral medium (CNM) structures under the effect of turbulence and thermal instability. Using low resolution simulations we explored the impact of the WNM initial density and properties of the turbulence (stirring in Fourier with a varying mix of solenoidal and compressive modes) on the cold gas formation. Two sets of initial conditions which match the observations were selected to produce high resolution simulations (1024^3) allowing to study in details the properties of the produced dense structures. For typical values of the density, pressure and velocity dispersion of the WNM in the solar neighborhood, the turbulent motions of the HI can not provoque the phase transition from WNM to CNM, whatever their amplitude and their distribution in solenoidal and compressive modes. On the other hand we show that a quasi-isothermal increase in WNM density of a factor of 2 to 4 is enough to induce the phase transition, leading to the transition of about 40 percent of the gas to the cold phase within 1 Myr. Given the observed properties of the HI in the local ISM, the WNM and individual CNM structures in the local ISM are sub or transsonic and their dynamics are tightly interwoven. The velocity field bears the evidence of subsonic turbulence with a 2D power spectrum following the Kolmogorov law as P(k) \sim k^{-8/3} while the density is highly contrasted with a singificantly shallower power spectrum, reminiscent of what is observed in the cold ISM. Supra-thermal line width observed for CNM might be the result of relative velocity between cold structures. Finally, the cold structures denser than 5 cm^{-3} reproduce well the laws M \sim L^{2.25-2.28} and sigma(v) \sim 0.5-0.8L^{1/3} generally observed in molecular clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

World-leading science with SPIRou - the nIR spectropolarimeter / high-precision velocimeter for CFHT

SPIRou is a near-infrared (nIR) spectropolarimeter / velocimeter proposed as a new-generation instrument for CFHT. SPIRou aims in particular at becoming world-leader on two forefront science topics, (i) the quest for habitable Earth-like planets around very- low-mass stars, and (ii) the study of low-mass star and planet formation in the presence of magnetic fields. In addition to these two main goals, SPIRou will be able to tackle many key programs, from weather patterns on brown dwarf to solar-system planet atmospheres, to dynamo processes in fully-convective bodies and planet habitability. The science programs that SPIRou proposes to tackle are forefront (identified as first priorities by most research agencies worldwide), ambitious (competitive and complementary with science programs carried out on much larger facilities, such as ALMA and JWST) and timely (ideally phased with complementary space missions like TESS and CHEOPS). SPIRou is designed to carry out its science mission with maximum efficiency and optimum precision. More specifically, SPIRou will be able to cover a very wide single-shot nIR spectral domain (0.98-2.35 μm) at a resolving power of 73.5K, providing unpolarized and polarized spectra of low-mass stars with a ~15% average throughput and a radial velocity (RV) precision of 1 m/s.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Adaptable Radiative Transfer Innovations for Submillimetre Telescopes (ARTIST) - Dust polarisation module (DustPol)

We present a new publicly available tool (DustPol) aimed to model the polarised thermal dust emission. The module DustPol, which is publicly available, is part of the ARTIST (Adaptable Radiative Transfer Innovations for Submillimetre Telescopes) package, which also offers tools for modelling the polarisation of line emission together with a model library and a Python-based user interface. DustPol can easily manage analytical as well as pre-gridded models to generate synthetic maps of the Stokes I, Q, and U parameters. These maps are stored in FITS format which is straightforwardly read by the data reduction software used, e.g., by the Atacama Large Millimeter Array (ALMA). This turns DustPol into a powerful engine for the prediction of the expected polarisation features of a source observed with ALMA or the Planck satellite as well as for the interpretation of existing submillimetre observations obtained with other telescopes. DustPol allows the parameterisation of the maximum degree of polarisation and we find that, in a prestellar core, if there is depolarisation, this effect should happen at densities of 10^6 cm-3 or larger. We compare a model generated by DustPol with the observational polarisation data of the low-mass Class 0 object NGC 1333 IRAS 4A, finding that the total and the polarised emission are consistent.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Collapse, outflows and fragmentation of massive, turbulent and magnetized prestellar barotropic cores

Stars and more particularly massive stars, have a drastic impact on galaxy evolution. Yet the conditions in which they form and collapse are still not fully understood. In particular, the influence of the magnetic field on the collapse of massive clumps is relatively unexplored, it is thus of great relevance in the context of the formation of massive stars to investigate its impact. We perform high resolution, MHD simulations of the collapse of hundred solar masses, turbulent and magnetized clouds, using the adaptive mesh refinement code RAMSES. We compute various quantities such as mass distribution, magnetic field and angular momentum within the collapsing core and study the episodic outflows and the fragmentation that occurs during the collapse. The magnetic field has a drastic impact on the cloud evolution. We find that magnetic braking is able to substantially reduce the angular momentum in the inner part of the collapsing cloud. Fast and episodic outflows are being launched with typical velocities of the order of 3-5 km s$^{-1}$ although the highest velocities can be as high as 30-40 km s$^{-1}$. The fragmentation in several objects, is reduced in substantially magnetized clouds with respect to hydrodynamical ones by a factor of the order of 1.5-2. We conclude that magnetic fields have a significant impact on the evolution of massive clumps. In combination with radiation, magnetic fields largely determine the outcome of massive core collapse. We stress that numerical convergence of MHD collapse is a challenging issue. In particular, numerical diffusion appears to be important at high density therefore possibly leading to an over-estimation of the number of fragments.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

CHESS, Chemical Herschel surveys of star forming regions:Peering into the protostellar shock L1157-B1

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Detection of hydrogen fluoride absorption in diffuse molecular clouds with Herschel/HIFI: a ubiquitous tracer of molecular gas

We discuss the detection of absorption by interstellar hydrogen fluoride (HF) along the sight line to the submillimeter continuum sources W49N and W51. We have used Herschel's HIFI instrument in dual beam switch mode to observe the 1232.4762 GHz J = 1 - 0 HF transition in the upper sideband of the band 5a receiver. We detected foreground absorption by HF toward both sources over a wide range of velocities. Optically thin absorption components were detected on both sight lines, allowing us to measure - as opposed to obtain a lower limit on - the column density of HF for the first time. As in previous observations of HF toward the source G10.6-0.4, the derived HF column density is typically comparable to that of water vapor, even though the elemental abundance of oxygen is greater than that of fluorine by four orders of magnitude. We used the rather uncertain N(CH)-N(H2) relationship derived previously toward diffuse molecular clouds to infer the molecular hydrogen column density in the clouds exhibiting HF absorption. Within the uncertainties, we find that the abundance of HF with respect to H2 is consistent with the theoretical prediction that HF is the main reservoir of gas-phase fluorine for these clouds. Thus, hydrogen fluoride has the potential to become an excellent tracer of molecular hydrogen, and provides a sensitive probe of clouds of small H2 column density. Indeed, the observations of hydrogen fluoride reported here reveal the presence of a low column density diffuse molecular cloud along the W51 sight line, at an LSR velocity of ~ 24kms-1, that had not been identified in molecular absorption line studies prior to the launch of Herschel.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamic star formation in the massive DR21 filament

The formation of massive stars is a highly complex process in which it is not clear whether the star-forming gas is in global gravitational collapse or in an equilibrium state, supported by turbulence. By studying one of the most massive and dense star-forming regions in the Galaxy at a distance of less than 3 kpc, the filament containing the well-known sources DR21 and DR21(OH), we expect to find observational signatures that allow to discriminate between the two views. We use molecular line data from our 13CO 1-0, CS 2-1, and N2H+ 1-0 survey of the Cygnus X region obtained with the FCRAO and high-angular resolution observations of CO, CS, HCO+, N2H+, and H2CO, obtained with the IRAM 30m telescope. We observe a complex velocity field and velocity dispersion in the DR21 filament in which regions of highest column-density, i.e. dense cores, have a lower velocity dispersion than the surrounding gas and velocity gradients that are not (only) due to rotation. Infall signatures in optically thick line profiles of HCO+ and 12CO are observed along and across the whole DR21 filament. From modelling the observed spectra, we obtain a typical infall speed of 0.6 km/s and mass accretion rates of the order of a few 10^-3 Msun/yr for the two main clumps constituting the filament. These massive (4900 and 3300 Msun) clumps are both gravitationally contracting. All observed kinematic features in the DR21 filament can be explained if it is formed by the convergence of flows at large scales and is now in a state of global gravitational collapse. Whether this convergence of flows originated from self-gravity at larger scales or from other processes can not be settled with the present study. The observed velocity field and velocity dispersion are consistent with results from (magneto)-hydrodynamic simulations where the cores lie at the stagnation points of convergent turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

First detection of ND in the solar-mass protostar IRAS16293-2422

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Herschel/HIFI discovery of interstellar chloronium (H$_2$Cl$^+$)

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Herschel/HIFI observations of interstellar OH+ and H2O+ towards W49N: a probe of diffuse clouds with a small molecular fraction

We report the detection of absorption by interstellar hydroxyl cations and water cations, along the sight-line to the bright continuum source W49N. We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 972 GHz N = 1 - 0 transition of OH+ and the 1115 GHz 1(11) - 0(00) transition of ortho-H2O+. The resultant spectra show absorption by ortho-H2O+, and strong absorption by OH+, in foreground material at velocities in the range 0 to 70 km/s with respect to the local standard of rest. The inferred OH+/H2O+ abundance ratio ranges from ~ 3 to ~ 15, implying that the observed OH+ arises in clouds of small molecular fraction, in the 2 - 8% range. This conclusion is confirmed by the distribution of OH+ and H2O+ in Doppler velocity space, which is similar to that of atomic hydrogen, as observed by means of 21 cm absorption measurements, and dissimilar from that typical of other molecular tracers. The observed OH+/H abundance ratio of a few E-8 suggests a cosmic ray ionization rate for atomic hydrogen of (0.6 - 2.4) E-16 s-1, in good agreement with estimates inferred previously for diffuse clouds in the Galactic disk from observations of interstellar H3+ and other species.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Herschel/HIFI observations of spectrally resolved methylidyne signatures toward the high-mass star-forming core NGC6334I

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ortho-to-para ratio of interstellar heavy water

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Strong absorption by interstellar hydrogen fluoride: Herschel/HIFI observations of the sight-line to G10.6-0.4 (W31C)

We report the detection of strong absorption by interstellar hydrogen fluoride along the sight-line to the submillimeter continuum source G10.6-0.4 (W31C). We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 1232.4763 GHz J=1-0 HF transition in the upper sideband of the Band 5a receiver. The resultant spectrum shows weak HF emission from G10.6-0.4 at LSR velocities in the range -10 to -3 km/s, accompanied by strong absorption by foreground material at LSR velocities in the range 15 to 50 km/s. The spectrum is similar to that of the 1113.3430 GHz 1(11)-0(00) transition of para-water, although at some frequencies the HF (hydrogen fluoride) optical depth clearly exceeds that of para-H2O. The optically-thick HF absorption that we have observed places a conservative lower limit of 1.6E+14 cm-2 on the HF column density along the sight-line to G10.6-0.4. Our lower limit on the HF abundance, 6E-9 relative to hydrogen nuclei, implies that hydrogen fluoride accounts for between ~ 30 and 100% of the fluorine nuclei in the gas phase along this sight-line. This observation corroborates theoretical predictions that - because the unique thermochemistry of fluorine permits the exothermic reaction of F atoms with molecular hydrogen - HF will be the dominant reservoir of interstellar fluorine under a wide range of conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The distribution of water in the high-mass star-forming region NGC 6334I

We present observations of twelve rotational transitions of H2O-16, H2O-18, and H2O-17 toward the massive star-forming region NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical conditions in these components. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of 200 K, whereas water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 4 10^-5 (outflow) and 10^-8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust grains. The H2O-18/H2O-17 abundance ratio is 3.2, which is consistent with the O-18/O-17 ratio determined from CO isotopologues. The ortho/para ratio in water appears to be relatively low 1.6(1) in the cold, quiescent gas, but close to the equilibrium value of three in the warmer outflow material (2.5(0.8)).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The methanol lines and hot core of OMC2-FIR4, an intermediate-mass protostar, with Herschel-HIFI

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10^-6 is placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display asymmetric profiles, which may arise from infall around the hot core.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Toward understanding the formation of multiple systems - A pilot IRAM-PdBI survey of Class 0 objects

The formation process of binary stars and multiple systems is poorly understood. Here, we seek to determine the typical outcome of protostellar collapse and to constrain models of binary formation by core fragmentation during collapse, using high-resolution millimeter continuum imaging of very young (Class 0) protostars observed at the beginning of the main accretion phase. We carried out a pilot high-resolution study of 5 Class 0 objects, using the most extended (A) configuration of the IRAM PdBI at 1.3 mm, which allow us to probe the multiplicity of Class 0 protostars down to separations a ~50 AU and circumstellar mass ratios q ~0.07. We show that our PdBI observations revealed only wide (>1500 AU) protobinary systems and/or outflow-generated features. When combined with previous millimeter interferometric observations of Class 0 protostars, our pilot PdBI study tentatively suggests that the binary fraction in the ~ 75-1000 AU range increases from the Class 0 to the Class I stage. It also seems to argue against purely hydrodynamic models of binary star formation. We briefly discuss possible alternative scenarios to reconcile the low multiplicity rate of Class 0 protostars on small scales with the higher binary fraction observed at later evolutionary stages.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Disk formation during collapse of magnetized protostellar cores

In the context of star and planet formation, understanding the formation of disks is of fundamental importance. Previous studies found that the magnetic field has a very strong impact on the collapse of a prestellar cloud, particularly in possibly suppressing the formation of a disk even for relatively modest values of the magnetic intensity. Since observations infer that cores have a substantial level of magnetization, this raises the question of how disks form. However, most studies have been restricted to the case in which the initial angle, $α$, between the magnetic field and the rotation axis equals 0$^\circ$. We explore and analyse the influence of non aligned configurations on disk formation. We perform 3D ideal MHD, AMR numerical simulations for various values of $μ$, the ratio of the mass-to-flux to the critical mass-to-flux, and various values of $α$. We find that disks form more easily as $α$ increases from 0 to 90$^\circ$. We propose that as the magnetized pseudo-disks become thicker with increasing $α$, the magnetic braking efficiency is lowered. We also find that even small values of $α$ ($\simeq$ 10-20$^\circ$) show significant differences with the alligned case. Within the framework of ideal MHD and for our choice of initial conditions, centrifugally supported disks cannot form for values of $μ$ smaller than $\simeq$3, if the magnetic field and the rotation axis are perpendicular, and smaller than about $\simeq$5-10 when they are perfectly aligned.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the Structure of the Turbulent Interstellar Clouds

{It is well established that the atomic interstellar hydrogen is filling the galaxies and constitutes the building blocks of molecular clouds.} {To understand the formation and the evolution of molecular clouds, it is necessary to investigate the dynamics of turbulent and thermally bistable as well as barotropic flows.} {We perform high resolution 3-dimensional hydrodynamical simulations of 2-phase, isothermal and polytropic flows.} {We compare the density PDF and Mach number density relation in the various simulations and conclude that 2-phase flows behave rather differently than polytropic flows. We also extract the clumps and study their statistical properties such as mass spectrum, mass-size relation and internal velocity dispersion. In each case, it is found that the behaviour is well represented by a simple powerlaw. While the various exponents inferred are very similar for the 2-phase, isothermal and polytropic flows, we nevertheless find significant differences, such as for example the internal velocity dispersion which is smaller for 2-phase flows.} {The structures statistics are very similar to what has been inferred from observations, in particular the mass spectrum, the mass-size relation and the velocity dispersion-size relation are all powerlaws whose indices well agree with the observed values. Our results suggest that in spite of various statistics being similar for 2-phase and polytropic flows, they nevertheless present significant differences, stressing the necessity to consider the proper thermal structure of the interstellar atomic hydrogen for computing its dynamics as well as the formation of molecular clouds.}

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Protostellar collapse: radiative and magnetic feedbacks on small scale fragmentation

It is established that both radiative transfer and magnetic field have a strong impact on the collapse and the fragmentation of prestellar dense cores, but no consistent calculation exists yet at such scales. We present original AMR calculations including magnetic field (in the ideal MHD limit) and radiative transfer, within the Flux Limited Diffusion approximation, of the collapse of a 1 solar mass dense core. We compare the results with calculations performed with a barotropic EOS. We show that radiative transfer has an important impact on the collapse and the fragmentation, through the cooling or heating of the gas, and is complementary of the magnetic field. A larger field yields a stronger magnetic braking, increasing the accretion rate and thus the effect of the radiative feedback. Even for a strongly magnetized core, where the dynamics of the collapse is dominated by the magnetic field, radiative transfer is crucial to determine the temperature and optical depth distributions, two potentially accessible observational diagnostics. A barotropic EOS cannot account for realistic fragmentation. The diffusivity of the numerical scheme, however, is found to strongly affect the output of the collapse, leading eventually to spurious fragmentation. Both radiative transfer and magnetic field must be included in numerical calculations of star formation to obtain realistic collapse configurations and observable signatures. Nevertheless, the numerical resolution and the robustness of the solver are of prime importance to obtain reliable results. When using an accurate solver, the fragmentation is found to always remain inhibited by the magnetic field, at least in the ideal MHD limit, even when radiative transfer is included.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Simulated CII observations for SPICA/SAFARI

We investigate the case of CII 158 micron observations for SPICA/SAFARI using a three-dimensional magnetohydrodynamical (MHD) simulation of the diffuse interstellar medium (ISM) and the Meudon PDR code. The MHD simulation consists of two converging flows of warm gas (10,000 K) within a cubic box 50 pc in length. The interplay of thermal instability, magnetic field and self-gravity leads to the formation of cold, dense clumps within a warm, turbulent interclump medium. We sample several clumps along a line of sight through the simulated cube and use them as input density profiles in the Meudon PDR code. This allows us to derive intensity predictions for the CII 158 micron line and provide time estimates for the mapping of a given sky area.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

On the structure of the turbulent interstellar atomic hydrogen. I- Physical characteristics

{We study in some details the statistical properties of the turbulent 2-phase interstellar atomic gas.{We present high resolution bidimensional numerical simulations of the interstellar atomic hydrogen which describe it over 3 to 4 orders of magnitude in spatial scales.}{The simulations produce naturally small scale structures having either large or small column density. It is tempting to propose that the former are connected to the tiny small scale structures observed in the ISM. We compute the mass spectrum of CNM structures and find that ${\cal N}(M) dM \propto M ^{-1.7} dM$, which is remarkably similar to the mass spectrum inferred for the CO clumps. We propose a theoretical explanation based on a formalism inspired from the Press & Schecter (1974) approach and used the fact that the turbulence within WNM is subsonic. This theory predicts ${\cal N}(M) \propto M ^{-5/3}$ in 2D and ${\cal N}(M) \propto M ^{-16/9}$ in 3D. We compute the velocity and the density power-spectra and conclude that, although the latter is rather flat, as observed in supersonic isothermal simulations, the former follows the Kolmogorov prediction and is dominated by its solenoidal component. This is due to the bistable nature of the flow which produces large density fluctuations even when the rms Mach number (of WNM) is not large. We also find that, whereas the energy at large scales is mainly in the WNM, at smaller scales, it is dominated by the kinetic energy of the CNM fragments.}

0 citations0 reviewsNo code linkedOpen access