Source author record

J. Graciá-Carpio

J. Graciá-Carpio 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

30works
2topics
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

30 published item(s)

preprint2023arXiv

Euclid preparation: XXIII. Derivation of galaxy physical properties with deep machine learning using mock fluxes and H-band images

Next generation telescopes, like Euclid, Rubin/LSST, and Roman, will open new windows on the Universe, allowing us to infer physical properties for tens of millions of galaxies. Machine learning methods are increasingly becoming the most efficient tools to handle this enormous amount of data, because they are often faster and more accurate than traditional methods. We investigate how well redshifts, stellar masses, and star-formation rates (SFR) can be measured with deep learning algorithms for observed galaxies within data mimicking the Euclid and Rubin/LSST surveys. We find that Deep Learning Neural Networks and Convolutional Neutral Networks (CNN), which are dependent on the parameter space of the training sample, perform well in measuring the properties of these galaxies and have a better accuracy than methods based on spectral energy distribution fitting. CNNs allow the processing of multi-band magnitudes together with $H_{\scriptscriptstyle\rm E}$-band images. We find that the estimates of stellar masses improve with the use of an image, but those of redshift and SFR do not. Our best results are deriving i) the redshift within a normalised error of less than 0.15 for 99.9$\%$ of the galaxies with S/N>3 in the $H_{\scriptscriptstyle\rm E}$-band; ii) the stellar mass within a factor of two ($\sim0.3 \rm dex$) for 99.5$\%$ of the considered galaxies; iii) the SFR within a factor of two ($\sim0.3 \rm dex$) for $\sim$70$\%$ of the sample. We discuss the implications of our work for application to surveys as well as how measurements of these galaxy parameters can be improved with deep learning.

preprint2023arXiv

Euclid preparation. XXX. Performance assessment of the NISP Red-Grism through spectroscopic simulations for the Wide and Deep surveys

This work focuses on the pilot run of a simulation campaign aimed at investigating the spectroscopic capabilities of the Euclid Near-Infrared Spectrometer and Photometer (NISP), in terms of continuum and emission line detection in the context of galaxy evolutionary studies. To this purpose we constructed, emulated, and analysed the spectra of 4992 star-forming galaxies at $0.3 \leq z \leq 2.5$ using the NISP pixel-level simulator. We built the spectral library starting from public multi-wavelength galaxy catalogues, with value-added information on spectral energy distribution (SED) fitting results, and from Bruzual and Charlot (2003) stellar population templates. Rest-frame optical and near-IR nebular emission lines were included using empirical and theoretical relations. We inferred the 3.5$σ$ NISP red grism spectroscopic detection limit of the continuum measured in the $H$ band for star-forming galaxies with a median disk half-light radius of \ang{;;0.4} at magnitude $H= 19.5\pm0.2\,$AB$\,$mag for the Euclid Wide Survey and at $H = 20.8\pm0.6\,$AB$\,$mag for the Euclid Deep Survey. We found a very good agreement with the red grism emission line detection limit requirement for the Wide and Deep surveys. We characterised the effect of the galaxy shape on the detection capability of the red grism and highlighted the degradation of the quality of the extracted spectra as the disk size increases. In particular, we found that the extracted emission line signal to noise ratio (SNR) drops by $\sim\,$45$\%$ when the disk size ranges from \ang{;;0.25} to \ang{;;1}. These trends lead to a correlation between the emission line SNR and the stellar mass of the galaxy and we demonstrate the effect in a stacking analysis unveiling emission lines otherwise too faint to detect.

preprint2022arXiv

Euclid Preparation: XIV. The Complete Calibration of the Color-Redshift Relation (C3R2) Survey: Data Release 3

The Complete Calibration of the Color-Redshift Relation (C3R2) survey is obtaining spectroscopic redshifts in order to map the relation between galaxy color and redshift to a depth of i ~ 24.5 (AB). The primary goal is to enable sufficiently accurate photometric redshifts for Stage IV dark energy projects, particularly Euclid and the Roman Space Telescope, which are designed to constrain cosmological parameters through weak lensing. We present 676 new high-confidence spectroscopic redshifts obtained by the C3R2 survey in the 2017B-2019B semesters using the DEIMOS, LRIS, and MOSFIRE multi-object spectrographs on the Keck telescopes. Combined with the 4454 redshifts previously published by this project, the C3R2 survey has now obtained and published 5130 high-quality galaxy spectra and redshifts. If we restrict consideration to only the 0.2 < z(phot) < 2.6 range of interest for the Euclid cosmological goals, then with the current data release C3R2 has increased the spectroscopic redshift coverage of the Euclid color space from 51% (as reported by Masters et al. 2015) to the current 91%. Once completed and combined with extensive data collected by other spectroscopic surveys, C3R2 should provide the spectroscopic calibration set needed to enable photometric redshifts to meet the cosmology requirements for Euclid, and make significant headway toward solving the problem for Roman.

preprint2022arXiv

Euclid preparation: XIX. Impact of magnification on photometric galaxy clustering

We investigate the importance of lensing magnification for estimates of galaxy clustering and its cross-correlation with shear for the photometric sample of Euclid. Using updated specifications, we study the impact of lensing magnification on the constraints and the shift in the estimation of the best fitting cosmological parameters that we expect if this effect is neglected. We follow the prescriptions of the official Euclid Fisher matrix forecast for the photometric galaxy clustering analysis and the combination of photometric clustering and cosmic shear. The slope of the luminosity function (local count slope), which regulates the amplitude of the lensing magnification, and the galaxy bias have been estimated from the Euclid Flagship simulation.We find that magnification significantly affects both the best-fit estimation of cosmological parameters and the constraints in the galaxy clustering analysis of the photometric sample. In particular, including magnification in the analysis reduces the 1$σ$ errors on $Ω_{\text{m},0}, w_{0}, w_a$ at the level of 20-35%, depending on how well we will be able to independently measure the local count slope. In addition, we find that neglecting magnification in the clustering analysis leads to shifts of up to 1.6$σ$ in the best-fit parameters. In the joint analysis of galaxy clustering, cosmic shear, and galaxy-galaxy lensing, magnification does not improve precision, but it leads to an up to 6$σ$ bias if neglected. Therefore, for all models considered in this work, magnification has to be included in the analysis of galaxy clustering and its cross-correlation with the shear signal ($3\times2$pt analysis) for an accurate parameter estimation.

preprint2022arXiv

Euclid preparation: XX. The Complete Calibration of the Color-Redshift Relation survey: LBT observations and data release

The Complete Calibration of the Color-Redshift Relation survey (C3R2) is a spectroscopic programme designed to empirically calibrate the galaxy color-redshift relation to the Euclid depth (I_E=24.5), a key ingredient for the success of Stage IV dark energy projects based on weak lensing cosmology. A spectroscopic calibration sample as representative as possible of the galaxies in the Euclid weak lensing sample is being collected, selecting galaxies from a self-organizing map (SOM) representation of the galaxy color space. Here, we present the results of a near-infrared H- and K-bands spectroscopic campaign carried out using the LUCI instruments at the LBT. For a total of 251 galaxies, we present new highly-reliable redshifts in the 1.3<= z <=1.7 and 2<= z<=2.7 ranges. The newly-determined redshifts populate 49 SOM cells which previously contained no spectroscopic measurements and almost double the occupation numbers of an additional 153 SOM cells. A final optical ground-based observational effort is needed to calibrate the missing cells in particular in the redshift range 1.7<= z<=2.7 that lack spectroscopic calibration. In the end, Euclid itself will deliver telluric-free NIR spectra that can complete the calibration.

preprint2022arXiv

Euclid preparation. XVIII. The NISP photometric system

Euclid will be the first space mission to survey most of the extragalactic sky in the 0.95-2.02 $μ$m range, to a 5$σ$ point-source median depth of 24.4 AB mag. This unique photometric data set will find wide use beyond Euclid's core science. In this paper, we present accurate computations of the Euclid Y_E, J_E and H_E passbands used by the Near-Infrared Spectrometer and Photometer (NISP), and the associated photometric system. We pay particular attention to passband variations in the field of view, accounting among others for spatially variable filter transmission, and variations of the angle of incidence on the filter substrate using optical ray tracing. The response curves' cut-on and cut-off wavelengths - and their variation in the field of view - are determined with 0.8 nm accuracy, essential for the photometric redshift accuracy required by Euclid. After computing the photometric zeropoints in the AB mag system, we present linear transformations from and to common ground-based near-infrared photometric systems, for normal stars, red and brown dwarfs, and galaxies separately. A Python tool to compute accurate magnitudes for arbitrary passbands and spectral energy distributions is provided. We discuss various factors from space weathering to material outgassing that may slowly alter Euclid's spectral response. At the absolute flux scale, the Euclid in-flight calibration program connects the NISP photometric system to Hubble Space Telescope spectrophotometric white dwarf standards; at the relative flux scale, the chromatic evolution of the response is tracked at the milli-mag level. In this way, we establish an accurate photometric system that is fully controlled throughout Euclid's lifetime.

preprint2021arXiv

$Euclid$ preparation: XV. Forecasting cosmological constraints for the $Euclid$ and CMB joint analysis

The combination and cross-correlation of the upcoming $Euclid$ data with cosmic microwave background (CMB) measurements is a source of great expectation since it will provide the largest lever arm of epochs, ranging from recombination to structure formation across the entire past light cone. In this work, we present forecasts for the joint analysis of $Euclid$ and CMB data on the cosmological parameters of the standard cosmological model and some of its extensions. This work expands and complements the recently published forecasts based on $Euclid$-specific probes, namely galaxy clustering, weak lensing, and their cross-correlation. With some assumptions on the specifications of current and future CMB experiments, the predicted constraints are obtained from both a standard Fisher formalism and a posterior-fitting approach based on actual CMB data. Compared to a $Euclid$-only analysis, the addition of CMB data leads to a substantial impact on constraints for all cosmological parameters of the standard $Λ$-cold-dark-matter model, with improvements reaching up to a factor of ten. For the parameters of extended models, which include a redshift-dependent dark energy equation of state, non-zero curvature, and a phenomenological modification of gravity, improvements can be of the order of two to three, reaching higher than ten in some cases. The results highlight the crucial importance for cosmological constraints of the combination and cross-correlation of $Euclid$ probes with CMB data.

preprint2021arXiv

Euclid preparation: I. The Euclid Wide Survey

Euclid is an ESA mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey (EWS) in visible and near-infrared, covering roughly 15,000 square degrees of extragalactic sky on six years. The wide-field telescope and instruments are optimized for pristine PSF and reduced straylight, producing very crisp images. This paper presents the building of the Euclid reference survey: the sequence of pointings of EWS, Deep fields, Auxiliary fields for calibrations, and spacecraft movements followed by Euclid as it operates in a step-and-stare mode from its orbit around the Lagrange point L2. Each EWS pointing has four dithered frames; we simulate the dither pattern at pixel level to analyse the effective coverage. We use up-to-date models for the sky background to define the Euclid region-of-interest (RoI). The building of the reference survey is highly constrained from calibration cadences, spacecraft constraints and background levels; synergies with ground-based coverage are also considered. Via purposely-built software optimized to prioritize best sky areas, produce a compact coverage, and ensure thermal stability, we generate a schedule for the Auxiliary and Deep fields observations and schedule the RoI with EWS transit observations. The resulting reference survey RSD_2021A fulfills all constraints and is a good proxy for the final solution. Its wide survey covers 14,500 square degrees. The limiting AB magnitudes ($5σ$ point-like source) achieved in its footprint are estimated to be 26.2 (visible) and 24.5 (near-infrared); for spectroscopy, the H$_α$ line flux limit is $2\times 10^{-16}$ erg cm$^{-2}$ s$^{-1}$ at 1600 nm; and for diffuse emission the surface brightness limits are 29.8 (visible) and 28.4 (near-infrared) mag arcsec$^{-2}$.

preprint2021arXiv

Euclid preparation: XVI. Exploring the ultra low-surface brightness Universe with Euclid/VIS

While Euclid is an ESA mission specifically designed to investigate the nature of Dark Energy and Dark Matter, the planned unprecedented combination of survey area ($\sim15\,000$ deg$^2$), spatial resolution, low sky-background, and depth also make Euclid an excellent space observatory for the study of the low surface brightness Universe. Scientific exploitation of the extended low surface brightness structures requires dedicated calibration procedures yet to be tested. We investigate the capabilities of Euclid to detect extended low surface brightness structure by identifying and quantifying sky background sources and stray-light contamination. We test the feasibility of generating sky flat-fields to reduce large-scale residual gradients in order to reveal the extended emission of galaxies observed in the Euclid Survey. We simulate a realistic set of Euclid/VIS observations, taking into account both instrumental and astronomical sources of contamination, including cosmic rays, stray-light, zodiacal light, ISM, and the CIB, while simulating the effects of the presence of background sources in the FOV. We demonstrate that a combination of calibration lamps, sky flats and self-calibration would enable recovery of emission at a limiting surface brightness magnitude of $μ=29.5^{+0.08}_{-0.27} $ mag arcsec$^{-2}$ ($3σ$, $10\times10$ arcsec$^2$) in the Wide Survey, reaching regions 2 magnitudes deeper in the Deep Surveys. Euclid/VIS has the potential to be an excellent low surface brightness observatory. Covering the gap between pixel-to-pixel calibration lamp flats and self-calibration observations for large scales, the application of sky flat-fielding will enhance the sensitivity of the VIS detector at scales of larger than 1 degree, up to the size of the FOV, enabling Euclid to detect extended surface brightness structures below $μ=31$ mag arcsec$^{-2}$ and beyond.

preprint2020arXiv

Euclid preparation: VIII. The Complete Calibration of the Colour-Redshift Relation survey: VLT/KMOS observations and data release

The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

preprint2019arXiv

AGN feedback in a galaxy merger: Multi-phase, galaxy-scale outflows including a fast molecular gas blob ~6 kpc away from IRAS F08572+3915

To understand the role that AGN feedback plays in galaxy evolution we need in-depth studies of the multi-phase structure and energetics of galaxy-wide outflows. In this work we present new, deep ($\sim$50 hr) NOEMA CO(1-0) line observations of the molecular gas in the powerful outflow driven by the AGN in the ultra-luminous infrared galaxy IRAS F08572+3915. We spatially resolve the outflow, finding that its most likely configuration is a wide-angle bicone aligned with the kinematic major axis of the rotation disk. The molecular gas in the wind reaches velocities up to approximately $\pm$1200 km s$^{-1}$ and transports nearly 20% of the molecular gas mass in the system. We detect a second outflow component located $\sim$6 kpc north-west from the galaxy moving away at $\sim$900 km s$^{-1}$, which could be the result of a previous episode of AGN activity. The total mass and energetics of the outflow, which includes contributions from the ionized, neutral, warm and cold molecular gas phases is strongly dominated by the cold molecular gas. In fact, the molecular mass outflow rate is higher than the star formation rate, even if we only consider the gas in the outflow that is fast enough to escape the galaxy, which accounts for about $\sim$40% of the total mass of the outflow. This results in an outflow depletion time for the molecular gas in the central $\sim$1.5 kpc region of only $\sim3$ Myr, a factor of $\sim2$ shorter than the depletion time by star formation activity.

preprint2019arXiv

Molecular Gas Inflows and Outflows in Ultraluminous Infrared Galaxies at $z\sim0.2$ and one QSO at $z=6.1$

Aims. We aim to search and characterize inflows and outflows of molecular gas in four ultraluminous infrared galaxies (ULIRGs) at $z\sim0.2-0.3$ and one distant QSO at $z=6.13$. Methods. We use Herschel PACS and ALMA Band 7 observations of the hydroxyl molecule (OH) line at rest-frame wavelength 119 $μ$m which in absorption can provide unambiguous evidence for inflows or outflows of molecular gas in nuclear regions of galaxies. Our study contributes to double the number of OH observations of luminous systems at $z\sim0.2-0.3$, and push the search for molecular outflows based on the OH transition to $z\sim6$. Results. We detect OH high-velocity absorption wings in three of the four ULIRGs. In two cases, IRAS F20036-1547 and IRAS F13352+6402, the blueshifted absorption profiles indicate the presence of powerful and fast molecular gas outflows. Consistent with an inside-out quenching scenario, these outflows are depleting the central reservoir of molecular gas at a similar rate than the intense star formation activity. In the case of the starburst-dominated system IRAS 10091+4704, we detect an inverted P-Cygni profile that is unique among ULIRGs and indicates the presence of a fast ($\sim400$ km s$^{-1}$) inflow of molecular gas at a rate of $\sim100~M_{\odot}~{\rm yr}^{-1}$ towards the central region. Finally, we tentatively detect ($\sim3σ$) the OH doublet in absorption in the $z=6.13$ QSO ULAS J131911+095051. The OH feature is blueshifted with a median velocity that suggests the presence of a molecular outflow, although characterized by a modest molecular mass loss rate of $\sim200~M_{\odot}~{\rm yr}^{-1}$. This value is comparable to the small mass outflow rates found in the stacking of the [CII] spectra of other $z\sim6$ QSOs and suggests that ejective feedback in this phase of the evolution of ULAS J131911+095051 has subsided.

preprint2016arXiv

Broad [CII] line wings as tracer of molecular and multi-phase outflows in infrared bright galaxies

We report a tentative correlation between the outflow characteristics derived from OH absorption at $119\,μ\text{m}$ and [CII] emission at $158\,μ\text{m}$ in a sample of 22 local and bright ultraluminous infrared galaxies (ULIRGs). For this sample we investigate whether [CII] broad wings are a good tracer of molecular outflows, and how the two tracers are connected. Fourteen objects in our sample have a broad wing component as traced by [CII], and all of these also show OH119 absorption indicative of an outflow (in 1 case an inflow). The other eight cases, where no broad [CII] component was found, are predominantly objects with no OH outflow or a low-velocity ($\leq 100\,\text{km s}^{-1}$) OH outflow. The full width at half maximum (FWHM) of the broad [CII] component shows a trend with the OH119 blue-shifted velocity, although with significant scatter. Moreover, and despite large uncertainties, the outflow masses derived from OH and broad [CII] show a 1:1 relation. The main conclusion is therefore that broad [CII] wings can be used to trace molecular outflows. This may be particularly relevant at high redshift, where the usual tracers of molecular gas (like low-J CO lines) become hard to observe. Additionally, observations of blue-shifted Na I D $λλ5890,5896$ absorption are available for ten of our sources. Outflow velocities of Na I D show a trend with OH velocity and broad [CII] FWHM. These observations suggest that the atomic and molecular gas phases of the outflow are connected.

preprint2016arXiv

Constraints on the Broad Line Region Properties and Extinction in Local Seyferts

We use high spectral resolution (R > 8000) data covering 3800-13000Å to study the physical conditions of the broad line region (BLR) of nine nearby Seyfert 1 galaxies. Up to six broad HI lines are present in each spectrum. A comparison - for the first time using simultaneous optical to near-infrared observations - to photoionisation calculations with our devised simple scheme yields the extinction to the BLR at the same time as determining the density and photon flux, and hence distance from the nucleus, of the emitting gas. This points to a typical density for the HI emitting gas of 10$^{11}$cm$^{-3}$ and shows that a significant amount of this gas lies at regions near the dust sublimation radius, consistent with theoretical predictions. We also confirm that in many objects the line ratios are far from case B, the best-fit intrinsic broad-line H$α$/H$β$ ratios being in the range 2.5-6.6 as derived with our photoionization modeling scheme. The extinction to the BLR, based on independent estimates from HI and HeII lines, is A$_V$ $\le$ 3 for Seyfert 1-1.5s, while Seyfert 1.8-1.9s have A$_V$ in the range 4-8. A comparison of the extinction towards the BLR and narrow line region (NLR) indicates that the structure obscuring the BLR exists on scales smaller than the NLR. This could be the dusty torus, but dusty nuclear spirals or filaments could also be responsible. The ratios between the X-ray absorbing column N$_H$ and the extinction to the BLR are consistent with the Galactic gas-to-dust ratio if N$_H$ variations are considered.

preprint2016arXiv

Searching for molecular outflows in Hyper-Luminous Infrared Galaxies

We present constraints on the molecular outflows in a sample of five Hyper-Luminous Infrared Galaxies using Herschel observations of the OH doublet at 119 μm. We have detected the OH doublet in three cases: one purely in emission and two purely in absorption. The observed emission profile has a significant blueshifted wing suggesting the possibility of tracing an outflow. Out of the two absorption profiles, one seems to be consistent with the systemic velocity while the other clearly indicates the presence of a molecular outflow whose maximum velocity is about ~1500 km/s. Our analysis shows that this system is in general agreement with previous results on Ultra-luminous Infrared Galaxies and QSOs, whose outflow velocities do not seem to correlate with stellar masses or starburst luminosities (star formation rates). Instead the galaxy outflow likely arises from an embedded AGN.

preprint2016arXiv

The far-infrared emitting region in local galaxies and QSOs: Size and scaling relations

We use Herschel 70 to 160um images to study the size of the far-infrared emitting region in about 400 local galaxies and quasar (QSO) hosts. The sample includes normal `main-sequence' star-forming galaxies, as well as infrared luminous galaxies and Palomar-Green QSOs, with different levels and structures of star formation. Assuming Gaussian spatial distribution of the far-infrared (FIR) emission, the excellent stability of the Herschel point spread function (PSF) enables us to measure sizes well below the PSF width, by subtracting widths in quadrature. We derive scalings of FIR size and surface brightness of local galaxies with FIR luminosity, with distance from the star-forming main-sequence, and with FIR color. Luminosities L_FIR ~ 10^11 L_Sun can be reached with a variety of structures spanning 2 dex in size. Ultraluminous L_FIR >~ 10^12 L_Sun galaxies far above the main-sequence inevitably have small Re_70um ~ 0.5 kpc FIR emitting regions with large surface brightness, and can be close to optically thick in the FIR on average over these regions. Compared to these local relations, first ALMA sizes for the dust emission regions in high redshift galaxies, measured at somewhat longer rest wavelengths, suggest larger sizes at the same IR luminosity. We report a remarkably tight relation with 0.15 dex scatter between FIR surface brightness and the ratio of [CII] 158um emission and FIR emission -- the so-called [CII]-deficit is more tightly linked to surface brightness than to FIR luminosity or FIR color. Among 33 z <~ 0.1 PG QSOs with typical L_FIR/L_Bol,AGN ~ 0.1, 19 have a measured 70um half light radius, with median Re_70um = 1.1kpc. This is consistent with the FIR size for galaxies with similar L_FIR but lacking a QSO, in accordance with a scenario where the rest FIR emission of these types of QSOs is, in most cases, due to host star formation.

preprint2016arXiv

Thick Disks, and an Outflow, of Dense Gas in the Nuclei of Nearby Seyfert Galaxies

We discuss the dense molecular gas in central regions of nearby Seyfert galaxies, and report new arcsec resolution observations of HCN(1-0) and HCO$^+$(1-0) for 3 objects. In NGC 3079 the lines show complex profiles as a result of self-absorption and saturated continuum absorption. H$^{13}$CN reveals the continuum absorption profile, with a peak close to the galaxy's systemic velocity that traces disk rotation, and a second feature with a blue wing extending to $-350$km s$^{-1}$ that most likely traces a nuclear outflow. The morphological and spectral properties of the emission lines allow us to constrain the dense gas dynamics. We combine our kinematic analysis for these 3 objects, as well as another with archival data, with a previous comparable analysis of 4 other objects, to create a sample of 8 Seyferts. In 7 of these, the emission line kinematics imply thick disk structures on radial scales of $\sim$100pc, suggesting such structures are a common occurrence. We find a relation between the circumnuclear LHCN and Mdyn that can be explained by a gas fraction of 10% and a conversion factor αHCN $\sim$ 10 between gas mass and HCN luminosity. Finally, adopting a different perspective to probe the physical properties of the gas around AGN, we report on an analysis of molecular line ratios which indicates that the clouds in this region are not self-gravitating.

preprint2015arXiv

A deep Herschel/PACS observation of CO(40-39) in NGC 1068: a search for the molecular torus

Emission from high-J CO lines in galaxies has long been proposed as a tracer of X-ray dominated regions (XDRs) produced by AGN. Of particular interest is the question of whether the obscuring torus, which is required by AGN unification models, can be observed via high-J CO cooling lines. Here we report on the analysis of a deep Herschel-PACS observation of an extremely high J CO transition (40-39) in the Seyfert 2 galaxy NGC 1068. The line was not detected, with a derived 3$σ$ upper limit of $2 \times 10^{-17}\,\text{W}\,\text{m}^{-2}$. We apply an XDR model in order to investigate whether the upper limit constrains the properties of a molecular torus in NGC 1068. The XDR model predicts the CO Spectral Line Energy Distributions for various gas densities and illuminating X-ray fluxes. In our model, the CO(40-39) upper limit is matched by gas with densities $\sim 10^{6}-10^{7}\,\text{cm}^{-3}$, located at $1.6-5\,\text{pc}$ from the AGN, with column densities of at least $10^{25}\,\text{cm}^{-2}$. At such high column densities, however, dust absorbs most of the CO(40-39) line emission at $λ= 65.69\, μ$m. Therefore, even if NGC 1068 has a molecular torus which radiates in the CO(40-39) line, the dust can attenuate the line emission to below the PACS detection limit. The upper limit is thus consistent with the existence of a molecular torus in NGC 1068. In general, we expect that the CO(40-39) is observable in only a few AGN nuclei (if at all), because of the required high gas column density, and absorption by dust.

preprint2015arXiv

Insights on the Dusty Torus and Neutral Torus from Optical and X-ray Obscuration in a Complete Volume Limited Hard X-ray AGN Sample

We describe a complete volume limited sample of nearby active galaxies selected by their 14-195keV luminosity, and outline its rationale for studying the mechanisms regulating gas inflow and outflow. We describe also a complementary sample of inactive galaxies, selected to match the AGN host galaxy properties. The active sample appears to have no bias in terms of AGN type, the only difference being the neutral absorbing column which is two orders of magnitude greater for the Seyfert 2s. In the luminosity range spanned by the sample, log L_{14-195keV} [erg/s] = 42.4-43.7, the optically obscured and X-ray absorbed fractions are 50-65%. The similarity of these fractions to more distant spectroscopic AGN samples, although over a limited luminosity range, suggests that the torus does not strongly evolve with redshift. Our sample confirms that X-ray unabsorbed Seyfert 2s are rare, comprising not more than a few percent of the Seyfert 2 population. At higher luminosities, the optically obscured fraction decreases (as expected for the increasing dust sublimation radius), but the X-ray absorbed fraction changes little. We argue that the cold X-ray absorption in these Seyfert 1s can be accounted for by neutral gas in clouds that also contribute to the broad line region (BLR) emission; and suggest that a geometrically thick neutral gas torus co-exists with the BLR and bridges the gap to the dusty torus.

preprint2014arXiv

High-lying OH absorption, [C II] deficits, and extreme $L_{\mathrm{FIR}}/M_{\mathrm{H2}}$ ratios in galaxies

Herschel/PACS observations of 29 local (Ultra-)Luminous Infrared Galaxies, including both starburst and AGN-dominated sources as diagnosed in the mid-infrared/optical, show that the equivalent width of the absorbing OH 65 um Pi_{3/2} J=9/2-7/2 line (W_{eq}(OH65)) with lower level energy E_{low}~300 K, is anticorrelated with the [C ii]158 um line to far-infrared luminosity ratio, and correlated with the far-infrared luminosity per unit gas mass and with the 60-to-100 um far-infrared color. While all sources are in the active L_{IR}/M_{H2}>50 Lsun/Msun mode as derived from previous CO line studies, the OH65 absorption shows a bimodal distribution with a discontinuity at L_{FIR}/M_{H2}~100 Lsun/Msun. In the most buried sources, OH65 probes material partially responsible for the silicate 9.7 um absorption. Combined with observations of the OH 71 um Pi_{1/2} J=7/2-5/2 doublet (E_{low}~415 K), radiative transfer models characterized by the equivalent dust temperature, Tdust, and the continuum optical depth at 100 um, tau_{100}, indicate that strong [C ii]158 um deficits are associated with far-IR thick (tau_{100}>~0.7, N_{H}>~10^{24} cm^{-2}), warm (Tdust>~60 K) structures where the OH 65 um absorption is produced, most likely in circumnuclear disks/tori/cocoons. With their high L_{FIR}/M_{H2} ratios and columns, the presence of these structures is expected to give rise to strong [C ii] deficits. W_{eq}(OH65) probes the fraction of infrared luminosity arising from these compact/warm environments, which is >~30-50% in sources with high W_{eq}({OH65}). Sources with high W_{eq}({OH65}) have surface densities of both L_{IR} and M_{H2} higher than inferred from the half-light (CO or UV/optical) radius, tracing coherent structures that represent the most buried/active stage of (circum)nuclear starburst-AGN co-evolution.

preprint2013arXiv

Massive Molecular Outflows and Evidence for AGN Feedback from CO Observations

We study the properties of massive, galactic-scale outflows of molecular gas and investigate their impact on galaxy evolution. We present new IRAM PdBI CO(1-0) observations of local ULIRGs and QSO hosts: clear signature of massive and energetic molecular outflows, extending on kpc scales, is found in the CO(1-0) kinematics of four out of seven sources, with measured outflow rates of several 100 Msun/yr. We combine these new observations with data from the literature, and explore the nature and origin of massive molecular outflows within an extended sample of 19 local galaxies. We find that starburst-dominated galaxies have an outflow rate comparable to their SFR, or even higher by a factor of ~ 2-4, implying that starbursts can indeed be effective in removing cold gas from galaxies. Nevertheless, our results suggest that the presence of an AGN can boost the outflow rate by a large factor, which is found to increase with the L_AGN/L_bol ratio. The gas depletion time-scales due to molecular outflows are anti-correlated with the presence and luminosity of an AGN in these galaxies, and range from a few hundred million years in starburst galaxies, down to just a few million years in galaxies hosting powerful AGNs. In quasar hosts the depletion time-scales due to the outflow are much shorter than the depletion time-scales due to star formation. We estimate the outflow kinetic power and find that, for galaxies hosting powerful AGNs, it corresponds to about 5% of the AGN luminosity, as expected by models of AGN feedback. Moreover, we find that momentum rates of about 20 L_AGN/c are common among the AGN-dominated sources in our sample. For "pure" starburst galaxies our data tentatively support models in which outflows are mostly momentum-driven by the radiation pressure from young stars onto dusty clouds.

preprint2013arXiv

The Mrk 231 molecular outflow as seen in OH

We report on the Herschel/PACS observations of OH in Mrk 231, with detections in 9 doublets observed within the PACS range, and present radiative transfer models for the outflowing OH. Signatures of outflowing gas are found in up to 6 OH doublets with different excitation requirements. At least two outflowing components are identified, one with OH radiatively excited, and the other with low excitation, presumably spatially extended. Particularly prominent, the blue wing of the absorption detected in the in-ladder 2Pi_{3/2} J=9/2-7/2 OH doublet at 65 um, with E_lower=290 K, indicates that the excited outflowing gas is generated in a compact and warm (circum)nuclear region. Because the excited, outflowing OH gas in Mrk 231 is associated with the warm, far-IR continuum source, it is likely more compact (diameter of 200-300 pc) than that probed by CO and HCN. Nevertheless, its mass-outflow rate per unit of solid angle as inferred from OH is similar to that previously derived from CO, >~70x(2.5x10^{-6}/X_{OH}) Msun yr^{-1} sr^{-1}, where X_{OH} is the OH abundance relative to H nuclei. In spherical symmetry, this would correspond to >~850x(2.5x10^{-6}/X_{OH}) Msun yr^{-1}, though significant collimation is inferred from the line profiles. The momentum flux of the excited component attains ~15 L_{AGN}/c, with an OH column density of (1.5-3)x10^{17} cm^-2 and a mechanical luminosity of ~10^{11} Lsun. The detection of very excited OH peaking at central velocities indicates the presence of a nuclear reservoir of gas rich in OH, plausibly the 130-pc scale circumnuclear torus previously detected in OH megamaser emission, that may be feeding the outflow. An exceptional ^{18}OH enhancement, with OH/^{18}OH<~30 at both central and blueshifted velocities, is likely the result of interstellar-medium processing by recent starburst/SNe activity.

preprint2012arXiv

Dust temperature and CO-to-H2 conversion factor variations in the SFR-M* plane

Deep Herschel imaging and 12CO(2-1) line luminosities from the IRAM PdBI are combined for a sample of 17 galaxies at z>1 from the GOODS-N field. The sample includes galaxies both on and above the main sequence (MS) traced by star-forming galaxies in the SFR-M* plane. The far-infrared data are used to derive dust masses, Mdust. Combined with an empirical prescription for the dependence of the gas-to-dust ratio on metallicity (GDR), the CO luminosities and Mdust values are used to derive for each galaxy the CO-to-H2 conversion factor, alpha_co. Like in the local Universe, the value of alpha_co is a factor of ~5 smaller in starbursts compared to normal star-forming galaxies (SFGs). We also uncover a relation between alpha_co and dust temperature (Tdust; alpha_co decreasing with increasing Tdust) as obtained from modified blackbody fits to the far-infrared data. While the absolute normalization of the alpha_co(Tdust) relation is uncertain, the global trend is robust against possible systematic biases in the determination of Mdust, GDR or metallicity. Although we cannot formally distinguish between a step and a smooth evolution of alpha_co with the dust temperature, we can conclude that in galaxies of near-solar metallicity, a critical value of Tdust=30K can be used to determine whether the appropriate alpha_co is closer to the starburst value (1.0 Msun(K kms pc^2)^-1, if Tdust>30K) or closer to the Galactic value (4.35 Msun (K kms pc^2)^-1, if Tdust<30K). This indicator has the great advantage of being less subjective than visual morphological classifications of mergers/SFGs, which can be difficult at high z because of the clumpy nature of SFGs. In the absence of far-infrared data, the offset of a galaxy from the main sequence (i.e., log[SSFR(galaxy)/SSFR_MS(M*,z)]) can be used to identify galaxies requiring the use of an alpha_co conversion factor lower than the Galactic value.

preprint2012arXiv

Excited OH+, H2O+, and H3O+ in NGC 4418 and Arp 220

We report on Herschel/PACS observations of absorption lines of OH+, H2O+ and H3O+ in NGC 4418 and Arp 220. Excited lines of OH+ and H2O+ with E_lower of at least 285 and \sim200 K, respectively, are detected in both sources, indicating radiative pumping and location in the high radiation density environment of the nuclear regions. Abundance ratios OH+/H2O+ of 1-2.5 are estimated in the nuclei of both sources. The inferred OH+ column and abundance relative to H nuclei are (0.5-1)x10^{16} cm-2 and \sim2x10^{-8}, respectively. Additionally, in Arp 220, an extended low excitation component around the nuclear region is found to have OH+/H2O+\sim5-10. H3O+ is detected in both sources with N(H3O+)\sim(0.5-2)x10^{16} cm-2, and in Arp 220 the pure inversion, metastable lines indicate a high rotational temperature of ~500 K, indicative of formation pumping and/or hot gas. Simple chemical models favor an ionization sequence dominated by H+ - O+ - OH+ - H2O+ - H3O+, and we also argue that the H+ production is most likely dominated by X-ray/cosmic ray ionization. The full set of observations and models leads us to propose that the molecular ions arise in a relatively low density (\gtrsim10^4 cm-3) interclump medium, in which case the ionization rate per H nucleus (including secondary ionizations) is zeta>10^{-13} s-1, a lower limit that is severalx10^2 times the highest rate estimates for Galactic regions. In Arp 220, our lower limit for zeta is compatible with estimates for the cosmic ray energy density inferred previously from the supernova rate and synchrotron radio emission, and also with the expected ionization rate produced by X-rays. In NGC 4418, we argue that X-ray ionization due to an AGN is responsible for the molecular ion production.

preprint2012arXiv

Herschel-PACS Observations of Far-IR CO Line Emission in NGC 1068: Highly Excited Molecular Gas in the Circumnuclear Disk

We report the detection of far-IR CO rotational emission from the prototypical Seyfert 2 galaxy NGC 1068. Using Herschel-PACS, we have detected 11 transitions in the J_upper=14-30 (E_upper/k_B = 580-2565 K) range, all of which are consistent with arising from within the central 10" (700 pc). The detected transitions are modeled as arising from 2 different components: a moderate excitation (ME) component close to the galaxy systemic velocity, and a high excitation (HE) component that is blueshifted by ~80 km s^{-1}. We employ a large velocity gradient (LVG) model and derive n_H2~10^{5.6} cm^{-3}, T_kin~170 K, and M_H2~10^{6.7} M_sun for the ME component, and n_H2~10^{6.4} cm^{-3}, T_kin~570 K, and M_H2~10^{5.6} M_sun for the HE component, although for both components the uncertainties in the density and mass are plus/minus (0.6-0.9) dex. We compare the CO line profiles with those of other molecular tracers observed at higher spatial and spectral resolution, and find that the ME transitions are consistent with these lines arising in the ~200 pc diameter ring of material traced by H_2 1-0 S(1) observations. The blueshift of the HE lines may also be consistent with the bluest regions of this H_2 ring, but a better kinematic match is found with a clump of infalling gas ~40 pc north of the AGN. We discuss the prospects of placing the HE component near the AGN, and conclude that while the moderate thermal pressure precludes an association with the ~1 pc radius H_2O maser disk, the HE component could potentially be located only a few parsecs more distant from the AGN, and might then provide the N_H~10^{25} cm^{-2} column obscuring the nuclear hard X-rays. Finally, we also report sensitive upper limits extending up to J_upper=50, which place constraints on a previous model prediction for the CO emission from the X-ray obscuring torus. [Abridged]

preprint2012arXiv

Herschel/PACS spectroscopy of NGC 4418 and Arp 220: H2O, H2^{18}O, OH, ^{18}OH, O I, HCN and NH3

Herschel/PACS spectroscopy of the luminous infrared galaxies NGC4418 and Arp220 reveals high excitation in H2O, OH, HCN, and NH3. In NGC4418, absorption lines were detected with E_low>800 K (H2O), 600 K (OH), 1075 K (HCN), and 600 K (NH3), while in Arp220 the excitation is somewhat lower. While outflow signatures in moderate excitation lines are seen in Arp220 as reported in previous studies, in NGC4418 the lines tracing its outer regions are redshifted relative to the nucleus, suggesting an inflow with Mdot<~12 Msun yr^{-1}. Both galaxies have warm (Tdust>~100 K) nuclear continuum components, together with a more extended component that is much more prominent and massive in Arp220. A chemical dichotomy is found in both sources: on the one hand, the nuclear regions have high H2O abundances, ~10^{-5}, and high HCN/H2O and HCN/NH3 column density ratios of 0.1-0.4 and 2-5, respectively, indicating a chemistry typical of evolved hot cores where grain mantle evaporation has occurred. On the other hand, the high OH abundance, with OH/H2O ratios of ~0.5, indicates the effects of X-rays and/or cosmic rays. The nuclear media have surface brightnesses >~10^{13} Lsun/kpc^2 and are estimated to be thick (N_H>~10^{25} cm^{-2}). While NGC4418 shows weak absorption in H2^{18}O and ^{18}OH, with a ^{16}O-to-^{18}O ratio of >~250-500, the strong absorption of the rare isotopologues in Arp220 indicates ^{16}O-to-^{18}O of 70-130. Further away from the nuclear regions, the H2O abundance decreases to <~10^{-7} and the OH/H2O ratio is 2.5-10. Despite the different scales of NGC4418, Arp220, and Mrk231, preliminary evidence is found for an evolutionary sequence from infall, hot-core like chemistry, and solar oxygen isotope ratio to high velocity outflow, disruption of the hot core chemistry and cumulative high mass stellar processing of 18O.

preprint2012arXiv

Spectroscopic FIR mapping of the disk and galactic wind of M82 with Herschel-PACS

[Abridged] We present maps of the main cooling lines of the neutral atomic gas ([OI] at 63 and 145 micron and [CII] at 158 micron) and in the [OIII] 88 micron line of the starburst galaxy M82, carried out with the PACS spectrometer on board the Herschel satellite. By applying PDR modeling we derive maps of the main ISM physical parameters, including the [CII] optical depth, at unprecedented spatial resolution (~300 pc). We can clearly kinematically separate the disk from the outflow in all lines. The [CII] and [OI] distributions are consistent with PDR emission both in the disk and in the outflow. Surprisingly, in the outflow, the atomic and the ionized gas traced by the [OIII] line both have a deprojected velocity of ~75 km/s, very similar to the average velocity of the outflowing cold molecular gas (~ 100 km/s) and several times smaller than the outflowing material detected in Halpha (~ 600 km/s). This suggests that the cold molecular and neutral atomic gas and the ionized gas traced by the [OIII] 88 micron line are dynamically coupled to each other but decoupled from the Halpha emitting gas. We propose a scenario where cold clouds from the disk are entrained into the outflow by the winds where they likely evaporate, surviving as small, fairly dense cloudlets (n_H\sim 500-1000 cm^-3, G_0\sim 500- 1000, T_gas\sim300 K). We show that the UV photons provided by the starburst are sufficient to excite the PDR shells around the molecular cores. The mass of the neutral atomic gas in the outflow is \gtrsim 5-12x 10^7 M_sun to be compared with that of the molecular gas (3.3 x 10^8 M_sun) and of the Halpha emitting gas (5.8 x 10^6 M_sun). The mass loading factor, (dM/dt)/SFR, of the molecular plus neutral atomic gas in the outflow is ~ 2. Energy and momentum driven outflow models can explain the data equally well, if all the outflowing gas components are taken into account.

preprint2011arXiv

Far-Infrared Line Deficits in Galaxies with Extreme Lfir/MH2 Ratios

We report initial results from the far-infrared fine structure line observations of a sample of 44 local starbursts, Seyfert galaxies and infrared luminous galaxies obtained with the PACS spectrometer on board Herschel. We show that the ratio between the far-infrared luminosity and the molecular gas mass, Lfir/MH2, is a much better proxy for the relative brightness of the far-infrared lines than Lfir alone. Galaxies with high Lfir/MH2 ratios tend to have weaker fine structure lines relative to their far-infrared continuum than galaxies with Lfir/MH2 < 80 Lsun/Msun. A deficit of the [C II] 158 micron line relative to Lfir was previously found with the ISO satellite, but now we show for the first time that this is a general aspect of all far-infrared fine structure lines, regardless of their origin in the ionized or neutral phase of the interstellar medium. The Lfir/MH2 value where these line deficits start to manifest is similar to the limit that separates between the two modes of star formation recently found in galaxies on the basis of studies of their gas-star formation relations. Our finding that the properties of the interstellar medium are also significantly different in these regimes provides independent support for the different star forming relations in normal disk galaxies and major merger systems. We use the spectral synthesis code Cloudy to model the emission of the lines. The expected increase of the ionization parameter with Lfir/MH2 can simultaneously explain the line deficits in the [C II], [N II] and [O I] lines.

preprint2011arXiv

Massive molecular outflows and negative feedback in ULIRGs observed by Herschel-PACS

Mass outflows driven by stars and active galactic nuclei are a key element in many current models of galaxy evolution. They may produce the observed black hole-galaxy mass relation and regulate and quench both star formation in the host galaxy and black hole accretion. However, observational evidence of such feedback processes through outflows of the bulk of the star forming molecular gas is still scarce. Here we report the detection of massive molecular outflows, traced by the hydroxyl molecule (OH), in far-infrared spectra of ULIRGs obtained with Herschel-PACS as part of the SHINING key project. In some of these objects the (terminal) outflow velocities exceed 1000 km/s, and their outflow rates (up to ~1200 M_sol/yr) are several times larger than their star formation rates. We compare the outflow signatures in different types of ULIRGs and in starburst galaxies to address the issue of the energy source (AGN or starburst) of these outflows. We report preliminary evidence that ULIRGs with a higher AGN luminosity (and higher AGN contribution to L_IR) have higher terminal velocities and shorter gas depletion time scales. The outflows in the observed ULIRGs are able to expel the cold gas reservoirs from the centres of these objects within ~1E6-1E8 years.

preprint2010arXiv

Herschel-PACS spectroscopy of IR-bright galaxies at high redshift

We present Herschel-PACS observations of rest-frame mid-infrared and far-infrared spectral line emissions from two lensed, ultra-luminous infrared galaxies at high redshift: MIPS J142824.0+352619 (MIPS J1428), a starburst-dominated system at z = 1.3, and IRAS F10214+4724 (F10214), a source at z = 2.3 hosting both star-formation and a luminous AGN. We have detected [OI]63 micron and [OIII]52 micron in MIPS J1428, and tentatively [OIII]52 micron in F10214. Together with the recent ZEUS-CSO [CII]158 micron detection in MIPS J1428 we can for the first time combine [OI], [CII] and far-IR (FIR) continuum measurements for photo-dissociation (PDR) modeling of an ultra-luminous (L_IR > 10^12 L_sun) star forming galaxy at the peak epoch of cosmic star formation. We find that MIPS J1428, contrary to average local ULIRGs, does not show a deficit in [OI] relative to FIR. The combination of far-UV flux G_0 and gas density n (derived from the PDR models), as well as the star formation efficiency (derived from CO and FIR) is similar to normal or starburst galaxies, despite the high infrared luminosity of this system. In contrast, F10214 has stringent upper limits on [OIV] and [SIII], and an [OIII]/FIR ratio at least an order of magnitude lower than local starbursts or AGN, similar to local ULIRGs.