Sebastien Muller

Molecules as tracers of galaxy evolution: an EMIR survey - I. Presentation of the data and first results

Aims. We investigate the molecular gas properties of a sample of 23 galaxies in order to find and test chemical signatures of galaxy evolution and to compare them to IR evolutionary tracers. Methods. Observation at 3 mm wavelengths were obtained with the EMIR broadband receiver, mounted on the IRAM 30 m telescope on Pico Veleta, Spain. We compare the emission of the main molecular species with existing models of chemical evolution by means of line intensity ratios diagrams and principal component analysis. Results. We detect molecular emission in 19 galaxies in two 8 GHz-wide bands centred at 88 and 112 GHz. The main detected molecules are CO, (CO)-C-13, HCN, HNC, HCO+, CN, and C2H. We also detect HC3N J = 10-9 in the galaxies IRAS 17208, IC 860, NGC 4418, NGC 7771, and NGC 1068. The only HC3N detections are in objects with HCO+/HCN 0.8). The brightest HC3N emission is found in IC 860, where we also detect the molecule in its vibrationally excited state. We find low HNC/HCN line ratios (<0.5), that cannot be explained by existing PDR or XDR chemical models. The intensities of HCO+ and HNC appear anti-correlated. No correlation is found between the HNC/HCN line ratio and dust temperature. All HNC-bright objects are either luminous IR galaxies (LIRG) or Seyferts. Galaxies with bright polycyclic aromatic hydrocarbons (PAH) emission show low HNC/HCO+ ratios. The CO/(CO)-C-13 ratio is positively correlated with the dust temperature and is generally higher than in our galaxy. The emission of CN and (CO)-O-18 is correlated. Conclusions. Bright HC3N emission in HCO+-faint objects may imply that these are not dominated by X-ray chemistry. Thus the HCN/HCO+ line ratio is not, by itself, a reliable tracer of XDRs. Bright HC3N and faint HCO+ could be signatures of embedded star-formation, instead of AGN activity. Mechanical heating caused by supernova explosions may be responsible for the low HNC/HCN and high HCO+/HCN ratios in some starbursts. We cannot exclude, however, that the discussed trends are largely caused by optical depth effects or excitation. Chemical models alone cannot explain all properties of the observed molecular emission. Better constraints to the gas spacial distribution and excitation are needed to distinguish abundance and excitation effects.

Molecules at z = 0.89 - A 4-mm-rest-frame absorption-line survey toward PKS 1830−211

We present the results of a 7 mm spectral survey of molecular absorption lines originating in the disk of a z = 0.89 spiral galaxy located in front of the quasar PKS 1830-211. Our survey was performed with the Australia Telescope Compact Array and covers the frequency interval 30-50 GHz, corresponding to the rest-frame frequency interval 57-94 GHz. A total of 28 different species, plus 8 isotopic variants, were detected toward the south-west absorption region, located about 2 kpc from the center of the z = 0.89 galaxy, which therefore has the largest number of detected molecular species of any extragalactic object so far. The results of our rotation diagram analysis show that the rotation temperatures are close to the cosmic microwave background temperature of 5.14 K that we expect to measure at z = 0.89, whereas the kinetic temperature is one order of magnitude higher, indicating that the gas is subthermally excited. The molecular fractional abundances are found to be in-between those in typical Galactic diffuse and translucent clouds, and clearly deviate from those observed in the dark cloud TMC 1 or in the Galactic center giant molecular cloud Sgr B2. The isotopic ratios of carbon, nitrogen, oxygen, and silicon deviate significantly from the solar values, which can be linked to the young age of the z = 0.89 galaxy and a release of nucleosynthesis products dominated by massive stars. Toward the north-east absorption region, where the extinction and column density of gas is roughly one order of magnitude lower than toward the SW absorption region, only a handful of molecules are detected. Their relative abundances are comparable to those in Galactic diffuse clouds. We also report the discovery of several new absorption components, with velocities spanning between -300 and +170 km s(-1). Finally, the line centroids of several species (e. g., CH3OH, NH3) are found to be significantly offset from the average velocity. If caused by a variation in the proton-to-electron mass ratio mu with redshift, these offsets yield an upper limit vertical bar Delta mu/mu vertical bar \textless 4x10(-6), which takes into account the kinematical noise produced by the velocity dispersion measured from a large number of molecular species.

The density, the cosmic microwave background, and the proton-to-electron mass ratio in a cloud at redshift 0.9

Based on measurements with the Effelsberg 100-m telescope, a multi-line study of molecular species is presented toward the gravitational lens system PKS 1830-211, which is by far the best known target to study dense cool gas in absorption at intermediate redshift. Determining average radial velocities and performing Large Velocity Gradient radiative transfer calculations, the aims of this study are (1) to determine the density of the gas, (2) to constrain the temperature of the cosmic microwave background (CMB), and (3) to evaluate the proton-to-electron mass ratio at redshift z ˜ 0.89. Analyzing data from six rotational HC3N transitions (this includes the J=7≤ftarrow6 line, which is likely detected for the first time in the interstellar medium) we obtain n(H2) ~ 2600 cm-3 for the gas density of the south-western absorption component, assuming a background source covering factor, which is independent of frequency. With a possibly more realistic frequency dependence proportional to ν0.5 (the maximal exponent permitted by observational boundary conditions), n(H2) ~ 1700 cm-3. Again toward the south-western source, excitation temperatures of molecular species with optically thin lines and higher rotational constants are, on average, consistent with the expected temperature of the cosmic microwave background, T_CMB = 5.14 K. However, individually, there is a surprisingly large scatter which far surpasses expected uncertainties. A comparison of CS J =1 ≤ftarrow0 and 4≤ftarrow3 optical depths toward the weaker north-western absorption component results in T_ex = 11 K and a 1-σ error of 3 K. For the main component, a comparison of velocities determined from ten optically thin NH3 inversion lines with those from five optically thin rotational transitions of HC3N, observed at similar frequencies, constrains potential variations of the proton-to-electron mass ratio μ to Δμ / μ < 1.4 × 10 -6 with 3-σ confidence. Also including optically thin rotational lines from other molecular species, it is emphasized that systematic errors are Δ V < 1 km s-1, corresponding to Δμ/μ< 1.0× 10-6.

UVMULTIFIT: A versatile tool for fitting astronomical radio interferometric data

Context. The analysis of astronomical interferometric data is often performed on the images obtained after deconvolving the interferometer’s point spread function. This strategy can be understood (especially for cases of sparse arrays) as fitting models to models, since the deconvolved images are already non-unique model representations of the actual data (i.e., the visibilities). Indeed, the interferometric images may be affected by visibility gridding, weighting schemes (e.g., natural vs. uniform), and the particulars of the (non-linear) deconvolution algorithms. Fitting models to the direct interferometric observables (i.e., the visibilities) is preferable in the cases of simple (analytical) sky intensity distributions. Aims. We present UVMULTIFIT, a versatile library for fitting visibility data, implemented in a Python-based framework. Our software is currently based on the CASA package, but can be easily adapted to other analysis packages, provided they have a Python API. Methods. The user can simultaneously fit an indefinite number of source components to the data, each of which depend on any algebraic combination of fitting parameters. Fits to individual spectral-line channels or simultaneous fits to all frequency channels are allowed. Results. We have tested the software with synthetic data and with real observations. In some cases (e.g., sources with sizes smaller than the diffraction limit of the interferometer), the results from the fit to the visibilities (e.g., spectra of close by sources) are far superior to the output obtained from the mere analysis of the deconvolved images. Conclusions. UVMULTIFIT is a powerful improvement of existing tasks to extract the maximum amount of information from visibility data, especially in cases close to the sensitivity/resolution limits of interferometric observations.

Detection of the SiNC radical in IRC+10216

Following discovery of the free radical SiCN in the C-star envelope IRC+10216, we report the detection in the same source of its isomer SiNC. The microwave spectra of SiNC and SiCN were studied in the laboratory and their rotational transition frequencies are accurately known. The ground fine structure state of SiNC, 2 Π1/2, gives rise to a series of rotational transitions, spaced by 12.8 GHz, each with Λ-doubling. Five weak lines are detected with the IRAM 30-m telescope at the frequencies of the J = 6.5 → 5. 5( e), 7.5 → 6.5 (f), 8.5 → 7.5 (f) and 10.5 → 9.5 (e) and (f) rotational transitions. Other SiNC lines from these or adjacent rotational transitions are found to be blended with stronger lines from known molecules. The lines assigned to SiNC have a cusped shape, characteristic of species confined to a hollow shell in the outer circumstellar envelope. They are twice weaker than their SiCN counterparts, which have the same shape, and presumably arise in the same region of the envelope. SiNC and SiCN have about the same abundance in IRC+10216, ∼4 × 10 −9 with respect to H2 .T his contrasts with HCN, HC3 Na nd HC 5N, for which the cyanide to isocyanide abundance ratio is >100.

Molecular line emission in NGC 1068 imaged with ALMA II. The chemistry of the dense molecular gas

Aims. We present a detailed analysis of Atacama Large Millimeter/submillimeter Array (ALMA) Bands 7 and 9 data of CO, HCO+, HCN, and CS, augmented with Plateau de Bure Interferometer (PdBI) data of the ~200 pc circumnuclear disc (CND) and the ~1.3 kpc starburst ring (SB ring) of NGC 1068, a nearby (D = 14 Mpc) Seyfert 2 barred galaxy. We aim to determine the physical characteristics of the dense gas present in the CND, and to establish whether the different line intensity ratios we find within the CND, as well as between the CND and the SB ring, are due to excitation effects (gas density and temperature differences) or to a different chemistry. Methods. We estimate the column densities of each species in local thermodynamic equilibrium (LTE). We then compute large one-dimensional, non-LTE radiative transfer grids (using RADEX) by using only the CO transitions first, and then all the available molecules to constrain the densities, temperatures, and column densities within the CND. We finally present a preliminary set of chemical models to determine the origin of the gas. Results. We find that, in general, the gas in the CND is very dense (>105 cm-3) and hot (T> 150 K), with differences especially in the temperature across the CND. The AGN position has the lowest CO/HCO+, CO/HCN, and CO/CS column density ratios. The RADEX analyses seem to indicate that there is chemical differentiation across the CND. We also find differences between the chemistry of the SB ring and some regions of the CND; the SB ring is also much colder and less dense than the CND. Chemical modelling does not succeed in reproducing all the molecular ratios with one model per region, suggesting the presence of multi-gas phase components. Conclusions. The LTE, RADEX, and chemical analyses all indicate that more than one gas-phase component is necessary to uniquely fit all the available molecular ratios within the CND. A higher number of molecular transitions at the ALMA resolution is necessary to determine quantitatively the physical and chemical characteristics of these components.

Probing isotopic ratios at z = 0.89: molecular line absorption in front of the quasar PKS 1830-211

With the Plateau de Bure interferometer, we have measured the C, N, O and S isotopic abundance ratios in the arm of a spiral galaxy with a redshift of 0.89. The galaxy is seen face-on according to HST images. Its bulge intercepts the line of sight to the radio-loud quasar PKS 1830−211, giving rise at mm wavelengths to two Einstein images located each behind a spiral arm. The arms appear in absorption in the lines of several molecules, giving the opportunity to study the chemical composition of a galaxy only a few Gyr old. The isotopic ratios in this spiral galaxy differ markedly from those observed in the Milky Way. The 17 O/ 18 Oa nd 14 N/ 15 N ratios are low, as one would expect from an object too young to let low mass stars play a major role in the regeneration of the gas.

A precise and accurate determination of the cosmic microwave background temperature at z =0 .89

Context. According to the Big Bang theory and as a consequence of adiabatic expansion of the Universe, the temperature of the cosmic microwave background (CMB) increases linearly with redshift. This relation is, however, poorly explored, and detection of any deviation would directly lead to (astro-)physics beyond the standard model. Aims. We aim to measure the temperature of the CMB with an accuracy of a few percent at z = 0.89 toward the molecular absorber in the galaxy lensing the quasar PKS 1830−211. Methods. We adopted a Monte-Carlo Markov chain approach, coupled with predictions from the non-LTE radiative transfer code RADEX, to solve the excitation conditions of a set of various molecular species directly from their spectra. Results. We determine TCMB = 5.08 ± 0.10 K at 68% confidence level. Our measurement is consistent with the value TCMB = 5.14 K predicted by the standard cosmological model with adiabatic expansion of the Universe. This is the most precise determination of TCMB at z > 0t o date.

Winds of change – a molecular outflow in NGC 1377? - The anatomy of an extreme FIR-excess galaxy

Aims. Our goal was to investigate the molecular gas distribution and kinematics in the extreme far-infrared (FIR) excess galaxy NGC 1377 and to address the nature and evolutionary status of the buried source. Methods. We used high- (0. �� 65 × 0. �� 52, (65 × 52 pc)) and low- (4. �� 88 × 2. �� 93) resolution SubMillimeter Array (SMA) observations to image the 12 CO and 13 CO 2−1 line emission. Results. We find bright, complex 12 CO 2−1 line emission in the inner 400 pc of NGC 1377. The 12 CO 2−1 line has wings that are tracing a kinematical component that appears to be perpendicular to the component traced by the line core. Together with an intrigu- ing X-shape of the integrated intensity and dispersion maps, this suggests that the molecular emission of NGC 1377 consists of a disk-outflow system. Lower limits to the molecular mass and outflow rate are Mout(H2) > 1 × 10 7 Mand u M > 8 Myr −1 .T he age of the proposed outflow is estimated to be 1.4 Myr, the extent to be 200 pc and the outflow speed to be Vout = 140 km s −1 . The total molecular mass in the SMA map is estimated to Mtot(H2) = 1.5 × 10 8 M� (on a scale of 400 pc) while in the inner r = 29 pc the molecular mass is Mcore(H2) = 1.7 × 10 7 Mwith a corresponding H2 column density of N(H2) = 3.4 × 10 23 cm −2 and an average 12 CO 2−1 brightness temperature of 19 K. 13 CO 2−1 emission is found at a factor 10 fainter than 12 CO in the low-resolution map while C 18 O2 −1 remains undetected. We find weak 1 mm continuum emission of 2.4 mJy with spatial extent less than 400 pc. Conclusions. Observing the molecular properties of the FIR-excess galaxy NGC 1377 allows us to probe the early stages of nuclear activity and the onset of feedback in active galaxies. The age of the outflow supports the notion that the current nuclear activity is young - a few Myr. The outflow may be powered by radiation pressure from a compact, dust enshrouded nucleus, but other driving mechanisms are possible. The buried source may be an active galactic nucleus (AGN) or an extremely young (1 Myr) compact star- burst. Limitations on size and mass lead us to favor the AGN scenario, but additional studies are required to settle this question. In either case, the wind with its implied mass outflow rate will quench the nuclear power source within the very short time of 5−25 Myr. It is possible, however, that the gas is unable to escape the galaxy and may eventually fall back onto NGC 1377 again.

A strong magnetic field in the jet base of a supermassive black hole

The polarized mark of magnetic fields Powerful twin jets of plasma often reach more than tens of thousands of light-years from their base in an active galactic nucleus (AGN). Astronomers are still at work investigating what can corral the jets so tightly and propel them so far. Martí-Vidal et al. may have found the answer hiding in polarized light signals that show evidence of a phenomenon called Faraday rotation. This measure can indicate the strength of the magnetic field present, which for the AGN PKS 1830-211 is as strong as a few Gauss. The knowledge that magnetic fields have a driving role brings us closer to understanding this phenomenon. Science, this issue p. 311 A polarized signal offers evidence for the agent that boosts and guides powerful jets in a distant active galaxy. Active galactic nuclei (AGN) host some of the most energetic phenomena in the universe. AGN are thought to be powered by accretion of matter onto a rotating disk that surrounds a supermassive black hole. Jet streams can be boosted in energy near the event horizon of the black hole and then flow outward along the rotation axis of the disk. The mechanism that forms such a jet and guides it over scales from a few light-days up to millions of light-years remains uncertain, but magnetic fields are thought to play a critical role. Using the Atacama Large Millimeter/submillimeter Array (ALMA), we have detected a polarization signal (Faraday rotation) related to the strong magnetic field at the jet base of a distant AGN, PKS 1830−211. The amount of Faraday rotation (rotation measure) is proportional to the integral of the magnetic field strength along the line of sight times the density of electrons. The high rotation measures derived suggest magnetic fields of at least tens of Gauss (and possibly considerably higher) on scales of the order of light-days (0.01 parsec) from the black hole.