Local luminous infrared galaxies. I. Spatially resolved observations with the Spitzer infrared spectrograph

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

Local luminous infrared galaxies: Spatially resolved mid-infrared observations with Spitzer/IRS

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

Determining star formation rates for infrared galaxies

We show that measures of star formation rates (SFRs) for infrared galaxies using either single-band 24 μm or extinction-corrected Paα luminosities are consistent in the total infrared luminosity = L(TIR) ~ 10^10 L_☉ range. MIPS 24 μm photometry can yield SFRs accurately from this luminosity upward: SFR(M_☉ yr^–1) = 7.8 × 10^–10 L(24 μm, L_☉) from L(TIR) = 5× 10^9 L_☉ to 10^11 L_☉ and SFR = 7.8 × 10^–10 L(24 μm, L_☉)(7.76 × 10^–11 L(24))^0.048 for higher L(TIR). For galaxies with L(TIR) ≥ 10^10 L_☉, these new expressions should provide SFRs to within 0.2 dex. For L(TIR) ≥ 10^11 L_☉, we find that the SFR of infrared galaxies is significantly underestimated using extinction-corrected Paα (and presumably using any other optical or near-infrared recombination lines). As a part of this work, we constructed spectral energy distribution templates for eleven luminous and ultraluminous purely star forming infrared galaxies and over the spectral range 0.4 μm to 30 cm. We use these templates and the SINGS data to construct average templates from 5 μm to 30 cm for infrared galaxies with L(TIR) = 5× 10^9 to 10^13 L_☉. All of these templates are made available online.

Near-infrared and star-forming properties of local luminous infrared galaxies

We use Hubble Space Telescope (HST) NICMOS continuum and Paα observations to study the near-infrared and star formation properties of a representative sample of 30 local (d ~ 35-75 Mpc) luminous infrared galaxies (LIRGs, infrared [8-1000 μm] luminosities of log L_IR = 11-11.9 L_☉). The data provide spatial resolutions of 25-50 pc and cover the central ~3.3-7.1 kpc regions of these galaxies. About half of the LIRGs show compact (~1-2 kpc) Paα emission with a high surface brightness in the form of nuclear emission, rings, and minispirals. The rest of the sample show Paα emission along the disk and the spiral arms extending over scales of 3-7 kpc and larger. About half of the sample contains H II regions with Hα luminosities significantly higher than those observed in normal galaxies. There is a linear empirical relationship between the mid-IR 24 μm and hydrogen recombination (extinction-corrected Paα) luminosity for these LIRGs, and the H II regions in the central part of M51. This relation holds over more than four decades in luminosity, suggesting that the mid-IR emission is a good tracer of the star formation rate (SFR). Analogous to the widely used relation between the SFR and total IR luminosity of R. Kennicutt, we derive an empirical calibration of the SFR in terms of the monochromatic 24 μm luminosity that can be used for luminous, dusty galaxies.

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

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

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

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

Lurking in ULIRGs: Molecular Gas in local Merging Galaxies

Mergers and interacting galaxies are pivotal to the evolution of galaxies in the universe. They are the sites of prodigious star formation and key to understanding the starburst processes: the physical and chemical properties and the dynamics of the molecular gas. ULIRGs or Ultraluminous Infrared Galaxies are a result of many of these mergers. They host extreme starbursts, AGNs, and mergers. They are the perfect laboratory to probe the connection between starbursts, black hole accretion and mergers and to further our understanding of star formation and merging.

NGC 6240 and Arp 220 can be considered the founding members of this very active class of objects. They are in different stages of merging and hence are excellent case studies to further our understanding about the merging process. We have imaged the dense star-forming regions of these galaxies at sub-arcsec resolution with CARMA C and B Configurations (2" and 0.5 - 0.8"). Multi-band imaging allows excitation analysis of HCN, HCO⁺, HNC, and CS along with CO transitions to constrain the properties of the gas. Our dataset is unique in that we have observed these lines at similar resolutions and high sensitivity which can be used to derive line ratios of faint high excitation lines.

Arp 220 has not had confirmed X-ray AGN detections for either nuclei. However, our observations indicate HCN/HNC ratios consistent with the chemistry of X-ray Dominated Regions (XDRs) -- a likely symptom of AGN. We calculated the molecular Hydrogen densities using each of the molecular species and conclude that assuming abundances of HNC and HCO⁺ similar to those in galactic sources are incorrect in the case of ULIRGs. The physical conditions in the dense molecular gas in ULIRGs alter these abundances. The derived H2 volume densities are ~ 5 x 10⁴ cm^-3 in both Arp 220 nuclei and ~ 10⁴ cm^-3 in NGC 6240.

Nuclear and extended spectra of NGC 1068-I. Hints from near-infrared spectroscopy

We report the first simultaneous zJHK spectroscopy on the archetypical Seyfert 2 galaxy NGC 1068 covering the wavelength region 0.9-2.4 mu m. The slit, aligned in the north-south direction and centred in the optical nucleus, maps a region 300 pc in radius at subarcsec resolution, with a spectral resolving power of 360 km s-1. This configuration allows us to study the physical properties of the nuclear gas including that of the north side of the ionization cone, map the strong excess of continuum emission in the K band and attributed to dust and study the variations, both in flux and profile, in the emission lines. Our results show the following. (1) Mid- to low-ionization emission lines are split into two components, whose relative strengths vary with the position along the slit and seem to be correlated with the jet. (2) The coronal lines are single-peaked and are detected only in the central few hundred of pc from the nucleus. (3) The absorption lines indicate the presence of intermediate age stellar population, which might be a significant contributor to the continuum in the near-IR spectra. (4) Through some simple photoionization models we find photoionization as the main mechanism powering the emitting gas. (5) Calculations using stellar features point to a mass concentration inside the 100-200 pc of about 1010 M(circle dot).

Nuclear and extended spectra of NGC 1068-II. Near-infrared stellar population synthesis

We performed stellar population synthesis on the nuclear and extended regions of NGC 1068 by means of near-infrared spectroscopy to disentangle their spectral energy distribution components. This is the first time that such a technique is applied to the whole 0.8-2.4 mu m wavelength interval in this galaxy. NGC 1068 is one of the nearest and probably the most studied Seyfert 2 galaxy, becoming an excellent laboratory to study the interaction between black holes, the jets that they can produce and the medium in which they propagate. Our main result is that traces of young stellar population are found at similar to 100 pc south of the nucleus. The contribution of a power-law continuum in the centre is about 25 per cent, which is expected if the light is scattered from a Seyfert 1 nucleus. We find peaks in the contribution of the featureless continuum about 100-150 pc from the nucleus on both sides. They might be associated with regions where the jet encounters dense clouds. Further support to this scenario is given by the peaks of hot dust distribution found around these same regions and the H(2) emission-line profile, leading us to propose that the peaks might be associated to regions where stars are being formed. Hot dust also has an important contribution to the nuclear region, reinforcing the idea of the presence of a dense, circumnuclear torus in this galaxy. Cold dust appears mostly in the south direction, which supports the view that the south-west emission is behind the plane of the galaxy and is extinguished very likely by dust in the plane. Intermediate-age stellar population contributes significantly to the continuum, especially in the inner 200 pc.

The effect of local galaxy surface density on star formation for HI selected galaxies

We present the result of investigations into two theories to explain the star formation rate (SFR)-density relationship. For regions of high galaxy density, either there are fewer star-forming galaxies or galaxies capable of forming stars are present but some physical process is suppressing their star formation. We use H I Parkes All-Sky Survey's (HIPASS) HI detected galaxies and infrared and radio fluxes to investigate SFRs and efficiencies with respect to local surface density. For nearby (vel < 10 000 km s(-1)) H I galaxies, we find a strong correlation between H I mass and SFR. The number of H I galaxies decreases with increasing local surface density. For H I galaxies (1000 < vel < 6000 km s(-1)), there is no significant change in the SFR or the efficiency of star formation with respect to local surface density. We conclude that the SFR-density relation is due to a decrease in the number of H I star-forming galaxies in regions of high galaxy density and not to the suppression of star formation.

IRAS 00183-7111: ALMA and X-ray view of an ultra luminous infrared galaxy

The aim of this Thesis work is to study the multi-frequency properties of the Ultra Luminous Infrared Galaxy (ULIRG) IRAS 00183-7111 (I00183) at z = 0.327, connecting ALMA sub-mm/mm observations with those at high energies in order to place constraints on the properties of its central power source and verify whether the gas traced by the CO may be responsible for the obscuration observed in X-rays. I00183 was selected from the so-called Spoon diagnostic diagram (Spoon et al. 2007) for mid-infrared spectra of infrared galaxies based on the equivalent width of the 6.2 μm Polycyclic Aromatic Hydrocarbon (PAH) emission feature versus the 9.7 μm silicate strength. Such features are a powerful tool to investigate the contribution of star formation and AGN activity in this class of objects. I00183 was selected from the top-left region of the plot where the most obscured sources, characterized by a strong Si absorption feature, are located. To link the sub-mm/mm to the X-ray properties of I00183, ALMA archival Cycle 0 data in Band 3 (87 GHz) and Band 6 (270 GHz) have been calibrated and analyzed, using CASA software. ALMA Cycle 0 was the Early Science program for which data reprocessing is strongly suggested. The main work of this Thesis consisted in reprocessing raw data to provide an improvement with respect to the available archival products and results, which were obtained using standard procedures. The high-energy data consists of Chandra, XMM-Newton and NuSTAR observations which provide a broad coverage of the spectrum in the energy range 0.5 − 30 keV. Chandra and XMM archival data were used, with an exposure time of 22 and 22.2 ks, respectively; their reduction was carried out using CIAO and SAS software. The 100 ks NuSTAR are still private and the spectra were obtained by courtesy of the PI (K. Iwasawa). A detailed spectral analysis was done using XSPEC software; the spectral shape was reproduced starting from simple phenomenological models, and then more physical models were introduced to account for the complex mechanisms that involve this source. In Chapter 1, an overview of the scientific background is discussed, with a focus on the target, I00183, and the Spoon diagnostic diagram, from which it was originally selected. In Chapter 2, the basic principles of interferometry are briefly introduced, with a description of the calibration theory applied to interferometric observations. In Chapter 3, ALMA and its capabilities, both current and future, are shown, explaining also the complex structure of the ALMA archive. In Chapter 4, the calibration of ALMA data is presented and discussed, showing also the obtained imaging products. In Chapter 5, the analysis and discussion of the main results obtained from ALMA data is presented. In Chapter 6, the X-ray observations, data reduction and spectral analysis are reported, with a brief introduction to the basic principle of X-ray astronomy and the instruments from which the observations were carried out. Finally, the overall work is summarized, with particular emphasis on the main obtained results and the possible future perspectives.

Search for sub-mm, mm and radio continuum emission from extremely red objects

We present the results of sub-mm, mm (850 mum, 450 mum and 1250 mum) and radio (1.4 and 4.8 GHz) continuum observations of a sample of 27 K-selected Extremely Red Objects, or EROs, (14 of which form a complete sample with K < 20 and I - K > 5) aimed at detecting dusty starbursts, deriving the fraction of UltraLuminous Infrared Galaxies (ULIGs) in ERO samples, and constraining their redshifts using the radio-FIR correlation. One ERO was tentatively detected at 1250 mum and two were detected at 1.4 GHz, one of which has a less secure identification as an ERO counterpart. Limits on their redshifts and their star forming properties are derived and discussed. We stacked the observations of the undetected objects at 850 mum, 1250 mum and 4.8 GHz in order to search for possible statistical emission from the ERO population as a whole, but no significant detections were derived either for the whole sample or as a function of the average NIR colours. These results strongly suggest that the dominant population of EROs with K < 20 is not comprised of ULIGs like HR 10, but is probably made of radio-quiet ellipticals and weaker starburst galaxies with L < 10(12) L . and SFR < few hundred M. yr(-1).

Supernovae and radio transients in M82

We present optical and near-infrared (NIR) photometry and NIR spectroscopy of SN 2004am, the only optically detected supernova (SN) in M82. These demonstrate that SN 2004am was a highly reddened Type II-P SN similar to the low-luminosity Type II-P events such as SNe 1997D and 2005cs. We show that SN 2004am was located coincident with the obscured super star cluster M82-L, and from the cluster age infer a progenitor mass of 12{^{+ 7}_{- 3}} M⊙. In addition to this, we present a high spatial resolution Gemini-North Telescope K-band adaptive optics image of the site of SN 2008iz and a second transient of uncertain nature, both detected so far only at radio wavelengths. Using image subtraction techniques together with archival data from the Hubble Space Telescope, we are able to recover a NIR transient source coincident with both objects. We find the likely extinction towards SN 2008iz to be not more than AV ˜ 10. The nature of the second transient remains elusive and we regard an extremely bright microquasar in M82 as the most plausible scenario.

The Mice at play in the CALIFA survey. A case study of a gas-rich major merger between first passage and coalescence

We present optical integral field spectroscopy (IFS) observations of the Mice, a major merger between two massive (≳10^11 M_⊙) gas-rich spirals NGC 4676A and B, observed between first passage and final coalescence. The spectra provide stellar and gas kinematics, ionised gas properties, and stellar population diagnostics, over the full optical extent of both galaxies with ~1.6 kpc spatial resolution. The Mice galaxies provide a perfect case study that highlights the importance of IFS data for improving our understanding of local galaxies. The impact of first passage on the kinematics of the stars and gas has been significant, with strong bars most likely induced in both galaxies. The barred spiral NGC 4676B exhibits a strong twist in both its stellar and ionised gas disk. The edge-on disk galaxy NGC 4676A appears to be bulge free, with a strong bar causing its “boxy” light profile. On the other hand, the impact of the merger on the stellar populations has been minimal thus far. By combining the IFS data with archival multiwavelength observations we show that star formation induced by the recent close passage has not contributed significantly to the total star formation rate or stellar mass of the galaxies. Both galaxies show bicones of high ionisation gas extending along their minor axes. In NGC 4676A the high gas velocity dispersion and Seyfert-like line ratios at large scaleheight indicate a powerful outflow. Fast shocks (vs ~ 350 km s^-1) extend to ~6.6 kpc above the disk plane. The measured ram pressure (P/k = 4.8 × 10^6 K cm^-3) and mass outflow rate (~8−20 M_⊙ yr^-1) are similar to superwinds from local ultra-luminous infrared galaxies, although NGC 4676A only has a moderate infrared luminosity of 3 × 10^10 L_⊙. Energy beyond what is provided by the mechanical energy of the starburst appears to be required to drive the outflow. Finally, we compare the observations to mock kinematic and stellar population maps extracted from a hydrodynamical merger simulation. The models show little enhancement in star formation during and following first passage, in agreement with the observations. We highlight areas where IFS data could help further constrain the models.

Role of galaxy mergers in cosmic star formation history

We present a morphology study of intermediate-redshift (0.2 < z < 1.2) luminous infrared galaxies (LIRGs) and general field galaxies in the GOODS fields using a revised asymmetry measurement method optimized for deep fields. By taking careful account of the importance of the underlying sky-background structures, our new method does not suffer from systematic bias and offers small uncertainties. By redshifting local LIRGs and low-redshift GOODS galaxies to different higher redshifts, we have found that the redshift dependence of the galaxy asymmetry due to surface-brightness dimming is a function of the asymmetry itself, with larger corrections for more asymmetric objects. By applying redshift-, infrared (IR)-luminosity- and optical-brightness-dependent asymmetry corrections, we have found that intermediate-redshift LIRGs generally show highly asymmetric morphologies, with implied merger fractions ~50% up to z = 1.2, although they are slightly more symmetric than local LIRGs. For general field galaxies, we find an almost constant relatively high merger fraction (20%-30%). The B-band luminosity functions (LFs) of galaxy mergers are derived at different redshifts up to z = 1.2 and confirm the weak evolution of the merger fraction after breaking the luminosity-density degeneracy. The IR LFs of galaxy mergers are also derived, indicating a larger merger fraction at higher IR luminosity. The integral of the merger IR LFs indicates a dramatic evolution of the merger-induced IR energy density [(1 + z)^~(5-6)], and that galaxy mergers start to dominate the cosmic IR energy density at z greater than or ~ 1.