A prominent jet-driven outflow of CO(2-1) molecular gas is found along the kinematic minor axis of the Seyfert 2 galaxy ESO 420-G13, at a distance of $340-600\, \rm{pc}$ from the nucleus. The wind morphology resembles a characteristic funnel shape, formed by a highly collimated filamentary emission at the base, likely tracing the jet propagation through a tenuous medium, until a bifurcation point at $440\, \rm{pc}$ where the jet hits a dense molecular core and shatters, dispersing the molecular gas into several clumps and filaments within the expansion cone. We also trace the jet in ionised gas within the inner $\lesssim 340\, \rm{pc}$ using the [NeII]$_{\rm 12.8 \mu m}$ line emission, where the molecular gas follows a circular rotation pattern. The wind outflow carries a mass of $\sim 8 \times 10^6\, \rm{M_\odot}$ at an average wind projected speed of $\sim 160\, \rm{km\,s^{-1}}$, which implies a mass outflow rate of $\sim 14\, \rm{M_\odot\,yr^{-1}}$. Based on the structure of the outflow and the budget of energy and momentum, we discard radiation pressure from the active nucleus, star formation, and supernovae as possible launching mechanisms. ESO 420-G13 is the second case after NGC 1377 where the presence of a previously unknown jet is revealed due to its interaction with the interstellar medium, suggesting that unknown jets in feeble radio nuclei might be more common than expected. Two possible jet-cloud configurations are discussed to explain the presence of an outflow at such distance from the AGN. The outflowing gas will likely not escape, thus a delay in the star formation rather than quenching is expected from this interaction, while the feedback effect would be confined within the central few hundred parsecs of the galaxy.

The nearby system 4C12.50, also known as IRAS 13451+1217 and PKS 1345+12, is a merger of gas-rich galaxies with infrared and radio activity. It has a perturbed interstellar medium (ISM) and a dense configuration of gas and dust around the nucleus. The radio emission at small ($\sim$100 pc) and large ($\sim$100 kpc) scales, as well as the large X-ray cavity in which the system is embedded, are indicative of a jet that could have affected the ISM. We carried out observations of the CO(1-0), (3-2), and (4-3) lines with the Atacama Large Millimeter Array (ALMA) to determine basic properties (i.e., extent, mass, and excitation) of the cold molecular gas in this system, including its already-known wind. The CO emission reveals the presence of gaseous streams related to the merger, which result in a small ($\sim$4kpc-wide) disk around the western nucleus. The disk reaches a rotational velocity of 200 $kms^{-1}$ , and has a mass of 3.8($\pm$0.4)$\times$10${^9}M_{\odot}$. It is truncated at a gaseous ridge north of the nucleus that is bright in [O III]. Regions with high-velocity CO emission are seen at signal-to-noise ratios of between 3 and 5 along filaments that radially extend from the nucleus to the ridge and that are bright in [O III] and stellar emission. A tentative wind detection is also reported in the nucleus and in the disk. The molecular gas speed could be as high as 2200 $kms^{-1}$ and the total wind mass could be as high as 1.5($\pm$0.1)$\times$10$^9M_{\odot}$. Energetically, it is possible that the jet, assisted by the radiation pressure of the active nucleus or the stars, accelerated clouds inside an expanding bubble.

We report ALMA observations of CO(3-2) emission in the Seyfert galaxy NGC613, at a spatial resolution of 17pc, as part of our NUclei of GAlaxies sample. Our aim is to investigate the morphology and dynamics of the gas inside the central kpc, and to probe nuclear fueling and feedback phenomena. The morphology of CO(3-2) line emission reveals a 2-arm trailing nuclear spiral at $r\lesssim$100pc and a circumnuclear ring at ~350pc radius, that is coincident with the star-forming ring seen in the optical images. The molecular gas in the galaxy disk is in a remarkably regular rotation, however, the kinematics in the nuclear region is very skewed. The nuclear spectrum of CO and dense gas tracers HCN(4-3), HCO+(4-3), and CS(7-6) show broad wings up to $\pm$300km/s, associated with a molecular outflow emanating from the nucleus (r~25pc). We derive a molecular outflow mass $M_{out}$=2x10$^6$M$_\odot$ and a mass outflow rate of $\dot{M}_{out}=$27$\rm M_\odot yr^{-1}$. The molecular outflow energetics exceed the values predicted by AGN feedback models: the kinetic power of the outflow corresponds to 20%L_AGN and the momentum rate is $\dot{M}_{out}v\sim400L_{AGN}/c$. The outflow is mainly boosted by the AGN through entrainment by the radio jet, but given the weak nuclear activity of NGC613, we might be witnessing a fossil outflow, resulted from a strong past AGN that now has already faded. The nuclear trailing spiral observed in CO emission is inside the ILR ring of the bar. We compute the gravitational torques exerted in the gas to estimate the efficiency of the angular momentum exchange. The gravity torques are negative from 25 to 100pc and the gas loses its angular momentum in a rotation period, providing evidence of a highly efficient inflow towards the center. This phenomenon shows that the massive central black hole has a significant dynamical influence on the gas, triggering its fueling.

We report ALMA observations of CO(3-2) emission in a sample of 7 Seyfert/LINER galaxies at the unprecedented spatial resolution of 0.1'' = 4-9 pc. Our aim is to explore the close environment of AGN, and the dynamical structures leading to its fueling. The selected galaxies host low-luminosity AGN, and have a wide range of activity types, and barred or ringed morphologies. The observed maps reveal the existence of circum-nuclear disk structures, defined by their morphology and decoupled kinematics. We call these structures "molecular tori", even though they appear often as disks, without holes in the center. They have varied orientations, unaligned with their host galaxy. The radius of the tori ranges from 6 to 27 pc, and their mass from 0.7 10$^7$ to 3.9 10$^7$ Msun. At larger scale, the gas is always piled up in a few 100~pc scale resonant rings, that play the role of a reservoir to fuel the nucleus. In some cases, a trailing spiral is observed inside the ring, providing evidence for feeding processes. Most frequently, the torus and the AGN are slightly off-centered, with respect to the bar-resonant ring position, implying that the black hole is wandering by a few 10~pc amplitude around the center of mass of the galaxy. Our spatial resolution allows us to measure gas velocities inside the sphere of influence of the central black holes. By fitting the observations with different simulated cubes, varying the torus inclination and the black hole mass, it is possible to estimate the mass of the central black hole, which is in general difficult for such late-type galaxies, with only a pseudo-bulge.

We present a detailed study of the molecular gas in the fast AGN-driven outflow in the nearby radio-loud Seyfert galaxy IC 5063. Using ALMA observations of a number of tracers (12CO(1-0), 12CO(2-1), 12CO(3-2), 13CO(2-1) and HCO+(4-3)), we map the differences in excitation, density and temperature of the gas. The results show that in the immediate vicinity of the radio jet, a fast outflow, with velocities up to 800 km/s, is occurring of which the gas has high excitation temperatures in the range 30-55 K, demonstrating the direct impact of the jet on the ISM. The relative brightness of the CO lines show that the outflow is optically thin. We estimate the mass of the molecular outflow to be 1.2 x 10^6 Msol and likely to be a factor 2-3 larger. This is similar to that of the outflow of atomic gas, but much larger than that of the ionised outflow, showing that the outflow is dominated by cold gas. The total mass outflow rate we estimate to be ~12 Msol/yr. The mass of the outflow is much smaller than the total gas mass of the ISM of IC 5063. Therefore, although the influence of the radio jet is very significant in the inner regions, globally speaking the impact will be very modest. We use RADEX modelling to explore the physical conditions of the molecular gas in the outflow. Models with the outflowing gas being quite clumpy give the most consistent results and our preferred solutions have kinetic temperatures in the range 20-100 K and densities between 10^5 and 10^6 cm^-3. The resulting pressures are 10^6-10^7.5 K cm^-3, about two orders of magnitude higher than in the outer quiescent disk. The results strongly suggest that the outflow is driven by the radio jet expanding into a clumpy medium, creating a cocoon of gas which is pushed away from the jet axis resulting in a lateral outflow, very similar to what is predicted by numerical simulations.

High resolution submm observations are important in probing the morphology, column density and dynamics of obscured active galactic nuclei (AGNs). With high resolution (0.06 x 0.05) ALMA 690 GHz observations we have found bright (TB >80 K) and compact (FWHM 10x7 pc) CO 6-5 line emission in the nucleus of the extremely radio-quiet galaxy NGC1377. The integrated CO 6-5 intensity is aligned with the previously discovered jet/outflow of NGC1377 and is tracing the dense (n>1e4 cm-3), hot gas at the base of the outflow. The velocity structure is complex and shifts across the jet/outflow are discussed in terms of jet-rotation or separate, overlapping kinematical components. High velocity gas (deltaV +-145 km/s) is detected inside r<2-3 pc and we suggest that it is emerging from an inclined rotating disk or torus of position angle PA=140+-20 deg with a dynamical mass of approx 3e6 Msun. This mass is consistent with that of a supermassive black hole (SMBH), as inferred from the M-sigma relation. The gas mass of the proposed disk/torus constitutes <3% of the nuclear dynamical mass. In contrast to the intense CO 6-5 line emission, we do not detect dust continuum with an upper limit of S(690GHz)<2mJy. The corresponding, 5 pc, H2 column density is estimated to N(H2)<3e23 cm-2, which is inconsistent with a Compton Thick (CT) source. We discuss the possibility that CT obscuration may be occuring on small (subparsec) or larger scales. From SED fitting we suggest that half of the IR emission of NGC1377 is nuclear and the rest (mostly the far-infrared (FIR)) is more extended. The extreme radio quietness, and the lack of emission from other star formation tracers, raise questions on the origin of the FIR emission. We discuss the possibility that it is arising from the dissipation of shocks in the molecular jet/outflow or from irradiation by the nuclear source along the poles.

Using CO (4-3) and (2-1) Atacama Large Millimeter Array (ALMA) data, we prove that the molecular gas in the jet-driven winds of the galaxy IC5063 is more highly excited than the rest of the molecular gas in the disk of the same galaxy. On average, the CO(4-3)/CO(2-1) flux ratio is 1 for the disk and 5 for the jet accelerated or impacted gas. Spatially-resolved maps reveal that in regions associated with winds, the CO(4-3)/CO(2-1) flux ratio significantly exceeds the upper limit of 4 for optically thick gas. It frequently takes values between 5 and 11, and it occasionally further approaches the upper limit of 16 for optically thin gas. Excitation temperatures of 30-100 K are common for the molecules in these regions. If all of the outflowing molecular gas is optically thin, at 30-50 K, then its mass is 2*10^6 M_sun. This lower mass limit is an order of magnitude below the mass derived from the CO(2-1) flux in the case of optically thick gas. Molecular winds can thus be less massive, but more easily detectable at high z than they were previously thought to be.

We present a coherent database of spectroscopic observations of far-IR fine-structure lines from the Herschel/PACS archive for a sample of 170 local AGN, plus a comparison sample of 20 starburst galaxies and 43 dwarf galaxies. Published Spitzer/IRS and Herschel/SPIRE line fluxes are included to extend our database to the full 10-600 $\mu m$ spectral range. The observations are compared to a set of CLOUDY photoionisation models to estimate the above physical quantities through different diagnostic diagrams. We confirm the presence of a stratification of gas density in the emission regions of the galaxies, which increases with the ionisation potential of the emission lines. The new [OIV]25.9$\mu m$/[OIII]88$\mu m$ vs [NeIII]15.6$\mu m$/[NeII]12.8$\mu m$ diagram is proposed as the best diagnostic to separate: $i)$ AGN activity from any kind of star formation; and $ii)$ low-metallicity dwarf galaxies from starburst galaxies. Current stellar atmosphere models fail to reproduce the observed [OIV]25.9$\mu m$/[OIII]88$\mu m$ ratios, which are much higher when compared to the predicted values. Finally, the ([NeIII]15.6$\mu m$ + [NeII]12.8$\mu m$)/([SIV]10.5$\mu m$ + [SIII]18.7$\mu m$) ratio is proposed as a promising metallicity tracer to be used in obscured objects, where optical lines fail to accurately measure the metallicity. The diagnostic power of mid- to far-infrared spectroscopy shown here for local galaxies will be of crucial importance to study galaxy evolution during the dust-obscured phase at the peak of the star formation and black-hole accretion activity ($1 < z < 4$). This study will be addressed by future deep spectroscopic surveys with present and forthcoming facilities such as JWST, ALMA, and SPICA.

With high resolution (0."25 x 0."18) ALMA CO 3-2 observations of the nearby (D=21 Mpc, 1"=102 pc), extremely radio-quiet galaxy NGC1377, we have discovered a high-velocity, very collimated nuclear outflow which we interpret as a molecular jet with a projected length of +-150 pc. Along the jet axis we find strong velocity reversals where the projected velocity swings from -150 km/s to +150 km/s. A simple model of a molecular jet precessing around an axis close to the plane of the sky can reproduce the observations. The velocity of the outflowing gas is difficult to constrain due to the velocity reversals but we estimate it to be between 240 and 850 km/s and the jet to precess with a period P=0.3-1.1 Myr. The CO emission is clumpy along the jet and the total molecular mass in the high-velocity (+-(60 to 150 km/s)) gas lies between 2e6 Msun (light jet) and 2e7 Msun (massive jet). There is also CO emission extending along the minor axis of NGC1377. It holds >40% of the flux in NGC1377 and may be a slower, wide-angle molecular outflow which is partially entrained by the molecular jet. We discuss the driving mechanism of the molecular jet and suggest that it is either powered by a very faint radio jet or by an accretion disk-wind similar to those found towards protostars. The nucleus of NGC1377 harbours intense embedded activity and we detect emission from vibrationally excited HCN J=4-3 v_2=1f which is consistent with hot gas and dust. We find large columns of H2 in the centre of NGC1377 which may be a sign of a high rate of recent gas infall. The dynamical age of the molecular jet is short (<1 Myr), which could imply that it is young and consistent with the notion that NGC1377 is caught in a transient phase of its evolution. However, further studies are required to determine the age of the molecular jet, its mass and the role it is playing in the growth of the nucleus of NGC1377.

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.

Gas can be violently stripped from their galaxy disks in rich clusters, and be dispersed over 100kpc-scale tails or plumes. Young stars have been observed in these tails, suggesting they are formed in situ. This will contribute to the intracluster light, in addition to tidal stripping of old stars. We want to quantify the efficiency of intracluster star formation. We present CO(1--0) and CO(2--1) observations, made with the IRAM-30m telescope, towards the ram-pressure stripped tail northeast of NGC4388 in Virgo. HII regions found all along the tails, together with dust patches have been targeted. We detect molecular gas in 4 positions along the tail, with masses between 7x10$^5$ to 2x10$^6$ M$_\odot$. Given the large distance from the NGC 4388 galaxy, the molecular clouds must have formed in situ, from the HI gas plume. We compute the relation between surface densities of star formation and molecular gas in these regions, and find that the star formation has very low efficiency. The corresponding depletion time of the molecular gas can be up to 500 Gyr and more. Since this value exceeds a by far Hubble time, this gas will not be converted into stars, and will stay in a gaseous phase to join the intracluster medium.

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

We present a comparison of the molecular gas properties in the outflow vs. in the ambient medium of the local prototype radio-loud and ultraluminous-infrared galaxy 4C12.50 (IRAS13451+1232), using new data from the IRAM Plateau de Bure interferometer and 30m telescope, and the Herschel space telescope. Previous H_2 (0-0) S(1) and S(2) observations with the Spitzer space telescope had indicated that the warm (~400K) molecular gas in 4C12.50 is made up of a 1.4(+-0.2)x10^8 M_sun ambient reservoir and a 5.2(+-1.7)x10^7 M_sun outflow. The new CO(1-0) data cube indicates that the corresponding cold (25K) H_2 gas mass is 1.0(+-0.1)x10^10 M_sun for the ambient medium and <1.3x10^8 M_sun for the outflow, when using a CO-intensity-to-H_2-mass conversion factor alpha of 0.8 M_sun /(K km/s pc^2). The combined mass outflow rate is high, 230-800 M_sun/yr, but the amount of gas that could escape the galaxy is low. A potential inflow of gas from a 3.3(+-0.3)x10^8 M_sun tidal tail could moderate any mass loss. The mass ratio of warm-to-cold molecular gas is >= 30 times higher in the outflow than in the ambient medium, indicating that a non-negligible fraction of the accelerated gas is heated to temperatures at which star formation is inefficient. This conclusion is robust against the use of different alpha factor values, and/or different warm gas tracers (H_2 vs. H_2 plus CO): with the CO-probed gas mass being at least 40 times lower at 400K than at 25K, the total warm-to-cold mass ratio is always lower in the ambient gas than in the entrained gas. Heating of the molecular gas could facilitate the detection of new outflows in distant galaxies by enhancing their emission in intermediate rotational number CO lines.

The main energy-generating mechanisms in galaxies are black hole (BH) accretion and star formation (SF) and the interplay of these processes is driving the evolution of galaxies. MIR/FIR spectroscopy are able to distinguish between BH accretion and SF, as it was shown in the past by infrared spectroscopy from the space by the Infrared Space Observatory and Spitzer. Spitzer and Herschel spectroscopy together can trace the AGN and the SF components in galaxies, with extinction free lines, almost only in the local Universe, except for a few distant objects. One of the major goals of the study of galaxy evolution is to understand the history of the luminosity source of galaxies along cosmic time. This goal can be achieved with far-IR spectroscopic cosmological surveys. SPICA in combination with ground based large single dish submillimeter telescopes, such as CCAT, will offer a unique opportunity to do this. We use galaxy evolution models linked to the observed MIR-FIR counts (including Herschel) to predict the number of sources and their IR lines fluxes, as derived from observations of local galaxies. A shallow survey in an area of 0.5 square degrees, with a typical integration time of 1 hour per pointing, will be able to detect thousands of galaxies in at least three emission lines, using SAFARI, the far-IR spectrometer onboard of SPICA.

We present new stellar velocity dispersion measurements for four luminous quasars with the NIFS instrument and the ALTAIR laser guide star adaptive optics system on the Gemini North 8-m telescope. Stellar velocity dispersion measurements and measurements of the supermassive black hole masses in luminous quasars are necessary to investigate the coevolution of black holes and galaxies, trace the details of accretion, and probe the nature of feedback. We find that higher-luminosity quasars with higher-mass black holes are not offset with respect to the MBH-sigma relation exhibited by lower-luminosity AGNs with lower-mass black holes, nor do we see correlations with galaxy morphology. As part of this analysis, we have recalculated the virial products for the entire sample of reverberation-mapped AGNs and used these data to redetermine the mean virial factor hfi that places the reverberation data on the quiescent M_BH-sigma relation. With our updated measurements and new additions to the AGN sample, we obtain <f> = 4.31 +/- 1.05, which is slightly lower than, but consistent with, most previous determinations.

We have compiled a large sample of 151 high redshift (z=0.5-4) galaxies selected at 24 microns (S24>100 uJy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the mid-infrared spectrum into contributions from star formation and activity in the galactic nuclei. In addition, we have a wealth of photometric data from Spitzer IRAC/MIPS and Herschel PACS/SPIRE. We explore how effective different infrared color combinations are at separating our mid-IR spectroscopically determined active galactic nuclei from our star forming galaxies. We look in depth at existing IRAC color diagnostics, and we explore new color-color diagnostics combining mid-IR, far-IR, and near-IR photometry, since these combinations provide the most detail about the shape of a source's IR spectrum. An added benefit of using a color that combines far-IR and mid-IR photometry is that it is indicative of the power source driving the IR luminosity. For our data set, the optimal color selections are S250/S24 vs. S8.0/S3.6 and S100/S24 vs. S8.0/S3.6; both diagnostics have ~10% contamination rate in the regions occupied primarily by star forming galaxies and active galactic nuclei, respectively. Based on the low contamination rate, these two new IR color-color diagnostics are ideal for estimating both the mid-IR power source of a galaxy when spectroscopy is unavailable and the dominant power source contributing to the IR luminosity. In the absence of far-IR data, we present color diagnostics using the WISE mid-IR bands which can efficiently select out high z (z~2) star forming galaxies.

We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J=4-3 through J=13-12, 6 H2O rotational lines, and [CI] and [NII] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We applied both C and J shock models to the H2 v=1-0 S(1) and v=2-1 S(1) lines and the CO rotational ladder. The CO ladder is best reproduced by a model with shock velocity v_s=10 km s^-1 and a pre-shock density n_H=5 * 10^4 cm^-3. We find that the solution best fitting the H2 lines is degenerate: The shock velocities and number densities range between v_s = 17 - 47 km s^-1 and n_H=10^7 - 5 * 10^4 cm^-3, respectively. The H2 lines thus need a much more powerful shock than the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 km s^-1) shocks while only a small fraction (< 1 percent) of the ISM is exposed to the high velocity shocks. This implies that the gas is rapidly loosing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks.

Many theoretical models require powerful active galactic nuclei (AGNs) to suppress star formation in distant galaxies and reproduce the observed properties of today's massive galaxies. A recent study based on Herschel-SPIRE submillimetre observations claimed to provide direct support for this picture, reporting a significant decrease in the mean star-formation rates (SFRs) of the most luminous AGNs (Lx>10^44 erg/s) at z=1-3 in the Chandra Deep Field-North (CDF-N). In this letter we extend these results using Herschel-SPIRE 250um data in the COSMOS and CDF-S fields to achieve an order of magnitude improvement in the number of sources at Lx>10^44 erg/s. On the basis of our analysis, we find no strong evidence for suppressed star formation in Lx>10^44 erg/s AGNs at z=1-3. The mean SFRs of the AGNs are constant over the broad X-ray luminosity range of Lx~10^43-10^45 erg/s (with mean SFRs consistent with typical star-forming galaxies at z~2; <SFRs>~100-200 Msol/yr). We suggest that the previous CDF-N results were likely due to low number statistics. We discuss our results in the context of current theoretical models.

We present here a new spectral energy distribution (SED) fitting approach that we adopt to select radio-excess sources amongst distant star-forming galaxies in the GOODS-Herschel (North) field and to reveal the presence of hidden, highly obscured AGN. Through extensive SED analysis of 458 galaxies with radio 1.4 GHz and mid-IR 24 um detections using some of the deepest Chandra X-ray, Spitzer and Herschel infrared, and VLA radio data available to date, we have robustly identified a sample of 51 radio-excess AGN (~1300 deg^-2) out to redshift z~3. These radio-excess AGN have a significantly lower far-IR/radio ratio (q<1.68) than the typical relation observed for star-forming galaxies (q~2.2). We find that ~45% of these radio-excess sources have a dominant AGN component in the mid-IR band, while for the remainders the excess radio emission is the only indicator of AGN activity. The fraction of radio-excess AGN increases with X-ray luminosity reaching ~60% at Lx~10^44-10^45 erg/s, making these sources an important part of the total AGN population. However, almost half (24/51) of these radio-excess AGN are not detected in the deep Chandra X-ray data, suggesting that some of these sources might be heavily obscured. We also find that the specific star formation rates (sSFRs) of the radio-excess AGN are on average lower that those observed for X-ray selected AGN hosts, indicating that our sources are forming stars more slowly than typical AGN hosts, and possibly their star formation is progressively quenching.

We explore the effects of active galactic nuclei (AGN) and star formation activity on the infrared (0.3-1000 microns) spectral energy distributions of luminous infrared galaxies from z = 0.5 to 4.0. We have compiled a large sample of 151 galaxies selected at 24 microns (S24 > 100 uJy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the mid-IR spectrum into contributions from star formation and AGN activity. A significant portion (~25%) of our sample is dominated by an AGN in the mid-IR. Based on the mid-IR classification, we divide our full sample into four sub-samples: z~1 star-forming (SF) sources; z~2 SF sources; AGN with clear 9.7 micron silicate absorption; and AGN with featureless mid-IR spectra. From our large spectroscopic sample and wealth of multi-wavelength data, including deep Herschel imaging at 100, 160, 250, 350, and 500 microns, we use 95 galaxies with complete spectral coverage to create a composite spectral energy distribution (SED) for each sub-sample. We then fit a two-temperature component modified blackbody to the SEDs. We find that the IR SEDs have similar cold dust temperatures, regardless of the mid-IR power source, but display a marked difference in the warmer dust temperatures. We calculate the average effective temperature of the dust in each sub-sample and find a significant (~20 K) difference between the SF and AGN systems. We compare our composite SEDs to local templates and find that local templates do not accurately reproduce the mid-IR features and dust temperatures of our high redshift systems. High redshift IR luminous galaxies contain significantly more cool dust than their local counterparts. We find that a full suite of photometry spanning the IR peak is necessary to accurately account for the dominant dust temperature components in high redshift IR luminous galaxies.

In this paper, we characterize the infrared spectral energy distributions (SEDs) of mid-IR selected z~0.3-3.0 and L_IR~10^11-10^13Lsun galaxies, and study how their SEDs differ from those of local and high-z analogs. Our mid-IR flux-limited sample of 191 sources is unique in size, and spectral coverage, including Spitzer mid-IR spectroscopy. Here we add Herschel photometry at 250um, 350um, and 500um, which allows us to obtain accurate total IR luminosities, as well as constrain the relative contributions of AGN and starbursts to those luminosities. Our sample constitutes ~23% AGN (i.e. where the AGN contributes >50% of L_IR), ~30% starbursts (where AGN contributes <20% of L_IR and the mid-IR spectra are starburst-like); and ~47% composites (which show both significant AGN and starburst activity). The AGN-dominated sources divide into ones that show a strong silicate 9.7um absorption feature, implying highly obscured systems, and ones that do not. The high-tau_9.7 sources are half of our z>1.2 AGN, but show SEDs that are extremely rare among local AGN. The SEDs of our z~2 starburst-dominated ULIRGs are much closer to those of local LIRGs than ULIRGs. This is consistent with our earlier finding that, unlike local ULIRGs, our high-z starbursts are typically only in the early stages of a merger. The SEDs of the composite sources are most similar to the local archetypal warm ULIRG, Mrk231. In summary, our results show that there is strong evolution in the SEDs between local and z~2 IR-luminous galaxies, as well as that there is a wide range of SEDs among high redshift IR-luminous sources. The publicly-available SED templates we derive from our sample will be particularly useful for infrared population synthesis models, as well as in the interpretation of other mid-IR high-z galaxies in particular those detected by the recent all sky WISE survey.

We carried out CO(1-0) and CO(2-1) observations of 21 different regions in the vicinity of M86, NGC4438, and along the 120 kpc-long, Ha-emitting filamentary trail that connects them, aiming to test whether molecular gas can survive to be transferred from a spiral to an elliptical galaxy in Virgo's 10^7K intracluster medium (ICM). We targeted Ha-emitting regions that could be associated with the interface between cold molecular clouds and the hot ionized ICM. The data, obtained with the 30m telescope of the Institut de Radioastronomie Millimetrique, led to the detection of molecular gas close to M86. CO gas with a recession velocity that is similar to that of the stars, -265 km/s, and with a corresponding H2 mass of 2*10^7 M_sun, was detected ~10 kpc southeast of the nucleus of M86, near the peak of its HI emission. We argue that it is possible for this molecular gas either to have formed in situ from HI, or to have been stripped from NGC4438 directly in molecular form. In situ formation is nonetheless negligible for the 7*10^6 M_sun of gas detected at 12:26:15.9+12:58:49, at ~10 kpc northeast of M86, where no (strong) HI emission is present. This detection provides evidence for the survival of molecular gas in filaments for timescales of ~100 Myr. An amount equivalent to 5*10^7 M_sun of H2 gas that could be lost to the ICM or to neighboring galaxies was also discovered in the tidal tail northwest of NGC4438. A scenario of gas being alternatively brought to M86 from NGC4388 on its south was also examined but it was considered unlikely due to the non detection of CO below or at the HI stream velocities, 2000-2700 km/s.

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.

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.

We queried the Spitzer archive for high-resolution observations with the Infrared Spectrograph of optically selected active galactic nuclei (AGN) for the purpose of identifying sources with resolved fine-structure lines that would enable studies of the narrow-line region (NLR) at mid-infrared wavelengths. By combining 298 Spitzer spectra with 6 Infrared Space Observatory spectra, we present kinematic information of the NLR for 81 z<=0.3 AGN. We used the [NeV], [OIV], [NeIII], and [SIV] lines, whose fluxes correlate well with each other, to probe gas photoionized by the AGN. We found that the widths of the lines are, on average, increasing with the ionization potential of the species that emit them. No correlation of the line width with the critical density of the corresponding transition was found. The velocity dispersion of the gas, sigma, is systematically higher than that of the stars, sigma_*, in the AGN host galaxy, and it scales with the mass of the central black hole, M_BH. Further correlations between the line widths and luminosities L, and between L and M_BH, are suggestive of a three dimensional plane connecting log(M_BH) to a linear combination of log(sigma) and log(L). Such a plane can be understood within the context of gas motions that are driven by AGN feedback mechanisms, or virialized gas motions with a power-law dependence of the NLR radius on the AGN luminosity. The M_BH estimates obtained for 35 type 2 AGN from this plane are consistent with those obtained from the M_BH-sigma_* relation.

To optimally identify candidates of the Compton-thick (CT) active galactic nuclei (AGN) that contribute to the unresolved X-ray background in infrared surveys, a tracer of column density is desirable in addition to an AGN indicator. In a recent study, we aimed to test whether the 9.7 micron silicate absorption feature can be used for this purpose when seen at high optical depths. We found that the extreme criterion of optical thickness at 9.7 micron is efficient in identifying CT objects among local AGN. Having identified six of the nine CT AGN in the 12 micron sample with Spitzer and X-ray spectra, we expanded this analysis at intermediate/high z, using all GOODS and FLS sources with Spitzer and X-ray observations. We found 12 sources with tau_9.7>1 that host an AGN between 0.8<z<2.7. Four of them are likely to be CT according to their low X-ray to 6 micron luminosity ratio. Surveys with complete coverage in both mid-infrared spectra and X-ray data can provide large populations of such sources, as at least 5-9% of all infrared bright galaxies in the GOODS and FLS samples are tau_9.7>1 AGN.

Querying the Spitzer archive for optically-selected active galactic nuclei (AGN) observed in high-resolution mode spectroscopy, we identified radio and/or interacting galaxies with highly turbulent motions of the H2 gas at a temperature of a few hundred Kelvin. Unlike all other AGN that have unresolved H2 line profiles at a spectral resolution of ~600, 3C236, 3C293, IRAS09039+0503, MCG-2-58-22, and Mrk463E have intrinsic velocity dispersions exceeding 200 km/s for at least two of the rotational S0, S1, S2, and S3 lines. In a sixth source, 4C12.50, a blue wing was detected in the S1 and S2 line profiles, indicating the presence of a warm molecular gas component moving at -640 km/s with respect to the bulk of the gas at systemic velocity. Its mass is 5.2*10^7 M_sun, accounting for more than one fourth of the H2 gas at 374K, but less than 1% of the cold H2 gas computed from CO observations. Because no diffuse gas component of 4C12.50 has been observed to date to be moving at more than 250 km/s from systemic velocity, the H2 line wings are unlikely to be tracing gas in shock regions along the tidal tails of this merging system. They can instead be tracing gas driven by a jet or entrained by a nuclear outflow, which is known to emerge from the west nucleus of 4C12.50. It is improbable that such an outflow, with an estimated mass loss rate of 130 M_sun/yr, entirely quenches the star formation around this nucleus.

Star formation in massive galaxies is quenched at some point during hierarchical mass assembly. To understand where and when the quenching processes takes place, we study the evolution of the total star formation rate per unit total halo mass (\Sigma(SFR/M)) in three different mass scales: low mass halos (field galaxies), groups, and clusters, up to a redshift ~1.6. We use deep far-infrared PACS data at 100 and 160 um to accurately estimate the total star formation rate of the Luminous Infrared Galaxy population of 9 clusters with mass ~10^{15} M_{\odot}, and 9 groups/poor clusters with mass ~ 5 x 10^{13} M_{\odot}. Estimates of the field \Sigma(SFR/M) are derived from the literature, by dividing the star formation rate density by the mean comoving matter density of the universe. The field \Sigma(SFR/M) increases with redshift up to z~1 and it is constant thereafter. The evolution of the \Sigma(SFR/M)-z relation in galaxy systems is much faster than in the field. Up to redshift z~0.2, the field has a higher \Sigma(SFR/M) than galaxy groups and galaxy clusters. At higher redshifts, galaxy groups and the field have similar \Sigma(SFR/M), while massive clusters have significantly lower \Sigma(SFR/M) than both groups and the field. There is a hint of a reversal of the SFR activity vs. environment at z~1.6, where the group \Sigma(SFR/M) lies above the field \Sigma(SFR/M)-z relation. We discuss possible interpretations of our results in terms of the processes of downsizing, and star-formation quenching.

We present the deepest far-IR observations obtained with Herschel and examine the 3-500um SEDs of galaxies at 0<z<2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer and AKARI data. We find that the ratio of total IR luminosity to rest-frame 8um luminosity, IR8 (=Lir/L8), follows a Gaussian distribution centered on IR8=4 and defines an IR main sequence (MS). A minority population (<20 %) of outliers producing a tail skewed toward higher values of IR8 consist of starbursts (SB) with compact projected star formation densities. IR8 can be used to separate galaxies with normal and extended modes of star formation from compact SBs with high-IR8, high projected IR surface brightness (>3x10^10 Lsun kpc^-2) and a high specific SFR (i.e., SBs). The rest-frame, UV-2700A size of these distant SBs is typically half that of MS galaxies, supporting the correlation between star formation density and SB activity that is measured for the local sample. Locally, (U)LIRGs are systematically in the SB mode, whereas most distant (U)LIRGs form stars in the "normal" MS mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z>1.5 by previous studies. MS galaxies have strong PAH emission line features, a broad far-IR bump resulting from a combination of dust temperatures (Tdust~15-50 K), and an effective Tdust~31 K, as derived from the peak wavelength of their IR SED. Galaxies in the SB regime instead exhibit weak PAH EW and a sharper far-IR bump with an effective Tdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray AGNs is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty SBs. After correcting for the effect of SBs on IR8, we identify new candidates for extremely obscured AGNs.

We study the impact of galaxy-galaxy interactions on the FIR properties of galaxies and its evolution at 0<z<1.2. Using the high-z galaxies in the fields of GOODS observed by Herschel in the framework of the GOODS-Herschel key program and the local IRAS or AKARI-selected galaxies in the field of SDSS DR7, we investigate the dependence of galaxy properties on the morphology of and the distance to the nearest neighbor galaxy. We find that the SFRs and the SSFRs of galaxies, on average, depend on the morphology of and the distance to the nearest neighbor galaxy in this redshift range. When a late-type galaxy has a close neighbor galaxy, the SFR and the SSFR increase as it approaches a late-type neighbor, which is supported by K-S and MC tests with a significance level of >99%. However, the SFR and the SSFR decrease or do not change much as it approaches an early-type neighbor. The bifurcations of SFRs and SSFRs depending on the neighbor's morphology seem to occur at R_n~0.5r_vir,nei, which is supported by K-S and MC tests with a significance level of >98%. For all redshift bins, the SSFRs of late-type galaxies interacting with late-type neighbors are increased by factors of about 1.8+/-0.7 and 4.0+/-1.2 compared to those of non-interacting galaxies when the pair separation is smaller than 0.5r_vir,nei and 0.1r_vir,nei, respectively. The Tdust of both local and high-z late-type galaxies strongly interacting with late-type neighbors (i.e. R_n < 0.1r_vir,nei) appears to be higher than that of non-interacting galaxies with a significance level of 96-99%. However, Tdust of local late-type galaxies strongly interacting with early-type neighbors seems to be lower than or similar to that of non-interacting galaxies. Our results suggest that galaxy-galaxy interactions and mergers have been strongly affecting the SFR and the dust properties of star-forming galaxies over at least 8 billion years.

We present a study of the infrared properties of X-ray selected, moderate luminosity (Lx=10^{42}-10^{44}ergs/s) active galactic nuclei (AGNs) up to z~3, to explore the links between star formation in galaxies and accretion onto their central black holes. We use 100um and 160um fluxes from GOODS-Herschel -the deepest survey yet undertaken by the Herschel telescope- and show that in >94 per cent of cases these fluxes are dominated by the host. We find no evidence of any correlation between the X-ray and infrared luminosities of moderate AGNs at any redshift, suggesting that star-formation is decoupled from nuclear (AGN) activity. The star formation rates of AGN hosts increase strongly with redshift; by a factor of 43 from z<0.1 to z=2-3 for AGNs with the same X-ray luminosities. This increase is consistent with the factor of 25-50 increase in the specific star formation rates (SSFRs) of normal, star-forming (main-sequence) galaxies. Indeed, the average SSFRs of AGN hosts are only marginally (20 per cent) lower than those of main-sequence galaxies, with this small deficit being due to a fraction of AGNs residing in quiescent (low-SSFR) galaxies. We estimate 79+/-10 per cent of moderate AGNs are hosted in main-sequence galaxies, 15+/-7 per cent in quiescent galaxies and <10 per cent in strongly starbursting galaxies. The fractions of all main sequence galaxies at z<2 experiencing a period of moderate nuclear activity is strongly dependent on galaxy stellar mass (Mstars); rising from a few per cent at Mstars~10^{10}Msun to >20 per cent at Mstars>10^{11}Msun. Our results indicate that it is galaxy stellar mass that is most important in dictating whether a galaxy hosts a moderate luminosity AGN. We argue that the majority of moderate nuclear activity is fuelled by internal mechanisms rather than violent mergers, suggesting that disk instabilities could be an important AGN feeding mechanism.

Using extremely deep PACS 100- and 160um Herschel data from the GOODS-Herschel program, we identify 21 infrared bright galaxies previously missed in the deepest 24um surveys performed by MIPS. These MIPS dropouts are predominantly found in two redshift bins, centred at z ~0.4 and ~1.3. Their S_100/S_24 flux density ratios are similar to those of local LIRGs and ULIRGs, whose silicate absorption features at 18um (at z ~ 0.4) and 9.7um (at z ~ 1.3) are shifted into the 24um MIPS band at these redshifts. The high-z sub-sample consists of 11 infrared luminous sources, accounting for ~2% of the whole GOODS-Herschel sample and putting strong upper limits on the fraction of LIRGs/ULIRGs at 1.0<z<1.7 that are missed by the 24um surveys. We find that a S_100/S_24 > 43 colour cut selects galaxies with a redshift distribution similar to that of the MIPS dropouts and when combined with a second colour cut, S_16/S_8 > 4, isolates sources at 1.0 < z < 1.7. We show that these sources have elevated specific star formation rates (sSFR) compared to main sequence galaxies at these redshifts and are likely to be compact starbursts with moderate/strong 9.7um silicate absorption features in their mid-IR spectra. Herschel data reveal that their infrared luminosities extrapolated from the 24um flux density are underestimated, on average, by a factor of ~3. These silicate break galaxies account for 16% (8%) of the ULIRG (LIRG) population in the GOODS fields, indicating a lower limit in their space density of 2.0 \times 10^(-5) Mpc^(-3). Finally, we provide estimates of the fraction of z < 2 MIPS dropout sources as a function of the 24-, 100-, 160-, 250- and 350um sensitivity limits, and conclude that previous predictions of a population of silicate break galaxies missed by the major 24um extragalactic surveys have been overestimated.

We take advantage of the sensitivity and resolution of Herschel at 100 and 160 micron to directly image the thermal dust emission and investigate the infrared luminosities, L(IR), and dust obscuration of typical star-forming (L*) galaxies at high redshift. Our sample consists of 146 UV-selected galaxies with spectroscopic redshifts 1.5<z<2.6 in the GOODS-North field. Supplemented with deep Very Large Array (VLA) and Spitzer imaging, we construct median stacks at the positions of these galaxies at 24, 100, and 160 micron, and 1.4 GHz. The comparison between these stacked fluxes and a variety of dust templates and calibrations implies that typical star-forming galaxies with UV luminosities L(UV)>1e10 Lsun at z~2 are luminous infrared galaxies (LIRGs) with a median L(IR)=(2.2+/-0.3)e11 Lsun. Typical galaxies at 1.5<z<2.6 have a median dust obscuration L(IR)/L(UV) = 7.1+/-1.1, which corresponds to a dust correction factor, required to recover the bolometric star formation rate (SFR) from the unobscured UV SFR, of 5.2+/-0.6. This result is similar to that inferred from previous investigations of the UV, H-alpha, 24 micron, radio, and X-ray properties of the same galaxies studied here. Stacking in bins of UV slope implies that L* galaxies with redder spectral slopes are also dustier, and that the correlation between UV slope and dustiness is similar to that found for local starburst galaxies. Hence, the rest-frame 30 and 50 micron fluxes validate on average the use of the local UV attenuation curve to recover the dust attenuation of typical star-forming galaxies at high redshift. In the simplest interpretation, the agreement between the local and high redshift UV attenuation curves suggests a similarity in the dust production and stellar and dust geometries of starburst galaxies over the last 10 billion years.

Dust attenuation curves in external galaxies are useful to study their dust properties as well as to interpret their intrinsic spectral energy distributions. In particular the presence or absence of a UV bump at 2175 A remains an open issue which has consequences on the interpretation of broad band colours of distant galaxies. We study the dust attenuation curve in the UV range at z >1. In particular we search for the presence of a UV bump. We use deep photometric data of the CDFS obtained with intermediate and broad band filters by the MUSYC project to sample the UV rest-frame of galaxies with 1<z <2. Herschel/PACS and Spitzer/MIPS data are used to measure the dust emission. 30 galaxies were selected with high S/N in all bands. Their SEDs from the UV to the far-IR are fitted using the CIGALE code and the characteristics of the dust attenuation curve are obtained. The mean dust attenuation curve we derive exhibits a UV bump at 2175A whose amplitude corresponds to 35 % (76%) that of the MW (LMC2 supershell) extinction curve. An analytical expression of the average attenuation curve is given, it is found slightly steeper than the Calzetti et al. one, although at a 1 sigma level. Our galaxy sample is used to study the derivation of the slopes of the UV continuum from broad band colours, including the GALEX FUV-NUV colour. Systematic errors induced by the presence of the bump are quantified. We compare dust attenuation factors measured with CIGALE to the slope of the UV continuum and find that there is a large scatter around the relation valid for local starbursts (0.7 mag). The uncertainties on the determination of the UV slope lead to an extra systematic error of the order of 0.3 to 0.7 mag on dust attenuation when a filter overlaps the UV bump.

Spitzer/IRS has revealed many sources with very deep Si features at 9.7micron (tau>1). We set out to investigate whether a strong Si absorption feature is a good indicator for the presence of a heavily obscured AGN. We compile X-ray spectroscopic observations available in the literature on the optically-thick,tau(9.7)>1 sources from the IRAS Seyfert sample. We find that the majority of the high-tau optically confirmed Seyferts (6/9) in this sample are probably CT. Thus we provide direct evidence for a connection between mid-IR optically-thick galaxies and CT AGN, with the success rate being close to 70% in the local Universe. This is at least comparable, if not better, than other rates obtained with photometric information in the mid to far-IR, or even mid-IR to Xray. However, this technique cannot provide complete CT AGN samples,ie there are many CT AGN which do not show significant Si absorption, with the most notable example being N1068. Having assessed the validity of the high 9.7micron technique locally, we attempt to construct a sample of candidate CT AGN at higher redshifts. We compile a sample of 7 high-tau sources in the GOODS and 5 in the Spitzer FLS. All these have been selected to have no PAH features EW(6.2)<0.3 in order to maximize the probability that they are AGN. 6 out of 7 sources in the GOODS have been detected in X-rays, while for the five FLS sources only X-ray flux upper limits are available. The high X-ray luminosities of the detected GOODS sources corroborates that these are AGN. For FLS, ancillary optical spectroscopy reveals hidden nuclei in two more sources. SED fitting can support the presence of an AGN in the vast majority of sources. We cannot derive useful X-ray spectroscopy constraints on whether these are CT. However, the low LX/L6 ratios, suggest that at least 4 out of the 6 detected sources in GOODS may be associated with CT AGN.

We present rest-frame 15 and 24 um luminosity functions and the corresponding star-forming luminosity functions at z<0.3 derived from the 5MUSES sample. Spectroscopic redshifts have been obtained for ~98% of the objects and the median redshift is ~0.12. The 5-35 um IRS spectra allow us to estimate accurately the luminosities and build the luminosity functions. Using a combination of starburst and quasar templates, we quantify the star-formation and AGN contributions in the mid-IR SED. We then compute the star-formation luminosity functions at 15 um and 24 um, and compare with the total 15 um and 24 um luminosity functions. When we remove the contribution of AGN, the bright end of the luminosity function exhibits a strong decline, consistent with the exponential cutoff of a Schechter function. Integrating the differential luminosity function, we find that the fractional contribution by star formation to the energy density is 58% at 15 um and 78% at 24 um, while it goes up to ~86% when we extrapolate our mid-IR results to the total IR luminosity density. We confirm that the active galactic nuclei play more important roles energetically at high luminosities. Finally, we compare our results with work at z~0.7 and confirm that evolution on both luminosity and density is required to explain the difference in the LFs at different redshifts.

We present new results on the physical nature of infrared-luminous sources at 0.5<z<2.8 as revealed by HST/NICMOS imaging and IRS mid-infrared spectroscopy. Our sample consists of 134 galaxies selected at 24\mum with a flux of S(24\mum) > 0.9 mJy. We find many (~60%) of our sources to possess an important bulge and/or central point source component, most of which reveal additional underlying structures after subtraction of a best-fit sersic (or sersic+PSF) profile. Based on visual inspection of the NIC2 images and their residuals, we estimate that ~80% of all our sources are mergers. We calculate lower and upper limits on the merger fraction to be 62% and 91% respectively. At z < 1.5, we observe objects in early (pre-coalescence) merging stages to be mostly disk and star formation dominated, while we find mergers to be mainly bulge-dominated and AGN-starburst composites during coalescence and then AGN-dominated in late stages. This is analogous to what is observed in local ULIRGs. At z \geq 1.5, we find a dramatic rise in the number of objects in pre-coalescence phases of merging, despite an increase in the preponderance of AGN signatures in their mid-IR spectra and luminosities above 10^12.5 L_sun. We further find the majority of mergers at those redshifts to retain a disk-dominated profile during coalescence. We conclude that, albeit still driven by mergers, these high-z ULIRGs are substantially different in nature from their local counterparts and speculate that this is likely due to their higher gas content. Finally, we observe obscured ({\tau}_{9.7\mum} > 3.36) quasars to live in faint and compact hosts and show that these are likely high-redshift analogs of local dense-core mergers. We find late-stage mergers to show predominantly unobscured AGN spectra, but do not observe other morphological classes to occupy any one specific region in the Spoon diagram. [abridged]

We study a 24\,$\mu$m selected sample of 330 galaxies observed with the Infrared Spectrograph for the 5\,mJy Unbiased Spitzer Extragalactic Survey. We estimate accurate total infrared luminosities by combining mid-IR spectroscopy and mid-to-far infrared photometry, and by utilizing new empirical spectral templates from {\em Spitzer} data. The infrared luminosities of this sample range mostly from 10$^9$L$_\odot$ to $10^{13.5}$L$_\odot$, with 83% in the range 10$^{10}$L$_\odot$$<$L$_{\rm IR}$$<10^{12}$L$_\odot$. The redshifts range from 0.008 to 4.27, with a median of 0.144. The equivalent widths of the 6.2\,$\mu$m aromatic feature have a bimodal distribution. We use the 6.2\,$\mu$m PAH EW to classify our objects as SB-dominated (44%), SB-AGN composite (22%), and AGN-dominated (34%). The high EW objects (SB-dominated) tend to have steeper mid-IR to far-IR spectral slopes and lower L$_{\rm IR}$ and redshifts. The low EW objects (AGN-dominated) tend to have less steep spectral slopes and higher L$_{\rm IR}$ and redshifts. This dichotomy leads to a gross correlation between EW and slope, which does not hold within either group. AGN dominated sources tend to have lower log(L$_{\rm PAH 7.7\mu m}$/L$_{\rm PAH 11.3\mu m}$) ratios than star-forming galaxies, possibly due to preferential destruction of the smaller aromatics by the AGN. The log(L$_{\rm PAH 7.7\mu m}$/L$_{\rm PAH 11.3\mu m}$) ratios for star-forming galaxies are lower in our sample than the ratios measured from the nuclear spectra of nearby normal galaxies, most probably indicating a difference in the ionization state or grain size distribution between the nuclear regions and the entire galaxy. Finally, we provide a calibration relating the monochromatic 5.8, 8, 14 and 24um continuum or Aromatic Feature luminosity to L$_{\rm IR}$ for different types of objects.

The Ultra Luminous InfraRed Galaxy Mrk 231 reveals up to seven rotational lines of water (H2O) in emission, including a very high-lying (E_{upper}=640 K) line detected at a 4sigma level, within the Herschel/SPIRE wavelength range, whereas PACS observations show one H2O line at 78 microns in absorption, as found for other H2O lines previously detected by ISO. The absorption/emission dichotomy is caused by the pumping of the rotational levels by far-infrared radiation emitted by dust, and subsequent relaxation through lines at longer wavelengths, which allows us to estimate both the column density of H2O and the general characteristics of the underlying far-infrared continuum source. Radiative transfer models including excitation through both absorption of far-infrared radiation emitted by dust and collisions are used to calculate the equilibrium level populations of H2O and the corresponding line fluxes. The highest-lying H2O lines detected in emission, with levels at 300-640 K above the ground state, indicate that the source of far-infrared radiation responsible for the pumping is compact (radius=110-180 pc) and warm (T_{dust}=85-95 K), accounting for at least 45% of the bolometric luminosity. The high column density, N(H2O)~5x10^{17} cm^{-2}, found in this nuclear component, is most probably the consequence of shocks/cosmic rays, an XDR chemistry, and/or an "undepleted chemistry" where grain mantles are evaporated. A more extended region, presumably the inner region of the 1-kpc disk observed in other molecular species, could contribute to the flux observed in low-lying H2O lines through dense hot cores, and/or shocks. The H2O 78 micron line observed with PACS shows hints of a blue-shifted wing seen in absorption, possibly indicating the occurrence of H2O in the prominent outflow detected in OH (Fischer et al., this volume).

We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk231. In total 25 lines are detected, including CO J=5-4 through J=13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J=8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J=8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J=8. X-rays from the accreting supermassive black hole in Mrk231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.

We present CO observations of 9 ULIRGs at z~2 with S(24\mu m)>1mJy, previously confirmed with the mid-IR spectra in the Spitzer First Look Survey. All targets are required to have accurate redshifts from Keck/GEMINI near-IR spectra. Using the Plateau de Bure millimeter-wave Interferometer (PdBI) at IRAM, we detect CO J(3-2) [7 objects] or J(2-1) [1 object] line emission from 8 sources with integrated intensities Ic ~(5-9)sigma. The CO detected sources have a variety of mid-IR spectra, including strong PAH, deep silicate absorption and power-law continuum, implying that these molecular gas rich objects at z~2 could be either starbursts or dust obscured AGNs. The measured line luminosity L'[CO] is (1.28-3.77)e+10[K km/s pc^2]. The averaged molecular gas mass M(H2) is 1.7e+10Msun, assuming CO-to-H2 conversion factor of 0.8Msun/[K km/s pc^2]. Three sources (33%) -- MIPS506, MIPS16144 & MIPS8342 -- have double peak velocity profiles. The CO double peaks in MIPS506 and MIPS16144 show spatial separations of 45kpc and 10.9kpc, allowing the estimates of the dynamical masses of 3.2e+11*sin^(-2)(i)Msun and 5.4e+11*sin^{-2}(i)Msun respectively. The implied gas fraction, M(gas)/M(dyn), is 3% and 4%, assuming an average inclination angle. Finally, the analysis of the HST/NIC2 images, mid-IR spectra and IR SED revealed that most of our sources are mergers, containing dust obscured AGNs dominating the luminosities at (3-6)um. Together, these results provide some evidence suggesting SMGs, bright 24um z~2 ULIRGs and QSOs could represent three different stages of a single evolutionary sequence, however, a complete physical model would require much more data, especially high spatial resolution spectroscopy.

We present results from infrared spectroscopic projects that aim to test the relation between the mass of a black hole, M_BH, and the velocity dispersion of the stars in its host-galaxy bulge. We demonstrate that near-infrared, high-resolution spectroscopy assisted by adaptive optics is key in populating the high-luminosity end of the relation. We show that the velocity dispersions of mid-infrared, high-ionization lines originating from gas in the narrow-line region of the active galactic nucleus follow the same relation. This result provides a way of inferring M_BH estimates for the cosmologically significant population of obscured, type-2 AGN that can be applicable to data from spectrographs on the next generation infrared telescopes.

We present a catalog of mid-infrared (MIR) spectra of 150 infrared (IR) luminous galaxies in the Spitzer extragalactic first look survey obtained with IRS on board Spitzer. The sample is selected to be brighter than ~0.9 mJy at 24 micron and it has a z distribution in the range [0.3,3.5] with a peak at z=1. It primarily comprises ultraluminous IR galaxies at z>1 and luminous IR galaxies at z<1, as estimated from their monochromatic 14 micron luminosities. The number of sources with spectra that are dominated by an active galactic nucleus (AGN) continuum is 49, while 39 sources have strong, star-formation related features. For this classification, we used the equivalent width (EW) of the 11.3 micron polycyclic aromatic hydrocarbon (PAH) feature. Several intermediate/high z starbursts have higher PAH EW than local ULIRGs. An increase in the AGN activity is observed with increasing z and luminosity, based on the decreasing EW of PAHs and the increasing [NeIII]/[NeII] ratio. Spectral stacking leads to the detection of the 3.3 micron PAH, the H2 0-0 S(1) and S(3) lines, and the [NeV] line. We observe differences in the flux ratios of PAHs in the stacked spectra of IR-luminous galaxies with z or luminosity, which are not due to extinction effects. When placing the observed galaxies on IR color-color diagrams, we find that the wedge defining AGN comprises most sources with continuum-dominated spectra, but also contains many starbursts. The comparison of the 11.3 micron PAH EW and the H-band effective radius, measured from HST data, indicates that sources with EW>2 micron, are typically more extended than ~3 kpc. However, there is no strong correlation between the MIR spectral type and the near-IR extent of the sources. [Abridged].

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)

We present the first use of the Gemini North laser guide star adaptive optics (LGS AO) system and an integral field unit (IFU) to measure the stellar velocity dispersion of the host of a luminous quasar. The quasar PG1426+015 (z=0.086) was observed with the Near-Infrared Integral Field Spectrometer (NIFS) on the 8m Gemini North telescope in the H-band as part of the Science Verification phase of the new ALTAIR LGS AO system. The NIFS IFU and LGS AO are well suited for host studies of luminous quasars because one can achieve a large ratio of host to quasar light. We have measured the stellar velocity dispersion of PG1426+015 from 0.1'' to 1'' (0.16 kpc to 1.6 kpc) to be 217+/-15 km/s based on high signal-to-noise ratio measurements of Si I, Mg I, and several CO bandheads. This new measurement is a factor of four more precise than a previous measurement obtained with long-slit spectroscopy and good, natural seeing, yet was obtained with a shorter net integration time. We find that PG1426+015 has a velocity dispersion that places it significantly above the M-sigma relation of quiescent galaxies and lower-luminosity active galactic nuclei with black hole masses estimated from reverberation mapping. We discuss several possible explanations for this discrepancy that could be addressed with similar observations of a larger sample of luminous quasars.

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.

We present the far-IR, millimeter, and radio photometry as well as optical and near-IR spectroscopy of a sample of 48 z~1-3 Spitzer-selected ULIRGs with IRS mid-IR spectra. Our goals are to compute their bolometric emission, and to determine both the presence and relative strength of their AGN and starburst components. We find that strong-PAH sources tend to have higher 160um and 1.2mm fluxes than weak-PAH sources. The depth of the 9.7um silicate feature does not affect MAMBO detectability. We fit the far-IR SEDs of our sample and find an average <L_{IR}>~7x10^{12}Lsun for our z>1.5 sources. Spectral decomposition suggests that strong-PAH sources typically have ~20-30% AGN fractions. Weak-PAH sources by contrast tend to have >~70% AGN fractions, with a few sources having comparable contributions of AGN and starbursts. The optical line diagnostics support the presence of AGN in the bulk of the weak-PAH sources. With one exception, our sources are narrow-line sources, show no obvious correspondence between the optical extinction and the silicate feature depth, and, in two cases, show evidence for outflows. Radio AGN are present in both strong-PAH and weak-PAH sources. This is supported by our sample's far-IR-to-radio ratios (q) being consistently below the average value of 2.34 for local star-forming galaxies. We use survival analysis to include the lower-limits given by the radio-undetected sources, arriving at <q>=2.07+/-0.01 for our z>1.5 sample. In total, radio and, where available, optical line diagnostics support the presence of AGN in 57% of the z>1.5 sources, independent of IR-based diagnostics. For higher-z sources, the AGN luminosities alone are estimated to be >10^{12}Lsun, which, supported by the [OIII] luminosities, implies that the bulk of our sources host obscured quasars.

We present relations of the black hole mass and the optical luminosity with the velocity dispersion and the luminosity of the [Ne V] and the [O IV] high-ionization lines in the mid-infrared (MIR) for 28 reverberation-mapped active galactic nuclei. We used high-resolution Spitzer Infrared Spectrograph and Infrared Space Observatory Short Wavelength Spectrometer data to fit the profiles of these MIR emission lines that originate from the narrow-line region of the nucleus. We find that the lines are often resolved and that the velocity dispersion of [Ne V] and [O IV] follows a relation similar to that between the black hole mass and the bulge stellar velocity dispersion found for local galaxies. The luminosity of the [Ne V] and the [O IV] lines in these sources is correlated with that of the optical 5100A continuum and with the black hole mass. Our results provide a means to derive black hole properties in various types of active galactic nuclei, including highly obscured systems.

This is the second paper studying the QSOs in the spitzer QUEST sample. Previously we presented new PAH measurements and argued that most of the observed far infrared (FIR) radiation is due to star-forming activity. Here we present spectral energy distributions (SEDs) by supplementing our data with optical, NIR and FIR observations. We define two sub-groups of ``weak FIR'' and ``strong FIR'' QSOs, and a third group of FIR non-detections. Assuming a starburst origin for the FIR, we obtain ``intrinsic'' AGN SEDs by subtracting a starburst template from the mean SEDs. The resulting SEDs are remarkably similar for all groups. They show three distinct peaks corresponding to two silicate emission features and a 3mic bump that we interpret as the signature of the hottest AGN dust. They also display drops beyond 20mic that we interpret as the signature of the minimum temperature (about 200K) dust. This component must be optically thin to explain the silicate emission and the slope of the long wavelength continuum. We discuss the merits of an alternative model where most of the FIR emission is due to AGN heating. Such models are unlikely to explain the properties of our QSOs but they cannot be ruled out for more luminous objects. We also find correlations between the luminosity at 5100A and two infrared starburst indicators: L(60mic) and L(PAH 7.7mic). The correlation of L(5100A) with L(60mic) can be used to measure the relative growth rates and lifetimes of the black hole and the new stars.

We present host-galaxy velocity dispersions of 12 local (mainly Palomar-Green) quasi-stellar objects (QSOs) measured directly from the stellar CO absorption features in the H band. The mean bulge velocity dispersion of the QSOs in our sample is 186 km/s with a standard deviation of 24 km/s. The measurement of the stellar velocity dispersion in QSOs enables us to place them on observational diagrams such as the local black-hole mass to bulge-velocity-dispersion relation and the fundamental plane of early-type galaxies. Concerning the former relation, these QSOs have higher black hole masses than most Seyfert 1 AGNs with similar velocity dispersions. On the fundamental plane, PG QSOs are located between the regions occupied by moderate-mass and giant ellipticals. The QSO bulge and black hole masses, computed from the stellar velocity dispersions, are of order 10^11 M_sun and 10^8 M_sun respectively. The Eddington efficiency of their black holes is on average 0.25, assuming that all of the bolometric luminosity originates from the active nucleus. Our data are consistent with other lines of evidence that Palomar-Green QSOs are related to galaxy mergers with gas-rich components and that they are formed in a manner similar to the most massive Ultraluminous Infrared Galaxies, regardless of their far-infrared emission. However, PG QSOs seem to have smaller host dispersions and different formation mechanisms than QSOs with supermassive black holes of 5x10^8-10^9 M_sun that accrete at low rates and reside in massive spheroids.

We present results from our Very Large Telescope large program to study the dynamical evolution of local Ultraluminous Infrared Galaxies (ULIRGs) and QSOs. This paper is the second in a series presenting the stellar kinematics of 54 ULIRGs, derived from high resolution, long-slit H- and K-band spectroscopy. The data presented here, including observations of 17 new targets, are mainly focused on sources that have coalesced into a single nucleus. The stellar kinematics, extracted from the CO ro-vibrational bandheads in our spectra, indicate that ULIRG remnants are dynamically heated systems with a mean dispersion of 161 km/s. The combination of kinematic, structural, and photometric properties of the remnants indicate that they mostly originate from major mergers and that they result in the formation of systems supported by random motions, therefore, elliptical galaxies. The peak of the velocity dispersion distribution and the locus of ULIRGs on the fundamental plane of early-type galaxies indicate that the end products of ultraluminous mergers are typically moderate-mass ellipticals (of stellar mass ~10^10 - 10^11 M_sun). Converting the host dispersion into black hole mass with the aid of the M_BH-sigma relation yields black hole mass estimates of the order 10^7 - 10^8 M_sun and high accretion rates with Eddington efficiencies often >0.5.

This paper addresses the origin of the far-infrared (FIR) continuum of QSOs, based on the Quasar and ULIRG Evolution Study (QUEST) of nearby QSOs and ULIRGs using observations with the Spitzer Space Telescope. For 27 Palomar-Green QSOs at z <~ 0.3, we derive luminosities of diagnostic lines ([NeII]12.8um, [NeV]14.3um, [OIV]25.9um) and emission features (PAH7.7um emission which is related to star formation), as well as continuum luminosities over a range of mid- to far-infrared wavelengths between 6 and 60um. We detect star-formation related PAH emission in 11/26 QSOs and fine-structure line emission in all of them, often in multiple lines. The detection of PAHs in the average spectrum of sources which lack individual PAH detections provides further evidence for the widespread presence of PAHs in QSOs. Similar PAH/FIR and [NeII]/FIR ratios are found in QSOs and in starburst-dominated ULIRGs and lower luminosity starbursts. We conclude that the typical QSO in our sample has at least 30% but likely most of the far-infrared luminosity (~ 10^(10...12)Lsun) arising from star formation, with a tendency for larger star formation contribution at the largest FIR luminosities. In the QSO sample, we find correlations between most of the quantities studied including combinations of AGN tracers and starburst tracers. The common scaling of AGN and starburst luminosities (and fluxes) is evidence for a starburst-AGN connection in luminous AGN. Strong correlations of far-infrared continuum and starburst related quantities (PAH, low excitation [NeII]) offer additional support for the starburst origin of far-infrared emission.

We report the results from a deep HST NICMOS H-band imaging survey of a carefully selected sample of 33 luminous, late-stage galactic mergers at z < 0.3. Signs of a recent galactic interaction are seen in all of the objects in the HST sample, including all 7 IR-excess Palomar-Green (PG) QSOs in the sample. Unsuspected double nuclei are detected in 5 ULIRGs. A detailed two-dimensional analysis of the surface brightness distributions in these objects indicates that the great majority (81%) of the single-nucleus systems show a prominent early-type morphology. However, low-surface-brightness exponential disks are detected on large scale in at least 4 of these sources. The hosts of 'warm' AGN-like systems are of early type and have less pronounced merger-induced morphological anomalies than the hosts of cool systems with LINER or HII region-like nuclear optical spectral types. The host sizes and luminosities of the 7 PG~QSOs in our sample are statistically indistinguishable from those of the ULIRG hosts. In comparison, highly luminous quasars, such as those studied by Dunlop et al. (2003), have hosts which are larger and more luminous. The hosts of ULIRGs and PG QSOs lie close to the locations of intermediate-size (about 1 -- 2 L*) spheroids in the photometric projection of the fundamental plane of ellipticals, although there is a tendency in our sample for the ULIRGs with small hosts to be brighter than normal spheroids. Excess emission from a young stellar population in the ULIRG/QSO hosts may be at the origin of this difference. Our results provide support for a possible merger-driven evolutionary connection between cool ULIRGs, warm ULIRGs, and PG~QSOs although this sequence may break down at low luminosity. (abridged)

We present a study of the dynamical evolution of Ultraluminous Infrared Galaxies (ULIRGs), merging galaxies of infrared luminosity >10^12 L_sun. During our Very Large Telescope large program, we have obtained ISAAC near-infrared, high-resolution spectra of 54 ULIRGs (at several merger phases) and 12 local Palomar-Green QSOs to investigate whether ULIRGs go through a QSO phase during their evolution. One possible evolutionary scenario is that after nuclear coalescence, the black hole radiates close to Eddington to produce QSO luminosities. The mean stellar velocity dispersion that we measure from our spectra is similar (~160 km/s) for 30 post-coalescence ULIRGs and 7 IR-bright QSOs. The black holes in both populations have masses of order 10^7-10^8 M_sun (calculated from the relation to the host dispersion) and accrete at rates >0.5 Eddington. Placing ULIRGs and IR-bright QSOs on the fundamental plane of early-type galaxies shows that they are located on a similar region (that of moderate-mass ellipticals), in contrast to giant ellipticals and radio-loud QSOs. While this preliminary comparison of the ULIRG and QSO host kinematical properties indicates that (some) ULIRGs may undergo a QSO phase in their evolutionary history before they settle down as ellipticals, further data on non-IR excess QSOs are necessary to test this scenario.

We present first results from our Very Large Telescope large program to study the dynamical evolution of Ultraluminous Infrared Galaxies (ULIRGs), which are the products of mergers of gas-rich galaxies. The full data set consists of high resolution, long-slit, H- and K-band spectra of 38 ULIRGs and 12 QSOs (between 0.042<z<0.268). In this paper, we present the sources that have not fully coalesced, and therefore have two distinct nuclei. This sub-sample consists of 21 ULIRGs, the nuclear separation of which varies between 1.6 and 23.3 kpc. From the CO bandheads that appear in our spectra, we extract the stellar velocity dispersion, sigma, and the rotational velocity, V_rot. The stellar dispersion equals 142 km/s on average, while V_rot is often of the same order. We combine our spectroscopic results with high-resolution infrared (IR) imaging data to study the conditions for ULIRG activity in interacting pairs. We find that the majority of ULIRGs are triggered by almost equal-mass major mergers of 1.5:1 average ratio. Less frequently, 3:1 encounters are also observed in our sample. However, less violent mergers of mass ratio >3:1 typically do not force enough gas into the center to generate ULIRG luminosities.

We present results from our VLT Large Program to study the dynamical evolution of Ultraluminous Infrared Galaxies (ULIRGs). We have so far obtained near-infrared high-resolution ISAAC spectra of 53 local ULIRGs at several merger timescales and 12 Palomar-Green QSOs (more than half of which are IR-bright sources). We have extracted the stellar velocity dispersion and rotational velocity along our slits to derive the kinematics of the merging galaxies. These quantities enable us to answer the following questions about the evolution of ULIRGs: 1) What are the progenitor mass ratios?, 2) How do the stellar kinematics evolve?, and, 3) Is there a connection between ULIRGs and QSOs? We find that the Ultraluminous phase is mainly triggered by mergers of approximately equal mass galaxies, however, less violent minor mergers (of progenitor mass ratio ~3:1) also exist in our sample. Dynamical heating of the merging hosts is observed as the stellar systematic rotation decreases with time in favour of the increase of random motions. The merger remnants, being dispersion-dominated systems with non-negligible rotation, resemble elliptical galaxies. Placing ULIRGs on the fundamental plane of early-type galaxies shows that they resemble intermediate-mass ellipticals. After the nuclear coalescence, the black hole masses of ULIRGs, calculated from their relation to the host dispersions, are of the order 10^7-10^8 M_sun. To investigate whether ULIRGs go through a QSO phase during their evolution, we perform a similar (preliminary) analysis of the IR-bright-QSO kinematics. We find that the average dispersion of the IR-bright QSOs is similar to that of the ULIRG remnants, indicating that evolutionary links between the two populations may exist.

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]

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.

We use the radiation transfer simulation of Xilouris et al. (1999) to constrain the quantity of dust in three nearby spiral galaxies (NGC 4013, NGC 5907 and NGC 4565). The predicted visual optical depth from the model is compared with the thermal continuum radiation detected from NGC 4013 and NGC 5907 at 850 microns and from NGC 4565 at 1.2 mm. The former is based on SCUBA images of NGC 4013 and NGC 5907, reduced and presented for the first time in this work. The comparison of visual optical depth and 850 microns (1.2 mm) emission yields the emissivity of dust grains in the submillimeter (millimeter) waveband. We infer a value of 1.2x10^{4} for the emissivity at 850 microns which is a factor 4 higher than the benchmark, semi-empirical model of Draine & Lee (1984). At 1.2 mm our values are a factor 1.5 higher than this model. Our estimates are more closely aligned with recent measurements carried out in the laboratory on amorphous carbon and silicate particulates. A comparison between the distribution of 850 microns (1.2 mm) surface brightness and the intensity levels in the ^{12}CO(1-0) and 21 cm lines underlines the spatial association between dust detected in the submillimeter/millimeter waveband and molecular gas clouds. We suggest that the relatively high emissivity values that we derive may be attributable to amorphous, fluffy grains situated in denser gas environments.