-
Shocks in the warm neutral medium II -- Origin of neutral carbon at high pressure
Authors:
Benjamin Godard,
Guillaume Pineau Des Forêts,
Jeremy La Porte,
Mona Merlin-Weck
Abstract:
Aims: Ultraviolet (UV) lines of neutral carbon observed in absorption in the local diffuse interstellar medium (ISM) have long revealed that a substantial fraction of the mass of the gas lies at a thermal pressure one to three orders of magnitude above that of the bulk of the ISM. In this paper, we propose that this enigmatic component originates from shocks propagating at intermediate (…
▽ More
Aims: Ultraviolet (UV) lines of neutral carbon observed in absorption in the local diffuse interstellar medium (ISM) have long revealed that a substantial fraction of the mass of the gas lies at a thermal pressure one to three orders of magnitude above that of the bulk of the ISM. In this paper, we propose that this enigmatic component originates from shocks propagating at intermediate ($V_S > 30$ km s$^{-1}$) and high velocities ($V_S \geqslant 100$ km s$^{-1}$) in the Warm Neutral Medium (WNM).
Methods: Shock waves irradiated by the standard interstellar radiation field (ISRF) are modelled using the Paris-Durham shock code designed to follow the dynamical, thermal, and chemical evolutions of shocks with velocities up to 500 km s$^{-1}$. Each observed line of sight is decomposed into a high pressure and a low pressure components. The column density of carbon at high pressure is confronted to the model predictions to derive the number of shocks along the line of sight and their total dissipation rate.
Results: Phase transition shocks spontaneously lead to the presence of high pressure gas in the diffuse ISM and are found to naturally produce neutral carbon with excitation conditions and linewidths in remarkable agreement with the observations. The amounts of neutral carbon at high pressure detected over a sample of 89 lines of sight imply a dissipation rate of mechanical energy with a median of $\sim 3x10^{-25}$ erg cm$^{-3}$ s$^{-1}$ and a dispersion of about a factor of three. This distribution of the dissipation rate weakly depends on the detailed characteristics of shocks as long as they propagate at velocities between 30 and 200 kms s$^{-1}$ in a medium with a preshock density $n_H^0 \ge 0.3$ cm s$^{-3}$ and a transverse magnetic field $B_0 \leqslant 3$ $μ$G. We not only show that this solution is consistent with a scenario of shocks driven by supernovae remnants (SNR) but also that this scenario is, in fact, unavoidable. Any line of sight in the observational sample is bound to intercept SNRs, mostly distributed in the spiral arms of the Milky Way, and expanding in the diffuse ionized and neutral phases of the Galaxy. Surprisingly, the range of dissipation rate derived here, in events that probably drive turbulence in the WNM, is found to be comparable to the distribution of the kinetic energy transfer rate of the turbulent cascade derived from the observations of CO in the Cold Neutral Medium (CNM).
Conclusions: This work reveals a possible direct tracer of the mechanisms by which mechanical energy is injected in the ISM. It also suggests that a still unknown connection exists between the amount of energy dissipated during the injection process in the WNM and that used to feed interstellar turbulence and the turbulent cascade observed in the CNM.
△ Less
Submitted 28 June, 2024;
originally announced June 2024.
-
Radiative and mechanical energies in galaxies I. Contributions of molecular shocks and PDRs in 3C 326 N
Authors:
J. A. Villa-Vélez,
B. Godard,
P. Guillard,
G. Pineau des Forêts
Abstract:
Context: Atomic and molecular lines in galaxies offer insights into energy budgets and feedback mechanisms. Aims: This study establishes a new framework for interpreting these lines and deducing energy budgets from observations. Methods: Atomic and molecular lines detected in a given object are assumed to result from the combination of distributions of shocks and photo-dissociation regions (PDR).…
▽ More
Context: Atomic and molecular lines in galaxies offer insights into energy budgets and feedback mechanisms. Aims: This study establishes a new framework for interpreting these lines and deducing energy budgets from observations. Methods: Atomic and molecular lines detected in a given object are assumed to result from the combination of distributions of shocks and photo-dissociation regions (PDR). Using the Paris-Durham shock code and the Meudon PDR code, emissions are computed over a wide range of parameters. Total emissions are calculated using probability distribution functions, with a defined distance metric based on observed and predicted intensity ratios. Results: We analyze the radio galaxy 3C 326 N, finding both shocks and PDRs necessary to explain the line fluxes. Viable solutions occur only at low densities ($\rm n_H < 100 cm^{-3}$), indicating emission from diffuse interstellar matter. The optimal solution involves low-velocity shocks (5-20 km/s) in PDRs illuminated by UV radiation ten times stronger than in the solar neighborhood. The H$2$ 0-0 S(0) $28 μ$m, [CII] $158 μ$m, and [OI] $63 μ$m lines originate from PDRs, while other H$2$ lines mostly come from shocks. The reprocessed radiative and mechanical energies are $\rm {L_{UV} = 6.3\times10^9 L\odot}$ and $\rm {L_K = 3.9\times10^8 L_\odot}$, respectively, in agreement with 3C 326 N's infrared luminosity, and consistent with 1% of the AGN jet kinetic power dissipated in the interstellar medium. Conclusions: This study demonstrates that the radiative and mechanical energy budgets of galaxies can be derived from observations of atomic and molecular lines alone. It highlights the unexpected importance of the diffuse medium for 3C 326 N. Comparison with new JWST data for 3C 326 N shows striking agreement, opening new prospects for predicting and interpreting extragalactic observations.
△ Less
Submitted 3 May, 2024;
originally announced May 2024.
-
Shocks in the warm neutral medium I -- Theoretical model
Authors:
Benjamin Godard,
Guillaume Pineau Des Forêts,
Shmuel Bialy
Abstract:
Context. Atomic and molecular line emissions from shocks may provide valuable information on the injection of mechanical energy in the interstellar medium (ISM), the generation of turbulence, and the processes of phase transition between the Warm Neutral Medium (WNM) and the Cold Neutral Medium (CNM).Aims. In this series of papers, we investigate the properties of shocks propagating in the WNM. Ou…
▽ More
Context. Atomic and molecular line emissions from shocks may provide valuable information on the injection of mechanical energy in the interstellar medium (ISM), the generation of turbulence, and the processes of phase transition between the Warm Neutral Medium (WNM) and the Cold Neutral Medium (CNM).Aims. In this series of papers, we investigate the properties of shocks propagating in the WNM. Our objective is to identify the tracers of these shocks, use them to interpret ancillary observations of the local diffuse matter, and provide predictions for future observations.Methods. Shocks propagating in the WNM are studied using the Paris-Durham shock code, a multi-fluid model built to follow the thermodynamical and chemical structures of shock waves, at steady-state, in a plane-parallel geometry. The code, already designed to take into account the impact of an external radiation field, is updated to treat self-irradiated shocks at intermediate (30
△ Less
Submitted 3 July, 2024; v1 submitted 30 April, 2024;
originally announced April 2024.
-
Shock excitation of H$_2$ in the James Webb Space Telescope era
Authors:
L. E. Kristensen,
B. Godard,
P. Guillard,
A. Gusdorf,
G. Pineau des Forets
Abstract:
(Abridged) H2 is the most abundant molecule in the Universe. Thanks to its widely spaced energy levels, it predominantly lights up in warm gas, T > 100 K, such as shocked regions, and it is one of the key targets of JWST observations. These include shocks from protostellar outflows, all the way up to starburst galaxies and AGN. Shock models are able to simulate H2 emission. We aim to explore H2 ex…
▽ More
(Abridged) H2 is the most abundant molecule in the Universe. Thanks to its widely spaced energy levels, it predominantly lights up in warm gas, T > 100 K, such as shocked regions, and it is one of the key targets of JWST observations. These include shocks from protostellar outflows, all the way up to starburst galaxies and AGN. Shock models are able to simulate H2 emission. We aim to explore H2 excitation using such models, and to test over which parameter space distinct signatures are produced in H2 emission. We present simulated H2 emission using the Paris-Durham shock code over an extensive grid of 14,000 plane-parallel stationary shock models, a large subset of which are exposed to an external UV radiation field. The grid samples 6 input parameters: preshock density, shock velocity, transverse magnetic field strength, UV radiation field strength, cosmic-ray-ionization rate, and PAH abundance. Physical quantities, such as temperature, density, and width, have been extracted along with H2 integrated line intensities. The strength of the transverse magnetic field, set by the scaling factor, b, plays a key role in the excitation of H2. At low values of b (<~ 0.3, J-type shocks), H2 excitation is dominated by vibrationally excited lines; at higher values (b >~ 1, C-type shocks), rotational lines dominate the spectrum for shocks with an external radiation field comparable to (or lower than) the solar neighborhood. Shocks with b >= 1 can be spatially resolved with JWST for nearby objects. When the input kinetic energy flux increases, the excitation and integrated intensity of H2 increases similarly. An external UV field mainly serves to increase the excitation, particularly for shocks where the input radiation energy is comparable to the input kinetic energy flux. These results provide an overview of the energetic reprocessing of input kinetic energy flux and the resulting H2 line emission.
△ Less
Submitted 9 July, 2023;
originally announced July 2023.
-
3D chemical structure of the diffuse turbulent ISM II -- Origin of CH$^+$, new solution to an 80 years mystery
Authors:
Benjamin Godard,
Guillaume Pineau Des Forêts,
Patrick Hennebelle,
Elena Bellomi,
Valeska Valdivia
Abstract:
Aims: The large abundances of CH$^+$ in the diffuse interstellar medium (ISM) are a long standing issue of our understanding of the thermodynamical and chemical states of the gas. We investigate, here, the formation of CH+ in turbulent and multiphase environments, where the heating of the gas is almost solely driven by the photoelectric effect. Methods: The diffuse ISM is simulated using the magne…
▽ More
Aims: The large abundances of CH$^+$ in the diffuse interstellar medium (ISM) are a long standing issue of our understanding of the thermodynamical and chemical states of the gas. We investigate, here, the formation of CH+ in turbulent and multiphase environments, where the heating of the gas is almost solely driven by the photoelectric effect. Methods: The diffuse ISM is simulated using the magnetohydrodynamic (MHD) code RAMSES which self-consistently computes the dynamical and thermal evolution of the gas along with the time-dependent evolutions of the abundances of H$^+$, H, and H$_2$. The rest of the chemistry, including the abundance of CH$^+$, is computed in post-processing, at equilibrium, under the constraint of out-ofequilibrium of H$^+$, H, and H$_2$. The comparison with the observations is performed taking into account an often neglected, yet paramount, piece of information, namely the length of the intercepted diffuse matter along the observed lines of sight. Results: The quasi totality of the mass of CH$^+$ originates from the unstable gas, in environments where the kinetic temperature is larger than 600 K, the density ranges between 0.6 and 10 cm$^{-3}$, the electronic fraction ranges between 3 x 10$^{-4}$ and 6 x 10$^{-3}$, and the molecular fraction is smaller than 0.4. Its formation is driven by warm and out-of-equilibrium H$_2$ initially formed in the cold neutral medium (CNM) and injected in more diffuse environments and even the warm neutral medium (WNM) through a combination of advection and thermal instability. The simulation which displays the tightest agreement with the HI-to-H$_2$ transition and the thermal pressure distribution observed in the Solar Neighborhood is found to naturally reproduce the observed abundances of CH$^+$, the dispersion of observations, the probability of occurrence of most of the lines of sight, the fraction of non-detections of CH$^+$, and the distribution of its line profiles. The amount of CH$^+$ and the statistical properties of the simulated lines of sight are set by the fraction of unstable gas rich in H$_2$ which is controlled, on Galactic scales, by the mean density of the diffuse ISM (or, equivalently, its total mass), the amplitude of the mean UV radiation field, and the strength of the turbulent forcing. Conclusions: This work offers a new and natural solution to an 80 years old chemical riddle. The almost ubiquitous presence of CH$^+$ in the diffuse ISM likely results from the exchanges of matter between the CNM and the WNM induced by the combination of turbulent advection and thermal instability, without the need to invoke ambipolar diffusion or regions of intermittent turbulent dissipation. Through two phase turbulent mixing, CH$^+$ might thus be a tracer of the H$_2$ mass loss rate of CNM clouds.
△ Less
Submitted 21 September, 2022;
originally announced September 2022.
-
Wide-angle protostellar outflows driven by narrow jets in stratified cores
Authors:
M. Rabenanahary,
S. Cabrit,
Z. Meliani,
G. Pineau des Forêts
Abstract:
Most simulations of outflow feedback on star formation are based on the assumption that outflows are driven by a wide angle "X-wind," rather than a narrow jet. However, the arguments initially raised against pure jet-driven flows were based on steady ejection in a uniform medium, a notion that is no longer supported based on recent observations. We aim to determine whether a pulsed narrow jet laun…
▽ More
Most simulations of outflow feedback on star formation are based on the assumption that outflows are driven by a wide angle "X-wind," rather than a narrow jet. However, the arguments initially raised against pure jet-driven flows were based on steady ejection in a uniform medium, a notion that is no longer supported based on recent observations. We aim to determine whether a pulsed narrow jet launched in a density-stratified, self-gravitating core could reproduce typical molecular outflow properties, without the help of a wide-angle wind component. We performed axisymmetric hydrodynamic simulations using the MPI-AMRVAC code with optically thin radiative cooling on timescales up to 10000 yrs. Then we computed and compared the predicted properties with observational data. First, the jet-driven shell expands faster and wider through a core with steeply decreasing density than through an uniform core. Second, when blown into the same singular flattened core, a jet-driven shell has a similar width as a wide-angle wind-driven shell in the first few hundred years, but a decelerating expansion on long timescales. The flow adopts a conical shape and a base opening angle reaching up to $90\unicode{xb0}$. Third, after $\sim$ 10000 yrs, a pulsed jet-driven shell shows fitting features and a qualitative resemblance with recent observations of protostellar outflows with the Atacama Large Millimeter Array (ALMA), such as HH46-47 and CARMA-7. In particular, similarities are seen in the shell widths, opening angles, position-velocity diagrams, and mass-velocity distribution, with some showing a closer resemblance than in simulations based on a wide-angle "X-wind" model. Therefore, a realistic ambient density stratification in addition to millenia-long integration times are equally essential to reliably predict the properties of outflows driven by a pulsed jet and to confront them with the observations.
△ Less
Submitted 12 April, 2022;
originally announced April 2022.
-
Gas condensation in Brightest Group Galaxies unveiled with MUSE
Authors:
V. Olivares,
P. Salome,
S. L. Hamer,
F. Combes,
M. Gaspari,
K. Kolokythas,
E. O'Sullivan,
R. S. Beckmann,
A. Babul,
F. L. Polles,
M. Lehnert,
S. I. Loubser,
M. Donahue,
M. -L. Gendron-Marsolais,
P. Lagos,
G. Pineau des Forets,
B. Godard,
T. Rose,
G. Tremblay,
G. Ferland,
P. Guillard
Abstract:
The origin of the cold gas in central galaxies in groups is still a matter of debate. We present Multi-Unit Spectroscopic Explorer (MUSE) observations of 18 optically selected local Brightest Group Galaxies (BGGs) to study the kinematics and distribution of the optical emission-line gas. MUSE observations reveal a distribution of gas morphologies including ten complex networks of filaments extendi…
▽ More
The origin of the cold gas in central galaxies in groups is still a matter of debate. We present Multi-Unit Spectroscopic Explorer (MUSE) observations of 18 optically selected local Brightest Group Galaxies (BGGs) to study the kinematics and distribution of the optical emission-line gas. MUSE observations reveal a distribution of gas morphologies including ten complex networks of filaments extending up to 10 kpc to two compact (<3 kpc) and five extended (>5 kpc) disk-dominated structures. Some rotating disks show rings and elongated structures arising from the central disk. The kinematics of the stellar component is mainly rotation-dominated, which is very different from the disturbed kinematics and distribution found in the filamentary sources. The ionized gas is kinematically decoupled from the stellar component for most systems, suggesting an external origin for the gas. We also find that the Halpha luminosity correlates with the cold molecular mass. By exploring the thermodynamical properties of the hot atmospheres, we find that the filamentary sources and compact disks are found in systems with small central entropy values and tcool/teddy ratios. This suggests that, like for Brightest Cluster Galaxies in cool core clusters, the ionized gas are likely formed from hot halo gas condensations, consistently with the Chaotic Cold Accretion simulations (as shown via the C-ratio, Tat, and k-plot). We note that gaseous rotating disks are more frequent than in BCGs. An explanation for the origin of the gas in those objects is a contribution to gas fueling by mergers or group satellites, as qualitatively hinted by some sources of the present sample. Nonetheless, we discuss the possibility that some extended disks could also be a transition stage in an evolutionary sequence including filaments, extended disks and compact disks, as described by hot gas condensation models of cooling flows.
△ Less
Submitted 19 January, 2022;
originally announced January 2022.
-
The ortho-to-para ratio of water in interstellar clouds
Authors:
A. Faure,
P. Hily-Blant,
C. Rist,
G. Pineau des Forêts,
A. Matthews,
D. R. Flower
Abstract:
The nuclear-spin chemistry of interstellar water is investigated using the University of Grenoble Alpes Astrochemical Network (UGAN). This network includes reactions involving the different nuclear-spin states of the hydrides of carbon, nitrogen, oxygen and sulphur, as well as their deuterated forms. Nuclear-spin selection rules are implemented within the scrambling hypothesis for reactions involv…
▽ More
The nuclear-spin chemistry of interstellar water is investigated using the University of Grenoble Alpes Astrochemical Network (UGAN). This network includes reactions involving the different nuclear-spin states of the hydrides of carbon, nitrogen, oxygen and sulphur, as well as their deuterated forms. Nuclear-spin selection rules are implemented within the scrambling hypothesis for reactions involving up to seven protons. The abundances and ortho-to-para ratios (OPRs) of gas-phase water and water ions (H$_2$O$^+$ and H$_3$O$^+$) are computed under the steady-state conditions representative of a dark molecular cloud and during the early phase of gravitational collapse of a prestellar core. The model incorporates the freezing of the molecules on to grains, simple grain surface chemistry and cosmic-ray induced and direct desorption of ices. The predicted OPRs are found to deviate significantly from both thermal and statistical values and to be independent of temperature below $\sim $30~K. The OPR of H$_2$O is shown to lie between 1.5 and 2.6, depending on the spin-state of H$_2$, in good agreement with values derived in translucent clouds with relatively high extinction. In the prestellar core collapse calculations, the OPR of H$_2$O is shown to reach the statistical value of 3 in regions with severe depletion ($n_{\rm H}> 10^7$~cm$^{-3}$). We conclude that a low water OPR ($\lesssim 2.5$) is consistent with gas-phase ion-neutral chemistry and reflects a gas with OPR(H$_2)\lesssim 1$. Available OPR measurements in protoplanetary disks and comets are finally discussed.
△ Less
Submitted 6 January, 2022;
originally announced January 2022.
-
Sulfur gas-phase abundance in dense cores
Authors:
P. Hily-Blant,
G. Pineau des Forêts,
A. Faure,
F. Lique
Abstract:
The abundance of volatile sulfur in dense clouds is long-standing problem in studies of the physics and chemistry of star-forming regions. Sulfur is an important species because its low ionization potential may possibly make it an important charge carrier. The observed sulfur-bearing species in the gas-phase of dense clouds represent only a minor fraction of the cosmic sulfur abundance, which has…
▽ More
The abundance of volatile sulfur in dense clouds is long-standing problem in studies of the physics and chemistry of star-forming regions. Sulfur is an important species because its low ionization potential may possibly make it an important charge carrier. The observed sulfur-bearing species in the gas-phase of dense clouds represent only a minor fraction of the cosmic sulfur abundance, which has been interpreted as a signature of sulfur depletion into ices at the surface of dust grains. However, atomic sulfur, which could be the main gas-phase carrier, cannot be observed directly in cold cores. We present measurements of the nitrogen sulfide (NS) radical toward four dense cores performed with the IRAM-30m telescope. Analytical chemical considerations and chemical models over a wide parameter space show that the NS:N2H+ abundance ratio provides a direct constraint on the abundance of gas-phase atomic sulfur. Toward early-type cores, we find that $n(\rm S)/n_{\rm H}$ is close, or even equal, to the cosmic abundance of sulfur, 14$\times 10^{-6}$, demonstrating that sulfur is not depleted and is atomic, which is in agreement with chemical models. More chemically evolved cores show sulfur depletion by factors up to 100 in their densest parts. In L1544, atomic sulfur depletion is shown to increase with increasing density. Future observations are needed to discover the solid-phase carrier of sulfur. The initial steps of the collapse of pre-stellar cores in the high sulfur abundance regime also need to be explored from their chemical and dynamical perspectives.
△ Less
Submitted 5 December, 2021; v1 submitted 2 December, 2021;
originally announced December 2021.
-
Self-generated ultraviolet radiation in molecular shock waves II. CH+ and the interpretation of emission from shock ensembles
Authors:
Andrew Lehmann,
Benjamin Godard,
Gillaume Pineau des Forêts,
Alba Vidal-García,
Edith Falgarone
Abstract:
Shocks, modelled over a broad range of parameters, are used to construct a new tool to deduce the mechanical energy and physical conditions from observed atomic or molecular emission lines. We compute magnetised, molecular shock models with velocities $V_s=5$-$80$ km s$^{-1}$, preshock proton densities $n_{\rm H}=10^2$-$10^6$ cm$^{-3}$, weak or moderate magnetic field strengths, and in the absence…
▽ More
Shocks, modelled over a broad range of parameters, are used to construct a new tool to deduce the mechanical energy and physical conditions from observed atomic or molecular emission lines. We compute magnetised, molecular shock models with velocities $V_s=5$-$80$ km s$^{-1}$, preshock proton densities $n_{\rm H}=10^2$-$10^6$ cm$^{-3}$, weak or moderate magnetic field strengths, and in the absence or presence of an external UV radiation field. We develop a simple emission model of an ensemble of shocks for connecting any observed emission lines to the mechanical energy and physical conditions of the system. For this range of parameters we find the full diversity (C-, C$^*$-, CJ-, and J-type) of magnetohydrodynamic shocks. H$_2$ and H are dominant coolants, with up to 30% of the shock kinetic flux escaping in Ly$α$ photons. The reformation of molecules in the cooling tail means H$_2$ is even a good tracer of dissociative shocks and shocks that were initially fully atomic. For each shock model we provide integrated intensities of rovibrational lines of H$_2$, CO, and CH$^+$, atomic H lines, and atomic fine-structure and metastable lines. We demonstrate how to use these shock models to deduce the mechanical energy and physical conditions of extragalactic environments. As a template example, we interpret the CH$^+$(1-0) emission from the Eyelash starburst galaxy. A mechanical energy injection rate of at least $10^{11}$ $L_\odot$ into molecular shocks is required to reproduce the observed line. The low-velocity, externally irradiated shocks are at least an order magnitude more efficient than the most efficient shocks with no external irradiation, in terms of the total mechanical energy required. We predict differences of more than 2 orders of magnitude in intensities of the pure rotational lines of CO, Ly$α$, metastable lines of O, S$^+$, and N, between representative models.
△ Less
Submitted 28 November, 2021;
originally announced November 2021.
-
Excitation mechanisms in the intracluster filaments surrounding Brightest Cluster Galaxies
Authors:
F. L. Polles,
P. Salomé,
P. Guillard,
B. Godard,
G. Pineau des Forêts,
V. Olivares,
R. S. Beckmann,
R. E. A. Canning,
F. Combes,
Y. Dubois,
A. C. Edge,
A. C. Fabian,
G. J. Ferland,
S. L. Hamer,
M. D. Lehnert
Abstract:
The excitation of the filamentary gas structures surrounding giant elliptical galaxies at the center of cool-core clusters, a.k.a BCGs (brightest cluster galaxies), is key to our understanding of active galactic nucleus feedback, and of the impact of environmental and local effects on star formation. We investigate the contribution of the thermal radiation from the cooling flow surrounding BCGs to…
▽ More
The excitation of the filamentary gas structures surrounding giant elliptical galaxies at the center of cool-core clusters, a.k.a BCGs (brightest cluster galaxies), is key to our understanding of active galactic nucleus feedback, and of the impact of environmental and local effects on star formation. We investigate the contribution of the thermal radiation from the cooling flow surrounding BCGs to the excitation of the filaments. We explore the effects of small levels of extra-heating (turbulence), and of metallicity, on the optical and infrared lines. Using the Cloudy code, we model the photoionization and photodissociation of a slab of gas of optical depth AV{\leq}30mag at constant pressure, in order to calculate self-consistently all of the gas phases, from ionized gas to molecular gas. The ionizing source is the EUV and soft X-ray radiation emitted by the cooling gas. We test these models comparing their predictions to the rich multi-wavelength observations, from optical to submillimeter. These models reproduce most of the multi-wavelength spectra observed in the nebulae surrounding the BCGs, not only the LINER-like optical diagnostics: [O iii]λ 5007 Å/H\b{eta}, [N ii]λ 6583 Å/Hα and ([S ii]λ 6716 Å+[S ii]λ 6731 Å)/Hα but also the infrared emission lines from the atomic gas. The modeled ro-vib H2 lines also match observations, which indicates that near and mid-IR H2 lines are mostly excited by collisions between H2 molecules and secondary electrons produced naturally inside the cloud by the interaction between the X-rays and the cold gas in the filament. However, there is still some tension between ionized and molecular line tracers (i.e. CO), which requires to optimize the cloud structure and the density of the molecular zone.
△ Less
Submitted 17 March, 2021;
originally announced March 2021.
-
Extremely broad Lyman-alpha line emission from the molecular intra-group medium in Stephan's Quintet: evidence for a turbulent cascade in a highly clumpy multi-phase medium?
Authors:
P. Guillard,
P. N Appleton,
F. Boulanger,
J. M. Shull,
M. D. Lehnert,
G. Pineau des Forets,
E. Falgarone,
M. E. Cluver,
C. K. Xu,
S. C. Gallagher,
P. A. Duc
Abstract:
We present Hubble Space Telescope Cosmic Origin Spectrograph (COS) UV line spectroscopy and integral-field unit (IFU) observations of the intra-group medium in Stephan's Quintet (SQ). SQ hosts a 30 kpc long shocked ridge triggered by a galaxy collision at a relative velocity of 1000 km/s, where large amounts of molecular gas coexist with a hot, X-ray emitting, plasma. COS spectroscopy at five posi…
▽ More
We present Hubble Space Telescope Cosmic Origin Spectrograph (COS) UV line spectroscopy and integral-field unit (IFU) observations of the intra-group medium in Stephan's Quintet (SQ). SQ hosts a 30 kpc long shocked ridge triggered by a galaxy collision at a relative velocity of 1000 km/s, where large amounts of molecular gas coexist with a hot, X-ray emitting, plasma. COS spectroscopy at five positions sampling the diverse environments of the SQ intra-group medium reveals very broad (2000 km/s) Ly$α$ line emission with complex line shapes. The Ly$α$ line profiles are similar to or much broader than those of H$β$, [CII]$\lambda157.7μ$m and CO~(1-0) emission. The extreme breadth of the Ly$α$ emission, compared with H$β$, implies resonance scattering within the observed structure. Scattering indicates that the neutral gas of the intra-group medium is clumpy, with a significant surface covering factor. We observe significant variations in the Ly$α$/H$β$ flux ratio between positions and velocity components. From the mean line ratio averaged over positions and velocities, we estimate the effective escape fraction of Ly$α$ photons to be 10-30%. Remarkably, over more than four orders of magnitude in temperature, the powers radiated by X-rays, Ly$α$, H$_2$, [CII] are comparable within a factor of a few, assuming that the ratio of the Ly$α$ to H$_2$ fluxes over the whole shocked intra-group medium stay in line with those observed at those five positions. Both shocks and mixing layers could contribute to the energy dissipation associated with a turbulent energy cascade. Our results may be relevant for the cooling of gas at high redshifts, where the metal content is lower than in this local system, and a high amplitude of turbulence is more common.
△ Less
Submitted 20 October, 2021; v1 submitted 12 February, 2021;
originally announced February 2021.
-
Self-generated ultraviolet radiation in molecular shock waves I. Effects of Lyman $α$, Lyman $β$, and two-photon continuum
Authors:
A. Lehmann,
B. Godard,
G. Pineau des Forêts,
E. Falgarone
Abstract:
Shocks are ubiquitous in the interstellar and intergalactic media, where their chemical and radiative signatures reveal the physical conditions in which they arise. Detailed astrochemical models of shocks at all velocities are necessary to understand the physics of many environments including protostellar outflows, supernova remnants, and galactic outflows.
We present an accurate treatment of th…
▽ More
Shocks are ubiquitous in the interstellar and intergalactic media, where their chemical and radiative signatures reveal the physical conditions in which they arise. Detailed astrochemical models of shocks at all velocities are necessary to understand the physics of many environments including protostellar outflows, supernova remnants, and galactic outflows.
We present an accurate treatment of the self-generated UV radiation in intermediate velocity, stationary, weakly magnetised, J-type, molecular shocks. Shock solutions computed with the Paris-Durham shock code are post-processed using a multi-level accelerated $Λ$-iteration radiative transfer algorithm to compute Ly$α$, Ly$β$, and 2-photon continuum emission. The subsequent impacts on the ionisation and dissociation of key atomic and molecular species as well as on the heating by the photoelectric effect take the wavelength dependent cross-sections and the fluid velocity profile into account. We analyse shock models with velocities $V=25-60$ km/s, propagating in dense ($n \geq 10^4$ ${\rm cm}^{-3}$), shielded gas.
Self-absorption traps Ly$α$ photons in a small region in the shock, though a large fraction escapes into the line wings. We find a critical velocity $V\sim 30$ km/s above which shocks produce a Ly$α$ photon flux exceeding that of the standard ISRF. The escaping photons generate a warm slab (T~100 K) ahead of the shock as well as pre-ionise C and S. These shocks are traced by bright atomic fine structure (e.g. O and S) and metastable (e.g. O and C) lines, substantive molecular emission (e.g. H2, OH, and CO), enhanced column densities of several species (e.g. CH+ and HCO+), as well as a severe destruction of H2O. As much as 13-21% of the initial kinetic energy of the shock escapes in Ly$α$ and Ly$β$ photons if the dust opacity in the radiative precursor allows it.
△ Less
Submitted 2 October, 2020;
originally announced October 2020.
-
Depletion and fractionation of nitrogen in collapsing cores
Authors:
P. Hily-Blant,
G. Pineau des Forêts,
A. Faure,
D. R. Flower
Abstract:
Measurements of the nitrogen isotopic ratio in Solar System comets show a constant value, ~140, which is three times lower than the protosolar ratio, a highly significant difference that remains unexplained. Observations of static starless cores at early stages of collapse confirm the theoretical expectation that nitrogen fractionation in interstellar conditions is marginal for most species. Yet,…
▽ More
Measurements of the nitrogen isotopic ratio in Solar System comets show a constant value, ~140, which is three times lower than the protosolar ratio, a highly significant difference that remains unexplained. Observations of static starless cores at early stages of collapse confirm the theoretical expectation that nitrogen fractionation in interstellar conditions is marginal for most species. Yet, observed isotopic ratios in N2H+ are at variance with model predictions. These gaps in our understanding of how the isotopic reservoirs of nitrogen evolve, from interstellar clouds to comets, and, more generally, to protosolar nebulae, may have their origin in missing processes or misconceptions in the chemistry of interstellar nitrogen. So far, theoretical studies of nitrogen fractionation in starless cores have addressed the quasi-static phase of their evolution such that the effect of dynamical collapse on the isotopic ratio is not known. In this paper, we investigate the fractionation of 14N and 15N during the gravitational collapse of a pre-stellar core through gas-phase and grain adsorption and desorption reactions. The initial chemical conditions, which are obtained in steady state after typically a few Myr, show low degrees of fractionation in the gas phase, in agreement with earlier studies. However, during collapse, the differential rate of adsorption of 14N- and 15N-containing species onto grains results in enhanced 15N:14N ratios, in better agreement with the observations. Furthermore, we find differences in the behavior, with increasing density, of the isotopic ratio in different species. We find that the collapse must take place on approximately one free-fall timescale, based on the CO abundance profile in L183 [see the end in the PDF file]
△ Less
Submitted 14 September, 2020;
originally announced September 2020.
-
3D chemical structure of diffuse turbulent ISM. I. Statistics of the HI-to-H$_2$ transition
Authors:
Elena Bellomi,
Benjamin Godard,
Patrick Hennebelle,
Valeska Valdivia,
Guillaume Pineau des Forêts,
Pierre Lesaffre,
Michel Pérault
Abstract:
We studied the statistical properties of the HI-to-H$_2$ transition observed in absorption in the local diffuse and multiphase ISM to identify the physical processes controlling the probability of occurrence of any line of sight. The turbulent diffuse ISM is modeled using the RAMSES code, which includes detailed treatments of the magnetohydrodynamics, the thermal evolution of the gas, and the chem…
▽ More
We studied the statistical properties of the HI-to-H$_2$ transition observed in absorption in the local diffuse and multiphase ISM to identify the physical processes controlling the probability of occurrence of any line of sight. The turbulent diffuse ISM is modeled using the RAMSES code, which includes detailed treatments of the magnetohydrodynamics, the thermal evolution of the gas, and the chemistry of H$_2$. The impacts of the UV radiation field, the mean density, the turbulent forcing, the integral scale, the magnetic field, and the gravity on the molecular content of the gas are explored through a parametric study covering a wide range of physical conditions. The statistics of the HI-to-H$_2$ transition are interpreted through analytical prescriptions and compared with the observations using a modified and robust version of the Kolmogorov-Smirnov test. The results of one simulation, convolved with the distribution of distances of the observational sample, are able to explain most of the statistical properties of the HI-to-H$_2$ transition observed in the local ISM. The tightest agreement is obtained for a neutral diffuse gas modeled over ~200 pc, with a mean density of $1-2$ cm$^{-3}$, illuminated by the standard interstellar UV radiation field, and stirred up by a large-scale compressive turbulent forcing. Within this configuration, the 2D probability histogram of the column densities of H and H$_2$ is remarkably stable and is almost unaltered by gravity, the strength of the turbulent forcing, the resolution of the simulation, or the strength of the magnetic field $B_x$. The weak effect of the resolution and our analytical prescription suggest that the column densities of HI are likely built up in large-scale WNM and CNM structures correlated in density over ~20 pc and ~10 pc, respectively, while those of H$_2$ are built up in CNM structures between ~3 pc and ~10 pc.
△ Less
Submitted 25 October, 2020; v1 submitted 11 September, 2020;
originally announced September 2020.
-
Dust coagulation feedback on magnetohydrodynamic resistivities in protostellar collapse
Authors:
V. Guillet,
P. Hennebelle,
G. Pineau des Forêts,
A. Marcowith,
B. Commerçon,
P. Marchand
Abstract:
The degree of coupling between the gas and the magnetic field during the collapse of a core and the subsequent formation of a disk depends on the assumed dust size distribution. We study the impact of grain-grain coagulation on the evolution of magnetohydrodynamic (MHD) resistivities during the collapse of a prestellar core. We use a 1-D model to follow the evolution of the dust size distribution,…
▽ More
The degree of coupling between the gas and the magnetic field during the collapse of a core and the subsequent formation of a disk depends on the assumed dust size distribution. We study the impact of grain-grain coagulation on the evolution of magnetohydrodynamic (MHD) resistivities during the collapse of a prestellar core. We use a 1-D model to follow the evolution of the dust size distribution, out-of-equilibrium ionization state and gas chemistry during the collapse of a prestellar core. To compute the grain-grain collisional rate, we consider models for both random and systematic, size-dependent, velocities. We include grain growth through grain-grain coagulation and ice accretion, but ignore grain fragmentation. Starting with a MRN (Mathis et al. 1977) size distribution, we find that coagulation in grain-grain collisions generated by hydrodynamical turbulence is not efficient at removing the smallest grains, and as a consequence does not affect much the evolution of the Hall and ambipolar diffusion MHD resistivities which still severly drop during the collapse like in models without coagulation. The inclusion of systematic velocities, possibly induced by the presence of ambipolar diffusion, increases the coagulation rate between small and large grains, removing small grains earlier in the collapse and therefore limiting the drop in the Hall and ambipolar diffusion resistivities. At intermediate densities ($n_{\rm H} \sim 10^8\,$cm$^{-3}$), the Hall and ambipolar diffusion resistivities are found to be higher by 1 to 2 orders of magnitude in models with coagulation than in models where coagulation is ignored, and also higher than in a toy model without coagulation where all grains smaller than $0.1\,μ$m would have been removed in the parent cloud before the collapse. When grain drift velocities induced by ambipolar diffusion are included, dust coagulation ... (abridged)
△ Less
Submitted 8 July, 2020;
originally announced July 2020.
-
Constraining MHD disk winds with ALMA. Apparent rotation signatures and application to HH212
Authors:
B. Tabone,
S. Cabrit,
G. Pineau des Forêts,
J. Ferreira,
A. Gusdorf,
L. Podio,
E. Bianchi,
E. Chapillon,
C. Codella,
F. Gueth
Abstract:
Large millimeter interferometers are revealing a growing number of rotating outflows, which are suggested to trace magneto-centrifugal disk winds (MHD DWs). However, their impact on disk accretion is not yet well quantified. Here we identify systematic biases in retrieving the true launch zone, magnetic lever arm, and angular momentum flux of an MHD DW from apparent rotation signatures. Synthetic…
▽ More
Large millimeter interferometers are revealing a growing number of rotating outflows, which are suggested to trace magneto-centrifugal disk winds (MHD DWs). However, their impact on disk accretion is not yet well quantified. Here we identify systematic biases in retrieving the true launch zone, magnetic lever arm, and angular momentum flux of an MHD DW from apparent rotation signatures. Synthetic position-velocity cuts are constructed from self-similar MHD DWs over a broad range of parameters, and three different methods are applied for estimating the specific angular momentum. We find that the launch radius inferred using the well-known relation from Anderson et al. (2006) can markedly differ from the true outermost launch radius $r_{out}$ of the DW. The "double-peak separation" and "flow width" methods provide only a strict lower limit to $r_{out}$. This bias is independent of angular resolution and can reach a factor ten. In contrast, the "rotation curve" method gives a good estimate of $r_{out}$ when the flow is well resolved, and an upper limit otherwise. The magnetic lever arm is always underestimated. Only comparison with synthetic predictions can take into account properly all observational effects. As an application, we present a comparison with ALMA observations of HH212 at resolutions from 250 au to 16 au, which represents the most stringent observational test of MHD DW to date. This comparison confirms our predicted biases for the double-peak separation method, and the large $r_{out}\sim40~$au and small magnetic lever arm first suggested by Tabone et al. (2017). We also derive the first accurate analytical expression for the fraction of disk angular momentum extracted by an MHD disk wind of given radial extent, magnetic lever arm, and mass flux. Application to HH212 confirms that MHD DWs are serious candidates for the steady angular momentum extraction process in young disks.
△ Less
Submitted 19 April, 2020;
originally announced April 2020.
-
Molecule formation in dust-poor irradiated jets I. Stationary disk winds
Authors:
B. Tabone,
B. Godard,
G. Pineau des Forêts,
S. Cabrit,
E. F. van Dishoeck
Abstract:
Recent ALMA observations suggest that the highest velocity part of molecular protostellar jets are launched from the dust-sublimation regions of the accretion disks (<0.3 au). However, formation and survival of molecules in inner protostellar disk winds, in the presence of a harsh FUV radiation field and the absence of dust, remain unexplored. We aim at determining if simple molecules can be synth…
▽ More
Recent ALMA observations suggest that the highest velocity part of molecular protostellar jets are launched from the dust-sublimation regions of the accretion disks (<0.3 au). However, formation and survival of molecules in inner protostellar disk winds, in the presence of a harsh FUV radiation field and the absence of dust, remain unexplored. We aim at determining if simple molecules can be synthesized and spared in fast and collimated dust-free disk winds or if a fraction of dust is necessary to explain the observed molecular abundances. This work is based on the Paris-Durham shock code designed to model irradiated environments. Fundamental properties of the dust-free chemistry are investigated from single point models. A laminar 1D disk wind model is then built using a parametric flow geometry. This model includes time-dependent chemistry and the attenuation of the radiation field by gas-phase photoprocesses. We show that a small fraction of H2 (< 1e-2), primarily formed through the H- route, can efficiently initiate molecule synthesis such as CO and SiO above TK ~ 800 K. The attenuation of the radiation field by atomic species (eg. C, Si, S) proceeds through continuum self-shielding. This process ensures efficient formation of CO, OH, SiO, H2O through neutral-neutral reactions, and the survival of these molecules. Class 0 dust-free winds with high mass-loss rates ($\dot{M}_w >$ 2e-6 Msun/yr) are predicted to be rich in molecules if warm (TK > 800 K). The molecular content of disk winds is very sensitive to the presence of dust and a mass-fraction of surviving dust as small as 1e-5 significantly increases the H2O and SiO abundances. Chemistry of high-velocity jets is a powerful tool to probe their content in dust and uncover their launching point. Models of internal shocks are required to fully exploit the current (sub-)millimeter observations and prepare future JWST observations.
△ Less
Submitted 3 March, 2020;
originally announced March 2020.
-
Dense gas formation and destruction in a simulated Perseus-like galaxy cluster with spin-driven black hole feedback
Authors:
R. S. Beckmann,
Y. Dubois,
P. Guillard,
P. Salome,
V. Olivares,
F. Polles,
C. Cadiou,
F. Combes,
S. Hamer,
M. D. Lehnert,
G. Pineau des Forets
Abstract:
Extended filamentary H$α$ emission nebulae are a striking feature of nearby galaxy clusters but the formation mechanism of the filaments, and the processes which shape their morphology remain unclear. We conduct an investigation into the formation, evolution and destruction of dense gas in the center of a simulated, Perseus-like, cluster under the influence of a spin-driven jet. We particularly st…
▽ More
Extended filamentary H$α$ emission nebulae are a striking feature of nearby galaxy clusters but the formation mechanism of the filaments, and the processes which shape their morphology remain unclear. We conduct an investigation into the formation, evolution and destruction of dense gas in the center of a simulated, Perseus-like, cluster under the influence of a spin-driven jet. We particularly study the role played by condensation of dense gas from the diffuse intracluster medium, and the impact of direct uplifting of existing dense gas by the jets, in determining the spatial distribution and kinematics of the dense gas. We present a hydrodynamical simulation of an idealised Perseus-like cluster using the adaptive mesh refinement code {\sc ramses}. Our simulation includes a supermassive black hole (SMBH) that self-consistently tracks its spin evolution via its local accretion, and in turn drives a large-scale jet whose direction is based on the black hole's spin evolution. We show that the formation and destruction of dense gas is closely linked to the SMBH's feedback cycle, and that its morphology is highly variable throughout the simulation. While extended filamentary structures readily condense from the hot intra-cluster medium, they are easily shattered into an overly clumpy distribution of gas during their interaction with the jet driven outflows. Condensation occurs predominantly onto infalling gas located 5 - 15 kpc from the center during quiescent phases of the central AGN, when the local ratio of the cooling time to free fall time falls below 20, i.e. when $t_{\rm cool}/t_{\rm ff} < 20$. We find evidence for both condensation and uplifting of dense gas, but caution that purely hydrodynamical simulations struggle to effectively regulate the cluster cooling cycle and produce overly clumpy distributions of dense gas morphologies, compared to observation.
△ Less
Submitted 17 September, 2019; v1 submitted 3 September, 2019;
originally announced September 2019.
-
Warm H$_2$ as a probe of massive accretion and feedback through shocks and turbulence across cosmic time
Authors:
Philip Appleton,
Lee Armus,
Francois Boulanger,
Charles M. Bradford,
Jonathan Braine,
Volker Bromm,
Peter Capak,
Michelle Cluver,
Asantha Cooray,
Tanio Diaz-Santos,
Eiichi Egami,
Bjorn Emonts,
Pierre Guillard,
George Helou,
Lauranne Lanz,
Susanne Madden,
Anne Medling,
Ewan O'Sullivan,
Patrick Ogle,
Alexandra Pope,
Guillaume Pineau des Forêts,
J. Michael Shull,
John-David Smith,
Aditya Togi,
C. Kevin Xu
Abstract:
Galaxy formation depends on a complex interplay between gravitational collapse, gas accretion, merging, and feedback processes. Yet, after many decades of investigation, these concepts are poorly understood. This paper presents the argument that warm H$_2$ can be used as a tool to unlock some of these mysteries. Turbulence, shocks and outflows, driven by star formation, AGN activity or inflows, ma…
▽ More
Galaxy formation depends on a complex interplay between gravitational collapse, gas accretion, merging, and feedback processes. Yet, after many decades of investigation, these concepts are poorly understood. This paper presents the argument that warm H$_2$ can be used as a tool to unlock some of these mysteries. Turbulence, shocks and outflows, driven by star formation, AGN activity or inflows, may prevent the rapid buildup of star formation in galaxies. Central to our understanding of how gas is converted into stars is the process by which gas can dissipate its mechanical energy through turbulence and shocks in order to cool. H$_2$ lines provide direct quantitative measurements of kinetic energy dissipation in molecular gas in galaxies throughout the Universe. Based on the detection of very powerful H$_2$ lines from z = 2 galaxies and proto-clusters at the detection limits of {\it Spitzer}, we are confident that future far-IR and UV H$_2$ observations will provide a wealth of new information and insight into galaxy evolution to high-z. Finally, at the very earliest epoch of star and galaxy formation, warm H$_2$ may also provide a unique glimpse of molecular gas collapse at 7 $<$ z $<$ 12 in massive dark matter (DM) halos on their way to forming the very first galaxies. Such measurements are beyond the reach of existing and planned observatories.
△ Less
Submitted 15 March, 2019;
originally announced March 2019.
-
Ubiquitous cold and massive filaments in cool core clusters
Authors:
V. Olivares,
P. Salomé,
F. Combes,
S. Hamer,
P. Guillard,
M. D. Lehnert,
F. Polles,
R. S. Beckmann,
Y. Dubois,
M. Donahue,
A. Edge,
A. C. Fabian,
B. McNamara,
T. Rose,
H. Russell,
G. Tremblay,
A. Vantyghem,
R. E. A. Canning,
G. Ferland,
B. Godard,
M. Hogan,
S. Peirani,
G. Pineau des Forets
Abstract:
Multi-phase filamentary structures around Brightest Cluster Galaxies are likely a key step of AGN-feedback. We observed molecular gas in 3 cool cluster cores: Centaurus, Abell S1101, and RXJ1539.5 and gathered ALMA and MUSE data for 12 other clusters. Those observations show clumpy, massive and long, 3--25 kpc, molecular filaments, preferentially located around the radio bubbles inflated by the AG…
▽ More
Multi-phase filamentary structures around Brightest Cluster Galaxies are likely a key step of AGN-feedback. We observed molecular gas in 3 cool cluster cores: Centaurus, Abell S1101, and RXJ1539.5 and gathered ALMA and MUSE data for 12 other clusters. Those observations show clumpy, massive and long, 3--25 kpc, molecular filaments, preferentially located around the radio bubbles inflated by the AGN (Active Galactic Nucleus). Two objects show nuclear molecular disks. The optical nebula is certainly tracing the warm envelopes of cold molecular filaments. Surprisingly, the radial profile of the H$α$/CO flux ratio is roughly constant for most of the objects, suggesting that (i) between 1.2 to 7 times more cold gas could be present and (ii) local processes must be responsible for the excitation. Projected velocities are between 100--400 km s$^{-1}$, with disturbed kinematics and sometimes coherent gradients. This is likely due to the mixing in projection of several thin unresolved filaments. The velocity fields may be stirred by turbulence induced by bubbles, jets or merger-induced sloshing. Velocity and dispersions are low, below the escape velocity. Cold clouds should eventually fall back and fuel the AGN. We compare the filament's radial extent, r$_{fil}$, with the region where the X-ray gas can become thermally unstable. The filaments are always inside the low-entropy and short cooling time region, where t$_{cool}$/t$_{ff}$<20 (9 of 13 sources). The range t$_{cool}$/t$_{ff}$, 8-23 at r$_{fil}$, is likely due to (i) a more complex gravitational potential affecting the free-fall time (e.g., sloshing, mergers); (ii) the presence of inhomogeneities or uplifted gas in the ICM, affecting the cooling time. For some of the sources, r$_{fil}$ lies where the ratio of the cooling time to the eddy-turnover time, t$_{cool}$/t$_{eddy}$, is approximately unity.
△ Less
Submitted 25 February, 2019;
originally announced February 2019.
-
Models of irradiated molecular shocks
Authors:
Benjamin Godard,
Guillaume Pineau des Forêts,
Pierre Lesaffre,
Andrew Lehmann,
Antoine Gusdorf,
Edith Falgarone
Abstract:
Aims. The goal of the paper is to present a detailed study of the propagation of low velocity (5 to 25 km s-1) stationary molecular shocks in environments illuminated by an external ultraviolet (UV) radiation field. In particular, we intend to show how the structure, dynamics, energetics, and chemical properties of shocks are modified by UV photons and to estimate how efficiently shocks can produc…
▽ More
Aims. The goal of the paper is to present a detailed study of the propagation of low velocity (5 to 25 km s-1) stationary molecular shocks in environments illuminated by an external ultraviolet (UV) radiation field. In particular, we intend to show how the structure, dynamics, energetics, and chemical properties of shocks are modified by UV photons and to estimate how efficiently shocks can produce line emission. Methods. We implemented several key physico-chemical processes in the Paris-Durham shock code to improve the treatment of the radiative transfer and its impact on dust and gas particles. We propose a new integration algorithm to find the steady-state solutions of magnetohydrodynamics equations in a range of parameters in which the fluid evolves from a supersonic to a subsonic regime. We explored the resulting code over a wide range of physical conditions, which encompass diffuse interstellar clouds and hot and dense photon-dominated regions (PDR). Results. We find that C-type shock conditions cease to exist as soon as G0 > 0.2 (nH/cm-3)^1/2. Such conditions trigger the emergence of another category of stationary solutions, called C*-type and CJ-type shocks, in which the shocked gas is momentarily subsonic along its trajectory. These solutions are shown to be unique for a given set of physical conditions and correspond to dissipative structures in which the gas is heated up to temperatures comprised between those found in C-type and adiabatic J-type shocks. High temperatures combined with the ambient UV field favour the production or excitation of a few molecular species to the detriment of others, hence leading to specific spectroscopic tracers such as rovibrational lines of H2 and rotational lines of CH+. Unexpectedly, the rotational lines of CH+ may carry as much as several percent of the shock kinetic energy.
△ Less
Submitted 14 January, 2019;
originally announced January 2019.
-
Water and interstellar complex organics associated with the HH 212 protostellar disc - On disc atmospheres, disc winds, and accretion shocks
Authors:
C. Codella,
E. Bianchi,
B. Tabone,
C. -F. Lee,
S. Cabrit,
C. Ceccarelli,
L. Podio,
F. Bacciotti,
R. Bachiller,
E. Chapillon,
F. Gueth,
A. Gusdorf,
B. Lefloch,
S. Leurini,
G. Pineau des Forets,
K. L. J. Rygl,
M. Tafalla
Abstract:
The HH 212 protostellar system, in Orion B, has been mapped thanks to ALMA-Band 7 Cycle 1 and Cycle 4 observations of dueterated water (HDO) and acetaldehyde (CH$_3$CHO) emission with an angular resolution down to $\sim$0.15 arcsec (60 au). Many emission lines due to 14 CH$_3$CHO and 1 HDO transitions at high excitation ($E_{\rm u}$ between 163 K and 335 K) have been imaged in the inner $\sim$ 70…
▽ More
The HH 212 protostellar system, in Orion B, has been mapped thanks to ALMA-Band 7 Cycle 1 and Cycle 4 observations of dueterated water (HDO) and acetaldehyde (CH$_3$CHO) emission with an angular resolution down to $\sim$0.15 arcsec (60 au). Many emission lines due to 14 CH$_3$CHO and 1 HDO transitions at high excitation ($E_{\rm u}$ between 163 K and 335 K) have been imaged in the inner $\sim$ 70 au region. The local thermal equilibrium analysis of the CH$_3$CHO emission leads to a temperature of 78$\pm$14 K and a column density of 7.6$\pm$3.2 $\times$ 10$^{15}$ cm$^{-2}$, which, when $N_{\rm H_2}$ of 10$^{24}$ cm$^{-2}$ is assumed, leads to an abundance of $X_{\rm CH_3CHO}$ $\simeq$ 8 $\times$ 10$^{-9}$. The large velocity gradient analysis of the HDO emission also places severe constraints on the volume density, n$_{\rm H_2}$ $\geq$ 10$^8$ cm$^{-3}$. The line profiles are 5--7 km s$^{-1}$ wide, and CH$_3$CHO and HDO both show a $\pm$ 2 km s$^{-1}$ velocity gradient over a size of $\sim$ 70 au (blue-shifted emission towards the north-west and red-shifted emission towards the south-east) along the disc equatorial plane, in agreement with what was found so far using other molecular tracers. The kinematics of CH$_3$CHO and HDO are consistent with the occurrence of a centrifugal barrier, that is, the infalling envelope-rotating disc ring, which is chemically enriched through low-velocity accretion shocks. The emission radius is $\sim$ 60 au, in good agreement with what was found before for another interstellar complex organic molecule such as NH$_2$CHO. We support a vertical structure for the centrifugal barrier, suggesting the occurrence of two outflowing, expanding, and rotating rings above and below (of about 40-45 au) the optically thick equatorial disc plane. It is tempting to speculate that these rings could probe the basis of a wind launched from this region.
△ Less
Submitted 20 June, 2018;
originally announced June 2018.
-
Herschel Spectroscopy of the Taffy Galaxies (UGC 12914/12915 = VV 254): Enhanced [C II] emission in the collisionally-formed bridge
Authors:
B. W. Peterson,
P. N. Appleton,
T. Bitsakis,
P. Guillard,
K. Alatalo,
F. Boulanger,
M. Cluver,
P. -A. Duc,
E. Falgarone,
S. Gallagher,
Y. Gao,
G. Helou,
T. H. Jarrett,
B. Joshi,
U. Lisenfeld,
N. Lu,
P. Ogle,
G. Pineau des Forêts,
P. van der Werf,
C. K. Xu
Abstract:
Using the PACS and SPIRE spectrometers on-board Herschel, we obtained observations of the Taffy galaxies (UGC 12914/12915) and bridge. The Taffy system is believed to be the result of a face-on collision between two gas-rich galaxies, in which the stellar disks passed through each other, but the gas was dispersed into a massive H I and molecular bridge between them. Emission is detected and mapped…
▽ More
Using the PACS and SPIRE spectrometers on-board Herschel, we obtained observations of the Taffy galaxies (UGC 12914/12915) and bridge. The Taffy system is believed to be the result of a face-on collision between two gas-rich galaxies, in which the stellar disks passed through each other, but the gas was dispersed into a massive H I and molecular bridge between them. Emission is detected and mapped in both galaxies and the bridge in the [C II]157.7 $μ$m and [O I]63.2 $μ$m fine-structure lines. Additionally, SPIRE FTS spectroscopy detects the [C I] $^3$P$_2$$\rightarrow$$^3$P$_1$(809.3 GHz) and [C I] $^3$P$_1$$\rightarrow$$3$P$_0$(492.2 GHz) neutral carbon lines, and weakly detects high-J CO transitions in the bridge. These results indicate that the bridge is composed of a warm multi-phase medium consistent with shock and turbulent heating. Despite low star formation rates in the bridge, the [C II] emission appears to be enhanced, reaching [C II]/FIR ratios of 3.3% in parts of the bridge. Both the [C II] and [O I] lines show broad intrinsic multi-component profiles, similar to those seen in previous CO 1-0 and H I observations. The [C II] emission shares similar line profiles with both the double-peaked H I profiles and shares a high-velocity component with single-peaked CO profiles in the bridge, suggesting that the [C II] emission originates in both the neutral and molecular phases. We show that it is feasible that a combination of turbulently heated H$_2$ and high column-density H I, resulting from the galaxy collision, is responsible for the enhanced [C II] emission.
△ Less
Submitted 30 January, 2018;
originally announced January 2018.
-
Interaction between a pulsating jet and a surrounding disk wind. A hydrodynamical perspective
Authors:
B. Tabone,
A. Raga,
S. Cabrit,
G. Pineau des Forêts
Abstract:
The molecular richness of fast protostellar jets within 20-100 au of their source, despite strong ultraviolet irradiation, remains a challenge for the models investigated so far. We aim to investigate the effect of interaction between a time-variable jet and a surrounding steady disk wind, to assess the possibility of jet chemical enrichement by the wind, and the characteristic signatures of such…
▽ More
The molecular richness of fast protostellar jets within 20-100 au of their source, despite strong ultraviolet irradiation, remains a challenge for the models investigated so far. We aim to investigate the effect of interaction between a time-variable jet and a surrounding steady disk wind, to assess the possibility of jet chemical enrichement by the wind, and the characteristic signatures of such a configuration. We have constructed an analytic model of a jet bow shock driven into a surrounding slower disk wind in the thin shell approximation. The refilling of the post bow shock cavity from below by the disk wind is also studied. An extension of the model to the case of two or more successive internal working surfaces (IWS) is made. We then compared this analytic model with numerical simulations with and without a surrounding disk wind. We find that at early times (of order the variability period), jet bow shocks travel in refilled pristine disk-wind material, before interacting with the cocoon of older bow shocks. This opens the possibility of bow shock chemical enrichment (if the disk wind is molecular and dusty) and of probing the unperturbed disk wind structure near the jet base. Several distinctive signatures of the presence of a surrounding disk wind are identified, in the bow shock morphology and kinematics. Numerical simulations validate our analytical approach and further show that at large scale, the passage of many jet IWS inside a disk wind produces a stationary V-shaped cavity, closing down onto the axis at a finite distance from the source.
△ Less
Submitted 10 January, 2018;
originally announced January 2018.
-
ALMA discovery of a rotating SO/SO$_2$ flow in HH212. A possible MHD disk wind?
Authors:
B. Tabone,
S. Cabrit,
E. Bianchi,
J. Ferreira,
G. Pineau des Forêts,
C. Codella,
A. Gusdorf,
F. Gueth,
L. Podio,
E. Chapillon
Abstract:
We wish to constrain the possible contribution of a magnetohydrodynamic disk wind (DW) to the HH212 molecular jet. We mapped the flow base with ALMA Cycle 4 at 0.13"$\sim$ 60 au resolution and compared these observations with synthetic DW predictions. We identified, in SO/SO$_2$, a rotating flow that is wider and slower than the axial SiO jet. The broad outflow cavity seen in C$^{34}$S is not carv…
▽ More
We wish to constrain the possible contribution of a magnetohydrodynamic disk wind (DW) to the HH212 molecular jet. We mapped the flow base with ALMA Cycle 4 at 0.13"$\sim$ 60 au resolution and compared these observations with synthetic DW predictions. We identified, in SO/SO$_2$, a rotating flow that is wider and slower than the axial SiO jet. The broad outflow cavity seen in C$^{34}$S is not carved by a fast wide-angle wind but by this slower agent. Rotation signatures may be fitted by a DW of a moderate lever arm launched out to $\sim$ 40 au with SiO tracing dust-free streamlines from 0.05-0.3 au. Such a DW could limit the core-to-star efficiency to $\leq$50%.
△ Less
Submitted 28 October, 2017; v1 submitted 3 October, 2017;
originally announced October 2017.
-
Nature of shocks revealed by SOFIA OI observations in the Cepheus E protostellar outflow
Authors:
A. Gusdorf,
S. Anderl,
B. Lefloch,
S. Leurini,
H. Wiesemeyer,
R. Guesten,
M. Benedettini,
C. Codella,
B. Godard,
A. I. Gomez-Ruiz,
K. Jacobs,
L. E. Kristensen,
P. Lesaffre,
G. Pineau des Forets,
D. C. Lis
Abstract:
Protostellar jets and outflows are key features of the star-formation process, and primary processes of the feedback of young stars on the interstellar medium. Understanding the underlying shocks is necessary to explain how jets and outflows are launched, and to quantify their chemical and energetic impacts on the surrounding medium. We performed a high-spectral resolution study of the [OI]…
▽ More
Protostellar jets and outflows are key features of the star-formation process, and primary processes of the feedback of young stars on the interstellar medium. Understanding the underlying shocks is necessary to explain how jets and outflows are launched, and to quantify their chemical and energetic impacts on the surrounding medium. We performed a high-spectral resolution study of the [OI]$_{\rm 63 μm}$ emission in the outflow of the intermediate-mass Class 0 protostar Cep E-mm. We present observations of the OI $^3$P$_1 \rightarrow$ $^3$P$_2$, OH between $^2Π_{1/2}$ $J = 3/2$ and $J = 1/2$ at 1837.8 GHz, and CO (16-15) lines with SOFIA-GREAT at three positions in the Cep E outflow: mm (the driving protostar), BI (in the southern lobe), and BII (the terminal position in the southern lobe). The CO line is detected at all three positions. The OI line is detected in BI and BII, whereas the OH line is not detected. In BII, we identify three kinematical components in OI and CO, already detected in CO: the jet, the HH377 terminal bow-shock, and the outflow cavity. The OI column density is higher in the outflow cavity than in the jet, which itself is higher than in the terminal shock. The terminal shock is where the abundance ratio of OI to CO is the lowest (about 0.2), whereas the jet component is atomic (ratio $\sim$2.7). In the jet, we compare the OI observations with shock models that successfully fit the integrated intensity of 10 CO lines: these models do not fit the OI data. The high intensity of OI emission points towards the propagation of additional dissociative or alternative FUV-irradiated shocks, where the illumination comes from the shock itself. From the sample of low-to-high mass protostellar outflows where similar observations have been performed, the effects of illumination seem to increase with the mass of the protostar.
△ Less
Submitted 12 April, 2017;
originally announced April 2017.
-
Powerful H2 Line Cooling in Stephan's Quintet II. Group-wide Gas and Shock Modeling of the Warm H2 and a Comparison with [CII]157.7um Emission and Kinematics
Authors:
P. N. Appleton,
P. Guillard,
A. Togi,
K. Alatalo,
F. Boulanger,
M. Cluver,
G. Pineau des Forets,
U. Lisenfeld,
P. Ogle,
C. K. Xu
Abstract:
We map for the first time the two-dimensional H_2 excitation of warm intergalactic gas in Stephan's Quintet on group-wide (50 x 35 kpc^2) scales to quantify the temperature, mass and warm-H_2 mass fraction as a function of position using Spitzer. Molecular gas temperatures are seen to rise (to T > 700K) and the slope of the power--law density-temperature relation flattens along the main ridge of t…
▽ More
We map for the first time the two-dimensional H_2 excitation of warm intergalactic gas in Stephan's Quintet on group-wide (50 x 35 kpc^2) scales to quantify the temperature, mass and warm-H_2 mass fraction as a function of position using Spitzer. Molecular gas temperatures are seen to rise (to T > 700K) and the slope of the power--law density-temperature relation flattens along the main ridge of the filament, defining the region of maximum heating. We also performed MHD modeling of the excitation properties of the warm gas, to map the velocity structure and energy deposition rate of slow and fast molecular shocks. Slow magnetic shocks were required to explain the power radiated from the lowest--lying rotational states of H_2, and strongly support the idea that energy cascades down to small scales and low velocities from the fast collision of NGC 7318b with group-wide gas. The highest levels of heating of the warm H_2 is strongly correlated with the large-scale stirring of the medium as measured by [CII] spectroscopy with Herschel. H_2 is also seen associated with a separate bridge which extends towards the Seyfert nucleus in NGC 7319, both from Spitzer and CARMA CO observations. This opens up the possibility that both galaxy collisions and AGN outflows can turbulently heat gas on large-scales in compact groups. The observations provide a laboratory for studying the effects of turbulent energy dissipation on group-wide scales that may provide clues about the heating and cooling of gas at high-z in early galaxy and protogalaxy formation.
△ Less
Submitted 17 January, 2017; v1 submitted 11 January, 2017;
originally announced January 2017.
-
Molecule survival in magnetized protostellar disk winds. II. Predicted H2O line profiles versus Herschel/HIFI observations
Authors:
W. Yvart,
S. Cabrit,
G. Pineau des Forets,
J. Ferreira
Abstract:
We investigate whether the broad wings of H2O emission identified with Herschel towards low-mass Class 0 and Class 1 protostars may be consistent with an origin in a dusty MHD disk wind, and the constraints it would set on the underlying disk properties. We present synthetic H2O line profiles predictions for a typical MHD disk wind solution with various values of disk accretion rate, stellar mass,…
▽ More
We investigate whether the broad wings of H2O emission identified with Herschel towards low-mass Class 0 and Class 1 protostars may be consistent with an origin in a dusty MHD disk wind, and the constraints it would set on the underlying disk properties. We present synthetic H2O line profiles predictions for a typical MHD disk wind solution with various values of disk accretion rate, stellar mass, extension of the launching area, and view angle. We compare them in terms of line shapes and intensities with the HIFI profiles observed by the WISH Key Program. We find that a dusty MHD disk wind launched from 0.2--0.6 AU AU to 3--25 AU can reproduce to a remarkable degree the observed shapes and intensities of the broad H2O component, both in the fundamental 557 GHz line and in more excited lines. Such a model also readily reproduces the observed correlation of 557 GHz line luminosity with envelope density, if the infall rate at 1000 AU is 1--3 times the disk accretion rate in the wind ejection region. It is also compatible with the typical disk size and bolometric luminosity in the observed targets. However, the narrower line profiles in Class 1 sources suggest that MHD disk winds in these sources, if present, would have to be slower and/or less water rich than in Class 0 sources. In conclusion, MHD disk winds appear as a valid (though not unique) option to consider for the origin of the broad H2O component in low-mass protostars. ALMA appears ideally suited to further test this model by searching for resolved signatures of the warm and slow wide-angle molecular wind that would be predicted.
△ Less
Submitted 12 November, 2015;
originally announced November 2015.
-
Warm molecular Hydrogen at high redshift with the James Webb Space Telescope
Authors:
P. Guillard,
F. Boulanger,
M. D. Lehnert,
P. N. Appleton,
G. Pineau des Forêts
Abstract:
The build-up of galaxies is regulated by a complex interplay between gravitational collapse, galaxy merging and feedback related to AGN and star formation. The energy released by these processes has to dissipate for gas to cool, condense, and form stars. How gas cools is thus a key to understand galaxy formation. \textit{Spitzer Space Telescope} infrared spectroscopy revealed a population of galax…
▽ More
The build-up of galaxies is regulated by a complex interplay between gravitational collapse, galaxy merging and feedback related to AGN and star formation. The energy released by these processes has to dissipate for gas to cool, condense, and form stars. How gas cools is thus a key to understand galaxy formation. \textit{Spitzer Space Telescope} infrared spectroscopy revealed a population of galaxies with weak star formation and unusually powerful H$_2$ line emission. This is a signature of turbulent dissipation, sustained by large-scale mechanical energy injection. The cooling of the multiphase interstellar medium is associated with emission in the H$_2$ lines. These results have profound consequences on our understanding of regulation of star formation, feedback and energetics of galaxy formation in general. The fact that H$_2$ lines can be strongly enhanced in high-redshift turbulent galaxies will be of great importance for the \textit{James Webb Space Telescope} observations which will unveil the role that H$_2$ plays as a cooling agent in the era of galaxy assembly.
△ Less
Submitted 7 October, 2015;
originally announced October 2015.
-
Challenging shock models with SOFIA OH observations in the high-mass star-forming region Cepheus A
Authors:
A. Gusdorf,
R. Guesten,
K. M. Menten,
D. R. Flower,
G. Pineau des Forets,
C. Codella,
T. Csengeri,
A. I. Gomez-Ruiz,
S. Heyminck,
K. Jacobs,
L. E. Kristensen,
S. Leurini,
M. A. Requena-Torres,
S. F. Wampfler,
H. Wiesemeyer,
F. Wyrowski
Abstract:
OH is a key molecule in H2O chemistry, a valuable tool for probing physical conditions, and an important contributor to the cooling of shock regions. OH participates in the re-distribution of energy from the protostar towards the surrounding ISM. Our aim is to assess the origin of the OH emission from the Cepheus A massive star-forming region and to constrain the physical conditions prevailing in…
▽ More
OH is a key molecule in H2O chemistry, a valuable tool for probing physical conditions, and an important contributor to the cooling of shock regions. OH participates in the re-distribution of energy from the protostar towards the surrounding ISM. Our aim is to assess the origin of the OH emission from the Cepheus A massive star-forming region and to constrain the physical conditions prevailing in the emitting gas. We thus want to probe the processes at work during the formation of massive stars. We present spectrally resolved observations of OH towards the outflows of Cepheus A with the GREAT spectrometer onboard the SOFIA telescope. Three triplets were observed at 1834.7 GHz, 1837.8 GHz, and 2514.3 GHz (163.4, 163.1, and 119.2 microns), at angular resolutions of 16.3", 16.3", and 11.9", respectively. We present the CO (16-15) spectrum at the same position. We compared the integrated intensities in the redshifted wings to shock models. The two triplets near 163 microns are detected in emission with blending hyperfine structure unresolved. Their profiles and that of CO can be fitted by a combination of 2 or 3 Gaussians. The observed 119.2 microns triplet is seen in absorption, since its blending hyperfine structure is unresolved, but with three line-of-sight components and a blueshifted emission wing consistent with that of the other lines. The OH line wings are similar to those of CO, suggesting that they emanate from the same shocked structure. Under this common origin assumption, the observations fall within the model predictions and within the range of use of our model only if we consider that four shock structures are caught in our beam. Our comparisons suggest that the observations might be consistently fitted by a J-type model with nH > 1e5 cm-3, v > 20 km/s, and with a filling factor of ~1. Such a high density is generally found in shocks associated to high-mass protostars.
△ Less
Submitted 28 September, 2015;
originally announced September 2015.
-
Sulphur-bearing molecules in diffuse molecular clouds: new results from SOFIA/GREAT and the IRAM 30 m telescope
Authors:
D. A. Neufeld,
B. Godard,
M. Gerin,
G. Pineau des Forêts,
C. Bernier,
E. Falgarone,
U. U. Graf,
R. Güsten,
E. Herbst,
P. Lesaffre,
P. Schilke,
P. Sonnentrucker,
H. Wiesemeyer
Abstract:
We have observed five sulphur-bearing molecules in foreground diffuse molecular clouds lying along the sight-lines to five bright continuum sources. We have used the GREAT instrument on SOFIA to observe the 1383 GHz $^2Π_{3/2} J=5/2-3/2$ transitions of SH towards the star-forming regions W31C, G29.96-0.02, G34.3+0.1, W49N and W51, detecting foreground absorption towards all five sources; and the E…
▽ More
We have observed five sulphur-bearing molecules in foreground diffuse molecular clouds lying along the sight-lines to five bright continuum sources. We have used the GREAT instrument on SOFIA to observe the 1383 GHz $^2Π_{3/2} J=5/2-3/2$ transitions of SH towards the star-forming regions W31C, G29.96-0.02, G34.3+0.1, W49N and W51, detecting foreground absorption towards all five sources; and the EMIR receivers on the IRAM 30m telescope at Pico Veleta to detect the H$_2$S 1(10)-1(01), CS J=2-1 and SO 3(2)-2(1) transitions. In nine foreground absorption components detected towards these sources, the inferred column densities of the four detected molecules showed relatively constant ratios, with N(SH)/N(H$_2$S) in the range 1.1 - 3.0, N(CS)/N(H$_2$S) in the range 0.32 - 0.61, and N(SO)/N(H$_2$S) in the range 0.08 - 0.30. The observed SH/H$_2$ ratios - in the range (0.5-2.6) $\times 10^{-8}$ - indicate that SH (and other sulphur-bearing molecules) account for << 1% of the gas-phase sulphur nuclei. The observed abundances of sulphur-bearing molecules, however, greatly exceed those predicted by standard models of cold diffuse molecular clouds, providing further evidence for the enhancement of endothermic reaction rates by elevated temperatures or ion-neutral drift. We have considered the observed abundance ratios in the context of shock and turbulent dissipation region (TDR) models. Using the TDR model, we find that the turbulent energy available at large scale in the diffuse ISM is sufficient to explain the observed column densities of SH and CS. Standard shock and TDR models, however, fail to reproduce the column densities of H$_2$S and SO by a factor of about 10; more elaborate shock models - in which account is taken of the velocity drift, relative to H$_2$, of SH molecules produced by the dissociative recombination of H$_3$S$^+$ - reduce this discrepancy to a factor ~ 3.
△ Less
Submitted 19 February, 2015;
originally announced February 2015.
-
Turbulent energy dissipation and intermittency in ambipolar diffusion magnetohydrodynamics
Authors:
Georgios Momferratos,
Pierre Lesaffre,
Edith Falgarone,
Guillaume Pineau des Forêts
Abstract:
The dissipation of kinetic and magnetic energy in the interstellar medium (ISM) can proceed through viscous, Ohmic or ambipolar diffusion (AD). It occurs at very small scales compared to the scales at which energy is presumed to be injected. This localized heating may impact the ISM evolution but also its chemistry, thus providing observable features. Here, we perform 3D spectral simulations of de…
▽ More
The dissipation of kinetic and magnetic energy in the interstellar medium (ISM) can proceed through viscous, Ohmic or ambipolar diffusion (AD). It occurs at very small scales compared to the scales at which energy is presumed to be injected. This localized heating may impact the ISM evolution but also its chemistry, thus providing observable features. Here, we perform 3D spectral simulations of decaying magnetohydrodynamic turbulence including the effects of AD. We find that the AD heating power spectrum peaks at scales in the inertial range, due to a strong alignment of the magnetic and current vectors in the dissipative range. AD affects much greater scales than the AD scale predicted by dimensional analysis. We find that energy dissipation is highly concentrated on thin sheets. Its probability density function follows a lognormal law with a power-law tail which hints at intermittency, a property which we quantify by use of structure function exponents. Finally, we extract structures of high dissipation, defined as connected sets of points where the total dissipation is most intense and we measure the scaling exponents of their geometric and dynamical characteristics: the inclusion of AD favours small sizes in the dissipative range.
△ Less
Submitted 12 February, 2015;
originally announced February 2015.
-
Collisional excitation of water by hydrogen atoms
Authors:
F. Daniel,
A. Faure,
P. J. Dagdigian,
M. -L. Dubernet,
F. lique,
G. Pineau des forêts
Abstract:
We present quantum dynamical calculations that describe the rotational excitation of H$_2$O due to collisions with H atoms. We used a recent, high accuracy potential energy surface, and solved the collisional dynamics with the close-coupling formalism, for total energies up to 12 000 cm$^{-1}$. From these calculations, we obtained collisional rate coefficients for the first 45 energy levels of bot…
▽ More
We present quantum dynamical calculations that describe the rotational excitation of H$_2$O due to collisions with H atoms. We used a recent, high accuracy potential energy surface, and solved the collisional dynamics with the close-coupling formalism, for total energies up to 12 000 cm$^{-1}$. From these calculations, we obtained collisional rate coefficients for the first 45 energy levels of both ortho- and para-H$_2$O and for temperatures in the range T = 5-1500 K. These rate coefficients are subsequently compared to the values previously published for the H$_2$O / He and H$_2$O / H$_2$ collisional systems. It is shown that no simple relation exists between the three systems and that specific calculations are thus mandatory.
△ Less
Submitted 29 October, 2014;
originally announced October 2014.
-
Exceptional AGN-driven turbulence inhibits star formation in the 3C 326N radio-galaxy
Authors:
P. Guillard,
F. Boulanger,
M. D. Lehnert,
G. Pineau des Forêts,
F. Combes,
E. Falgarone,
J. Bernard-Salas
Abstract:
We detect bright [CII]158$μ$m line emission from the radio galaxy 3C 326N at z=0.09, which shows weak star formation ($SFR<0.07$M$_{\odot}$~yr$^{-1}$) despite having strong H$_2$ line emission and $2\times 10^9$M$_{\odot}$ of molecular gas. The [CII] line is twice as strong as the 0-0S(1) 17$μ$m H$_2$ line, and both lines are much in excess what is expected from UV heating. We combine infrared Spi…
▽ More
We detect bright [CII]158$μ$m line emission from the radio galaxy 3C 326N at z=0.09, which shows weak star formation ($SFR<0.07$M$_{\odot}$~yr$^{-1}$) despite having strong H$_2$ line emission and $2\times 10^9$M$_{\odot}$ of molecular gas. The [CII] line is twice as strong as the 0-0S(1) 17$μ$m H$_2$ line, and both lines are much in excess what is expected from UV heating. We combine infrared Spitzer and Herschel data with gas and dust modeling to infer the gas physical conditions. The [CII] line traces 30 to 50% of the molecular gas mass, which is warm (70<T<100K) and at moderate densities $700<n_{H}<3000$cm$^{-3}$. The [CII] line is broad with a blue-shifted wing, and likely to be shaped by a combination of rotation, outflowing gas, and turbulence. It matches the near-infrared H$_2$ and the Na D optical absorption lines. If the wing is interpreted as an outflow, the mass loss rate would be larger than 20M$_{\odot}$/yr, and the depletion timescale shorter than the orbital timescale ($10^8$yr). These outflow rates may be over-estimated because the stochastic injection of turbulence on galactic scales can contribute to the skewness of the line profile and mimic outflowing gas. We argue that the dissipation of turbulence is the main heating process of this gas. Cosmic rays can also contribute to the heating but they require an average gas density larger than the observational constraints. We show that strong turbulent support maintains a high gas vertical scale height (0.3-4kpc) in the disk and can inhibit the formation of gravitationally-bound structures at all scales, offering a natural explanation for the weakness of star formation in 3C 326N. To conclude, the bright [CII] line indicates that strong AGN jet-driven turbulence may play a key role in enhancing the amount of molecular gas (positive feedback) but yet can prevent star formation on galactic scales (negative feedback).
△ Less
Submitted 7 November, 2014; v1 submitted 22 October, 2014;
originally announced October 2014.
-
Chemical probes of turbulence in the diffuse medium: the TDR model
Authors:
Benjamin Godard,
Edith Falgarone,
Guillaume Pineau Des Forêts
Abstract:
Context. Tens of light hydrides and small molecules have now been detected over several hundreds sight lines sampling the diffuse interstellar medium (ISM) in both the Solar neighbourhood and the inner Galactic disk. They provide unprecedented statistics on the first steps of chemistry in the diffuse gas. Aims. These new data confirm the limitations of the traditional chemical pathways driven by t…
▽ More
Context. Tens of light hydrides and small molecules have now been detected over several hundreds sight lines sampling the diffuse interstellar medium (ISM) in both the Solar neighbourhood and the inner Galactic disk. They provide unprecedented statistics on the first steps of chemistry in the diffuse gas. Aims. These new data confirm the limitations of the traditional chemical pathways driven by the UV photons and the cosmic rays (CR) and the need for additional energy sources, such as turbulent dissipation, to open highly endoenergetic formation routes. The goal of the present paper is to further investigate the link between specific species and the properties of the turbulent cascade in particular its space-time intermittency. Methods. We have analysed ten different atomic and molecular species in the framework of the updated model of turbulent dissipation regions (TDR). We study the influence on the abundances of these species of parameters specific to chemistry (density, UV field, and CR ionisation rate) and those linked to turbulence (the average turbulent dissipation rate, the dissipation timescale, and the ion neutral velocity drift in the regions of dissipation). Results. The most sensitive tracers of turbulent dissipation are the abundances of CH+ and SH+, and the column densities of the J = 3, 4, 5 rotational levels of H2 . The abundances of CO, HCO+, and the intensity of the 158 $μ$m [CII] emission line are significantly enhanced by turbulent dissipation. The vast diversity of chemical pathways allows the independent determinations of free parameters never estimated before: an upper limit to the average turbulent dissipation rate, $\overline{\varepsilon}$ < 10$^{-23}$ erg cm$^{-3}$ s$^{-1}$ for $n_H$=20 cm$^{-3}$, from the CH+ abundance; an upper limit to the ion-neutral velocity drift, $u_{in}$ < 3.5 km s$^{-1}$, from the SH+ to CH+ abundance ratio; and a range of dissipation timescales, 100 < $τ_V$ < 1000 yr, from the CO to HCO+ abundance ratio. For the first time, we reproduce the large abundances of CO observed on diffuse lines of sight, and we show that CO may be abundant even in regions with UV-shieldings as low as $5 \times 10^{-3}$ mag. The best range of parameters also reproduces the abundance ratios of OH, C2H, and H2O to HCO+ and are consistent with the known properties of the turbulent cascade in the Galactic diffuse ISM. Conclusions. Our results disclose an unexpected link between the dissipation of turbulence and the emergence of molecular richness in the diffuse ISM. Some species, such as CH+ or SH+, turn out to be unique tracers of the energy trail in the ISM. In spite of some degeneracy, the properties of the turbulent cascade, down to dissipation, can be captured through specific molecular abundances.
△ Less
Submitted 16 August, 2014;
originally announced August 2014.
-
Shocks in dense clouds. IV. Effects of grain-grain processing on molecular line emission
Authors:
S. Anderl,
V. Guillet,
G. Pineau des Forêts,
D. R. Flower
Abstract:
Grain-grain processing has been shown to be an indispensable ingredient of shock modelling in high density environments. For densities higher than \sim10^5 cm-3, shattering becomes a self-enhanced process that imposes severe chemical and dynamical consequences on the shock characteristics. Shattering is accompanied by the vaporization of grains, which can directly release SiO to the gas phase. Giv…
▽ More
Grain-grain processing has been shown to be an indispensable ingredient of shock modelling in high density environments. For densities higher than \sim10^5 cm-3, shattering becomes a self-enhanced process that imposes severe chemical and dynamical consequences on the shock characteristics. Shattering is accompanied by the vaporization of grains, which can directly release SiO to the gas phase. Given that SiO rotational line radiation is used as a major tracer of shocks in dense clouds, it is crucial to understand the influence of vaporization on SiO line emission. We have developed a recipe for implementing the effects of shattering and vaporization into a 2-fluid shock model, resulting in a reduction of computation time by a factor \sim100 compared to a multi-fluid modelling approach. This implementation was combined with an LVG-based modelling of molecular line radiation transport. Using this model we calculated grids of shock models to explore the consequences of different dust-processing scenarios. Grain-grain processing is shown to have a strong influence on C-type shocks for a broad range of magnetic fields: they become hotter and thinner. The reduction in column density of shocked gas lowers the intensity of molecular lines, at the same time as higher peak temperatures increase the intensity of highly excited transitions compared to shocks without grain-grain processing. For OH the net effect is an increase in line intensities, while for CO and H2O it is the contrary. The intensity of H2 emission is decreased in low transitions and increased for highly excited lines. For all molecules, the highly excited lines become sensitive to the value of the magnetic field. Although vaporization increases the intensity of SiO rotational lines, this effect is weakened by the reduced shock width. The release of SiO early in the hot shock changes the excitation characteristics of SiO radiation.
△ Less
Submitted 1 August, 2014;
originally announced August 2014.
-
Dense molecular globulettes and the dust arc towards the runaway O star AE Aur (HD 34078)
Authors:
P. Gratier,
J. Pety,
P. Boissé,
S. Cabrit,
P. Lesaffre,
M. Gerin,
G. Pineau des Forêts
Abstract:
Some runaway stars are known to display IR arc-like structures around them, resulting from their interaction with surrounding interstellar material. The properties of these features as well as the processes involved in their formation are still poorly understood. We aim at understanding the physical mechanisms that shapes the dust arc observed near the runaway O star AEAur (HD34078). We obtained a…
▽ More
Some runaway stars are known to display IR arc-like structures around them, resulting from their interaction with surrounding interstellar material. The properties of these features as well as the processes involved in their formation are still poorly understood. We aim at understanding the physical mechanisms that shapes the dust arc observed near the runaway O star AEAur (HD34078). We obtained and analyzed a high spatial resolution map of the CO(1-0) emission that is centered on HD34078, and that combines data from both the IRAM interferometer and 30m single-dish antenna. The line of sight towards HD34078 intersects the outer part of one of the detected globulettes, which accounts for both the properties of diffuse UV light observed in the field and the numerous molecular absorption lines detected in HD34078's spectra, including those from highly excited H2 . Their modeled distance from the star is compatible with the fact that they lie on the 3D paraboloid which fits the arc detected in the 24 μm Spitzer image. Four other compact CO globulettes are detected in the mapped area. These globulettes have a high density and linewidth, and are strongly pressure-confined or transient. The good spatial correlation between the CO globulettes and the IR arc suggests that they result from the interaction of the radiation and wind emitted by HD 34078 with the ambient gas. However, the details of this interaction remain unclear. A wind mass loss rate significantly larger than the value inferred from UV lines is favored by the large IR arc size, but does not easily explain the low velocity of the CO globulettes. The effect of radiation pressure on dust grains also meets several issues in explaining the observations. Further observational and theoretical work is needed to fully elucidate the processes shaping the gas and dust in bow shocks around runaway O stars. (Abridged)
△ Less
Submitted 26 March, 2014;
originally announced March 2014.
-
Interstellar chemistry of nitrogen hydrides in dark clouds
Authors:
Romane Le Gal,
Pierre Hily-Blant,
Alexandre Faure,
Guillaume Pineau des Forêts,
Claire Rist,
Sébastien Maret
Abstract:
The aim of the present work is to perform a comprehensive analysis of the interstellar chemistry of nitrogen, focussing on the gas-phase formation of the smallest polyatomic species and in particular nitrogen hydrides. We present a new chemical network in which the kinetic rates of critical reactions have been updated based on recent experimental and theoretical studies, including nuclear spin bra…
▽ More
The aim of the present work is to perform a comprehensive analysis of the interstellar chemistry of nitrogen, focussing on the gas-phase formation of the smallest polyatomic species and in particular nitrogen hydrides. We present a new chemical network in which the kinetic rates of critical reactions have been updated based on recent experimental and theoretical studies, including nuclear spin branching ratios. Our network thus treats the different spin symmetries of the nitrogen hydrides self-consistently together with the ortho and para forms of molecular hydrogen. This new network is used to model the time evolution of the chemical abundances in dark cloud conditions. The steady-state results are analysed, with special emphasis on the influence of the overall amounts of carbon, oxygen, and sulphur. Our calculations are also compared with Herschel/HIFI observations of NH, NH$_2$, and NH$_3$ detected towards the external envelope of the protostar IRAS 16293-2422. The observed abundances and abundance ratios are reproduced for a C/O gas-phase elemental abundance ratio of $\sim0.8$, provided that the sulphur abundance is depleted by a factor larger than 2. The ortho-to-para ratio of H$_2$ in these models is $\sim10^{-3}$. Our models also provide predictions for the ortho-to-para ratios of NH$_2$ and NH$_3$ of $\sim2.3$ and $\sim0.7$ respectively. We conclude that the abundances of nitrogen hydrides in dark cloud conditions are consistent with the gas-phase synthesis predicted with our new chemical network.
△ Less
Submitted 3 December, 2013; v1 submitted 20 November, 2013;
originally announced November 2013.
-
Shock-Enhanced C+ Emission and the Detection of H2O from Stephan's Quintet's Group-Wide Shock using Herschel
Authors:
P. N. Appleton,
P. Guillard,
F. Boulanger,
M. E. Cluver,
P. Ogle,
E. Falgarone,
G. Pineau Des Forets,
E. O'Sullivan,
P. -A. Duc,
S. Gallagher,
Y. Gao,
T. Jarrett,
I. Konstantopoulos,
U. Lisenfeld,
S. Lord,
N. Lu,
B. W. Peterson,
C. Struck,
E. Sturm,
R. Tuffs,
I. Valchanov,
P. van der Werf,
K. C. Xu
Abstract:
We present the first Herschel spectroscopic detections of the [OI]63 and [CII]158 micron fine-structure transitions, and a single para-H2O line from the 35 x 15 kpc^2 shocked intergalactic filament in Stephan's Quintet. The filament is believed to have been formed when a high-speed intruder to the group collided with clumpy intergroup gas. Observations with the PACS spectrometer provide evidence f…
▽ More
We present the first Herschel spectroscopic detections of the [OI]63 and [CII]158 micron fine-structure transitions, and a single para-H2O line from the 35 x 15 kpc^2 shocked intergalactic filament in Stephan's Quintet. The filament is believed to have been formed when a high-speed intruder to the group collided with clumpy intergroup gas. Observations with the PACS spectrometer provide evidence for broad (> 1000 km s^-1) luminous [CII] line profiles, as well as fainter [OI]63micron emission. SPIRE FTS observations reveal water emission from the p-H2O (111-000) transition at several positions in the filament, but no other molecular lines. The H2O line is narrow, and may be associated with denser intermediate-velocity gas experiencing the strongest shock-heating. The [CII]/PAH{tot) and [CII]/FIR ratios are too large to be explained by normal photo-electric heating in PDRs. HII region excitation or X-ray/Cosmic Ray heating can also be ruled out. The observations lead to the conclusion that a large fraction the molecular gas is diffuse and warm. We propose that the [CII], [OI] and warm H2 line emission is powered by a turbulent cascade in which kinetic energy from the galaxy collision with the IGM is dissipated to small scales and low-velocities, via shocks and turbulent eddies. Low-velocity magnetic shocks can help explain both the [CII]/[OI] ratio, and the relatively high [CII]/H2 ratios observed. The discovery that [CII] emission can be enhanced, in large-scale turbulent regions in collisional environments has implications for the interpretation of [CII] emission in high-z galaxies.
△ Less
Submitted 5 September, 2013;
originally announced September 2013.
-
The CN/C15 N isotopic ratio towards dark clouds
Authors:
Pierre Hily-Blant,
G. Pineau Des Forêts,
Alexandre Faure,
Romane Le Gal,
Marco Padovani
Abstract:
Understanding the origin of the composition of solar system cosmomaterials is a central question, not only in the cosmochemistry and astrochemistry fields, and requires various approaches to be combined. Measurements of isotopic ratios in cometary materials provide strong constraints on the content of the protosolar nebula. Their relation with the composition of the parental dark clouds is, howeve…
▽ More
Understanding the origin of the composition of solar system cosmomaterials is a central question, not only in the cosmochemistry and astrochemistry fields, and requires various approaches to be combined. Measurements of isotopic ratios in cometary materials provide strong constraints on the content of the protosolar nebula. Their relation with the composition of the parental dark clouds is, however, still very elusive. In this paper, we bring new constraints based on the isotopic composition of nitrogen in dark clouds, with the aim of understanding the chemical processes that are responsible for the observed isotopic ratios. We have observed and detected the fundamental rotational transition of C$^{15}$N towards two starless dark clouds, L1544 and L1498. We were able to derive the column density ratio of C$^{15}$N over $^{13}$CN towards the same clouds, and obtain the CN/C$^{15}$N isotopic ratios, which were found to be $500\pm75$ for both L1544 and L1498. These values are therefore marginally consistent with the protosolar value of 441. Moreover, this ratio is larger than the isotopic ratio of nitrogen measured in HCN. In addition, we present model calculations of the chemical fractionation of nitrogen in dark clouds, which make it possible to understand how CN can be deprived of $^{15}$N and HCN can simultaneously be enriched in heavy nitrogen. The non-fractionation of N2H+, however, remains an open issue and we propose some chemical way of alleviating the discrepancy between model predictions and the observed ratios.
△ Less
Submitted 20 August, 2013;
originally announced August 2013.
-
Low velocity shocks: signatures of turbulent dissipation in diffuse irradiated gas
Authors:
Pierre Lesaffre,
Guillaume Pineau des Forêts,
Benjamin Godard,
Pierre Guillard,
François Boulanger,
Edith Falgarone
Abstract:
We examine the chemical and emission properties of mildly irradiated (G0=1) magnetised shocks in diffuse media (nH=10^2 to 10^4 /cm3) at low to moderate velocities (from 3 to 40 km/s). Results: The formation of some molecules relies on endoergic reactions. In J-shocks, their abundances are enhanced by several orders of magnitude for shock velocities as low as 7 km/s. Otherwise most chemical proper…
▽ More
We examine the chemical and emission properties of mildly irradiated (G0=1) magnetised shocks in diffuse media (nH=10^2 to 10^4 /cm3) at low to moderate velocities (from 3 to 40 km/s). Results: The formation of some molecules relies on endoergic reactions. In J-shocks, their abundances are enhanced by several orders of magnitude for shock velocities as low as 7 km/s. Otherwise most chemical properties of J-type shocks vary over less than an order of magnitude between velocities from about 7 to about 30 km/s, where H2 dissociation sets in. C-type shocks display a more gradual molecular enhancement as the shock velocity increases. We quantify the energy flux budget (fluxes of kinetic, radiated and magnetic energies) with emphasis on the main cooling lines of the cold interstellar medium. Their sensitivity to shock velocity is such that it allows observations to constrain statistical distributions of shock velocities. We fit various probability distribution functions (PDFs) of shock velocities to spectroscopic observations of the galaxy-wide shock in Stephan's Quintet (SQ) and of a Galactic line of sight sampling diffuse molecular gas in Chamaeleon. In both cases, low velocities bear the greatest statistical weight and the PDF is consistent with a bimodal distribution. In the very low velocity shocks (below 5 km/s), dissipation is due to ion-neutral friction which powers H2 low energy transitions and atomic lines. In moderate velocity shocks (20 km/s and above), the dissipation is due to viscous heating and accounts for most of the molecular emission. In our interpretation a significant fraction of the gas on the line of sight is shocked (from 4% to 66%). For example, C+ emission may trace shocks in UV irradiated gas where C+ is the dominant carbon species.
△ Less
Submitted 31 January, 2013;
originally announced January 2013.
-
Discovery of interstellar mercapto radicals (SH) with the GREAT instrument on SOFIA
Authors:
D. A. Neufeld,
E. Falgarone,
M. Gerin,
B. Godard,
E. Herbst,
G. Pineau des Forêts,
A. I. Vasyunin,
R. Güsten,
H. Wiesemeyer,
O. Ricken
Abstract:
We report the first detection of interstellar mercapto radicals, obtained along the sight-line to the submillimeter continuum source W49N. We have used the GREAT instrument on SOFIA to observe the 1383 GHz Doublet Pi 3/2 J = 5/2 - 3/2 lambda doublet in the upper sideband of the L1 receiver. The resultant spectrum reveals SH absorption in material local to W49N, as well as in foreground gas, unasso…
▽ More
We report the first detection of interstellar mercapto radicals, obtained along the sight-line to the submillimeter continuum source W49N. We have used the GREAT instrument on SOFIA to observe the 1383 GHz Doublet Pi 3/2 J = 5/2 - 3/2 lambda doublet in the upper sideband of the L1 receiver. The resultant spectrum reveals SH absorption in material local to W49N, as well as in foreground gas, unassociated with W49N, that is located along the sight-line. For the foreground material at velocities in the range 37 - 44 km/s with respect to the local standard of rest, we infer a total SH column density ~ 2.6 E+12 cm-2, corresponding to an abundance of ~ 7 E-9 relative to H2, and yielding an SH/H2S abundance ratio ~ 0.13. The observed SH/H2S abundance ratio is much smaller than that predicted by standard models for the production of SH and H2S in turbulent dissipation regions and shocks, and suggests that the endothermic neutral-neutral reaction SH + H2 -> H2S + H must be enhanced along with the ion-neutral reactions believed to produce CH+ and SH+ in diffuse molecular clouds.
△ Less
Submitted 14 February, 2012;
originally announced February 2012.
-
Turbulent molecular gas and star formation in the shocked intergalactic medium of Stephan's Quintet
Authors:
Pierre Guillard,
François Boulanger,
Guillaume Pineau des Forêts,
Edith Falgarone,
Antoine Gusdorf,
Michelle Cluver,
Philip Appleton,
Ute Lisenfeld,
Pierre-Alain Duc,
Patrick Ogle,
Kevin Xu
Abstract:
We report on single-dish radio CO observations towards the inter-galactic medium (IGM) of the Stephan's Quintet (SQ) group of galaxies. Extremely bright mid-IR H2 rotational line emission from warm molecular gas has been detected by Spitzer in the kpc-scale shock created by a galaxy collision. We detect in the IGM CO(1-0), (2-1) and (3-2) line emission with complex profiles, spanning a velocity ra…
▽ More
We report on single-dish radio CO observations towards the inter-galactic medium (IGM) of the Stephan's Quintet (SQ) group of galaxies. Extremely bright mid-IR H2 rotational line emission from warm molecular gas has been detected by Spitzer in the kpc-scale shock created by a galaxy collision. We detect in the IGM CO(1-0), (2-1) and (3-2) line emission with complex profiles, spanning a velocity range of 1000 km/s. The spectra exhibit the pre-shock recession velocities of the two colliding gas systems (5700 and 6700 km/s), but also intermediate velocities. This shows that much of the molecular gas has formed out of diffuse gas accelerated by the galaxy-tidal arm collision. A total H2 mass of 5x10^9 Msun is detected in the shock. The molecular gas carries a large fraction of the gas kinetic energy involved in the collision, meaning that this energy has not been thermalized yet. The turbulent kinetic energy of the H2 gas is at least a factor of 5 greater than the thermal energy of the hot plasma heated by the collision. The ratio between the warm H2 mass derived from Spitzer IRS spectroscopy and the H2 mass derived from CO fluxes is ~0.3 in the IGM of SQ, which is 10-100 times higher than in star-forming galaxies. In the shocked region, the ratio of the PAH-to-CO surface luminosities, commonly used to measure the star formation efficiency of the H2 gas, is lower (up to a factor 75) than the observed values in star-forming galaxies. We suggest that turbulence fed by the galaxy-tidal arm collision maintains a high heating rate within the H2 gas. This interpretation implies that the velocity dispersion on the scale of giant molecular clouds in SQ is one order of magnitude larger than the Galactic value. The high amplitude of turbulence may explain why this gas is not forming stars efficiently. [abridged version]
△ Less
Submitted 30 March, 2012; v1 submitted 13 February, 2012;
originally announced February 2012.
-
Molecule survival in magnetized protostellar disk winds. I. Chemical model and first results
Authors:
D. Panoglou,
S. Cabrit,
G. Pineau des Forets,
P. J. V. Garcia,
J. Ferreira,
F. Casse
Abstract:
Molecular counterparts to atomic jets have been detected within 1000 AU of young stars. Reproducing them is a challenge for proposed ejection models. We explore whether molecules may survive in an MHD disk wind invoked to reproduce the kinematics and tentative rotation signatures of atomic jets in T Tauri stars. The coupled ionization, chemical and thermal evolution along dusty flow streamlines is…
▽ More
Molecular counterparts to atomic jets have been detected within 1000 AU of young stars. Reproducing them is a challenge for proposed ejection models. We explore whether molecules may survive in an MHD disk wind invoked to reproduce the kinematics and tentative rotation signatures of atomic jets in T Tauri stars. The coupled ionization, chemical and thermal evolution along dusty flow streamlines is computed for a prescribed MHD disk wind solution, using a method developed for magnetized shocks in the interstellar medium. Irradiation by wind-attenuated coronal X-rays and FUV photons from accretion hot spots is included, with self-shielding of H2 and CO. Disk accretion rates of 5e-6, 1e-6 and 1e-7 solar masses per year are considered, representative of low-mass young protostars (Class 0), evolved protostars (Class I) and very active T Tauri stars (Class II). The disk wind has an onion-like thermo-chemical structure, with streamlines launched from larger radii having lower temperature and ionisation, and higher H2 abundance. The coupling between charged and neutral fluids is sufficient to eject molecules from the disk out to 9 AU. The launch radius beyond which most H2 survives moves outward with evolutionary stage. CO survives in the Class 0 but is significantly photodissociated in the Class I/II. Balance between ambipolar heating and molecular cooling establishes an asymptotic temperature 700-3000 K, with cooler jets at earlier protostellar stages. Endothermic formation of H2O is efficient with abundances up to 1e-4, while CH+ and SH+ can exceed 1e-6 in the Class I/II winds. A centrifugal MHD disk wind launched from beyond 0.2-1 AU can produce molecular jets/winds up to speeds 100 km/s in young low-mass stars. The model predicts a high ratio H2/CO and an increase of molecular launch radius, temperature, and flow width as the source evolves, in agreement with current observed trends.
△ Less
Submitted 14 December, 2011;
originally announced December 2011.
-
Nitrogen hydrides and the H2 ortho-to-para ratio in dark clouds
Authors:
Vincent Dislaire,
Pierre Hily-Blant,
Alexandre Faure,
Sébastien Maret,
Aurore Bacmann,
Guillaume Pineau Des Forêts
Abstract:
Nitrogen bearing species are common tracers of the physical conditions in a wide variety of objects, and most remarkably in dark clouds. The reservoir of gaseous nitrogen is expected to be atomic or molecular, but none of the two species are observable in the dark gas. Their abundances therefore derive indirectly from those of N-bearing species through chemical modelling. The recent years have acc…
▽ More
Nitrogen bearing species are common tracers of the physical conditions in a wide variety of objects, and most remarkably in dark clouds. The reservoir of gaseous nitrogen is expected to be atomic or molecular, but none of the two species are observable in the dark gas. Their abundances therefore derive indirectly from those of N-bearing species through chemical modelling. The recent years have accumulated data which stress our incomplete understanding of the nitrogen chemistry in dark cloud conditions. To tackle this problem of the nitrogen chemistry in cold gas, we have revised the formation of nitrogen hydrides, which is initiated by the key reaction \ce{N+ + H2 -> NH+ + H}. We propose a new rate for this reaction which depends on the ortho-to-para ratio of H2. This new rate allows to reproduce the abundance ratios of the three nitrogen hydrides, NH, \ce{NH2}, and \ce{NH3}, observed towards IRAS16293-2422, provided that the channel leading to NH from the dissociative recombination of \ce{N2H+} is not closed at low temperature. The ortho-to-para ratio of H2 is constrained to O/P=$10^{-3}$ by the abundance ratio NH:NH2, which provides a new method to measure O/P. This work stresses the need for reaction rates at the low temperatures of dark clouds, and for branching ratios of critical dissociative recombination reactions.
△ Less
Submitted 7 November, 2011;
originally announced November 2011.
-
Hydrogen cyanide and isocyanide in prestellar cores
Authors:
M. Padovani,
C. M. Walmsley,
M. Tafalla,
P. Hily-Blant,
G. Pineau des Forêts
Abstract:
We studied the abundance of HCN, H13CN, and HN13C in a sample of prestellar cores, in order to search for species associated with high density gas. We used the IRAM 30m radiotelescope to observe along the major and the minor axes of L1498, L1521E, and TMC 2, three cores chosen on the basis of their CO depletion properties. We mapped the J=1-0 transition of HCN, H13CN, and HN13C towards the source…
▽ More
We studied the abundance of HCN, H13CN, and HN13C in a sample of prestellar cores, in order to search for species associated with high density gas. We used the IRAM 30m radiotelescope to observe along the major and the minor axes of L1498, L1521E, and TMC 2, three cores chosen on the basis of their CO depletion properties. We mapped the J=1-0 transition of HCN, H13CN, and HN13C towards the source sample plus the J=1-0 transition of N2H+ and the J=2-1 transition of C18O in TMC 2. We used two different radiative transfer codes, making use of recent collisional rate calculations, in order to determine more accurately the excitation temperature, leading to a more exact evaluation of the column densities and abundances. We find that the optical depths of both H13CN(1-0) and HN13C(1-0) are non-negligible, allowing us to estimate excitation temperatures for these transitions in many positions in the three sources. The observed excitation temperatures are consistent with recent computations of the collisional rates for these species and they correlate with hydrogen column density inferred from dust emission. We conclude that HCN and HNC are relatively abundant in the high density zone, n(H2) about 10^5 cm-3, where CO is depleted. The relative abundance [HNC]/[HCN] differs from unity by at most 30 per cent consistent with chemical expectations. The three hyperfine satellites of HCN(1-0) are optically thick in the regions mapped, but the profiles become increasingly skewed to the blue (L1498 and TMC 2) or red (L1521E) with increasing optical depth suggesting absorption by foreground layers.
△ Less
Submitted 25 August, 2011;
originally announced August 2011.
-
Excitation of H$_2$ in photodissociation regions as seen by Spitzer
Authors:
Emilie Habart,
Alain Abergel,
Francois Boulanger,
Christine Joblin,
Laurent Verstraete,
Mathieu Compiègne,
Guillaume Pineau des Forêts,
Jacques Le Bourlot
Abstract:
We present spectroscopic observations obtained with the infrared Spitzer Space Telescope, which provide insight into the H$_2$ physics and gas energetics in photodissociation Regions (PDRs) of low to moderate far-ultraviolet (FUV) fields and densities. We analyze data on six well known Galactic PDRs (L1721, California, N7023E, Horsehead, rho Oph, N2023N), sampling a poorly explored range of excita…
▽ More
We present spectroscopic observations obtained with the infrared Spitzer Space Telescope, which provide insight into the H$_2$ physics and gas energetics in photodissociation Regions (PDRs) of low to moderate far-ultraviolet (FUV) fields and densities. We analyze data on six well known Galactic PDRs (L1721, California, N7023E, Horsehead, rho Oph, N2023N), sampling a poorly explored range of excitation conditions ($χ\sim 5-10^3$), relevant to the bulk of molecular clouds in galaxies. Spitzer observations of H$_2$ rotational lines are complemented with H$_2$ data, including ro-vibrational line measurements, obtained with ground-based telescopes and ISO, to constrain the relative contributions of ultraviolet pumping and collisions to the H$_2$ excitation. The data analysis is supported by model calculations with the Meudon PDR code. The observed column densities of rotationally excited H$_2$ are observed to be much higher than PDR model predictions. In the lowest excitation PDRs, the discrepancy between the model and the data is about one order of magnitude for rotational levels $J \ge $3. We discuss whether an enhancement in the H$_2$ formation rate or a local increase in photoelectric heating, as proposed for brighter PDRs in former ISO studies, may improve the data-model comparison. We find that an enhancement in the H$_2$ formation rates reduces the discrepancy, but the models still fall short of the data. This large disagreement suggests that our understanding of the formation and excitation of H$_2$ and/or of PDRs energetics is still incomplete. We discuss several explanations, which could be further tested using the Herschel Space Telescope
△ Less
Submitted 23 December, 2010;
originally announced December 2010.
-
Spitzer spectral line mapping of the HH211 outflow
Authors:
O. Dionatos,
B. Nisini,
S. Cabrit,
L. Kristensen,
G. Pineau des Forets
Abstract:
Aims: We employ archival Spitzer slit-scan observations of the HH211 outflow in order to investigate its warm gas content, assess the jet mass flux in the form of H2 and probe for the existence of an embedded atomic jet. Methods: Detected molecular and atomic lines are interpreted by means of emission line diagnostics and an existing grid of molecular shock models. The physical properties of the w…
▽ More
Aims: We employ archival Spitzer slit-scan observations of the HH211 outflow in order to investigate its warm gas content, assess the jet mass flux in the form of H2 and probe for the existence of an embedded atomic jet. Methods: Detected molecular and atomic lines are interpreted by means of emission line diagnostics and an existing grid of molecular shock models. The physical properties of the warm gas are compared against other molecular jet tracers and to the results of a similar study towards the L1448-C outflow. Results: We have detected and mapped the v=0-0 S(0) - S(7) H2 lines and fine-structure lines of S, Fe+, and Si+. H2 is detected down to 5" from the source and is characterized by a "cool" T~300K and a "warm" T~1000 K component, with an extinction Av ~ 8 mag. The amount of cool H2 towards the jet agrees with that estimated from CO assuming fully molecular gas. The warm component is well fitted by C-type shocks with a low beam filling factor ~ 0.01-0.04 and a mass-flux similar to the cool H2. The fine-structure line emission arises from dense gas with ionization fraction ~0.5 - 5 x 10e-3, suggestive of dissociative shocks. Line ratios to sulfur indicate that iron and silicon are depleted compared to solar abundances by a factor ~10-50. Conclusions: Spitzer spectral mapping observations reveal for the first time a cool H$_2$ component towards the CO jet of HH211 consistent with the CO material being fully molecular and warm at ~ 300 K. The maps also reveal for the first time the existence of an embedded atomic jet in the HH211 outflow that can be traced down to the central source position. Its significant iron and silicon depletion excludes an origin from within the dust sublimation zone around the protostar. The momentum-flux seems insufficient to entrain the CO jet, although current uncertainties on jet speed and shock conditions are too large for a definite conclusion.
△ Less
Submitted 4 June, 2010;
originally announced June 2010.
-
Observations and modeling of the dust emission from the H2-bright galaxy-wide shock in Stephan's Quintet
Authors:
P. Guillard,
F. Boulanger,
M. E. Cluver,
P. N. Appleton,
G. Pineau des Forets,
P. Ogle
Abstract:
Spitzer Space Telescope observations revealed powerful mid-infrared (mid-IR) H2 rotational line emission from the Stephan's Quintet (SQ) X-ray emitting large scale shock associated with a collision between two galaxies. Because H2 forms on dust grains, the presence of H2 is physically linked to the survival of dust, and we expect some dust emission to come from the molecular gas. To test this inte…
▽ More
Spitzer Space Telescope observations revealed powerful mid-infrared (mid-IR) H2 rotational line emission from the Stephan's Quintet (SQ) X-ray emitting large scale shock associated with a collision between two galaxies. Because H2 forms on dust grains, the presence of H2 is physically linked to the survival of dust, and we expect some dust emission to come from the molecular gas. To test this interpretation, IR observations and dust modeling are used to identify and characterize the thermal dust emission from the shocked molecular gas. The spatial distribution of the IR emission allows us to isolate the faint PAH and dust continuum emission associated with the molecular gas in the SQ shock. We model the spectral energy distribution (SED) of this emission, and fit it to Spitzer observations. Faint PAH and dust continuum emission are detected in the SQ shock, outside star-forming regions. The 12/24um flux ratio in the shock is remarkably close to that of the diffuse Galactic interstellar medium, leading to a Galactic PAH/VSG abundance ratio. However, the properties of the PAH emission spectrum in the shock differ from that of the Galaxy, which may suggest an enhanced fraction of large and neutrals PAHs. The IR SED is consistent with the expected emission from dust associated with the warm (>150K) H2 gas, heated by a UV radiation field of intensity comparable to that of the solar neighborhood, in agreement with GALEX UV observations. The presence of PAHs and dust grains in the high-speed (1000km/s) galaxy collision suggests that dust survives. We propose that the dust that survived destruction was in pre-shock gas at densites larger than a few 0.1cm-3, which was not shocked at velocities larger than 200km/s. [abridged]
△ Less
Submitted 5 April, 2010;
originally announced April 2010.