H2 formation and excitation in the Stephan's Quintet galaxy-wide collision

Context. The Spitzer Space Telescope has detected a powerful (L(H2)~10^41 erg s-1) mid-infrared H2 emission towards the galaxy-wide collision in the Stephan's Quintet (SQ) galaxy group. This discovery was followed by the detection of more distant H2-luminous extragalactic sources, with almost no spectroscopic signatures of star formation. These observations set molecular gas in a new context where one has to describe its role as a cooling agent of energetic phases of galaxy evolution. Aims. The SQ postshock medium is observed to be multiphase, with H2 gas coexisting with a hot (~ 5 10^6 K), X-ray emitting plasma. The surface brightness of H2 lines exceeds that of the X-rays and the 0-0 S(1) H2 linewidth is ~ 900 km s-1, of the same order of the collision velocity. These observations raise three questions we propose to answer: (i) Why H2 is present in the postshock gas ? (ii) How can we account for the H2 excitation ? (iii) Why H2 is a dominant coolant ? Methods. We consider the collision of two flows of multiphase dusty gas. Our model quantifies the gas cooling, dust destruction, H2 formation and excitation in the postshock medium. Results. (i) The shock velocity, the post-shock temperature and the gas cooling timescale depend on the preshock gas density. The collision velocity is the shock velocity in the low density volume filling intercloud gas. This produces a ~ 5 10^6 K, dust-free, X-ray emitting plasma. The shock velocity is smaller in clouds. We show that gas heated to temperatures less than 10^6 K cools, keeps its dust content and becomes H2 within the SQ collision age (~ 5 10^6 years). (ii) Since the bulk kinetic energy of the H2 gas is the dominant energy reservoir, we consider that the H2 emission is powered by the dissipation of kinetic turbulent energy. (Abridged)Comment: 19 pages, 12 figures. Accepted for publication in Astronomy & Astrophysics Minor editing and typo

Temporal evolution of magnetic molecular shocks I. Moving grid simulations

We present time-dependent 1D simulations of multifluid magnetic shocks with chemistry resolved down to the mean free path. They are obtained with an adaptive moving grid implemented with an implicit scheme. We examine a broad range of parameters relevant to conditions in dense molecular clouds, with preshock densities between 10^3 and 10^5 cm-3, velocities between 10 and 40 km/s, and three different scalings for the transverse magnetic field: B=0,0.1,1 \mu G \sqrt{n.cm3}. We first use this study to validate the results of Chi\`eze, Pineau des For\^ets & Flower (1998), in particular the long delays necessary to obtain steady C-type shocks, and we provide evolutionary time-scales for a much greater range of parameters. We also present the first time-dependent models of dissociative shocks with a magnetic precursor, including the first models of stationary CJ shocks in molecular conditions. We find that the maximum speed for steady C-type shocks is reached before the occurrence of a sonic point in the neutral fluid, unlike previously thought. As a result, the maximum speed for C-shocks is lower than previously believed. Finally, we find a large amplitude bouncing instability in J-type fronts near the H2 dissociation limit (u ~ 25-30 km/s), driven by H2 dissociation/reformation. At higher speeds, we find an oscillatory behaviour of short period and small amplitude linked to collisional ionisation of H. Both instabilities are suppressed after some time when a magnetic field is present. In a companion paper, we use the present simulations to validate a new semi-analytical construction method for young low-velocity magnetic shocks based on truncated steady-state models.Comment: A&A in pres

Dissipative structures of diffuse molecular gas: I - Broad HCO+^+(1-0) emission

Results: We report the detection of broad HCO+(1-0) lines (10 mK < T < 0.5 K). The interpretation of 10 of the HCO+ velocity components is conducted in conjunction with that of the associated optically thin 13CO emission. The derived HCO+ column densities span a broad range, $10^{11}< N(HCO+)/\Delta v <4 \times 10^{12} \rm cm^2/(km/s^{-1}$, and the inferred HCO+ abundances, $2 \times 10^{-10}<X(HCO+) < 10^{-8}$, are more than one order of magnitude above those produced by steady-state chemistry in gas weakly shielded from UV photons, even at large densities. We compare our results with the predictions of non-equilibrium chemistry, swiftly triggered in bursts of turbulence dissipation and followed by a slow thermal and chemical relaxation phase, assumed isobaric. The set of values derived from the observations, i.e. large HCO+ abundances, temperatures in the range of 100--200 K and densities in the range 100--1000 cm3, unambiguously belongs to the relaxation phase. The kinematic properties of the gas suggest in turn that the observed HCO+ line emission results from a space-time average in the beam of the whole cycle followed by the gas and that the chemical enrichment is made at the expense of the non-thermal energy. Last, we show that the "warm chemistry" signature (i.e large abundances of HCO+, CH+, H20 and OH) acquired by the gas within a few hundred years, the duration of the impulsive chemical enrichment, is kept over more than thousand years. During the relaxation phase, the \wat/OH abundance ratio stays close to the value measured in diffuse gas by the SWAS satellite, while the OH/HCO+ ratio increases by more than one order of magnitude.Comment: 14 page

H_2 formation and excitation in the Stephan's Quintet galaxy-wide collision

Context. The Spitzer Space Telescope has detected a powerful (L_(H_2) ~ 10^(41) erg s^(-1)) mid-infrared H_2 emission towards the galaxy-wide collision in the Stephan's Quintet (henceforth SQ) galaxy group. This discovery was followed by the detection of more distant H_2-luminous extragalactic sources, with almost no spectroscopic signatures of star formation. These observations place molecular gas in a new context where one has to describe its role as a cooling agent of energetic phases of galaxy evolution. Aims. The SQ postshock medium is observed to be multiphase, with H_2 gas coexisting with a hot (~5 × 10^6 K), X-ray emitting plasma. The surface brightness of H_2 lines exceeds that of the X-rays and the 0-0 S(1)H_2 linewidth is ~900 km  s^(-1), of the order of the collision velocity. These observations raise three questions we propose to answer: (i) why is H_2 present in the postshock gas? (ii) How can we account for the H_2 excitation? (iii) Why is H_2 a dominant coolant? Methods. We consider the collision of two flows of multiphase dusty gas. Our model quantifies the gas cooling, dust destruction, H_2 formation and excitation in the postshock medium. Results. (i) The shock velocity, the post-shock temperature and the gas cooling timescale depend on the preshock gas density. The collision velocity is the shock velocity in the low density volume-filling intercloud gas. This produces a ~5 × 10^6 K, dust-free, X-ray emitting plasma. The shock velocity is lower in clouds. We show that gas heated to temperatures of less than 10^6 K cools, keeps its dust content and becomes H_2 within the SQ collision age (~5 × 10^6 years). (ii) Since the bulk kinetic energy of the H_2 gas is the dominant energy reservoir, we consider that the H_2 emission is powered by the dissipation of kinetic turbulent energy. We model this dissipation with non-dissociative MHD shocks and show that the H_2 excitation can be reproduced by a combination of low velocities shocks (5-20 km s^(-1)) within dense (n_H > 10^3 cm^(-3)) H_2 gas. (iii) An efficient transfer of the bulk kinetic energy to turbulent motion of much lower velocities within molecular gas is required to make H_2 a dominant coolant of the postshock gas. We argue that this transfer is mediated by the dynamic interaction between gas phases and the thermal instability of the cooling gas. We quantify the mass and energy cycling between gas phases required to balance the dissipation of energy through the H_2 emission lines. Conclusions. This study provides a physical framework to interpret H_2 emission from H_2-luminous galaxies. It highlights the role that H_2 formation and cooling play in dissipating mechanical energy released in galaxy collisions. This physical framework is of general relevance for the interpretation of observational signatures, in particular H_2 emission, of mechanical energy dissipation in multiphase gas

Observations and modeling of the dust emission from the H_2-bright galaxy-wide shock in Stephan's Quintet

Context. Spitzer Space Telescope observations have detected powerful mid-infrared (mid-IR) H_2 rotational line emission from the X-ray emitting large-scale shock (~15 × 35 kpc^2) associated with a galaxy collision in Stephan's Quintet (SQ). Because H_2 forms on dust grains, the presence of H_2 is physically linked to the survival of dust, and we expect some dust emission to originate in the molecular gas. Aims. To test this interpretation, IR observations and dust modeling are used to identify and characterize the thermal dust emission from the shocked molecular gas. Methods. 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. The radiation field is determined with GALEX UV, HST V-band, and ground-based near-IR observations. We consider two limiting cases for the structure of the H_2 gas: it is either diffuse and penetrated by UV radiation, or fragmented into clouds that are optically thick to UV. Results. Faint PAH and dust continuum emission are detected in the SQ shock, outside star-forming regions. The 12/24 μm 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 shock inferred from the PAH emission spectrum differ from those of the Galaxy, which may be indicative of an enhanced fraction of large and neutrals PAHs. In both models (diffuse or clumpy H_2 gas), the IR SED is consistent with the expected emission from dust associated with the warm (> 150 K) H_2 gas, heated by a UV radiation field of intensity comparable to that of the solar neighborhood. This is in agreement with GALEX UV observations that show that the intensity of the radiation field in the shock is GUV = 1.4±0.2 [Habing units]. Conclusions. The presence of PAHs and dust grains in the high-speed (~1000 km s^(-1)) galaxy collision suggests that dust survives. We propose that the dust that survived destruction was in pre-shock gas at densites higher than a few 0.1 cm^(-3), which was not shocked at velocities larger than ~200 km s^(-1). Our model assumes a Galactic dust-to-gas mass ratio and size distribution, and current data do not allow us to identify any significant deviations of the abundances and size distribution of dust grains from those of the Galaxy. Our model calculations show that far-IR Herschel observations will help in constraining the structure of the molecular gas, and the dust size distribution, and thereby to look for signatures of dust processing in the SQ shock

Dense molecular globulettes and the dust arc towards the runaway O star AE Aur (HD 34078)

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 {\mu}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)Comment: Accepted for publication in Astronomy & Astrophysic

Energetics of the molecular gas in the H_2 luminous radio galaxy 3C 326: Evidence for negative AGN feedback

We present a detailed analysis of the gas conditions in the H_2 luminous radio galaxy 3C 326 N at z ~ 0.1, which has a low star-formation rate (SFR ~ 0.07 M_⊙ yr^(−1)) in spite of a gas surface density similar to those in starburst galaxies. Its star-formation efficiency is likely a factor ~ 10−50 lower than those of ordinary star-forming galaxies. Combining new IRAM CO emission-line interferometry with existing Spitzer mid-infrared spectroscopy, we find that the luminosity ratio of CO and pure rotational H_2 line emission is factors 10−100 lower than what is usually found. This suggests that most of the molecular gas is warm. The Na D absorption-line profile of 3C 326 N in the optical suggests an outflow with a terminal velocity of ~−1800 km s^(−1) and a mass outflow rate of 30−40 M_⊙ yr^(−1), which cannot be explained by star formation. The mechanical power implied by the wind, of order 10^(43) erg s^(−1), is comparable to the bolometric luminosity of the emission lines of ionized and molecular gas. To explain these observations, we propose a scenario where a small fraction of the mechanical energy of the radio jet is deposited in the interstellar medium of 3C 326 N, which powers the outflow, and the line emission through a mass, momentum and energy exchange between the different gas phases of the ISM. Dissipation times are of order 10^(7−8) yrs, similar or greater than the typical jet lifetime. Small ratios of CO and PAH surface brightnesses in another 7 H_2 luminous radio galaxies suggest that a similar form of AGN feedback could be lowering star-formation efficiencies in these galaxies in a similar way. The local demographics of radio-loud AGN suggests that secular gas cooling in massive early-type galaxies of ≥ 10^(11) M_⊙ could generally be regulated through a fundamentally similar form of “maintenance-phase” AGN feedback

Collisional excitation of water by hydrogen atoms

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

Nitrogen superfractionation in dense cloud cores

We report new calculations of interstellar 15N fractionation. Previously, we have shown that large enhancements of 15N/14N can occur in cold, dense gas where CO is frozen out, but that the existence of an NH + N channel in the dissociative recombination of N2H+ severely curtails the fractionation. In the light of recent experimental evidence that this channel is in fact negligible, we have reassessed the 15N chemistry in dense cloud cores. We consider the effects of temperatures below 10 K, and of the presence of large amounts of atomic nitrogen. We also show how the temporal evolution of gas-phase isotope ratios is preserved as spatial heterogeneity in ammonia ice mantles, as monolayers deposited at different times have different isotopic compositions. We demonstrate that the upper layers of this ice may have 15N/14N ratios an order of magnitude larger than the underlying elemental value. Converting our ratios to delta-values, we obtain delta(15N) > 3,000 per mil in the uppermost layer, with values as high as 10,000 per mil in some models. We suggest that this material is the precursor to the 15N `hotspots' recently discovered in meteorites and IDPsComment: accepted by MNRA

ISOCAM spectro-imaging of the H2 rotational lines in the supernova remnant IC443

We report spectro-imaging observations of the bright western ridge of the supernova remnant IC 443 obtained with the ISOCAM circular variable filter (CVF) on board the Infrared Space Observatory (ISO). This ridge corresponds to a location where the interaction between the blast wave of the supernova and ambient molecular gas is amongst the strongest. The CVF data show that the 5 to 14 micron spectrum is dominated by the pure rotational lines of molecular hydrogen (v = 0--0, S(2) to S(8) transitions). At all positions along the ridge, the H2 rotational lines are very strong with typical line fluxes of 10^{-4} to 10^{-3} erg/sec/cm2/sr. We compare the data to a new time-dependent shock model; the rotational line fluxes in IC 443 are reproduced within factors of 2 for evolutionary times between 1,000 and 2,000 years with a shock velocity of 30 km/sec and a pre-shock density of 10^4 /cm3.Comment: To appear in Astronomy and Astrophysic