Proyecto de Investigación: PID2019-105423GA-I00
Cargando...
Colaboradores
Financiadores
ID
PID2019-105423GA-I00
Autores
Publicaciones
Physics of ULIRGs with MUSE and ALMA: The PUMA project II. Are local ULIRGs powered by AGN? The subkiloparsec view of the 220 GHz continuum
(EDP Sciences, 2021-07-12) Pereira Santaella, M.; Colina, L.; García Burillo, S.; Lamperti, I.; González Alfonso, E.; Perna, M.; Arribas, S.; Alonso Herrero, A.; Aalto, S.; Combes, F.; Labiano, Á.; Piqueras López, J.; Rigopoulou, D.; Van der Werf, P. P.; Comunidad de Madrid; Agencia Estatal de Investigación (AEI); Science and Technology Facilities Council (STFC); Pereira Santaella, M. [0000-0002-4005-9619]; Unidad de Excelencia Científica María de Maeztu Centro de Astrobiología del Instituto Nacional de Técnica Aeroespacial y CSIC, MDM-2017-0737
We analyze new high-resolution (400 pc) ∼220 GHz continuum and CO(2–1) Atacama Large Millimeter Array (ALMA) observations of a representative sample of 23 local (z < 0.165) ultra-luminous infrared systems (ULIRGs; 34 individual nuclei) as part of the “Physics of ULIRGs with MUSE and ALMA” (PUMA) project. The deconvolved half-light radii of the ∼220 GHz continuum sources, rcont, are between < 60 pc and 350 pc (median 80–100 pc). We associate these regions with the regions emitting the bulk of the infrared luminosity (LIR). The good agreement, within a factor of 2, between the observed ∼220 GHz fluxes and the extrapolation of the infrared gray-body as well as the small contributions from synchrotron and free–free emission support this assumption. The cold molecular gas emission sizes, rCO, are between 60 and 700 pc and are similar in advanced mergers and early interacting systems. On average, rCO are ∼2.5 times larger than rcont. Using these measurements, we derived the nuclear LIR and cold molecular gas surface densities (ΣLIR = 1011.5 − 1014.3 L⊙ kpc−2 and ΣH2 = 102.9 − 104.2 M⊙ pc−2, respectively). Assuming that the LIR is produced by star formation, the median ΣLIR corresponds to ΣSFR = 2500 M⊙ yr−1 kpc−2. This ΣSFR implies extremely short depletion times, ΣH2/ΣSFR < 1–15 Myr, and unphysical star formation efficiencies > 1 for 70% of the sample. Therefore, this favors the presence of an obscured active galactic nucleus (AGN) in these objects that could dominate the LIR. We also classify the ULIRG nuclei in two groups: (a) compact nuclei (rcont < 120 pc) with high mid-infrared excess emission (ΔL6−20 μm/LIR) found in optically classified AGN; and (b) nuclei following a relation with decreasing ΔL6−20 μm/LIR for decreasing rcont. The majority, 60%, of the nuclei in interacting systems lie in the low-rcont end (<120 pc) of this relation, while this is the case for only 30% of the mergers. This suggests that in the early stages of the interaction, the activity occurs in a very compact and dust-obscured region while, in more advanced merger stages, the activity is more extended, unless an optically detected AGN is present. Approximately two-thirds of the nuclei have nuclear radiation pressures above the Eddington limit. This is consistent with the ubiquitous detection of massive outflows in local ULIRGs and supports the importance of the radiation pressure in the outflow launching process.
Multiphase feedback processes in the Sy2 galaxy NGC 5643
(EDP Sciences, 2021-01-12) García Bernete, I.; Alonso Herrero, A.; García Burillo, S.; Pereira Santaella, M.; García Lorenzo, B.; Carrera, F. J.; Rigopoulou, D.; Ramos Almeida, C.; Villar Martín, M.; González Martín, O.; Hicks, E. K. S.; Labiano, Á.; Ricci, C.; Mateos, S.; Ministerio de Economía y Competitividad (MINECO); Science and Technology Facilities Council (STFC); Agencia Estatal de Investigación (AEI); Comunidad de Madrid; Fondo Nacional de Desarrollo Científico y Tecnológico (FONDECYT); García Lorenzo, M. B. [0000-0002-7228-7173]; Ramos Almeida, C. [0000-0001-8353-649X]; Carrero, F. J. [0000-0003-2135-9023]; Unidad de Excelencia Científica María de Maeztu Centro de Astrobiología del Instituto Nacional de Técnica Aeroespacial y CSIC, MDM-2017-0737; Unidad de Excelencia Científica María de Maeztu AGENCIA ESTATAL CONSEJO SUPERIOR DE INVESTIGACIONES CIENTIFICAS (CSIC), MDM-2017-0765
We study the multiphase feedback processes in the central ∼3 kpc of the barred Seyfert 2 galaxy NGC 5643. We used observations of the cold molecular gas (ALMA CO(2−1) transition) and ionized gas (MUSE IFU optical emission lines). We studied different regions along the outflow zone, which extends out to ∼2.3 kpc in the same direction (east-west) as the radio jet, as well as nuclear and circumnuclear regions in the host galaxy disk. The CO(2−1) line profiles of regions in the outflow and spiral arms show two or more different velocity components: one associated with the host galaxy rotation, and the others with out- or inflowing material. In the outflow region, the [O III]λ5007 Å emission lines have two or more components: the narrow component traces rotation of the gas in the disk, and the others are related to the ionized outflow. The deprojected outflowing velocities of the cold molecular gas (median Vcentral ∼ 189 km s−1) are generally lower than those of the outflowing ionized gas, which reach deprojected velocities of up to 750 km s−1 close to the active galactic nucleus (AGN), and their spatial profiles follow those of the ionized phase. This suggests that the outflowing molecular gas in the galaxy disk is being entrained by the AGN wind. We derive molecular and ionized outflow masses of ∼5.2 × 107 M⊙ (αCOGalactic) and 8.5 × 104 M⊙ and molecular and ionized outflow mass rates of ∼51 M⊙ yr−1 (αCOGalactic) and 0.14 M⊙ yr−1, respectively. This means that the molecular phase dominates the outflow mass and outflow mass rate, while the kinetic power and momentum of the outflow are similar in both phases. However, the wind momentum loads (Ṗout/ṖAGN) for the molecular and ionized outflow phases are ∼27−5 (αCOGalactic and αCOULIRGs) and < 1, which suggests that the molecular phase is not momentum conserving, but the ionized phase most certainly is. The molecular gas content (Meast ∼ 1.5 × 107 M⊙; αCOGalactic) of the eastern spiral arm is approximately 50−70% of the content of the western one. We interpret this as destruction or clearing of the molecular gas produced by the AGN wind impacting in the eastern side of the host galaxy (negative feedback process). The increase in molecular phase momentum implies that part of the kinetic energy from the AGN wind is transmitted to the molecular outflow. This suggests that in Seyfert-like AGN such as NGC 5643, the radiative or quasar and the kinetic or radio AGN feedback modes coexist and may shape the host galaxies even at kiloparsec scales through both positive and (mild) negative feedback.
A proto-pseudobulge in ESO 320-G030 fed by a massive molecular inflow driven by a nuclear bar
(EDP Sciences, 2021-01-07) González Alfonso, E.; Pereira Santaella, M.; Fischer, J.; García Burillo, S.; Yang, C.; Alonso Herrero, A.; Colina, L.; Ashby, M. L. N.; Smith, H. A.; Rico Villas, F.; Martín Pintado, J.; Cazzoli, S.; Stewart, F. P.; National Aeronautics and Space Administration (NASA); Agencia Estatal de Investigación (AEI); Comunidad de Madrid; European Commission (EC); 0000-0001-5285-8517; 0000-0001-6697-7808; 0000-0003-0444-6897; 0000-0002-8117-9991; 0000-0001-6794-2519; 0000-0002-7705-2525; 0000-0001-8266-8298; Centros de Excelencia Severo Ochoa, INSTITUTO DE ASTROFISICA DE ANDALUCIA (IAA), SEV-2017-0709; Unidad de Excelencia Científica María de Maeztu Centro de Astrobiología del Instituto Nacional de Técnica Aeroespacial y CSIC, MDM-2017-0737
Galaxies with nuclear bars are believed to efficiently drive gas inward, generating a nuclear starburst and possibly an active galactic nucleus. 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 or emission in 18 lines of H2O, which we combine with the ALMA H2O 4(23)-3(30) 448 GHz line (E-upper similar to 400 K) and continuum images to study the physical properties of the nuclear region. Radiative transfer models indicate that three nuclear components are required to account for the multi-transition H2O and continuum data. An envelope, with radius R similar to 130-150 pc, dust temperature T-dust approximate to 50 K, and N-H2 similar to 2x10(23) cm(-2), surrounds a nuclear disk with R similar to 40 pc that is optically thick in the far-infrared (tau (100 mu m)similar to 1.5-3, N-H2 similar to 2x10(24) cm(-2)). In addition, an extremely compact (R similar to 12 pc), warm (approximate to 100 K), and buried (tau (100 mu m)> 5, N-H2 greater than or similar to 5x10(24) cm(-2)) core component is required to account for the very high-lying H2O absorption lines. The three nuclear components account for 70% of the galaxy luminosity (SFR similar to 16-18 M-circle dot 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 M-inf similar to 20 M yr(-1), which is similar to the nuclear star formation rate (SFR), indicating that the starburst is sustained by the inflow. At lower r, similar to 100-150 pc, the inflow is best probed by the far-infrared OH ground-state doublets, with an estimated M-inf similar to 30 M yr(-1). The inferred short timescale of similar to 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 H218O, OH, 18OH, OH+, H2O+, H3O+, NH, NH2, NH3, CH, CH+, 13CH+, HF, SH, and C3, and we estimate their abundances.
Physics of ULIRGs with MUSE and ALMA: The PUMA project I. Properties of the survey and first MUSE data results
(EDP Sciences, 2021-02-16) Perna, M.; Arribas, S.; Pereira Santaella, M.; Colina, L.; Bellocchi, E.; Catalán Torrecilla, C.; Cazzoli, S.; Crespo Gómez, A.; Maiolino, R.; Piqueras López, J.; Rodríguez del Pino, B.; Comunidad de Madrid; Agencia Estatal de Investigación (AEI); European Research Council (ERC); Science and Technology Facilities Council (STFC); Perna, M. [0000-0002-0362-5941]; Arribas, S. [0000-0001-7997-1640]; Colina, L. [0000-0002-9090-4227]; Bellocchi, E. [0000-0001-9791-4228]; Catalán Torrecilla, C. [0000-0002-8067-0164]; Cazzoli, S. [0000-0002-7705-2525]; Maiolino, R. [0000-0002-4985-3819]; Piqueras López, J. [0000-0003-1580-1188]; Rodríguez del Pino, B. [0000-0001-5171-3930]; Centros de Excelencia Severo Ochoa, INSTITUTO DE ASTROFISICA DE ANDALUCIA (IAA), SEV-2017-0709
Context. Ultraluminous infrared galaxies (ULIRGs) are characterised by extreme starburst (SB) and active galactic nucleus (AGN) activity, and are therefore ideal laboratories for studying the outflow phenomena and their feedback effects. We have recently started a project called Physics of ULIRGs with MUSE and ALMA (PUMA), which is a survey of 25 nearby (z < 0.165) ULIRGs observed with the integral field spectrograph MUSE and the interferometer ALMA. This sample includes systems with both AGN and SB nuclear activity in the pre- and post-coalescence phases of major mergers.
Aims. The main goals of the project are (i) to study the prevalence of (ionised, neutral, and molecular) outflows as a function of the galaxy properties, (ii) to constrain the driving mechanisms of the outflows (e.g. distinguish between SB and AGN winds), and (iii) to identify and characterise feedback effects on the host galaxy. In this first paper, we present details on the sample selection, MUSE observations, and data reduction, and derive first high-level data products.
Methods. MUSE data cubes were analysed to study the dynamical status of each of the 21 ULIRGs observed so far, taking the stellar kinematics and the morphological properties inferred from MUSE narrow-band images into account. We also located the ULIRG nuclei, taking advantage of near-infrared (HST) and millimeter (ALMA) data, and studied their optical spectra to infer (i) the ionisation state through standard optical line ratio diagnostics, and (ii) outflows in both atomic ionised ([O III], Hα) and neutral (Na ID) gas.
Results. We show that the morphological and stellar kinematic classifications are consistent: post-coalescence systems are more likely associated with ordered motions, while interacting (binary) systems are dominated by non-ordered and streaming motions. We also find broad and asymmetric [O III] and Na ID profiles in almost all nuclear spectra, with line widths in the range [300 − 2000] km s−1, possibly associated with AGN- and SB-driven winds. This result reinforces previous findings that indicated that outflows are ubiquitous during the pre- and post-coalescence phases of major mergers.
Stellar kinematics in the nuclear regions of nearby LIRGs with VLT-SINFONI Comparison with gas phases and implications for dynamical mass estimations
(EDP Sciences, 2021-06-22) Crespo Gómez, A.; Piqueras López, J.; Arribas, S.; Pereira Santaella, M.; Colina, L.; Rodríguez del Pino, B.; National Aeronautics and Space Administration (NASA); Agencia Estatal de Investigación (AEI); Ministerio de Economía y Competitividad (MINECO); 0000-0003-2119-277X; 0000-0003-1580-1188; 0000-0001-7997-1640; 0000-0002-9090-4227; 0000-0001-5171-3930
Context. Nearby luminous infrared galaxies (LIRGs) are often considered to be the local counterpart of the star forming galaxy (SFG) population at z > 1. Therefore, local LIRGs are ideal systems with which to perform spatially resolved studies on the physical processes that govern these objects and to validate assumptions made in high-z studies because of a lack of sensitivity and/or spatial resolution.
Aims. In this work we analyse the spatially resolved kinematics of the stellar component in the inner r < 1–2 kpc of ten nearby (mean z = 0.014) LIRGs, establishing the dynamical state of the stars and estimating their dynamical masses (Mdyn). We compare the stellar kinematics with those for different gas phases, and analyse the relative effects of using different tracers when estimating dynamical masses.
Methods. We use seeing-limited SINFONI H- and K-band spectroscopy in combination with ancillary infrared (IR) imaging from various instruments (NICMOS/F160W, NACO/Ks and IRAC/3.6 μm). The stellar kinematics are extracted in both near-IR bands by fitting the continuum emission using pPXF. The velocity maps are then modelled as rotating discs and used to extract the geometrical parameters (i.e. centre, PA, and inclination), which are compared with their photometric counterparts extracted from the near-IR images. We use the stellar and the previously extracted gas velocity and velocity dispersion maps to estimate the dynamical mass using the different tracers.
Results. We find that the different gas phases have similar kinematics, whereas the stellar component is rotating with slightly lower velocities (i.e. V* ∼ 0.8Vg) but in significantly warmer orbits (i.e. σ* ∼ 2σg) than the gas phases, resulting in significantly lower V/σ for the stars (i.e. ∼1.5–2) than for the gas (i.e. ∼4–6). These ratios can be understood if the stars are rotating in thick discs while the gas phases are confined in dynamically cooler (i.e. thinner) rotating discs. However, these differences do not lead to significant discrepancies between the dynamical mass estimations based on the stellar and gas kinematics. This result suggests that the gas kinematics can be used to estimate Mdyn also in z ∼ 2 SFGs, a galaxy population that shares many structural and kinematic properties with local LIRGs.
