Source author record

F. Combes

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

199works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

199 published item(s)

preprint2025arXiv

High resolution mapping of molecular tori with ALMA

Recent high resolution mapping of the circum-nuclear regions of Active Galactic Nuclei (AGN) has revealed the existence of geometrically thin nuclear disks, in general randomly oriented with respect to their galaxy hosts. These molecular tori have typical radii of 10~pc, and contain a few 10$^7$ M$_\odot$ of H$_2$, with H$_2$ column densities between 10$^{23}$ and 10$^{25}$ cm$^{-2}$. We mapped two of the most massive of these molecular tori with higher resolution, in order to unveil their morphology and kinematics, their possible warp and clumpiness, and derive their stability and life-time. We used the highest resolution possible with ALMA (16~km baseline) in Band 7, taking into account for mapping CO(3-2) and HCO$^+$(4-3) the compromise between sensitivity and resolution. New features are discovered at the high resolution, obtained with a beam of 0.015\arcsec, equivalent to $\sim$ 1~pc scale, at their $\sim$ 15~Mpc distance. The molecular torus in NGC~613 appears like a ring, depleted in molecular gas near the center. The depletion region is displaced by 3~pc toward the NW from the AGN position, meaning some $m=1$ asymmetry in the torus. The molecular torus in NGC~1672 has a different position angle from previous lower-resolution observations, and is edge-on, revealing a geometrically very thin torus (axis ratio 6.5 to 10), with a clear warp. This confirms that the classical model of a simple geometrically thick dusty torus is challenged by high resolution observations. The nuclear disks appear clumpy, and slightly lopsided. The molecular outflow in NGC~613 is now resolved out. Well inside the sphere of influence of the black holes (BH), we are now able to determine more accurately their mass, for those Seyfert spiral galaxies, in a region of the M-sigma relation where the scatter is maximum.

preprint2022arXiv

First discoveries and localisations of Fast Radio Bursts with MeerTRAP: a real-time, commensal MeerKAT survey

We report on the discovery and localization of fast radio bursts (FRBs) from the MeerTRAP project, a commensal fast radio transient-detection programme at MeerKAT in South Africa. Our hybrid approach combines a coherent search with an average field-of-view of 0.4 $\rm deg^{2}$ with an incoherent search utilizing a field-of-view of $\sim$1.27 $\rm deg^{2}$ (both at 1284~MHz). Here, we present results on the first three FRBs: FRB 20200413A (DM=1990.05 pc cm$^{-3}$), FRB 20200915A (DM=740.65 pc cm$^{-3}$), and FRB 20201123A (DM=433.55 pc cm$^{-3}$). FRB 20200413A was discovered only in the incoherent beam. FRB 20200915A (also discovered only in the incoherent beam) shows speckled emission in the dynamic spectrum which cannot be explained by interstellar scintillation in our Galaxy or plasma lensing, and might be intrinsic to the source. FRB 20201123A shows a faint post-cursor burst about 200 ms after the main burst and warrants further follow-up to confirm whether it is a repeating FRB. FRB 20201123A also exhibits significant temporal broadening consistent with scattering by a turbulent medium. The broadening exceeds that predicted for medium along the sightline through our Galaxy. We associate this scattering with the turbulent medium in the environment of the FRB in the host galaxy. Within the approximately $1'$ localization region of FRB 20201123A, we identify one luminous galaxy ($r \approx 15.67$; J173438.35$-$504550.4) that dominates the posterior probability for a host association. The galaxy's measured properties are consistent with other FRB hosts with secure associations.

preprint2022arXiv

HI gas playing hide-and-seek around a powerful FRI-type quasar at z$\sim$2.1

We present optical spectroscopic and milli-arcsecond scale radio continuum observations of the quasar M1540-1453 ($z_{em}$ = 2.104$\pm$0.002) that shows associated HI 21-cm absorption at $z_{abs}$ = 2.1139. At sub-kpc scales, the powerful radio source with 1.4 GHz luminosity of $5.9\times10^{27}$ WHz$^{-1}$ shows Fanaroff-Riley (FR) class I morphology caused by the interaction with dense gas within 70 pc from the AGN. Interestingly, while there are indications for the presence of absorption from low-ionization species like FeII, SiII and SiIII in the optical spectrum, the expected strong damped Ly$α$ absorption is not detected at the redshift of the HI 21-cm absorber. In comparison to typical high-$z$ quasars, the Ly$α$ emission line is much narrower. The `ghostly' nature of the HI Ly$α$ absorber partially covering the broad line region of extent 0.05 pc and the detection of widespread HI 21-cm absorption covering the diffuse radio source (extent $>$425 pc) imply the presence of a large clumpy HI halo -- which may have been blown by the jet-ISM interaction. Further observations are needed to confirm the `ghostly' nature of the Ly$α$ absorber, and obtain a better understanding of the role played by the jet-ISM interaction in shaping the radio morphology of this powerful AGN. The study showcases how joint radio and optical analysis can shed light on gaseous environment and origin of radio morphology in AGN at high redshifts, when these are still the assembly sites of giant galaxies.

preprint2022arXiv

Linking pressure gradients with the stability of molecular clouds in galactic outflows

The jets launched by actively accreting black holes are capable of launching several of the massive (million or billion solar mass) molecular outflows observed in galaxies. These outflows could suppress or enhance star formation in galaxies. To investigate the stability of clouds capable to form stars in outflows, we modeled CO and HCO+ ALMA data of the galaxy IC5063, in which black-hole jets impact molecular clouds. Using a radiative transfer code that self-consistently performs astrochemical and thermal balance calculations based on the available gas heating sources, we found that mechanical heating and cosmic ray (CR) heating are fully capable of individually reproducing the data. In our best-fit model, CRs provide about 1/3rd of the dense gas heating at the radio lobes, emphasizing the role of this often neglected mechanism in heating the gas and potentially generating outflows. The gas temperature and density indicate that the jet passage leads to an increase of about 1 order of magnitude in the internal pressure Pi of molecular clouds (with Pi/k from 8*10^5 up to 7*10^6 K cm^-3), irrespective of the excitation mechanism. From the fluxes of [S II] and [N II] lines in VLT MUSE data, the external pressure Pe of molecular clouds increases in several regions enough to exceed Pi. This result leads us to conclude that we are observing the expansion of an ionized overpressurized cocoon that compresses molecular clouds and that could lead to their collapse. Some jet-impacted clouds, nonetheless, near pathways that the jet cleared have increased Pi and decreased Pe. They are likely to undergo evaporation of their outer layers. Part of the evaporated layers could mass load the outflow thanks to ram pressure from co-spatial ionized gas flows. The observed pressure changes thus suggest that both star formation enhancement and suppression could simultaneously occur.

preprint2022arXiv

MALS SALT-NOT survey of MIR-selected powerful radio-bright AGN at 0<z<3.5

We present results of an optical spectroscopic survey using SALT and NOT to build a WISE mid-infrared color-based, dust-unbiased sample of powerful radio-bright ($>$200 mJy at 1.4 GHz) AGN for the MeerKAT Absorption Line Survey (MALS). Our sample has 250 AGN (median $z=1.8$) showing emission lines, 26 with no emission lines, and 27 without optical counterparts. Overall, our sample is fainter ($Δi$=0.6 mag) and redder ($Δ(g-i)$=0.2 mag) than radio-selected quasars, and representative of fainter quasar population detected in optical surveys. About 20% of the sources are narrow line AGN (NLAGN) $-$ 65% of these, at $z < 0.5$ are galaxies without strong nuclear emission, and 10% at $z>1.9$, have emission line ratios similar to radio galaxies. The farthest NLAGN in our sample is M1513$-$2524 ($z_{em}=3.132$), and the largest (size$\sim$330 kpc) is M0909$-$3133 ($z_{em}=0.884$). We discuss in detail 110 AGN at $1.9 < z < 3.5$. Despite representing the radio loudest quasars (median $R$=3685), their Eddington ratios are similar to the SDSS quasars having lower $R$. We detect 4 CIV BALQSOs, all among AGN with least $R$, and highest black hole masses and Eddington ratios. The BAL detection rate ($4^{+3}_{-2}$%) is consistent with that seen in extremely powerful ($L_{1.4GHz}>10^{25}$ WHz$^{-1}$) quasars. Using optical light-curves, radio polarization and $γ$-ray detections, we identify 7 high-probability BL Lacs. We also summarize the full MALS footprint to search for HI 21-cm and OH 18-cm lines at $z<2$.

preprint2022arXiv

Mapping HI 21-cm in the Klemola 31 group at z = 0.029: emission and absorption towards PKS2020-370

We present MeerKAT Absorption Line Survey (MALS) observations of the HI gas in the Klemola31 galaxy group ($z=0.029$), located along the line of sight to the radio-loud quasar PKS2020-370 ($z=1.048$). Four galaxies of the group are detected in HI emission, and HI absorption is also detected in front of PKS2020-370 in Klemola31A. The emission and absorption are somewhat compensating on the line of sight of the quasar, and the derived column density of the absorption appears under-estimated, with respect to the neighbouring emission. A symmetric tilted-ring model of Klemola31A, assuming the absorbing gas in regular rotation in the plane, yields a rather high spin temperature of 530K. An alternative interpretation is that the absorbing gas is extra-planar, which will also account for its non-circular motion. The NaI/CaII ratio also suggests that the absorbing gas is unrelated to cold HI disk. Two of the galaxies in the Klemola group are interacting with a small companion, and reveal typical tidal tails, and velocity perturbations. Only one of the galaxies, ESO400-13, reveals a strong HI deficiency, and a characteristic ram-pressure stripping, with a total asymmetry in the distribution of its gas. Since a small galaxy group as Klemola31 is not expected to host a dense intra-group gas, this galaxy must be crossing the group at a very high velocity, mostly in the sky plane.

preprint2022arXiv

NOEMA observations support a recoiling black hole in 3C186

3C186 is a powerful radio loud quasar (QSO) at the center of a cool-core cluster at z=1.06. Previous studies reported evidence for a projected spatial offset of ~1'' between the isophotal center of the galaxy and the point-source QSO as well as a spectral shift of ~2000 km/s between the narrow and broad line region of the system. In this work we report high-resolution molecular gas CO(4-3) observations of the system taken with NOEMA interferometer. We clearly detect a large reservoir of molecular gas, $M_{H_2}\sim8\times10^{10}~M_\odot$, that is cospatial with the host galaxy and likely associated with a rotating disk-like structure. We firmly confirm both the spatial offset of the galaxy's gas reservoir with respect to the continuum emission of the QSO and the spectral offset with respect to the redshift of the broad line region. Our morphological and kinematical analysis confirms that the most likely scenario to explain the 3C186 system is that the QSO is a kicked super-massive black hole (SMBH), which we believe may have resulted from a strong gravitational wave recoil as two SMBHs coalesced after the merger of their host galaxies.

preprint2022arXiv

The Close AGN Reference Survey (CARS): Tracing the circumnuclear star formation in the super-Eddington NLS1 Mrk 1044

The host galaxy conditions for rapid supermassive black hole growth are poorly understood. Narrow-line Seyfert 1 (NLS1) galaxies often exhibit high accretion rates and are hypothesized to be prototypes of active galactic nuclei (AGN) at an early stage of their evolution. We present VLT MUSE NFM-AO observations of Mrk 1044, the nearest super-Eddington accreting NLS1. Together with archival MUSE WFM data we aim to understand the host galaxy processes that drive Mrk 1044's black hole accretion. We extract the faint stellar continuum emission from the AGN-deblended host and perform spatially resolved emission line diagnostics with an unprecedented resolution. Combining both MUSE WFM and NFM-AO observations, we use a kinematic model of a thin rotating disk to trace the stellar and ionized gas motion from 10$\,$kpc down to 30$\,$pc around the nucleus. Mrk 1044's stellar kinematics follow circular rotation, whereas the ionized gas shows tenuous spiral features in the center. We resolve a compact star forming circumnuclear ellipse (CNE) that has a semi-minor axis of 306$\,$pc. Within this CNE, the gas is metal rich and its line ratios are entirely consistent with excitation by star formation. With an integrated SFR of $0.19 \pm 0.05 \,{\rm M}_\odot\,{\rm yr}^{-1}$ the CNE contributes 27% of the galaxy-wide star formation. We conclude that Mrk 1044's nuclear activity has not yet affected the circumnuclear star formation. Thus, Mrk 1044 is consistent with the idea that NLS1s are young AGN. A simple mass budget consideration suggests that the circumnuclear star formation and AGN phase are connected and the patterns in the ionized gas velocity field are a signature of the ongoing AGN feeding.

preprint2022arXiv

Unveiling the Main Sequence to Starburst Transition Region with a Sample of Intermediate Redshift Luminous Infrared Galaxies

We present a CO(3-2) study of four systems composed of six (ultra) luminous infrared galaxies (U/LIRGs), located at 0.28 < z < 0.44, that straddle the transition region between regular star forming galaxies and starbursts. These galaxies benefit from previous multi-wavelength analysis allowing in depth exploration of an understudied population of U/LIRGs at a time when the universe is experiencing a rapid decline in star formation rate density. We detect CO(3-2) emission in four targets and these galaxies fall between the loci of regular star forming galaxies and starbursts on the Kennicutt-Schmidtt relation. Compared to low luminosity LIRGs and high luminosity ULIRGs at similar redshifts, we find they all have similar molecular gas budgets with the difference in their star formation rates (SFR) driven by the star formation efficiency (SFE). This suggests that at these redshifts large molecular gas reservoirs must coincide with an increased SFE to transition a galaxy into the starburst regime. We studied the structure and kinematics and found our four detections are either interacting or have disturbed morphology which may be driving the SFE. One of the CO(3-2) non-detections has a strong continuum detection, and has been previously observed in H$α$, suggesting an unusual interstellar medium for a ULIRG. We conclude that our sample of transitioning U/LIRGs fill the gap between regular star forming galaxies and starbursts, suggest a continuous change in SFE between these two populations and the increased SFE may be driven by morphology and differing stages of interaction.

preprint2022arXiv

Virgo Filaments II: Catalog and First Results on the Effect of Filaments on galaxy properties

Virgo is the nearest galaxy cluster; it is thus ideal for studies of galaxy evolution in dense environments in the local Universe. It is embedded in a complex filamentary network of galaxies and groups, which represents the skeleton of the large scale Laniakea supercluster. Here we assemble a comprehensive catalog of galaxies extending up to ~12 virial radii in projection from Virgo to revisit the Cosmic Web structure around it. This work is the foundation of a series of papers that will investigate the multi-wavelength properties of galaxies in the Cosmic Web around Virgo. We match spectroscopically confirmed sources from several databases and surveys including HyperLeda, NASA Sloan Atlas, NED, and ALFALFA. The sample consists of ~7000 galaxies. By exploiting a tomographic approach, we identify 13 filaments, spanning several Mpc in length. Long > 17 Mpc/h filaments, tend to be thin (< 1 Mpc/h in radius) and with a low density contrast (< 5), while shorter filaments show a larger scatter in their structural properties. Overall, we find that filaments are a transitioning environment between the field and cluster in terms of local densities, galaxy morphologies, and fraction of barred galaxies. Denser filaments have a higher fraction of early type galaxies, suggesting that the morphology-density relation is already in place in the filaments, before galaxies fall into the cluster itself. We release the full catalog of galaxies around Virgo and their associated properties.

preprint2021arXiv

Dense Molecular Clouds in the Crab Supernova Remnant

Molecular emission was imaged with ALMA from numerous components near and within bright H2-emitting knots and absorbing dust globules in the Crab Nebula. These observations provide a critical test of how energetic photons and particles produced in a young supernova remnant interact with gas, cleanly differentiating between competing models. The four fields targeted show contrasting properties but within them, seventeen distinct molecular clouds are identified with CO emission; a few also show emission from HCO+, SiO and/or SO. These observations are compared with Cloudy models of these knots. It has been suggested that the Crab filaments present an exotic environment in which H2 emission comes from a mostly-neutral zone probably heated by cosmic rays produced in the supernova surrounding a cool core of molecular gas. Our model is consistent with the observed CO J=3-2 line strength. These molecular line emitting knots in the Crab present a novel phase of the ISM representative of many important astrophysical environments.

preprint2021arXiv

The Close AGN Reference Survey (CARS): IFU survey data and the BH mass dependence of long-term AGN variability

[Abridged] AGN are thought to be intimately connected with their host galaxies through feeding and feedback processes. A spatially resolved multiwavelength survey is required to map the interaction of AGN with their host galaxies on different spatial scales and different phases of the ISM. The goal of CARS is to obtain the necessary spatially resolved multiwavelength observations for an unbiased sample of local unobscured luminous AGN. We present the overall CARS survey design and the associated wide-field optical IFU spectroscopy for all 41 CARS targets at z<0.06 randomly selected from the Hamburg/ESO survey of luminous unobscured AGN. This data set provides the backbone of CARS and allows us to characterize host galaxy morphologies, AGN parameters, precise systemic redshifts, and ionized gas distributions including excitation conditions, kinematics, and metallicities in unprecedented detail. We focus our study on the size of the ENLR which has been traditionally connected to AGN luminosity. Given the large scatter in the ENLR size-luminosity relation, we performed a large parameter search to identify potentially more fundamental relations. Remarkably, we identified the strongest correlation between the maximum projected ENLR size and the black hole mass, consistent with an $R_\mathrm{ENLR,max}\sim M_\mathrm{BH}^{0.5}$ relationship. We interpret the maximum ENLR size as a timescale indicator of a single BH radiative-efficient accretion episode for which we inferred log(t_AGN) = (0.45+- 0.08)log(M_BH)+1.78 using forward modeling. The extrapolation of our inferred relation toward higher BH masses is consistent with an independent lifetime estimate from the HeII proximity zones around luminous AGN at z~3. While our proposed link between the BH mass and AGN lifetime might be a secondary correlation itself or impacted by unknown biases, it has a few relevant implications if confirmed.

preprint2021arXiv

The Close AGN Reference Survey (CARS): Locating the [O III] wing component in luminous local Type 1 AGN

[Abridged]The strong asymmetry in the optical [O III]$λ$5007 emission line is one of the best signatures of AGN-driven warm (~10$^4$ K) ionized gas outflows on host galaxy scales. While large spectroscopic surveys like SDSS have characterized the kinematics of [O III] for large samples of AGN, estimating the associated energetics requires spatially resolving these outflows with, for example, IFU studies. As part of CARS we obtained spatially-resolved IFU spectroscopy for a representative sample of 39 luminous type 1 AGN at 0.01<z<0.06 with MUSE and VIMOS IFUs at the VLT to infer the spatial location of the ionized gas outflows. We compare the light distributions of the [O III] wing to that of the H$β$ broad emission line region, a classical point source (PSF). We then use the PSF to distinguish between the unresolved and resolved [O III] wing emission. We further determine its location using spectro-astrometry for the point-like sources. The [O III] wing is spatially unresolved in 23 out of the 36 AGN with >80 % of the flux associated with a point-like source. We measure <100 pc offsets in the spatial location of the outflow from the AGN nucleus using the spectro-astrometry technique for these sources. For the other 13 AGN, the [O III] wing emission is resolved and possibly extended on kpc scale. We conclude that [O III] wing emission can be compact or extended in an unbiased luminous AGN sample, where both cases are likely to appear. Electron density in the compact [O III] wing regions (median $n_e$~1900 cm$^{-3}$) is nearly a magnitude higher than in the extended ones (median $n_e$~500 cm$^{-3}$). The presence of spatially extended and compact [O III] wing emission is unrelated to the AGN bolometric luminosity and to inclination effects, which means other features such as time delays, or mechanical feedback/radio jets may shape the ionized gas outflow properties.

preprint2021arXiv

The Close AGN Reference Survey (CARS): No obvious signature of AGN feedback on star formation, but subtle trends

[Abridged] Active Galactic Nuclei (AGN) are thought to be responsible for the suppression of star formation in massive ~10$^{10}$ M$_\odot$ galaxies. While this process is a key feature in numerical simulations, it is not yet unambiguously confirmed in observational studies. Characterization of the star formation rate (SFR) in AGN host galaxies is challenging as AGN light contaminates most SFR tracers. We aim to obtain and compare SFR estimates from different tracers for AGN host galaxies in the Close AGN Reference Survey (CARS) to provide new observational insights. We construct integrated panchromatic spectral energy distributions (SED) to measure the FIR luminosity as a tracer for the recent (< 100 Myr) SFR. In addition, we use integral-field unit observation of the CARS targets to employ the H$α$ luminosity decontaminated by AGN excitation as a proxy for the current (< 5 Myr) SFR. We find that significant differences in specific SFR of the AGN host galaxies as compared with the larger galaxy population disappear once cold gas mass, in addition to stellar mass, is used to predict the SFR. We identify individual galaxies with a significant difference in their SFR which can be related to a recent enhancement or decline in their SFR history that might be related to various processes including interactions, gas consumption, outflows and AGN feedback. AGN can occur in various stages of galaxy evolution which makes it difficult to relate the SFR solely to the impact of the AGN. We do not find any strong evidence for global positive or negative AGN feedback in the CARS sample. However, there is tentative evidence that 1) the relative orientation of the AGN engine with respect to the host galaxies might alter the efficiency of AGN feedback and 2) the recent SFH is an additional tool to identify rapid changes in galaxy growth driven by the AGN or other processes.

preprint2021arXiv

Virgo Filaments I: Processing of gas in cosmological filaments around the Virgo cluster

Galaxies have different morphology, gas content and star formation rate (SFR) in dense environments like galaxy clusters. The impact of environmental density extends to several virial radii, and galaxies are pre-processed in filaments and groups, before falling into the cluster. Our goal is to quantify this pre-processing, in terms of gas content and SFR, as a function of density in cosmic filaments. We have observed the two first CO transitions in 163 galaxies with the IRAM-30m telescope, and added 82 more measurements from the literature, for a sample of 245 galaxies in the filaments around Virgo. We gathered HI-21cm measurements from the literature, and observed 69 galaxies with the Nançay telescope, to complete our sample. We compare our filament galaxies with comparable samples from the Virgo cluster and with the isolated galaxies of the AMIGA sample. We find a progression from field galaxies to filament and cluster ones for decreasing SFR, increasing fraction of galaxies in the quenching phase, increasing proportion of early-type galaxies and decreasing gas content. Galaxies in the quenching phase, defined as having SFR below 1/3 of the main sequence rate, are only between 0-20\% in the isolated sample, while they are 20-60\% in the filaments and 30-80\% in the Virgo cluster. Processes that lead to star formation quenching are already at play in filaments. They depend mostly on the local galaxy density, while the distance to the filament spine is a secondary parameter. While the HI to stellar mass ratio decreases with local density by ~1 dex in the filaments, and ~2 dex in Virgo with respect to the field, the decrease is much less for the H$_2$ to stellar mass ratio. As the environmental density increases, the gas depletion time decreases, since the gas content decreases faster than the SFR. This suggests that gas depletion significantly precedes star formation quenching.

preprint2020arXiv

A molecular absorption line survey toward the AGN of Hydra-A

We present Atacama Large Millimeter/submillimeter Array observations of the brightest cluster galaxy Hydra-A, a nearby ($z=0.054$) giant elliptical galaxy with powerful and extended radio jets. The observations reveal CO(1-0), CO(2-1), $^{13}$CO(2-1), CN(2-1), SiO(5-4), HCO$^{+}$(1-0), HCO$^{+}$(2-1), HCN(1-0), HCN(2-1), HNC(1-0) and H$_{2}$CO(3-2) absorption lines against the galaxy's bright and compact active galactic nucleus. These absorption features are due to at least 12 individual molecular clouds which lie close to the centre of the galaxy and have velocities of approximately $-50$ to $+10$ km/s relative to its recession velocity, where positive values correspond to inward motion. The absorption profiles are evidence of a clumpy interstellar medium within brightest cluster galaxies composed of clouds with similar column densities, velocity dispersions and excitation temperatures to those found at radii of several kpc in the Milky Way. We also show potential variation in a $\sim 10$ km/s wide section of the absorption profile over a two year timescale, most likely caused by relativistic motions in the hot spots of the continuum source which change the background illumination of the absorbing clouds.

preprint2020arXiv

A Virgo Environmental Survey Tracing Ionised Gas Emission. VESTIGE VIII. Bridging the cluster-ICM-galaxy evolution at small scales

We measure FIR emission from tails of stripped dust following the ionised and atomic gas components in galaxies undergoing ram pressure stripping. We study the dust-to-gas relative distribution and mass ratio in the stripped interstellar medium and relate them to those of the intra-cluster medium, thus linking the cluster-ICM-galaxy evolution at small-scales. The galaxy sample consists of three Scd Virgo galaxies with stellar masses in the range $10^9\lesssim \mathrm{M_{*}} \lesssim 10^{10}\, \mathrm{M_{\odot}}$, and within 1 Mpc from the cluster centre, namely NGC 4330, NGC 4522, and NGC 4654. Through the analysis of VESTIGE H$α$, $Herschel$ SPIRE far-infrared, and VIVA HI data, we trace the spatial distribution of the tails and infer the dust and gas masses from the measured far-infrared 250 $μ$m and HI flux densities. Dust-to-gas mass ratios (DGRs) in the tails are analysed as a function of the galaxy mass, metallicity, and dust temperature. Along the stripped component, the dust distribution closely follows the HI and H$α$ emitting gas, all extending beyond the optical disc. In these regions, the DGRs are $2.0\pm0.6\times10^{-3}$, $0.7\pm0.1\times10^{-3}$, and $0.4\pm0.03\times10^{-3}$, for NGC 4330, NGC 4522, and NGC 4654, respectively, i.e. up to a factor of 15 less than the values measured in the main body of nearby galaxies. We also find a negative trend in the DGR as a function of the metallicity that can be explained in terms of a dust component more centrally concentrated in more metal-rich systems. Together with the finding that the stripped dust is cold, $T_{d} \lesssim 25\, K$, our results support an outside-in stripping scenario of the galaxy interstellar medium. This study shows that ram pressure stripping is a key mechanism in the building up of the Virgo intra-cluster component injecting dust grains into the ICM, thus contributing to its metal enrichment.

preprint2020arXiv

ALMA observations of CS in NGC 1068: chemistry and excitation

We present results from Atacama Large Millimeter/submillimeter Array (ALMA) observations of CS from the nearby galaxy NGC 1068 ($\sim14$ Mpc). This Seyfert 2 barred galaxy possesses a circumnuclear disc (CND, $r\sim200$ pc) and a starburst ring (SB ring, $r\sim1.3$ kpc). These high-resolution maps ($\sim0.5$", $\sim35$ pc) allow us to analyse specific sub-regions in the galaxy and investigate differences in line intensity ratios and physical conditions, particularly those between the CND and SB ring. Local thermodynamic equilibrium (LTE) analysis of the gas is used to calculate CS densities in each sub-region, followed by non-LTE analysis conducted using the radiative transfer code RADEX to fit observations and constrain gas temperature, CS column density and hydrogen density. Finally, the chemical code UCLCHEM is used to reconstruct the gas, allowing an insight into its origin and chemical history. The density of hydrogen in the CND is found to be $\geq10^5$ cm$^{-2}$, although exact values vary, reaching $10^6$ cm$^{-2}$ at the AGN. The conditions in the two arms of the SB ring appear similar to one another, though the density found ($\sim10^4$ cm$^{-2}$) is lower than in the CND. The temperature in the CND increases from east to west, and is also overall greater than found in the SB ring. These modelling methods indicate the requirement for multi-phase gas components in order to fit the observed emission over the galaxy. A larger number of high resolution transitions across the SLED may allow for further constraining of the conditions, particularly in the SB ring.

preprint2020arXiv

ALMA resolves the remarkable molecular jet and rotating wind in the extremely radio-quiet galaxy NGC 1377

Submillimetre and millimetre observations are important in probing the properties of the molecular gas and dust around obscured active galactic nuclei (AGNs) and their feedback. With very high-resolution (0."02x0."03 (2x3 pc)) ALMA 345 GHz observations of CO 3-2, HCO$^+$ 4-3, HCN 4-3 $ν_2$=1$f$, and continuum we have studied the molecular outflow and nucleus of the extremely radio-quiet lenticular galaxy NGC1377. The outflow is resolved, revealing a 150 pc long, clumpy, high-velocity, collimated molecular jet. The molecular emission is emerging from the spine of the jet with an average diameter of 3-7 pc. A narrow-angle, rotating molecular wind surrounds the jet and is enveloped by a larger-scale, slower CO-emitting structure. The jet and narrow wind are turbulent ($σ>$40 kms$^{-1}$) and have steep radial gas excitation gradients. The jet shows velocity reversals that we propose are caused by precession, or episodic directional changes. We suggest that an important process powering the outflow is magneto-centrifugal driving. In contrast, the large-scale CO-envelope may be a slow wind, or cocoon that stems from jet-wind interactions. An asymmetric, nuclear r$\sim$2 pc and hot (>180 K) dust structure with a high molecular column density, N(H$_2$)$\sim1.8 \times 10^{24}$ cm$^{-2}$, is detected in continuum and vibrationally excited HCN. Its luminosity is likely powered by a buried AGN. The mass of the supermassive black hole (SMBH) is estimated to $\sim9\times10^6$ M$_\odot$ and the SMBH of NGC1377 appears to be at the end of an intense phase of accretion. The nuclear growth may be fuelled by low-angular momentum gas inflowing from gas ejected in the molecular jet and wind. Such a feedback-loop of cyclic accretion and outflows would be an effective process in growing the nuclear SMBH. This result invites new questions as to SMBH growth processes in obscured, dusty galaxies.

preprint2020arXiv

Blind HI and OH absorption line search: first results with MALS and uGMRT processed using ARTIP

We present details of the Automated Radio Telescope Imaging Pipeline (ARTIP) and results of a sensitive blind search for HI and OH absorbers at $z<0.4$ and $z<0.7$, respectively. ARTIP is written in Python 3.6, extensively uses the Common Astronomy Software Application (CASA) tools and tasks, and is designed to enable the geographically-distributed MeerKAT Absorption Line Survey (MALS) team to collaboratively process large volumes of radio interferometric data. We apply it to the first MALS dataset obtained using the 64-dish MeerKAT radio telescope and 32K channel mode of the correlator. With merely 40 minutes on target, we present the most sensitive spectrum of PKS1830-211 ever obtained and characterize the known HI ($z=0.19$) and OH ($z=0.89$) absorbers. We further demonstrate ARTIP's capabilities to handle realistic observing scenarios by applying it to a sample of 72 bright radio sources observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) to blindly search for HI and OH absorbers. We estimate the numbers of HI and OH absorbers per unit redshift to be $n_{21}(z\sim0.18)<$0.14 and $n_{\rm OH}(z\sim0.40)<$0.12, respectively, and constrain the cold gas covering factor of galaxies at large impact parameters (50 kpc $<ρ<$ 150 kpc) to be less than 0.022. Due to the small redshift path, $Δz\sim$13 for HI with column density$>5.4\times10^{19}$ cm$^{-2}$, the survey has probed only the outskirts of star-forming galaxies at $ρ>30$ kpc. MALS with the expected $Δz\sim10^{3-4}$ will overcome this limitation and provide stringent constraints on the cold gas fraction of galaxies in diverse environments over $0<z<1.5$.

preprint2020arXiv

Cold molecular gas and PAH emission in Seyfert galaxies

We investigate the relation between the detection of the $11.3\,μ$m PAH feature in the nuclear ($\sim 24-230\,$pc) regions of 22 nearby Seyfert galaxies and the properties of the cold molecular gas. For the former we use ground-based (0.3-0.6" resolution) mid-infrared (mid-IR) spectroscopy. The cold molecular gas is traced by ALMA and NOEMA high (0.2-1.1") angular resolution observations of the CO(2-1) transition. Galaxies with a nuclear detection of the $11.3\,μ$m PAH feature contain more cold molecular gas (median $1.6\times 10^7\,M_\odot$) and have higher column densities ($N({\rm H}_2) = 2 \times 10^{23}\,{\rm cm}^{-2}$) over the regions sampled by the mid-IR slits than those without a detection. This suggests that molecular gas plays a role in shielding the PAH molecules in the harsh environments of Seyfert nuclei. Choosing the PAH molecule naphthalene as an illustration, we compute its half-life in the nuclear regions of our sample when exposed to 2.5keV hard X-ray photons. We estimate shorter half-lives for naphthalene in nuclei without a $11.3\,μ$m PAH detection than in those with a detection. The Spitzer/IRS PAH ratios on circumnuclear scales ($\sim$ 4" $\sim$ 0.25-1.3kpc) are in between model predictions for neutral and partly ionized PAHs. However, Seyfert galaxies in our sample with the highest nuclear H$_2$ column densities are not generally closer to the neutral PAH tracks. This is because in the majority of our sample galaxies, the CO(2-1) emission in the inner $\sim$ 4" is not centrally peaked and in some galaxies traces circumnuclear sites of strong star formation activity. Spatially resolved observations with the MIRI medium-resolution spectrograph (MRS) on the James Webb Space Telescope will be able to distinguish the effects of an active galactic nucleus (AGN) and star formation on the PAH emission in nearby AGN.

preprint2020arXiv

Detection of deuterated molecules, but not of lithium hydride, in the z=0.89 absorber toward PKS1830-211

Deuterium and lithium are light elements of high cosmological and astrophysical importance. In this work we report the first detection of deuterated molecules and a search for lithium hydride, 7LiH, at redshift z=0.89 in the spiral galaxy intercepting the line of sight to the quasar PKS1830-211. We used ALMA to observe several submillimeter lines of ND, NH2D, and HDO, and their related isotopomers NH2, NH3, and H2^18O, in absorption against the southwest image of the quasar, allowing us to derive XD/XH abundance ratios. The absorption spectra mainly consist of two distinct narrow velocity components for which we find remarkable differences. One velocity component shows XD/XH abundances that is about 10 times larger than the primordial elemental D/H ratio, and no variability of the absorption profile during the time span of our observations. [...] The second component has XD/XH abundances that are 100 times larger than the primordial D/H ratio, a deepening of the absorption by a factor of two within a few months, and a rich chemical composition, with relative enhancements of N2H+, CH3OH, SO2, and complex organic molecules. We therefore speculate that this component is associated with the analog of a Galactic dark cloud, while the first component is likely more diffuse. Our search for the 7LiH (1--0) line was unsuccessful and we derive an upper limit 7LiH/H2 = 4 x 10^-13 (3sigma) in the z=0.89 absorber toward PKS1830-211. Besides, with ALMA archival data, we could not confirm the previous tentative detections of this line in the z=0.68 absorber toward B0218+357; we derive an upper limit 7LiH/H2 = 5 x 10^-11 (3sigma), although this is less constraining than our limit toward PKS1830-211. We conclude that, as in the Milky Way, only a tiny fraction of lithium nuclei are possibly bound in LiH in these absorbers at intermediate redshift.

preprint2020arXiv

Environmental processing in cluster core galaxies at z=1.7

Today, the brightest cluster galaxies (BCGs) are passive and very massive galaxies at the center of their clusters, and they still accrete mass through swallowing companions and gas from cooling flows. However their formation history is not well known. We report CO(4$\rightarrow$3) and continuum map observations of the SpARCS1049+56 BCG at $z=1.709$, one of the most distant known BCGs. Our observations yield $M_{{\rm H}_2}<1.1\times10^{10}M_\odot$ for the BCG; while in CO(4$\rightarrow$3), we detect two gas-rich companions at the northeast and southeast of the BCG, within 20 kpc, with ${L^\prime_{\rm CO(4\rightarrow3)}=(5.8\pm0.6)\times10^{9}}$ K km s$^{-1}$ pc$^2$ and $(7.4\pm0.7)\times10^{9}$ K km s$^{-1}$ pc$^2$, respectively. The northern companion is associated with a pair of merging cluster galaxies, while the southern one shows a southern tail in CO(4$\rightarrow$3), which was also detected in continuum, and we suggest it to be the most distant jellyfish galaxy for which ram pressure stripping is effectively able to strip off its dense molecular gas. This study probes the presence of rare gas-rich systems in the very central region of a distant cluster core, which will potentially merge into the BCG itself. Currently, we may thus be seeing the reversal of the star formation versus density relation at play in the distant universe. This is the first time the assembly of high-$z$ progenitors of our local BCGs can be studied in such great detail.

preprint2020arXiv

Molecular gas and star formation activity in LIRGs in clusters at intermediate redshifts

We investigate the role of dense Mpc-scale environments in processing molecular gas of cluster galaxies as they fall into the cluster cores. We consider $\sim20$ luminous infrared galaxies (LIRGs) in intermediate-$z$ clusters, from the Hershel Lensing Survey and the Local Cluster Substructure Survey. They include MACS J0717.5+3745 at $z=0.546$ and Abell 697, 963, 1763, and 2219 at $z=0.2-0.3$. We have performed far infrared to ultraviolet spectral energy distribution modeling of the LIRGs, which span cluster-centric distances within $r/r_{200}\simeq0.2-1.6$. We have observed the LIRGs in CO(1$\rightarrow$0) or CO(2$\rightarrow$1) with the Plateau de Bure interferometer and its successor NOEMA, as part of five observational programs carried out between 2012 and 2017. We have compared the molecular gas to stellar mass ratio $M(H_2)/M_\star$, star formation rate (SFR), and depletion time ($τ_{\rm dep}$) of the LIRGs with those of a compilation of cluster and field star forming galaxies. The targeted LIRGs have SFR, $M(H_2)/M_\star$, and $τ_{\rm dep}$ that are consistent with those of both main sequence (MS) field galaxies and star forming galaxies from the comparison sample. However we find that the depletion time, normalized to the MS value, increases with increasing $r/r_{200}$, with a significance of $2.8σ$, which is ultimately due to a deficit of cluster core LIRGs with $τ_{\rm dep}\gtrsimτ_{\rm dep,MS}$. We suggest that a rapid exhaustion of the molecular gas reservoirs occurs in the cluster LIRGs and is effective in suppressing their star formation. This mechanism may explain the exponential decrease of the fraction of cluster LIRGs with cosmic time. The compression of the gas in LIRGs, possibly induced by intra-cluster medium shocks, may be responsible for the short depletion timescales, observed in a large fraction of cluster core LIRGs.

preprint2020arXiv

Molecular gas in CLASH brightest cluster galaxies at $z\sim0.2-0.9$

Brightest cluster galaxies (BCGs) are excellent laboratories to study galaxy evolution in dense Mpc-scale environments. We have observed in CO(1-0), CO(2-1), CO(3-2), or CO(4-3), with the IRAM-30m, 18 BCGs at $z\sim0.2-0.9$ that are drawn from the CLASH survey. Our sample includes RX1532, which is our primary target, being among the BCGs with the highest star formation rate (SFR$\gtrsim100~M_\odot$/yr) in the CLASH sample. We unambiguously detected both CO(1-0) and CO(3-2) in RX1532, yielding a large reservoir of molecular gas, $M_{H_2}=(8.7\pm1.1)\times10^{10}~M_\odot$, and a high level of excitation $r_{31}=0.75\pm0.12$. A morphological analysis of the HST I-band image of RX1532 reveals the presence of clumpy substructures both within and outside the half-light radius $r_e=(11.6\pm0.3)$ kpc, similarly to those found independently both in ultraviolet and in H$_α$ in previous work. We tentatively detected CO(1-0) or CO(2-1) in four other BCGs, with molecular gas reservoirs in the range $M_{H_2}=2\times10^{10-11} M_\odot$. For the remaining 13 BCGs we set robust upper limits of $M_{H_2}/M_\star\lesssim0.1$, which are among the lowest molecular gas to stellar mass ratios found for distant ellipticals and BCGs. By comparison with distant cluster galaxies observed in CO our study shows that RX1532 ($M_{H_2}/M_\star = 0.40\pm0.05$) belongs to the rare population of star forming and gas-rich BCGs in the distant universe. By using available X-ray based estimates of the central intra-cluster medium entropy, we show that the detection of large reservoirs of molecular gas $M_{H_2}\gtrsim10^{10}~M_\odot$ in distant BCGs is possible when the two conditions are met: i) high SFR and ii) low central entropy, which favors the condensation and the inflow of gas onto the BCGs themselves, similarly to what has been previously found for some local BCGs.

preprint2020arXiv

Molecular gas in distant brightest cluster galaxies

The stellar mass assembly of brightest cluster galaxies (BCGs) is still debated. We have observed in CO with the IRAM-30m two star forming BCGs belonging to SpARCS clusters, 3C 244.1 ($z=0.4$) and SDSS J161112.65+550823.5 ($z=0.9$), and compared their molecular gas and star formation properties with those of a compilation of $\sim100$ distant cluster galaxies, including nine additional distant BCGs at ${z\sim0.4-3.5}$. We have set robust upper limits of $M_{{\rm H}_2}<1.0\times10^{10}~M_\odot$ and $<2.8\times10^{10}~M_\odot$ to their molecular gas content, respectively, as well as to the molecular gas to stellar mass ratio $M({\rm H}_2)/M_\star\lesssim0.2$ and depletion time $τ_{\rm dep}\lesssim40$ Myr of the two targeted BCGs. They are thus among the distant cluster galaxies with the lowest gas fractions and shortest depletion times. The majority, $64\%\pm15\%$ and $73\%\pm18\%$, of the 11 BCGs with observations in CO have lower $M({\rm H}_2)/M_\star$ and $τ_{\rm dep}$, respectively, than those estimated for main sequence galaxies. Statistical analysis tentatively suggests that the values of $M({\rm H}_2)/M_\star$ and $τ_{\rm dep}$ for the 11 BCGs deviates, with a significance of $\sim2σ$, from those of the comparison sample of cluster galaxies. A morphological analysis for a subsample of seven BCGs with archival HST observations reveals that $71\%\pm17\%$ of the BCGs are compact or show star-forming substructures/components. Our results suggest a scenario where distant star forming BCGs assemble a significant fraction $\sim16\%$ of their stellar mass in the relatively short timescale $τ_{\rm dep}$, while environmental mechanisms might prevent the replenishment of gas feeding the star formation. We speculate that compact components favor the rapid exhaustion of molecular gas and help to quench the BCGs. Distant star forming BCGs are excellent targets for ALMA and JWST.

preprint2019arXiv

A CO molecular gas wind 340 pc away from the Seyfert 2 nucleus in ESO 420-G13 probes an elusive radio jet

A prominent jet-driven outflow of CO(2-1) molecular gas is found along the kinematic minor axis of the Seyfert 2 galaxy ESO 420-G13, at a distance of $340-600\, \rm{pc}$ from the nucleus. The wind morphology resembles a characteristic funnel shape, formed by a highly collimated filamentary emission at the base, likely tracing the jet propagation through a tenuous medium, until a bifurcation point at $440\, \rm{pc}$ where the jet hits a dense molecular core and shatters, dispersing the molecular gas into several clumps and filaments within the expansion cone. We also trace the jet in ionised gas within the inner $\lesssim 340\, \rm{pc}$ using the [NeII]$_{\rm 12.8 μm}$ line emission, where the molecular gas follows a circular rotation pattern. The wind outflow carries a mass of $\sim 8 \times 10^6\, \rm{M_\odot}$ at an average wind projected speed of $\sim 160\, \rm{km\,s^{-1}}$, which implies a mass outflow rate of $\sim 14\, \rm{M_\odot\,yr^{-1}}$. Based on the structure of the outflow and the budget of energy and momentum, we discard radiation pressure from the active nucleus, star formation, and supernovae as possible launching mechanisms. ESO 420-G13 is the second case after NGC 1377 where the presence of a previously unknown jet is revealed due to its interaction with the interstellar medium, suggesting that unknown jets in feeble radio nuclei might be more common than expected. Two possible jet-cloud configurations are discussed to explain the presence of an outflow at such distance from the AGN. The outflowing gas will likely not escape, thus a delay in the star formation rather than quenching is expected from this interaction, while the feedback effect would be confined within the central few hundred parsecs of the galaxy.

preprint2016arXiv

A new look at the kinematics of the bulge from an N-body model

(Abridged) By using an N-body simulation of a bulge that was formed via a bar instability mechanism, we analyse the imprints of the initial (i.e. before bar formation) location of stars on the bulge kinematics, in particular on the heliocentric radial velocity distribution of bulge stars. Four different latitudes were considered: $b=-4^\circ$, $-6^\circ$, $-8^\circ$, and $-10^\circ$, along the bulge minor axis as well as outside it, at $l=\pm5^\circ$ and $l=\pm10^\circ$. The bulge X-shaped structure comprises stars that formed in the disk at different locations. Stars formed in the outer disk, beyond the end of the bar, which are part of the boxy peanut-bulge structure may show peaks in the velocity distributions at positive and negative heliocentric radial velocities with high absolute values that can be larger than 100 $\rm km$ $\rm s^{-1}$, depending on the observed direction. In some cases the structure of the velocity field is more complex and several peaks are observed. Stars formed in the inner disk, the most numerous, contribute predominantly to the X-shaped structure and present different kinematic characteristics. Our results may enable us to interpret the cold high-velocity peak observed in the APOGEE commissioning data, as well as the excess of high-velocity stars in the near and far arms of the X-shaped structure at $l$=$0^\circ$ and $b$=$-6^\circ$. When compared with real data, the kinematic picture becomes more complex due to the possible presence in the observed samples of classical bulge and/or thick disk stars. Overall, our results point to the existence of complex patterns and structures in the bulge velocity fields, which are generated by the bar. This suggests that caution should be used when interpreting the bulge kinematics: the presence of substructures, peaks and clumps in the velocity fields is not necessarily a sign of past accretion events.

preprint2016arXiv

A precessing molecular jet signaling an obscured, growing supermassive black hole in NGC1377?

With high resolution (0."25 x 0."18) ALMA CO 3-2 observations of the nearby (D=21 Mpc, 1"=102 pc), extremely radio-quiet galaxy NGC1377, we have discovered a high-velocity, very collimated nuclear outflow which we interpret as a molecular jet with a projected length of +-150 pc. Along the jet axis we find strong velocity reversals where the projected velocity swings from -150 km/s to +150 km/s. A simple model of a molecular jet precessing around an axis close to the plane of the sky can reproduce the observations. The velocity of the outflowing gas is difficult to constrain due to the velocity reversals but we estimate it to be between 240 and 850 km/s and the jet to precess with a period P=0.3-1.1 Myr. The CO emission is clumpy along the jet and the total molecular mass in the high-velocity (+-(60 to 150 km/s)) gas lies between 2e6 Msun (light jet) and 2e7 Msun (massive jet). There is also CO emission extending along the minor axis of NGC1377. It holds >40% of the flux in NGC1377 and may be a slower, wide-angle molecular outflow which is partially entrained by the molecular jet. We discuss the driving mechanism of the molecular jet and suggest that it is either powered by a very faint radio jet or by an accretion disk-wind similar to those found towards protostars. The nucleus of NGC1377 harbours intense embedded activity and we detect emission from vibrationally excited HCN J=4-3 v_2=1f which is consistent with hot gas and dust. We find large columns of H2 in the centre of NGC1377 which may be a sign of a high rate of recent gas infall. The dynamical age of the molecular jet is short (<1 Myr), which could imply that it is young and consistent with the notion that NGC1377 is caught in a transient phase of its evolution. However, further studies are required to determine the age of the molecular jet, its mass and the role it is playing in the growth of the nucleus of NGC1377.

preprint2016arXiv

A search for HI and OH absorption in z > 3 CO emitters

We present the results of a survey for HI 21-cm and OH 18-cm absorption in seven strong CO emitters at z > 3. Despite reaching limits comparable to those required to detect 21-cm absorption at lower redshifts, we do not detect either transition in any of the objects searched. We believe that this is due to the high redshift selection causing all of our targets to have ultra-violet luminosities above the critical value, where all of the atomic gas in the host galaxy disk is suspected to be ionised. However, not only are all of our targets bright in CO emission, but detection of CO above the critical UV luminosity is generally not uncommon. This suggests that the molecular gas is shielded from the radiation or is physically remote from the source of the continuum emission, as it appears to be from CO observations of high redshift radio galaxies.

preprint2016arXiv

ALMA observations of cold molecular gas filaments trailing rising radio bubbles in PKS0745-191

We present ALMA observations of the CO(1-0) and CO(3-2) line emission tracing filaments of cold molecular gas in the central galaxy of the cluster PKS0745-191. The total molecular gas mass of 4.6 +/- 0.3 x 10^9 solar masses, assuming a Galactic X_{CO} factor, is divided roughly equally between three filaments each extending radially 3-5 kpc from the galaxy centre. The emission peak is located in the SE filament roughly 1 arcsec (2 kpc) from the nucleus. The velocities of the molecular clouds in the filaments are low, lying within +/-100 km/s of the galaxy's systemic velocity. Their FWHMs are less than 150 km/s, which is significantly below the stellar velocity dispersion. Although the molecular mass of each filament is comparable to a rich spiral galaxy, such low velocities show that the filaments are transient and the clouds would disperse on <10^7 yr timescales unless supported, likely by the indirect effect of magnetic fields. The velocity structure is inconsistent with a merger origin or gravitational free-fall of cooling gas in this massive central galaxy. If the molecular clouds originated in gas cooling even a few kpc from their current locations their velocities would exceed those observed. Instead, the projection of the N and SE filaments underneath X-ray cavities suggests they formed in the updraft behind bubbles buoyantly rising through the cluster atmosphere. Direct uplift of the dense gas by the radio bubbles appears to require an implausibly high coupling efficiency. The filaments are coincident with low temperature X-ray gas, bright optical line emission and dust lanes indicating that the molecular gas could have formed from lifted warmer gas that cooled in situ.

preprint2016arXiv

ALMA resolves the torus of NGC 1068: continuum and molecular line emission

We have used the Atacama Large Millimeter Array (ALMA) to map the emission of the CO(6-5) molecular line and the 432 μm continuum emission from the 300 pc-sized circumnuclear disk (CND) of the nearby Seyfert 2 galaxy NGC 1068 with a spatial resolution of ~4 pc. These observations spatially resolve the CND and, for the first time, image the dust emission, the molecular gas distribution, and the kinematics from a 7-10 pc-diameter disk that represents the submillimeter counterpart of the putative torus of NGC 1068. We fitted the nuclear spectral energy distribution of the torus using ALMA and near and mid-infrared (NIR/MIR) data with CLUMPY models. The mass and radius of the best-fit solution for the torus are both consistent with the values derived from the ALMA data alone: Mgas_torus=(1+-0.3)x10^5 Msun and Rtorus=3.5+-0.5 pc. The dynamics of the molecular gas in the torus show non-circular motions and enhanced turbulence superposed on the rotating pattern of the disk. The kinematic major axis of the CO torus is tilted relative to its morphological major axis. By contrast with the nearly edge-on orientation of the H2O megamaser disk, we have found evidence suggesting that the molecular torus is less inclined (i=34deg-66deg) at larger radii. The lopsided morphology and complex kinematics of the torus could be the signature of the Papaloizou-Pringle instability, long predicted to likely drive the dynamical evolution of active galactic nuclei (AGN) tori.

preprint2016arXiv

ALMA reveals optically thin, highly excited CO gas in the jet-driven winds of the galaxy IC5063

Using CO (4-3) and (2-1) Atacama Large Millimeter Array (ALMA) data, we prove that the molecular gas in the jet-driven winds of the galaxy IC5063 is more highly excited than the rest of the molecular gas in the disk of the same galaxy. On average, the CO(4-3)/CO(2-1) flux ratio is 1 for the disk and 5 for the jet accelerated or impacted gas. Spatially-resolved maps reveal that in regions associated with winds, the CO(4-3)/CO(2-1) flux ratio significantly exceeds the upper limit of 4 for optically thick gas. It frequently takes values between 5 and 11, and it occasionally further approaches the upper limit of 16 for optically thin gas. Excitation temperatures of 30-100 K are common for the molecules in these regions. If all of the outflowing molecular gas is optically thin, at 30-50 K, then its mass is 2*10^6 M_sun. This lower mass limit is an order of magnitude below the mass derived from the CO(2-1) flux in the case of optically thick gas. Molecular winds can thus be less massive, but more easily detectable at high z than they were previously thought to be.

preprint2016arXiv

Exploring the GalMer database: bar properties and non-circular motions

We use Tree-SPH simulations from the GalMer database by Chilingarian et al. to characterize and quantify the non-circular motions induced by the presence of bar-like structures on the observed rotation curve of barred galaxies derived from empirical models of their line-of-sight velocity maps. The GalMer database consists of SPH simulations of galaxies spanning a wide range of morphological types and sizes. The aim is to compare the intrinsic velocities and bar properties from the simulations with those derived from pseudo-observations. This allows us to estimate the amount of non-circularity and to test the various methods used to derive the bar properties and rotation curves. The intrinsic velocities in the simulations are calculated from the gravitational forces whereas the observed rotation velocities are derived by applying the ROTCUR and DiskFit algorithms to well-resolved observations of intermediate-inclination, strongly barred galaxies. Our results confirm that the tilted ring method implemented in ROTCUR systematically underestimates/overestimates the rotational velocities by up to 40 percent in the inner part of the galaxy when the bar is aligned with one of the symmetry axes for all the models. For the DiskFit analysis, we find that it produces unrealistic values for all the models used in this work when the bar is within $\sim$10 degrees of the major or minor axis.

preprint2016arXiv

High-redshift major mergers weakly enhance star formation

Galaxy mergers are believed to trigger strong starbursts. This is well assessed by observations in the local Universe. However the efficiency of this mechanism has poorly been tested so far for high redshift, actively star forming, galaxies. We present a suite of pc-resolution hydrodynamical numerical simulations to compare the star formation process along a merging sequence of high and low z galaxies, by varying the gas mass fraction between the two models. We show that, for the same orbit, high-redshift gas-rich mergers are less efficient than low-redshift ones at producing starbursts: the star formation rate excess induced by the merger and its duration are both around 10 times lower than in the low gas fraction case. The mechanisms that account for the star formation triggering at low redshift - the increased compressive turbulence, gas fragmentation, and central gas inflows - are only mildly, if not at all, enhanced for high gas fraction galaxy encounters. Furthermore, we show that the strong stellar feedback from the initially high star formation rate in high redshift galaxies does not prevent an increase of the star formation during the merger. Our results are consistent with the observed increase of the number of major mergers with increasing redshift being faster than the respective increase in the number of starburst galaxies.

preprint2016arXiv

HST imaging of the dusty filaments and nucleus swirl in NGC4696 at the centre of the Centaurus Cluster

Narrow-band HST imaging has resolved the detailed internal structure of the 10 kpc diameter H alpha+[NII] emission line nebulosity in NGC4696, the central galaxy in the nearby Centaurus cluster, showing that the dusty, molecular, filaments have a width of about 60pc. Optical morphology and velocity measurements indicate that the filaments are dragged out by the bubbling action of the radio source as part of the AGN feedback cycle. Using the drag force we find that the magnetic field in the filaments is in approximate pressure equipartition with the hot gas. The filamentary nature of the cold gas continues inward, swirling around and within the Bondi accretion radius of the central black hole, revealing the magnetic nature of the gas flows in massive elliptical galaxies. HST imaging resolves the magnetic, dusty, molecular filaments at the centre of the Centaurus cluster to a swirl around and within the Bondi radius.

preprint2016arXiv

Molecular Gas Along a Bright H-alpha Filament in 2A 0335+096 Revealed by ALMA

We present ALMA CO(1-0) and CO(3-2) observations of the brightest cluster galaxy (BCG) in the 2A 0335+096 galaxy cluster (z = 0.0346). The total molecular gas mass of (1.13+/-0.15) x 10^9 M_sun is divided into two components: a nuclear region and a 7 kpc long dusty filament. The central molecular gas component accounts for (3.2+/-0.4) x 10^8 M_sun of the total supply of cold gas. Instead of forming a rotationally-supported ring or disk, it is composed of two distinct, blueshifted clumps south of the nucleus and a series of low-significance redshifted clumps extending toward a nearby companion galaxy. The velocity of the redshifted clouds increases with radius to a value consistent with the companion galaxy, suggesting that an interaction between these galaxies <20 Myr ago disrupted a pre-existing molecular gas reservoir within the BCG. Most of the molecular gas, (7.8+/-0.9) x 10^8 M_sun, is located in the filament. The CO emission is co-spatial with a 10^4 K emission-line nebula and soft X-rays from 0.5 keV gas, indicating that the molecular gas has cooled out of the intracluster medium over a period of 25-100 Myr. The filament trails an X-ray cavity, suggesting that the gas has cooled from low entropy gas that has been lifted out of the cluster core and become thermally unstable. We are unable to distinguish between inflow and outflow along the filament with the present data. Cloud velocities along the filament are consistent with gravitational free-fall near the plane of the sky, although their increasing blueshifts with radius are consistent with outflow.

preprint2016arXiv

Non-thermal emission and dynamical state of massive galaxy clusters from CLASH sample

Massive galaxy clusters are the most violent large scale structures undergoing merger events in the Universe. Based upon their morphological properties in X-rays, they are classified as un-relaxed and relaxed clusters and often host (a fraction of them) different types of non-thermal radio emitting components, viz., haloes, mini-haloes, relics and phoenix within their Intra Cluster Medium (ICM). The radio haloes show steep (alpha = -1.2) and ultra steep (alpha < -1.5) spectral properties at low radio frequencies, giving important insights on the merger (pre or post) state of the cluster. Ultra steep spectrum radio halo emissions are rare and expected to be the dominating population to be discovered via LOFAR and SKA in the future. Further, the distribution of matter (morphological information), alignment of hot X-ray emitting gas from the ICM with the total mass (dark + baryonic matter) and the bright cluster galaxy (BCG) is generally used to study the dynamical state of the cluster. We present here a multi wavelength study on 14 massive clusters from the CLASH survey and show the correlation between the state of their merger in X-ray and spectral properties (1.4 GHz - 150 MHz) at radio wavelengths. Using the optical data we also discuss about the gas-mass alignment, in order to understand the interplay between dark and baryonic matter in massive galaxy clusters.

preprint2016arXiv

OH+ and H2O+ absorption toward PKS1830-211

We report the detection of OH+ and H2O+ in the z=0.89 absorber toward the lensed quasar PKS1830-211. The abundance ratio of OH+ and H2O+ is used to quantify the molecular hydrogen fraction (fH2) and the cosmic-ray ionization rate of atomic hydrogen (zH) along two lines of sight, located at ~2 kpc and ~4 kpc to either side of the absorber's center. The molecular fraction decreases outwards, from ~0.04 to ~0.02, comparable to values measured in the Milky Way at similar galactocentric radii. For zH, we find values of ~2x10^-14 s^-1 and ~3x10^-15 s^-1, respectively, which are slightly higher than in the Milky Way at comparable galactocentric radii, possibly due to a higher average star formation activity in the z=0.89 absorber. The ALMA observations of OH+, H2O+, and other hydrides toward PKS1830-211 reveal the multi-phase composition of the absorbing gas. Taking the column density ratios along the southwest and northeast lines of sight as a proxy of molecular fraction, we classify the species ArH+, OH+, H2Cl+, H2O+, CH, and HF as tracing gases increasingly more molecular. Incidentally, our data allow us to improve the accuracy of H2O+ rest frequencies and thus refine the spectroscopic parameters.

preprint2016arXiv

Optical Emission Line Nebulae in Galaxy Cluster Cores 1: The Morphological, Kinematic and Spectral Properties of the Sample

We present an Integral Field Unit survey of 73 galaxy clusters and groups with the VIsible Multi Object Spectrograph (VIMOS) on VLT. We exploit the data to determine the H$α$ gas dynamics on kpc-scales to study the feedback processes occurring within the dense cluster cores. We determine the kinematic state of the ionised gas and show that the majority of systems ($\sim$ 2/3) have relatively ordered velocity fields on kpc scales that are similar to the kinematics of rotating discs and are decoupled from the stellar kinematics of the Brightest Cluster Galaxy. The majority of the H$α$ flux ($>$ 50%) is typically associated with these ordered kinematics and most systems show relatively simple morphologies suggesting they have not been disturbed by a recent merger or interaction. Approximately 20% of the sample (13/73) have disturbed morphologies which can typically be attributed to AGN activity disrupting the gas. Only one system shows any evidence of an interaction with another cluster member. A spectral analysis of the gas suggests that the ionisation of the gas within cluster cores is dominated by non stellar processes, possibly originating from the intracluster medium itself.

preprint2016arXiv

The Close AGN Reference Survey (CARS) - Mrk 1018 returns to the shadows after 30 years as a Seyfert 1

We report the discovery that the known `changing look' AGN Mrk 1018 has changed spectral type for a second time. New VLT-MUSE data taken in 2015 as part of the Close AGN Reference Survey (CARS) shows that the AGN has returned to its original Seyfert 1.9 classification. The CARS sample is selected to contain only bright type 1 AGN, but Mrk 1018's broad emission lines and continuum, typical of type 1 AGN, have almost entirely disappeared. We use spectral fitting of the MUSE spectrum and previously available spectra to determine the drop in broad line flux and the Balmer decrement. We find that the broad line flux has decreased by a factor of 4.75+/- 0.5 in Hα since an SDSS spectrum was taken in 2000. The Balmer decrement has not changed significantly implying no enhanced reddening with time, but the remaining broad lines are more asymmetric than those present in the type 1 phase. We posit that the change is due to an intrinsic drop in flux from the accretion disk rather than variable extinction or a tidal disruption event.

preprint2016arXiv

The Close AGN Reference Survey (CARS) - What is causing Mrk1018's return to the shadows after 30 years?

We recently discovered that the active galactic nucleus (AGN) of Mrk 1018 has changed optical type again after 30 years as a type 1 AGN. Here we combine Chandra, NuStar, Swift, Hubble Space Telescope and ground-based observations to explore the cause of this change. The 2-10keV flux declines by a factor of ~8 between 2010 and 2016. We show with our X-ray observation that this is not caused by varying neutral hydrogen absorption along the line-of-sight up to the Compton-thick level. The optical-UV spectral energy distributions are well fit with a standard geometrically thin optically thick accretion disc model that seems to obey the expected $L\sim T^4$ relation. It confirms that a decline in accretion disc luminosity is the primary origin for the type change. We detect a new narrow-line absorber in Lya blue-shifted by ~700km/s with respect to the systemic velocity of the galaxy. This new Lya absorber could be evidence for the onset of an outflow or a companion black hole with associated gas that could be related to the accretion rate change. However, the low column density of the absorber means that it is not the direct cause for Mrk 1018's changing-look nature.

preprint2016arXiv

The low-frequency radio emission in blazar PKS2155-304

We report radio imaging and monitoring observations in the frequency range 0.235 - 2.7 GHz during the flaring mode of PKS 2155-304, one of the brightest BL Lac objects. The high sensitivity GMRT observations not only reveal extended kpc-scale jet and FRI type lobe morphology in this erstwhile `extended-core' blazar but also delineate the morphological details, thanks to its arcsec scale resolution. The radio light curve during the end phase of the outburst measured in 2008 shows high variability (8.5%) in the jet emission in the GHz range, compared to the lower core variability (3.2%) seen at the lowest frequencies. The excess of flux density with a very steep spectral index in the MHz range supports the presence of extra diffuse emission at low frequencies. The analysis of multi wavelength (radio/ optical/ gamma-ray) light curves at different radio frequencies confirms the variability of the core region and agrees with the scenario of high energy emission in gamma-rays due to inverse Compton emission from a collimated relativistic plasma jet followed by synchrotron emission in radio. Clearly, these results give an interesting insight of the jet emission mechanisms in blazars and highlight the importance of studying such objects with low frequency radio interferometers like LOFAR and the SKA and its precursor instruments.

preprint2016arXiv

The overmassive black hole in NGC 1277: new constraints from molecular gas kinematics

We report the detection of CO(1-0) emission from NGC 1277, a lenticular galaxy in the Perseus Cluster. NGC 1277 has previously been proposed to host an overmassive black hole (BH) compared to the galaxy bulge luminosity (mass), based on stellar-kinematic measurements. The CO(1-0) emission, observed with the IRAM Plateau de Bure Interferometer (PdBI) using both, a more compact (2.9-arcsec resolution) and a more extended (1-arcsec resolution) configuration, is likely to originate from the dust lane encompassing the galaxy nucleus at a distance of 0.9 arcsec (~320 pc). The double-horned CO(1-0) profile found at 2.9-arcsec resolution traces $1.5\times 10^8\ M_\odot$ of molecular gas, likely orbiting in the dust lane at $\sim 550\ \mathrm{km\ s^{-1}}$, which suggests a total enclosed mass of $\sim 2\times 10^{10}\ M_\odot$. At 1-arcsec resolution, the CO(1-0) emission appears spatially resolved along the dust lane in east-west direction, though at a low signal-to-noise ratio. In agreement with the previous stellar-kinematic measurements, the CO(1-0) kinematics is found to be consistent with an $\sim 1.7\times 10^{10}\ M_\odot$ BH for a stellar mass-to-light ratio of $M/L_V=6.3$, while a less massive BH of $\sim 5\times 10^{9}\ M_\odot$ is possible when assuming a larger $M/L_V=10$. While the molecular gas reservoir may be associated with a low level of star formation activity, the extended 2.6-mm continuum emission is likely to originate from a weak AGN, possibly characterized by an inverted radio-to-millimetre spectral energy distribution. Literature radio and X-ray data indicate that the BH in NGC 1277 is also overmassive with respect to the Fundamental Plane of BH activity.

preprint2015arXiv

A radio jet drives a molecular and atomic gas outflow in multiple regions within one square kiloparsec of the nucleus of the nearby galaxy IC5063

We analyzed near-infrared data of the nearby galaxy IC5063 taken with the Very Large Telescope SINFONI instrument. IC5063 is an elliptical galaxy that has a radio jet nearly aligned with the major axis of a gas disk in its center. The data reveal multiple signatures of molecular and atomic gas that has been kinematically distorted by the passage of the jet plasma or cocoon within an area of ~1 kpc^2. Concrete evidence that the interaction of the jet with the gas causes the gas to accelerate comes from the detection of outflows in four different regions along the jet trail: near the two radio lobes, between the radio emission tip and the optical narrow-line-region cone, and at a region with diffuse 17.8 GHz emission midway between the nucleus and the north radio lobe. The outflow in the latter region is biconical, centered 240 pc away from the nucleus, and oriented perpendicularly to the jet trail. The diffuse emission that is observed as a result of the gas entrainment or scattering unfolds around the trail and away from the nucleus with increasing velocity. It overall extends for >700 pc parallel and perpendicular to the trail. Near the outflow starting points, the gas has a velocity excess of 600 km/s to 1200 km/s with respect to ordered motions, as seen in [FeII], Pa alpha, or H2 lines. High H2 (1-0) S(3)/S(1) flux ratios indicate non-thermal excitation of gas in the diffuse outflow.

preprint2015arXiv

ALMA detection of [CII] 158 micron emission from a strongly lensed z=2 star-forming galaxy

Our objectives are to determine the properties of the interstellar medium (ISM) and of star-formation in typical star-forming galaxies at high redshift. Following up on our previous multi-wavelength observations with HST, Spitzer, Herschel, and the Plateau de Bure Interferometer (PdBI), we have studied a strongly lensed z=2.013 galaxy, the arc behind the galaxy cluster MACS J0451+0006, with ALMA to measure the [CII] 158 micron emission line, one of the main coolants of the ISM. [CII] emission from the southern part of this galaxy is detected at 10 $σ$. Taking into account strong gravitational lensing, which provides a magnification of $μ=49$, the intrinsic lensing-corrected [CII]158 micron luminosity is $L(CII)=1.2 \times 10^8 L_\odot$. The observed ratio of [CII]-to-IR emission, $L(CII)/L(FIR) \approx (1.2-2.4) \times 10^{-3}$, is found to be similar to that in nearby galaxies. The same also holds for the observed ratio $L(CII)/L(CO)=2.3 \times 10^3$, which is comparable to that of star-forming galaxies and active galaxy nuclei (AGN) at low redshift. We utilize strong gravitational lensing to extend diagnostic studies of the cold ISM to an order of magnitude lower luminosity ($L(IR) \sim (1.1-1.3) \times 10^{11} L_\odot$) and SFR than previous work at high redshift. While larger samples are needed, our results provide evidence that the cold ISM of typical high redshift galaxies has physical characteristics similar to normal star forming galaxies in the local Universe.

preprint2015arXiv

Dense gas tracing the collisional past of Andromeda. An atypical inner region?

The central kiloparsec region of the Andromeda galaxy is relatively gas poor, while the interstellar medium appears to be concentrated in a ring-like structure at about 10 kpc radius. The central gas depletion has been attributed to a possible head-on collision 200 Myr ago, supported by the existence of an offset inner ring of warm dust. We present new IRAM 30m radio telescope observations of the molecular gas in the central region, and the detection of CO and its isotopes $^{13}$CO(2-1) and C$^{18}$O(2-1), together with the dense gas tracers, HCN(1-0) and HCO+(1-0). A systematic study of the observed peak temperatures with non-LTE equilibrium simulations shows that the detected lines trace dense regions with n$_{H_2}$ in the range 2.5 $10^4 - 5.6 10^5$ cm$^{-3}$, while the gas is very clumpy with a beam filling factor of 0.5-2 10$^{-2}$. This is compatible with the dust mass derived from the far-infrared emission, assuming a dust-to-gas mass ratio of 0.01 with a typical clump size of 2 pc. We also show that the gas is optically thin in all lines except for $^{12}$CO(1-0) and $^{12}$CO(2-1), CO lines are close to their thermal equilibrium condition at 17.5-20 K, the molecular hydrogen density is larger than critical and HCN and HCO+ lines have a subthermal excitation temperature of 9 K with a density smaller than critical. The average $^{12}$CO/$^{13}$CO line ratio is high (~21), and close to the $^{12}$CO/C$^{18}$O ratio (~30) that was measured in the north-western region and estimated in the south-east stacking. The fact that the optically thin $^{13}$CO and C$^{18}$O lines have comparable intensities means that the secondary element $^{13}$C is depleted with respect to the primary $^{12}$C, as is expected just after a recent star formation. This suggests that there has been a recent starburst in the central region, supporting the head-on collision scenario.

preprint2015arXiv

Environmental regulation of cloud and star formation in galactic bars

The strong time-dependence of the dynamics of galactic bars yields a complex and rapidly evolving distribution of dense gas and star forming regions. Although bars mainly host regions void of any star formation activity, their extremities can gather the physical conditions for the formation of molecular complexes and mini-starbursts. Using a sub-parsec resolution hydrodynamical simulation of a Milky Way-like galaxy, we probe these conditions to explore how and where bar (hydro-)dynamics favours the formation or destruction of molecular clouds and stars. The interplay between the kpc-scale dynamics (gas flows, shear) and the parsec-scale (turbulence) is key to this problem. We find a strong dichotomy between the leading and trailing sides of the bar, in term of cloud fragmentation and in the age distribution of the young stars. After orbiting along the bar edge, these young structures slow down at the extremities of the bar, where orbital crowding increases the probability of cloud-cloud collision. We find that such events increase the Mach number of the cloud, leading to an enhanced star formation efficiency and finally the formation of massive stellar associations, in a fashion similar to galaxy-galaxy interactions. We highlight the role of bar dynamics in decoupling young stars from the clouds in which they form, and discuss the implications on the injection of feedback into the interstellar medium, in particular in the context of galaxy formation.

preprint2015arXiv

Explaining the formation of bulges with MOND

In the cold dark matter (CDM) paradigm, bulges easily form through galaxy mergers, either major or minor, or through clumpy disks in the early universe, where clumps are driven to the center by dynamical friction. Also pseudo-bulges, with a more disky morphology and kinematics, can form more slowly through secular evolution of a bar, where resonant stars are elevated out of the plane, in a peanut/box shape. As a result, in CDM cosmological simulations, it is very difficult to find a bulgeless galaxy, while they are observed very frequently in the local universe. A different picture emerges in alternative models of the missing mass problem. In MOND (MOdified Newtonian Dynamics), galaxy mergers are much less frequent, since the absence of dark matter halos reduces the dynamical friction between two galaxies. Also, while clumpy galaxies lead to rapid classical bulge formation in CDM, the inefficient dynamical friction with MOND in the early-universe galaxies prevents the clumps to coalesce together in the center to form spheroids. This leads to less frequent and less massive classical bulges. Bars in MOND are more frequent and stronger, and have a more constant pattern speed, which modifies significantly the pseudo-bulge morphology. The fraction of pseudo-bulges is expected to be dominant in MOND.

preprint2015arXiv

High Resolution Imaging of PHIBSS z~2 Main Sequence Galaxies in CO J=1-0

We present Karl G. Jansky Very Large Array observations of the CO J=1-0 transition in a sample of four $z\sim2$ main sequence galaxies. These galaxies are in the blue sequence of star-forming galaxies at their redshift, and are part of the IRAM Plateau de Bure HIgh-$z$ Blue Sequence Survey (PHIBSS) which imaged them in CO J=3-2. Two galaxies are imaged here at high signal-to-noise, allowing determinations of their disk sizes, line profiles, molecular surface densities, and excitation. Using these and published measurements, we show that the CO and optical disks have similar sizes in main-sequence galaxies, and in the galaxy where we can compare CO J=1-0 and J=3-2 sizes we find these are also very similar. Assuming a Galactic CO-to-H$_2$ conversion, we measure surface densities of $Σ_{mol}\sim1200$ M$_\odot$pc$^{-2}$ in projection and estimate $Σ_{mol}\sim500-900$ M$_\odot$pc$^{-2}$ deprojected. Finally, our data yields velocity-integrated Rayleigh-Jeans brightness temperature line ratios $r_{31}$ that are approximately unity. In addition to the similar disk sizes, the very similar line profiles in J=1-0 and J=3-2 indicate that both transitions sample the same kinematics, implying that their emission is coextensive. We conclude that in these two main sequence galaxies there is no evidence for significant excitation gradients or a large molecular reservoir that is diffuse or cold and not involved in active star-formation. We suggest that $r_{31}$ in very actively star-forming galaxies is likely an indicator of how well mixed the star formation activity and the molecular reservoir are.

preprint2015arXiv

High-resolution imaging of the molecular outflows in two mergers: IRAS17208-0014 and NGC1614

Galaxy evolution scenarios predict that the feedback of star formation and nuclear activity (AGN) can drive the transformation of gas-rich spiral mergers into ULIRGs, and, eventually, lead to the build-up of QSO/elliptical hosts. We study the role that star formation and AGN feedback have in launching and maintaining the molecular outflows in two starburst-dominated advanced mergers, NGC1614 and IRAS17208-0014, by analyzing the distribution and kinematics of their molecular gas reservoirs. We have used the PdBI array to image with high spatial resolution (0.5"-1.2") the CO(1-0) and CO(2-1) line emissions in NGC1614 and IRAS17208-0014, respectively. The velocity fields of the gas are analyzed and modeled to find the evidence of molecular outflows in these sources and characterize the mass, momentum and energy of these components. While most (>95%) of the CO emission stems from spatially-resolved (~2-3kpc-diameter) rotating disks, we also detect in both mergers the emission from high-velocity line wings that extend up to +-500-700km/s, well beyond the estimated virial range associated with rotation and turbulence. The kinematic major axis of the line wing emission is tilted by ~90deg in NGC1614 and by ~180deg in IRAS17208-0014 relative to their respective rotating disk major axes. These results can be explained by the existence of non-coplanar molecular outflows in both systems. In stark contrast with NGC1614, where star formation alone can drive its molecular outflow, the mass, energy and momentum budget requirements of the molecular outflow in IRAS17208-0014 can be best accounted for by the existence of a so far undetected (hidden) AGN of L_AGN~7x10^11 L_sun. The geometry of the molecular outflow in IRAS17208-0014 suggests that the outflow is launched by a non-coplanar disk that may be associated with a buried AGN in the western nucleus.

preprint2015arXiv

Mapping the inner regions of the polar disk galaxy NGC4650A with MUSE

[abridged] The polar disk galaxy NGC4650A was observed during the commissioning of the MUSE at the ESO VLT to obtain the first 2D map of the velocity and velocity dispersion for both stars and gas. The new MUSE data allow the analysis of the structure and kinematics towards the central regions of NGC4650A, where the two components co-exist. These regions were unexplored by the previous long-slit literature data available for this galaxy. The extended view of NGC~4650A given by the MUSE data is a galaxy made of two perpendicular disks that remain distinct and drive the kinematics right into the very centre of this object. In order to match this observed structure for NGC4650A, we constructed a multicomponent mass model made by the combined projection of two disks. By comparing the observations with the 2D kinematics derived from the model, we found that the modelled mass distribution in these two disks can, on average, account for the complex kinematics revealed by the MUSE data, also in the central regions of the galaxy where the two components coexist. This result is a strong constraint on the dynamics and formation history of this galaxy; it further supports the idea that polar disk galaxies like NGC~4650A were formed through the accretion of material that has different angular momentum.

preprint2015arXiv

Measuring star formation with resolved observations: the test case of M33

Context. Measuring star formation at a local scale is important to constrain star formation laws. Yet, it is not clear whether and how the measure of star formation is affected by the spatial scale at which a galaxy is observed. Aims. We want to understand the impact of the resolution on the determination of the spatially resolved star formation rate (SFR) and other directly associated physical parameters such as the attenuation. Methods. We have carried out a multi-scale, pixel-by-pixel study of the nearby galaxy M33. Assembling FUV, Halpha, 8, 24, 70, and 100 micron maps, we have systematically compared the emission in individual bands with various SFR estimators from a resolution of 33 pc to 2084 pc. Results. We have found that there are strong, scale-dependent, discrepancies up to a factor 3 between monochromatic SFR estimators and Halpha+24 micron. The scaling factors between individual IR bands and the SFR show a strong dependence on the spatial scale and on the intensity of star formation. Finally, strong variations of the differential reddening between the nebular emission and the stellar continuum are seen, depending on the specific SFR (sSFR) and on the resolution. At the finest spatial scales, there is little differential reddening at high sSFR. The differential reddening increases with decreasing sSFR. At the coarsest spatial scales the differential reddening is compatible with the canonical value found for starburst galaxies. Conclusions. Our results confirm that monochromatic estimators of the SFR are unreliable at scales smaller than 1 kpc. Furthermore, the extension of local calibrations to high redshift galaxies presents non-trivial challenges as the properties of these systems may be poorly known.

preprint2015arXiv

Molecular depletion times and the CO-to-H2 conversion factor in metal-poor galaxies

Tracing molecular hydrogen content with carbon monoxide in low-metallicity galaxies has been exceedingly difficult. Here we present a new effort, with IRAM 30-m observations of 12CO(1-0) of a sample of 8 dwarf galaxies having oxygen abundances ranging from 12+logO/H=7.7 to 8.4. CO emission is detected in all galaxies, including the most metal-poor galaxy of our sample (0.1 Zsun); to our knowledge this is the largest number of 12CO(1-0) detections ever reported for galaxies with 12+logO/H<=8 (0.2 Zsun) outside the Local Group. We calculate stellar masses (Mstar) and star-formation rates (SFRs), and analyze our results by combining our observations with galaxy samples from the literature. Extending previous results for a correlation of the molecular gas depletion time, tau(dep), with Mstar and specific SFR (sSFR), we find a variation in tau(dep) of a factor of 200 or more (from <50 Myr to 10 Gyr) over a spread of 1000 in sSFR and Mstar. We exploit the variation of tau(dep) to constrain the CO-to-H2 mass conversion factor alpha(CO) at low metallicity, and assuming a power-law variation find alpha(CO) \propto (Z/Zsun)^1.9, similar to results based on dust continuum measurements compared with gas mass. By including HI measurements, we show that the fraction of total gas mass relative to the baryonic mass is higher in galaxies that are metal poor, of low mass, and of high sSFR. Finally, comparisons of the data with star-formation models of the molecular gas phases suggest that, at metallicities Z/Zsun<=0.2, there are some discrepancies with model predictions.

preprint2015arXiv

Molecular gas content in strongly-lensed z~1.5-3 star-forming galaxies with low IR luminosities

To extend the molecular gas measurements to typical star-forming galaxies (SFGs) with SFR < 40 Msun yr^{-1} and M* < 2.5x10^{10} Msun at z~1.5-3, we have observed CO emission with the IRAM Plateau de Bure Interferometer and 30m telescope for five strongly-lensed galaxies selected from the Herschel Lensing Survey. These observations are combined with a compilation of CO measurements from the literature. We infer the luminosity correction factors r2,1 = 0.81+/-0.20 and r3,1 = 0.57+/-0.15 for the J=2 and J=3 CO transitions, respectively, valid for SFGs at z>1. The combined sample of CO-detected SFGs at z>1 shows a large spread in star formation efficiency (SFE), such that SFE extend beyond the low values of local spirals and overlap the distribution of z>1 sub-mm galaxies. We find that the spread in SFE (or equivalently in molecular gas depletion timescale) is due to primarily the specific star formation rate, but also stellar mass and redshift. Correlations of SFE with the offset from the main-sequence and the compactness of the starburst are less clear. The increase of the molecular gas depletion timescale with M* now revealed by low M* SFGs at z>1 and observed at z=0 is in contrast to the admitted constant molecular gas depletion timescale and the linear Kennicutt-Schmidt relation. We confirm an increase of the molecular gas fraction (fgas) from z~0.2 to z~1.2, followed by a very mild increase toward higher redshifts. At each redshift fgas shows a large dispersion due to the dependence of fgas on M*, producing a gradient of increasing fgas with decreasing M*. We provide the first measure of fgas of z>1 SFGs at the low-M* end (10^{9.4} < M*/Msun < 10^{9.9}), reaching a mean fgas = 0.69+/-0.18, which shows a clear fgas upturn. Finally, we find evidence for a non-universal dust-to-gas ratio among high-redshift SFGs and sub-mm galaxies, local spirals and ULIRGs with near-solar metallicities.

preprint2015arXiv

Probing highly obscured, self-absorbed galaxy nuclei with vibrationally excited HCN

We present high resolution (0."4) IRAM PdBI and ALMA mm and submm observations of the (ultra) luminous infrared galaxies ((U)LIRGs) IRAS17208-0014, Arp220, IC860 and Zw049.057 that reveal intense line emission from vibrationally excited ($ν_2$=1) J=3-2 and 4-3 HCN. The emission is emerging from buried, compact (r<17-70 pc) nuclei that have very high implied mid-infrared surface brightness $>$$5\times 10^{13}$ L$_{\odot}$ kpc$^{-2}$. These nuclei are likely powered by accreting supermassive black holes (SMBHs) and/or hot (>200 K) extreme starbursts. Vibrational, $ν_2$=1, lines of HCN are excited by intense 14 micron mid-infrared emission and are excellent probes of the dynamics, masses, and physical conditions of (U)LIRG nuclei when H$_2$ column densities exceed $10^{24}$ cm$^{-2}$. It is clear that these lines open up a new interesting avenue to gain access to the most obscured AGNs and starbursts. Vibrationally excited HCN acts as a proxy for the absorbed mid-infrared emission from the embedded nuclei, which allows for reconstruction of the intrinsic, hotter dust SED. In contrast, we show strong evidence that the ground vibrational state ($ν$=0), J=3-2 and 4-3 rotational lines of HCN and HCO$^+$ fail to probe the highly enshrouded, compact nuclear regions owing to strong self- and continuum absorption. The HCN and HCO$^+$ line profiles are double-peaked because of the absorption and show evidence of non-circular motions - possibly in the form of in- or outflows. Detections of vibrationally excited HCN in external galaxies are so far limited to ULIRGs and early-type spiral LIRGs, and we discuss possible causes for this. We tentatively suggest that the peak of vibrationally excited HCN emission is connected to a rapid stage of nuclear growth, before the phase of strong feedback.

preprint2015arXiv

Quantifying stellar radial migration in an N-body simulation: blurring, churning, and the outer regions of galaxy discs

Radial stellar migration in galactic discs has received much attention in studies of galactic dynamics and chemical evolution, but remains a dynamical phenomenon that needs to be fully quantified. In this work, using a Tree-SPH simulation of an Sb-type disc galaxy, we quantify the effects of blurring (epicyclic excursions) and churning (change of guiding radius). We quantify migration (either blurring or churning) both in terms of flux (the number of migrators passing at a given radius), and by estimating the population of migrators at a given radius at the end of the simulation compared to non-migrators, but also by giving the distance over which the migration is effective at all radii. We confirm that the corotation of the bar is the main source of migrators by churning in a bar-dominated galaxy, its intensity being directly linked to the episode of a strong bar, in the first 1-3 Gyr of the simulation. We show that within the outer Lindblad resonance (OLR), migration is strongly dominated by churning, while blurring gains progressively more importance towards the outer disc and at later times. Most importantly, we show that the OLR limits the exchange of angular momentum, separating the disc in two distinct parts with minimal or null exchange, except in the transition zone, which is delimited by the position of the OLR at the epoch of the formation of the bar, and at the final epoch. We discuss the consequences of these findings for our understanding of the structure of the Milky Way disc. Because the Sun is situated slightly outside the OLR, we suggest that the solar vicinity may have experienced very limited churning from the inner disc.

preprint2015arXiv

The nuclear gas disk of NGC 1566 dissected by SINFONI and ALMA

We present the results of near-infrared H- and K-band European Southern Observatory SINFONI integral field spectroscopy of the Seyfert galaxy NGC 1566. We investigate the central kpc of this nearby galaxy, concentrating on excitation conditions, morphology, and stellar content. NGC 1566 was selected from our NUGA (-south) sample and is a ringed, spiral galaxy with a stellar bar. We present emission and absorption line measurements in the central kpc of NGC 1566. Broad and narrow Brγ lines were detected. The detection of a broad Brγ component is a clear sign of a super-massive black hole in the center. Blackbody emission temperatures of ~1000 K are indicative of a hot dust component, the torus, in the nuclear region. The molecular hydrogen lines, hydrogen recombination lines, and [FeII] indicate that the excitation at the center is coming from an AGN. The central region is predominantly inhabited by molecular gas, dust, and an old K-M type giant stellar population. The molecular gas and stellar velocity maps both show a rotation pattern. The molecular gas velocity field shows a perturbation toward the center that is typical for bars or spiral density waves. The molecular gas species of warm H_2(1-0)S(1) and cold ^{12}CO(3-2) gas trace a nuclear gas disk of about 3" in radius with a nuclear spiral reaching toward the nucleus. From the equivalent width of H_2(1-0)S(1) a molecular ring with r<~3" can be inferred. This spiral seems to be an instrument that allows gas to fall toward the nucleus down to <50 pc scales. The excitation of molecular hydrogen in the nuclear gas disk is not clear but diagnostic diagrams show a distinction between the nuclear region and a <9 Myr old star forming region at the southwestern spiral arm. Possibly shocked gas is detected ~2" from the center, which is visible in dispersion maps of H$_2$(1-0)S(1) and ^{12}CO(3-2) and in the 0.87 mm continuum.

preprint2014arXiv

[CII] and 12CO(1-0) Emission Maps in HLSJ091828.6+514223: A Strongly Lensed Interacting System at z=5.24

We present Submillimeter Array (SMA) [CII] 158um and Jansky Very Large Array (JVLA) $^{12}$CO(1-0) line emission maps for the bright, lensed, submillimeter source at $z=5.2430$ behind Abell 773: HLSJ091828.6+514223 (HLS0918). We combine these measurements with previously reported line profiles, including multiple $^{12}$CO rotational transitions, [CI], water and [NII], providing some of the best constraints on the properties of the interstellar medium (ISM) in a galaxy at $z>5$. HLS0918 has a total far-infrared (FIR) luminosity L_FIR(8-1000um) = (1.6$\pm$0.1)x10^14 L_sun/mu, where the total magnification mu_total = 8.9$\pm$1.9, via a new lens model from the [CII] and continuum maps. Despite a HyLIRG luminosity, the FIR continuum shape resembles that of a local LIRG. We simultaneously fit all of the observed spectral line profiles, finding four components which correspond cleanly to discrete spatial structures identified in the maps. The two most redshifted spectral components occupy the nucleus of a massive galaxy, with a source plane separation <1 kpc. The reddest dominates the continuum map (de-magnified L_FIR = (1.1$\pm$0.2)x10^13 L_sun), and excites strong water emission in both nuclear components via a powerful FIR radiation field from the intense star formation. A third star-forming component is most likely a region of a merging companion (dV ~ 500 km/s) exhibiting generally similar gas properties. The bluest component originates from a spatially distinct region, and photo-dissociation region (PDR) analysis suggests that it is lower density, cooler and forming stars less vigorously than the other components. Strikingly, it has very strong [NII] emission which may suggest an ionized, molecular outflow. This comprehensive view of gas properties and morphology in HLS0918 previews the science possible for a large sample of high-redshift galaxies once ALMA attains full sensitivity.

preprint2014arXiv

A $10^{10}$ Solar Mass Flow of Molecular Gas in the Abell 1835 Brightest Cluster Galaxy

We report ALMA Early Science observations of the Abell 1835 brightest cluster galaxy (BCG) in the CO (3-2) and CO (1-0) emission lines. We detect $5\times 10^{10}~\rm M_\odot$ of molecular gas within 10 kpc of the BCG. Its ensemble velocity profile width of $\sim 130 ~\rm km~s^{-1}$ FWHM is too narrow for the molecular cloud sto be supported in the galaxy by dynamic pressure. The gas may instead be supported in a rotating, turbulent disk oriented nearly face-on. Roughly $10^{10}~\rm M_\odot$ of molecular gas is projected $3-10 ~\rm kpc$ to the north-west and to the east of the nucleus with line of sight velocities lying between $-250 ~\rm km~s^{-1}$ to $+480 ~\rm km~s^{-1}$ with respect to the systemic velocity. The high velocity gas may be either inflowing or outflowing. However, the absence of high velocity gas toward the nucleus that would be expected in a steady inflow, and its bipolar distribution on either side of the nucleus, are more naturally explained as outflow. Star formation and radiation from the AGN are both incapable of driving an outflow of this magnitude. If so, the molecular outflow may be associated a hot outflow on larger scales reported by Kirkpatrick and colleagues. The molecular gas flow rate of approximately $200~\rm M_\odot ~yr^{-1}$ is comparable to the star formation rate of $100-180~\rm M_\odot ~yr^{-1}$ in the central disk. How radio bubbles would lift dense molecular gas in their updrafts, how much gas will be lost to the BCG, and how much will return to fuel future star formation and AGN activity are poorly understood. Our results imply that radio-mechanical (radio mode) feedback not only heats hot atmospheres surrounding elliptical galaxies and BCGs, it is able to sweep higher density molecular gas away from their centers.

preprint2014arXiv

ALMA observations of feeding and feedback in nearby Seyfert galaxies: an AGN-driven outflow in NGC 1433

We report ALMA observations of CO(3-2) emission in the Seyfert 2 double-barred galaxy NGC1433, at the unprecedented spatial resolution of 0.5"=24 pc. Our aim is to probe AGN feeding and feedback phenomena through the morphology and dynamics of the gas inside the central kpc. The CO map, which covers the whole nuclear region (nuclear bar and ring), reveals a nuclear gaseous spiral structure, inside the nuclear ring encircling the nuclear stellar bar. This gaseous spiral is well correlated with the dusty spiral seen in Hubble Space Telescope images. The nuclear spiral winds up in a pseudo-ring at 200 pc radius, which might correspond to the inner ILR. Continuum emission is detected at 0.87 mm only at the very centre, and its origin is more likely thermal dust emission than non-thermal emission from the AGN. It might correspond to the molecular torus expected to exist in this Seyfert 2 galaxy. The HCN(4-3) and HCO+(4-3) lines were observed simultaneously, but only upper limits are derived, with a ratio to the CO(3-2) line lower than 1/60 at 3sigma, indicating a relatively low abundance of very dense gas. The kinematics of the gas over the nuclear disk reveal rather regular rotation only slightly perturbed by streaming motions due to the spiral; the primary and secondary bars are too closely aligned with the galaxy major or minor axis to leave a signature in the projected velocities. Near the nucleus, there is an intense high-velocity CO emission feature redshifted to 200 km/s (if located in the plane), with a blue-shifted counterpart, at 2" (100 pc) from the centre. While the CO spectra are quite narrow in the centre, this wide component is interpreted as an outflow involving a molecular mass of 3.6 10^6 Mo and a flow rate 7 Mo/yr. The flow could be in part driven by the central star formation, but is mainly boosted by the AGN through its wind or radio jets.

preprint2014arXiv

ALMA reveals the feeding of the Seyfert 1 nucleus in NGC 1566

We report ALMA observations of CO(3-2) emission in the Seyfert 1 galaxy NGC 1566, at a spatial resolution of 25 pc. Our aim is to investigate the morphology and dynamics of the gas inside the central kpc, and to probe nuclear fueling and feedback phenomena. NGC 1566 has a nuclear bar of 1.7 kpc radius and a conspicuous grand design spiral starting from this radius. The ALMA field of view, of diameter 0.9 kpc, lies well inside the nuclear bar and reveals a molecular trailing spiral structure from 50 to 300~pc in size, which is contributing to fuel the nucleus, according to its negative gravity torques. The spiral starts with a large pitch angle from the center and then winds up in a pseudo-ring at the inner Lindblad resonance (ILR) of the nuclear bar. This is the first time that a trailing spiral structure is clearly seen driving the gas inwards inside the ILR ring of the nuclear bar. This phenomenon shows that the massive central black hole has a significant dynamical influence on the gas, triggering its fueling. The gaseous spiral is well correlated with the dusty spiral seen through extinction in HST images, and also with a spiral feature emitting 0.87mm continuum. This continuum emission must come essentially from cold dust heated by the interstellar radiation field. The HCN(4-3) and HCO+(4-3) lines were simultaneously mapped and detected in the nuclear spiral. The HCO+(4-3) line is 3 times stronger than the HCN(4-3), as expected when star formation excitation dominates over active galactic nucleus (AGN) heating. The CO(3-2)/HCO+(4-3) integrated intensity ratio is \sim 100. The molecular gas is in remarkably regular rotation, with only slight non-circular motions at the periphery of the nuclear spiral arms. These perturbations are quite small, and no outflow nor AGN feedback is detected.

preprint2014arXiv

An ALMA Early Science survey of molecular absorption lines toward PKS1830-211 -- Analysis of the absorption profiles

We present the first results of an ALMA spectral survey of strong absorption lines for common interstellar species in the z=0.89 molecular absorber toward the lensed blazar PKS1830-211. The dataset brings essential information on the structure and composition of the absorbing gas in the foreground galaxy. In particular, we find absorption over large velocity intervals (gtrsim 100 km/s) toward both lensed images of the blazar. This suggests either that the galaxy inclination is intermediate and that we sample velocity gradients or streaming motions in the disk plane, that the molecular gas has a large vertical distribution or extraplanar components, or that the absorber is not a simple spiral galaxy but might be a merger system. The number of detected species is now reaching a total of 42 different species plus 14 different rare isotopologues toward the SW image, and 14 species toward the NE line-of-sight. The abundances of CH, H2O, HCO+, HCN, and NH3 relative to H2 are found to be comparable to those in the Galactic diffuse medium. Of all the lines detected so far toward PKS1830-211, the ground-state line of ortho-water has the deepest absorption. We argue that ground-state lines of water have the best potential for detecting diffuse molecular gas in absorption at high redshift.

preprint2014arXiv

Bulge formation in disk galaxies with MOND

The formation of galaxies and their various components can be stringent tests of dark matter models and of gravity theories. In the standard cold dark matter (CDM) model, spheroids are formed through mergers in a strongly hierarchical scenario, and also in the early universe through dynamical friction in clumpy galaxies. More secularly, pseudo-bulges are formed by the inner vertical resonance with bars. The high efficiency of bulge formation is in tension with observations in the local universe of a large amount of bulge-less spiral galaxies. In the present work, the formation of bulges in very gas-rich galaxies, as those in the early universe, is studied in the Milgrom's MOdified Newtonian Dynamics (MOND), through multi-grid simulations of the non-linear gravity, including the gas dissipation, star formation and feedback. Clumpy disks are rapidly formed, as in their Newtonian equivalent systems. However, the dynamical friction is not as efficient, in the absence of dark matter halos, and the clumps have no time to coalesce into the center to form bulges, before they are eroded by stellar feedback and shear forces. Previous work has established that mergers are less frequent in MOND, and classical bulges are expected less massive. It is now shown that gas-rich clumpy galaxies in the early universe do not form bulges. Since pseudo-bulges are formed with a similar rate as in the Newtonian equivalent systems, it can be expected that the contribution of pseudo-bulges is significantly higher in MOND.

preprint2014arXiv

CO map and steep Kennicutt-Schmidt relation in the extended UV disk of M63

Results from the UV satellite GALEX revealed large extensions of disks in some nearby spiral galaxies, extending out to 3 to 4 times the isophotal radius, r25. M63 is a remarkable example of a spiral galaxy with one of the most extended UV disks, so it offers the opportunity to search for the molecular gas and characterize the star formation in outer disk regions as revealed by the UV emission. We obtained deep CO(1-0) and CO(2-1) observations on the IRAM 30 m telescope along the major axis of the M63 disk from the center out to the galactocentric radius rgal = 1.6 r25 and over a bright UV region at rgal = 1.36 r25. CO(1-0) is detected all along the M63 major axis out to r25, and CO(2-1) is confined to rgal = 0.68 r25, which may betray lower excitation temperatures in the outer disk. CO(1-0) is also detected in the external bright UV region of M63. The radial profiles of the CO emission and of the Halpha, 24 micron, NUV and FUV star formation tracers and HI taken from the literature show a severe drop with the galactocentric radius, such that beyond r25 they are all absent with the exception of a faint UV emission and HI. The CO emission detection in the external UV region, where the UV flux is higher than the UV flux observed beyond r25, highlights a tight correlation between the CO and UV fluxes, namely the amount of molecular gas and the intensity of star formation. This external UV region is dominated by the atomic gas, suggesting that HI is more likely the precursor of H2 rather than the product of UV photodissociation. A broken power law needs to be invoked to describe the Kennicutt-Schmidt (K-S) relation of M63 from the center of the galaxy out to rgal = 1.36 r25. While all along the major axis out to r25 the K-S relation is almost linear, in the external UV region the SFR regime is highly nonlinear and characterized by a steep K-S relation and very low star formation efficiency.

preprint2014arXiv

Cold gas in group-dominant elliptical galaxies

We present IRAM 30m telescope observations of the CO(1-0) and (2-1) lines in a sample of 11 group-dominant elliptical galaxies selected from the CLoGS nearby groups sample. Our observations confirm the presence of molecular gas in 4 of the 11 galaxies at >4 sigma significance, and combining these with data from the literature we find a detection rate of 43+-14%, comparable to the detection rate for nearby radio galaxies, suggesting that group-dominant ellipticals may be more likely to contain molecular gas than their non-central counterparts. Those group-dominant galaxies which are detected typically contain ~2x10^8 Msol of molecular gas, and although most have low star formation rates (<1 Msol/yr) they have short depletion times, indicating that the gas must be replenished on timescales ~100 Myr. Almost all of the galaxies contain active nuclei, and we note while the data suggest that CO may be more common in the most radio-loud galaxies, the mass of molecular gas required to power the active nuclei through accretion is small compared to the masses observed. We consider possible origin mechanisms for the gas, through cooling of stellar ejecta within the galaxies, group-scale cooling flows, and gas-rich mergers, and find probable examples of each type within our sample, confirming that a variety of processes act to drive the build up of molecular gas in group-dominant ellipticals.

preprint2014arXiv

Combined CO & Dust Scaling Relations of Depletion Time and Molecular Gas Fractions with Cosmic Time, Specific Star Formation Rate and Stellar Mass

We combine molecular gas masses inferred from CO emission in 500 star forming galaxies (SFGs) between z=0 and 3, from the IRAM-COLDGASS, PHIBSS1/2 and other surveys, with gas masses derived from Herschel far-IR dust measurements in 512 galaxy stacks over the same stellar mass/redshift range. We constrain the scaling relations of molecular gas depletion time scale (tdepl) and gas to stellar mass ratio (Mmolgas/M*) of SFGs near the star formation main-sequence with redshift, specific star formation rate (sSFR) and stellar mass (M*). The CO- and dust-based scaling relations agree remarkably well. This suggests that the CO-H2 mass conversion factor varies little within 0.6dex of the main sequence (sSFR(ms,z,M*)), and less than 0.3dex throughout this redshift range. This study builds on and strengthens the results of earlier work. We find that tdepl scales as (1+z)^-0.3 *(sSFR/sSFR(ms,z,M*))^-0.5, with little dependence on M*. The resulting steep redshift dependence of Mmolgas/M* ~(1+z)^3 mirrors that of the sSFR and probably reflects the gas supply rate. The decreasing gas fractions at high M* are driven by the flattening of the SFR-M* relation. Throughout the redshift range probed a larger sSFR at constant M* is due to a combination of an increasing gas fraction and a decreasing depletion time scale. As a result galaxy integrated samples of the Mmolgas-SFR rate relation exhibit a super-linear slope, which increases with the range of sSFR. With these new relations it is now possible to determine Mmolgas with an accuracy of 0.1dex in relative terms, and 0.2dex including systematic uncertainties.

preprint2014arXiv

Detection of chloronium and measurement of the 35Cl/37Cl isotopic ratio at z=0.89 toward PKS1830-211

We report the first extragalactic detection of chloronium (H2Cl+), in the z=0.89 absorber in front of the lensed blazar PKS1830-211. The ion is detected through its 1_11-0_00 line along two independent lines of sight toward the North-East and South-West images of the blazar. The relative abundance of H2Cl+ is significantly higher (by a factor ~7) in the NE line of sight, which has a lower H2/H fraction, indicating that H2Cl+ preferably traces the diffuse gas component. From the ratio of the H2^35Cl+ and H2^37Cl+ absorptions toward the SW image, we measure a 35Cl/37Cl isotopic ratio of 3.1 (-0.2; +0.3) at z=0.89, similar to that observed in the Galaxy and the solar system.

preprint2014arXiv

Exceptional AGN-driven turbulence inhibits star formation in the 3C 326N radio-galaxy

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).

preprint2014arXiv

Fueling the central engine of radio galaxies. III. Molecular gas and star formation efficiency of 3C 293

Aims. We investigate the different manifestations of AGN feedback in the evolved, powerful radio source 3C293 and their impact on the molecular gas of its host galaxy, which harbors young star-forming regions and fast outflows of HI and ionized gas. Methods. We study the distribution and kinematics of the molecular gas of 3C293 using high spatial resolution observations of the CO(1-0) and CO(2-1) lines, and the 3 and 1mm continuum taken with the IRAM PdBI. We mapped the molecular gas of 3C293 and compared it with the dust and star-formation images of the host. We searched for signatures of outflow motions in the CO kinematics, and reexamined the evidence of outflowing gas in the HI spectra. We also derived the star formation rate (SFR) and efficiency (SFE) of the host with all available SFR tracers from the literature, and compared them with the SFE of young and evolved radio galaxies and normal star-forming galaxies. Results. The CO(1-0) emission line shows that the molecular gas in 3C293 is distributed along a massive (2.2E10 Msun) warped disk with diameter of 21 kpc that rotates around the AGN. Our data show that the dust and the star formation are clearly associated with the CO disk. The CO(2-1) emission is located in the inner 7 kpc (diameter) region around the AGN, coincident with the inner part of the CO(1-0) disk. Both the CO(1-0) and CO(2-1) spectra reveal the presence of an absorber against the central regions of 3C293 that is associated with the disk. We do not detect any fast (>500 km/s) outflow motions in the cold molecular gas. The host of 3C293 shows an SFE consistent with the Kennicutt-Schmidt law. The apparently low SFE of evolved radio galaxies may be caused by an underestimation of the SFR and/or an overestimation of the molecular gas densities in these sources. We find no signatures of AGN feedback in the molecular gas of 3C293.

preprint2014arXiv

Future mmVLBI Research with ALMA: A European vision

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

preprint2014arXiv

Galaxy size trends as a consequence of cosmology

We show that recently documented trends in galaxy sizes with mass and redshift can be understood in terms of the influence of underlying cosmic evolution; a holistic view which is complimentary to interpretations involving the accumulation of discreet evolutionary processes acting on individual objects. Using standard cosmology theory, supported with results from the Millennium simulations, we derive expected size trends for collapsed cosmic structures, emphasising the important distinction between these trends and the assembly paths of individual regions. We then argue that the observed variation in the stellar mass content of these structures can be understood to first order in terms of natural limitations of cooling and feedback. But whilst these relative masses vary by orders of magnitude, galaxy and host radii have been found to correlate linearly. We explain how these two aspects will lead to galaxy sizes that closely follow observed trends and their evolution, comparing directly with the COSMOS and SDSS surveys. Thus we conclude that the observed minimum radius for galaxies, the evolving trend in size as a function of mass for intermediate systems, and the observed increase in the sizes of massive galaxies, may all be considered an emergent consequence of the cosmic expansion.

preprint2014arXiv

Heating of the molecular gas in the massive outflow of the local ultraluminous-infrared and radio-loud galaxy 4C12.50

We present a comparison of the molecular gas properties in the outflow vs. in the ambient medium of the local prototype radio-loud and ultraluminous-infrared galaxy 4C12.50 (IRAS13451+1232), using new data from the IRAM Plateau de Bure interferometer and 30m telescope, and the Herschel space telescope. Previous H_2 (0-0) S(1) and S(2) observations with the Spitzer space telescope had indicated that the warm (~400K) molecular gas in 4C12.50 is made up of a 1.4(+-0.2)x10^8 M_sun ambient reservoir and a 5.2(+-1.7)x10^7 M_sun outflow. The new CO(1-0) data cube indicates that the corresponding cold (25K) H_2 gas mass is 1.0(+-0.1)x10^10 M_sun for the ambient medium and <1.3x10^8 M_sun for the outflow, when using a CO-intensity-to-H_2-mass conversion factor alpha of 0.8 M_sun /(K km/s pc^2). The combined mass outflow rate is high, 230-800 M_sun/yr, but the amount of gas that could escape the galaxy is low. A potential inflow of gas from a 3.3(+-0.3)x10^8 M_sun tidal tail could moderate any mass loss. The mass ratio of warm-to-cold molecular gas is >= 30 times higher in the outflow than in the ambient medium, indicating that a non-negligible fraction of the accelerated gas is heated to temperatures at which star formation is inefficient. This conclusion is robust against the use of different alpha factor values, and/or different warm gas tracers (H_2 vs. H_2 plus CO): with the CO-probed gas mass being at least 40 times lower at 400K than at 25K, the total warm-to-cold mass ratio is always lower in the ambient gas than in the entrained gas. Heating of the molecular gas could facilitate the detection of new outflows in distant galaxies by enhancing their emission in intermediate rotational number CO lines.

preprint2014arXiv

HI study of the environment around ESO 243-49, the host galaxy of an intermediate mass black hole

The lenticular galaxy ESO 243-49 hosts the ultraluminous X-ray source HLX-1, the best candidate intermediate mass black hole (IMBH) currently known. The environments of IMBHs remain unknown, however the proposed candidates include the nuclei of dwarf galaxies or globular clusters. Evidence at optical wavelengths points at HLX-1 being the remnant of an accreted dwarf galaxy. Here we report the Australia Telescope Compact Array radio observations of HI emission in and around ESO 243-49 searching for signatures of a recent merger event. No HI line emission is detected in ESO 243-49 with a 5$σ$ upper limit on the HI gas mass of a few $10^8 M_{\odot}$. A likely reason for this non-detection is the cluster environment depleting ESO 243-49's HI gas reservoir. The upper limit is consistent with an interpretation of HLX-1 as a dwarf satellite of ESO 243-49, however more sensitive observations are required for a detection. We detect ~$5 \times 10^8 M_{\odot}$ of HI gas in the peculiar spiral galaxy AM 0108-462, located at a projected distance of ~170 kpc from ESO 243-49. This amount of HI gas is ~10 times less than in spiral galaxies with similar optical and near-infrared properties in the field, strengthening the conclusion that the cluster environment indeed depletes the HI gas reservoir of these two galaxies. Here we also report observations of AM 0108-462 in several optical and near-infrared bands using the Magellan 6.5 m telescopes, and archival X-ray and ultraviolet observations with XMM-Newton and Swift. These data combined with the HI line data suggest it is likely that AM 0108-462 is experiencing a merger event.

preprint2014arXiv

Jet-induced star formation in 3C 285 and Minkowski Object

How efficiently star formation proceeds in galaxies is still an open question. Recent studies suggest that AGN can regulate the gas accretion and thus slow down star formation (negative feedback). However, evidence of AGN positive feedback has also been observed in a few radio galaxies (eg. Centaurus A). Here we present CO observations of 3C 285 and Minkowski Object (MO), which are examples of jet-induced star formation. A spot (named 09.6) aligned with the 3C 285 radio jet, at a projected distance of ~70 kpc from the galaxy centre, shows star formation, detected in optical emission. MO is located along the jet of NGC 541 and also shows star formation. To know the distribution of molecular gas along the jets is a way to study the physical processes at play in the AGN interaction with the intergalactic medium. We observed CO lines in 3C 285, NGC 541, 09.6 and MO with the IRAM-30m telescope. In the central galaxies, the spectra present a double-horn profile, typical of a rotation pattern, from which we are able to estimate the molecular gas density profile of the galaxy. The molecular gas appears to be in a compact reservoir. In addition, no kinematic signature of a molecular outflow is detected by the 30m-telescope. Interestingly, 09.6 and MO are not detected in CO. The cold gas mass upper limits are consistent with a star formation induced by the compression of dense ambient material by the jet. The depletion time scales are of the order of and even smaller than what is found in 3C 285, NGC 541 and local spiral galaxies (10^9 yr). The molecular gas surface density in 09.6 follows a Schmidt-Kennicutt law if the emitting region is very compact, while MO is found to have a much higher SFE (very short depletion time). Higher sensitivity and spatial resolution are necessary to detect CO in the spots of star formation, and map the emission in these jet-induced star forming regions.

preprint2014arXiv

Mapping a stellar disk into a boxy bulge: The outside-in part of the Milky Way bulge formation

By means of idealized, dissipationless N-body simulations which follow the formation and subsequent buckling of a stellar bar, we study the characteristics of boxy/peanut-shaped bulges and compare them with the properties of the stellar populations in the Milky Way bulge. The main results of our modeling, valid for the general family of boxy/peanut shaped bulges, are the following: (i) because of the spatial redistribution in the disk initiated at the epoch of bar formation, stars from the innermost regions to the outer Lindblad resonance of the stellar bar are mapped into a boxy bulge; (ii) the contribution of stars to the local bulge density depends on their birth radius: stars born in the innermost disk tend to dominate the innermost regions of the boxy bulge, while stars originating closer to the OLR are preferentially found in the outer regions of the boxy/peanut structure; (iii) stellar birth radii are imprinted in the bulge kinematics, the larger the birth radii of stars ending up in the bulge, the greater their rotational support and the higher their line-of- sight velocity dispersions (but note that this last trend depends on the bar viewing angle); (iv) the higher the classical bulge-over-disk ratio, the larger its fractional contribution of stars at large vertical distance from the galaxy mid-plane. (ABRIDGED) On the basis of their chemical and kinematic characteristics, the results of our modeling suggests that the populations A, B and C, as defined by the ARGOS survey, can be associated, respectively, with the inner thin disk, to the young thick and to the old thick disk, following the nomenclature recently suggested for stars in the solar neighborhood by Haywood et al. (2013).

preprint2014arXiv

Metallicity and kinematics of the bar in-situ

Constraints on the Galactic bulge/bar structure and formation history from stellar kinematics and metallicities mainly come from relatively high-latitude fields (|b|>4) where a complex mix of stellar population is seen. We aim here to constrain the formation history of the Galactic bar by studying the radial velocity and metallicity distributions of stars in-situ (|b|<1). We observed red clump stars in four fields along the bar's major axis (l=10,6,-6 and b=0 plus a field at l=0,b=1) with low-resolution spectroscopy from VLT/FLAMES, observing around the CaII triplet. We developed robust methods for extracting radial velocity and metallicity estimates from these low signal-to-noise spectra. We derived distance probability distributions using Bayesian methods rigorously handling the extinction law. We present radial velocities and metallicity distributions, as well as radial velocity trends with distance. We observe an increase in the radial velocity dispersion near the Galactic plane. We detect the streaming motion of the stars induced by the bar in fields at l=+/-6, the highest velocity components of this bar stream being metal-rich ([Fe/H]~0.2 dex). Our data is consistent with a bar inclined at 26+/-3 from the Sun-Galactic centre line. We observe a significant fraction of metal-poor stars, in particular in the field at l=0,b=1. We confirm the flattening of the metallicity gradient along the minor axis when getting closer to the plane, with a hint that it could actually be inverted. Our stellar kinematics corresponds to the expected behaviour of a bar issued from the secular evolution of the Galactic disc. The mix of several populations, seen further away from the plane, is also seen in the bar in-situ since our metallicity distributions highlight a different spatial distribution between metal-poor and metal-rich stars, the more metal-poor stars being more centrally concentrated.

preprint2014arXiv

Models of AGN feedback

The physical processes responsible of sweeping up the surrounding gas in the host galaxy of an AGN, and able in some circumstances to expel it from the galaxy, are not yet well known. The various mechanisms are briefly reviewed: quasar or radio modes, either momentum-conserving outflows, energy-conserving outflows, or intermediate. They are confronted to observations, to know whether they can explain the M-sigma relation, quench the star formation or whether they can also provide some positive feedback and how the black hole accretion history is related to that of star formation.

preprint2014arXiv

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

We present a detailed analysis of 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 disk (CND) and the ~ 1.3 kpc starburst ring (SB ring) of NGC~1068, a nearby (D = 14 Mpc) Seyfert 2 barred galaxy. We aim at determining the physical characteristics of the dense gas present in the CND and 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. 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 first only the CO transitions, and then all the available molecules in order 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. We find that in general the gas in the CND is very dense (> 10^5 cm^-3) and hot (T> 150K), 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. 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. 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.

preprint2014arXiv

Molecular line emission in NGC1068 imaged with ALMA. I An AGN-driven outflow in the dense molecular gas

We investigate the fueling and the feedback of star formation and nuclear activity in NGC1068, a nearby (D=14Mpc) Seyfert 2 barred galaxy, by analyzing the distribution and kinematics of the molecular gas in the disk. We have used ALMA to map the emission of a set of dense molecular gas tracers (CO(3-2), CO(6-5), HCN(4-3), HCO+(4-3) and CS(7-6)) and their underlying continuum emission in the central r ~ 2kpc of NGC1068 with spatial resolutions ~ 0.3"-0.5" (~ 20-35pc). Molecular line and dust continuum emissions are detected from a r ~ 200pc off-centered circumnuclear disk (CND), from the 2.6kpc-diameter bar region, and from the r ~ 1.3kpc starburst (SB) ring. Most of the emission in HCO+, HCN and CS stems from the CND. Molecular line ratios show dramatic order-of-magnitude changes inside the CND that are correlated with the UV/X-ray illumination by the AGN, betraying ongoing feedback. The gas kinematics from r ~ 50pc out to r ~ 400pc reveal a massive (M_mol ~ 2.7 (+0.9, -1.2) x 10^7 Msun) outflow in all molecular tracers. The tight correlation between the ionized gas outflow, the radio jet and the occurrence of outward motions in the disk suggests that the outflow is AGN-driven. The outflow rate estimated in the CND, dM/dt ~ 63 (+21, -37) Msun yr^-1, is an order of magnitude higher than the star formation rate at these radii, confirming that the outflow is AGN-driven. The power of the AGN is able to account for the estimated momentum and kinetic luminosity of the outflow. The CND mass load rate of the CND outflow implies a very short gas depletion time scale of <=1 Myr.

preprint2014arXiv

New constraints on dust emission and UV attenuation of z=6.5-7.5 galaxies from millimeter observations

We have targeted two recently discovered Lyman break galaxies (LBGs) to search for dust continuum and [CII] 158 micron line emission. The strongly lensed z~6.8 LBG A1703-zD1 behind the galaxy cluster Abell 1703, and the spectroscopically confirmed z=7.508 LBG z8-GND-5296 in the GOODS-N field have been observed with the Plateau de Bure interferometer (PdBI) at 1.2mm. These observations have been combined with those of three z>6.5 Lya emitters (named HCM6A, Himiko, and IOK-1), for which deep measurements were recently obtained with the PdBI and ALMA. [CII] is undetected in both galaxies, providing a deep upper limit for Abell1703-zD1, comparable to recent ALMA non-detections. Dust continuum emission from Abell1703-zD1 and z8-GND-5296 is not detected with an rms of 0.12 and 0.16 mJy/beam. From these non-detections we derive upper limits on their IR luminosity and star formation rate, dust mass, and UV attenuation. Thanks to strong gravitational lensing the limit for Abell1703-zD1 is probing the sub-LIRG regime ($L_{IR} <8.1 \times 10^{10}$ Lsun) and very low dust masses ($M_d<1.6 \times 10^7$ Msun). We find that all five galaxies are compatible with the Calzetti IRX-$β$ relation, their UV attenuation is compatible with several indirect estimates from other methods (the UV slope, extrapolation of the attenuation measured from the IR/UV ratio at lower redshift, and SED fits), and the dust-to-stellar mass ratio is not incompatible with that of galaxies from z=0 to 3. For their stellar mass the high-z galaxies studied here have an attenuation below the one expected from the mean relation of low redshift (z<1.5) galaxies. More and deeper (sub)-mm data are clearly needed to directly determine the UV attenuation and dust content of the dominant population of high-z star-forming galaxies and to establish more firmly their dependence on stellar mass, redshift, and other properties.

preprint2014arXiv

Polar-ring galaxies: the SDSS view on the symbiotic galaxies

Polar-ring galaxies are multi-spin systems, showing star formation in a blue late-type component, perpendicular to a red early-type one, revealing how galaxy formation can sometimes occur in successive steps. We perform two-dimensional decomposition in the $g$, $r$, $i$ bandpasses of 50 polar-ring galaxies (PRGs) from the Sloan Digital Sky Survey. Each object was fit with a Sérsic host galaxy and a Sérsic ring. Our general results are: (i) The central (host) galaxies of the PRGs are non-dwarf sub-$L^{\ast}$ galaxies with colors typical for early-type galaxies. (ii) Polar structures in our sample are, on average, fainter and bluer than their host galaxies. (iii) In most galaxies, the stellar mass M$_*$ of the polar component is not negligible in comparison with that of the host. (iv) The distributions of the host galaxies on the size -- luminosity and Kormendy diagrams are shifted by $\sim 1^m$ to fainter magnitudes in comparison with E/S0 galaxies. It means that the PRGs hosts are more similar to quenched disks than to ordinary early-type galaxies. (v) All the PRGs in our sample are detected in mid-infrared by WISE, and we derive from the 22$μ$m luminosity their star formation rate (SFR). Their SFR/M$_*$ ratio is larger than for the early-type galaxy sample of Atlas$^{\rm 3D}$, showing that the star forming disk brings a significant contribution to the new stars. Globally, PRGs appear frequently on the green valley in the mass-color diagram, revealing the symbiotic character between a red-sequence host and a blue cloud ring.

preprint2014arXiv

Reconstructing the star formation history of the Milky Way disc(s) from chemical abundances

We develop a chemical evolution model in order to study the star formation history of the Milky Way. Our model assumes that the Milky Way is formed from a closed box-like system in the inner regions, while the outer parts of the disc experience some accretion. Unlike the usual procedure, we do not fix the star formation prescription (e.g. Kennicutt law) in order to reproduce the chemical abundance trends. Instead, we fit the abundance trends with age in order to recover the star formation history of the Galaxy. Our method enables one to recover with unprecedented accuracy the star formation history of the Milky Way in the first Gyrs, in both the inner (R<7-8kpc) and outer (R>9-10kpc) discs as sampled in the solar vicinity. We show that, in the inner disc, half of the stellar mass formed during the thick disc phase, in the first 4-5 Gyr. This phase was followed by a significant dip in the star formation activity (at 8-9 Gyr) and a period of roughly constant lower level star formation for the remaining 8 Gyr. The thick disc phase has produced as many metals in 4 Gyr as the thin disc in the remaining 8 Gyr. Our results suggest that a closed box model is able to fit all the available constraints in the inner disc. A closed box system is qualitatively equivalent to a regime where the accretion rate, at high redshift, maintains a high gas fraction in the inner disc. In such conditions, the SFR is mainly governed by the high turbulence of the ISM. By z~1 it is possible that most of the accretion takes place in the outer disc, while the star formation activity in the inner disc is mostly sustained by the gas not consumed during the thick disc phase, and the continuous ejecta from earlier generations of stars. The outer disc follows a star formation history very similar to that of the inner disc, although initiated at z~2, about 2 Gyr before the onset of the thin disc formation in the inner disc.

preprint2014arXiv

Science Case and Requirements for the MOSAIC Concept for a Multi-Object Spectrograph for the European Extremely Large Telescope

Over the past 18 months we have revisited the science requirements for a multi-object spectrograph (MOS) for the European Extremely Large Telescope (E-ELT). These efforts span the full range of E-ELT science and include input from a broad cross-section of astronomers across the ESO partner countries. In this contribution we summarise the key cases relating to studies of high-redshift galaxies, galaxy evolution, and stellar populations, with a more expansive presentation of a new case relating to detection of exoplanets in stellar clusters. A general requirement is the need for two observational modes to best exploit the large (>40 sq. arcmin) patrol field of the E-ELT. The first mode ('high multiplex') requires integrated-light (or coarsely resolved) optical/near-IR spectroscopy of >100 objects simultaneously. The second ('high definition'), enabled by wide-field adaptive optics, requires spatially-resolved, near-IR of >10 objects/sub-fields. Within the context of the conceptual study for an ELT-MOS called MOSAIC, we summarise the top-level requirements from each case and introduce the next steps in the design process.

preprint2014arXiv

The IRAM M33 CO(2-1) Survey - A complete census of the molecular gas out to 7 kpc

In order to study the ISM and the interplay between the atomic and molecular components in a low-metallicity environment, we present a complete high angular and spectral resolution map and data cube of the 12CO(2-1) emission from the Local Group galaxy M33. Its metallicity is roughly half-solar, such that we can compare its ISM with that of the Milky Way with the main changes being the metallicity and the gas mass fraction. The data have a 12" angular resolution (50pc) with a spectral resolution of 2.6 km/s and a mean noise level of 20 mK per channel in antenna temperature. A radial cut along the major axis was also observed in the 12CO(1-0) line. The CO data cube and integrated intensity map are optimal when using HI data to define the baseline window and the velocities over which the CO emission is integrated. Great care was taken when building these maps, testing different windowing and baseline options and investigating the effect of error beam pickup. The total CO(2-1) luminosity is 2.8e7 K km/s pc2, following the spiral arms in the inner disk. There is no clear variation in the CO(2-1/1-0) intensity ratio with radius and the average value is roughly 0.8. The total molecular gas mass is estimated, using a N(H2)/Ico(1-0)=4e20 cm-2/(K km/s) conversion factor, to be 3.1e8 Msol. The CO spectra in the cube were shifted to zero velocity by subtracting the velocity of the HI peak from the CO spectra. Hence, the velocity dispersion between the atomic and molecular components is extremely low, independently justifying the use of the HI line in building our maps. Stacking the spectra in concentric rings shows that the CO linewidth and possibly the CO-HI velocity dispersion decrease in the outer disk. Using the CO(2-1) emission to trace the molecular gas, the probability distribution function of the H2 column density shows an excess at high column density above a log normal distribution.

preprint2014arXiv

Why the Milky Way's bulge is not only a bar formed from a cold thin disk

By analyzing a N-body simulation of a bulge formed simply via a bar instability mechanism operating on a kinematically cold stellar disk, and by comparing the results of this analysis with the structural and kinematic properties of the main stellar populations of the Milky Way bulge, we conclude that the bulge of our Galaxy is not a pure stellar bar formed from a pre-existing thin stellar disk, as some studies have recently suggested. On the basis of several arguments emphasized in this paper, we propose that the bulge population which, in the Milky Way, is observed not to be part of the peanut structure corresponds to the old galactic thick disk, thus implying that the Milky Way is a pure thin+thick disk galaxy, with only a possible limited contribution of a classical bulge.

preprint2013arXiv

A gas-rich AGN near the centre of a galaxy cluster at z ~ 1.4

The formation of the first virialized structures in overdensities dates back to ~9 Gyr ago, i.e. in the redshift range z ~ 1.4 - 1.6. Some models of structure formation predict that the star formation activity in clusters was high at that epoch, implying large reservoirs of cold molecular gas. Aiming at finding a trace of this expected high molecular gas content in primeval clusters, we searched for the 12CO(2-1) line emission in the most luminous active galactic nucleus (AGN) of the cluster around the radio galaxy 7C 1756+6520 at z ~ 1.4, one of the farthest spectroscopic confirmed clusters. This AGN, called AGN.1317, is located in the neighbourhood of the central radio galaxy at a projected distance of ~780 kpc. The IRAM Plateau de Bure Interferometer was used to investigate the molecular gas quantity in AGN.1317, observing the 12CO(2-1) emission line. We detect CO emission in an AGN belonging to a galaxy cluster at z ~ 1.4. We measured a molecular gas mass of 1.1 x 10^10 Msun, comparable to that found in submillimeter galaxies. In optical images, AGN.1317 does not seem to be part of a galaxy interaction or merger.We also derived the nearly instantaneous star formation rate (SFR) from Halpha flux obtaining a SFR ~65 Msun/yr. This suggests that AGN.1317 is actively forming stars and will exhaust its reservoir of cold gas in ~0.2-1.0 Gyr.

preprint2013arXiv

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

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

preprint2013arXiv

A Ten Billion Solar Mass Outflow of Molecular Gas Launched by Radio Bubbles in the Abell 1835 Brightest Cluster Galaxy

We report ALMA Early Science observations of the Abell 1835 brightest cluster galaxy (BCG) in the CO (3-2) and CO (1-0) emission lines. We detect 5E10 solar masses of molecular gas within 10 kpc of the BCG. Its velocity width of ~130 km/s FWHM is too narrow to be supported by dynamical pressure. The gas may instead be supported in a rotating, turbulent disk oriented nearly face-on. The disk is forming stars at a rate of 100-180 solar masses per year. Roughly 1E10 solar masses of molecular gas is projected 3-10 kpc to the north-west and to the east of the nucleus with line of sight velocities lying between -250 km/s to +480 km/s with respect to the systemic velocity. Although inflow cannot be ruled out, the rising velocity gradient with radius is consistent with a broad, bipolar outflow driven by radio jets or buoyantly rising X-ray cavities. The molecular outflow may be associated with an outflow of hot gas in Abell 1835 seen on larger scales. Molecular gas is flowing out of the BCG at a rate of approximately 200 solar masses per year, which is comparable to its star formation rate. How radio bubbles lift dense molecular gas in their updrafts, how much gas will be lost to the BCG, and how much will return to fuel future star formation and AGN activity are poorly understood. Our results imply that radio-mechanical (radio mode) feedback not only heats hot atmospheres surrounding elliptical galaxies and BCGs, it is able to sweep higher density molecular gas away from their centers.

preprint2013arXiv

Acetone in Orion BN/KL - High-resolution maps of a special oxygen-bearing molecule

As one of the prime targets of interstellar chemistry study, Orion BN/KL clearly shows different molecular distributions between large nitrogen- (e.g., C2H5CN) and oxygen-bearing (e.g., HCOOCH3) molecules. However, acetone (CH3)2CO, a special complex O-bearing molecule, has been shown to have a very different distribution from other typical O-bearing molecules in the BN/KL region. We searched for acetone within our IRAM Plateau de Bure Interferometer 3 mm and 1.3 mm data sets. Twenty-two acetone lines were searched within these data sets. The angular resolution ranged from 1.8 X 0.8 to 6.0 X 2.3 arcsec^2, and the spectral resolution ranged from 0.4 to 1.9 km s-1. Nine of the acetone lines appear free of contamination. Three main acetone peaks (Ace-1, 2, and 3) are identified in Orion BN/KL. The new acetone source Ace-3 and the extended emission in the north of the hot core region have been found for the first time. An excitation temperature of about 150 K is determined toward Ace-1 and Ace-2, and the acetone column density is estimated to be 2-4 X 10^16 cm-2 with a relative abundance of 1-6 X 10^-8 toward these two peaks. Acetone is a few times less abundant toward the hot core and Ace-3 compared with Ace-1 and Ace-2. We find that the overall distribution of acetone in BN/KL is similar to that of N-bearing molecules, e.g., NH3 and C2H5CN, and very different from those of large O-bearing molecules, e.g., HCOOCH3 and (CH3)2O. Our findings show the acetone distribution is more extended than in previous studies and does not originate only in those areas where both N-bearing and O-bearing species are present. Moreover, because the N-bearing molecules may be associated with shocked gas in Orion BN/KL, this suggests that the formation and/or destruction of acetone may involve ammonia or large N-bearing molecules in a shocked-gas environment.

preprint2013arXiv

ALMA observations of cool dust in a low-metallicity starburst, SBS0335-052

We present Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 0 Band 7 observations of an extremely metal-poor dwarf starburst galaxy in the Local Universe, SBS0335-052 (12+log(O/H)~7.2). With these observations, dust is detected at 870micron (ALMA Band 7), but 87% of the flux in this band is due to free-free emission from the starburst. We have compiled a spectral energy distribution (SED) of SBS0335-052 that spans almost 6 orders of magnitude in wavelength and fit it with a spherical dust shell heated by a single-age stellar population; the best-fit model gives a dust mass of (3.8+/-0.6)x10^4 Msun. We have also constructed a SED including Herschel archival data for IZw18, another low-metallicity dwarf starburst (12+log(O/H)=7.17), and fit it with a similar model to obtain a dust mass of (3.4+/-1.0)x10^2 Msun. Compared with their atomic gas mass, the dust mass of SBS0335-052 far exceeds the prediction of a linear trend of dust-to-gas mass ratio with metallicity, while IZw18 falls far below. We use gas scaling relations to assess a putative missing gas component in both galaxies and find that the missing, possibly molecular, gas in SBS0335-052 is a factor of 6 times higher than the value inferred from the observed HI column density; in IZw18 the missing component is 4 times smaller. Ultimately, despite their similarly low metallicity, the differences in gas and dust column densities in SBS0335-052 and IZw18 suggest that metal abundance does not uniquely define star-formation processes. At some level, self-shielding and the survival of molecules may depend just as much on gas and dust column density as on metallicity. The effects of low metallicity may at least be partially compensated for by large column densities in the interstellar medium.

preprint2013arXiv

An extended Herschel drop-out source in the center of AS1063, a 'normal' dusty galaxy at z=6.1 or SZ substructures?

In the course of our 870um APEX/LABOCA follow up of the Herschel Lensing Survey we have detected a source in AS1063 (RXC J2248.7-4431), that has no counterparts in any of the Herschel PACS/SPIRE bands, it is a Herschel 'drop-out' with S_870/S_500>0.5. The 870um emission is extended and centered on the brightest cluster galaxy suggesting either a multiply imaged background source or substructure in the Sunyaev-Zel'dovich (SZ) increment due to inhomogeneities in the hot cluster gas of this merging cluster. We discuss both interpretations with emphasis on the putative lensed source. Based on the observed properties and on our lens model we find that this source could be the first SMG with a moderate far infrared luminosity (L_FIR<10^12 L_sol) detected so far at z>4. In deep HST observations we identified a multiply imaged z~6 source and we measured its spectroscopic redshift z=6.107 with VLT/FORS. This source could be associated with the putative SMG but it is most likely offset spatially by 10-30kpc and they could be interacting galaxies. With a FIR luminosity in the range [5-15]x10^{11} L_sol corresponding to a star formation rate in the range [80-260]M_sol/yr, this SMG would be more representative than the extreme starbursts usually detected at z>4. With a total magnification of ~25 it would open a unique window to the 'normal' dusty galaxies at the end of the epoch of reionization.

preprint2013arXiv

Cold dust in giant barred galaxy NGC1365

Observations of galaxies at sub-millimeter wavelengths, where the emission is mainly due to cold dust, are required to constrain the dust physical properties and provide important insight on the gas content of galaxies. We mapped NGC1365 at 870um with LABOCA, the Large APEX Bolometer Camera, allowing us to probe the central mass concentration as well as the rate at which the gas flows to the center. We obtained the dust physical properties both globally and locally for different locations in the galaxy. A 20 K modified black body represents about 98% of the total dust content of the galaxy, the rest can be represented by a warmer dust component of 40 K. The bar exhibits an east-west asymmetry in the dust distribution: The eastern bar is heavier than the western bar by more than a factor of 4. Integrating the dust SED, we derive a total infrared (IR) luminosity of 9.8 x 10^{10} L_{\odot} leading to a dust-enshrouded star formation rate of ~16.7 M_{\odot} per year in NGC1365. We derive the gas mass from the measurements of the dust emission resulting in a CO-to-H2 conversion factor of X_{CO}~1.2x10^{20} mol cm^{-2} (K km s^{-1})^{-1} in the central disk including the bar. Taking into account the metallicity variation, the central gas mass concentration is only ~20% at R<40" (3.6 kpc). On the other hand, the time-scale with which the gas flows into the center, ~300 Myr, is rather short. This indicates that the current central mass in NGC1365 is evolving fast due to the strong bar.

preprint2013arXiv

Cold gas dynamics in Hydra-A: evidence for a rotating disk

We present multi-frequency observations of the radio galaxy Hydra-A (3C218) located in the core of a massive, X-ray luminous galaxy cluster. IFU spectroscopy is used to trace the kinematics of the ionised and warm molecular hydrogen which are consistent with a ~ 5 kpc rotating disc. Broad, double-peaked lines of CO(2-1), [CII]157 $μ$m and [OI]63 $μ$m are detected. We estimate the mass of the cold gas within the disc to be M$_{gas}$ = 2.3 $\pm$ 0.3 x 10$^9$ M$_{\odot}$. These observations demonstrate that the complex line profiles found in the cold atomic and molecular gas are related to the rotating disc or ring of gas. Finally, an HST image of the galaxy shows that this gas disc contains a substantial mass of dust. The large gas mass, SFR and kinematics are consistent with the levels of gas cooling from the ICM. We conclude that the cold gas originates from the continual quiescent accumulation of cooled ICM gas. The rotation is in a plane perpendicular to the projected orientation of the radio jets and ICM cavities hinting at a possible connection between the kpc-scale cooling gas and the accretion of material onto the black hole. We discuss the implications of these observations for models of cold accretion, AGN feedback and cooling flows.

preprint2013arXiv

Cold gas in the inner regions of intermediate redshift clusters

Determining gas content and star formation rate has known remarkable progress in field galaxies, but has been much less investigated in galaxies inside clusters. We present the first CO observations of luminous infrared galaxies (LIRGs) inside the virial radii of two intermediate redshift clusters, CL1416+4446 (z=0.397) and CL0926+1242 (z=0.489). We detect three galaxies at high significance (5 to 10 sigma), and provide robust estimates of their CO luminosities, L'CO. In order to put our results into a general context, we revisit the relation between cold and hot gas and stellar mass in nearby field and cluster galaxies. We find evidence that at fixed LIR (or fixed stellar mass), the frequency of high L'CO galaxies is lower in clusters than in the field, suggesting environmental depletion of the reservoir of cold gas. The level of star formation activity in a galaxy is primarily linked to the amount of cold gas, rather than to the galaxy mass or the lookback time. In clusters, just as in the field, the conversion between gas and stars seems universal. The relation between LIR and L'CO for distant cluster galaxies extends the relation of nearby galaxies to higher IR luminosities. Nevertheless, the intermediate redshift galaxies fall well within the dispersion of the trend defined by local systems. Considering that L'CO is generally derived from the CO(1-0) line and sensitive to the vast majority of the molecular gas in the cold interstellar medium of galaxies, but less to the part which will actually be used to form stars, we suggest that molecular gas can be stripped before the star formation rate is affected. Combining the sample of Geach et al. (2009, 2011) and ours, we find evidence for a decrease in CO towards the cluster centers. This is the first hint of an environmental impact on cold gas at intermediate redshift.

preprint2013arXiv

Gas accretion in disk galaxies

Gas accretion is necessary to maintain star formation, spiral and bar structure, and secular evolution in galaxies. This can occur through tidal interaction, or mass accretion from cosmic filaments. Different processes will be reviewed to drive gas towards galaxy centers and trigger starbursts and AGN. The efficiency of these dynamical processes can be estimated through simulations and checked by observations at different redshift, across the Hubble time. Large progress has been made on galaxies at moderate and high redshifts, allowing to interpret the star formation history and star formation efficiency as a function of gas content, dynamical state and galaxy evolution.

preprint2013arXiv

Influence of baryonic physics in galaxy simulations: a semi-analytic treatment of the molecular component

Recent work in galaxy formation has enlightened the important role of baryon physics, to solve the main problems encountered by the standard theory at the galactic scale, such as the galaxy stellar mass functions, or the missing satellites problem. The present work aims at investigating in particular the role of the cold and dense molecular phase, which could play a role of gas reservoir in the outer galaxy discs, with low star formation efficiency. Through TreeSPH simulations, implementing the cooling to low temperatures, and the inclusion of the molecular hydrogen component, several feedback efficiencies are studied, and results on the gas morphology and star formation are obtained. It is shown that molecular hydrogen allows some slow star formation (with gas depletion times of about 5 Gyr) to occur in the outer parts of the discs. This dense and quiescent phase might be a way to store a significant fraction of dark baryons, in a relatively long time-scale, in the complete baryonic cycle, connecting the galaxy discs to hot gaseous haloes and to the cosmic filaments.

preprint2013arXiv

Massive molecular gas flows in the Abell 1664 brightest cluster galaxy

We report ALMA Early Science CO(1-0) and CO(3-2) observations of the brightest cluster galaxy (BCG) in Abell 1664. The BCG contains 1.1x10^{10} solar masses of molecular gas divided roughly equally between two distinct velocity systems: one from -250 to +250 km/s centred on the BCG's systemic velocity and a high velocity system blueshifted by 570 km/s with respect to the systemic velocity. The BCG's systemic component shows a smooth velocity gradient across the BCG center with velocity proportional to radius suggestive of solid body rotation about the nucleus. However, the mass and velocity structure are highly asymmetric and there is little star formation coincident with a putative disk. It may be an inflow of gas that will settle into a disk over several 10^8 yr. The high velocity system consists of two gas clumps, each ~2 kpc across, located to the north and southeast of the nucleus. Each has a line of sight velocity spread of 250-300 km/s. The velocity of the gas in the high velocity system tends to increase towards the BCG center and could signify a massive high velocity flow onto the nucleus. However, the velocity gradient is not smooth and these structures are also coincident with low optical-UV surface brightness regions, which could indicate dust extinction associated with each clump. If so, the high velocity gas would be projected in front of the BCG and moving toward us along the line of sight in a massive outflow most likely driven by the AGN. A merger origin is unlikely but cannot be ruled out.

preprint2013arXiv

PHIBSS: molecular gas, extinction, star formation and kinematics in the z=1.5 star forming galaxy EGS13011166

We report matched resolution, imaging spectroscopy of the CO J=3-2 line (with the IRAM Plateau de Bure millimeter interferometer) and of the H-alpha line (with LUCI at the Large Binocular Telescope)in the massive z=1.53 main-sequence galaxy EGS 13011166, as part of the "Plateau de Bure high-z, blue sequence survey (PHIBSS). We combine these data with HST V-J-J-H-band maps to derive spatially resolved distributions of stellar surface density, star formation rate, molecular gas surface density, optical extinction and gas kinematics. The spatial distribution and kinematics of the ionized and molecular gas are remarkably similar and are well modeled by a turbulent, globally Toomre unstable rotating disk. The stellar surface density distribution is smoother than the clumpy rest-frame UV/optical light distribution, and peaks in an obscured, star forming massive bulge near the dynamical center. The molecular gas surface density and the effective optical screen extinction track each other and are well modeled by a 'mixed' extinction model. The inferred slope of the spatially resolved molecular gas to star formation rate relation depends strongly on the adopted extinction model and can vary from 0.8 to 1.7. For the preferred mixed dust-gas model we find a near linear slope.

preprint2013arXiv

Polar ring galaxies as tests of gravity

Polar ring galaxies are ideal objects with which to study the three-dimensional shapes of galactic gravitational potentials since two rotation curves can be measured in two perpendicular planes. Observational studies have uncovered systematically larger rotation velocities in the extended polar rings than in the associated host galaxies. In the dark matter context, this can only be explained through dark halos that are systematically flattened along the polar rings. Here, we point out that these objects can also be used as very effective tests of gravity theories, such as those based on Milgromian dynamics (MOND). We run a set of polar ring models using both Milgromian and Newtonian dynamics to predict the expected shapes of the rotation curves in both planes, varying the total mass of the system, the mass of the ring with respect to the host, as well as the size of the hole at the center of the ring. We find that Milgromian dynamics not only naturally leads to rotation velocities being typically higher in the extended polar rings than in the hosts, as would be the case in Newtonian dynamics without dark matter, but that it also gets the shape and amplitude of velocities correct. Milgromian dynamics thus adequately explains this particular property of polar ring galaxies.

preprint2013arXiv

Probing the jet base of the blazar PKS1830-211 from the chromatic variability of its lensed images. Serendipitous ALMA observations of a strong gamma-ray flare

The launching mechanism of the jets of active galactic nuclei is observationally poorly constrained, due to the large distances to these objects and the very small scales (sub-parsec) involved. In order to better constrain theoretical models, it is especially important to get information from the region close to the physical base of the jet, where the plasma acceleration takes place. In this paper, we report multi-epoch and multi-frequency continuum observations of the z=2.5 blazar PKS1830-211 with ALMA, serendipitously coincident with a strong $γ$-ray flare reported by Fermi-LAT. The blazar is lensed by a foreground z=0.89 galaxy, with two bright images of the compact core separated by 1". Our ALMA observations individually resolve these two images (although not any of their substructures), and we study the change of their relative flux ratio with time (four epochs spread over nearly three times the time delay between the two lensed images) and frequency (between 350 and 1050 GHz, rest-frame of the blazar), during the $γ$-ray flare. In particular, we detect a remarkable frequency-dependent behaviour of the flux ratio, which implies the presence of a chromatic structure in the blazar (i.e., a core-shift effect). We rule out the possiblity of micro- and milli-lensing effects and propose instead a simple model of plasmon ejection in the blazar's jet to explain the time and frequency variability of the flux ratio. We suggest that PKS1830-211 is likely one of the best sources to probe the activity at the base of a blazar's jet at submillimeter wavelengths, thanks to the peculiar geometry of the system. The implications of the core-shift in absorption studies of the foreground z=0.89 galaxy (e.g., constraints on the cosmological variations of fundamental constants) are discussed.

preprint2013arXiv

Search for cold and hot gas in the ram pressure stripped Virgo dwarf galaxy IC3418

We present IRAM 30m sensitive upper limits on CO emission in the ram pressure stripped dwarf Virgo galaxy IC3418 and in a few positions covering HII regions in its prominent 17 kpc UV/Ha gas-stripped tail. In the central few arcseconds of the galaxy, we report a possible marginal detection of about 1x10^6 M_sun of molecular gas (assuming a Galactic CO-to-H_2 conversion factor) that could correspond to a surviving nuclear gas reservoir. We estimate that there is less molecular gas in the main body of IC3418, by at least a factor of 20, than would be expected from the pre-quenching UV-based star formation rate assuming the typical gas depletion timescale of 2 Gyr. Given the lack of star formation in the main body, we think the H_2-deficiency is real, although some of it may also arise from a higher CO-to-H_2 factor typical in low-metallicity, low-mass galaxies. The presence of HII regions in the tail of IC3418 suggests that there must be some dense gas; however, only upper limits of < 1x10^6 M_sun were found in the three observed points in the outer tail. This yields an upper limit on the molecular gas content of the whole tail < 1x10^7 M_sun, which is an amount similar to the estimates from the observed star formation rate over the tail. We also present strong upper limits on the X-ray emission of the stripped gas in IC3418 from a new Chandra observation. The measured X-ray luminosity of the IC3418 tail is about 280 times lower than that of ESO 137-001, a spiral galaxy in a more distant cluster with a prominent ram pressure stripped tail. Non-detection of any diffuse X-ray emission in the IC3418 tail may be due to a low gas content in the tail associated with its advanced evolutionary state and/or due to a rather low thermal pressure of the surrounding intra-cluster medium.

preprint2013arXiv

Signatures of radial migration in barred galaxies: Azimuthal variations in the metallicity distribution of old stars

By means of N-body simulations, we show that radial migration in galaxy disks, induced by bar and spiral arms, leads to significant azimuthal variations in the metallicity distribution of old stars at a given distance from the galaxy center. Metals do not show an axisymmetric distribution during phases of strong migration. Azimuthal variations are visible during the whole phase of strong bar phase, and tend to disappear as the effect of radial migration diminishes, together with a reduction in the bar strength. These results suggest that the presence of inhomogeneities in the metallicity distribution of old stars in a galaxy disk can be a probe of ongoing strong migration. Such signatures may be detected in the Milky Way by Gaia (and complementary spectroscopic data), as well as in external galaxies, by IFU surveys like CALIFA and ATLAS3D. Mixing - defined as the tendency toward a homogeneous, azimuthally symmetric, stellar distribution in the disk - and migration turns out to be two distinct processes, the effects of mixing starting to be visible when strong migration is over.

preprint2013arXiv

Spectral Energy Distributions of HII regions in M33 (HerM33es)

Within the framework of the Herschel M 33 extended survey HerM33es we study the Spectral Energy Distribution (SED) of a set of HII regions in M 33 as a function of the morphology. We present a catalogue of 119 HII regions morphologically classified: 9 filled, 47 mixed, 36 shell, and 27 clear shell HII regions. For each object we extract the photometry at twelve available wavelength bands (from FUV-1516A to IR-250mi) and obtain the SED. We also obtain emission line profiles across the regions to study the location of the stellar, ionised gas, and dust components. We find trends for the SEDs related to the morphology, showing that the star and gas-dust configuration affects the ratios of the emission in different bands. The mixed and filled regions show higher emission at 24mi than the shells and clear shells, which could be due to the proximity of the dust to the stellar clusters in the case of mixed and filled regions. The FIR peak for shells and clear shells seems to be located towards longer wavelengths, indicating that the dust is colder for this type of objects.The logarithmic 100/70mi ratio for filled and mixed regions remains constant over one order of magnitude in Halpha and FUV surface brightness, while the shells and clear shells exhibit a wider range of values of almost two orders of magnitude. We derive dust masses and temperatures fitting the individual SEDs with dust models proposed in the literature. The derived dust mass range is between 10^2-10^4 Msun and the cold dust temperature spans T(cold)~12-27 K. The spherical geometrical model proposed for the Halpha clear shells is confirmed by the emission profile obtained from the observations and is used to infer the electron density within the envelope: the typical electron density is 0.7+-0.3 cm^-3, while filled regions can reach values two to five times higher.

preprint2013arXiv

The Central 3 kpc of NGC 5850

NGC 5850 is a nearby (z=0.0085) early type spiral galaxy classified as LINER. It is considered as a prototype double-barred system. Our optical Integral Field Spectroscopic (IFS) data of the central 21x19 arcsec^2 of NGC 5850 show extended LINER-like emission which we ascribe to the presence of a hot and evolved stellar population, possibly together with a faint AGN. Additionally NGC 5850 shows extended `composite' ionization patterns, likely to stem from a mixture of LINER-like ionization and photoionization by star formation. The kinematics of the gas deviates strongly from a simple rotational structure.

preprint2013arXiv

Towards a resolved Kennicutt-Schmidt law at high redshift

Massive galaxies in the distant Universe form stars at much higher rates than today. Although direct resolution of the star forming regions of these galaxies is still a challenge, recent molecular gas observations at the IRAM Plateau de Bure interferometer enable us to study the star formation efficiency on subgalactic scales around redshift z = 1.2. We present a method for obtaining the gas and star formation rate (SFR) surface densities of ensembles of clumps composing galaxies at this redshift, even though the corresponding scales are not resolved. This method is based on identifying these structures in position-velocity diagrams corresponding to slices within the galaxies. We use unique IRAM observations of the CO(3-2) rotational line and DEEP2 spectra of four massive star forming distant galaxies - EGS13003805, EGS13004291, EGS12007881, and EGS13019128 in the AEGIS terminology - to determine the gas and SFR surface densities of the identifiable ensembles of clumps that constitute them. The integrated CO line luminosity is assumed to be directly proportional to the total gas mass, and the SFR is deduced from the [OII] line. We identify the ensembles of clumps with the angular resolution available in both CO and [OII] spectroscopy; i.e., 1-1.5". SFR and gas surface densities are averaged in areas of this size, which is also the thickness of the DEEP2 slits and of the extracted IRAM slices, and we derive a spatially resolved Kennicutt-Schmidt (KS) relation on a scale of ~8 kpc. The data generally indicates an average depletion time of 1.9 Gyr, but with significant variations from point to point within the galaxies.

preprint2013arXiv

Variation in the dust emissivity index across M33 with Herschel and Spitzer (HerM33es)

We study the wavelength dependence of the dust emission as a function of position and environment across the disk of M33 at a linear resolution of 160 pc using Spitzer and Herschel photometric data. Expressing the emissivity of the dust as a power law, the power-law exponent (beta) is estimated from two independent approaches designed to properly treat the degeneracy between beta and the dust temperature. Both beta and the dust temperature are higher in the inner disk than in the outer disk, contrary to reported beta-T anti-correlations found in other sources. In the cold + warm dust model, the warm component and the ionized gas (Halpha) have a very similar distribution across the galaxy, demonstrating that the model separates the components in an appropriate fashion. The flocculent spiral arms and the dust lanes are evident in the map of the cold component. Both cold and warm dust column densities are high in star forming regions and reach their maxima toward the giant star forming complexes NGC604 and NGC595. beta declines from close to 2 in the center to about 1.3 in the outer disk. beta is positively correlated with star formation and with molecular gas column, as traced by Halpha and CO emission. The lower dust emissivity index in the outer parts of M33 is likely related to the reduced metallicity (different grain composition) and possibly different size distribution. It is not due to the decrease in stellar radiation field or temperature in a simple way because the FIR-bright regions in the outer disk also have a low beta. Like most spirals, M33 has a (decreasing) radial gradient in star formation and molecular-to-atomic gas ratio such that the regions bright in Halpha or CO tend to trace the inner disk, making it difficult to distinguish between their effects on the dust.

preprint2012arXiv

A bright z=5.2 lensed submillimeter galaxy in the field of Abell 773: HLSJ091828.6+514223

During our Herschel Lensing Survey (HLS) of massive galaxy clusters, we have discovered an exceptionally bright source behind the z=0.22 cluster Abell 773, which appears to be a strongly lensed submillimeter galaxy (SMG) at z=5.2429. This source is unusual compared to most other lensed sources discovered by Herschel so far, because of its higher submm flux (\sim 200mJy at 500\micron) and its high redshift. The dominant lens is a foreground z=0.63 galaxy, not the cluster itself. The source has a far-infrared (FIR) luminosity of L_FIR= 1.1 10^{14}/μLo, where μis the magnification factor, likely \sim 11. We report here the redshift identification through CO lines with the IRAM-30m, and the analysis of the gas excitation, based on CO(7-6), CO(6-5), CO(5-4) detected at IRAM and the CO(2-1) at the EVLA. All lines decompose into a wide and strong red component, and a narrower and weaker blue component, 540\kms apart. Assuming the ultraluminous galaxy (ULIRG) CO-to-H2 conversion ratio, the H2 mass is 5.8 10^{11}/μMo, of which one third is in a cool component. From the CI line we derive a CI/H2 number abundance of 6 10^{-5} similar to that in other ULIRGs. The H2O line is strong only in the red velocity component, with an intensity ratio I(H_2O)/I(CO) \sim 0.5, suggesting a strong local FIR radiation field, possibly from an active nucleus (AGN) component. We detect the [NII]205\mics line for the first time at high-z. It shows comparable blue and red components, with a strikingly broad blue one, suggesting strong ionized gas flows.

preprint2012arXiv

A cold gas reservoir to fuel M31 nuclear black hole and stellar cluster

With IRAM-30m/HERA, we have detected CO(2-1) gas complexes within 30 arcsec (~100 pc) from the center of M31 amounting to a minimum total mass of 4.2 x 10^4 Msol (one third of the positions are detected). Averaging the whole HERA field, we have shown that there is no additional undetected diffuse component. We show that the above gas detection is associated with gas lying on the far side as no extinction is observed in the optical, but some emission is present on infra-red Spitzer maps. The kinematics is complex. (1) The velocity pattern is mainly redshifted: the dynamical center of the gas differs from the black hole position and the maximum of optical emission, and only the red-shifted side is seen in our data. (2) Several velocity components are detected in some lines of sight. Our interpretation is supported by the reanalysis of the effect of dust on a complete planetary nebula sample. Two dust components are detected with respective position angles of 37 deg and -66 deg. This is compatible with the superposition of the (PA=37 deg) disk dominated by the 10-kpc ring and the inner 0.7-kpc ring detected in infrared data, which position angle (-66 deg) can be measured for the first time. The large scale disk, which dominates the HI data, is very inclined (i=77 deg), warped and superposed on the line of sight on the less inclined inner ring. The detected CO emission might come from both components.

preprint2012arXiv

A precise and accurate determination of the cosmic microwave background temperature at z=0.89

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. We aim at measuring 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 PKS1830-211. We adopt a Monte-Carlo Markov Chain approach, coupled with predictions from the non-LTE radiative transfer code RADEX, to solve the excitation of a set of various molecular species directly from their spectra. We determine Tcmb=5.08 pm 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>0 to date.

preprint2012arXiv

An upper limit to the variation in the fundamental constants at redshift z = 5.2

Aims. We constrain a hypothetical variation in the fundamental physical constants over the course of cosmic time. Methods. We use unique observations of the CO(7-6) rotational line and the [CI] 3P_2 - 3P_1 fine-structure line towards a lensed galaxy at redshift z = 5.2 to constrain temporal variations in the constant F = alpha^2/mu, where mu is the electron-to-proton mass ratio and alpha is the fine-structure constant. The relative change in F between z = 0 and z = 5.2, dFF = (F_obs - F_lab)/F_lab, is estimated from the radial velocity offset, dV = V_rot - V_fs, between the rotational transitions in carbon monoxide and the fine-structure transition in atomic carbon. Results. We find a conservative value dV = 1 +/- 5 km/s (1sigma C.L.), which when interpreted in terms of dFF gives dFF < 2x10^-5. Independent methods restrict the mu-variations at the level of dmm < 1x10^-7 at z = 0.7 (look-back time t_z0.7 = 6.4 Gyr). Assuming that temporal variations in mu, if any, are linear, this leads to an upper limit on dmm < 2x10^-7 at z = 5.2 (t_z5.2 = 12.9 Gyr). From both constraints on dFF and dmm, one obtains for the relative change in alpha the estimate daa < 8x10^-6, which is at present the tightest limit on daa at early cosmological epochs.

preprint2012arXiv

CH3OCH3 in Orion-KL: a striking similarity with HCOOCH3

We used several data sets from the Plateau de Bure Interferometer to map the dimethyl ether emission in Orion-KL with different arcsec spatial resolutions and different energy levels to compare with our previous methyl formate maps. Our data show remarkable similarity between the dimethyl ether (CH3OCH3) and the methyl formate (HCOOCH3) distributions even on a small scale (1.8"x0.8" or about 500 AU). This long suspected similarity, seen from both observational and theoretical arguments, is demonstrated with unprecedented confidence, with a correlation coefficient of maps of 0.8. A common precursor is the simplest explanation of our correlation. Comparisons with previous laboratory work and chemical models suggest the major role of grain surface chemistry and a recent release, probably with little processing, of mantle molecules by shocks. In this case the CH3O radical produced from methanol ice would be the common precursor (whereas ethanol, C2H5OH, is produced from the radical CH2OH). The alternative gas phase scheme, where protonated methanol CH3OH2+ is the common precursor to produce methyl formate and dimethyl ether through reactions with HCOOH and CH3OH, is also compatible with our data. Our observations cannot yet definitely allow a choice between the different chemical processes, but the tight correlation between the distributions of HCOOCH3 and CH3OCH3 strongly contrasts with the different behavior we observe for the distributions of ethanol and formic acid. This provides a very significant constraint on models.

preprint2012arXiv

Cold and warm molecular gas in the outflow of 4C12.50

We present deep observations of the 12CO(1-0) and (3-2) lines in the ultra-luminous infrared and radio galaxy 4C12.50, carried out with the 30m telescope of the Institut de Radioastronomie Millimetrique. Our observations reveal the cold molecular gas component of a warm molecular gas outflow that was previously known from Spitzer Space Telescope data. The 12CO(3-2) profile indicates the presence of absorption at -950 km/s from systemic velocity with a central optical depth of 0.22. Its profile is similar to that of the HI absorption that was seen in radio data of this source. A potential detection of the (0-1) absorption enabled us to place an upper limit of 0.03 on its central optical depth, and to constrain the excitation temperature of the outflowing CO gas to >=65K assuming that the gas is thermalized. If the molecular clouds fully obscure the background millimeter continuum that is emitted by the radio core, the H2 column density is >=1.8*10^22 /cm^2. The outflow then carries an estimated cold H2 mass of at least 4.2*10^3 M_sun along the nuclear line of sight. This mass will be even higher when integrated over several lines of sight, but if it were to exceed 3*10^9 M_sun, the outflow would most likely be seen in emission. Since the ambient cold gas reservoir of 4C12.50 is 1.0*10^10 M_sun, the outflowing-to-ambient mass ratio of the warm gas (37%) could be elevated with respect to that of the cold gas.

preprint2012arXiv

Cool and warm dust emission from M33 (HerM33es)

We study the far-infrared emission from the nearby spiral galaxy M33 in order to investigate the dust physical properties such as the temperature and the luminosity density across the galaxy. Taking advantage of the unique wavelength coverage (100, 160, 250, 350 and 500 micron) of the Herschel Space Observatory and complementing our dataset with Spitzer-IRAC 5.8 and 8 micron and Spitzer-MIPS 24 and 70 micron data, we construct temperature and luminosity density maps by fitting two modified blackbodies of a fixed emissivity index of 1.5. We find that the 'cool' dust grains are heated at temperatures between 11 and 28 K with the lowest temperatures found in the outskirts of the galaxy and the highest ones in the center and in the bright HII regions. The infrared/submillimeter total luminosity (5 - 1000 micron) is estimated to be 1.9x10^9 Lsun. 59% of the total luminosity of the galaxy is produced by the 'cool' dust grains (~15 K) while the rest 41% is produced by 'warm' dust grains (~55 K). The ratio of the cool-to-warm dust luminosity is close to unity (within the computed uncertainties), throughout the galaxy, with the luminosity of the cool dust being slightly enhanced in the center of the galaxy. Decomposing the emission of the dust into two components (one emitted by the diffuse disk of the galaxy and one emitted by the spiral arms) we find that the fraction of the emission in the disk in the mid-infrared (24 micron) is 21%, while it gradually rises up to 57% in the submillimeter (500 micron). We find that the bulk of the luminosity comes from the spiral arm network that produces 70% of the total luminosity of the galaxy with the rest coming from the diffuse dust disk. The 'cool' dust inside the disk is heated at a narrow range of temperatures between 18 and 15 K (going from the center to the outer parts of the galaxy).

preprint2012arXiv

Dark matter distribution and its impact on the evolution of galaxy disks

The role of dark matter halos in galaxy disk evolution is reviewed, in particular the stabilisation of disks through self-gravity reduction, or the bar development through angular momentum exchange. Triaxial dark halos tend to weaken bars. When the dark mass inside the bar region is negligible, the bar develops through angular momentum exchange between inner and outer disk, and between stars and gas. Self-regulating cycles on the bar strength may develop in the presence of external gas accretion. Dynamical friction on dark halos slows down bars, which puts constraints on the dark matter amount inside the inner disk. During galaxy formation, baryons can lose most of their angular momentum if the infall is misaligned with the dark matter axes. Stable disks can form aligned with the minor axis of the dark halo. A sudden change in the infall direction, otherwise steady, can produce the peculiar polar ring galaxies. The dark matter halo can then be aligned along the polar disk. Misaligned infall can also maintain lopsidedness, which is only rarely due to galaxy interactions and mergers.

preprint2012arXiv

Dust and gas power-spectrum in M33 (HERM33ES)

Power spectra of de-projected images of late-type galaxies in gas and/or dust emission are very useful diagnostics of the dynamics and stability of their interstellar medium. Previous studies have shown that the power spectra can be approximated as two power-laws, a shallow one at large scales (larger than 500 pc) and a steeper one at small scales, with the break between the two corresponding to the line-of-sight thickness of the galaxy disk. We present a thorough analysis of the power spectra of the dust and gas emission at several wavelengths in the nearby galaxy M33. In particular, we use the recently obtained images at five wavelengths by PACS and SPIRE onboard Herschel. The large dynamical range (2-3 dex in scale) of most images allow us to determine clearly the change in slopes from -1.5 to -4, with some variations with wavelength. The break scale is increasing with wavelength, from 100 pc at 24 and 100micron to 350 pc at 500micron, suggesting that the cool dust lies in a thicker disk than the warm dust, may be due to star formation more confined to the plane. The slope at small scale tends to be steeper at longer wavelength, meaning that the warmer dust is more concentrated in clumps. Numerical simulations of an isolated late-type galaxy, rich in gas and with no bulge, like M33, are carried out, in order to better interpret these observed results. Varying the star formation and feedback parameters, it is possible to obtain a range of power-spectra, with two power-law slopes and breaks, which nicely bracket the data. The small-scale power-law is indeed reflecting the 3D behaviour of the gas layer, steepening strongly while the feedback smoothes the structures, by increasing the gas turbulence. M33 appears to correspond to a fiducial model with an SFR of $\sim$ 0.7 Mo/yr, with 10% supernovae energy coupled to the gas kinematics.

preprint2012arXiv

Evolution of Galactic Discs: Multiple Patterns, Radial Migration and Disc Outskirts

We investigate the evolution of galactic disks in N-body Tree-SPH simulations. We find that disks, initially truncated at three scale-lengths, can triple their radial extent, solely driven by secular evolution. Both Type I (single exponential) and Type II (down-turning) observed disk surface-brightness profiles can be explained by our findings. We relate these results to the strong angular momentum outward transfer, resulting from torques and radial migration associated with multiple patterns, such as central bars and spiral waves of different multiplicity. We show that even for stars ending up on cold orbits, the changes in angular momentum exhibit complex structure as a function of radius, unlike the expected effect of transient spirals alone. Focussing on one of our models, we find evidence for non-linear coupling among m=1, 2, 3 and 4 density waves, where m is the pattern multiplicity. We suggest that the naturally occurring larger resonance widths at galactic radii beyond four scale-lengths may have profound consequences on the formation and location of breaks in disk density profiles, provided spirals are present at such large distances. We also consider the effect of gas inflow and show that when in-plane smooth gas accretion of ~5 M_sun/yr is included, the outer disks become more unstable, leading to a strong increase in the stellar velocity dispersion. This, in turn, causes the formation of a Type III (up-turning) profile in the old stellar population. We propose that observations of Type III surface brightness profiles, combined with an up-turn in the stellar velocity dispersions beyond the disk break, could be a signature of ongoing gas-accretion. The results of this study suggest that disk outskirts comprised of stars migrated from the inner disk would have relatively large radial velocity dispersions, and significant thickness when seen edge-on. [Abridged]

preprint2012arXiv

Fueling the central engine of radio galaxies. II. The footprints of AGN feedback on the ISM of 3C 236

Aims: We study the emission of molecular gas in 3C236, a FR II radio source at z~0.1, and search for the footprints of AGN feedback. 3C236 shows signs of a reactivation of its AGN triggered by a recent minor merger episode. Observations have also previously identified an extreme HI outflow in this source. Methods: The IRAM PdBI has been used to study the distribution and kinematics of molecular gas in 3C236 by imaging with high spatial resolution the emission of the 12CO(2-1) line in the nucleus of the galaxy. We have searched for outflow signatures in the CO map. We have also derived the SFR in 3C236 using data available from the literature at UV, optical and IR wavelengths, to determine the star-formation efficiency of molecular gas. Results: The CO emission in 3C236 comes from a spatially resolved 2.6 kpc disk with a regular rotating pattern. Within the limits imposed by the sensitivity and velocity coverage of the CO data, we do not detect any outflow signatures in the cold molecular gas. The disk has a cold gas mass M(H2)~2.1x10^9 Msun. We determine a new value for the redshift of the source zCO=0.09927. The similarity between the CO and HI profiles indicates that the deep HI absorption in 3C236 can be accounted for by a rotating HI structure, restricting the evidence of HI outflow to the most extreme velocities. In the light of the new redshift, the analysis of the ionized gas kinematics reveals a 1000 km/s outflow. As for the CO emitting gas, outflow signatures are nevertheless absent in the warm molecular gas emission traced by infrared H2 lines. The star-formation efficiency in 3C236 is consistent with the value measured in normal galaxies, which follow the canonical KS-law. This result, confirmed to hold in other young radio sources examined in this work, is in stark contrast with the factor of 10-50 lower SFE that seems to characterize evolved powerful radio galaxies.

preprint2012arXiv

Gas fraction and star formation efficiency at z < 1.0

After new observations of 39 galaxies at z = 0.6-1.0 obtained at the IRAM 30m telescope, we present our full CO line survey covering the redshift range 0.2 < z < 1. Our aim is to determine the driving factors accounting for the steep decline in the star formation rate during this epoch. We study both the gas fraction, defined as Mgas/(Mgas+Mstar), and the star formation efficiency (SFE) defined by the ratio between far-infrared luminosity and molecular gas mass (LFIR/M(H2), i.e. a measure for the inverse of the gas depletion time. The sources are selected to be ultra-luminous infrared galaxies (ULIRGs), with LFIR greater than 10^12 Lo and experiencing starbursts. When we adopt a standard ULIRG CO-to-H2 conversion factor, their molecular gas depletion time is less than 100 Myr. Our full survey has now filled the gap of CO observations in the 0.2<z<1 range covering almost half of cosmic history. The detection rate in the 0.6 < z < 1 interval is 38% (15 galaxies out of 39), compared to 60% for the 0.2<z<0.6 interval. The average CO luminosity is L'CO = 1.8 10^10 K km/s pc^2, corresponding to an average H2 mass of 1.45 10^10 Mo. From observation of 7 galaxies in both CO(2-1) and CO(4-3), a high gas excitation has been derived; together with the dust mass estimation, this supports the choice of our low ULIRG conversion factor between CO luminosity and H2, for our sample sources. We find that both the gas fraction and the SFE significantly increase with redshift, by factors of 3 +-1 from z=0 to 1, and therefore both quantities play an important role and complement each other in cosmic star formation evolution.

preprint2012arXiv

Multiphase Signatures of AGN Feedback in Abell 2597

We present new Chandra X-ray observations of the brightest cluster galaxy (BCG) in the cool core cluster Abell 2597. The data reveal an extensive kpc-scale X-ray cavity network as well as a 15 kpc filament of soft-excess gas exhibiting strong spatial correlation with archival VLA radio data. In addition to several possible scenarios, multiwavelength evidence may suggest that the filament is associated with multiphase (10^3 - 10^7 K) gas that has been entrained and dredged-up by the propagating radio source. Stemming from a full spectral analysis, we also present profiles and 2D spectral maps of modeled X-ray temperature, entropy, pressure, and metal abundance. The maps reveal an arc of hot gas which in projection borders the inner edge of a large X-ray cavity. Although limited by strong caveats, we suggest that the hot arc may be (a) due to a compressed rim of cold gas pushed outward by the radio bubble or (b) morphologically and energetically consistent with cavity-driven active galactic nucleus (AGN) heating models invoked to quench cooling flows, in which the enthalpy of a buoyant X-ray cavity is locally thermalized as ambient gas rushes to refill its wake. If confirmed, this would be the first observational evidence for this model.

preprint2012arXiv

PHIBSS: molecular gas content and scaling relations in z~1-3 normal star forming galaxies

We present PHIBSS, the IRAM Plateau de Bure high-z blue sequence CO 3-2 survey of the molecular gas properties in normal star forming galaxies (SFGs) near the cosmic star formation peak. PHIBSS provides 52 CO detections in two redshift slices at z~1.2 and 2.2, with log(M*(M_solar))>10.4 and log(SFR(M_solar/yr))>1.5. Including a correction for the incomplete coverage of the M*-SFR plane, we infer average gas fractions of ~0.33 at z~1.2 and ~0.47 at z~2.2. Gas fractions drop with stellar mass, in agreement with cosmological simulations including strong star formation feedback. Most of the z~1-3 SFGs are rotationally supported turbulent disks. The sizes of CO and UV/optical emission are comparable. The molecular gas - star formation relation for the z=1-3 SFGs is near-linear, with a ~0.7 Gyrs gas depletion timescale; changes in depletion time are only a secondary effect. Since this timescale is much less than the Hubble time in all SFGs between z~0 and 2, fresh gas must be supplied with a fairly high duty cycle over several billion years. At given z and M*, gas fractions correlate strongly with the specific star formation rate. The variation of specific star formation rate between z~0 and 3 is mainly controlled by the fraction of baryonic mass that resides in cold gas.

preprint2012arXiv

Residual Cooling and Persistent Star Formation amid AGN Feedback in Abell 2597

New Chandra X-ray and Herschel FIR observations enable a multiwavelength study of active galactic nucleus (AGN) heating and intracluster medium (ICM) cooling in the brightest cluster galaxy of Abell 2597. The new Chandra observations reveal the central < 30 kiloparsec X-ray cavity network to be more extensive than previously thought, and associated with enough enthalpy to theoretically inhibit the inferred classical cooling flow. Nevertheless, we present new evidence, consistent with previous results, that a moderately strong residual cooling flow is persisting at 4%-8% of the classically predicted rates in a spatially structured manner amid the feedback-driven excavation of the X-ray cavity network. New Herschel observations are used to estimate warm and cold dust masses, a lower-limit gas-to-dust ratio, and a star formation rate consistent with previous measurements. The cooling time profile of the ambient X-ray atmosphere is used to map the locations of the observational star formation entropy threshold as well as the theoretical thermal instability threshold. Both lie just outside the < 30 kpc central region permeated by X-ray cavities, and star formation as well as ionized and molecular gas lie interior to both. The young stars are distributed in an elongated region that is aligned with the radio lobes, and their estimated ages are both younger and older than the X-ray cavity network, suggesting both jet-triggered as well as persistent star formation over the current AGN feedback episode. Bright X-ray knots that are coincident with extended Ly-alpha and FUV continuum filaments motivate a discussion of structured cooling from the ambient hot atmosphere along a projected axis that is perpendicular to X-ray cavity and radio axis. We conclude that the cooling ICM is the dominant contributor of the cold gas reservoir fueling star formation and AGN activity in the Abell 2597 BCG.

preprint2012arXiv

Star Formation Efficiency at Intermediate Redshift

Star formation is evolving very fast in the second half of the Universe, and it is yet unclear whether this is due to evolving gas content, or evolving star formation efficiency (SFE). We have carried out a survey of ultra-luminous galaxies (ULIRG) between z=0.2 and 1, to check the gas fraction in this domain of redshift which is still poorly known. Our survey with the IRAM-30m detected 33 galaxies out of 69, and we derive a significant evolution of both the gas fraction and SFE of ULIRGs over the whole period, and in particular a turning point around z=0.35. The result is sensitive to the CO-to-H2, conversion factor adopted, and both gas fraction and SFE have comparable evolution, when we adopt the low starburst conversion factor of α=0.8 Mo/(K km/s pc^2). Adopting a higher αwill increase the role of the gas fraction. Using α=0.8, the SFE and the gas fraction for z=0.2-1.0 ULIRGs are found to be significantly higher, by a factor 3, than for local ULIRGs, and are comparable to high redshift ones. We compare this evolution to the expected cosmic H2 abundance and the cosmic star formation history.

preprint2012arXiv

Survival of molecular gas in Virgo's hot intracluster medium: CO near M86

We carried out CO(1-0) and CO(2-1) observations of 21 different regions in the vicinity of M86, NGC4438, and along the 120 kpc-long, Ha-emitting filamentary trail that connects them, aiming to test whether molecular gas can survive to be transferred from a spiral to an elliptical galaxy in Virgo's 10^7K intracluster medium (ICM). We targeted Ha-emitting regions that could be associated with the interface between cold molecular clouds and the hot ionized ICM. The data, obtained with the 30m telescope of the Institut de Radioastronomie Millimetrique, led to the detection of molecular gas close to M86. CO gas with a recession velocity that is similar to that of the stars, -265 km/s, and with a corresponding H2 mass of 2*10^7 M_sun, was detected ~10 kpc southeast of the nucleus of M86, near the peak of its HI emission. We argue that it is possible for this molecular gas either to have formed in situ from HI, or to have been stripped from NGC4438 directly in molecular form. In situ formation is nonetheless negligible for the 7*10^6 M_sun of gas detected at 12:26:15.9+12:58:49, at ~10 kpc northeast of M86, where no (strong) HI emission is present. This detection provides evidence for the survival of molecular gas in filaments for timescales of ~100 Myr. An amount equivalent to 5*10^7 M_sun of H2 gas that could be lost to the ICM or to neighboring galaxies was also discovered in the tidal tail northwest of NGC4438. A scenario of gas being alternatively brought to M86 from NGC4388 on its south was also examined but it was considered unlikely due to the non detection of CO below or at the HI stream velocities, 2000-2700 km/s.

preprint2012arXiv

The bright-end of the luminosity function at z~9

We report new constraints on the galaxy luminosity function at z~9 based on observations carried out with ESO/VLT FORS2, HAWK-I and X-Shooter around the lensing cluster A2667, as part of our project aimed at selecting z~7-10 candidates accessible to spectroscopy. Only one J-dropout source was selected in this field fulfilling the color and magnitude criteria. This source was recently confirmed as a mid-z interloper based on X-Shooter spectroscopy. The depth and the area covered by our survey are well suited to set strong constraints on the bright-end of the galaxy luminosity function and hence on the star formation history at very high redshift. The non-detection of reliable J-dropout sources over the ~36arcmin2 field of view towards A2667 was used to carefully determine the lens-corrected effective volume and the corresponding upper-limit on the density of sources. The strongest limit is obtained for Phi(M_{1500}=-21.4+/-0.50)<6.70x10^{-6}Mpc^{-3}mag^{-1} at z~9. A maximum-likelihood fit of the luminosity function using all available data points including the present new result yields M*>-19.7 with fixed alpha=-1.74 and Phi*=1.10x10^{-3}Mpc^{-3}. The corresponding star formation rate density should be rho_{SFR}<5.97x10^{-3}M_{solar}/yr/Mpc^{3} at z~9. These results are in good agreement with the most recent estimates already published in this range of redshift and for this luminosity domain. This new result confirms the decrease in the density of luminous galaxies at very high-redshift, hence providing strong constraints for the design of future surveys aiming to explore the very high-redshift Universe.

preprint2011arXiv

A very extended molecular web around NGC 1275

We present the first detection of CO emission lines in the Halpha filaments at distances as far as 50 kpc from the centre of the galaxy NGC 1275. This gas is probably dense (>=10E3 cm-3). However, it is not possible to accurately determine the density and the kinetic temperature of this relatively warm gas (Tkin~20-500K) with the current data only. The amount of molecular gas in the filaments is large 10E9 Msun (assuming a Galactic N(H2)/Ico ratio). This is 10% of the total mass of molecular gas detected in this cD galaxy. This gas has large-scale velocities comparable to those seen in Halpha. The origin of the filaments is still unclear, but their formation is very likely linked to the AGN positive feedback (Revaz et al., 2008) that regulates the cooling of the surrounding X-ray-emitting gas as suggested by numerical simulations. We also present high-resolution spectra of the galaxy core. The spatial characteristics of the double-peaked profile suggest that the molecular web of filaments and streamers penetrates down to radii of less than 2 kpc from the central AGN and eventually feed the galaxy nucleus. The mass of gas inside the very central region is ~10E^9 Msun, and is similar to the mass of molecular gas found in the filaments.

preprint2011arXiv

A view of the narrow-line region in the infrared: active galactic nuclei with resolved fine-structure lines in the Spitzer archive

We queried the Spitzer archive for high-resolution observations with the Infrared Spectrograph of optically selected active galactic nuclei (AGN) for the purpose of identifying sources with resolved fine-structure lines that would enable studies of the narrow-line region (NLR) at mid-infrared wavelengths. By combining 298 Spitzer spectra with 6 Infrared Space Observatory spectra, we present kinematic information of the NLR for 81 z<=0.3 AGN. We used the [NeV], [OIV], [NeIII], and [SIV] lines, whose fluxes correlate well with each other, to probe gas photoionized by the AGN. We found that the widths of the lines are, on average, increasing with the ionization potential of the species that emit them. No correlation of the line width with the critical density of the corresponding transition was found. The velocity dispersion of the gas, sigma, is systematically higher than that of the stars, sigma_*, in the AGN host galaxy, and it scales with the mass of the central black hole, M_BH. Further correlations between the line widths and luminosities L, and between L and M_BH, are suggestive of a three dimensional plane connecting log(M_BH) to a linear combination of log(sigma) and log(L). Such a plane can be understood within the context of gas motions that are driven by AGN feedback mechanisms, or virialized gas motions with a power-law dependence of the NLR radius on the AGN luminosity. The M_BH estimates obtained for 35 type 2 AGN from this plane are consistent with those obtained from the M_BH-sigma_* relation.

preprint2011arXiv

Dust heating sources in galaxies: the case of M33 (HERM33ES)

Dust emission is one of the main windows to the physics of galaxies and to star formation as the radiation from young, hot stars is absorbed by the dust and reemitted at longer wavelengths. The recently launched Herschel satellite now provides a view of dust emission in the far-infrared at an unequaled resolution and quality up to 500 μm. In the context of the Herschel HERM33ES open time key project, we are studying the moderately inclined Scd local group galaxy M33 which is located only 840 kpc away. In this article, using Spitzer and Herschel data ranging from 3.6 μm to 500 μm, along with HI, Hα maps, and GALEX ultraviolet data we have studied the emission of the dust at the high spatial resolution of 150 pc. Combining Spitzer and Herschel bands, we have provided new, inclination corrected, resolved estimators of the total infrared brightness and of the star formation rate from any combination of these bands. The study of the colors of the warm and cold dust populations shows that the temperature of the former is, at high brightness, dictated by young massive stars but, at lower brightness, heating is taken over by the evolved populations. Conversely, the temperature of the cold dust is tightly driven by the evolved stellar populations.

preprint2011arXiv

Far infrared constraints on the contamination by dust obscured galaxies of high-z dropout searches

The spectral energy distributions (SED) of dusty galaxies at intermediate redshift may look similar to very high redshift galaxies in the optical/near infrared (NIR) domain. This can lead to the contamination of high redshift galaxy searches based on broad band optical/NIR photometry by lower redshift dusty galaxies as both kind of galaxies cannot be distinguished. The contamination rate could be as high as 50%. {This work shows how the far infrared (FIR) domain can help to recognize likely low-z interlopers in an optical/NIR search for high-z galaxies.} We analyse the FIR SEDs of two galaxies proposed as very high redshift ($z>7$) dropout candidates based on deep Hawk-I/VLT observations. The FIR SEDs are sampled with PACS/Herschel at 100 and 160\,$μ$m, with SPIRE/Herschel at 250, 350 and 500\,$μ$m and with LABOCA/APEX at 870\,$μ$m. We find that redshifts $>7$ would imply extreme FIR SEDs (with dust temperatures $>100$\,K and FIR luminosities $>10^{13}$\,$L_{\odot}$). At z$\sim$2, instead, the SEDs of both sources would be compatible with that of typical ULIRGs/SMGs. Considering all the data available for these sources from visible to FIR we re-estimate the redshifts and we find $z\sim$1.6--2.5. Due to the strong spectral breaks observed in these galaxies, standard templates from the literature fail to reproduce the visible-near IR part of the SEDs even when additional extinction is included. These sources resemble strongly dust obscured galaxies selected in Spitzer observations with extreme visible-to-FIR colors, and the galaxy GN10 at $z=4$. Galaxies with similar SEDs could contaminate other high redshift surveys.

preprint2011arXiv

Galaxy Evolution and Star Formation Efficiency at 0.2 < z < 0.6

We present the results of a CO line survey of 30 galaxies at moderate redshift (z \sim 0.2-0.6), with the IRAM 30m telescope, with the goal to follow galaxy evolution and in particular the star formation efficiency (SFE) as defined by the ratio between far-infrared luminosity and molecular gas mass (LFIR/M(H2)). The sources are selected to be ultra-luminous infrared galaxies (ULIRGs), with LFIR larger than 2.8 10^{12} Lsol, experiencing starbursts; adopting a low ULIRG CO-to-H2 conversion factor, their gas consumption time-scale is lower than 10^8 yr. To date only very few CO observations exist in this redshift range that spans nearly 25% of the universe's age. Considerable evolution of the star formation rate is already observed during this period. 18 galaxies out of our sample of 30 are detected (of which 16 are new detections), corresponding to a detection rate of 60%. The average CO luminosity for the 18 galaxies detected is L'CO = 2 10^{10} K km/s pc^2, corresponding to an average H2 mass of 1.6 10^{10} Msol. The FIR luminosity correlates well with the CO luminosity, in agreement with the correlation found for low and high redshift ULIRGs. Although the conversion factor between CO luminosity and H2 mass is uncertain, we find that the maximum amount of gas available for a single galaxy is quickly increasing as a function of redshift. Using the same conversion factor, the SFEs for z\sim 0.2-0.6 ULIRGs are found to be significantly higher, by a factor 3, than for local ULIRGs, and are comparable to high redshift ones. We compare this evolution to the expected cosmic H2 abundance and the cosmic star formation history.

preprint2011arXiv

Gas inflows, star formation and metallicity evolution in galaxy pairs

It has been known since many decades that galaxy interactions can induce star formation (hereafter SF) enhancements and that one of the driving mechanisms of this enhancement is related to gas inflows into the central galaxy regions, induced by asymmetries in the stellar component, like bars. In the last years many evidences have been accumulating, showing that interacting pairs have central gas-phase metallicities lower than those of field galaxies, by {\sim} 0.2-0.3 dex on average. These diluted ISM metallicities have been explained as the result of inflows of metal-poor gas from the outer disk to the galaxy central regions. A number of questions arises: What's the timing and the duration of this dilution? How and when does the SF induced by the gas inflow enrich the circumnuclear gas with re-processed material? Is there any correlation between the timing and strength of the dilution and the timing and intensity of the SF? By means of Tree-SPH simulations of galaxy major interactions, we have studied the effect that gas inflows have on the ISM dilution, and the effect that the induced SF has, subsequently, in re-enriching the nuclear gas. In this contribution, we present the main results of this study.

preprint2011arXiv

Giant Molecular Clouds in the Local Group Galaxy M33

We present an analysis of the systematic CO(2-1) survey at 12" resolution covering most of the local group spiral M 33 which, at a distance of 840 kpc, is close enough that individual giant molecular clouds (GMCs) can be identified. The goal of this work is to study the properties of the GMCs in this subsolar metallicity galaxy. The CPROPS (Cloud Properties) algorithm (Rosolowsky & Leroy 2006) was used to identify 337 GMCs in M 33, the largest sample to date in an external galaxy. The sample is used to study the GMC luminosity function, or mass spectrum under the assumption of a constant N(H2)/ICO ratio. We find that n(L)dL = K*L^(-2.0\pm0.1) for the entire sample. However, when the sample is divided into inner and outer disk samples, the exponent changes from 1.6 \pm 0.2 for the centre 2 kpc to 2.3 \pm 0.2 for galactocentric distances larger than 2 kpc. Based on the emission in the FUV, Halpha, 8mu, and 24mu bands, each cloud was classified in terms of its star forming activity - no star formation, embedded, or exposed star formation (visible in FUV and Halpha). At least one sixth of the clouds had no (massive) star formation, suggesting that the average time required for star formation to start is about one sixth of the total time for which the object is identifiable as a GMC. The clouds without star formation have significantly lower CO luminosities than those with star formation, whether embedded or exposed, presumably related to the lack of heating sources. Taking the cloud sample as a whole, the main non-trivial correlation was the decrease in cloud CO brightness (or luminosity) with galactocentric radius. The complete cloud catalog, including CO and HI spectra and the CO contours on the FUV, Halpha, 8mu, and 24mu images is presented in the appendix.

preprint2011arXiv

HCOOCH3 as a probe of temperature and structure of Orion-KL

We studied the O-bearing molecule HCOOCH3 to characterize the physical conditions of the different molecular source components in Orion-KL. We identify 28 methyl formate emission peaks throughout the 50" field of observations. The two strongest peaks are in the Compact Ridge (MF1) and in the SouthWest of the Hot Core (MF2). Spectral confusion is still prevailing as half of the expected transitions are blended over the region. Assuming that the transitions are thermalized, we derive the temperature at the five main emission peaks. At the MF1 position we find a temperature of 80K in a 1.8"x0.8" beam size and 120K on a larger scale (3.6" x2.2"), suggesting an external source of heating, whereas the temperature is about 130K at the MF2 position on both scales. Transitions of HCOOCH3 in vt=1 are detected as well and the good agreement of the positions on the rotational diagrams between the vt=0 and the vt=1 transitions suggests a similar temperature. The velocity of the gas is between 7.5 and 8.0km/s depending on the positions and column density peaks vary from 1.6x10^16 to 1.6x10^17cm^-2. A second velocity component is observed around 9-10 km/s in a North-South structure stretching from the Compact Ridge up to the BN object; this component is warmer at the MF1 peak. The two other C2H4O2 isomers are not detected and the derived upper limit for the column density is <3x10^14cm^-2 for glycolaldehyde and <2x10^15cm^-2 for acetic acid. From the 223GHz continuum map, we identify several dust clumps with associated gas masses in the range 0.8 to 5.8Msun. Assuming that the HCOOCH3 is spatially distributed as the dust, we find relative abundances of HCOOCH3 in the range <0.1x10^-8 to 5.2x10^-8. We suggest a relation between the methyl formate distribution and shocks as traced by 2.12 mum H2 emission.

preprint2011arXiv

Molecular Gas and Star Formation in Local Early-Type Galaxies

The molecular gas content of local early-type galaxies is constrained and discussed in relation to their evolution. First, as part of the Atlas3D survey, we present the first complete, large (260 objects), volume-limited single-dish survey of CO in normal local early-type galaxies. We find a surprisingly high detection rate of 22%, independent of luminosity and at best weakly dependent on environment. Second, the extent of the molecular gas is constrained with CO synthesis imaging, and a variety of morphologies is revealed. The kinematics of the molecular gas and stars are often misaligned, implying an external gas origin in over a third of the systems, although this behaviour is drastically diffferent between field and cluster environments. Third, many objects appear to be in the process of forming regular kpc-size decoupled disks, and a star formation sequence can be sketched by piecing together multi-wavelength information on the molecular gas, current star formation, and young stars. Last, early-type galaxies do not seem to systematically obey all our usual prejudices regarding star formation, following the standard Schmidt-Kennicutt law but not the far infrared-radio correlation. This may suggest a greater diversity in star formation processes than observed in disk galaxies. Using multiple molecular tracers, we are thus starting to probe the physical conditions of the cold gas in early-types.

preprint2011arXiv

Molecular gas in galaxies at all redshifts

I review some recent results about the molecular content of galaxies, obtained essentially from the CO lines, but also dense tracers, or the dust continuum emission. New results have been obtained on molecular cloud physics, and their efficiency to form stars, shedding light on the Kennicutt-Schmidt law as a function of surface density and galaxy type. Large progress has been made on galaxy at moderate and high redshifts, allowing to interprete the star formation history and star formation efficiency as a function of gas content, or galaxy evolution. In massive galaxies, the gas fraction was higher in the past, and galaxy disks were more unstable and more turbulent. ALMA observations will allow the study of more normal galaxies at high z with higher spatial resolution and sensitivity.

preprint2011arXiv

Molecular Gas in NUclei of GAlaxies (NUGA) XIV. The barred LINER/Seyfert 2 galaxy NGC 3627

We present CO(1-0) and CO(2-1) maps of the interacting barred LINER/Seyfert 2 galaxy NGC 3627 obtained with the IRAM interferometer at resolutions of 2.1" x 1.3" and 0.9" x 0.6", respectively. The molecular gas emission shows a nuclear peak, an elongated bar-like structure of ~18" (~900 pc) diameter in both CO maps and, in CO(1-0), a two-arm spiral feature from r~9" (~450 pc) to r~16" (~800 pc). The inner ~18" bar-like structure, with a north/south orientation (PA = 14°), forms two peaks at the extremes of this elongated emission region. The kinematics of the inner molecular gas shows signatures of non-circular motions associated both with the 18" bar-like structure and the spiral feature detected beyond it. The 1.6 micron H-band 2MASS image of NGC 3627 shows a stellar bar with a PA = -21°, different from the PA (= 14°) of the CO bar-like structure, indicating that the gas is leading the stellar bar. The torques computed with the HST-NICMOS F160W image and our PdBI maps are negative down to the resolution limit of our images, ~60 pc in CO(2-1). If the bar ends at ~3 kpc, coincident with corotation (CR), the torques are negative between the CR of the bar and the nucleus, down to the resolution limit of our observations. This scenario is compatible with a recently-formed rapidly rotating bar which has had insufficient time to slow down because of secular evolution, and thus has not yet formed an inner Lindblad resonance (ILR). The presence of molecular gas inside the CR of the primary bar, where we expect that the ILR will form, makes NGC 3627 a potential smoking gun of inner gas inflow. The gas is fueling the central region, and in a second step could fuel directly the active nucleus.

preprint2011arXiv

Molecular gas in NUclei of GAlaxies (NUGA) XV. Molecular gas kinematics in the inner 3kpc of NGC6951

Within the NUclei of GAlaxies project we have obtained IRAM PdBI and 30m 12CO(1-0) and 12CO(2-1) observations of the spiral galaxy NGC 6951. Previous work shows that there is indirect evidence of gas inflow from 3 kpc down to small radii: a large-scale stellar bar, a prominent starburst ring (r~580 pc) and a LINER/Seyfert 2 nucleus. In this paper we study the gas kinematics as traced by the CO line emission in detail. We quantify the influence of the large-scale stellar bar by constructing an analytical model of the evolution of gas particles in a barred potential. From this model gravitational torques and mass accumulation rates are computed. We compare our model-based gravitational torque results with previous observationally-based ones. The model also shows that the large-scale stellar bar is indeed the dominant force for driving the gas inward, to the starburst ring. Inside the ring itself a nuclear stellar oval might play an important role. Detailed analysis of the CO gas kinematics there shows that emission arises from two co-spatial, but kinematically distinct components at several locations. The main emission component can always be related to the overall bar-driven gas kinematics. The second component exhibits velocities that are larger than expected for gas on stable orbits, has a molecular gas mass of 1.8x10^6Msun, is very likely connected to the nuclear stellar oval, and is consistent with inflowing motion towards the very center. This may form the last link in the chain of gas inflow towards the active galactic nucleus in NGC 6951.

preprint2011arXiv

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. [...] 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. [...] We also report the discovery of several new absorption components, with velocities spanning between -300 and +170 km/s. 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 |Delta_mu/mu|<4e-6, which takes into account the kinematical noise produced by the velocity dispersion measured from a large number of molecular species.

preprint2011arXiv

On the absence of molecular absorption in high redshift millimetre-band searches

We have undertaken a search for millimetre-wave band absorption (through the CO and HCO+ rotational transitions) in the host galaxies of reddened radio sources (z = 0.405-1.802). Despite the colour selection (optical-near infrared colours of V - K > 5 in all but one source), no absorption was found in any of the eight quasars for which the background continuum flux was detected. On the basis of the previous (mostly intervening) H2 and OH detections, the limits reached here and in some previous surveys should be deep enough to detect molecular absorption according to their V - K colours. However, our survey makes the assumption that the reddening is associated with dust close to the emission redshift of the quasar and that the narrow millimetre component of this emission is intercepted by the compact molecular cores. By using the known millimetre absorbers to define the colour depth and comparing this with the ultra-violet luminosities of the sources, we find that, even if these assumptions are valid, only twelve of the forty objects (mainly from this work) are potentially detectable. This is assuming an excitation temperature of 10 K at z=0, with the number decreasing with increasing temperatures (to zero detectable at 100 K).

preprint2011arXiv

Radial migration in galactic disks caused by resonance overlap of multiple patterns: Self-consistent simulations

We have recently identified a new radial migration mechanism resulting from the overlap of spiral and bar resonances in galactic disks. Here we confirm the efficiency of this mechanism in fully self-consistent, Tree-SPH simulations, as well as high-resolution pure N-body simulations. In all barred cases we clearly identify the effect of spiral-bar resonance overlap by measuring a bimodality in the changes of angular momentum in the disk, dL, whose maxima are near the bar's corotation and outer Lindblad resonance. This contrasts with the smooth distribution of dL for a simulation with no stable bar present, where strong radial migration is induced by multiple spirals. The presence of a disk gaseous component appears to increase the rate of angular momentum exchange by about 20%. The efficiency of this mechanism is such that galactic stellar disks can extend to over 10 scale-lengths within 1-3 Gyr in both Milky Way size and low-mass galaxies (circular velocity ~100 km/s). We also show that metallicity gradients can flatten in less than 1 Gyr rendering mixing in barred galaxies an order of magnitude more efficient than previously thought.

preprint2011arXiv

SMA/PdBI multiple line observations of the nearby Seyfert2 galaxy NGC 1068: Shock related gas kinematics and heating in the central 100pc?

We present high angular resolution (0.5-2.0") observations of the mm continuum and the 12CO(J=3-2), 13CO(J=3-2), 13CO(J=2-1), C18O(J=2-1), HCN(J=3-2), HCO+(J=4-3) and HCO+(J=3-2) line emission in the circumnuclear disk (r=100pc) of the proto-typical Seyfert type-2 galaxy NGC1068, carried out with the Submillimeter Array. We further include in our analysis new 13CO(J=1-0) and improved 12CO(J=2-1) observations of NGC1068 at high angular resolution (1.0-2.0") and sensitivity, conducted with the IRAM Plateau de Bure Interferometer. Based on the complex dynamics of the molecular gas emission indicating non-circular motions in the central ~100pc, we propose a scenario in which part of the molecular gas in the circumnuclear disk of NGC1068 is radially blown outwards as a result of shocks. This shock scenario is further supported by quite warm (Tkin>=200K) and dense (nH2=10^4cm^-3) gas constrained from the observed molecular line ratios. The HCN abundance in the circumnuclear disk is found to be [HCN]/[12CO]=10^-3.5. This is slightly higher than the abundances derived for galactic and extragalactic starforming/starbursting regions. This results lends further support to X-ray enhanced HCN formation in the circumnuclear disk of NGC1068, as suggested by earlier studies. The HCO+ abundance ([HCO+]/[12CO]=10^-5) appears to be somewhat lower than that of galactic and extragalactic starforming/starbursting regions. When trying to fit the cm to mm continuum emission by different thermal and non-thermal processes, it appears that electron-scattered synchrotron emission yields the best results while thermal free-free emission seems to over-predict the mm continuum emission.

preprint2011arXiv

The Herschel M33 extended survey (HerM33es): PACS spectroscopy of the star forming region BCLMP 302

Context: The emission line of [CII] at 158 micron is one of the strongest cooling lines of the interstellar medium (ISM) in galaxies. Aims: Disentangling the relative contributions of the different ISM phases to [CII] emission, is a major topic of the HerM33es program, a Herschel key project to study the ISM in the nearby spiral galaxy M33. Methods: Using PACS, we have mapped the emission of [CII] 158 micron, [OI] 63 micron, and other FIR lines in a 2'x2' region of the northern spiral arm of M33, centered on the HII region BCLMP302. At the peak of H-alpha emission, we have observed in addition a velocity resolved [CII] spectrum using HIFI. We use scatterplots to compare these data with PACS 160 micron continuum maps, and with maps of CO and HI data, at a common resolution of 12 arcsec or 50 pc. Maps of H-alpha and 24 micron emission observed with Spitzer are used to estimate the SFR. We have created maps of the [CII] and [OI] 63 micron emission and detected [NII] 122 micron and NIII 57 micron at individual positions. Results: The [CII] line observed with HIFI is significantly broader than that of CO, and slightly blue-shifted. In addition, there is little spatial correlation between [CII] observed with PACS and CO over the mapped region. There is even less spatial correlation between [CII] and the atomic gas traced by HI. Detailed comparison of the observed intensities towards the HII region with models of photo ionization and photon dominated regions, confirms that a significant fraction, 20--30%, of the observed [CII] emission stems from the ionized gas and not from the molecular cloud. The gas heating efficiency, using the ratio between [CII] and the TIR as a proxy, varies between 0.07 and 1.5%, with the largest variations found outside the HII region.

preprint2011arXiv

The metallicity dependence of the CO {\rightarrow} H_2 conversion factor in z>1 star forming galaxies

We use the first systematic samples of CO millimeter emission in z>1 'main-sequence' star forming galaxies (SFGs) to study the metallicity dependence of the conversion factor αCO, from CO line luminosity to molecular gas mass. The molecular gas depletion rate inferred from the ratio of the star formation rate (SFR) to CO luminosity, is ~1 Gyr-1 for near-solar metallicity galaxies with stellar masses above M_S~1e11 M_sun. In this regime the depletion rate does not vary more than a factor of two to three as a function of molecular gas surface density, or redshift between z~0 and 2. Below M_S the depletion rate increases rapidly with decreasing metallicity. We argue that this trend is not caused by starburst events, by changes in the physical parameters of the molecular clouds, or by the impact of the fundamental metallicity-SFR-stellar mass relation. A more probable explanation is that the conversion factor is metallicity dependent and that star formation can occur in 'CO-dark' gas. The trend is also expected theoretically from the effect of enhanced photodissociation of CO by ultraviolet radiation at low metallicity. From the available z~0 and z~1-3 samples we constrain the slope of the log(αCO) -log (metallicity) relation to range between -1 and -2, fairly insensitive to the assumed slope of the gas-star formation rate relation. Because of the lower metallicities near the peak of the galaxy formation activity at z~1-2 compared to z~0, we suggest that molecular gas masses estimated from CO luminosities have to be substantially corrected upward for galaxies below M_S.

preprint2011arXiv

Turbulent and fast motions of H2 gas in active galactic nuclei

Querying the Spitzer archive for optically-selected active galactic nuclei (AGN) observed in high-resolution mode spectroscopy, we identified radio and/or interacting galaxies with highly turbulent motions of the H2 gas at a temperature of a few hundred Kelvin. Unlike all other AGN that have unresolved H2 line profiles at a spectral resolution of ~600, 3C236, 3C293, IRAS09039+0503, MCG-2-58-22, and Mrk463E have intrinsic velocity dispersions exceeding 200 km/s for at least two of the rotational S0, S1, S2, and S3 lines. In a sixth source, 4C12.50, a blue wing was detected in the S1 and S2 line profiles, indicating the presence of a warm molecular gas component moving at -640 km/s with respect to the bulk of the gas at systemic velocity. Its mass is 5.2*10^7 M_sun, accounting for more than one fourth of the H2 gas at 374K, but less than 1% of the cold H2 gas computed from CO observations. Because no diffuse gas component of 4C12.50 has been observed to date to be moving at more than 250 km/s from systemic velocity, the H2 line wings are unlikely to be tracing gas in shock regions along the tidal tails of this merging system. They can instead be tracing gas driven by a jet or entrained by a nuclear outflow, which is known to emerge from the west nucleus of 4C12.50. It is improbable that such an outflow, with an estimated mass loss rate of 130 M_sun/yr, entirely quenches the star formation around this nucleus.

preprint2011arXiv

Variation in the dust spectral index across M33

Using the Herschel PACS and SPIRE FIR/submm data, we investigate variations in the dust spectral index $β$ in the nearby spiral galaxy M33 at a linear resolution of 160 pc. We use an iteration method in two different approaches, single and two-component modified black body models. In both approaches, $β$ is higher in the central disk than in the outer disk similar to the dust temperature. There is a positive correlation between $β$ and Halpha as well as with the molecular gas traced by CO(2-1). A Monte-Carlo simulation shows that the physical parameters are better constrained when using the two-component model.

preprint2010arXiv

100 mum and 160 mum emission as resolved star-formation rate estimators in M33 (HERM33ES)

Over the past few years several studies have provided estimates of the SFR (star-formation rate) or the total infrared luminosity from just one infrared band. However these relations are generally derived for entire galaxies, which are known to contain a large scale diffuse emission that is not necessarily related to the latest star-formation episode. We provide new relations to estimate the SFR from resolved star-forming regions at 100 mum and 160 mum. We select individual star-forming regions in the nearby (840 kpc) galaxy M33. We estimate the SFR combining the emission in Halpha and at 24 mum to calibrate the emission at 100 mum and 160 mum as SFR estimators, as mapped with PACS/Herschel. The data are obtained in the framework of the HERM33ES open time key project. There is less emission in the HII regions at 160 mum than at 100 mum. Over a dynamic range of almost 2 dex in Sigma(SFR) we find that the 100 mum emission is a nearly linear estimator of the SFR, whereas that at 160 mum is slightly superlinear. The behaviour of individual star-forming regions is surprisingly similar to that of entire galaxies. At high Sigma(SFR), star formation drives the dust temperature, whereas uncertainties and variations in radiation-transfer and dust-heated processes dominate at low Sigma(SFR). Detailed modelling of both galaxies and individual star forming regions will be needed to interpret similarities and differences between the two and assess the fraction of diffuse emission in galaxies.

preprint2010arXiv

A Study of the Gas-Star Formation Relation over Cosmic Time

We use the first systematic data sets of CO molecular line emission in z~1-3 normal star forming galaxies for a comparison of the dependence of galaxy-averaged star formation rates on molecular gas masses at low and high redshifts, and in different galactic environments. Although the current high-z samples are still small and biased toward the luminous and massive tail of the actively star-forming 'main-sequence', a fairly clear picture is emerging. Independent of whether galaxy integrated quantities or surface densities are considered, low- and high-z SFG galaxy populations appear to follow similar molecular gas-star formation relations with slopes 1.1 to 1.2. The gas-depletion time scale in these SFGs grows from 0.5 Gyrs at z~2 to 1.5 Gyrs at z~0. Because star formation depletion times are significantly smaller than the Hubble time at all redshifts sampled, star formation rates and gas fractions are set by the balance between gas accretion from the halo and stellar feedback. In contrast, very luminous gas rich major mergers at both low-z and high-z produce on average 4 to10 times more far-infrared luminosity per unit gas mass. Only some fraction of this difference can be explained by uncertainties in gas-mass or luminosity estimators; much of it must be intrinsic. The most likely interpretation is that the star formation relation is driven by global dynamical effects. For a given mass, the more compact merger systems produce stars more rapidly because their gas clouds are more compressed with shorter dynamical times, so that they churn more quickly through the available gas reservoir than the typical normal disk galaxies. When the dependence on galactic dynamical time scale is explicitly included, disk galaxies and mergers appear to follow similar gas to star-formation relations. The mergers may be forming stars at slightly higher efficiencies than the disks.

preprint2010arXiv

ALMA and the First Galaxies

ALMA will become fully operational in a few years and open a new window on primordial galaxies. The mm and submm domain is privileged, since the peak of dust emission between 60 and 100 microns is redshifted there for z= 5-10, and the continuum benefits from a negative K-correction. At least 100 times more sources than with present instruments could be discovered, so that more normal galaxies, with lower luminosities than huge starbursts and quasars will be surveyed. The high spatial resolution will suppress the confusion, which plagues today single dish bolometer surveys. Several CO lines detected in broad-band receivers will determine the redshift of objects too obscured to be seen in the optical. With the present instrumentation, only the most massive and gas rich objects have been detected in CO at high z, most of them being ultra-luminous starbursts with an extremely high star formation efficiency. However, selection biases are omni-present in this domain, and ALMA will statistically clarify the evolution of star formation efficiency, being fully complementary to JWST and ELTs.

preprint2010arXiv

Cool gas and dust in M33: Results from the Herschel M33 extended survey (HERM33ES)

We present an analysis of the first space-based far-IR-submm observations of M 33, which measure the emission from the cool dust and resolve the giant molecular cloud complexes. With roughly half-solar abundances, M33 is a first step towards young low-metallicity galaxies where the submm may be able to provide an alternative to CO mapping to measure their H$_2$ content. In this Letter, we measure the dust emission cross-section $σ$ using SPIRE and recent CO and \HI\ observations; a variation in $σ$ is present from a near-solar neighborhood cross-section to about half-solar with the maximum being south of the nucleus. Calculating the total H column density from the measured dust temperature and cross-section, and then subtracting the \HI\ column, yields a morphology similar to that observed in CO. The H$_2$/\HI\ mass ratio decreases from about unity to well below 10% and is about 15% averaged over the optical disk. The single most important observation to reduce the potentially large systematic errors is to complete the CO mapping of M 33.

preprint2010arXiv

Deep Herschel view of obscured star formation in the Bullet cluster

We use deep, five band (100-500um) data from the Herschel Lensing Survey (HLS) to fully constrain the obscured star formation rate, SFR_FIR, of galaxies in the Bullet cluster (z=0.296), and a smaller background system (z=0.35) in the same field. Herschel detects 23 Bullet cluster members with a total SFR_FIR = 144 +/- 14 M_sun yr^-1. On average, the background system contains brighter far-infrared (FIR) galaxies, with ~50% higher SFR_FIR (21 galaxies; 207 +/- 9 M_sun yr^-1). SFRs extrapolated from 24um flux via recent templates (SFR_24) agree well with SFR_FIR for ~60% of the cluster galaxies. In the remaining ~40%, SFR_24 underestimates SFR_FIR due to a significant excess in observed S_100/S_24 (rest frame S_75/S_18) compared to templates of the same FIR luminosity.

preprint2010arXiv

Discovery of an extremely bright submillimeter galaxy at z=3.93

Serendipitously we have discovered a rare, bright submillimeter galaxy (SMG) with a flux density of 30 +/- 2 mJy at lambda=1.2mm, using MAMBO2 at the IRAM 30-meter millimeter telescope. Although no optical counterpart is known for MM18423+5938, we were able to measure the redshift z=3.92960 +/- 0.00013 from the detection of CO lines using the IRAM Eight MIxer Receiver (EMIR). In addition, by collecting all available photometric data in the far-infrared and radio to constrain its spectral energy distribution, we derive the FIR luminosity 4.8 10^14/m Lsol and mass 6.0 10^9/m Msol for its dust, allowing for a magnification factor m caused by a probable gravitational lens. The corresponding star-formation rate is 8.3 10^4/m Msol/yr. The detection of three lines of the CO rotational ladder, and a significant upper limit for a fourth CO line, allow us to estimate an H2 mass of between 1.9 10^11/m Msol and 1.1 10^12/m Msol. The two lines CI(3p1-3p0) and CI(3p2-3p1) were clearly detected and yield a [CI]/[H2] number abundance between 1.4 10^-5 and 8.0 10^-5. Upper limits are presented for emission lines of HCN, HCO^+, HNC, H_2O and other molecules observed. The moderate excitation of the CO lines is indicative of an extended starburst, and excludes the dominance of an AGN in heating this high-redshift SMG.

preprint2010arXiv

First detection of the Sunyaev Zel'dovich effect increment at lambda < 650 um

The Sunyaev-Zel'dovich (SZ) effect is a spectral distortion of the cosmic microwave background as observed through the hot plasma in galaxy clusters. This distortion is a decrement in the CMB intensity for lambda > 1.3 mm, an increment at shorter wavelengths, and small again by lambda ~250 um. As part of the Herschel Lensing Survey (HLS) we have mapped 1E0657-56 (the Bullet cluster) with SPIRE with bands centered at 250, 350 and 500 um and have detected the SZ effect at the two longest wavelengths. The measured SZ effect increment central intensities are Delta I_{0} = 0.097 +- 0.019 MJy sr^{-1} at 350 um and Delta I_{0} = 0.268 +- 0.031 MJy sr^{-1} at 500 um, consistent with the SZ effect spectrum derived from previous measurements at 2 mm. No other diffuse emission is detected. The presence of the finite temperature SZ effect correction is preferred by the SPIRE data at a significance of 2.1 sigma, opening the possibility that the relativistic SZ effect correction can be constrained by SPIRE in a sample of clusters. The results presented here have important ramifications for both sub-mm measurements of galaxy clusters and blank field surveys with SPIRE.

preprint2010arXiv

Herschel observations of FIR emission lines in brightest cluster galaxies

The question of how much gas cools in the cores of clusters of galaxies has been the focus of many, multiwavelength studies in the past 30 years. In this letter we present the first detections of the strongest atomic cooling lines, [C II], [O I] and [N I] in two strong cooling flow clusters, A1068 and A2597, using Herschel PACS. These spectra indicate that the substantial mass of cold molecular gas (>10^9 Mo) known to be present in these systems is being irradiated by intense UV radiation, most probably from young stars. The line widths of these FIR lines indicate that they share dynamics similar but not identical to other ionised and molecular gas traced by optical, near-infrared and CO lines. The relative brightness of the FIR lines compared to CO and FIR luminosity is consistent with other star-forming galaxies indicating that the properties of the molecular gas clouds in cluster cores and the stars they form are not unusual. These results provide additional evidence for a reservoir of cold gas that is fed by the cooling of gas in the cores of the most compact clusters and provide important diagnostics of the temperature and density of the dense clouds this gas resides in.

preprint2010arXiv

Herschel photometry of brightest cluster galaxies in cooling flow clusters

The dust destruction timescales in the cores of clusters of galaxies are relatively short given their high central gas densities. However, substantial mid-infrared and sub-mm emission has been detected in many brightest cluster galaxies. In this letter we present Herschel PACS and SPIRE photometry of the brightest cluster galaxy in three strong cooling flow clusters, A1068, A2597 and Zw3146. This photometry indicates that a substantial mass of cold dust is present (>3 x 10^7 Mo) at temperatures significantly lower (20-28K) than previously thought based on limited MIR and/or sub-mm results. The mass and temperature of the dust appear to match those of the cold gas traced by CO with a gas-to-dust ratio of 80-120.

preprint2010arXiv

High molecular gas fractions in normal massive star forming galaxies in the young Universe

Stars form from cold molecular interstellar gas. Since this is relatively rare in the local Universe, galaxies like the Milky Way form only a few new stars per year. Typical massive galaxies in the distant Universe formed stars an order of magnitude more rapidly. Unless star formation was significantly more efficient, this difference suggests that young galaxies were much more gas rich. Molecular gas observations in the distant Universe have so far been largely restricted to very luminous, rare objects, including mergers and quasars. Here we report the results of a systematic survey of molecular gas in samples of typical massive star forming galaxies at <z>~1.2 and 2.3, when the Universe was 40% and 24% of its current age. Our measurements provide empirical evidence that distant star forming galaxies indeed were gas rich, and that the star formation efficiency is not strongly dependent on cosmic epoch. The average fraction of cold gas relative to total galaxy baryonic mass at z= 2.3 and z=1.2 is ~44% and 34%, three to ten times higher than in today's massive spiral galaxies. The slow decrease between z~2 and 1 probably requires a mechanism of semi-continuous replenishment of fresh gas to the young galaxies.

preprint2010arXiv

High-resolution mapping of the physical conditions in two nearby active galaxies based on 12CO(1-0), (2-1) and (3-2) lines

We present a detailed analysis of high resolution observations of the three lowest CO transitions in two nearby active galaxies, NGC4569 and NGC4826. The CO(1-0) and (2-1) lines were observed with the Plateau de Bure Interferometer and the CO(3-2) line with the Submillimeter Array. Combining these data allows us to compare the emission in the three lines and to map the line ratios, R21=I_{CO(2-1)}/I_{CO(1-0)} and R32=I_{CO(3-2)}/I_{CO(1-0)} at a resolution of ~2", i.e., a linear resolution of 160 pc for NGC4569 and 40 pc for NGC4826. In both galaxies the emission in the three lines is similarly distributed spatially and in velocity, and CO is less excited (R32<0.6) than in the Galactic Center or the centers of other active galaxies studied so far. According to a pseudo-LTE model the molecular gas in NGC4569 is cold and mainly optically thick in the CO(1-0) and (2-1) lines; less than 50% of the gas is optically thin in the CO(3-2) line. LVG modeling suggests the presence of an elongated ring of cold and dense gas coinciding with the ILR of the stellar bar. More excited gas is resolved in the circumnuclear disk of NGC4826. According to our pseudo-LTE model this corresponds to warmer gas with a ~50% of the CO(3-2) emission being optically thin. LVG modeling indicates the presence of a semicircular arc of dense and cold gas centered on the dynamical center and ~70 pc in radius. The gas temperature increases and its density decreases toward the center. A near side/far side asymmetry noticeable in the CO, R32 and Pa-alpha maps suggests that opacity effects play a role. Examining published CO maps of nearby active galaxies we find similar asymmetries suggesting that this could be a common phenomenon in active galaxies. These mainly qualitative results open new perspectives for the study of active galaxies with the future Atacama Large Millimeter/submillimeter Array.

preprint2010arXiv

Improving the identification of high-z Herschel sources with position priors and optical/NIR and FIR/mm photometric redshifts

We present preliminary results about the detection of high redshift (U)LIRGs in the Bullet cluster field by the PACS and SPIRE instruments within the Herschel Lensing Survey (HLS) Program. We describe in detail a photometric procedure designed to recover robust fluxes and deblend faint Herschel sources near the confusion noise. The method is based on the use of the positions of Spitzer/MIPS 24 um sources as priors. Our catalogs are able to reliably (5 sigma) recover galaxies with fluxes above 6 and 10 mJy in the PACS 100 and 160 um channels, respectively, and 12 to 18 mJy in the SPIRE bands. We also obtain spectral energy distributions covering the optical through the far-infrared/millimeter spectral ranges of all the Herschel detected sources, and analyze them to obtain independent estimations of the photometric redshift based on either stellar population or dust emission models. We exemplify the potential of the combined use of Spitzer position priors plus independent optical and IR photometric redshifts to robustly assign optical/NIR counterparts to the sources detected by Herschel and other (sub-)mm instruments.

preprint2010arXiv

Molecular and Atomic Gas in the Local Group Galaxy M33

We present high resolution large scale observations of the molecular and atomic gas in the Local Group Galaxy M33. The observations were carried out using the HERA at the 30m IRAM telescope in the CO(2-1) line achieving a resolution of 12"x2.6 km/s, enabling individual GMCs to be resolved. The observed region mainly along the major axis out to a radius of 8.5 kpc, and covers the strip observed with HIFI/PACS Spectrometers as part of the HERM33ES Herschel key program. The achieved sensitivity in main beam temperature is 20-50 mK at 2.6 km/s velocity resolution. The CO(2-1) luminosity of the observed region is 1.7\pm0.1x10^7 Kkm/s pc^2, corresponding to H2 masses of 1.9x10^8 Msun (including He), calculated with a NH2/ICO twice the Galactic value due to the half-solar metallicity of M33. HI 21 cm VLA archive observations were reduced and the mosaic was imaged and cleaned using the multi-scale task in CASA, yielding a series of datacubes with resolutions ranging from 5" to 25". The HI mass within a radius of 8.5 kpc is estimated to be 1.4x10^9 Msun. The azimuthally averaged CO surface brightness decreases exponentially with a scale length of 1.9\pm0.1 kpc whereas the atomic gas surface density is constant at Sigma_HI=6\pm2 Msun/pc^2 deprojected to face-on. The central kiloparsec H_2 surface density is Sigma_H2=8.5\pm0.2 Msun/pc^2. The star formation rate per unit molecular gas (SF Efficiency, the rate of transformation of molecular gas into stars), as traced by the ratio of CO to Halpha and FIR brightness, is constant with radius. The SFE appears 2-4 times greater than of large spiral galaxies. A morphological comparison of molecular and atomic gas with tracers of star formation shows good agreement between these maps both in terms of peaks and holes. A few exceptions are noted. Several spectra, including those of a molecular cloud situated more than 8 kpc from the galaxy center, are presented.

preprint2010arXiv

Molecular Gas in SAURON Early-Type Galaxies: Detection of 13CO and HCN Emission

In a pilot project to study the relationship between star formation and molecular gas properties in nearby normal early-type galaxies, we used the IRAM 30m telescope to observe the 13CO(J=1-0), 13CO(J=2-1), HCN(J=1-0) and HCO+(J=1-0) line emission in the four galaxies of the SAURON sample with the strongest 12CO emission. We report the detection of 13CO emission in all four SAURON sources and HCN emission in three sources, while no HCO+ emission was found to our detection limits in any of the four galaxies. We find that the 13CO/12CO ratios of three SAURON galaxies are somewhat higher than those in galaxies of different Hubble types. The HCN/12CO and HCN/13CO ratios of all four SAURON galaxies resemble those of nearby Seyfert and dwarf galaxies with normal star formation rates, rather than those of starburst galaxies. The HCN/HCO+ ratio is found to be relatively high (i.e., >1) in the three SAURON galaxies with detected HCN emission, mimicking the behaviour in other star-forming galaxies but being higher than in starburst galaxies. When compared to most galaxies, it thus appears that 13CO is enhanced (relative to 12CO) in three out of four SAURON galaxies and HCO+ is weak (relative to HCN) in three out of three galaxies. All three galaxies detected in HCN follow the standard HCN-infrared luminosity and dense gas fraction-star formation efficiency correlations. As already suggested by 12CO observations, when traced by infrared radiation, star formation in the three SAURON galaxies thus appears to follow the same physical laws as in galaxies of different Hubble types. The star formation rate and fraction of dense molecular gas however do not reach the high values found in nearby starburst galaxies, but rather resemble those of nearby normal star-forming galaxies.

preprint2010arXiv

PACS and SPIRE photometer maps of M33: First results of the Herschel M33 extended survey (HERM33ES)

Within the framework of the HERM33ES key project, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 micron wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 and 500 micron using also MIPS/Spitzer data, to derive first estimates of the dust physical conditions. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M33. An underlying diffuse disk is seen in all SPIRE maps (250-500 micron). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged flux densities. The two component model, with beta fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is found to be about 60K plus/minus 10K. The temperature of the cold component drops significantly from about 24K in the inner 2 kpc radius to 13K beyond 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for beta=1.5, averaged over the galaxy, is higher than the solar value by a factor of 1.5 and is roughly in agreement with the subsolar metallicity of M33.

preprint2010arXiv

Properties of compact 250 μm emission and HII regions in M33 (HERM33ES)

Within the framework of the HERM33ES Key Project, using the high resolution and sensitivity of the Herschel photometric data, we study the compact emission in the Local Group spiral galaxy M33 to investigate the nature of the compact SPIRE emission sources. We extracted a catalogue of sources at 250um in order to investigate the nature of this compact emission. Taking advantage of the unprecedented Herschel resolution at these wavelengths, we also focus on a more precise study of some striking Halpha shells in the northern part of the galaxy. We present a catalogue of 159 compact emission sources in M33 identified by SExtractor in the 250um SPIRE band that is the one that provides the best spatial resolution. We also measured fluxes at 24um and Halpha for those 159 extracted sources. The morphological study of the shells also benefits from a multiwavelength approach including Halpha, far-UV from GALEX, and infrared from both Spitzer IRAC 8um and MIPS 24um in order to make comparisons. For the 159 compact sources selected at 250um, we find a very strong Pearson correlation coefficient with the MIPS 24um emission (r24 = 0.94) and a rather strong correlation with the Halpha emission, although with more scatter (rHa = 0.83). The morphological study of the Halpha shells shows a displacement between far-ultraviolet, Halpha, and the SPIRE bands. The cool dust emission from SPIRE clearly delineates the Halpha shell structures. The very strong link between the 250um compact emission and the 24um and Halpha emissions, by recovering the star formation rate from standard recipes for HII regions, allows us to provide star formation rate calibrations based on the 250um compact emission alone. The different locations of the Halpha and far-ultraviolet emissions with respect to the SPIRE cool dust emission leads to a dynamical age of a few Myr for the Halpha shells and the associated cool dust.

preprint2010arXiv

Reionization by UV or X-ray sources

We present simulations of the 21-cm signal during the epoch of reionization. We focus on properly modeling the absorption regime in the presence of inhomogeneous Wouthuysen-Field effect and X-ray heating. We ran radiative transfer simulations for three bands in the source spectrum (Lyman, UV, and X-ray) to fully account for these processes. We find that the brightness temperature fluctuation of the 21 cm signal has an amplitude greater than 100 mK during the early reionization, up to 10 times greater than the typical amplitude of a few 10 mK obtained during the later emission phase. More importantly, we find that even a rather high contribution from QSO-like sources only damps the absorption regime without erasing it. Heating the IGM with X-ray takes time. Our results show that observations of the early reionization will probably benefit from a higher signal-to-noise value than during later stages. After analyzing the statistical properties of the signal (power spectrum and PDF) we find three diagnostics to constrain the level of X-ray, hence the nature of the first sources.

preprint2010arXiv

The dilution peak, metallicity evolution, and dating of galaxy interactions and mergers

Strong inflows of gas from the outer disk to the inner kiloparsecs are induced during the interaction of disk galaxies. This inflow of relatively low-metallicity gas dilutes the metallicity of the circumnuclear gas. We have investigated several aspects of the process as the timing and duration of the dilution and its correlation with the induced star formation. We analysed major (1:1) gas-rich interactions and mergers, spanning a range of initial orbital characteristics. Star formation and metal enrichment from SNe are included in our model. Our results show that the strongest trend is between the star formation rate and the dilution of the metals in the nuclear region; i.e., the more intense the central burst of star formation, the more the gas is diluted. This trend comes from strong inflows of relatively metal-poor gas from the outer regions of both disks, which fuels the intense star formation and lowers the overall metallicity for a time. The strong inflows happen on timescales of about 10^8 years or less, and the most intense star formation and lowest gas phase metallicities are seen generally after the first pericentre passage. As the star formation proceeds and the merger advances, the dilution reduces and enrichment becomes dominant - ultimately increasing the metallicity of the circumnuclear gas to a level higher than the initial metallicities of the merging galaxies. The "fly-bys" - pairs that interact but do not merge - also cause some dilution. We even see some dilution early in the merger or in the "fly-bys" and thus do not observe a strong trend between the nuclear metallicities and separation in our simulations until the merger is well advanced. We also analyse the O and Fe enrichment of the ISM, and show that the evolution of the alpha/Fe ratios, as well as the dilution of the central gas metallicity, can be used as a clock for "dating" the interaction.

preprint2010arXiv

The Energetics of Molecular Gas in NGC 891 from H2 and FIR Spectroscopy

We have studied the molecular hydrogen energetics of the edge-on spiral galaxy NGC\,891, using a 34-position map in the lowest three pure rotational H$_2$ lines observed with the Spitzer Infrared Spectrograph. The S(0), S(1), and S(2) lines are bright with an extinction corrected total luminosity of $\sim2.8 \times 10^{7}$ L$_{\odot}$, or 0.09\% of the total-infrared luminosity of NGC\,891. The H$_2$ line ratios are nearly constant along the plane of the galaxy -- we do not observe the previously reported strong drop-off in the S(1)/S(0) line intensity ratio in the outer regions of the galaxy, so we find no evidence for the very massive cold CO-free molecular clouds invoked to explain the past observations. The H$_2$ level excitation temperatures increase monotonically indicating more than one component to the emitting gas. More than 99\% of the mass is in the lowest excitation (T$_{ex}$ $\sim$125 K) ``warm'' component. In the inner galaxy, the warm H$_2$ emitting gas is $\sim$15\% of the CO(1-0)-traced cool molecular gas, while in the outer regions the fraction is twice as high. This large mass of warm gas is heated by a combination of the far-UV photons from stars in photo-dissociation regions (PDRs) and the dissipation of turbulent kinetic energy. Including the observed far-infrared [OI] and [CII] fine-structure line emission and far-infrared continuum emission in a self-consistent manner to constrain the PDR models, we find essentially all of the S(0) and most (70\%) of the S(1) line arises from low excitation PDRs, while most (80\%) of the S(2) and the remainder of the S(1) line emission arises from low velocity microturbulent dissipation.

preprint2010arXiv

The far-infrared/submillimeter properties of galaxies located behind the Bullet cluster

The Herschel Lensing Survey (HLS) takes advantage of gravitational lensing by massive galaxy clusters to sample a population of high-redshift galaxies which are too faint to be detected above the confusion limit of current far-infrared/submillimeter telescopes. Measurements from 100-500 micron bracket the peaks of the far-infrared spectral energy distributions of these galaxies, characterizing their infrared luminosities and star formation rates. We introduce initial results from our science demonstration phase observations, directed toward the Bullet cluster (1E0657-56). By combining our observations with LABOCA 870 micron and AzTEC 1.1 mm data we fully constrain the spectral energy distributions of 19 MIPS 24 micron selected galaxies which are located behind the cluster. We find that their colors are best fit using templates based on local galaxies with systematically lower infrared luminosities.This suggests that our sources are not like local ultra-luminous infrared galaxies in which vigorous star formation is contained in a compact highly dust-obscured region. Instead, they appear to be scaled up versions of lower luminosity local galaxies with star formation occurring on larger physical scales.

preprint2010arXiv

The Herschel Lensing Survey (HLS): overview

The Herschel Lensing Survey (HLS) will conduct deep PACS and SPIRE imaging of ~40 massive clusters of galaxies. The strong gravitational lensing power of these clusters will enable us to penetrate through the confusion noise, which sets the ultimate limit on our ability to probe the Universe with Herschel. Here, we present an overview of our survey and a summary of the major results from our Science Demonstration Phase (SDP) observations of the Bullet Cluster (z=0.297). The SDP data are rich, allowing us to study not only the background high-redshift galaxies (e.g., strongly lensed and distorted galaxies at z=2.8 and 3.2) but also the properties of cluster-member galaxies. Our preliminary analysis shows a great diversity of far-infrared/submillimeter spectral energy distributions (SEDs), indicating that we have much to learn with Herschel about the properties of galaxy SEDs. We have also detected the Sunyaev-Zel'dovich (SZ) effect increment with the SPIRE data. The success of this SDP program demonstrates the great potential of the Herschel Lensing Survey to produce exciting results in a variety of science areas.

preprint2009arXiv

Dark Matter tested with satellites

Recently, the distribution of velocity dispersion as far as 400kpc around red isolated galaxies was derived from statistical studies of satellites in the SDSS (Klypin & Prada 2009). This could help to constrain dark matter models at intermediate scales. We compare the predictions of different DM distributions, LCDM with NFW or cored profiles, and also modified gravity models, with observations. It is shown how the freedom in the various parameters (radial distribution of satellites, velocity anisotropy, external field effect), prevents to disentangle the models, which all can give pretty good fits to the data. In all cases, realistic radial variations of velocity anisotropy are used for the satellites, and a constant stellar-mass to light ratio for the host galaxies.

preprint2009arXiv

Generation of rotationally dominated galaxies by mergers of pressure-supported progenitors

Through the analysis of a set of numerical simulations of major mergers between initially non-rotating, pressure supported progenitor galaxies with a range of central mass concentrations, we have shown that: (1) it is possible to generate elliptical-like galaxies, with v/sigma > 1 outside one effective radius, as a result of the conversion of orbital- into internal-angular momentum; (2) the outer regions acquire part of the angular momentum first; (3) both the baryonic and the dark matter components of the remnant galaxy acquire part of the angular momentum, the relative fractions depend on the initial concentration of the merging galaxies. For this conversion to occur the initial baryonic component must be sufficiently dense and/or the encounter should take place on a orbit with high angular momentum. Systems with these hybrid properties have been recently observed through a combination of stellar absorption lines and planetary nebulae for kinematic studies of early-type galaxies. Our results are in qualitative agreement with such observations and demonstrate that even mergers composed of non-rotating, pressure-supported progenitor galaxies can produce early-type galaxies with significant rotation at large radii.

preprint2009arXiv

Molecular Gas in NUclei of GAlaxies (NUGA) XIII. The interacting Seyfert 2/LINER galaxy NGC 5953

We present $^{12}$CO(1--0) and $^{12}$CO(2--1) maps of the interacting Seyfert 2/LINER galaxy NGC 5953 obtained with the IRAM interferometer at resolutions of 2\farcs1 $\times$ 1\farcs4 and 1\farcs1 $\times$ 0\farcs7, respectively. The CO emission is distributed over a disk of diameter $\sim$16\arcsec ($\sim$2.2 kpc), within which are several, randomly distributed peaks. The strongest peak does not coincide with the nucleus, but is instead offset from the center, $\sim2-3^{\prime\prime}$ ($\sim$340 pc) toward the west/southwest. The kinematics of the molecular component are quite regular, as is typical of a rotating disk. We also compared the $^{12}$CO distribution of NGC 5953 with observations at other wavelengths in order to study correlations between different tracers of the interstellar medium. Using NIR images, we computed the gravity torques exerted by the stellar potential on the gas. The torques are predominantly positive in both $^{12}$CO(1--0) and $^{12}$CO(2--1), suggesting that gas is not flowing into the center, and less than 5% of the gas angular momentum is exchanged in each rotation. This comes from the regular and almost axisymmetric total mass and gas distributions in the center of the galaxy. In NGC 5953, the AGN is apparently not being actively fueled in the current epoch.

preprint2009arXiv

Molecular Gas in NUclei of GAlaxies (NUGA): XII. The head-on collision in NGC1961

We present high-resolution CO(1-0) and CO(2-1) maps of the LINER 2 galaxy NGC1961. This galaxy is unusual among late-type (Sc) disk galaxies in having a very large radial extent and inferred dynamical mass. We propose a head-on collision scenario to explain the perturbed morphology of this galaxy- both the off-centered rings and the inflated radius. This scenario is supported by the detection of a steep velocity gradient in the CO(1-0) map at the position of a southwest peak in radio continuum and near-infrared emission. This peak would represent the remnant of the disrupting companion. We use numerical models to demonstrate the plausibility of the scenario. While ram pressure stripping could in principle be important for shocking the atomic gas and produce the striking head-tail morphology, the non detection of this small galaxy group in X-ray emission suggests that any hot intragroup medium has too low a density. The collision with the companion of mass ratio 1:4 produces a lopsided wave, which accounts for the sharp boundary observed in the atomic gas on the southern side. We argue that NGC1961 has a strongly warped disk, which gives the false impression of a nearly face-on system; the main disk is actually more edge-on, and this error in the true inclination has led to the surprisingly high dynamical mass for a morphologically late-type galaxy. In addition, the outward propagating ring artificially enlarges the disk. The collision de-stabilizes the inner disk and can provide gas inflow to the active nucleus.

preprint2009arXiv

MOND and the Galaxies

We review galaxy formation and dynamics under the MOND hypothesis of modified gravity, and compare to similar galaxies in Newtonian dynamics with dark matter. The aim is to find peculiar predictions both to discriminate between various hypotheses, and to make the theory progress through different constraints, touching the interpolation function, or the fundamental acceleration scale. Galaxy instabilities, forming bars and bulges at longer term, evolve differently in the various theories, and help to bring constraints, together with the observations of bar frequency. Dynamical friction and the predicted merger rate could be a sensitive test of theories. The different scenarios of galaxy formation are compared within the various theories and observations.

preprint2009arXiv

On the survival of metallicity gradients to major dry-mergers

Using a large suite of galaxies with a variety of concentrations and metallicity gradients, we study the evolution of non-dissipative ("dry") equal mass mergers. Our purpose in generating these simulations is to explore how the metallicity gradient in dry mergers depends on the structure and metallicity gradients of the galaxies involved in the merger. Specifically, we would like to answer: Could dry mergers lead to metallicity gradients as observed in elliptical galaxies in the local Universe? Do dry mergers always lead to a flattening of the initial (i.e., pre-merger) gradient? From this modeling, we conclude that: The ratio of the remnant and the initial galaxy slopes span a wide range of values, up to values greater than 1 (with values greater than one resulting only when companions have gradients twice the progenitor). For a merger between two ellipticals having identical initial metallicity slopes (i.e., equal companion and galaxy slopes), the metallicity profile of the remnant flattens, with a final gradient about 0.6 times the initial one. Ellipticals can maintain their original pre-merger metallicity gradient if the companion slope is sufficiently steep. The amount of flattening neither depends on the characteristics of the orbit of the progenitors or on their initial concentration. Given the diversity in outcomes of the mergers, we conclude that dry mergers do not violate any observational constraints on the systematic characteristics of metallicity gradients in local ellipticals. In fact, dry mergers may be important within the context of the results of our simulations and may explain the large scatter in gradients for massive ellipticals and the relative paucity of massive ellipticals with no or shallow metallicity gradients.

preprint2009arXiv

Ram pressure stripping of tilted galaxies

Ram pressure stripping of galaxies in clusters can yield gas deficient disks. Previous numerical simulations based on various approaches suggested that, except for near edge-on disk orientations, the amount of stripping depends very little on the inclination angle. Following our previous study of face-on stripping, we extend the set of parameters with the disk tilt angle and explore in detail the effects of the ram pressure on the interstellar content (ISM) of tilted galaxies that orbit in various environments of clusters, with compact or extended distributions of the intra-cluster medium (ICM). We further study how results of numerical simulations could be estimated analytically. A grid of numerical simulations with varying parameters is produced using the tree/SPH code GADGET with a modified method for calculating the ISM-ICM interaction. These SPH calculations extend the set of existing results obtained from different codes using various numerical techniques. The simulations confirm the general trend of less stripping at orientations close to edge-on. The dependence on the disk tilt angle is more pronounced for compact ICM distributions, however it almost vanishes for strong ram pressure pulses. Although various hydrodynamical effects are present in the ISM-ICM interaction, the main quantitative stripping results appear to be roughly consistent with a simple scenario of momentum transfer from the encountered ICM. This behavior can also be found in previous simulations. To reproduce the numerical results we propose a fitting formula depending on the disk tilt angle and on the column density of the encountered ICM. Such a dependence is superior to that on the peak ram pressure used in previous simple estimates.

preprint2008arXiv

Simulations of galactic disks including a dark baryonic component

$Λ$CDM numerical simulations predict that the "missing baryons" reside in a Warm-Hot gas phase in the over-dense cosmic filaments. However, there are now several theoretical and observational arguments that support the fact that galactic disks may be more massive than usually thought, containing a substantial fraction of the "missing baryons". Hereafter, we present new N-body simulations of galactic disks, where the gas content has been multiplied by a factor 5. The stability of the disk is ensured by assuming that the ISM is composed out of two partially coupled phases, a warm phase, corresponding the observed CO and HI gas and a cold collisionless phase corresponding to the unseen baryons.

preprint2007arXiv

NUclei of GAlaxies: V. Radio emission in 7 NUGA sources

We present high angular resolution radio snap-shot observations of seven nearby low-luminosity active galaxies (LLAGN) from the NUclei of GAlaxies (NUGA) survey. The observations were conducted with MERLIN and EVN/VLBI at 18cm and 6cm. At all observed angular resolutions and frequencies, we find indications for extended emission in about ~40% of the sources, consistent with the decrease of flux with increasing angular resolution. The extended components resemble jet emission in a majority of cases, consistent with the optically thin synchrotron emission implied by their steep spectra. We consider the compact 6cm EVN/VLBI radio emission of our sources in the context of the ``fundamental plane'' that previous LLAGN studies identified within the three-dimensional parameter space of radio luminosity, X-ray luminosity, and black hole mass. We demonstrate, using NGC7217 and NGC1068 as particular examples, that high-resolution, multi-epoch radio observations offer useful information about the origin of offsets from the fundamental plane.

preprint2006arXiv

Secular evolution in galaxies

New observations in favour of a significant role of secular evolution are reviewed: central star formation boosted in pseudo-bulge barred galaxies, relations between bulge and disk, evidence for rejuvenated bulges. Numerical simulations have shown that secular evolution can occur through a cycle of bar formation and destruction, in which the gas plays a major role. Since bars are weakened or destroyed in gaseous disks, the high frequency of bars observed today requires external cold gas accretion, to replenish the disk and allow a new bar formation. The rate of gas accretion from external filaments is compatible with what is observed in cosmological simulations.

preprint2006arXiv

The Building of Galactic Disks: Insights from the Triangulum Spiral Galaxy Messier 33

The Triangulum Spiral Galaxy Messier 33 offers unique insights into the building of a galactic disk. We identify spectacular arcs of intermediate age (0.6 Gyr - 2 Gyr) stars in the low-metallicity outer disk. The northern arc spans approx. 120 degrees in azimuth and up to 5 arcmin in width. The arcs are located 2-3 disk scale lengths from the galaxy centre (where 1 disk scale length is equivalent to 0.1 degrees in the V-band) and lie precisely where there is a warp in the HI profile of M33. Warps and infall are inextricably linked (Binney, 1992). We present spectroscopy of candidate stars in the outer northern arc, secured using the Keck I telescope in Hawaii. The target stars have estimated visual magnitudes as faint as V ~ 25m. Absorption bands of CN are seen in all spectra reported in this review talk, confirming their carbon star status. Also presented are PAH emissivity radial profiles generated from IRAC observations of M33 using the Spitzer Space Telescope. A dramatic change of phase in the m=2 Fourier component is detected at the domain of the arcs. M33 serves as an excellent example how the disks of spiral galaxies in our Universe are built: as dynamically open systems, growing from the inward, outward.

preprint2005arXiv

Continuum emission in NGC1068 and NGC3147: Indications for a turnover in the core spectra

We present new interferometric observations of the continuum emission at mm wavelengths in the Seyfert galaxies NGC1068 and NGC3147. Three mm continuum peaks are detected in NGC1068, one centered on the core, one associated with the jet and the third one with the counter-jet. This is the first significant detection of the radio jet and counter-jet at mm wavelengths in NGC1068. While the fluxes of the jet components agree with a steep spectral index extrapolated from cm-wavelengths, the core fluxes indicate a turnover of the inverted cm- into a steep mm-spectrum at roughly 50GHz which is most likely caused by electron-scattered synchrotron emission. As in NGC 1068, the spectrum of the pointlike continuum source in NGC3147 also shows a turnover between cm and mm-wavelengths at 25GHz resulting from synchrotron self-absorption different to NGC1068. This strongly resembles the spectrum of SgrA*, the weakly active nucleus of our own galaxy, and M81*, a link between SgrA* and Seyfert galaxies in terms of activity sequence, which may display a similar turnover.

preprint2004arXiv

Efficiency of the dynamical mechanism

The most extreme starbursts occur in galaxy mergers, and it is now acknowledged that dynamical triggering has a primary importance in star formation. This triggering is due partly to the enhanced velocity dispersion provided by gravitational instabilities, such as density waves and bars, but mainly to the radial gas flows they drive, allowing large amounts of gas to condense towards nuclear regions in a small time scale. Numerical simulations with several gas phases, taking into account the feedback to regulate star formation, have explored the various processes, using recipes like the Schmidt law, moderated by the gas instability criterion. May be the most fundamental parameter in starbursts is the availability of gas: this sheds light on the amount of external gas accretion in galaxy evolution. The detailed mechanisms governing gas infall in the inner parts of galaxy disks are discussed.

preprint2002arXiv

Properties of Dark Matter Haloes

An overview is presented of the main properties of dark matter haloes, as we know them from observations, essentially from rotation curves around spiral and dwarf galaxies. Detailed rotation curves are now known for more than a thousand galaxies, revealing that they are not so flat in the outer parts, but rising for late-types, and falling for early-types. A well established result now is that most bright galaxies are not dominated by dark matter inside their optical disks. Only for dwarfs and LSB (Low Surface Brightness galaxies) dark matter plays a dominant role in the visible regions. The 3D-shape of haloes are investigated through several methods, that will be discussed: polar rings, flaring of HI planes, X-ray isophotes. It is not yet possible with rotation curves to know how far haloes extend, but tentatives have been made. It will be shown that the dark matter appears to be coupled to the gas in spirals and dwarfs, suggesting that dark baryons could play the major role in rotation curves. Theories proposing to replace the non-baryonic dark matter by a different dynamical or gravity law, such as MOND, have to take into account the dark baryons, especially since their spatial distribution is likely to be quite different from the visible matter.

preprint2001arXiv

MIFS: A Mega Integral Field Spectrograph for the VLT

We describe MIFS, a second generation integral-field spectrograph for the VLT, operating in the visible wavelength range. It combines a 1'x1' field of view with the improved spatial resolution provided by multi-conjugate adaptive optics and covers a large simultaneous spectral range (0.6-1.0 microns). A separate mode exploits the highest spatial resolution provided by adaptive optics. With this unique combination of capabilities, MIFS has a wide domain of application and a large discovery potential. The MIFS low-spatial resolution mode (sampled at 0.2") combined with the initial MCAO capabilities planned for the VLT will provide ultra-deep fields with a limiting magnitude for spectroscopy of R~28. MIFS will improve the present day detection limit of Lyman alpha emitters by a factor of 100, and will detect low-mass star-forming galaxies to z~7. The MIFS high-spatial resolution mode (3"x3" field sampled at 0.01") is optimized for the next step in (MC)AO. It will probe, e.g., the relationship between supermassive central black holes and their host galaxy and the physics of winds from accretion disks in young stellar objects at unprecedented spatial resolution. MIFS will extend Europe's lead in integral-field spectroscopy. It capitalizes on new developments in adaptive optics, and is a key step towards instrumentation for OWL.

preprint2000arXiv

Fueling the AGN

Active Galactic Nuclei are fueled from material (gas or stars) that are in general far away from the gravitational influence of the central black hole, the engine thought to be responsible for their activity. The required material has a lot of angular momentum, that is a priori quite difficult to evacuate. The various dynamical mechanisms that may play a role in this game are reviewed, including m=2 perturbations (bars and spirals), m=1 perturbations (spirals, warps, lopsidedness), and tidal interactions between galaxies and mergers. In the latest stages of the merger, a binary black hole could be formed, and its influence on the dynamics and fueling is discussed. Starbursts are often associated with AGN, and the nature of their particular connection, and their role in the nuclear fueling is described. Evolution of the fueling efficiency with redshift is addressed.

preprint1999arXiv

CMB and Molecules at High Redshift

It becomes possible now to detect cold molecules at high redshift in the millimeter domain. Since the first discovery in 1992 by Brown and van den Bout of CO lines at z=2.28 in a gravitationally lensed starburst galaxy, nearly ten objects are now known to possess large quantities of molecular gas beyond z=1 and up to z = 5, through millimeter and sub-millimeter emission lines. The continuum dust emission is the most easily detected: in the mm domain, the emission is stronger for the more redshifted objects. For the CO lines, the situation is less favorable, and the reported detections are helped by gravitational amplification. The increase of the CMB temperature T_{bg} with redshift helps the rotational line excitation (especially at high z), but not its detection. Absorption in front of quasars is a more sensitive probe of cold gas at high redshift, able to detect individual clouds of a few solar masses (instead of 10^{10} Mo for emission). From the diffuse components, one can measure the cosmic black body temperature as a function of redshift. The high column densities component allow to observe important molecules not observable from the ground, like O2, H2O and LiH for example.

preprint1999arXiv

Molecular gas and the dynamics of galaxies

In this review, I discuss some highlights of recent research on molecular gas in galaxies; large-scale CO maps of nearby galaxies are being made, which extend our knowledge on global properties, radial gradients, and spiral structure of the molecular ISM. Very high resolution are provided by the interferometers, that reveal high velocity gradients in galaxy nuclei, and formation of embedded structures, like bars within bars. Observation of the CO and other lines in starburst galaxies have questioned the H2-to-CO conversion factor. Surveys of dwarfs have shown how the conversion factor depends on metallicity. The molecular content is not deficient in galaxy clusters, as is the atomic gas. Galaxy interactions are very effective to enhance gas concentrations and trigger starbursts. Nuclear disks or rings are frequently observed, that concentrate the star formation activity. Since the density of starbursting galaxies is strongly increasing with redshift, the CO lines are a privileged tool to follow evolution of galaxies and observe the ISM dynamics at high redshift: due to the high excitation of the molecular gas, the stronger high-$J$ CO lines are redshifted into the observable band, which facilitates the detection.

preprint1999arXiv

N-body simulations of globular cluster tides

We present N-body simulations of globular clusters, in orbits around the Galaxy, in order to study quantitatively and geometrically the tidal effects they encounter. The clusters are modelised with multi-mass King-Michie models (Michie 1963), including mass segregation at initial conditions. The Galaxy is modelled as realistic as possible, with three components: bulge, disk and dark halo. The main finding is that there exist two giant tidal tails around the globuler cluster in permanence along its orbit, whatever this orbit. The length of these tails is of the order of 5 tidal radii, or greater. The escaped stars are distributed radially as a power law in density, with a slope of -4. The tails present substructures, or clumps, that are the relics of the strongest shocks. Due to the compressive disk-shocking, the clusters display a prolate shape which major axis is precessing around the z axis. The tails are preferentially formed by the lowest mass stars, as expected, so that the tidal truncation increases mass segregation. Internal rotation of the cluster increases the mass loss. The flattening of dark matter cannot influence significantly the dynamics of the clusters. The orientation and the strength of the tidal tails are signatures of the last disk crossing, so that observed tidal tails can constrain strongly the cluster orbit and the galactic model (vertical scale of the disc).

preprint1999arXiv

Small Group Dynamics and Extended Gas

Interactions between galaxies have spectacular effects on gas dynamics, and small galaxy groups are a privileged place to investigate them. In particular, they could test the existence of cold H2 gas as dark matter in the outer parts of galaxies. HI observations have revealed that galaxies in small groups are deficient in atomic gas, like in richer galaxy clusters such as Virgo, although in a lesser extent. Galaxy interactions could be the cause of this deficiency, stripping the gas out of galaxies and enriching the inter-cluster medium (ICM) in hot gas, which in turn strips gas through ram pressure. Alternatively, the gas present at the formation of the group could have been heated to its virial temperature, and be observed now as X-rays. The dynamical processes related to this extended gas in small galaxy groups are reviewed.

preprint1998arXiv

Dust and Molecules at High Redshift

In the last years, progress has been very rapid in the domain of molecules at high redshift, and we know in better detail now the molecular and dust content in several systems beyond z=1 and up to z = 5. The first discovery in 1992 by Brown and van den Bout of CO lines at z=2.28 in a gravitationally lensed starburst galaxy, strongly stimulated searches of other systems, but these were harder than foreseen, and less than 10 other systems have been discovered in CO emission. Redshifts range between 2 and 5, the largest being BR1202-0725 at z=4.69. Most of these systems, if not all, are gravitationally amplified objects. Some have been discovered first through their dust emission, relatively easy to detect because of the negative K-correction effect. The detection of all these systems could give an answer about the debated question of the star-formation rate as a function of redshift. The maximum of star-formation rate, found around z=2 from optical studies, could shift to higher z if the most remote objects are hidden by dust. Absorption in front of quasars can also probe cold gas at high redshift, taking advantage of very high spatial (milli arcsec) and spectral (30m/s) resolutions. From the diffuse components, one can measure the cosmic black body temperature as a function of redshift. All these preliminary studies will be carried out at large scales with future millimeter instruments, and some perspectives are given.

preprint1998arXiv

Fractal Dimensions and Scaling Laws in the Interstellar Medium and Galaxy Distributions: a new Field Theory Approach

We develop a field theoretical approach to the cold interstellar medium (ISM) and large structure of the universe. We show that a non-relativistic self- gravitating gas in thermal equilibrium with variable number of atoms or fragments is exactly equivalent to a field theory of a scalar field phi(x) with exponential self-interaction. We analyze this field theory perturbatively and non-perturbatively through the renormalization group(RG).We show scaling behaviour (critical) for a continuous range of the physical parameters as the temperature. We derive in this framework the scaling relation M(R) \sim R^{d_H} for the mass on a region of size R, and Delta v \sim R^\frac12(d_H -1) for the velocity dispersion. For the density-density correlations we find a power-law behaviour for large distances \sim |r_1 - r_2|^{2D - 6}.The fractal dimension D turns to be related with the critical exponent νby D = 1/ ν. Mean field theory yields ν= 1/2, D = 2. Both the Ising and the mean field values are compatible with the present ISM observational data:1.4\leq D \leq 2. We develop a field theoretical approach to the galaxy distribution considering a gas of self-gravitating masses on the FRW background, in quasi-thermal equi- librium. We show that it exhibits scaling behaviour by RG methods. The galaxy correlations are computed without assuming homogeneity. We find <ρ({\vec r_0})ρ({\vec r_0} + {\vec r}) > \sim r^{D-3} $. The theory allows to compute the three and higher density correlators without any assumption.We find that the connected N-points density scales as r_1^{N(D-3)}, when $ r_1 >> r_i, 2\leq i \leq N $. There are no free parameters in this theory.

preprint1998arXiv

Fractal Structures and Scaling Laws in the Universe: Statistical Mechanics of the Self-Gravitating Gas

Fractal structures are observed in the universe in two very different ways. Firstly, in the gas forming the cold interstellar medium in scales from 10^{-4} pc till 100 pc. Secondly, the galaxy distribution has been observed to be fractal in scales up to hundreds of Mpc. We give here a short review of the statistical mechanical (and field theoretical) approach developed by us. We consider a non-relativistic self-gravitating gas in thermal equilibrium at temperature T inside a volume V. The statistical mechanics of such system has special features and, as is known, the thermodynamical limit does not exist in its customary form. Moreover, the treatments through microcanonical, canonical and grand canonical ensembles yield different results.We present here for the first time the equation of state for the self-gravitating gas in the canonical ensemble. We find that it has the form p = [N T/ V] f(eta), where p is the pressure, N is the number of particles and η\equiv {G m^2 N \over V^{1/3} T}. The N \to\infty and V \to\infty limit exists keeping ηfixed. We compute the function f(η) using Monte Carlo simulations and for small eta analytically. We compute the thermodynamic quantities of the system as free energy, entropy, chemical potential, specific heat, compressibility and speed of sound. We reproduce the well-known gravitational phase transition associated to the Jeans' instability. Namely, a gaseous phase for eta < eta_c and a condensed phase for eta > eta_c. Moreover, we derive the precise behaviour of the physical quantities near the transition. In particular, the pressure vanishes as p \sim(eta_c-eta)^B with B \sim 0.2 and eta_c \sim 1.6 and the energy fluctuations diverge as \sim(eta_c-eta)^{B-1}. The speed of sound decreases monotonically and approaches the value sqrt{T/6} at the transition.

preprint1998arXiv

Molecular gas in the barred spiral M100 - II. 12 CO(1-0) interferometer observations and numerical simulations

Using the IRAM interferometer we have mapped at high resolution (2"x1") the 12CO(1-0) emission in the nucleus of the doubled barred SABbc spiral M100. Molecular gas is distributed in a two spiral arm structure starting from the end points of the nuclear bar (r=600 pc) and a central source. The kinematics of the gas indicates the existence of a steep rotation curve (180 km/s at 100 pc) and strong streaming motions characteristic of a trailing spiral wave inside corotation. Gas flow simulations analyse the gas response to a gravitational potential derived from the K-band plate, including the two nested bars. We develop two families of models: first, a single pattern speed solution shared by the outer bar+spiral and by the nuclear bar, and secondly, a two independent bars solution, where the nuclear bar is dynamically decoupled and rotates faster than the primary bar. We found the best fit solution consisting of a fast pattern (160 km/s/kpc) for the nuclear bar (with corotation at 1.2 kpc) decoupled from the slow pattern of the outer bar+spiral (23 km/s/kpc, with corotation at 8-9 kpc).

preprint1998arXiv

The fractal structure of the universe : a new field theory approach

While the universe becomes more and more homogeneous at large scales, statistical analysis of galaxy catalogs have revealed a fractal structure at small-scales (λ< 100 h^{-1} Mpc), with a fractal dimension D=1.5-2 (Sylos Labini et al 1996). We study the thermodynamics of a self-gravitating system with the theory of critical phenomena and finite-size scaling and show that gravity provides a dynamical mechanism to produce this fractal structure. We develop a field theoretical approach to compute the galaxy distribution, assuming them to be in quasi-isothermal equilibrium. Only a limited, (although large), range of scales is involved, between a short-distance cut-off below which other physics intervene, and a large-distance cut-off, where the thermo- dynamic equilibrium is not satisfied. The galaxy ensemble can be considered at critical conditions, with large density fluctuations developping at any scale. From the theory of critical phenomena, we derive the two independent critical exponents nu and eta and predict the fractal dimension D = 1/nu to be either 1.585 or 2, depending on whether the long-range behaviour is governed by the Ising or the mean field fixed points, respectively. Both set of values are compatible with present observations. In addition, we predict the scaling behaviour of the gravitational potential to be r^{-(1 + eta)/2}. That is, r^{-0.5} for mean field or r^{- 0.519} for the Ising fixed point. The theory allows to compute the three and higher density correlators without any assumption or Ansatz. We find that the N-points density scales as r_1^{(N-1)(D-3)}, when r_1 >> r_i, 2 leq i leq N . There are no free parameters in this theory.

preprint1998arXiv

The Role of Interactions and Mergers

Galaxy evolution depends strongly on the environment. Dynamical interactions and subsequent evolution make galaxies more concentrated, with higher surface densities, and also trigger star-formation, which consumes the available cold gas. Already at z=0, a large variety of galaxy types are observed, with different evolution stages, from the unevolved gas rich dwarf irregulars, or low surface brightness galaxies, to early-type concentrated galaxies, with no remaining gas. The dynamical processes of galaxy interactions, including internal evolution, are reviewed and evidence is shown of much larger interaction/merging rate at high redshift.

preprint1997arXiv

Detection of water at z = 0.685 towards B0218+357

We report the detection of the H_2O molecule in absorption at a redshift z = 0.68466 in front of the gravitationally lensed quasar B0218+357. We detect the fundamental transition of ortho-water at 556.93 GHz (redshifted to 330.59 GHz). The line is highly optically thick and relatively wide (15 km/s FWHM), with a profile that is similar to that of the previously detected CO(2--1) and HCO^+(2--1) optically thick absorption lines toward this quasar. From the measured level of the continuum at 330.59 GHz, which corresponds to the level expected from the power-law spectrum $S(ν) \propto ν^{-0.25}$ already observed at lower frequencies, we deduce that the filling factor of the H_2O absorption is large. It was already known from the high optical thickness of the CO, ^{13}CO and C^{18}O lines that the molecular clouds entirely cover one of the two lensed images of the quasar (all its continuum is absorbed); our present results indicate that the H_2O clouds are covering a comparable surface. The H_2O molecules are therefore not confined to small cores with a tiny filling factor, but are extended over parsec scales. The H_2O line has a very large optical depth, and only isotopic lines could give us the water abundance. We have also searched for the 183 GHz line in absorption, obtaining only an upper limit; this yields constraints on the excitation temperature.

preprint1996arXiv

Search for molecular absorption in the tori of active galactic nuclei

We describe a search for molecular absorption at millimetre wavelengths associated with dusty molecular tori in active galactic nuclei (AGN). The sample observed consists of 11 flat-spectrum radio sources known to have red optical to infra-red colours plus two steep-spectrum narrow-line radio galaxies. Spectra of the sources were obtained in the 3-, 2- and 1.3-millimetre bands at frequencies corresponding to common molecular transitions of CO, HCO+, HCN and CS at the AGN redshift. No absorptions were detected in any of the sources. We calculated upper limits to the column density in molecular absorption, using an excitation temperature of 10 K, to be N(CO) < 10^{15} - 10^{16} cm^-2, equivalent to hydrogen columns of order N(H) < 10^{19} - 10^{20} cm^-2. These limits are significantly lower than the values N(H) \approx (2 - 6) 10^{21} cm^-2 that might be expected if the red colours of these sources were due to dust absorption at the quasar redshift as suggested by Webster et al. (1995). Should the excitation temperature of the molecular transitions be higher than 100K, the upper limits to the H2 column densities would be greater than those derived from the red colours. To explain the lack of molecular absorption we conclude that either the optical extinction takes place outside the host galaxy (along the line of sight), or the excitation temperature of the molecular transitions is very high, or the obscuration is not associated with significant amounts of cold molecular gas. It is quite possible that the hard X-ray flux from the central source of these AGN is strong enough to photo-dissociate the molecules.

preprint1996arXiv

Secular Evolution of Galaxy Morphologies

Today we have numerous evidences that spirals evolve dynamically through various secular or episodic processes, such as bar formation and destruction, bulge growth and mergers, sometimes over much shorter periods than the standard galaxy age of 10-15 Gyr. This, coupled to the known properties of the Hubble sequence, leads to a unique sense of evolution: from Sm to Sa. Linking this to the known mass components provides new indications on the nature of dark matter in galaxies. The existence of large amounts of yet undetected dark gas appears as the most natural option. Bounds on the amount of dark stars can be given since their formation is mostly irreversible and requires obviously a same amount of gas.

preprint1996arXiv

Self-gravity as an explanation of the fractal structure of the interstellar medium

The gas clouds of the interstellar medium have a fractal structure, the origin of which has generally been thought to lie in turbulence. The energy of turbulence could come from galactic rotation at large-scale, then cascade down to be dissipated on small-scales by viscosity; it has been suggested that such turbulence helps to prevent massive molecular clouds from collapsing in response to their own gravity. Here we show that, on the contrary, self-gravity itself may be the dominant factor in making clouds fractal. We develop a field-theory approach to the structure of clouds, assuming them to be isothermal, and with only gravitational interactions; we find that the observed fractal dimension of the clouds arise naturally from this approach. Although this result does not imply that turbulence is not important, it does demonstrate that the fractal structure can be understood without it.

preprint1996arXiv

Very Cold Gas and Dark Matter

We have recently proposed a new candidate for baryonic dark matter: very cold molecular gas, in near-isothermal equilibrium with the cosmic background radiation at 2.73 K. The cold gas, of quasi-primordial abundances, is condensed in a fractal structure, resembling the hierarchical structure of the detected interstellar medium. We present some perspectives of detecting this very cold gas, either directly or indirectly. The H$_2$ molecule has an "ultrafine" structure, due to the interaction between the rotation-induced magnetic moment and the nuclear spins. But the lines fall in the km domain, and are very weak. The best opportunity might be the UV absorption of H$_2$ in front of quasars. The unexpected cold dust component, revealed by the COBE/FIRAS submillimetric results, could also be due to this very cold H$_2$ gas, through collision-induced radiation, or solid H$_2$ grains or snowflakes. The $γ$-ray distribution, much more radially extended than the supernovae at the origin of cosmic rays acceleration, also points towards and extended gas distribution.

preprint1993arXiv

Is Dark Matter in Spiral Galaxies Cold Gas? I. Observational Constraints and Dynamical Clues About Galaxy Evolution

Based on dynamical constraints about the Hubble sequence evolution, observational data and a number of "conspiracies", we propose that the dark matter around spiral galaxies is in the form of cold gas, essentially in molecular form and rotationally supported. (full A&A paper (in press) available by anonymous ftp at obssd8.unige.ch in /pub/fractal as postscript file: dm_paper_I.ps (170k), or papers I & II + figures as a compressed tar file dm_papers.Z.tar (2.1 Mb)).

preprint1993arXiv

Is Dark Matter in Spiral Galaxies Cold Gas? II. Fractal Models and Star Non-Formation

Gas cloud models taking into account the recently disclosed fractal structure of cold gas are set up, showing that large errors in the classical gas mass determination based on smooth cloud models can easily follow if the gas is in reality fractal. Fractal clouds must present both optically thin and optically thick clumps in any single wavelength observations. The observed fractal dimension of the cold ISM suggests that mass underestimates by a factor 10 or more are typical. Due to its low temperature (around 3 K), and its condensed fractal structure, together with its low metallicity, the outer gas would be almost invisible for usual detectors. (A&A paper (in press) and figures available by anonymous ftp at obssd8.unige.ch in /pub/fractal as postscript file: dm_paper_II.ps (251k) and FIG*.ps, or papers I & II + figures as a compressed tar file dm_papers.Z.tar (2.1 Mb)).