Researcher profile

Miguel Pereira-Santaella

Miguel Pereira-Santaella contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
30works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

30 published item(s)

preprint2026arXiv

Evidence of Feedback Effects in Low-luminosity Active Galactic Nuclei Revealed by JWST Spectroscopy

This letter presents an analysis of the infrared ($\sim 3-28\,μm$) spectra extracted from the nuclear ($r < 150$ pc) regions of four low-luminosity active galactic nuclei (AGN), observed by JWST NIRSpec/IFU and MIRI/MRS as an extension of the Galaxy Activity, Torus, and Outflow Survey (GATOS). We find that, compared to higher-luminosity AGN, these low-luminosity AGN exhibit distinct properties in their emission of ionized gas, polycyclic aromatic hydrocarbons (PAHs), and molecular hydrogen (H$_2$). Specifically, the low-luminosity AGN exhibit relatively weak high-ionization potential lines (e.g., [Ne V] and [O IV]), and the line ratios suggest that fast radiative shocks (with $v_{\rm s}$ of $\sim \rm 100s\,km\,s^{-1}$) are the primary excitation source of ionized gas therein. Under the low-excitation conditions of their nuclear regions, these low-luminosity AGN generally exhibit a higher fraction of PAHs with large size ($N_{\rm C} \gtrsim 200$), reflecting the preferential destruction of smaller PAH molecules by AGN feedback. Furthermore, the H$_2$ transitions in these low-luminosity AGN are not fully thermalized, with slow, plausibly jet-driven molecular shocks (with $v_{\rm s} \leq \rm 10\,km\,s^{-1}$) likely being the extra excitation source. Taken together with results from the literature, these findings indicate that feedback operates in both low- and high-luminosity AGN, albeit its impact varies with AGN luminosity. In particular, systematic variations in PAH band ratios are found across AGN, demonstrating the differing influence of feedback in AGN of varying luminosities and highlighting the potential of PAH band ratios as diagnostics for distinguishing kinetic- and radiative-mode AGN feedback.

preprint2026arXiv

Measuring the Central Dark Mass in NGC 4258 with JWST/NIRSpec Stellar Kinematics

We present a new stellar dynamical measurement of the supermassive black hole (SMBH) mass in the nearby spiral galaxy NGC 4258, a critical benchmark for extragalactic mass measurements. We use archival JWST/NIRSpec IFU data (G235H/F170LP grating) to extract high-resolution two-dimensional stellar kinematics from the CO bandhead absorption features within the central $3&#39;&#39; \times 3&#39;&#39;$. We extract the stellar kinematics after correcting for instrumental artifacts and separating the stellar light from the non-thermal AGN continuum. We employ Jeans Anisotropic Models (JAM) to fit the observed kinematics, exploring a grid of 12 models to systematically test the impact of different assumptions for the point-spread function, stellar mass-to-light ratio ($M/L$) profile, and orbital anisotropy. All 12 models provide broadly acceptable fits, albeit with minor differences. The ensemble median and 68% (1$σ$) bootstrap confidence intervals of our 12 models yield a black hole mass of $M_{\rm BH} = (4.08^{+0.19}_{-0.33}) \times 10^7$ M$_\odot$. This paper showcases the utility of using the full model ensemble to robustly account for systematic uncertainties, rather than relying on formal errors from a single preferred model, as has been common practice. Our result is just 5% larger than, and consistent with, the benchmark SMBH mass derived from water maser dynamics, validating the use of NIRSpec stellar kinematics for robust SMBH mass determination. Our analysis demonstrates JWST&#39;s capability to resolve the SMBH&#39;s sphere of influence and deliver precise dynamical masses, even in the presence of significant AGN continuum emission.

preprint2025arXiv

Galaxy Activity, Torus and Outflow Survey (GATOS) X: Molecular gas clumpiness under the influence of AGN

The distribution of molecular gas on small scales regulates star formation and the growth of supermassive black holes in galaxy centers, yet the role of active galactic nuclei (AGN) feedback in shaping this distribution remains poorly constrained. We investigate how AGN influence the small-scale structure of molecular gas in galaxy centers, by measuring the clumpiness of CO(3 - 2) emission observed with the Atacama Large Millimeter/submillimeter Array (ALMA) in the nuclear regions (50 - 200 pc from the AGN) of 16 nearby Seyfert galaxies from the Galaxy Activity, Torus, and Outflow Survey (GATOS). To quantify clumpiness, we apply three different methods: (1) the median of the pixel-by-pixel contrast between the original and smoothed maps; (2) the ratio of the total excess flux to the total flux, after substracting the background smoothed emission; and (3) the fraction of total flux coming from clumpy regions, interpreted as the mass fraction in clumps. We find a negative correlation between molecular gas clumpiness and AGN X-ray luminosity (L_X), suggesting that higher AGN activity is associated with smoother gas distributions. All methods reveal a turnover in this relation around L_X = 10^{42} erg/s, possibly indicating a threshold above which AGN feedback becomes efficient at dispersing dense molecular structures and suppressing future star formation. Our findings provide new observational evidence that AGN feedback can smooth out dense gas structures in galaxy centers.

preprint2025arXiv

Impact of AGN and nuclear star formation on the ISM turbulence of galaxies: Insights from JWST/MIRI spectroscopy

Active galactic nuclei (AGN), star formation (SF), and galaxy interactions can drive turbulence in the gas of the ISM, which in turn plays a role in the SF within galaxies. The impact on molecular gas is of particular importance, as it serves as the primary fuel for SF. Our goal is to investigate the origin of turbulence and the emission of molecular gas, as well as low- and intermediate-ionization gas, in the inner few kpc of both AGN hosts and SF galaxies. We use JWST MIRI/MRS observations of a sample consisting of 54 galaxies at z<0.1. We present fluxes of the H2 S(5)6.9091, [Ar II]6.9853, [FeII]5.3403, and [Ar III]8.9914 lines, along with velocity dispersion from W80. For galaxies with coronal emission, [Mg V]5.6098 is also included. Line ratios are compared to photoionization and shock models to explore the origin of the gas emission. AGNs exhibit broader emission lines than SFGs, with the largest velocity dispersions observed in radio-strong (RS) AGNs. H2 gas is less turbulent compared to ionized gas, while coronal gas presents higher velocity dispersions. The W80 values for the ionized gas exhibits a decrease from the nucleus out to radii of approximately 0.5--1 kpc, followed by an outward increase up to 2-3 kpc. In contrast, the H2 line widths generally display increasing profiles with distance from the center. Correlations W80 and line ratios such as H2 S(5)/[ArII] and [FeII]/[ArII] indicate that the most turbulent gas is associated with shocks, enhancing H2 and [FeII] emissions. We speculate that these shocked gas regions are produced by AGN outflows and jet-cloud interactions in AGN-dominated sources, while in SFGs, they may be created by stellar winds and mergers. This shock-induced gas heating may be an important mechanism of AGN (or stellar) feedback, preventing the gas from cooling and forming new stars.

preprint2022arXiv

On the Viability of Determining Galaxy Properties from Observations I: Star Formation Rates and Kinematics

We explore how observations relate to the physical properties of the emitting galaxies by post-processing a pair of merging $z\sim2$ galaxies from the cosmological, hydrodynamical simulation NewHorizon using LCARS (Light from Cloudy Added to RAMSES) to encode the physical properties of the simulated galaxy into H$α$ emission line. By carrying out mock observations and analysis on these data cubes we ascertain which physical properties of the galaxy will be recoverable with the HARMONI spectrograph on the European Extremely Large Telescope (ELT). We are able to estimate the galaxy&#39;s star formation rate and dynamical mass to a reasonable degree of accuracy, with values within a factor of $1.81$ and $1.38$ of the true value. The kinematic structure of the galaxy is also recovered in mock observations. Furthermore, we are able to recover radial profiles of the velocity dispersion and are therefore able to calculate how the dynamical ratio varies as a function of distance from the galaxy centre. Finally, we show that when calculated on galaxy scales the dynamical ratio does not always provide a reliable measure of a galaxy&#39;s stability against gravity or act as an indicator of a minor merger.

preprint2021arXiv

HARMONI view of the host galaxies of active galactic nuclei around cosmic noon: Resolved stellar morpho-kinematics and the M$_{BH}$-$σ_{\star}$ relation

Context. The formation and evolution of galaxies appear linked to the growth of supermassive black holes, as evidenced by empirical scaling relations in nearby galaxies. Understanding this co-evolution over cosmic time requires the revelation of the dynamical state of galaxies and the measurement of the mass of their central black holes (MBH) at a range of cosmic distances. Bright active galactic nuclei (AGNs) are ideal for this purpose. Aims. The High Angular Resolution Monolithic Optical and Near-infrared Integral field spectrograph (HARMONI), the first light integral-field spectrograph for the Extremely Large Telescope (ELT), will transform visible and near-infrared ground-based astrophysics thanks to its advances in sensitivity and angular resolution. We aim to analyse the capabilities of HARMONI to reveal the stellar morpho-kinematic properties of the host galaxies of AGNs at about cosmic noon. Methods. We made use of the simulation pipeline for HARMONI (HSIM) to create mock observations of representative AGN host galaxies at redshifts around cosmic noon. We used observations taken with the Multi-Unit Spectroscopic Explorer (MUSE) of nearby galaxies showing different morphologies and dynamical stages combined with theoretical AGN spectra to create the target inputs for HSIM. Results. According to our simulations, an on-source integration time of three hours should be enough to measure the MBH and to trace the morphology and stellar kinematics of the brightest host galaxies of AGNs beyond cosmic noon. For host galaxies with stellar masses $<$10$^{11}$ solar masses, longer exposure times are mandatory to spatially resolve the stellar kinematics.

preprint2021arXiv

Predicting the Observability of Population III Stars with ELT-HARMONI via the Helium $1640{\rmÅ}$ emission line

Population III (Pop. III) stars, as of yet, have not been detected, however as we move into the era of extremely large telescopes this is likely to change. One likely tracer for Pop. III stars is the HeII$\lambda1640$ emission line, which will be detectable by the HARMONI spectrograph on the European Extremely Large Telescope (ELT) over a broad range of redshifts ($2\leq z\leq14$). By post-processing galaxies from the cosmological, AMR-hydrodynamical simulation NewHorizon with theoretical spectral energy distributions (SED) for Pop. III stars and radiative transfer (i.e. the Yggdrasil Models and CLOUDY look-up tables respectively) we are able to compute the flux of HeII$\lambda1640$ for individual galaxies. From mock 10 hour observations of these galaxies we show that HARMONI will be able to detect Pop. III stars in galaxies up to $z\sim10$ provided Pop. III stars have a top heavy Initial Mass Function (IMF). Furthermore, we find that should Pop. III stars instead have an IMF similar to those of the Pop. I stars, the HeII$\lambda1640$ line would only be observable for galaxies with Pop. III stellar masses in excess of $10^{7}{\,M}_\odot\,{\rm yr}^{-1}$, average stellar age $<1{\rm Myr}$ at $z=4$. Finally, we are able to determine the minimal intrinsic flux required for HARMONI to detect Pop. III stars in a galaxy up to $z=10$.

preprint2021arXiv

SOFIA observations of far-IR fine-structure lines in galaxies to measure metallicity

We present new and archival SOFIA FIFI-LS far-IR spectroscopic observations of 25 local galaxies of either the [OIII]52um and/or the [NIII]57um lines. Including other 31 galaxies from Herschel-PACS, we discuss a local sample of 47 galaxies, including HII region, luminous IR, low-metallicity dwarf and Seyfert galaxies. Analyzing the mid- to far-IR fine-structure lines of this sample, we assess the metallicity and compare with the optical spectroscopy estimates. Using the IR, we find a similar O/H--N/O relation to that known in the optical. As opposite, we find systematically lower N/O IR abundances when compared to the optical determinations, especially at high values of N/O (log(N/O) > -0.8). We explore various hypotheses to account for this difference: (i) difference in ionization structure traced by optical (O+, N+ regions) versus IR lines (O++, N++ regions); (ii) contamination of diffuse ionized gas affecting the optical lines used to compute the N/O abundance; (iii) dust obscuration affecting the optical-based determinations. However, we have not found any correlation of the Delta(N/O)= (N/O)_OPT-(N/O)_IR with either ionization, or electron density, or optical extinction. We speculatively suggest that accretion of metal-poor gas from the circumgalactic medium could provide an explanation for this difference, because the rapid decrease of total abundances during infall is followed by a N/O ratio decrease due to primary production of young - possibly embedded - massive stars, are preferentially traced by the IR diagnostics, while optical diagnostics would better trace the secondary production, when both N/O and O/H abundance ratios will increase.

preprint2021arXiv

Subarcsecond LOFAR imaging of Arp299 at 150 MHz. Tracing the nuclear and diffuse extended emission of a bright LIRG

We study for the first time the low-frequency ($\sim$150 MHz) radio brightness distribution of Arp~299 at subarcsecond resolution, tracing in both compact and extended emission regions the local spectral energy distribution (SED) in order to characterize the dominant emission and absorption processes. We analysed the spatially resolved emission of Arp 299 revealed by 150 MHz international baseline Low-Frequency Array (LOFAR) and 1.4, 5.0, and 8.4 GHz Very Large Array (VLA) observations. We present the first subarcsecond (0.4&#34;$\sim$100~pc) image of the whole Arp~299 system at 150~MHz. The high surface brightness sensitivity of our LOFAR observations ($\sim$100 $μ$Jy/beam) allowed us to detect all of the nuclear components detected at higher frequencies, as well as the extended steep-spectrum emission surrounding the nuclei. We obtained spatially resolved, two-point spectral index maps for the whole galaxy: the compact nuclei show relatively flat spectra, while the extended, diffuse component shows a steep spectrum. We fitted the radio SED of the nuclear regions using two different models: a continuous free-free medium model and a clumpy model. The continuous model can explain the SED of the nuclei assuming a population of relativistic electrons subjected to synchrotron, bremsstrahlung, and ionization losses. The clumpy model fits assuming relativistic electrons with negligible energy losses, and thermal fractions that are more typical of star-forming galaxies than those required for the continuous model. Our results confirm the usefulness of combining spatially resolved radio imaging at both MHz and GHz frequencies to characterize in detail the radio emission properties of LIRGs from the central 100 pc out to the kiloparsec galaxy-wide scales.

preprint2020arXiv

A proto-pseudobulge in ESO 320-G030 fed by a massive molecular inflow driven by a nuclear bar

Galaxies with nuclear bars are believed to efficiently drive gas inward, generating a nuclear starburst and possibly an active galactic nucleus (AGN). We confirm this scenario for the isolated, double-barred, luminous infrared galaxy ESO 320-G030 based on an analysis of Herschel and ALMA spectroscopic observations. Herschel/PACS and SPIRE observations of ESO 320-G030 show absorption/emission in 18 lines of H2O, which we combine with the ALMA H2O 423-330 448 GHz line (Eupper~400 K) and continuum images to study the nuclear region. Radiative transfer models indicate that 3 nuclear components are required to account for the H2O and continuum data. An envelope, with R~130-150 pc, T_dust~50 K, and N_H2~2x10^{23} cm^{-2}, surrounds a nuclear disk with R~40 pc and tau_100um~1.5-3 (N_H2~2x10^{24} cm^{-2}) and an extremely compact (R~12 pc), warm (~100 K), and buried (tau_100um>5, N_H2>~5x10^{24} cm^{-2}) core component. The three nuclear components account for 70% of the galaxy L_IR (SFR~16-18 Msun yr^{-1}). The nucleus is fed by a molecular inflow observed in CO 2-1 with ALMA, which is associated with the nuclear bar. With decreasing radius (r=450-225 pc), the mass inflow rate increases up to ~20 Msun yr^{-1}, which is similar to the nuclear SFR, indicating that the starburst is sustained by the inflow. At lower r, the inflow is best probed by the far-infrared OH ground-state doublets, with an estimated inflow rate of ~30 Msun yr^{-1}. The short timescale of ~20 Myr for nuclear gas replenishment indicates quick secular evolution, and indicates that we are witnessing an intermediate stage (<100 Myr) proto-pseudobulge fed by a massive inflow that is driven by a strong nuclear bar. We also apply the H2O model to the Herschel far-infrared spectroscopic observations of H2^{18}O, OH, $^{18}OH, OH+, H2O^+, H3O^+, NH, NH2, NH3, CH, CH^+, ^{13}CH^+, HF, SH, and C3, and estimate their abundances.

preprint2020arXiv

Extinction in the 11.2 micron PAH band and the low L_11.2/L_IR in ULIRGs

We present a method for recovering the intrinsic (extinction-corrected) luminosity of the 11.2 micron PAH band in galaxy spectra. Using 105 high S/N Spitzer/IRS spectra of star-forming galaxies, we show that the equivalent width ratio of the 12.7 and 11.2 micron PAH bands is independent on the optical depth, with small dispersion of ~5% indicative of a nearly constant intrinsic flux ratio R_int = (f_12.7/f_11.2)_int = 0.377 +/- 0.020. Conversely, the observed flux ratio, R_obs = (f_12.7/f_11.2)_obs strongly correlates with the silicate strength (S_sil) confirming that differences in R_obs reflect variation in the optical depth. The relation between R_obs and S_sil reproduces predictions for the Galactic Centre extinction law but disagrees with other laws. We calibrate the total extinction affecting the 11.2 micron PAH from R_obs, which we apply to another sample of 215 galaxies with accurate measurements of the total infrared luminosity (L_IR) to investigate the impact of extinction on L_11.2/L_IR. Correlation between L_11.2/L_IR and R_obs independently on L_IR suggests that increased extinction explains the well known decrease in the average L_11.2/L_IR at high L_IR. The extinction-corrected L_11.2 is proportional to L_IR in the range L_IR/L_sun = 10^9--10^13. These results consolidate L_11.2 as a robust tracer of star formation in galaxies.

preprint2020arXiv

Simulating gas kinematic studies of high-redshift galaxies with the HARMONI Integral Field Spectrograph

We present simulated observations of gas kinematics in galaxies formed in 10 pc resolution cosmological simulations with the hydrodynamical + N-body code RAMSES, using the new RAMSES2HSIM pipeline with the simulated observing pipeline (HSIM) for the ELT HARMONI IFU spectrograph. We post-process the galaxy&#39;s gas kinematics and Hα line emission for each simulation cell, and integrate the emission to produce an extinction-corrected input cube. We then simulate observations of the input cube with HARMONI, for a range of exposure times, spatial sampling, and spectral resolution. We analyze the mock observations to recover galaxy properties such as its kinematics and compare with the known simulation values. We investigate the cause of biases between the &#39;real&#39; and &#39;observed&#39; kinematic values, demonstrating the sensitivity of the inferred rotation curve to knowledge of the instrument&#39;s point spread function.

preprint2016arXiv

VLT-SINFONI sub-kpc study of the star formation in local LIRGs and ULIRGs: Analysis of the global $Σ_{\rm SFR}$ structure and characterisation of individual star-forming clumps

We present a 2D study of star formation at kpc and sub-kpc scales of a sample of local (z<0.1) U/LIRGs, based on near-IR VLT-SINFONI observations. We obtained integrated measurements of the star formation rate (SFR) and star formation rate surface density, together with their 2D distributions, based on Br_gamma and Pa_alpha emission. We observe a tight linear correlation between the SFR derived from our extinction-corrected measurements and that derived from 24 micron data, and a reasonable agreement with SFR derived from total IR luminosity. Our near-IR measurements are on average a factor 3 larger than optical SFR, even when extinction corrections are applied. We found that LIRGs have a median-observed star formation rate surface density of 1.72 Msun/yr/kpc^2 for the extinction-corrected distribution, whilst ULIRGs have 0.23 Msun/yr/kpc^2, respectively. These median values for ULIRGs increase up to 2.90 Msun/yr/kpc^2, when only their inner regions, covering the same size as the average FoV of LIRGs, are considered. We identified a total of 95 individual SF clumps in our sample, with sizes within 60-1500pc, and extinction-corrected Pa_alpha luminosities of 10^5-10^8 Lsun. Star-forming clumps in LIRGs are about ten times larger and thousands of times more luminous than typical clumps in spiral galaxies. Clumps in ULIRGs have sizes similar (x0.5-1) to those of high-z clumps, having Pa_alpha luminosities similar to some high-z clumps, and about 10 times less luminous than the most luminous high-z clumps identified so far. We also observed a change in the slope of the L-r relation. A likely explanation is that most luminous galaxies are interacting and merging, and therefore their size represents a combination of the distribution of the star-forming clumps within each galaxy in the system plus the effect of the projected distance.

preprint2015arXiv

Star-formation histories of local luminous infrared galaxies

We present the analysis of the integrated spectral energy distribution (SED) from the ultraviolet (UV) to the far-infrared and H$α$ of a sample of 29 local systems and individual galaxies with infrared (IR) luminosities between 10^11 Lsun and 10^11.8 Lsun. We have combined new narrow-band H$α$+[NII] and broad-band g, r optical imaging taken with the Nordic Optical Telescope (NOT), with archival GALEX, 2MASS, Spitzer, and Herschel data. The SEDs (photometry and integrated H$α$ flux) have been fitted with a modified version of the MAGPHYS code using stellar population synthesis models for the UV-near-IR range and thermal emission models for the IR emission taking into account the energy balance between the absorbed and re-emitted radiation. From the SED fits we derive the star-formation histories (SFH) of these galaxies. For nearly half of them the star-formation rate appears to be approximately constant during the last few Gyrs. In the other half, the current star-formation rate seems to be enhanced by a factor of 3-20 with respect to that occured ~1 Gyr ago. Objects with constant SFH tend to be more massive than starbursts and they are compatible with the expected properties of a main-sequence (M-S) galaxy. Likewise, the derived SFHs show that all our objects were M-S galaxies ~1 Gyr ago with stellar masses between 10^10.1 and 10^11.5 Msun. We also derived from our fits the average extinction (A_v=0.6-3 mag) and the polycyclic aromatic hydrocarbons (PAH) luminosity to L(IR) ratio (0.03-0.16). We combined the A_v with the total IR and H$α$ luminosities into a diagram which can be used to identify objects with rapidly changing (increasing or decreasing) SFR during the last 100 Myr.

preprint2014arXiv

Far-infrared line spectra of Seyfert galaxies from the Herschel-PACS Spectrometer

We present spectroscopic observations of FIR fine-structure lines of 26 Seyfert galaxies obtained with the Herschel-PACS spectrometer. These observations are complemented by spectroscopy with Spitzer-IRS and Herschel-SPIRE. The ratios of the OIII, NII, SIII and NeV lines have been used to determine electron densities in the ionised gas regions. The CI lines, observed with SPIRE, have been used to measure the densities in the neutral gas, while the OI lines provide a measure of the gas temperature, at densities below 10000 cm-3. Using the OI145/63um and SIII33/18um line ratios we find an anti-correlation of the temperature with the gas density. Using various fine-structure line ratios, we find that density stratification is common in these active galaxies. On average, the electron densities increase with the ionisation potential of the ions producing the NII, SIII and NeV emission. The infrared emission lines arise partly in the Narrow Line Region (NLR) photoionised by the AGN central engine, partly in HII regions photo ionised by hot stars and partly in neutral gas in photo-dissociated regions (PDRs). We attempt to separate the contributions to the line emission produced in these different regions by comparing our emission line ratios to empirical and theoretical values. In particular, we tried to separate the contribution of AGN and star formation by using a combination of Spitzer and Herschel lines, and we found that, besides the well known mid-IR line ratios, the mixed mid-IR/far-IR line ratio of OIII88um/OIV26um can reliably discriminate the two emission regimes, while the far-IR line ratio of CII157um/OI63um is only able to mildly separate the two regimes. By comparing the observed CII157um/NII205um ratio with photoionisation models, we also found that most of the CII emission in the galaxies we examined is due to PDRs.

preprint2014arXiv

Outflow of hot and cold molecular gas from the obscured secondary nucleus of NGC3256: closing in on feedback physics

The nuclei of merging galaxies are often deeply buried in dense layers of gas and dust. In these regions, gas outflows driven by starburst and AGN activity are believed to play a crucial role in the evolution of these galaxies. However, to fully understand this process it is essential to resolve the morphology and kinematics of such outflows. Using near-IR integral-field spectroscopy obtained with VLT/SINFONI, we detect a kpc-scale structure of high-velocity molecular hydrogen (H2) gas associated with the deeply buried secondary nucleus of the IR-luminous merger NGC3256. We show that this structure is likely the hot component of a molecular outflow, which is detected also in the cold molecular gas by Sakamoto et al. This outflow, with a molecular gas mass of M(H2)~2x10^7 Msun, is among the first to be spatially resolved in both the hot H2 gas with VLT/SINFONI and the cold CO-emitting gas with ALMA. The hot and cold components share a similar morphology and kinematics, with a hot-to-cold molecular gas mass ratio of ~6x10^-5. The high (~100 pc) resolution at which we map the geometry and velocity structure of the hot outflow reveals a biconical morphology with opening angle ~40 deg and gas spread across a FWZI~1200 km/s. Because this collimated outflow is oriented close to the plane of the sky, the molecular gas may reach maximum intrinsic outflow velocities of ~1800 km/s, with an average mass outflow rate of at least ~20 Msun/yr. By modeling the line-ratios of various near-IR H2 transitions, we show that the H2 gas in the outflow is heated through shocks or X-rays to a temperature of ~1900K. The energy needed to drive the outflow is likely provided by a hidden Compton-thick AGN or by the nuclear starburst. We show that the global kinematics of the molecular outflow in NGC3256 mimic those of CO-outflows that have been observed at low spatial resolution in starburst- and active galaxies.

preprint2014arXiv

Spatially resolved physical conditions of molecular gas and potential star formation tracers in M83, revealed by the Herschel SPIRE FTS

Since the launch of the Herschel Space Observatory, our understanding about the photo-dissociation regions (PDR) has taken a step forward. In the bandwidth of the Fourier Transform Spectrometer (FTS) of the Spectral and Photometric Imaging REceiver (SPIRE) on board Herschel, ten CO rotational transitions, including J=4-3 to J=13-12, and three fine structure lines, including [CI] 609, [CI] 370, and [NII] 250 micron, are covered. In this paper, we present our findings from the FTS observations at the nuclear region of M83, based on the spatially resolved physical parameters derived from the CO spectral line energy distribution (SLED) map and the comparisons with the dust properties and star-formation tracers. We discuss (1) the potential of using [NII] 250 and [CI] 370 micron as star-formation tracers; (2) the reliability of tracing molecular gas with CO; (3) the excitation mechanisms of warm CO; (4) the possibility of studying stellar feedback by tracing the thermal pressure of molecular gas in the nuclear region of M83.

preprint2014arXiv

Warm molecular gas temperature distribution in six local infrared bright Seyfert galaxies

We simultaneously analyze the spectral line energy distributions (SLEDs) of CO and H2 of six local luminous infrared (IR) Seyfert galaxies. For the CO SLEDs, we used new Herschel/SPIRE FTS data (from J=4-3 to J=13-12) and ground-based observations for the lower-J CO transitions. The H2 SLEDs were constructed using archival mid-IR Spitzer/IRS and near-IR VLT/SINFONI data for the rotational and ro-vibrational H2 transitions, respectively. In total, the SLEDs contain 26 transitions with upper level energies between 5 and 15000 K. A single, constant density, model (n$_{H_2}$ ~ 10$^{4.5-6}$ cm$^{-3}$) with a broken power-law temperature distribution reproduces well both the CO and H2 SLEDs. The power-law indices are $β_1$ ~ 1-3 for warm molecular gas (20 K < T < 100 K) and $β_2$ ~ 4-5 for hot molecular gas (T > 100 K). We show that the steeper temperature distribution (higher $β$) for hot molecular gas can be explained by shocks and photodissociation region (PDR) models, however, the exact $β$ values are not reproduced by PDR or shock models alone and a combination of both is needed. We find that the three major mergers among our targets have shallower temperature distributions for warm molecular gas than the other three spiral galaxies. This can be explained by a higher relative contribution of shock excitation, with respect to PDR excitation, for the warm molecular gas in these mergers. For only one of the mergers, IRASF 05189-2524, the shallower H2 temperature distribution differs from that of the spiral galaxies. The presence of a bright active galactic nucleus in this source might explain the warmer molecular gas observed.

preprint2013arXiv

Herschel SPIRE-FTS Observations of Excited CO and [CI] in the Antennae (NGC 4038/39): Warm and Cold Molecular Gas

We present Herschel SPIRE-FTS observations of the Antennae (NGC 4038/39), a well studied, nearby ($22$ Mpc) ongoing merger between two gas rich spiral galaxies. We detect 5 CO transitions ($J=4-3$ to $J=8-7$), both [CI] transitions and the [NII]$205μm$ transition across the entire system, which we supplement with ground based observations of the CO $J=1-0$, $J=2-1$ and $J=3-2$ transitions, and Herschel PACS observations of [CII] and [OI]$63μm$. Using the CO and [CI] transitions, we perform both a LTE analysis of [CI], and a non-LTE radiative transfer analysis of CO and [CI] using the radiative transfer code RADEX along with a Bayesian likelihood analysis. We find that there are two components to the molecular gas: a cold ($T_{kin}\sim 10-30$ K) and a warm ($T_{kin} \gtrsim 100$ K) component. By comparing the warm gas mass to previously observed values, we determine a CO abundance in the warm gas of $x_{CO} \sim 5\times 10^{-5}$. If the CO abundance is the same in the warm and cold gas phases, this abundance corresponds to a CO $J=1-0$ luminosity-to-mass conversion factor of $α_{CO} \sim 7 \ M_{\odot}{pc^{-2} \ (K \ km \ s^{-1})^{-1}}$ in the cold component, similar to the value for normal spiral galaxies. We estimate the cooling from H$_2$, [CII], CO and [OI]$63μm$ to be $\sim 0.01 L_{\odot}/M_{\odot}$. We compare PDR models to the ratio of the flux of various CO transitions, along with the ratio of the CO flux to the far-infrared flux in NGC 4038, NGC 4039 and the overlap region. We find that the densities recovered from our non-LTE analysis are consistent with a background far-ultraviolet field of strength $G_0\sim 1000$. Finally, we find that a combination of turbulent heating, due to the ongoing merger, and supernova and stellar winds are sufficient to heat the molecular gas.

preprint2013arXiv

Herschel/SPIRE Sub-millimeter Spectra of Local Active Galaxies

We present the sub-millimeter spectra from 450 GHz to 1550 GHz of eleven nearby active galaxies observed with the SPIRE Fourier Transform Spectrometer (SPIRE/FTS) onboard Herschel. We detect CO transitions from J_up = 4 to 12, as well as the two [CI] fine structure lines at 492 and 809 GHz and the [NII] 461 GHz line. We used radiative transfer models to analyze the observed CO spectral line energy distributions (SLEDs). The FTS CO data were complemented with ground-based observations of the low-J CO lines. We found that the warm molecular gas traced by the mid-J CO transitions has similar physical conditions (n_H2 ~ 10^3.2 - 10^3.9 cm^-3 and T_kin ~ 300 - 800 K) in most of our galaxies. Furthermore, we found that this warm gas is likely producing the mid-IR rotational H2 emission. We could not determine the specific heating mechanism of the warm gas, however it is possibly related to the star-formation activity in these galaxies. Our modeling of the [CI] emission suggests that it is produced in cold (T_kin < 30 K) and dense (n_H2 > 10^3 cm^-3) molecular gas. Transitions of other molecules are often detected in our SPIRE/FTS spectra. The HF J=1-0 transition at 1232 GHz is detected in absorption in UGC05101 and in emission in NGC7130. In the latter, near-infrared pumping, chemical pumping, or collisional excitation with electrons are plausible excitation mechanisms likely related to the AGN of this galaxy. In some galaxies few H2O emission lines are present. Additionally, three OH+ lines at 909, 971, and 1033 GHz are identified in NGC7130.

preprint2013arXiv

Local Luminous Infrared Galaxies. III. Co-evolution of Black Hole Growth and Star Formation Activity?

Local luminous infrared (IR) galaxies (LIRGs) have both high star formation rates (SFR) and a high AGN (Seyfert and AGN/starburst composite) incidence. Therefore, they are ideal candidates to explore the co-evolution of black hole (BH) growth and star formation (SF) activity, not necessarily associated with major mergers. Here, we use Spitzer/IRS spectroscopy of a complete volume-limited sample of local LIRGs (distances of <78Mpc). We estimate typical BH masses of 3x10^7 M_sun using [NeIII]15.56micron and optical [OIII]5007A gas velocity dispersions and literature stellar velocity dispersions. We find that in a large fraction of local LIRGs the current SFR is taking place not only in the inner nuclear ~1.5kpc region, as estimated from the nuclear 11.3micron PAH luminosities, but also in the host galaxy. We next use the ratios between the SFRs and BH accretion rates (BHAR) to study whether the SF activity and BH growth are contemporaneous in local LIRGs. On average, local LIRGs have SFR to BHAR ratios higher than those of optically selected Seyferts of similar AGN luminosities. However, the majority of the IR-bright galaxies in the RSA Seyfert sample behave like local LIRGs. Moreover, the AGN incidence tends to be higher in local LIRGs with the lowest SFRs. All this suggests that in local LIRGs there is a distinct IR-bright star forming phase taking place prior to the bulk of the current BH growth (i.e., AGN phase). The latter is reflected first as a composite and then as a Seyfert, and later as a non-LIRG optically identified Seyfert nucleus with moderate SF in its host galaxy.

preprint2012arXiv

Submillimetre line spectrum of the Seyfert galaxy NGC1068 from the Herschel-SPIRE Fourier Transform Spectrometer

The first complete submillimetre spectrum (190-670um) of the Seyfert 2 galaxy NGC1068 has been observed with the SPIRE Fourier Transform Spectrometer onboard the {\it Herschel} Space Observatory. The sequence of CO lines (Jup=4-13), lines from water, the fundamental rotational transition of HF, two o-H_2O+ lines and one line each from CH+ and OH+ have been detected, together with the two [CI] lines and the [NII]205um line. The observations in both single pointing mode with sparse image sampling and in mapping mode with full image sampling allow us to disentangle two molecular emission components, one due to the compact circum-nuclear disk (CND) and one from the extended region encompassing the star forming ring (SF-ring). Radiative transfer models show that the two CO components are characterized by density of n(H_2)=10^4.5 and 10^2.9 cm^-3 and temperature of T=100K and 127K, respectively. The comparison of the CO line intensities with photodissociation region (PDR) and X-ray dominated region (XDR) models, together with other observational constraints, such as the observed CO surface brightness and the radiation field, indicate that the best explanation for the CO excitation of the CND is an XDR with density of n(H_2) 10^4 cm^-3 and X-ray flux of 9 erg s^-1 cm^-2, consistent with illumination by the active galactic nucleus, while the CO lines in the SF-ring are better modeled by a PDR. The detected water transitions, together with those observed with the \her \sim PACS Spectrometer, can be modeled by an LVG model with low temperature (T_kin \sim 40K) and high density (n(H_2) in the range 10^6.7-10^7.9 cm^-3).

preprint2012arXiv

The NGC7771+NGC7770 Minor Merger: Harassing the Little One?

Numerical simulations of minor mergers, typically having mass ratios greater than 3:1, predict little enhancement in the global star formation activity. However, these models also predict that the satellite galaxy is more susceptible to the effects of the interaction than the primary. We use optical integral field spectroscopy and deep optical imaging to study the NGC7771+NGC7770 interacting system (~10:1 stellar mass ratio) to test these predictions. We find that the satellite galaxy NGC7770 is currently experiencing a galaxy-wide starburst with most of the optical light being from young and post-starburst stellar populations(<1Gyr). This galaxy lies off of the local star-forming sequence for composite galaxies with an enhanced integrated specific star formation rate. We also detect in the outskirts of NGC7770 Halpha emitting gas filaments. This gas appears to have been stripped from one of the two galaxies and is being excited by shocks. All these results are consistent with a minor-merger induced episode(s) of star formation in NGC7770 after the first close passage. Such effects are not observed on the primary galaxy NGC7771.

preprint2011arXiv

Local Luminous Infrared Galaxies. II. AGN Activity from Spitzer/IRS spectra

We quantify the active galactic nucleus (AGN) contribution to the mid-infrared (mid-IR) and the total infrared (IR, 8-1000micron) emission in a complete volume-limited sample of 53 local luminous infrared galaxies (LIRGs). We decompose the Spitzer Infrared Spectrograph (IRS) low-resolution 5-38micron spectra of the LIRGs into AGN and starburst components using clumpy torus models and star-forming galaxy templates, respectively. We find that 50% (25/50) of local LIRGs have an AGN component detected with this method. There is good agreement between these AGN detections through mid-IR spectral decomposition and other AGN indicators, such as the optical spectral class, mid-IR spectral features and X-ray properties. Taking all the AGN indicators together, the AGN detection rate in the individual nuclei of LIRGs is ~62%. The derived AGN bolometric luminosities are in the range L_bol(AGN)=0.4 -50x10^{43} erg/s. The AGN bolometric contribution to the IR luminosities of the galaxies is generally small, with 70% of LIRGs having L_bol(AGN)/L_IR<0.05. Only ~8% of local LIRGs have a significant AGN bolometric contribution L_bol(AGN)/L_IR > 0.25. From the comparison of our results with literature results of ultraluminous infrared galaxies, we confirm that in the local universe the AGN bolometric contribution to the IR luminosity increases with the IR luminosity of the galaxy/system. If we add up the AGN bolometric luminosities we find that AGNs only account for 5%^{+8%}_{-3%} of the total IR luminosity produced by local LIRGs (with and without AGN detections). This proves that the bulk of the IR luminosity of local LIRGs is due to star formation activity. Taking the newly determined IR luminosity density of LIRGs in the local universe, we then estimate an AGN IR luminosity density of Omega_IR(AGN) = 3x10^5 L_sun Mpc^{-3}$ in LIRGs.

preprint2011arXiv

The X-ray emission of local luminous infrared galaxies

We study the X-ray emission of a representative sample of 27 local luminous infrared galaxies (LIRGs). The median IR luminosity of our sample is log L_IR/L_sun = 11.2, thus the low-luminosity end of the LIRG class is well represented. We used new XMM-Newton data as well as Chandra and XMM-Newton archive data. The soft X-ray (0.5-2 keV) emission of most of the galaxies (>80%), including LIRGs hosting a Seyfert 2 nucleus, is dominated by star-formation related processes. These LIRGs follow the star-formation rate (SFR) versus soft X-ray luminosity correlation observed in local starbursts. We find that ~15% of the non-Seyfert LIRGs (3 out of 20) have an excess hard X-ray emission relative to that expected from star-formation that might indicate the presence of an obscured AGN. The rest of the non-Seyfert LIRGs follow the SFR versus hard X-ray (2-10 keV) luminosity correlation of local starbursts. The non-detection of the 6.4 keV Fe K alpha emission line in the non-Seyfert LIRGs allows us to put an upper limit to the bolometric luminosity of an obscured AGN, L_bol <1043 erg s-1 . That is, in these galaxies, if they hosted a low luminosity AGN, its contribution to total luminosity would be less than 10%. Finally we estimate that the AGN contribution to the total luminosity for our sample of local LIRGs is between 7% and 10%.

preprint2010arXiv

An X-ray Study of Local Infrared Bright Galaxies

We are carrying out detailed study of the X-ray and infrared (IR) properties of a sample of local (d < 70 Mpc) luminous infrared galaxies (LIRGs) using XMM-Newton and Spitzer (imaging and spectroscopy). The main goal is to study the extreme processes of star formation and/or active galactic nuclei (AGN) taking place in this cosmologically important class of galaxies. In this proceedings we present the preliminary results obtained from the analysis of the XMM-Newton X-ray images and the X-ray spectral modeling.

preprint2010arXiv

Local Luminous Infrared Galaxies. I. Spatially resolved observations with Spitzer/IRS

We present results from the Spitzer/IRS spectral mapping observations of 15 local luminous infrared galaxies (LIRGs). In this paper we investigate the spatial variations of the mid-IR emission which includes: fine structure lines, molecular hydrogen lines, polycyclic aromatic features (PAHs), continuum emission and the 9.7um silicate feature. We also compare the nuclear and integrated spectra. We find that the star formation takes place in extended regions (several kpc) as probed by the PAH emission as well as the [NeII] and [NeIII] emissions. The behavior of the integrated PAH emission and 9.7um silicate feature is similar to that of local starburst galaxies. We also find that the minima of the [NeIII]/[NeII] ratio tends to be located at the nuclei and its value is lower than that of HII regions in our LIRGs and nearby galaxies. It is likely that increased densities in the nuclei of LIRGs are responsible for the smaller nuclear [NeIII]/[NeII] ratios. This includes the possibility that some of the most massive stars in the nuclei are still embedded in ultracompact HII regions. In a large fraction of our sample the 11.3um PAH emission appears more extended than the dust 5.5um continuum emission. We find a dependency of the 11.3um PAH/7.7 um PAH and [NeII]/11.3um PAH ratios with the age of the stellar populations. Smaller and larger ratios respectively indicate recent star formation. The estimated warm (300 K < T < 1000 K) molecular hydrogen masses are of the order of 10^8 M_Sun, which are similar to those found in ULIRGs, local starbursts and Seyfert galaxies. Finally we find that the [NeII] velocity fields for most of the LIRGs in our sample are compatible with a rotating disk at ~kpc scales, and they are in a good agreement with H-alpha velocity fields.

preprint2010arXiv

The Mid-Infrared High-Ionization Lines from Active Galactic Nuclei and Star-Forming Galaxies

We used Spitzer/IRS spectroscopic data on 426 galaxies including quasars, Seyferts, LINER and HII galaxies to investigate the relationship among the mid-IR emission lines. There is a tight linear correlation between the [Ne V]14.3 um and 24.3 um (97.1 eV) and the [O IV]25.9 um (54.9 eV) high-ionization emission lines. The correlation also holds for these high-ionization emission lines and the [Ne III]15.56 um (41 eV) emission line, although only for active galaxies. We used these correlations to calculate the [Ne III] excess due to star formation in Seyfert galaxies. We also estimated the [O IV] luminosity due to star formation in active galaxies and determined that it dominates the [O IV] emission only if the contribution of the active nucleus to the total luminosity is below 5%. We find that the AGN dominates the [O IV] emission in most Seyfert galaxies, whereas star-formation adequately explains the observed [O IV] emission in optically classified HII galaxies. Finally we computed photoionization models to determine the physical conditions of the narrow line region where these high-ionization lines originate. The estimated ionization parameter range is -2.8 < log U < -2.5 and the total hydrogen column density range is 20 < log nH (cm-2) < 21.

preprint2009arXiv

Local Luminous Infrared Galaxies: Spatially resolved mid-infrared observations with Spitzer/IRS

Luminous Infrared (IR) Galaxies (LIRGs) are an important cosmological class of galaxies as they are the main contributors to the co-moving star formation rate density of the universe at z=1. In this paper we present a GTO Spitzer IRS program aimed to obtain spectral mapping of a sample of 14 local (d<76Mpc) LIRGs. The data cubes map, at least, the central 20arcsec x 20arcsec to 30arcsec x 30arcsec regions of the galaxies, and use all four IRS modules covering the full 5-38micron spectral range. The final goal of this project is to characterize fully the mid-IR properties of local LIRGs as a first step to understanding their more distant counterparts. In this paper we present the first results of this GTO program. The IRS spectral mapping data allow us to build spectral maps of the bright mid-IR emission lines (e.g., [NeII], [NeIII], [SIII], H_2), continuum, the 6.2 and 11.3micron PAH features, and the 9.7micron silicate feature, as well as to extract 1D spectra for regions of interest in each galaxy. The IRS data are used to obtain spatially resolved measurements of the extinction using the 9.7micron silicate feature, and to trace star forming regions using the neon lines and the PAH features. We also investigate a number of AGN indicators, including the presence of high excitation emission lines and a strong dust continuum emission at around 6micron. We finally use the integrated Spitzer/IRS spectra as templates of local LIRGs. We discuss several possible uses for these templates, including the calibration of the star formation rate of IR-bright galaxies at high redshift. We also predict the intensities of the brightest mid-IR emission lines for LIRGs as a function of redshift, and compare them with the expected sensitivities of future space IR missions.

preprint2009arXiv

The Extreme Star Formation Activity of Arp299 Revealed by Spitzer IRS Spectral Mapping

We present Spitzer/IRS spectral mapping observations of the luminous infrared galaxy (LIRG) Arp299 (IC694 + NGC3690) covering the central 45arcsec ~ 9kpc. The integrated mid-IR spectrum of Arp299 is similar to that of local starbursts despite its strongly interacting nature and high infrared luminosity, L_IR ~ 6x10^11 Lsun. This is explained because the star formation (probed by e.g. high [NeIII]15.56micron/[NeII]micron line ratios) is spread across at least 6-8kpc. Moreover, a large fraction of this star formation is taking place in young regions of moderate mid-IR optical depths such as the C+C&#39; complex in the overlap region between the two galaxies and in HII regions in the disks of the galaxies. It is only source A, the nuclear region of IC694, that shows the typical mid-IR characteristics of ultraluminous infrared galaxies (ULIRGs, L_IR > 10^12 Lsun), that is, very compact (less than 1kpc) and dust-enshrouded star formation resulting in a deep silicate feature and moderate equivalent widths of the PAHs. The nuclear region of NGC3690, known as source B1, hosts a low-luminosity AGN and is surrounded by regions of star formation. Although the high excitation [NeV]14.32micron line typical of AGN is not detected in B1, its upper limit is consistent with the value expected from the X-ray luminosity. The AGN emission is detected in the form of a strong hot dust component that accounts for 80-90% of the 6micron luminosity of B1. The similarity between the Arp299 integrated mid-IR spectrum and those of high-z ULIRGs suggests that Arp299 may represent a local example, albeit with lower IR luminosity and possibly higher metallicity, of the star-formation processes occurring at high-z.