Modelling the dust content of spiral galaxies: More dust mass vs. enhanced dust grain emissivity

We present detailed modelling of the spectral energy distribution (SED) of the spiral galaxies NGC 891, NGC 4013, and NGC 5907 in the far-infrared (FIR) and sub-millimeter (submm) wavelengths. The model takes into account the emission of the diffuse dust component, which is heated by the UV and optical radiation fields produced by the stars, as well as the emission produced locally in star forming HII complexes. The radiative transfer simulations of Xilouris et al. (1999) in the optical bands are used to constrain the stellar and dust geometrical parameters, as well as the total amount of dust. We find that the submm emission predicted by our model can not account for the observed fluxes at these wavelengths. We examine two cases, one having more dust embedded in a second thin disk and another allowing for an enhanced submillimeter emissivity of the dust grains. We argue that both cases can equally well reproduce the observed SED. The case of having more dust embedded in a second disk though, is not supported by the near-infrared observations and thus more realistic distributions of the dust (i.e., in spiral arms and clumps) have to be examined in order to better fit the surface brightness of each galaxy.Comment: To appear in the proceedings of: "The Spectral Energy Distribution of Gas-Rich Galaxies: Confronting Models with Data", Heidelberg, 4-8 Oct. 2004, eds. C.C. Popescu and R.J. Tuffs, AIP Conf. Ser., in pres

Is the Galactic submillimeter dust emissivity underestimated?

We present detailed modeling of the spectral energy distribution (SED) of the spiral galaxies NGC 891, NGC 4013, and NGC 5907 in the far-infrared and submm wavelengths. The model takes into account the emission produced by the diffuse dust and the star forming HII complexes. The dust mass is constrained by radiative transfer simulations in the optical (Xilouris et al. 1999). We find that the submm emission predicted by our model cannot account for the observed fluxes. Two scenarios may account for the "missing" submm flux. In the first scenario (Popescu et al. 2000), additional dust (to that derived from the optical, and associated with young stars) is embedded in the galaxy in the form of a thin disk and gives rise to additional submm emission. The other scenario investigates whether the average submm emissivity of the dust grains is higher than the values widely used in Galactic environments. In this case, the dust mass is equal to that derived from the optical observations, and the submm emissivity is treated as a free parameter calculated by fitting our model to the observed SED. We find the submm emissivity value to be ~3 times that often used for our Galaxy. While both scenarios reproduce the observed 850 micron surface brightness, the extra embedded dust model is not supported by the near infrared observations. We, thus, find that the enhanced dust submm emissivity scenario is the most plausible. [abridged]Comment: 12 pages, 10 figures, accepted for publication in Astronomy and Astrophysic

A radio jet drives a molecular and atomic gas outflow in multiple regions within one square kiloparsec of the nucleus of the nearby galaxy IC5063

We analyzed near-infrared data of the nearby galaxy IC5063 taken with the Very Large Telescope SINFONI instrument. IC5063 is an elliptical galaxy that has a radio jet nearly aligned with the major axis of a gas disk in its center. The data reveal multiple signatures of molecular and atomic gas that has been kinematically distorted by the passage of the jet plasma or cocoon within an area of ~1 kpc^2. Concrete evidence that the interaction of the jet with the gas causes the gas to accelerate comes from the detection of outflows in four different regions along the jet trail: near the two radio lobes, between the radio emission tip and the optical narrow-line-region cone, and at a region with diffuse 17.8 GHz emission midway between the nucleus and the north radio lobe. The outflow in the latter region is biconical, centered 240 pc away from the nucleus, and oriented perpendicularly to the jet trail. The diffuse emission that is observed as a result of the gas entrainment or scattering unfolds around the trail and away from the nucleus with increasing velocity. It overall extends for >700 pc parallel and perpendicular to the trail. Near the outflow starting points, the gas has a velocity excess of 600 km/s to 1200 km/s with respect to ordered motions, as seen in [FeII], Pa alpha, or H2 lines. High H2 (1-0) S(3)/S(1) flux ratios indicate non-thermal excitation of gas in the diffuse outflow.Comment: Accepted for publication in Ap

Ultra-luminous Infrared Galaxies in Sloan Digital Sky Survey Data Release 6

Ultra-luminous infrared galaxies (ULIRGs) are interesting objects with dramatic properties. Many efforts have been made to understand the physics of their luminous infrared emission and evolutionary stages. However, a large ULIRG sample is still needed to study the properties of their central black holes (BHs), the BH-host galaxy relation, and their evolution. We identified 308 ULIRGs from the Sloan Digital Sky Survey Data Release 6, and classified them into the NL ULIRGs (with only narrow emission lines) and the Type I ULIRGs (with broad emission lines). About 56% of ULIRGs in our total sample show interaction features, and this percentage is 79% for redshift z < 0.2. Optical identifications of these ULIRGs show that the AGN percentage is at least 49%, and the percentage increases with the infrared luminosity. We found 62 Type I ULIRGs, and estimated their BH masses and velocity dispersions from their optical spectra. Together with known Type I ULIRGs in the literature, a sample of 90 Type I ULIRGs enables us to make a statistical study. We found that the BH masses of Type I ULIRGs are typically smaller than those of PG QSOs, and most Type I ULIRGs follow the M_{BH}-sigma relation. However, some ULIRGs with larger Eddington ratio deviate from this relation, even the line width of the [OIII] narrow line (NL) core or the [SII] line was used as the surrogate of velocity dispersion. This implies that at least some ULIRGs are probably still in the early evolution stage toward QSOs. The anti-correlation between the mass deviation from the M_{BH}-sigma relation and the Eddington ratio supports that the evolution of Type I ULIRGs is probably followed by the building up of the M_{BH}-sigma relation and the evolution to the QSO phase.Comment: 46 pages, 15 figures, 2 tables; published by Ap

X-QUEST: A Comprehensive X-ray Study of Local ULIRGs and QSOs

We present results from the X-ray portion of a multi-wavelength study of local ULIRGs and QSOs called QUEST (Quasar-ULIRG Evolution STudy). The data consist of new and archival X-ray data on 40 ULIRGs and 26 PG QSOs taken with Chandra and XMM-Newton. A combination of traditional and hardness ratio spectral fitting methods is used to characterize the X-ray properties of these objects. The absorption-corrected 2-10 keV to bolometric luminosity ratios of the ULIRGs and PG QSOs suggest that the likelihood for dominant nuclear activity increases along the merger sequence from "cool" ULIRGs, "warm" ULIRGs, infrared-bright QSOs, and infrared-faint QSOs. The starburst dominates the total power in ULIRGs prior to the merger, and this is followed by rapid black hole growth during and after coalescence. These results are in general agreement with those obtained in the mid-infrared with Spitzer and recent numerical simulations.Comment: 63 pages preprint style including 16 figures and 11 tables; accepted by Ap

Cold and warm molecular gas in the outflow of 4C12.50

We present deep observations of the 12CO(1-0) and (3-2) lines in the ultra-luminous infrared and radio galaxy 4C12.50, carried out with the 30m telescope of the Institut de Radioastronomie Millimetrique. Our observations reveal the cold molecular gas component of a warm molecular gas outflow that was previously known from Spitzer Space Telescope data. The 12CO(3-2) profile indicates the presence of absorption at -950 km/s from systemic velocity with a central optical depth of 0.22. Its profile is similar to that of the HI absorption that was seen in radio data of this source. A potential detection of the (0-1) absorption enabled us to place an upper limit of 0.03 on its central optical depth, and to constrain the excitation temperature of the outflowing CO gas to >=65K assuming that the gas is thermalized. If the molecular clouds fully obscure the background millimeter continuum that is emitted by the radio core, the H2 column density is >=1.8*10^22 /cm^2. The outflow then carries an estimated cold H2 mass of at least 4.2*10^3 M_sun along the nuclear line of sight. This mass will be even higher when integrated over several lines of sight, but if it were to exceed 3*10^9 M_sun, the outflow would most likely be seen in emission. Since the ambient cold gas reservoir of 4C12.50 is 1.0*10^10 M_sun, the outflowing-to-ambient mass ratio of the warm gas (37%) could be elevated with respect to that of the cold gas.Comment: A&A letters, in pres

HST NICMOS imaging of z~2, 24 micron-selected Ultraluminous Infrared Galaxies

We present Hubble Space Telescope NICMOS H-band imaging of 33 Ultraluminous Infrared Galaxies (ULIRGs) at z~2 that were selected from the 24 micron catalog of the Spitzer Extragalactic First Look Survey. The images reveal that at least 17 of the 33 objects are associated with interactions. Up to one fifth of the sources in our sample could be minor mergers whereas only 2 systems are merging binaries with luminosity ratio <=3:1, which is characteristic of local ULIRGs. The rest-frame optical luminosities of the sources are of the order 10^10-10^11 L_sun and their effective radii range from 1.4 to 4.9 kpc. The most compact sources are either those with a strong active nucleus continuum or those with a heavy obscuration in the mid-infrared regime, as determined from Spitzer Infra-Red Spectrograph data. The luminosity of the 7.7 micron feature produced by polycyclic aromatic hydrocarbon molecules varies significantly among compact systems whereas it is typically large for extended systems. A bulge-to-disk decomposition performed for the 6 brightest (m_H<20) sources in our sample indicates that they are best fit by disk-like profiles with small or negligible bulges, unlike the bulge-dominated remnants of local ULIRGs. Our results provide evidence that the interactions associated with ultraluminous infrared activity at z~2 can differ from those at z~0.Comment: ApJ, in press. Document revised to match the journal versio

A Deep HST H-Band Imaging Survey of Massive Gas-Rich Mergers. II. The QUEST PG QSOs

We report the results from a deep HST NICMOS H-band imaging survey of 28 z < 0.3 QSOs from the Palomar-Green (PG) sample. This program is part of QUEST (Quasar / ULIRG Evolution STudy) and complements a similar set of data on 26 highly-nucleated ULIRGs presented in Paper I. Our analysis indicates that the fraction of QSOs with elliptical hosts is higher among QSOs with undetected far-infrared (FIR) emission, small infrared excess, and luminous hosts. The hosts of FIR-faint QSOs show a tendency to have less pronounced merger-induced morphological anomalies and larger QSO-to-host luminosity ratios on average than the hosts of FIR-bright QSOs, consistent with late-merger evolution from FIR-bright to FIR-faint QSOs. The spheroid sizes and total host luminosities of the radio-quiet PG QSOs in our sample are statistically indistinguishable from the ULIRG hosts presented in Paper I, while those of radio-loud PG QSOs are systematically larger and more luminous. ULIRGs and PG QSOs with elliptical hosts fall near, but not exactly on, the fundamental plane of inactive spheroids. We confirm the systematic trend noted in Paper I for objects with small (< 2 kpc) spheroids to be up to ~1 mag. brighter than inactive spheroids. The host colors and wavelength dependence of their sizes support the idea that these deviations are due at least in part to non-nuclear star formation. However, the amplitudes of these deviations does not depend on host R-H colors. Taken at face value (i.e., no correction for extinction or the presence of a young stellar population), the H-band spheroid-host luminosities imply BH masses ~5 -- 200 x 10^7 M_sun and sub-Eddington mass accretion rates for both QSOs and ULIRGs. These results are compared with published BH mass estimates derived from other methods. (abridged)Comment: Accepted for publication in the Astrophysical Journal, Vol. 701, August 20 issue. Paper with high-resolution figures can be downloaded at http://www.astro.umd.edu/~veilleux/pubs/nicmos2.pd

Far-IR/Submillimeter Spectroscopic Cosmological Surveys: Predictions of Infrared Line Luminosity Functions for z<4 Galaxies

Star formation and accretion onto supermassive black holes in the nuclei of galaxies are the two most energetic processes in the Universe, producing the bulk of the observed emission throughout its history. We simulated the luminosity functions of star-forming and active galaxies for spectral lines that are thought to be good spectroscopic tracers of either phenomenon, as a function of redshift. We focused on the infrared (IR) and sub-millimeter domains, where the effects of dust obscuration are minimal. Using three different and independent theoretical models for galaxy formation and evolution, constrained by multi-wavelength luminosity functions, we computed the number of star-forming and active galaxies per IR luminosity and redshift bin. We converted the continuum luminosity counts into spectral line counts using relationships that we calibrated on mid- and far-IR spectroscopic surveys of galaxies in the local universe. Our results demonstrate that future facilities optimized for survey-mode observations, i.e., the Space Infrared Telescope for Cosmology and Astrophysics (SPICA) and the Cerro Chajnantor Atacama Telescope (CCAT), will be able to observe thousands of z>1 galaxies in key fine-structure lines, e.g., [SiII], [OI], [OIII], [CII], in a half-square-degree survey, with one hour integration time per field of view. Fainter lines such as [OIV], [NeV] and H_2 (0-0)S1 will be observed in several tens of bright galaxies at 1<z<2, while diagnostic diagrams of active-nucleus vs star-formation activity will be feasible even for normal z~1 galaxies. We discuss the new parameter space that these future telescopes will cover and that strongly motivate their construction.Comment: Accepted for publication in The Astrophysical Journal on 20/10/2011, 17 pages, 13 figure