Source author record

L. Colina

L. Colina 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

55works
6topics
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

55 published item(s)

preprint2025arXiv

JWST/MIRI observations of the young TWA 27 system: hydrocarbon disk chemistry, silicate clouds, evidence for a CPD

The Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST) enables the characterisation of young self-luminous gas giants at previously inaccessible wavelengths, revealing physical processes in gas, dust, and clouds. We characterise the young planetary system TWA 27 (2M1207) in the mid-infrared (MIR), studying the atmosphere and disk spectra of the M9 brown dwarf TWA 27A and its L6 planetary-mass companion TWA 27b. We obtained data with the MIRI Medium Resolution Spectrometer (MRS) from 4.9 to 20 um, and MIRI Imaging in the F1000W and F1500W filters. We applied high-contrast imaging methods to extract the companion's spectral energy distribution up to 15 um at 0.78 arcsec separation and a contrast of 60. Combining these with published JWST/NIRSpec spectra, we analysed the 1-20 um range using self-consistent atmospheric grids and 0D slab models for molecular disk emission. The atmosphere of TWA 27A is well fitted by a BT-SETTL model with Teff 2780 K, log g 4.3, plus a 740 K blackbody for the inner disk rim. The disk shows at least 11 organic molecules, with no water or silicate dust emission detected. The atmosphere of TWA 27b is matched by a Teff 1400 K low-gravity model with extinction, best fit by the ExoREM grid. MIRI spectra and photometry for TWA 27b reveal a silicate cloud absorption feature between 8-10 um and significant (>5 sigma) infrared excess at 15 um consistent with circumplanetary disk emission. These MIRI observations provide new insights into TWA 27, revealing diverse features to study the formation and evolution of circumplanetary disks and young dusty atmospheres.

preprint2022arXiv

A dusty compact object bridging galaxies and quasars at cosmic dawn

Understanding how super-massive black holes form and grow in the early Universe has become a major challenge since the discovery of luminous quasars only 700 million years after the Big Bang. Simulations indicate an evolutionary sequence of dust-reddened quasars emerging from heavily dust-obscured starbursts that then transition to unobscured luminous quasars by expelling gas and dust. Although the last phase has been identified out to a redshift of 7.6, a transitioning quasar has not been found at similar redshifts owing to their faintness at optical and near-infrared wavelengths. Here we report observations of an ultraviolet compact object, GNz7q, associated with a dust-enshrouded starburst at a redshift of z=7.1899+/-0.0005. The host galaxy is more luminous in dust emission than any other known object at this epoch, forming 1,600 solar masses of stars per year within a central radius of 480 parsec. A red point source in the far-ultraviolet is identified in deep, high-resolution imaging and slitless spectroscopy. GNz7q is extremely faint in X-rays, which indicates the emergence of a uniquely ultraviolet compact star-forming region or a Compton-thick super-Eddington black-hole accretion disk at the dusty starburst core. In the latter case, the observed properties are consistent with predictions from cosmological simulations and suggest that GNz7q is an antecedent to unobscured luminous quasars at later epochs.

preprint2022arXiv

Compact molecular gas emission in local LIRGs among low- and high-z galaxies

We present new CO(2-1) observations of a representative sample of 24 local (z$<$0.02) luminous infrared galaxies (LIRGs) obtained at high spatial resolution ($<$100 pc) from ALMA. We derive the effective radii of the CO(2-1) and the 1.3 mm continuum emissions using the curve-of-growth method. LIRGs show an extremely compact cold molecular gas distribution (median R$_{CO}$ $\sim$0.7 kpc), which is a factor 2 smaller than the ionized gas, and 3.5 times smaller than the stellar size. The molecular size of LIRGs is similar to that of early-type galaxies (R$_{CO}\sim$1 kpc) and about a factor of 6 more compact than local Spirals of similar stellar mass. Only the CO emission in low-z ULIRGs is more compact than these local LIRGs by a factor of 2. Compared to high-z (1$<$z$<$6) systems, the stellar sizes and masses of local LIRGs are similar to those of high-z MS star-forming galaxies (SFG) and about a factor of 2-3 lower than sub-mm galaxies (SMG). The molecular sizes of high-z MS SFGs and SMGs are larger than those derived for LIRGs by a factor of $\sim$3 and $\sim$8, respectively. These results indicate that while low-z LIRGs and high-z MS-SFGs have similar stellar masses and sizes, the regions of current star formation (traced by the ionized gas) and of potential star-formation (traced by the molecular gas) are substantially smaller in LIRGs, and constrained to the central kpc region. High-z galaxies represent a wider population but their star-forming regions are more extended, even covering the overall size of the host galaxy. High-z galaxies have larger fractions of gas than low-z LIRGs, and therefore the formation of stars could be induced by interactions and mergers in extended disks or filaments with large enough molecular gas surface density involving physical mechanisms similar to those identified in the central kpc of LIRGs.

preprint2022arXiv

Spatially resolved star-formation relations of dense molecular gas in NGC 1068

We analyse the influence of the dynamical environment on the star formation (SF) relations of the dense molecular gas in the starburst (SB) ring of the Seyfert 2 galaxy NGC 1068. We used ALMA to image the emission of the 1-0 transitions of HCN and HCO+ with a resolution of 56 pc. We also used ancillary data of CO(1-0) at a resolution of ~100 pc, and CO(3-2) and its underlying continuum emission at ~40 pc. These observations allow us to probe a range of molecular gas densities (n(H2)~10$^{3-5}cm^{-3}$). The SF rate (SFR) is derived from Pa$α$ line emission imaged by HST/NICMOS. We analysed how SF relations change depending on the choice of aperture sizes and molecular gas tracer. The scatter in the Kennicutt-Schmidt relation is about a factor of two to three lower for the HCN and HCO+ lines compared to CO(1-0) for a common aperture. Correlations lose statistical significance below a critical spatial scale $\approx$300-400 pc. The SF efficiency of the dense molecular gas (SFEdense) shows a scattered distribution as a function of the HCN luminosity (L'(HCN)) around a mean value of $\simeq0.01$Myr$^{-1}$. An alternative prescription for SF relations, linking the SFEdense and the boundedness of the gas measured by the parameter b$\equivΣ$dense/$σ^2$, where $Σ$dense is the dense molecular gas surface density and $σ$ the velocity dispersion, resolves the degeneracy associated with the SFEdense-L'(HCN) plot. We identify two branches in the SFEdense-b plot that correspond to two dynamical environments in the SB ring, which are defined by their proximity to the bar-ring interface region. This region corresponds to the crossing of two density wave resonances, where an increased rate of cloud-cloud collisions would favour an enhanced compression of molecular gas. Our results suggest that galactic dynamics plays a major role in the efficiency of the gas conversion into stars.

preprint2022arXiv

The PUMA project. III. Incidence and properties of ionised gas disks in ULIRGs, associated velocity dispersion and its dependence on starburstiness

A classical scenario suggests that ULIRGs transform colliding spiral galaxies into a spheroid dominated early-type galaxy. Recent high-resolution simulations have instead shown that, under some circumstances, rotation disks can be preserved during the merging process or rapidly regrown after coalescence. Our goal is to analyze in detail the ionised gas kinematics in a sample of ULIRGs to infer the incidence of gas rotational dynamics in late-stage interacting galaxies and merger remnants. We analysed MUSE data of a sample of 20 nearby (z<0.165) ULIRGs, as part of the "Physics of ULIRGs with MUSE and ALMA" (PUMA) project. We found that 27% individual nuclei are associated with kpc-scale disk-like gas motions. The rest of the sample displays a plethora of gas kinematics, dominated by winds and merger-induced flows, which make the detection of rotation signatures difficult. On the other hand, the incidence of stellar disk-like motions is ~2 times larger than gaseous disks, as the former are probably less affected by winds and streams. The eight galaxies with a gaseous disk present relatively high intrinsic gas velocity dispersion (sigma = 30-85 km/s), rotationally-supported motions (with gas rotation velocity over velocity dispersion vrot/sigma > 1-8), and dynamical masses in the range (2-7)x1e10 Msun. By combining our results with those of local and high-z disk galaxies from the literature, we found a significant correlation between sigma and the offset from the main sequence (MS), after correcting for their evolutionary trends. Our results confirm the presence of kpc-scale rotating disks in interacting galaxies and merger remnants, with an incidence going from 27% (gas) to ~50% (stars). The ULIRGs gas velocity dispersion is up to a factor of ~4 higher than in local normal MS galaxies, similar to high-z starbursts as presented in the literature.

preprint2021arXiv

Duality in spatially resolved star-formation relations in local LIRGs

We analyse the star formation (SF) relations in a sample of 16 nearby luminous infrared galaxies (LIRGs) with more than 2800 regions defined on scales of 90 to 500 pc. We used ALMA to map the distribution of the cold molecular gas traced by the J = 2-1 line of CO and archival PaαHST/NICMOS imaging to trace the recent SF. In four objects, we find two different branches in the Kennicutt-Schmidt relation at 90 pc scales, suggesting the existence of a duality in this relation. The two branches correspond to two different dynamical environments within each galaxy. One branch, which corresponds to the central region of these galaxies 90% of the regions are located at radii < 0.85 kpc), shows higher gas and star formation rate surface densities with higher velocity dispersion. The other branch, which shows lower molecular gas and SF rate surface densities, corresponds to the more external disk regions (r $\sim$ 1 kpc). Despite the scatter, the SF efficiency of the galaxies with a dual behaviour increases with increasing boundedness as measured by the b parameter (b $\equiv$ $Σ_{H_{2}}$/$σ^{2}$ $\propto$ $α_{vir}^{-1}$). At larger spatial scales (250 and 500 pc), the duality disappears. The rest of the sample does not show evidence of this dual behaviour at any scale.

preprint2020arXiv

Physics of ULIRGs with MUSE and ALMA: The PUMA project I. Properties of the survey and first MUSE data results

Ultraluminous infrared galaxies (ULIRGs) are characterised by extreme starburst (SB) and AGN activity, and are therefore ideal laboratories for studying the outflow phenomena. 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. The main goals of the project are to study the prevalence of multi-phase outflows as a function of the galaxy properties, to constrain the driving mechanisms of the outflows (e.g. distinguish between SB and AGN winds), and to identify feedback effects on the host galaxy. In this first paper, we present details on the sample selection, MUSE observations, and derive first data products. MUSE data 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, using near-IR (HST) and mm (ALMA) data, and studied their optical spectra to infer the ionisation state through BPT diagnostics, and outflows in both ionised and neutral gas. 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 [OIII] and NaID profiles in almost all nuclear spectra, with line widths in the range 300-2000 km/s, 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.

preprint2020arXiv

Simulating JWST deep extragalactic imaging surveys and physical parameter recovery

We present a new prospective analysis of deep multi-band imaging with the James Webb Space Telescope (JWST). In this work, we investigate the recovery of high-redshift $5<z<12$ galaxies through extensive image simulations of accepted JWST programs such as CEERS in the EGS field and HUDF GTO. We introduce complete samples of $\sim300,000$ galaxies with stellar masses $\log(M_*/M_\odot)>6$ and redshifts $0<z<15$, as well as galactic stars, into realistic mock NIRCam, MIRI and HST images to properly describe the impact of source blending. We extract the photometry of the detected sources as in real images and estimate the physical properties of galaxies through spectral energy distribution fitting. We find that the photometric redshifts are primarily limited by the availability of blue-band and near-infrared medium-band imaging. The stellar masses and star-formation rates are recovered within $0.25$ and $0.3$ dex respectively, for galaxies with accurate photometric redshifts. Brown dwarfs contaminating the $z>5$ galaxy samples can be reduced to $<0.01$ arcmin$^{-2}$ with a limited impact on galaxy completeness. We investigate multiple high-redshift galaxy selection techniques and find the best compromise between completeness and purity at $5<z<10$ using the full redshift posterior probability distributions. In the EGS field, the galaxy completeness remains higher than $50\%$ for $m_\text{UV}<27.5$ sources at all redshifts, and the purity is maintained above $80$ and $60\%$ at $z\leq7$ and $10$ respectively. The faint-end slope of the galaxy UV luminosity function is recovered with a precision of $0.1-0.25$, and the cosmic star-formation rate density within $0.1$ dex. We argue in favor of additional observing programs covering larger areas to better constrain the bright end.

preprint2019arXiv

Rest-frame UV Properties of Luminous Strong Gravitationally Lensed Ly$α$ Emitters from the BELLS GALLERY Survey

We present deep rest-frame UV spectroscopic observations using the Gran Telescopio Canarias of six gravitationally lensed Lya emitters (LAEs) at $2.36<z<2.82$ selected from the BELLS GALLERY survey. By taking the magnifications into account, we show that LAEs can be as luminous as L(Lya) = 30x10$^{42}$ erg s-1 and M(UV) = -23 (AB) without invoking an AGN component, in contrast with previous findings. We measure Lya rest-frame equivalent widths, EW(Lya), ranging from 16Åto 50Åand Lya escape fractions, fesc(Lya), from 10% to 40%. Large EW(Lya) and fesc(Lya) are found predominantly in LAEs showing weak low-ionization ISM absorption (EW < 1Å) and narrow Lya profiles (< 300 km s-1 FWHM) with their peak close (< 80 km s-1) to their systemic redshifts, suggestive of less scatter from low HI column densities that favours the escape of Lya photons. We infer stellar metallicities of Z/Zsun ~ 0.2 in almost all LAEs by comparing the P-Cygni profiles of the wind lines NV1240Åand CIV1549Åwith those from stellar synthesis models. We also find a trend between M(UV) and the velocity offset of ISM absorption lines, such as the most luminous LAEs experience stronger outflows. The most luminous LAEs show star formation rates up to 180 Msun yr-1, yet they appear relatively blue ($β$(UV) ~ -1.8 to -2.0) showing evidence of little dust attenuation (E(B-V) = 0.10-0.14). These luminous LAEs may be particular cases of young starburst galaxies that have had no time to form large amounts of dust. If so, they are ideal laboratories to study the early phase of massive star formation, stellar and dust mass growth, and chemical enrichment histories of starburst galaxies at high-z.

preprint2016arXiv

A SINFONI view of the nuclear activity and circum-nuclear star formation in NGC 4303

We present new maps of emission-line flux distributions and kinematics in both ionized (traced by HI and [FeII] lines) and molecular (H2) gas of the inner 0.7x0.7kpc2 of the galaxy NGC4303, with a spatial resolution 40-80pc and velocity resolution 90-150 km/s obtained from near-IR integral field specroscopy using the VLT instrument SINFONI. The most promiment feature is a 200-250pc ring of circum-nuclear star-forming regions. The emission from ionized and molecular gas shows distinct flux distributions: while the strongest HI and [FeII] emission comes from regions in the west side of the ring (ages~4Myr), the H2 emission is strongest at the nucleus and in the east side of the ring (ages>10Myr). We find that regions of enhanced hot H2 emission are anti-correlated with those of enhanced [FeII] and HI emission, which can be attributed to post starburst regions that do not have ionizing photons anymore but still are hot enough (~2000K) to excite the H2 molecule. The line ratios are consistent with the presence of an AGN at the nucleus. The youngest regions have stellar masses in the range 0.3-1.5E5 MSun and ionized and hot molecular gas masses of ~0.25-1.2E4 Msun and 2.5-5 Msun, respectively. The stellar and gas velocity fields show a rotation pattern, with the gas presenting larger velocity amplitudes than the stars, with a deviation observed for the H2 along the nuclear bar, where increased velocity dispersion is also observed, possibly associated with non circular motions along the bar. The stars in the ring show smaller velocity dispersion than the surroundings, that can be attributed to a cooler dynamics due to their recent formation from cool gas.

preprint2016arXiv

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

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

preprint2016arXiv

High-velocity extended molecular outflow in the star-formation dominated luminous infrared galaxy ESO 320-G030

We analyze new high spatial resolution (~60 pc) ALMA CO(2-1) observations of the isolated luminous infrared galaxy ESO 320-G030 (d=48 Mpc) in combination with ancillary HST optical and near-IR imaging as well as VLT/SINFONI near-IR integral field spectroscopy. We detect a high-velocity (~450 km/s) spatially resolved (size~2.5 kpc; dynamical time ~3 Myr) massive (~10^7 Msun; mass rate~2-8 Msun/yr) molecular outflow originated in the central ~250 pc. We observe a clumpy structure in the outflowing cold molecular gas with clump sizes between 60 and 150 pc and masses between 10^5.5 and 10^6.4 Msun. The mass of the clumps decreases with increasing distance, while the velocity is approximately constant. Therefore, both the momentum and kinetic energy of the clumps decrease outwards. In the innermost (~100 pc) part of the outflow, we measure a hot-to-cold molecular gas ratio of 7x10^-5, which is similar to that measured in other resolved molecular outflows. We do not find evidence of an ionized phase in this outflow. The nuclear IR and radio properties are compatible with strong and highly obscured star-formation (A_k ~ 4.6 mag; SFR~15 Msun/yr). We do not find any evidence for the presence of an active galactic nucleus. We estimate that supernova explosions in the nuclear starburst (ν(SN) ~ 0.2 yr^-1) can power the observed molecular outflow. The kinetic energy and radial momentum of the cold molecular phase of the outflow correspond to about 2% and 20%, respectively, of the supernovae output. The cold molecular outflow velocity is lower than the escape velocity, so the gas will likely return to the galaxy disk. The mass loading factor is ~0.1-0.5, so the negative feedback due to this star-formation powered molecular outflow is probably limited.

preprint2016arXiv

Neutral gas outflows in nearby [U]LIRGs via optical NaD feature

We studied the properties of the neutral gas in a sample of 38 local [U]LIRGs, which mainly covers the less explored LIRG luminosity range. This study is based on the analysis of the spatially integrated and spatially resolved spectra of the NaD feature obtained with the VIMOS/VLT-IFS. Analyzing spatially integrated spectra, we find that the contribution of stars to the observed NaD-EWs is generally small (<35%), and therefore this feature is dominated by ISM absorption. After subtracting the stellar contribution, we find that the pure-ISM integrated spectra generally show blueshifted NaD profiles, indicating neutral gas outflow velocities V in the range 65-260 km/s. Excluding galaxies with AGNs, V shows a dependency with the SFR of the type V \propto SFR^0.15 (in rather good agreement with previous results). The spatially resolved analysis was performed for 40 galaxies, 22 of which have neutral gas velocity fields dominated by noncircular motions with signatures of cone-like winds. However, a modest number of targets (11/40) show disk rotation. We found that the wind masses are in the range (0.4-7.5)x10^8 Msun, reaching up to 3% of the dynamical mass of the host. The mass rates are typically only 0.2-0.4 times the SFR indicating that, the mass loss is too small to slow down the SF significantly. In the majority of cases, the velocity of the outflowing gas is not sufficient to escape the host potential well and, therefore, most of the gas rains back into the disk. On average V/v(esc) is higher in less massive galaxies, confirming that the galaxy mass has a primary role in shaping the recycling of gas and metals. The comparison between the wind power and kinetic power of the starburst associated with SNe indicates that the starburst could drive the outflows in nearly all the [U]LIRGs galaxies, as the wind power is generally lower than 20% of the kinetic power supplied by the starburst.

preprint2016arXiv

Sub-kpc star-formation law in the local luminous infrared galaxy IC 4687 as seen by ALMA

We analyze the spatially resolved (250 pc scales) and integrated star-formation (SF) law in the local luminous infrared galaxy (LIRG) IC4687. This is one of the first studies of the SF law on a starburst LIRG at these small spatial scales. We combined new interferometric ALMA CO(2-1) data with existing HST/NICMOS Pa$α$ narrow-band imaging and VLT/SINFONI near-IR integral field spectroscopy to obtain accurate extinction corrected SF rate (SFR) and cold molecular gas surface densities ($Σ_{gas}$ and $Σ_{SFR}$). We find that IC4687 forms stars very efficiently with an average depletion time ($t_{dep}$) of 160 Myr for the individual 250 pc regions. This is approximately one order of magnitude shorter than the $t_{dep}$ of local normal spirals and also shorter than that of main-sequence high-z objects, even when we use a Galactic $α_{CO}$ conversion factor. This result suggests a bimodal SF law in the $Σ_{SFR} \propto Σ_{gas}^{N}$ representation. A universal SF law is recovered if we normalize the $Σ_{gas}$ by the global dynamical time. However, at the spatial scales studied here, we find that the SF efficiency (or $t_{dep}$) does not depend on the local dynamical time for this object. Therefore, an alternative normalization (e.g., free-fall time) should be found if a universal SF law exists at these scales.

preprint2016arXiv

The complex evolutionary paths of local infrared bright galaxies: a high angular resolution mid-infrared view

We investigate the evolutionary connection between local IR-bright galaxies ($\log L_{\rm IR}\ge 11.4\,L_\odot$) and quasars. We use high angular resolution ($\sim$ 0.3-0.4 arcsec $\sim$ few hundred parsecs) $8-13\,μ$m ground-based spectroscopy to disentangle the AGN mid-IR properties from those of star formation. The comparison between the nuclear $11.3\,μ$m PAH feature emission and that measured with Spitzer/IRS indicates that the star formation is extended over a few kpc in the IR-bright galaxies. The AGN contribution to the total IR luminosity of IR-bright galaxies is lower than in quasars. Although the dust distribution is predicted to change as IR-bright galaxies evolve to IR-bright quasars and then to optical quasars, we show that the AGN mid-IR emission of all the quasars in our sample is not significantly different. In contrast, the nuclear emission of IR-bright galaxies with low AGN contributions appears more heavily embedded in dust although there is no clear trend with the interaction stage or projected nuclear separation. This suggests that the changes in the distribution of the nuclear obscuring material may be taking place rapidly and at different interaction stages washing out the evidence of an evolutionary path. When compared to normal AGN, the nuclear star formation activity of quasars appears to be dimming whereas it is enhanced in some IR-bright nuclei, suggesting that the latter are in an earlier star-formation dominated phase.

preprint2016arXiv

The complex, dusty narrow-line region of NGC 4388: Gas-jet interactions, outflows, and extinction revealed by near-IR spectroscopy

We present Gemini/GNIRS spectroscopy of the Seyfert 2 galaxy NGC 4388, with simultaneous coverage from 0.85 - 2.5 $μ$m. Several spatially-extended emission lines are detected for the first time, both in the obscured and unobscured portion of the optical narrow line region (NLR), allowing us to assess the combined effects of the central continuum source, outflowing gas and shocks generated by the radio jet on the central 280 pc gas. The HI and [FeII] lines allow us to map the extinction affecting the NLR. We found that the nuclear region is heavily obscured, with E(B-V) ~1.9 mag. To the NE of the nucleus and up to ~150 pc, the extinction remains large, ~1 mag or larger, consistent with the system of dust lanes seen in optical imaging. We derived position-velocity diagrams for the most prominent lines as well as for the stellar component. Only the molecular gas and the stellar component display a well-organized pattern consistent with disk rotation. Other emission lines are kinematically perturbed or show little evidence of rotation. Extended high-ionization emission of sulfur, silicon and calcium is observed to distances of at least 200 pc both NE and SW of the nucleus. We compared flux ratios between these lines with photoionization models and conclude that radiation from the central source alone cannot explain the observed high-ionization spectrum. Shocks between the radio-jet and the ambient gas are very likely an additional source of excitation. We conclude that NGC 4388 is a prime laboratory to study the interplay between all these mechanisms.

preprint2016arXiv

The impact of JWST broad-band filter choice on photometric redshift estimation

The determination of galaxy redshifts in James Webb Space Telescope (JWST)'s blank-field surveys will mostly rely on photometric estimates, based on the data provided by JWST's Near-Infrared Camera (NIRCam) at 0.6-5.0 μm and Mid Infrared Instrument (MIRI) at λ>5.0 μm. In this work we analyse the impact of choosing different combinations of NIRCam and MIRI broad-band filters (F070W to F770W), as well as having ancillary data at λ<0.6 μm, on the derived photometric redshifts (zphot) of a total of 5921 real and simulated galaxies, with known input redshifts z=0-10. We found that observations at λ<0.6 μm are necessary to control the contamination of high-z samples by low-z interlopers. Adding MIRI (F560W and F770W) photometry to the NIRCam data mitigates the absence of ancillary observations at λ<0.6 μm and improves the redshift estimation. At z=7-10, accurate zphot can be obtained with the NIRCam broad bands alone when S/N>=10, but the zphot quality significantly degrades at S/N<=5. Adding MIRI photometry with one magnitude brighter depth than the NIRCam depth allows for a redshift recovery of 83-99%, depending on SED type, and its effect is particularly noteworthy for galaxies with nebular emission. The vast majority of NIRCam galaxies with [F150W]=29 AB mag at z=7-10 will be detected with MIRI at [F560W, F770W]<28 mag if these sources are at least mildly evolved or have spectra with emission lines boosting the mid-infrared fluxes.

preprint2015arXiv

A CO-rich merger shaping a powerful and hyper-luminous infrared radio galaxy at z=2: the Dragonfly Galaxy

In the low-redshift Universe, the most powerful radio sources are often associated with gas-rich galaxy mergers or interactions. We here present evidence for an advanced, gas-rich (`wet') merger associated with a powerful radio galaxy at a redshift of z~2. This radio galaxy, MRC 0152-209, is the most infrared-luminous high-redshift radio galaxy known in the southern hemisphere. Using the Australia Telescope Compact Array, we obtained high-resolution CO(1-0) data of cold molecular gas, which we complement with HST/WFPC2 imaging and WHT long-slit spectroscopy. We find that, while roughly M(H2) ~ 2 x 10$^{10}$ M$_{\odot}$ of molecular gas coincides with the central host galaxy, another M(H2) ~ 3 x 10$^{10}$ M$_{\odot}$ is spread across a total extent of ~60 kpc. Most of this widespread CO(1-0) appears to follow prominent tidal features visible in the rest-frame near-UV HST/WFPC2 imaging. Ly$α$ emission shows an excess over HeII, but a deficiency over L(IR), which is likely the result of photo-ionisation by enhanced but very obscured star formation that was triggered by the merger. In terms of feedback, the radio source is aligned with widespread CO(1-0) emission, which suggests that there is a physical link between the propagating radio jets and the presence of cold molecular gas on scales of the galaxy's halo. Its optical appearance, combined with the transformational stage at which we witness the evolution of MRC 0152-209, leads us to adopt the name `Dragonfly Galaxy'.

preprint2015arXiv

A mid-infrared spectroscopic atlas of local active galactic nuclei on sub-arcsecond resolution using GTC/CanariCam

We present an atlas of mid-infrared (mid-IR) ~7.5-13micron spectra of 45 local active galactic nuclei (AGN) obtained with CanariCam on the 10.4m Gran Telescopio CANARIAS (GTC) as part of an ESO/GTC large program. The sample includes Seyferts and other low luminosity AGN (LLAGN) at a median distance of 35Mpc and luminous AGN, namely PG quasars, (U)LIRGs, and radio galaxies (RG) at a median distance of 254Mpc. To date, this is the largest mid-IR spectroscopic catalog of local AGN at sub-arcsecond resolution (median 0.3arcsec). The goal of this work is to give an overview of the spectroscopic properties of the sample. The nuclear 12micron luminosities of the AGN span more than four orders of magnitude, nu*Lnu(12micron)~ 3e41-1e46erg/s. In a simple mid-IR spectral index vs. strength of the 9.7micron silicate feature diagram most LLAGN, Seyfert nuclei, PG quasars, and RGs lie in the region occupied by clumpy torus model tracks. However, the mid-IR spectra of some might include contributions from other mechanisms. Most (U)LIRG nuclei in our sample have deeper silicate features and flatter spectral indices than predicted by these models suggesting deeply embedded dust heating sources and/or contribution from star formation. The 11.3micron PAH feature is clearly detected in approximately half of the Seyfert nuclei, LLAGN, and (U)LIRGs. While the RG, PG quasars, and (U)LIRGs in our sample have similar nuclear 12micron luminosities, we do not detect nuclear PAH emission in the RGs and PG quasars.

preprint2015arXiv

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

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

preprint2015arXiv

Sub-arcsec mid-IR observations of NGC 1614: Nuclear star-formation or an intrinsically X-ray weak AGN?

We present new mid-infrared N-band spectroscopy and Q-band photometry of the local luminous infrared galaxy NGC1614, one of the most extreme nearby starbursts. We analyze the mid-IR properties of the nucleus (central 150 pc) and four regions of the bright circumnuclear (diameter~600 pc) star-forming (SF) ring of this object. The nucleus differs from the circumnuclear SF ring by having a strong 8-12 micron continuum (low 11.3 micron PAH equivalent width). These characteristics, together with the nuclear X-ray and sub-mm properties, can be explained by an X-ray weak active galactic nucleus (AGN), or by peculiar SF with a short molecular gas depletion time and producing an enhanced radiation field density. In either case, the nuclear luminosity (L(IR) < 6e43 erg/s) is only <5% of the total bolometric luminosity of NGC1614. So this possible AGN does not dominate the energy output in this object. We also compare three star-formation rate (SFR) tracers (Pa$α$, 11.3 micron PAH, and 24 micron emissions) at 150 pc scales in the circumnuclear ring. In general, we find that the SFR is underestimated (overestimated) by a factor of 2-4 (2-3) using the 11.3 micron PAH (24 micron) emission with respect to the extinction corrected Pa$α$ SFR. The former can be explained because we do not include diffuse PAH emission in our measurements, while the latter might indicate that the dust temperature is particularly warmer in the central regions of NGC1614.

preprint2015arXiv

The Mid-Infrared Instrument for JWST, II: Design and Build

The Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST) provides measurements over the wavelength range 5 to 28.5 microns. MIRI has, within a single 'package', four key scientific functions: photometric imaging, coronagraphy, single-source low-spectral resolving power (R ~ 100) spectroscopy, and medium-resolving power (R ~ 1500 to 3500) integral field spectroscopy. An associated cooler system maintains MIRI at its operating temperature of < 6.7 K. This paper describes the driving principles behind the design of MIRI, the primary design parameters, and their realization in terms of the 'as-built' instrument. It also describes the test program that led to delivery of the tested and calibrated Flight Model to NASA in 2012, and the confirmation after delivery of the key interface requirements.

preprint2015arXiv

The Mid-Infrared Instrument for the James Webb Space Telescope I: Introduction

MIRI (the Mid-Infrared Instrument for the James Webb Space Telescope (JWST)) operates from 5 to 28.5 microns and combines over this range: 1.) unprecedented sensitivity levels; 2.) sub-arcsec angular resolution; 3.) freedom from atmospheric interference; 4.) the inherent stability of observing in space; and 5.) a suite of versatile capabilities including imaging, low and medium resolution spectroscopy (with an integral field unit), and coronagraphy. We illustrate the potential uses of this unique combination of capabilities with various science examples: 1.) imaging exoplanets; 2.) transit and eclipse spectroscopy of exoplanets; 3.) probing the first stages of star and planet formation, including identifying bioactive molecules; 4.) determining star formation rates and mass growth as galaxies are assembled; and 5.) characterizing the youngest massive galaxies. This paper is the introduction to a series of ten covering all aspects of the instrument.

preprint2015arXiv

The Spanish Square Kilometre Array White Book

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

Nuclear 11.3$μ$m PAH emission in local active galactic nuclei

We present Gran Telescopio CANARIAS CanariCam 8.7$μ$m imaging and 7.5-13$μ$m spectroscopy of six local systems known to host an active galactic nucleus (AGN) and have nuclear star formation. Our main goal is to investigate whether the molecules responsible for the 11.3$μ$m polyclyclic aromatic hydrocarbon (PAH) feature are destroyed in the close vicinity of an AGN. We detect 11.3$μ$m PAH feature emission in the nuclear regions of the galaxies as well as extended PAH emission over a few hundred parsecs. The equivalent width (EW) of the feature shows a minimum at the nucleus but increases with increasing radial distances, reaching typical star-forming values a few hundred parsecs away from the nucleus. The reduced nuclear EW are interpreted as due to increased dilution from the AGN continuum rather than destruction of the PAH molecules. We conclude that at least those molecules responsible for the 11.3$μ$m PAH feature survive in the nuclear environments as close as 10pc from the AGN and for Seyfert-like AGN luminosities. We propose that material in the dusty tori, nuclear gas disks, and/or host galaxies of AGN is likely to provide the column densities necessary to protect the PAH molecules from the AGN radiation field.

preprint2014arXiv

Spatially resolved kinematics, galactic wind, and quenching of star formation in the luminous infrared galaxy IRAS F11506-3851

We present a multi-wavelength integral field spectroscopic study of the low-z LIRG IRAS F11506-3851, on the basis of VIMOS and SINFONI (ESO-VLT) observations. The morphology and the 2D kinematics of the gaseous (neutral and ionized) and stellar components have been mapped using the NaD doublet, the H$α$ line, and the near-IR CO(2-0) and CO(3-1) bands. The kinematics of the ionized gas and the stars are dominated by rotation, with large observed velocity amplitudes and centrally peaked velocity dispersion maps. The stars lag behind the warm gas and represent a dynamically hotter system, as indicated by the observed dynamical ratios. Thanks to these IFS data we have disentangled the contribution of the stars and the ISM to the NaD feature, finding that it is dominated by the absorption of neutral gas clouds in the ISM. The neutral gas 2D kinematics shows a complex structure dominated by two components. On the one hand, the thick slowly rotating disk lags significantly compared to the ionized gas and the stars, with an irregular and off-center velocity dispersion map. On the other hand, a kpc-scale neutral gas outflow is observed along the semi-minor axis of the galaxy, as revealed by large blueshifted velocities (30-154 km/s). We derive an outflowing mass rate in neutral gas of about 48 $\dot{M_{\rm w}}$/yr. Although this implies a global mass loading factor of 1.4, the 2D distribution of the ongoing SF suggests a much larger value of mass loading factor associated with the inner regions (R$<$200 pc), where the current SF represents only 3 percent of the total. All together these results strongly suggest that we are witnessing (nuclear) quenching due to SF feedback in IRAS F11506-3851. However, the relatively large mass of molecular gas detected in the nuclear region via the H2 1-0 S(1) line suggests that further episodes of SF may take place again.

preprint2013arXiv

Uncovering the deeply embedded AGN activity in the nuclear regions of the interacting galaxy Arp299

We present mid-infrared (MIR) 8-13micron spectroscopy of the nuclear regions of the interacting galaxy Arp299 (IC694+NGC3690) obtained with CanariCam (CC) on the 10.4m Gran Telescopio Canarias (GTC). The high angular resolution (~0.3-0.6arcsec) of the data allows us to probe nuclear physical scales between 60 and 120pc, which is a factor of 10 improvement over previous MIR spectroscopic observations of this system. The GTC/CC spectroscopy displays evidence of deeply embedded Active Galactic Nucleus (AGN) activity in both nuclei. The GTC/CC nuclear spectrum of NGC3690/Arp299-B1 can be explained as emission from AGN-heated dust in a clumpy torus with both a high covering factor and high extinction along the line of sight. The estimated bolometric luminosity of the AGN in NGC3690 is 3.2(+/-0.6)x10^44 erg/s. The nuclear GTC/CC spectrum of IC694/Arp299-A shows 11.3micron polycyclic aromatic hydrocarbon (PAH) emission stemming from a deeply embedded (A_V~24mag) region of less than 120pc in size. There is also a continuum-emitting dust component. If associated with th putative AGN in IC694, we estimate that it would be approximately 5 times less luminous than the AGN in NGC3690. The presence of dual AGN activity makes Arp299 a good example to study such phenomenon in the early coalescence phase of interacting galaxies.

preprint2013arXiv

VLT-SINFONI integral field spectroscopy of low-z luminous and ultraluminous infrared galaxies II. 2D extinction structure and distance effects

We present a 2D study of the internal extinction on (sub)kiloparsec scales of a sample of local (z<0.1) (U)LIRGs, based on near-infrared Pa_a, Br_d, and Br_g line ratios, obtained with VLT-SINFONI integral-field spectroscopy. The 2D extinction (Av) distributions of the objects, map regions of ~3x3 kpc (LIRGs) and ~12x12 kpc (ULIRGs), with average angular resolutions (FWHM) of ~0.2 kpc and ~0.9 kpc, respectively. The individual Av galaxy distributions indicate a very clumpy dust structure already on sub-kiloparsec scales, with values (per spaxel) ranging from Av ~ 1 to 20 mag in LIRGs, and from Av ~ 2 to 15 mag in ULIRGs. As a class, the median values of the distributions are Av=5.3 mag and Av=6.5 mag for the LIRG and ULIRG subsamples, respectively. We evaluated the effects of the galaxy distance in the measurements of the extinction as a function of the linear scale (in kpc) of the spaxel (i.e. due to the limited angular resolution of the observations). If the distribution of the gas/dust and star-forming regions in local LIRGs (63 Mpc, 40 pc/spaxel on average) is the same for galaxies at greater distances, the observed median Av values based on emission line ratios would be a factor ~0.8 lower at the average distance of our ULIRG sample (328 Mpc, 0.2 kpc/spaxel), and a factor ~0.67 for galaxies located at distances of more than 800 Mpc (0.4 kpc/spaxel). This distance effect would have implications for deriving the intrinsic extinction in high-z star-forming galaxies and for subsequent properties such as star formation rate, star formation surface density, and KS- law, based on H_alpha line fluxes. If local LIRGs are analogues of the main-sequence star-forming galaxies at cosmological distances, the extinction values (Av) derived from the observed emission lines in these high-z sources would need to be increased by a factor 1.4 on average.

preprint2012arXiv

Gas fraction and star formation efficiency at z < 1.0

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

preprint2012arXiv

Integral Field Spectroscopy based Hα sizes of local Luminous and Ultraluminous Infrared Galaxies. A Direct Comparison with high-z Massive Star Forming Galaxies

Aims. We study the analogy between local U/LIRGs and high-z massive SFGs by comparing basic Hα structural characteristics, such as size, and luminosity (and SFR) surface density, in an homogeneous way (i.e. same tracer and size definition, similar physical scales). Methods. We use Integral Field Spectroscopy based Hα emission maps for a representative sample of 54 local U/LIRGs (66 galaxies). From this initial sample we select 26 objects with Hα luminosities (L(Hα)) similar to those of massive (i.e. M\ast \sim 10^10 M\odot or larger) SFGs at z \sim 2, and observed on similar physical scales. Results. The sizes of the Hα emitting region in the sample of local U/LIRGs span a large range, with r1/2(Hα) from 0.2 to 7 kpc. However, about 2/3 of local U/LIRGs with Lir > 10^11.4 L\odot have compact Hα emission (i.e. r1/2 < 2 kpc). The comparison sample of local U/LIRGs also shows a higher fraction (59%) of objects with compact Hα emission than the high-z sample (25%). This gives further support to the idea that for this luminosity range the size of the star forming region is a distinctive factor between local and distant galaxies of similar SF rates. However, when using Hα as a tracer for both local and high-z samples, the differences are smaller than the ones recently reported using a variety of other tracers. Despite of the higher fraction of galaxies with compact Hα emission, a sizable group (\sim 1/3) of local U/LIRGs are large (i.e. r1/2 > 2 kpc). These are systems showing pre-coalescence merger activity and they are indistinguishable from the massive high-z SFGs galaxies in terms of their Hα sizes, and luminosity and SFR surface densities.

preprint2012arXiv

Spatially Resolved [FeII] 1.64 μm Emission in NGC 5135. Clues for Understanding the Origin of the Hard X-rays in Luminous Infrared Galaxies

Spatially resolved near-IR and X-ray imaging of the central region of the Luminous Infrared Galaxy NGC 5135 is presented. The kinematical signatures of strong outflows are detected in the [FeII]1.64 μm emission line in a compact region at 0.9 kpc from the nucleus. The derived mechanical energy release is consistent with a supernova rate of 0.05-0.1 yr$^{-1}$. The apex of the outflowing gas spatially coincides with the strongest [FeII] emission peak and with the dominant component of the extranuclear hard X-ray emission. All these features provide evidence for a plausible direct physical link between supernova-driven outflows and the hard X-ray emitting gas in a LIRG. This result is consistent with model predictions of starbursts concentrated in small volumes and with high thermalization efficiencies. A single high-mass X-ray binary (HMXB) as the major source of the hard X-ray emission although not favoured, cannot be ruled out. Outside the AGN, the hard X-ray emission in NGC 5135 appears to be dominated by the hot ISM produced by supernova explosions in a compact star-forming region, and not by the emission due to HMXB. If this scenario is common to U/LIRGs, the hard X-rays would only trace the most compact (< 100 pc) regions with high supernova and star formation densities, therefore a lower limit to their integrated star formation. The SFR derived in NGC 5135 based on its hard X-ray luminosity is a factor of two and four lower than the values obtained from the 24 μm and soft X-ray luminosities, respectively.

preprint2012arXiv

Spatially resolved kinematics of the central regions of M83: hidden mass signatures and the role of supernovae

The barred grand-design spiral M83 (NGC 5236) is one of the most studied galaxies given its proximity, orientation, and particular complexity. Nonetheless, many aspects of the central regions remain controversial conveying our limited understanding of the inner gas and stellar kinematics, and ultimately of the nucleus evolution. In this work, we present AO VLT-SINFONI data of its central ~235x140 pc with an unprecedented spatial resolution of ~0.2 arcsec, corresponding to ~4 pc. We have focused our study on the distribution and kinematics of the stars and the ionised and molecular gas by studying in detail the Pa_alpha and Br_gamma emission, the H_2 1-0S(1) line at 2.122 micron and the [FeII] line at 1.644 micron, together with the CO absorption bands at 2.293 micron and 2.323 micron. Our results reveal a complex situation where the gas and stellar kinematics are totally unrelated. Supernova explosions play an important role in shaping the gas kinematics, dominated by shocks and inflows at scales of tens of parsecs that make them unsuitable to derive general dynamical properties. We propose that the location of the nucleus of M83 is unlikely to be related to the off-centre 'optical nucleus'. The study of the stellar kinematics reveals that the optical nucleus is a gravitationally bound massive star cluster with M_dyn = (1.1 \pm 0.4)x10^7 M_sun, formed by a past starburst. The kinematic and photometric analysis of the cluster yield that the stellar content of the cluster is well described by an intermediate age population of log T(yr) = 8.0\pm0.4, with a mass of M \simeq (7.8\pm2.4)x10^6 M_sun.

preprint2012arXiv

Star Formation Efficiency at Intermediate Redshift

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

preprint2012arXiv

Studying the kinematic asymmetries of disks and post-coalescence mergers using a new `kinemetry' criterion

We have obtained VIMOS/VLT optical integral field spectroscopy (IFS) data for a sample of 4 LIRGs which have been selected at a similar distance ($\sim$ 70 Mpc) to avoid relative resolution effects. They have been classified in two groups (isolated disk and post-coalescence mergers) according to their morphology. The $kinemetry$ method (developed by Krajnovic and coworkers) is used to characterize the kinematic properties of these galaxies and to discuss new criteria for distinguishing their status. We present and discuss new kinematic maps (i.e., velocity field and velocity dispersion) for these four galaxies. The morphological and kinematic classifications of these systems are consistent, with disks having lower kinematic asymmetries than post-coalescence mergers. We then propose and discuss a new kinematic criterion to differentiate these two groups. This criterion distinguishes better these two categories and has the advantage of being less sensitive to angular resolution effects. According to the previous criteria,the present post-coalescence systems would have been classified as disks. This indicates that the separation of disks from mergers is subjective to the definition of `merger'. It also suggests that previous estimates of the merger/disk ratio could have been underestimated, but larger samples are necessary to establish a firmer conclusion.

preprint2012arXiv

The Nuclear Infrared Emission of Low-Luminosity Active Galactic Nuclei

We present high-resolution mid-infrared (MIR) imaging, nuclear spectral energy distributions (SEDs) and archival Spitzer spectra for 22 low-luminosity active galactic nuclei (LLAGN; Lbol \lesssim 10^42 erg/sec). Infrared (IR) observations may advance our understanding of the accretion flows in LLAGN, the fate of the obscuring torus at low accretion rates, and, perhaps, the star formation histories of these objects. However, while comprehensively studied in higher-luminosity Seyferts and quasars, the nuclear IR properties of LLAGN have not yet been well-determined. We separate the present LLAGN sample into three categories depending on their Eddington ratio and radio emission, finding different IR characteristics for each class. (I) At the low-luminosity, low-Eddington ratio (log Lbol/LEdd < -4.6) end of the sample, we identify "host-dominated" galaxies with strong polycyclic aromatic hydrocarbon bands that may indicate active (circum-)nuclear star formation. (II) Some very radio-loud objects are also present at these low Eddington ratios. The IR emission in these nuclei is dominated by synchrotron radiation, and some are likely to be unobscured type 2 AGN that genuinely lack a broad line region. (III) At higher Eddington ratios, strong, compact nuclear sources are visible in the MIR images. The nuclear SEDs of these galaxies are diverse; some resemble typical Seyfert nuclei, while others lack a well-defined MIR "dust bump". Strong silicate emission is present in many of these objects. We speculate that this, together with high ratios of silicate strength to hydrogen column density, could suggest optically thin dust and low dust-to-gas ratios, in accordance with model predictions that LLAGN do not host a Seyfert-like obscuring torus.

preprint2012arXiv

VLT-SINFONI integral field spectroscopy of low-z luminous and ultraluminous infrared galaxies I. Atlas of the 2D gas structure

We present an atlas of a sample of local (U)LIRGs covering the luminosity range log(L_IR/L_sun)=11.1-12.4. The atlas is based on near-infrared H and K-band VLT-SINFONI IFS, and presents the ionised, partially ionised, and warm molecular gas 2D flux distributions and kinematics over a FoV of 3x3 kpc (LIRGs) and 12x12kpc (ULIRGs) and with average linear resolutions of 0.2kpc and 0.9kpc, respectively. The different phases of the gas show a wide morphological variety with the nucleus as the brightest Br_g source for 33% of the LIRGs and 71% of the ULIRGs, whereas all the (U)LIRGs have their maximum H_2 emission in their nuclear regions. In LIRGs, the ionised gas distribution is dominated by the emission from the star-forming rings or giant HII regions in the spiral arms. The Br_g and [FeII] line at 1.644 micron trace the same structures, although the emission peaks at different locations in some of the objects, and the [FeII] seems to be more extended and diffuse. The ULIRG subsample contains mainly pre-coalescence interacting systems. Although the peaks of the molecular gas emission and the continuum coincide in 71% of the ULIRGs, regions with intense Pa_a (Br_g) emission tracing luminous star-forming regions located at distances of 2-4kpc away from the nucleus are also detected, usually associated with secondary nuclei or tidal tails. The gas kinematics in LIRGs are primarily due to rotational motions around the centre of the galaxy, although local deviations associated with radial flows and/or regions of higher velocity dispersions are present. The ionised and molecular gas share the same kinematics to first order, showing slight differences in the velocity amplitudes in some cases, whereas the average velocity dispersions are compatible within uncertainties. As expected, the kinematics of the ULIRG subsample is more complex, owing to the interacting nature of the objects of the sample.

preprint2011arXiv

Discovery of Two Supernovae in the Nuclear Regions of the Luminous Infrared Galaxy IC 883

We report the discovery of two consecutive supernovae (SNe), 2010cu and 2011hi, located at 0.37" (180 pc) and 0.79" (380 pc) projected distance respectively from the centre of the K-band nucleus of the luminous infrared galaxy IC 883. The SNe were discovered in an ongoing near-infrared K-band search for core-collapse SNe in such galaxies using the ALTAIR/NIRI adaptive optics system with laser guide star at the Gemini-North Telescope. These are thus the closest SNe yet discovered to a LIRG nucleus in optical or near-infrared wavelengths. The near-infrared light curves and colours of both SNe are consistent with core-collapse events. Both SNe seem to suffer from relatively low host galaxy extinction suggesting that regardless of their low projected galactocentric distances, they are not deeply buried in the nuclear regions of the host galaxy.

preprint2011arXiv

Extranuclear Halpha-emitting complexes in low-z (U)LIRGs: Precursors of tidal dwarf galaxies?

(Abridged)This paper characterizes the physical and kinematic properties of external massive star-forming regions in a sample of (U)LIRGs. We use high angular resolution ACS images from the HST B and I bands, as well as Halpha-line emission maps obtained with IFS. We find 31 external Halpha-emitting (young star-forming) complexes in 11 (U)LIRGs. These complexes have in general similar sizes, luminosities, and metallicities to extragalactic giant HII regions and TDG candidates found in less luminous mergers and compact groups of galaxies. We assess the mass content and the likelihood of survival as TDGs of the 22 complexes with simple structures in the HST images based on their photometric, structural, and kinematic properties. The dynamical tracers used (radius-sigma and luminosity-sigma diagrams) indicate that most of the complexes might be self-gravitating entities. The resistance to forces from the parent galaxy is studied by considering the tidal mass of the candidate and its relative velocity with respect to the parent galaxy. After combining the results of previous studies of TDG searches in ULIRGs a total of 9 complexes satisfy most of the applied criteria and thus show a high-medium or high likelihood of survival, their total mass likely being compatible with that of dwarf galaxies. They are defined as TDG candidates. We propose that they probably formed more often during the early phases of the interaction. Combining all data for complexes with IFS data where a significant fraction of the system is covered, we infer a TDG production rate of 0.3 candidates with the highest probabilities of survival per system for the (U)LIRGs class. This rate, though, might decrease to 0.1 after the systems in (U)LIRGs have evolved for 10 Gyr, for long-lived TDGs, which would imply that no more than 5-10 % of the overall dwarf population could be of tidal origin.

preprint2011arXiv

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

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

preprint2011arXiv

Integrated spectra extraction based on signal-to-noise optimization using Integral Field Spectroscopy

We propose and explore the potential of a method to extract high signal-to-noise (S/N) integrated spectra related to physical and/or morphological regions on a 2-dimensional field using Integral Field Spectroscopy (IFS) observations by employing an optimization procedure based on either continuum (stellar) or line (nebular) emission features. The optimization method is applied to a set of IFS VLT-VIMOS observations of (U)LIRG galaxies, describing the advantages of the optimization by comparing the results with a fixed-aperture, single spectrum case, and by implementing some statistical tests. We demonstrate that the S/N of the IFS optimized integrated spectra is significantly enhanced when compared with the single aperture unprocessed case. We provide an iterative user-friendly and versatile IDL algorithm that allows the user to spatially integrate spectra following more standard procedures. This is made available to the community as part of the PINGSoft IFS software package.

preprint2011arXiv

Ionized outflows in SDSS type 2 quasars at z~0.3-0.6

We have analyzed the spatially integrated kinematic properties of the ionized gas within the inner r<few kpc in 13 optically selected SDSS type 2 quasars at z~0.3-0.6, using the [OIII]4959,5007 lines. The line profiles show a significant asymmetry in 11 objects. There is a clear preference for blue asymmetries, which are found in 9/13 quasars at 10% intensity level. In coherence with studies on other types of active and non active galaxies, we propose that the asymmetries are produced by outflows where differential dust extinction is at work. This scenario is favoured by other results we find: in addition to quiescent ambient gas, whose kinematic properties are consistent with gravitational motions, we have discovered highly perturbed gas in all objects. This gas emits very broad lines (R=FWHM[OIII]/FWHM(stars)>~2). While the quiescent gas shows small or null velocity shifts relative to the systemic velocity, the highly perturbed gas trends to show larger shifts which, moreover, are blueshifts in general. Within a given object, the most perturbed gas trends to have the largest blueshift as well. All together support that the perturbed gas, which is responsible for the blue asymmetries of the line profiles, is outflowing. Although some bias affects the sample, we argue that ionized gas outflows are a common phenomenon in optically selected type 2 quasars at 0.3< z < 0.6.

preprint2011arXiv

Star formation laws in Luminous Infrared Galaxies. New observational constraints on models

The observational study of star formation relations in galaxies is central to unraveling the physical processes at work on local and global scales. We wish to expand the sample of extreme starbursts, represented by local LIRGs and ULIRGs, with high quality observations in the 1-0 line of HCN. We study if a universal law can account for the star formation relations observed for the dense molecular gas in normal star forming galaxies and extreme starbursts. We have used the IRAM 30m telescope to observe a sample of 19 LIRGs in the 1-0 lines of CO, HCN and HCO+. The analysis of the new data proves that the efficiency of star formation in the dense molecular gas (SFE-dense) of extreme starbursts is a factor 3-4 higher compared to normal galaxies. We find a duality in Kennicutt-Schmidt (KS) laws that is reinforced if we account for the different conversion factor for HCN (alpha-HCN) in extreme starbursts and for the unobscured star formation rate in normal galaxies. This result extends to the higher molecular densities probed by HCN lines the more extreme bimodal behavior of star formation laws, derived from CO molecular lines by two recent surveys. We have confronted our observations with the predictions of theoretical models in which the efficiency of star formation is determined by the ratio of a constant star formation rate per free-fall time (SFR-ff) to the local free-fall time. We find that it is possible to fit the observed differences in the SFE-dense between normal galaxies and LIRGs/ULIRGs using a common constant SFR-ff and a set of physically acceptable HCN densities, but only if SFR-ff~0.005-0.01 and/or if alpha-HCN is a factor of a few lower than our favored values. Star formation recipes that explicitly depend on the galaxy global dynamical time scales do not significantly improve the fit to the new HCN data presented in this work.

preprint2011arXiv

The local LIRG NGC5135: AGN and SN activity traced by NIR IFU spectroscopy

By observing the local luminous infrared galaxy NGC5135 with the near-IR IFU spectrograph SINFONI (ESO-VLT), we report an approx. 600 pc (in projection) AGN outflow traced by [SiVI]1.96mum emission. This is the largest outflow traced by a coronal line ever reported. Its large spatial scale suggests that shocks, in addition to AGN continuum emission, are needed to locally produce [SiVI] emission. We also show, for the first time, clear kinematical evidence of the AGN-outflow vs. ISM interaction through variations in the 2D velocity fields of different gas phases. Such local perturbations in the kinematics clearly match the outflow structure. We use the [FeII]1.64mum emission, a supernovae tracer, to estimate the supernovae rate in different star-forming knots (approx. 250 pc across) within the central 2.3 kpc of NGC5135. The estimated supernovae rates go from 0.02 to 0.08 yr^{-1} being in excellent agreement with predictions from 6 cm radio emission in the same areas.

preprint2010arXiv

A High Spatial Resolution Mid-Infrared Spectroscopic Study of the Nuclei and Star-Forming Regions in Luminous Infrared Galaxies

We present a high spatial (diffraction-limited) resolution (~0.3") mid-infrared (MIR) spectroscopic study of the nuclei and star-forming regions of 4 local luminous infrared galaxies (LIRGs) using T-ReCS on the Gemini South telescope. We investigate the spatial variations of the features seen in the N-band spectra of LIRGs on scales of ~100 pc, which allow us to separate the AGN emission from that of the star formation (SF). We compare our Gemini T-ReCS nuclear and integrated spectra of LIRGs with those obtained with Spitzer IRS. The 9.7um silicate absorption feature is weaker in the nuclei of the LIRGs than in the surrounding regions. This is probably due to the either clumpy or compact environment of the central AGN or young, nuclear starburst. We find that the [NeII] luminosity surface density is tightly and directly correlated with that of Pa-alpha for the LIRG star-forming regions (slope of 1.00+-0.02). Although the 11.3um PAH feature shows also a trend with Pa-alpha, this is not common for all the regions. We also find that the [NeII]\Pa-alpha ratio does not depend on the Pa-alpha equivalent width (EW), i.e., on the age of the ionizing stellar populations, suggesting that, on the scales probed here, the [NeII] emission line is a good tracer of the SF activity in LIRGs. On the other hand, the 11.3um PAH\Pa-alpha ratio increases for smaller values of the Pa-alpha EW (increasing ages), indicating that the 11.3um PAH feature can also be excited by older stars than those responsible for the Pa-alpha emission. Additional high spatial resolution observations are essential to investigate, in a statistical way, the star formation in local LIRGs at the smallest scales and to probe ultimately whether they share the same physical properties as high-z LIRGs, ULIRGs and submillimiter galaxies.

preprint2010arXiv

A Spitzer/IRS Study of Local Luminous Infrared Galaxies

We present the first results of our program to study a sample of local luminous infrared galaxies (LIRGs, L_IR = 10^11-10^12 L_sun) with the Spitzer infrared spectrograph (IRS). In these proceedings we investigate the behavior of the 9.7 um silicate feature in LIRGs. As opposed to the extreme silicate absorptions observed in ultraluminous infrared galaxies (ULIRGs, L_IR = 10^12-10^13 L_sun), LIRGs exhibit intermediate silicate absorption features, comparable to those of starburst galaxies. We also find that most of the LIRGs have the minima of the [NeIII]/[NeII] ratio located at their nuclei. It is likely that increased densities in the nuclei are responsible for the smaller nuclear ratios. In the nuclei, it is also possible that the most massive stars are either absent, or still embedded in ultracompact HII regions. Finally we discuss the possible contribution of an AGN to the nuclear mid-IR emission of the galaxy, which in general is low in these local LIRGs.

preprint2010arXiv

Extreme Starbursts in the Local Universe

The "Extreme starbursts in the local universe" workshop was held at the Insituto de Astrofisica de Andalucia in Granada, Spain on 21-25 June 2010. Bearing in mind the advent of a new generation of facilities such as JWST, Herschel, ALMA, eVLA and eMerlin, the aim of the workshop was to bring together observers and theorists to review the latest results. The purpose of the workshop was to address the following issues: what are the main modes of triggering extreme starbursts in the local Universe? How efficiently are stars formed in extreme starbursts? What are the star formation histories of local starburst galaxies? How well do the theoretical simulations model the observations? What can we learn about starbursts in the distant Universe through studies of their local counterparts? How important is the role of extreme starbursts in the hierarchical assembly of galaxies? How are extreme starbursts related to the triggering of AGN in the nuclei of galaxies? Overall, 41 talks and 4 posters with their corresponding 10 minutes short talks were presented during the workshop. In addition, the workshop was designed with emphasis on discussions, and therefore, there were 6 discussion sessions of up to one hour during the workshop. Here is presented a summary of the purposes of the workshop as well as a compilation of the abstracts corresponding to each of the presentations. The summary and conclusions of the workshop along with a description of the future prospects by Sylvain Veilleux can be found in the last section of this document. A photo of the assistants is included.

preprint2009arXiv

Integral field optical spectroscopy of a representative sample of ULIRGs: I. The Data

We present a project aimed at studying the structure, dust distribution, ionization state, and kinematics of a representative sample of 22 ULIRGs. The galaxies in the sample undergo different merger phases (they are evenly divided between pre- and post-coalescence systems) and ionization stages (27% HII, 32% LINER, 18% Seyfert, and 23% mixed classifications) over a wide infrared luminosity range (11.8<Lir/Lsun<12.6), which also includes some galaxies of low luminosity. The main aims of this paper are to present the sample and discuss the structure of the stellar and ionized gas components. Our results imply that evolution caused by a merger is occurring in the ionized gas structure of ULIRGs. The present study relies on the use of integral field optical spectroscopy data obtained with the INTEGRAL instrument at the 4.2 m William Herschel Telescope.

preprint2009arXiv

Integral field optical spectroscopy of a representative sample of ULIRGs: II. Two-dimensional kpc-scale extinction structure

We investigate the two-dimensional kpc-scale structure of the extinction in a representative sample of local ULIRGs using the Halpha/Hbeta line ratio.We use optical integral field spectroscopy obtained with the INTEGRAL instrument at the William Herschel Telescope. Complementary optical and near-IR high angular resolution HST images have also been used. The extinction exhibits a very complex and patchy structure in ULIRGs on kpc scales, from basically transparent regions to others deeply embedded in dust (Av~0.0 to Av~8.0 mag). Nuclear extinction covers a broad range in Av from 0.6 to 6 mag, 69% of the nuclei having Av>2.0 mag. Extinction in the external regions is substantially lower than in the nuclei with 64% of the ULIRGs in the sample having median Av of less than 2 mag for the entire galaxy. While post-coalescence nuclei tend to cluster around Av values of 2 to 3 mag, pre-coalescence nuclei appear more homogeneously distributed over the entire 0.4 mag <Av< 7.7 mag range. For the average extinction (Av~2.0 derived for the ULIRGs of the sample, the ratio of the de-reddened to observed SFR values is 6. The extinction-corrected, Halpha-based SFR ranges from 10 to 300 Msun/yr. For only 28% of the cases the de-reddened SFR is <20 Msun/yr, whereas for the observed SFR this percentage increases to 72%. The IR-based SFR is always higher than the optical-based one, with differences ranging from about 2 to up to 30. The nuclear observed SFR has an average contribution to the total one of 16% for the entire sample. Once corrected for extinction, the average value becomes 31%. Because of mostly extinction effects, the optical (I-band) half-light radius in the sample galaxies is on average a factor 2.3 larger than the corresponding near-IR (H-band) value.

preprint2009arXiv

Near-IR Integral Field Spectroscopy study of the Star Formation and AGN of the LIRG NGC 5135

We present a study of the central 2.3 kpc of NGC 5135, a nearby Luminous Infrared Galaxy (LIRG) with an AGN and circumnuclear starburst. Our main results are based on intermediate spectral resolution (~ 3000-4000) near infrared data taken with the SINFONI integral field spectrograph at the ESO VLT. The ionization of the different phases of the interstellar gas and the complex structures of the star formation have been mapped. Individual regions of interest have been identified and studied in detail. For the first time in this galaxy, we have detected the presence of a high excitation ionization cone centered on the AGN by using the [SiVI] (1.96 micron) line. So far, this structure is the largest reported in the literature for this coronal line, extending (in projection) as far as ~ 600 pc from the galaxy nucleus. In a complex spatial distribution, a variety of mechanisms are driving the gas ionization, including SNe remnant shocks, young stars and AGN photoionization. The excitation of the molecular gas, however, is mainly produced by X-rays and SNe remnant shocks. UV-mechanisms like fluorescence represent a marginal overall contribution to this process, contrary to the expectations we might have for a galaxy with a recent and strong star formation. Our SNe rate estimations from [FeII] (1.64 micron) are in excellent agreement with 6 cm radio emission predictions. Typical SNe rates between 0.01-0.04 yr^-1 were found for individual ~ 200 pc-scale regions, with an overall SNe rate of 0.4-0.5 yr^-1. Even though NGC 5135 has suffered a recent starburst (6-7 Myr ago), the data strongly suggest the presence of a second, older stellar population dominated by red giant/supergiant stars.

preprint2009arXiv

Radio monitoring of NGC 7469: Late time radio evolution of SN 2000ft and the circumnuclear starburst in NGC 7469

We present the results of an eight-year long monitoring of the radio emission from the Luminous Infrared Galaxy (LIRG) NGC 7469, using 8.4 GHz Very Large Array (VLA) observations at 0.3'' resolution. Our monitoring shows that the late time evolution of the radio supernova SN 2000ft follows a decline very similar to that displayed at earlier times of its optically thin phase. The late time radio emission of SN 2000ft is therefore still being powered by its interaction with the presupernova stellar wind, and not with the interstellar medium (ISM). Indeed, the ram pressure of the presupernova wind is ρ_w v_w^2 \approx 7.6E-9 dyn/cm^2, at a supernova age of approximately 2127 days, which is significantly larger than the expected pressure of the ISM around SN 2000ft. At this age, the SN shock has reached a distance r_{sh \approx 0.06 pc, and our observations are probing the interaction of the SN with dense material that was ejected by the presupernova star about 5820 years prior to its explosion. From our VLA monitoring, we estimate that the swept-up mass by the supernova shock after about six years of expansion is \approx 0.29 M_sun, assuming an average expansion speed of the supernova of 10000 km/s. We also searched for recently exploded core-collapse supernovae in our VLA images. Apart from SN 2000ft (S_ν\approx 1760 microJy at its peak, corresponding to 1.1E28 erg/s/Hz, we found no evidence for any other radio supernova (RSN) more luminous than \approx 6.0E26 erg/s/Hz, which suggests that no other Type IIn SN has exploded since 2000 in the circumnuclear starburst of NGC 7469.

preprint2005arXiv

LINER-like Extended Nebulae in ULIRGs: Shocks Generated by Merger Induced Flows

In this work we studied the two-dimensional ionization structure of the circumnuclear and extranuclear regions in a sample of six low-z Ultraluminous Infrared Galaxies using Integral Field Spectroscopy. The ionization conditions in the extranuclear regions of these galaxies (~5-15 kpc) are typical of LINERs as obtained from the Veilleux-Osterbrock line ratio diagnostic diagrams. The range of observed line ratios is best explained by the presence of fast shocks with velocities of 150 to 500 km s^{-1}, while the ionization by an AGN or nuclear starburst is in general less likely. The comparison of the two-dimensional ionization level and velocity dispersion in the extranuclear regions of these galaxies shows a positive correlation, further supporting the idea that shocks are indeed the main cause of ionization. The origin of these shocks is also investigated. Despite the likely presence of superwinds in the circumnuclear regions of these systems, no evidence for signatures of superwinds such as double velocity components are found in the extended extranuclear regions. We consider a more likely explanation for the presence of shocks, the existence of tidally induced large scale gas flows caused by the merging process itself, as evidenced by the observed velocity fields characterized by peak-to-peak velocities of 400 km s^{-1}, and velocity dispersions of up to 200 km s^{-1}.

preprint2001arXiv

HST/WFPC2 Imaging of the QDOT Ultraluminous Infrared Galaxy Sample

We present HST WFPC2 V band imaging for 23 Ultraluminous Infrared Galaxies taken from the QDOT redshift survey. The fraction of sources observed to be interacting is 87%. Most of the merging systems show a number of compact `knots', whose colour and brightness differ substantially from their immediate surroundings. Colour maps for nine of the objects show a non-uniform colour structure. Features include blue regions located towards the centres of merging systems which are likely to be areas of enhanced star formation, and compact red regions which are likely to be dust shrouded starbursts or AGN. The host galaxies of the QSOs in the sample were found to be either interacting systems or ellipticals. Our data shows no evidence that ULIRGs are a simple transition stage between galaxy mergers and QSOs. We propose an alternative model for ULIRGs based on the morphologies in our sample and previous N-body simulations. Under this model ULIRGs as a class are much more diverse than a simple transition between galaxy merger and QSO. The evolution of IR power source and merger morphology in ULIRGs is driven solely by the local environment and the morphologies of the merger progenitors.

preprint1999arXiv

Nuclear Bar, Star Formation and Gas Fueling in the Active Galaxy NGC 4303

A combination of Hubble Space Telescope (HST) WFPC2 and NICMOS images are used to investigate the gas/dust and stellar structure inside the central 300 pc of the nearby active galaxy NGC 4303. The NICMOS H-band (F160W) image reveals a bright core and a nuclear elongated bar-like structure of 250 pc in diameter. The bar is centered on the bright core, and its major axis is oriented in proyection along the spin axis of the nuclear gaseous rotating disk recently detected (Colina & Arribas 1999). The V-H (F606W - F160W) image reveals a complex gas/dust distribution with a two-arm spiral structure of about 225 pc in radius. The southwestern arm is traced by young star-forming knots while the northeastern arm is detected by the presence of dust lanes. These spirals do not have a smooth structure but rather they are made of smaller flocculent spirals or filament-like structures. The magnitudes and colors of the star-forming knots are typical of clusters of young stars with masses of 0.5 to 1 x $10^5 M_{solar}, and ages of 5 to 25 million years. The overall structure of the nuclear spirals as well as the size, number and masses of the star-forming knots are explained in the context of a massive gaseous nuclear disk subject to self-gravitational instabilities and to the gravitational field created by the nuclear bar. According to the model, the gaseous disk has a mass of about 5 x 10^7 M_{solar} inside a radius of 400 pc, the bar has a radius of 150 pc and a pattern speed of about 0.5 Myr^{-1}, and the average mass accretion rate into the core (R < 8 pc) is about 0.01 M_{solar}$ yr^{-1} for about 80 Myr.