Galactic outflows are known to consist of several gas phases, however, so far the connection between these multiple phases has been investigated little and only in a few objects. In this paper, we analyse MUSE/VLT data of 26 local (U)LIRGs and study their ionised and neutral atomic phases. We also include objects from the literature to obtain a total sample of 31 galaxies with spatially resolved multi-phase outflow information. We find that the ionized phase of the outflows has on average an electron density three times higher than the disc ($n_{\rm e, disc}$ $\sim$ 150 cm$^{-3}$ vs $n_{\rm e, outflow}$ $\sim$ 500 cm$^{-3}$), suggesting that cloud compression in the outflow is more important that cloud dissipation. We find that the difference in extinction between outflow and disc correlates with the outflow gas mass. Together with the analysis of the outflow velocities, this suggests that at least some of the outflows are associated with the ejection of dusty clouds from the disc. The presence of dust in outflows is relevant for potential formation of molecules inside them. We combine our data with millimetre data to investigate the molecular phase. We find that the molecular phase accounts for more than 60 $\%$ of the total mass outflow rate in most objects and this fraction is higher in AGN-dominated systems. The neutral atomic phase contributes of the order of 10 $\%$, while the ionized phase is negligible. The ionized-to-molecular mass outflow rate declines slightly with AGN luminosity, although with a large scatter.

The Galactic Center black hole Sagittarius A* is a variable NIR source that exhibits bright flux excursions called flares. The low-flux density turnover of the flux distribution is below the sensitivity of current single-aperture telescopes. We use the unprecedented resolution of the GRAVITY instrument at the VLTI. Our light curves are unconfused, overcoming the confusion limit of previous photometric studies. We analyze the light curves using standard statistical methods and obtain the flux distribution. We find that the flux distribution of SgrA* turns over at a median flux density of (1.1\pm0.3)mJy. We measure the percentiles of the flux distribution and use them to constrain the NIR K-band SED. Furthermore, we find that the flux distribution is intrinsically right-skewed to higher flux density in log space. Flux densities below 0.1mJy are hardly ever observed. In consequence, a single powerlaw or lognormal distribution does not suffice to describe the observed flux distribution in its entirety. However, if one takes into account a power law component at high flux densities, a lognormal distribution can describe the lower end of the observed flux distribution. We confirm the RMS-flux relation for Sgr~A* and find it to be linear for all flux densities in our observation. We conclude that Sgr~A* has two states: the bulk of the emission is generated in a lognormal process with a well-defined median flux density and this quiescent emission is supplemented by sporadic flares that create the observed power law extension of the flux distribution.

The star S2 orbiting the compact radio source Sgr A* is a precision probe of the gravitational field around the closest massive black hole (candidate). Over the last 2.7 decades we have monitored the star's radial velocity and motion on the sky, mainly with the SINFONI and NACO adaptive optics (AO) instruments on the ESO VLT, and since 2017, with the four-telescope interferometric beam combiner instrument GRAVITY. In this paper we report the first detection of the General Relativity (GR) Schwarzschild Precession (SP) in S2's orbit. Owing to its highly elliptical orbit (e = 0.88), S2's SP is mainly a kink between the pre-and post-pericentre directions of motion ~ +- 1 year around pericentre passage, relative to the corresponding Kepler orbit. The superb 2017-2019 astrometry of GRAVITY defines the pericentre passage and outgoing direction. The incoming direction is anchored by 118 NACO-AO measurements of S2's position in the infrared reference frame, with an additional 75 direct measurements of the S2-Sgr A* separation during bright states ('flares') of Sgr A*. Our 14-parameter model fits for the distance, central mass, the position and motion of the reference frame of the AO astrometry relative to the mass, the six parameters of the orbit, as well as a dimensionless parameter f_SP for the SP (f_SP = 0 for Newton and 1 for GR). From data up to the end of 2019 we robustly detect the SP of S2, del phi = 12' per orbital period. From posterior fitting and MCMC Bayesian analysis with different weighting schemes and bootstrapping we find f_SP = 1.10 +- 0.19. The S2 data are fully consistent with GR. Any extended mass inside S2's orbit cannot exceed ~ 0.1% of the central mass. Any compact third mass inside the central arcsecond must be less than about 1000 M_sun.

We report on the determination of electron densities, and their impact on the outflow masses and rates, measured in the central few hundred parsecs of 11 local luminous active galaxies. We show that the peak of the integrated line emission in the AGN is significantly offset from the systemic velocity as traced by the stellar absorption features, indicating that the profiles are dominated by outflow. In contrast, matched inactive galaxies are characterised by a systemic peak and weaker outflow wing. We present three independent estimates of the electron density in these AGN, discussing the merits of the different methods. The electron density derived from the [SII] doublet is significantly lower than than that found with a method developed in the last decade using auroral and transauroral lines, as well as a recently introduced method based on the ionization parameter. The reason is that, for gas photoionized by an AGN, much of the [SII] emission arises in an extended partially ionized zone where the implicit assumption that the electron density traces the hydrogen density is invalid. We propose ways to deal with this situation and we derive the associated outflow rates for ionized gas, which are in the range 0.001--0.5 M$_{\odot}$yr$^{-1}$ for our AGN sample. We compare these outflow rates to the relation between $\dot{M}_{out}$ and $L_{AGN}$ in the literature, and argue that it may need to be modified and rescaled towards lower mass outflow rates.

The inner regions of the discs of high-mass young stellar objects (HMYSOs) are still poorly known due to the small angular scales and the high visual extinction involved. We deploy near-infrared (NIR) spectro-interferometry to probe the inner gaseous disc in HMYSOs and investigate the origin and physical characteristics of the CO bandhead emission (2.3-2.4 $\mu$m). We present the first GRAVITY/VLTI observations at high spectral (R=4000) and spatial (mas) resolution of the CO overtone transitions in NGC 2024 IRS2. The continuum emission is resolved in all baselines and is slightly asymmetric, displaying small closure phases ($\leq$8$^{\circ}$). Our best ellipsoid model provides a disc inclination of 34$^{\circ}$$\pm$1$^{\circ}$, a disc major axis position angle of 166$^{\circ}$$\pm$1$^{\circ}$, and a disc diameter of 3.99$\pm$0.09 mas (or 1.69$\pm$0.04 au, at a distance of 423 pc). The small closure phase signals in the continuum are modelled with a skewed rim, originating from a pure inclination effect. For the first time, our observations spatially and spectrally resolve the first four CO bandheads. Changes in visibility, as well as differential and closure phases across the bandheads are detected. Both the size and geometry of the CO-emitting region are determined by fitting a bidimensional Gaussian to the continuum-compensated CO bandhead visibilities. The CO-emitting region has a diameter of 2.74$\pm^{0.08}_{0.07}$ mas (1.16$\pm$0.03 au), and is located in the inner gaseous disc, well within the dusty rim, with inclination and $PA$ matching the dusty disc geometry, which indicates that both dusty and gaseous discs are coplanar. Physical and dynamical gas conditions are inferred by modelling the CO spectrum. Finally, we derive a direct measurement of the stellar mass of $M_*\sim$14.7$^{+2}_{-3.6}$ M$_{\odot}$ by combining our interferometric and spectral modelling results.

Infrared observations of Sgr A* probe the region close to the event horizon of the black hole at the Galactic center. These observations can constrain the properties of low-luminosity accretion as well as that of the black hole itself. The GRAVITY instrument at the ESO VLTI has recently detected continuous circular relativistic motion during infrared flares which has been interpreted as orbital motion near the event horizon. Here we analyze the astrometric data from these flares, taking into account the effects of out-of-plane motion and orbital shear of material near the event horizon of the black hole. We have developed a new code to predict astrometric motion and flux variability from compact emission regions following particle orbits. Our code combines semi-analytic calculations of timelike geodesics that allow for out-of-plane or elliptical motions with ray tracing of photon trajectories to compute time-dependent images and light curves. We apply our code to the three flares observed with GRAVITY in 2018. We show that all flares are consistent with a hotspot orbiting at R$\sim$9 gravitational radii with an inclination of $i\sim140^\circ$. The emitting region must be compact and less than $\sim5$ gravitational radii in diameter. We place a further limit on the out-of-plane motion during the flare.

Aims. We aim to search and characterize inflows and outflows of molecular gas in four ultraluminous infrared galaxies (ULIRGs) at $z\sim0.2-0.3$ and one distant QSO at $z=6.13$. Methods. We use Herschel PACS and ALMA Band 7 observations of the hydroxyl molecule (OH) line at rest-frame wavelength 119 $\mu$m which in absorption can provide unambiguous evidence for inflows or outflows of molecular gas in nuclear regions of galaxies. Our study contributes to double the number of OH observations of luminous systems at $z\sim0.2-0.3$, and push the search for molecular outflows based on the OH transition to $z\sim6$. Results. We detect OH high-velocity absorption wings in three of the four ULIRGs. In two cases, IRAS F20036-1547 and IRAS F13352+6402, the blueshifted absorption profiles indicate the presence of powerful and fast molecular gas outflows. Consistent with an inside-out quenching scenario, these outflows are depleting the central reservoir of molecular gas at a similar rate than the intense star formation activity. In the case of the starburst-dominated system IRAS 10091+4704, we detect an inverted P-Cygni profile that is unique among ULIRGs and indicates the presence of a fast ($\sim400$ km s$^{-1}$) inflow of molecular gas at a rate of $\sim100~M_{\odot}~{\rm yr}^{-1}$ towards the central region. Finally, we tentatively detect ($\sim3\sigma$) the OH doublet in absorption in the $z=6.13$ QSO ULAS J131911+095051. The OH feature is blueshifted with a median velocity that suggests the presence of a molecular outflow, although characterized by a modest molecular mass loss rate of $\sim200~M_{\odot}~{\rm yr}^{-1}$. This value is comparable to the small mass outflow rates found in the stacking of the [CII] spectra of other $z\sim6$ QSOs and suggests that ejective feedback in this phase of the evolution of ULAS J131911+095051 has subsided.

Our objective is to estimate the C/O ratio in the atmosphere of beta Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry. We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain K-band spectro-interferometric data on $\beta$ Pic b. We extracted a medium resolution (R=500) K-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio. Our new astrometry disfavors a circular orbit for $\beta$ Pic b ($e=0.15^{+0.05}_{-0.04}$). Combined with previous results and with Hipparcos/GAIA measurements, this astrometry points to a planet mass of $M = 12.7\pm{}2.2\,M_\mathrm{Jup}$. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of beta Pic b. In particular, the C/O ratios derived with the two codes are identical ($0.43\pm{}0.05$ vs $0.43^{+0.04}_{-0.03}$). We argue that if the stellar C/O in $\beta$ Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.

We present near-infrared interferometric data on the Seyfert 2 galaxy NGC 1068, obtained with the GRAVITY instrument on the European Southern Observatory Very Large Telescope Interferometer. The extensive baseline coverage from 5 to 60 M\lambda allowed us to reconstruct a continuum image of the nucleus with an unrivaled 0.2 pc resolution in the K-band. We find a thin ring-like structure of emission with a radius r = 0.24+/-0.03 pc, inclination i = 70+/-5 deg, position angle PA = -50+/-4 deg, and h/r < 0.14, which we associate with the dust sublimation region. The observed morphology is inconsistent with the expected signatures of a geometrically and optically thick torus. Instead, the infrared emission shows a striking resemblance to the 22 GHz maser disc, which suggests they share a common region of origin. The near-infrared spectral energy distribution indicates a bolometric luminosity of (0.4-4.7) x 10^45 erg/s, behind a large A_K ~ 5.5 (A_V ~ 90) screen of extinction that also appears to contribute significantly to obscuring the broad line region.

To understand the role that AGN feedback plays in galaxy evolution we need in-depth studies of the multi-phase structure and energetics of galaxy-wide outflows. In this work we present new, deep ($\sim$50 hr) NOEMA CO(1-0) line observations of the molecular gas in the powerful outflow driven by the AGN in the ultra-luminous infrared galaxy IRAS F08572+3915. We spatially resolve the outflow, finding that its most likely configuration is a wide-angle bicone aligned with the kinematic major axis of the rotation disk. The molecular gas in the wind reaches velocities up to approximately $\pm$1200 km s$^{-1}$ and transports nearly 20% of the molecular gas mass in the system. We detect a second outflow component located $\sim$6 kpc north-west from the galaxy moving away at $\sim$900 km s$^{-1}$, which could be the result of a previous episode of AGN activity. The total mass and energetics of the outflow, which includes contributions from the ionized, neutral, warm and cold molecular gas phases is strongly dominated by the cold molecular gas. In fact, the molecular mass outflow rate is higher than the star formation rate, even if we only consider the gas in the outflow that is fast enough to escape the galaxy, which accounts for about $\sim$40% of the total mass of the outflow. This results in an outflow depletion time for the molecular gas in the central $\sim$1.5 kpc region of only $\sim3$ Myr, a factor of $\sim2$ shorter than the depletion time by star formation activity.

We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high FIR surface brightness, and combine with results from the literature. CO line ratios and outflow structure provide some constraints on the conversion from observables to quantities such as molecular mass outflow rates. Ratios between outflow emission in higher J CO transitions and in CO(1-0) typically are consistent with excitation Ri1<~1. For IRAS 13120-5453, however, R31=2.10 indicates optically thin CO in the outflow. Like much of the outflow literature, we use alpha(CO) = 0.8, and we present arguments for using C=1 in deriving molecular mass outflow rates Mdot = C*M*v/R. We compare the two main methods for molecular outflow detection: CO mm interferometry and Herschel OH spectroscopy. For 26 sources studied with both methods, we find 80% agreement in detecting vout>~150km/s outflows, and non-matches can be plausibly ascribed to outflow geometry and SNR. For 12 bright ULIRGs with detailed OH-based outflow modeling, CO outflows are detected in all but one. Outflow masses, velocities, and sizes for these 11 sources agree well between the two methods, and modest remaining differences may relate to the different but overlapping regions sampled by CO emission and OH absorption. Outflow properties correlate better with AGN luminosity and with bolometric luminosity than with FIR surface brightness. The most massive outflows are found for systems with current AGN activity, but significant outflows in non-AGN systems must relate to star formation or to AGN activity in the recent past. We report scaling relations for the increase of outflow mass, rate, momentum rate, and kinetic power with bolometric luminosity. Short ~10^6yr flow times and some sources with resolved multiple outflow episodes support a role of intermittent driving, likely by AGN. (abridged)

The formation and the evolution of protoplanetary disks are important stages in the lifetime of stars. The processes of disk evolution and planet formation are intrinsically linked. We spatially resolve with GRAVITY/VLTI in the K-band the sub au-scale region of 27 stars to gain statistical understanding of their properties. We look for correlations with stellar parameters, such as luminosity, mass, temperature and age. Our sample also cover a range of various properties in terms of reprocessed flux, flared or flat morphology, and gaps. We developed semi-physical geometrical models to fit our interferometric data. Our best models correspond to smooth and wide rings, implying that wedge-shaped rims at the dust sublimation edge are favored, as found in the H-band. The closure phases are generally non-null with a median value of ~10 deg, indicating spatial asymmetries of the intensity distributions. Multi-size grain populations could explain the closure phase ranges below 20-25 deg but other scenarios should be invoked to explain the largest ones. Our measurements extend the Radius-Luminosity relation to ~1e4 Lsun and confirm the significant spread around the mean relation observed in the H-band. Gapped sources exhibit a large N-to-K band size ratio and large values of this ratio are only observed for the members of our sample that would be older than 1 Ma, less massive, and with lower luminosity. In the 2 Ms mass range, we observe a correlation in the increase of the relative age with the transition from group II to group I, and an increase of the N-to-K size ratio. However, the size of the current sample does not yet permit us to invoke a clear universal evolution mechanism across the HAeBe mass range. The measured locations of the K-band emission suggest that these disks might be structured by forming young planets, rather than by depletion due to EUV, FUV, and X-ray photo-evaporation.

We report on our combined analysis of HST, VLT/MUSE, VLT/SINFONI, and ALMA observations of the local Seyfert 2 galaxy, NGC 5728 to investigate in detail the feeding and feedback of the AGN. The datasets simultaneously probe the morphology, excitation, and kinematics of the stars, ionized gas, and molecular gas over a large range of spatial scales (10 pc--10 kpc). NGC 5728 contains a large stellar bar which is driving gas along prominent dust lanes to the inner 1 kpc where the gas settles into a circumnuclear ring. The ring is strongly star forming and contains a substantial population of young stars as indicated by the lowered stellar velocity dispersion and gas excitation consistent with HII regions. We model the kinematics of the ring using the velocity field of the CO (2--1) emission and stars and find it is consistent with a rotating disk. The outer regions of the disk, where the dust lanes meet the ring, show signatures of inflow at a rate of 1 M$_{\sun}$ yr$^{-1}$. Inside the ring, we observe three molecular gas components corresponding to the circular rotation of the outer ring, a warped disk, and the nuclear stellar bar. The AGN is driving an ionized gas outflow that reaches a radius of 250 pc with a mass outflow rate of 0.08 M$_{\sun}$ yr$^{-1}$ consistent with its luminosity and scaling relations from previous studies. While we observe distinct holes in CO emission which could be signs of molecular gas removal, we find that largely the AGN is not disrupting the structure of the circumnuclear region.

We use VLTI/GRAVITY near-infrared interferometry measurements of 8 bright, Type 1 AGN to study the size and structure of hot dust heated by the central engine. We partially resolve each source, and report Gaussian FWHM sizes in the range 0.3-0.8 milliarcseconds. In all but one object, we find no evidence for significant elongation or asymmetry (closure phases < 1 deg). The effective physical radius increases with bolometric luminosity as found from past reverberation and interferometry measurements. The measured sizes for Seyfert galaxies are systematically larger than for the two quasars in our sample when measured relative to the previously reported R ~ L^1/2 relationship explained by emission at the sublimation radius. This could be evidence of evolving near-infrared emission region structure as a function of central luminosity.

Precise measurements of the S-stars orbiting SgrA* have set strong constraints on the nature of the compact object at the centre of the Milky Way. The presence of a black hole in that region is well established, but its neighboring environment is still an open debate. In that respect, the existence of dark matter in that central region may be detectable due to its strong signatures on the orbits of stars: the main effect is a Newtonian precession which will affect the overall pericentre shift of S2, the latter being a target measurement of the GRAVITY instrument. The exact nature of this dark matter (e.g., stellar dark remnants or diffuse dark matter) is unknown. This article assumes it to be an scalar field of toroidal distribution, associated with ultra-light dark matter particles, surrounding the Kerr black hole. Such a field is a form of "hair" expected in the context of superradiance, a mechanism that extracts rotational energy from the black hole. Orbital signatures for the S2 star are computed and shown to be detectable by GRAVITY. The scalar field can be constrained because the variation of orbital elements depends both on the relative mass of the scalar field to the black hole and on the field mass coupling parameter.

We present a 0.16% precise and 0.27% accurate determination of R0, the distance to the Galactic Center. Our measurement uses the star S2 on its 16-year orbit around the massive black hole Sgr A* that we followed astrometrically and spectroscopically for 27 years. Since 2017, we added near-infrared interferometry with the VLTI beam combiner GRAVITY, yielding a direct measurement of the separation vector between S2 and Sgr A* with an accuracy as good as 20 micro-arcsec in the best cases. S2 passed the pericenter of its highly eccentric orbit in May 2018, and we followed the passage with dense sampling throughout the year. Together with our spectroscopy, in the best cases with an error of 7 km/s, this yields a geometric distance estimate: R0 = 8178 +- 13(stat.) +- 22(sys.) pc. This work updates our previous publication in which we reported the first detection of the gravitational redshift in the S2 data. The redshift term is now detected with a significance level of 20 sigma with f_redshift = 1.04 +- 0.05.

To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes on the central star, and integrated off-axis on the HR8799e planet situated at 390 mas from the star. Data reduction included post-processing to remove the flux leaking from the central star and to extract the coherent flux of the planet. The inferred K band spectrum of the planet has a spectral resolution of 500. We also derive the astrometric position of the planet relative to the star with a precision on the order of 100$\,\mu$as. The GRAVITY astrometric measurement disfavors perfectly coplanar stable orbital solutions. A small adjustment of a few degrees to the orbital inclination of HR 8799 e can resolve the tension, implying that the orbits are close to, but not strictly coplanar. The spectrum, with a signal-to-noise ratio of $\approx 5$ per spectral channel, is compatible with a late-type L brown dwarf. Using Exo-REM synthetic spectra, we derive a temperature of $1150\pm50$\,K and a surface gravity of $10^{4.3\pm0.3}\,$cm/s$^{2}$. This corresponds to a radius of $1.17^{+0.13}_{-0.11}\,R_{\rm Jup}$ and a mass of $10^{+7}_{-4}\,M_{\rm Jup}$, which is an independent confirmation of mass estimates from evolutionary models. Our results demonstrate the power of interferometry for the direct detection and spectroscopic study of exoplanets at close angular separations from their stars.

During its orbit around the four million solar mass black hole Sagittarius A* the star S2 experiences significant changes in gravitational potential. We use this change of potential to test one part of the Einstein equivalence principle: the local position invariance (LPI). We study the dependency of different atomic transitions on the gravitational potential to give an upper limit on violations of the LPI. This is done by separately measuring the redshift from hydrogen and helium absorption lines in the stellar spectrum during its closest approach to the black hole. For this measurement we use radial velocity data from 2015 to 2018 and combine it with the gravitational potential at the position of S2, which is calculated from the precisely known orbit of S2 around the black hole. This results in a limit on a violation of the LPI of $|\beta_{He}-\beta_{H}| = (2.4 \pm 5.1) \cdot 10^{-2}$. The variation in potential that we probe with this measurement is six magnitudes larger than possible for measurements on Earth, and a factor ten larger than in experiments using white dwarfs. We are therefore testing the LPI in a regime where it has not been tested before.

The broadening of atomic emission lines by high-velocity motion of gas near accreting supermassive black holes is an observational hallmark of quasars. Observations of broad emission lines could potentially constrain the mechanism for transporting gas inwards through accretion disks or outwards through winds. The size of this broad-line region has been estimated by measuring the light travel time delay between the variable nuclear continuum and the emission lines - a method known as reverberation mapping. In some models the emission lines arise from a continuous outflow, whereas in others they are produced by orbiting gas clouds. Directly imaging such regions has not hitherto been possible because of their small angular sizes (< 0.1 milli-arcseconds). Here we report a spatial offset (with a spatial resolution of ten micro-arcseconds or about 0.03 parsecs for a distance of 550 million parsecs) between the red and blue photo-centres of the broad Paschen-{\alpha} line of the quasar 3C 273 perpendicular to the direction of its radio jet. This spatial offset corresponds to a gradient in the velocity of the gas and thus implies that the gas is orbiting the central supermassive black hole. The data are well fitted by a broad-line-region model of a thick disk of gravitationally bound material orbiting a black hole of 300 million solar masses. We infer a disk radius of 150 light days; a radius of 100-400 light days was found previously using reverberation mapping. The rotation axis of the disk aligns in inclination and position angle with the radio jet. Our results support the methods that are often used to estimate the masses of accreting supermassive black holes and to study their evolution over cosmic time.

We report the detection of continuous positional and polarization changes of the compact source SgrA* in high states ('flares') of its variable near- infrared emission with the near-infrared GRAVITY-Very Large Telescope Interferometer (VLTI) beam-combining instrument. In three prominent bright flares, the position centroids exhibit clockwise looped motion on the sky, on scales of typically 150 micro-arcseconds over a few tens of minutes, corresponding to about 30% the speed of light. At the same time, the flares exhibit continuous rotation of the polarization angle, with about the same 45(+/-15)-minute period as that of the centroid motions. Modelling with relativistic ray tracing shows that these findings are all consistent with a near face-on, circular orbit of a compact polarized 'hot spot' of infrared synchrotron emission at approximately six to ten times the gravitational radius of a black hole of 4 million solar masses. This corresponds to the region just outside the innermost, stable, prograde circular orbit (ISCO) of a Schwarzschild-Kerr black hole, or near the retrograde ISCO of a highly spun-up Kerr hole. The polarization signature is consistent with orbital motion in a strong poloidal magnetic field.

This work presents an interferometric study of the massive-binary fraction in the Orion Trapezium Cluster with the recently comissioned GRAVITY instrument. We observe a total of 16 stars of mainly OB spectral type. We find three previously unknown companions for $\theta ^1$ Ori B, $\theta ^2$ Ori B, and $\theta ^2$ Ori C. We determine a separation for the previously suspected companion of NU Ori. We confirm four companions for $\theta ^1$ Ori A, $\theta ^1$ Ori C, $\theta ^1$ Ori D, and $\theta ^2$ Ori A, all with substantially improved astrometry and photometric mass estimates. We refine the orbit of the eccentric high-mass binary $\theta ^1$ Ori C and we are able to derive a new orbit for $\theta ^1$ Ori D. We find a system mass of 21.7 $M_{\odot}$ and a period of $53$ days. Together with other previously detected companions seen in spectroscopy or direct imaging, eleven of the 16 high-mass stars are multiple systems. We obtain a total number of 22 companions with separations up to 600 AU. The companion fraction of the early B and O stars in our sample is about 2, significantly higher than in earlier studies of mostly OB associations. The separation distribution hints towards a bimodality. Such a bimodality has been previously found in A stars, but rarely in OB binaries, which up to this point have been assumed to be mostly compact with a tail of wider companions. We also do not find a substantial population of equal-mass binaries. The observed distribution of mass ratios declines steeply with mass, and like the direct star counts, indicates that our companions follow a standard power law initial mass function. Again, this is in contrast to earlier findings of flat mass ratio distributions in OB associations. We exclude collision as a dominant formation mechanism but find no clear preference for core accretion or competitive accretion.

Eta Car is one of the most intriguing luminous blue variables in the Galaxy. Observations and models at different wavelengths suggest a central binary with a 5.54 yr period residing in its core. 2D and 3D radiative transfer and hydrodynamic simulations predict a primary with a dense and slow stellar wind that interacts with the faster and lower density wind of the secondary. The wind-wind collision scenario suggests that the secondary's wind penetrates the primary's wind creating a low-density cavity in it, with dense walls where the two winds interact. We aim to trace the inner ~5-50 au structure of Eta Car's wind-wind interaction, as seen through BrG and, for the first time, through the He I 2s-2p line. We have used spectro-interferometric observations with GRAVITY at the VLTI. Our modeling of the continuum allows us to estimate its FWHM angular size close to 2 mas and an elongation ratio of 1.06 +/- 0.05 over a PA = 130 +/- 20 deg. Our CMFGEN modeling helped us to confirm that the role of the secondary should be taken into account to properly reproduce the observed BrG and He I lines. Chromatic images across BrG reveal a southeast arc-like feature, possibly associated to the hot post-shocked winds flowing along the cavity wall. The images of He I 2s-2p served to constrain the 20 mas structure of the line-emitting region. The observed morphology of He I suggests that the secondary is responsible for the ionized material that produces the line profile. Both the BrG and the He I 2s-2p maps are consistent with previous hydrodynamical models of the colliding wind scenario. Future dedicated simulations together with an extensive interferometric campaign are necessary to refine our constraints on the wind and stellar parameters of the binary, which finally will help us predict the evolutionary path of Eta Car.

MICADO will enable the ELT to perform diffraction limited near-infrared observations at first light. The instrument's capabilities focus on imaging (including astrometric and high contrast) as well as single object spectroscopy. This contribution looks at how requirements from the observing modes have driven the instrument design and functionality. Using examples from specific science cases, and making use of the data simulation tool, an outline is presented of what we can expect the instrument to achieve.

The highly elliptical, 16-year-period orbit of the star S2 around the massive black hole candidate Sgr A* is a sensitive probe of the gravitational field in the Galactic centre. Near pericentre at 120 AU, ~1400 Schwarzschild radii, the star has an orbital speed of ~7650 km/s, such that the first-order effects of Special and General Relativity have now become detectable with current capabilities. Over the past 26 years, we have monitored the radial velocity and motion on the sky of S2, mainly with the SINFONI and NACO adaptive optics instruments on the ESO Very Large Telescope, and since 2016 and leading up to the pericentre approach in May 2018, with the four-telescope interferometric beam-combiner instrument GRAVITY. From data up to and including pericentre, we robustly detect the combined gravitational redshift and relativistic transverse Doppler effect for S2 of z ~ 200 km/s / c with different statistical analysis methods. When parameterising the post-Newtonian contribution from these effects by a factor f, with f = 0 and f = 1 corresponding to the Newtonian and general relativistic limits, respectively, we find from posterior fitting with different weighting schemes f = 0.90 +/- 0.09 (stat) +\- 0.15 (sys). The S2 data are inconsistent with pure Newtonian dynamics.

Nuclear outflows driven by accreting massive black holes are one of the main feedback mechanisms invoked at high-z to reproduce the distinct separation between star-forming, disk galaxies and quiescent spheroidal systems. Yet, our knowledge of feedback at high-z remains limited by the lack of observations of the multiple gas phases in galaxy outflows. In this work we use new deep, high-spatial resolution ALMA CO(3-2) and archival VLT/SINFONI H$\alpha$ observations to study the molecular and ionized components of the AGN-driven outflow in zC400528 ---a massive, main sequence galaxy at z=2.3 in the process of quenching. We detect a powerful molecular outflow that shows a positive velocity gradient and extends for at least ~10 kpc from the nuclear region, about three times the projected size of the ionized wind. The molecular gas in the outflow does not reach velocities high enough to escape the galaxy and is therefore expected to be reaccreted. Keeping in mind the various assumptions involved in the analysis, we find that the mass and energetics of the outflow are dominated by the molecular phase. The AGN-driven outflow in zC400528 is powerful enough to deplete the molecular gas reservoir on a timescale at least twice shorter than that needed to exhaust it by star formation. This suggests that the nuclear outflow is one of the main quenching engines at work in the observed suppression of the central star-formation activity in zC400528.

ERIS is an instrument that will both extend and enhance the fundamental diffraction limited imaging and spectroscopy capability for the VLT. It will replace two instruments that are now being maintained beyond their operational lifetimes, combine their functionality on a single focus, provide a new wavefront sensing module that makes use of the facility Adaptive Optics System, and considerably improve their performance. The instrument will be competitive with respect to JWST in several regimes, and has outstanding potential for studies of the Galactic Center, exoplanets, and high redshift galaxies. ERIS had its final design review in 2017, and is expected to be on sky in 2020. This contribution describes the instrument concept, outlines its expected performance, and highlights where it will most excel.

The oxygen-bearing molecular ions OH+, H2O+, and H3O+ are key species that probe the ionization rate of (partially) molecular gas that is ionized by X-rays and cosmic rays permeating the interstellar medium. We report Herschel far-infrared and submillimeter spectroscopic observations of OH+ in Mrk 231, showing both ground-state P-Cygni profiles, and excited line profiles with blueshifted absorption wings extending up to ~1000 km s^{-1}. In addition, OH+ probes an excited component peaking at central velocities, likely arising from the torus probed by the OH centimeter-wave megamaser. Four lines of H2O+ are also detected at systemic velocities, but H3O+ is undetected. Based on our earlier OH studies, we estimate an abundance ratio of OH/OH+~5-10 for the outflowing components and ~20 for the torus, and an OH+ abundance relative to H nuclei of ~>10^{-7}. We also find high OH+/H2O+ and OH+/H3O+ ratios, both are ~>4 in the torus and ~>10-20 in the outflowing gas components. Chemical models indicate that these high OH+ abundances relative to OH, H2O+, and H3O+ are characteristic of gas with a high ionization rate per unit density, \zeta/n_H~(1-5)x10^{-17} cm^3 s^{-1} and ~(1-2)x10^{-16} cm^3 s^{-1} for the above components, respectively, and an ionization rate of \zeta~(0.5-2)x10^{-12} s^{-1}. X-rays appear to be unable to explain the inferred ionization rate, and thus we suggest that low-energy (10-400 MeV) cosmic-rays are primarily responsible for the ionization with \dot{M}_{CR}~0.01 M_{sun} yr^{-1} and \dot{E}_{CR}~10^{44} erg s^{-1}, the latter corresponding to 1% of the AGN luminosity and similar to the energetics of the molecular outflow. We suggest that cosmic-rays accelerated in the forward shock associated with the molecular outflow are responsible for the ionization, as they diffuse through the outflowing molecular phase downstream.

We use the Herschel/PACS spectrometer to study the global and spatially resolved far-infrared (FIR) fine-structure line emission in a sample of 52 galaxies that constitute the SHINING survey. These galaxies include star-forming, active-galactic nuclei (AGN), and luminous infrared galaxies (LIRGs). We find an increasing number of galaxies (and kiloparsec size regions within galaxies) with low line-to-FIR continuum ratios as a function of increasing FIR luminosity ($L_{\mathrm{FIR}}$), dust infrared color, $L_{\mathrm{FIR}}$ to molecular gas mass ratio ($L_{\mathrm{FIR}}/M_{\mathrm{mol}}$), and FIR surface brightness ($\Sigma_{\mathrm{FIR}}$). The correlations between the [CII]/FIR or [OI]/FIR ratios with $\Sigma_{\mathrm{FIR}}$ are remarkably tight ($\sim0.3$ dex scatter over almost four orders of magnitude in $\Sigma_{\mathrm{FIR}}$). We observe that galaxies with $L_{\mathrm{FIR}}/M_{\mathrm{mol}} \gtrsim 80\,L_{\odot}\,M_{\odot}^{-1}$ and $\Sigma_{\mathrm{FIR}}\gtrsim10^{11}$ $L_{\odot}$ kpc$^{-2}$ tend to have weak fine-structure line-to-FIR continuum ratios, and that LIRGs with infrared sizes $\gtrsim1$ kpc have line-to-FIR ratios comparable to those observed in typical star-forming galaxies. We analyze the physical mechanisms driving these trends in Paper II (Herrera-Camus et al. 2018). The combined analysis of the [CII], [NII], and [OIII] lines reveals that the fraction of the [CII] line emission that arises from neutral gas increases from 60% to 90% in the most active star-forming regions and that the emission originating in the ionized gas is associated with low-ionization, diffuse gas rather than with dense gas in HII regions. Finally, we report the global and spatially resolved line fluxes of the SHINING galaxies to enable the comparison and planning of future local and high-$z$ studies.

The SHINING survey (Paper I; Herrera-Camus et al. 2018) offers a great opportunity to study the properties of the ionized and neutral media of galaxies from prototypical starbursts and active galactic nuclei (AGN) to heavily obscured objects. Based on Herschel/PACS observations of the main far-infrared (FIR) fine-structure lines, in this paper we analyze the physical mechanisms behind the observed line deficits in galaxies, the apparent offset of luminous infrared galaxies (LIRGs) from the mass-metallicity relation, and the scaling relations between [CII] 158 $\mu$m line emission and star formation rate (SFR). Based on a toy model and the Cloudy code, we conclude that the increase in the ionization parameter with FIR surface brightness can explain the observed decrease in the line-to-FIR continuum ratio of galaxies. In the case of the [CII] line, the increase in the ionization parameter is accompanied by a reduction in the photoelectric heating efficiency and the inability of the line to track the increase in the FUV radiation field as galaxies become more compact and luminous. In the central $\sim$kiloparsec regions of AGN galaxies we observe a significant increase in the [OI] 63 $\mu$m/[CII] line ratio; the AGN impact on the line-to-FIR ratios fades on global scales. Based on extinction-insensitive metallicity measurements of LIRGs we confirm that they lie below the mass-metallicity relation, but the offset is smaller than those reported in studies that use optical-based metal abundances. Finally, we present scaling relations between [CII] emission and SFR in the context of the main-sequence of star-forming galaxies.

We explore the relationship between X-ray absorption and optical obscuration within the BAT AGN Spectroscopic Survey (BASS) which has been collecting and analyzing the optical and X-ray spectra for 641 hard X-ray selected ($E>14$ keV) active galactic nuclei (AGN). We use the deviation from a linear broad H$\alpha$-to-X-ray relationship as an estimate of the maximum optical obscuration towards the broad line region and compare the $A_{\rm V}$ to the hydrogen column densities ($N_{\rm H}$) found through systematic modeling of their X-ray spectra. We find that the inferred columns implied by $A_{\rm V}$ towards the broad line region (BLR) are often orders of magnitude less than the columns measured towards the X-ray emitting region indicating a small scale origin for the X-ray absorbing gas. After removing 30\% of Sy 1.9s that potentially have been misclassified due to outflows, we find that 86\% (164/190) of the Type 1 population (Sy 1--1.9) are X-ray unabsorbed as expected based on a single obscuring structure. However, 14\% (26/190), of which 70\% (18/26) are classified as Sy 1.9, are X-ray absorbed, suggesting the broad line region itself is providing extra obscuration towards the X-ray corona. The fraction of X-ray absorbed Type 1 AGN remains relatively constant with AGN luminosity and Eddington ratio, indicating a stable broad line region covering fraction.

Using new APEX and JCMT spectroscopy of the CO 2-1 line, we undertake a controlled study of cold molecular gas in moderately luminous Active Galactic Nuclei (AGN) and inactive galaxies from the Luminous Local AGN with Matched Analogs (LLAMA) survey. We use spatially resolved infrared photometry of the LLAMA galaxies from 2MASS, WISE, IRAS & Herschel, corrected for nuclear emission using multi-component spectral energy distribution (SED) fits, to examine the dust-reprocessed star-formation rates (SFRs), molecular gas fractions and star formation efficiencies (SFEs) over their central 1 - 3 kpc. We find that the gas fractions and central SFEs of both active and inactive galaxies are similar when controlling for host stellar mass and morphology (Hubble type). The equivalent central molecular gas depletion times are consistent with the discs of normal spiral galaxies in the local Universe. Despite energetic arguments that the AGN in LLAMA should be capable of disrupting the observable cold molecular gas in their central environments, our results indicate that nuclear radiation only couples weakly with this phase. We find a mild preference for obscured AGN to contain higher amounts of central molecular gas, which suggests a connection between AGN obscuration and the gaseous environment of the nucleus. Systems with depressed SFEs are not found among the LLAMA AGN. We speculate that the processes that sustain the collapse of molecular gas into dense pre-stellar cores may also be a prerequisite for the inflow of material on to AGN accretion disks.

In a complete sample of local 14-195 keV selected AGNs and inactive galaxies, matched by their host galaxy properties, we study the spatially resolved stellar kinematics and luminosity distributions at near-infrared wavelengths on scales of 10-150 pc, using SINFONI on the VLT. In this paper, we present the first half of the sample, which comprises 13 galaxies, 8 AGNs and 5 inactive galaxies. The stellar velocity fields show a disk-like rotating pattern, for which the kinematic position angle is in agreement with the photometric position angle obtained from large scale images. For this set of galaxies, the stellar surface brightness of the inactive galaxy sample is generally comparable to the matched sample of AGN but extends to lower surface brightness. After removal of the bulge contribution, we find a nuclear stellar light excess with an extended nuclear disk structure, and which exhibits a size-luminosity relation. While we expect the excess luminosity to be associated with a dynamically cooler young stellar population, we do not typically see a matching drop in dispersion. This may be because these galaxies have pseudo-bulges in which the intrinsic dispersion increases towards the centre. And although the young stars may have an impact in the observed kinematics, their fraction is too small to dominate over the bulge and compensate the increase in dispersion at small radii, so no dispersion drop is seen. Finally, we find no evidence for a difference in the stellar kinematics and nuclear stellar luminosity excess between these active and inactive galaxies.

A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas. Outflowing gas is identified in the far-infrared through P-Cygni line shapes and absorption blueshifted wings in molecular lines with high dipolar moments, and through emission line wings of fine-structure lines of ionized gas. We quantify the detectability of galaxy-scale massive molecular and ionized outflows as a function of redshift in AGN-dominated, starburst-dominated, and main-sequence galaxies, explore the detectability of metal-rich inflows in the local Universe, and describe the most significant synergies with other current and future observatories that will measure feedback in galaxies via complementary tracers at other wavelengths.

IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features.

To investigate the inner regions of protoplanetary disks, we performed near-infrared interferometric observations of the classical TTauri binary system S CrA. We present the first VLTI-GRAVITY high spectral resolution ($R\sim$4000) observations of a classical TTauri binary, S CrA (composed of S CrA N and S CrA S and separated by $\sim$1.4"), combining the four 8-m telescopes in dual-field mode. Our observations in the near-infrared K-band continuum reveal a disk around each binary component, with similar half-flux radii of about 0.1 au at d$\sim$130 pc, inclinations ($i=$28$\pm$3$^o$\ and $i=$22$\pm$6$^o$), and position angles (PA=0$^o\pm$6$^o$ and PA=-2$^o\pm$12$^o$), suggesting that they formed from the fragmentation of a common disk. The S CrA N spectrum shows bright HeI and Br$\gamma$ line emission exhibiting inverse P-Cygni profiles, typically associated with infalling gas. The continuum-compensated Br$\gamma$ line visibilities of S CrA N show the presence of a compact Br$\gamma$ emitting region the radius of which is about $\sim$0.06 au, which is twice as big as the truncation radius. This component is mostly tracing a wind. Moreover, a slight radius change between the blue- and red-shifted Br$\gamma$ line components is marginally detected. The presence of an inverse P-Cygni profile in the HeI and Br$\gamma$ lines, along with the tentative detection of a slightly larger size of the blue-shifted Br$\gamma$ line component, hint at the simultaneous presence of a wind and magnetospheric accretion in S CrA N.

We use Herschel data to analyze the size of the far-infrared 70micron emission for z<0.06 local samples of 277 hosts of Swift-BAT selected active galactic nuclei (AGN), and 515 comparison galaxies that are not detected by BAT. For modest far-infrared luminosities 8.5<log(LFIR)<10.5, we find large scatter of half light radii Re70 for both populations, but a typical Re70 <~ 1 kpc for the BAT hosts that is only half that of comparison galaxies of same far-infrared luminosity. The result mostly reflects a more compact distribution of star formation (and hence gas) in the AGN hosts, but compact AGN heated dust may contribute in some extremely AGN-dominated systems. Our findings are in support of an AGN-host coevolution where accretion onto the central black hole and star formation are fed from the same gas reservoir, with more efficient black hole feeding if that reservoir is more concentrated. The significant scatter in the far-infrared sizes emphasizes that we are mostly probing spatial scales much larger than those of actual accretion, and that rapid accretion variations can smear the distinction between the AGN and comparison categories. Large samples are hence needed to detect structural differences that favour feeding of the black hole. No size difference AGN host vs. comparison galaxies is observed at higher far-infrared luminosities log(LFIR)>10.5 (star formation rates >~ 6 Msun/yr), possibly because these are typically reached in more compact regions in the first place.

We present an analysis of far--infrared (FIR) [CII] and [OI] fine structure line and continuum observations obtained with $Herschel$/PACS, and CO(1-0) observations obtained with the IRAM Plateau de Bure Interferometer, of Lyman Break Analogs (LBAs) at $z\sim 0.2$. The principal aim of this work is to determine the typical ISM properties of $z\sim 1-2$ Main Sequence (MS) galaxies, with stellar masses between $10^{9.5}$ and $10^{11}$ $M_{\odot}$, which are currently not easily detectable in all these lines even with ALMA and NOEMA. We perform PDR modeling and apply different IR diagnostics to derive the main physical parameters of the FIR emitting gas and dust and we compare the derived ISM properties to those of galaxies on and above the MS at different redshifts. We find that the ISM properties of LBAs are quite extreme (low gas temperature, high density and thermal pressure) with respect to those found in local normal spirals and more active local galaxies. LBAs have no [CII] deficit despite having the high specific star formation rates (sSFRs) typical of starbursts. Although LBAs lie above the local MS, we show that their ISM properties are more similar to those of high-redshift MS galaxies than of local galaxies above the main sequence. This data set represents an important reference for planning future ALMA [CII] observations of relatively low-mass MS galaxies at the epoch of the peak of the cosmic star formation.

In Tombesi et al. (2015), we reported the first direct evidence for a quasar accretion disk wind driving a massive molecular outflow. The target was F11119+3257, an ultraluminous infrared galaxy (ULIRG) with unambiguous type-1 quasar optical broad emission lines. The energetics of the accretion disk wind and molecular outflow were found to be consistent with the predictions of quasar feedback models where the molecular outflow is driven by a hot energy-conserving bubble inflated by the inner quasar accretion disk wind. However, this conclusion was uncertain because the energetics were estimated from the optically thick OH 119 um transition profile observed with Herschel. Here, we independently confirm the presence of the molecular outflow in F11119+3257, based on the detection of broad wings in the CO(1-0) profile derived from ALMA observations. The broad CO(1-0) line emission appears to be spatially extended on a scale of at least ~7 kpc from the center. Mass outflow rate, momentum flux, and mechanical power of (80-200) R_7^{-1} M_sun/yr, (1.5-3.0) R_7^{-1} L_AGN/c, and (0.15-0.40)% R_7^{-1} L_AGN are inferred from these data, assuming a CO-to-H_2 conversion factor appropriate for a ULIRG (R_7 is the radius of the outflow normalized to 7 kpc and L_AGN is the AGN luminosity). These rates are time-averaged over a flow time scale of 7x10^6 yrs. They are similar to the OH-based rates time-averaged over a flow time scale of 4x10^5 yrs, but about a factor 4 smaller than the local ("instantaneous"; <10^5 yrs) OH-based estimates cited in Tombesi et al. The implications of these new results are discussed in the context of time-variable quasar-mode feedback and galaxy evolution. The need for an energy-conserving bubble to explain the molecular outflow is also re-examined.

We present the first optical observation at sub-milliarcsecond (mas) scale of the microquasar SS 433 obtained with the GRAVITY instrument on the VLT interferometer. The 3.5 hour exposure reveals a rich K-band spectrum dominated by hydrogen Br$\gamma\ $ and \ion{He}{i} lines, as well as (red-shifted) emission lines coming from the jets. The K-band continuum emitting region is dominated by a marginally resolved point source ($<$ 1 mas) embedded inside a diffuse background accounting for 10\% of the total flux. The jet line positions agree well with the ones expected from the jet kinematic model, an interpretation also supported by the consistent sign (i.e. negative/positive for the receding/approaching jet component) of the phase shifts observed in the lines. The significant visibility drop across the jet lines, together with the small and nearly identical phases for all baselines, point toward a jet that is offset by less than 0.5 mas from the continuum source and resolved in the direction of propagation, with a typical size of 2 mas. The jet position angle of $\sim$80$^{\circ}$ is consistent with the expected one at the observation date. Jet emission so close to the central binary system would suggest that line locking, if relevant to explain the amplitude and stability of the 0.26c jet velocity, operates on elements heavier than hydrogen. The Br$\gamma\ $ profile is broad and double peaked. It is better resolved than the continuum and the change of the phase signal sign across the line on all baselines suggests an East-West oriented geometry alike the jet direction and supporting a (polar) disk wind origin.

We observe the HMXB BP Cru using interferometry in the near-infrared K band with VLTI/GRAVITY. Continuum visibilities are at most partially resolved, consistent with the predicted size of the hypergiant. Differential visibility amplitude ($\Delta |V| \sim 5\%$) and phase ($\Delta \phi \sim 2 \deg$) signatures are observed across the HeI $2.059 \mu$m and Br$\gamma$ lines, the latter seen strongly in emission, unusual for the donor star's spectral type. For a baseline $B \sim 100$m, the differential phase RMS $\sim 0.2 \deg$ corresponds to an astrometric precision of $\sim 2 \mu$as. A model-independent analysis in the marginally resolved limit of interferometry reveals asymmetric and extended emission with a strong wavelength dependence. We propose geometric models based on an extended and distorted wind and/or a high density gas stream, which has long been predicted to be present in this system. The observations show that optical interferometry is now able to resolve HMXBs at the spatial scale at which accretion takes place, and therefore probe the effects of the gravitational and radiation fields of the compact object on its environment.

GRAVITY is a new instrument to coherently combine the light of the European Southern Observatory Very Large Telescope Interferometer to form a telescope with an equivalent 130 m diameter angular resolution and a collecting area of 200 m$^2$. The instrument comprises fiber fed integrated optics beam combination, high resolution spectroscopy, built-in beam analysis and control, near-infrared wavefront sensing, phase-tracking, dual beam operation and laser metrology [...]. This article gives an overview of GRAVITY and reports on the performance and the first astronomical observations during commissioning in 2015/16. We demonstrate phase tracking on stars as faint as m$_K$ ~ 10 mag, phase-referenced interferometry of objects fainter than m$_K$ ~ 15 mag with a limiting magnitude of m$_K$ ~ 17 mag, minute long coherent integrations, a visibility accuracy of better than 0.25 %, and spectro-differential phase and closure phase accuracy better than 0.5{\deg}, corresponding to a differential astrometric precision of better than 10 microarcseconds ({\mu}as). The dual-beam astrometry, measuring the phase difference of two objects with laser metrology, is still under commissioning. First observations show residuals as low as 50 {\mu}as when following objects over several months. We illustrate the instrument performance with the observations of archetypical objects for the different instrument modes. Examples include the Galactic Center supermassive black hole and its fast orbiting star S2 for phase referenced dual beam observations and infrared wavefront sensing, the High Mass X-Ray Binary BP Cru and the Active Galactic Nucleus of PDS 456 for few {\mu}as spectro-differential astrometry, the T Tauri star S CrA for a spectro-differential visibility analysis, {\xi} Tel and 24 Cap for high accuracy visibility observations, and {\eta} Car for interferometric imaging with GRAVITY.

Recent observations have revealed massive galactic molecular outflows that may have physical conditions (high gas densities) required to form stars. Indeed, several recent models predict that such massive galactic outflows may ignite star formation within the outflow itself. This star-formation mode, in which stars form with high radial velocities, could contribute to the morphological evolution of galaxies, to the evolution in size and velocity dispersion of the spheroidal component of galaxies, and would contribute to the population of high-velocity stars, which could even escape the galaxy. Such star formation could provide in-situ chemical enrichment of the circumgalactic and intergalactic medium (through supernova explosions of young stars on large orbits), and some models also predict that it may contribute substantially to the global star formation rate observed in distant galaxies. Although there exists observational evidence for star formation triggered by outflows or jets into their host galaxy, as a consequence of gas compression, evidence for star formation occurring within galactic outflows is still missing. Here we report new spectroscopic observations that unambiguously reveal star formation occurring in a galactic outflow at a redshift of 0.0448. The inferred star formation rate in the outflow is larger than 15 Msun/yr. Star formation may also be occurring in other galactic outflows, but may have been missed by previous observations owing to the lack of adequate diagnostics.

We report on the energetics of molecular outflows in 14 local Ultraluminous Infrared Galaxies (ULIRGs) that show unambiguous outflow signatures (P-Cygni profiles or high-velocity absorption wings) in the far-infrared lines of OH measured with the Herschel/PACS spectrometer. Detection of both ground-state (at 119 and 79 um) and one or more radiatively-excited (at 65 and 84 um) lines allows us to model the nuclear gas (<~300 pc) as well as the more extended components using spherically symmetric radiative transfer models. The highest molecular outflow velocities are found in buried sources, in which slower but massive expansion of the nuclear gas is also observed. With the exception of a few outliers, the outflows have momentum fluxes of (2-5)xL_IR/c and mechanical luminosities of (0.1-0.3)% of L_IR. The moderate momentum boosts in these sources (<~3) suggest that the outflows are mostly momentum-driven by the combined effects of AGN and nuclear starbursts, as a result of radiation pressure, winds, and supernovae remnants. In some sources (~20%), however, powerful (10^{10.5-11} Lsun) AGN feedback and (partially) energy-conserving phases are required, with momentum boosts in the range 3-20. These outflows appear to be stochastic strong-AGN feedback events that occur throughout the merging process. In a few sources, the outflow activity in the innermost regions has subsided in the last ~1 Myr. While OH traces the molecular outflows at sub-kpc scales, comparison of the masses traced by OH with those previously inferred from tracers of more extended outflowing gas suggests that most mass is loaded (with loading factors of Mdot/SFR=1-10) from the central galactic cores (a few x 100 pc). Outflow depletion timescales are <10^8 yr, shorter than the gas consumption timescales by factors of 1.1-15, and are anti-correlated with the AGN luminosity.

We present Herschel PACS observations of the [CII] 158 micron emission line in a sample of 24 intermediate mass (9<logM$_\ast$/M$_\odot$<10) and low metallicity (0.4< Z/Z$_\odot$<1.0) galaxies from the xCOLD GASS survey. Combining them with IRAM CO(1-0) measurements, we establish scaling relations between integrated and molecular region [CII]/CO(1-0) luminosity ratios as a function of integrated galaxy properties. A Bayesian analysis reveals that only two parameters, metallicity and offset from the star formation main sequence, $\Delta$MS, are needed to quantify variations in the luminosity ratio; metallicity describes the total dust content available to shield CO from UV radiation, while $\Delta$MS describes the strength of this radiation field. We connect the [CII]/CO luminosity ratio to the CO-to-H$_2$ conversion factor and find a multivariate conversion function $\alpha_{CO}$, which can be used up to z~2.5. This function depends primarily on metallicity, with a second order dependence on $\Delta$MS. We apply this to the full xCOLD GASS and PHIBSS1 surveys and investigate molecular gas scaling relations. We find a flattening of the relation between gas mass fraction and stellar mass at logM$_\ast$/M$_\odot$<10. While the molecular gas depletion time varies with sSFR, it is mostly independent of mass, indicating that the low L$_{CO}$/SFR ratios long observed in low mass galaxies are entirely due to photodissociation of CO, and not to an enhanced star formation efficiency.

We discuss the environment of local hard X-ray selected active galaxies, with reference to two independent group catalogues. We find that the fraction of these AGN in S0 host galaxies decreases strongly as a function of galaxy group size (halo mass) - which contrasts with the increasing fraction of galaxies of S0 type in denser environments. However, there is no evidence for an environmental dependence of AGN in spiral galaxies. Because most AGN are found in spiral galaxies, this dilutes the signature of environmental dependence for the population as a whole. We argue that the differing results for AGN in disk-dominated and bulge-dominated galaxies is related to the source of the gas fuelling the AGN, and so may also impact the luminosity function, duty cycle, and obscuration. We find that there is a significant difference in the luminosity function for AGN in spiral and S0 galaxies, and tentative evidence for some difference in the fraction of obscured AGN.

The VLTI instrument GRAVITY combines the beams from four telescopes and provides phase-referenced imaging as well as precision-astrometry of order 10 microarcseconds by observing two celestial objects in dual-field mode. Their angular separation can be determined from their differential OPD (dOPD) when the internal dOPDs in the interferometer are known. Here, we present the general overview of the novel metrology system which performs these measurements. The metrology consists of a three-beam laser system and a homodyne detection scheme for three-beam interference using phase-shifting interferometry in combination with lock-in amplifiers. Via this approach the metrology system measures dOPDs on a nanometer-level.

We use high spectral resolution (R > 8000) data covering 3800-13000\r{A} to study the physical conditions of the broad line region (BLR) of nine nearby Seyfert 1 galaxies. Up to six broad HI lines are present in each spectrum. A comparison - for the first time using simultaneous optical to near-infrared observations - to photoionisation calculations with our devised simple scheme yields the extinction to the BLR at the same time as determining the density and photon flux, and hence distance from the nucleus, of the emitting gas. This points to a typical density for the HI emitting gas of 10$^{11}$cm$^{-3}$ and shows that a significant amount of this gas lies at regions near the dust sublimation radius, consistent with theoretical predictions. We also confirm that in many objects the line ratios are far from case B, the best-fit intrinsic broad-line H$\alpha$/H$\beta$ ratios being in the range 2.5-6.6 as derived with our photoionization modeling scheme. The extinction to the BLR, based on independent estimates from HI and HeII lines, is A$_V$ $\le$ 3 for Seyfert 1-1.5s, while Seyfert 1.8-1.9s have A$_V$ in the range 4-8. A comparison of the extinction towards the BLR and narrow line region (NLR) indicates that the structure obscuring the BLR exists on scales smaller than the NLR. This could be the dusty torus, but dusty nuclear spirals or filaments could also be responsible. The ratios between the X-ray absorbing column N$_H$ and the extinction to the BLR are consistent with the Galactic gas-to-dust ratio if N$_H$ variations are considered.

MICADO will equip the E-ELT with a first light capability for diffraction limited imaging at near-infrared wavelengths. The instrument's observing modes focus on various flavours of imaging, including astrometric, high contrast, and time resolved. There is also a single object spectroscopic mode optimised for wavelength coverage at moderately high resolution. This contribution provides an overview of the key functionality of the instrument, outlining the scientific rationale for its observing modes. The interface between MICADO and the adaptive optics system MAORY that feeds it is summarised. The design of the instrument is discussed, focussing on the optics and mechanisms inside the cryostat, together with a brief overview of the other key sub-systems.

We present the results from a systematic search for galactic-scale, molecular (OH 119 $\mu$m) outflows in a sample of 52 Local Volume ($d < 50$ Mpc) Burst Alert Telescope detected active galactic nuclei (BAT AGN) with \emph{Herschel}-PACS. We combine the results from our analysis of the BAT AGN with the published \emph{Herschel}/PACS data of 43 nearby ($z<0.3$) galaxy mergers, mostly ultraluminous infrared galaxies (ULIRGs) and QSOs. The objects in our sample of BAT AGN have, on average, $\sim 10-100$ times lower AGN luminosities, star formation rates (SFRs), and stellar masses than those of the ULIRG and QSO sample. OH 119 $\mu$m is detected in 42 of our BAT AGN targets. Evidence for molecular outflows (i.e. OH absorption profiles with median velocities more blueshifted than $-$50 km s$^{-1}$ and/or blueshifted wings with 84-percentile velocities less than $-$300 km s$^{-1}$) is seen in only four BAT AGN (NGC~7479 is the most convincing case). Evidence for molecular inflows (i.e. OH absorption profiles with median velocities more redshifted than 50 km s$^{-1}$) is seen in seven objects, although an inverted P-Cygni profile is detected unambiguously in only one object (Circinus). Our data show that both the starburst and AGN contribute to driving OH outflows, but the fastest OH winds require AGN with quasar-like luminosities. We also confirm that the total absorption strength of OH 119 $\mu$m is a good proxy for dust optical depth as it correlates strongly with the 9.7 $\mu$m silicate absorption feature, a measure of obscuration originating in both the nuclear torus and host galaxy disk.

We use Herschel 70 to 160um images to study the size of the far-infrared emitting region in about 400 local galaxies and quasar (QSO) hosts. The sample includes normal `main-sequence' star-forming galaxies, as well as infrared luminous galaxies and Palomar-Green QSOs, with different levels and structures of star formation. Assuming Gaussian spatial distribution of the far-infrared (FIR) emission, the excellent stability of the Herschel point spread function (PSF) enables us to measure sizes well below the PSF width, by subtracting widths in quadrature. We derive scalings of FIR size and surface brightness of local galaxies with FIR luminosity, with distance from the star-forming main-sequence, and with FIR color. Luminosities L_FIR ~ 10^11 L_Sun can be reached with a variety of structures spanning 2 dex in size. Ultraluminous L_FIR >~ 10^12 L_Sun galaxies far above the main-sequence inevitably have small Re_70um ~ 0.5 kpc FIR emitting regions with large surface brightness, and can be close to optically thick in the FIR on average over these regions. Compared to these local relations, first ALMA sizes for the dust emission regions in high redshift galaxies, measured at somewhat longer rest wavelengths, suggest larger sizes at the same IR luminosity. We report a remarkably tight relation with 0.15 dex scatter between FIR surface brightness and the ratio of [CII] 158um emission and FIR emission -- the so-called [CII]-deficit is more tightly linked to surface brightness than to FIR luminosity or FIR color. Among 33 z <~ 0.1 PG QSOs with typical L_FIR/L_Bol,AGN ~ 0.1, 19 have a measured 70um half light radius, with median Re_70um = 1.1kpc. This is consistent with the FIR size for galaxies with similar L_FIR but lacking a QSO, in accordance with a scenario where the rest FIR emission of these types of QSOs is, in most cases, due to host star formation.

We present constraints on the molecular outflows in a sample of five Hyper-Luminous Infrared Galaxies using Herschel observations of the OH doublet at 119 {\mu}m. We have detected the OH doublet in three cases: one purely in emission and two purely in absorption. The observed emission profile has a significant blueshifted wing suggesting the possibility of tracing an outflow. Out of the two absorption profiles, one seems to be consistent with the systemic velocity while the other clearly indicates the presence of a molecular outflow whose maximum velocity is about ~1500 km/s. Our analysis shows that this system is in general agreement with previous results on Ultra-luminous Infrared Galaxies and QSOs, whose outflow velocities do not seem to correlate with stellar masses or starburst luminosities (star formation rates). Instead the galaxy outflow likely arises from an embedded AGN.

We report a tentative correlation between the outflow characteristics derived from OH absorption at $119\,\mu\text{m}$ and [CII] emission at $158\,\mu\text{m}$ in a sample of 22 local and bright ultraluminous infrared galaxies (ULIRGs). For this sample we investigate whether [CII] broad wings are a good tracer of molecular outflows, and how the two tracers are connected. Fourteen objects in our sample have a broad wing component as traced by [CII], and all of these also show OH119 absorption indicative of an outflow (in 1 case an inflow). The other eight cases, where no broad [CII] component was found, are predominantly objects with no OH outflow or a low-velocity ($\leq 100\,\text{km s}^{-1}$) OH outflow. The full width at half maximum (FWHM) of the broad [CII] component shows a trend with the OH119 blue-shifted velocity, although with significant scatter. Moreover, and despite large uncertainties, the outflow masses derived from OH and broad [CII] show a 1:1 relation. The main conclusion is therefore that broad [CII] wings can be used to trace molecular outflows. This may be particularly relevant at high redshift, where the usual tracers of molecular gas (like low-J CO lines) become hard to observe. Additionally, observations of blue-shifted Na I D $\lambda\lambda 5890,5896$ absorption are available for ten of our sources. Outflow velocities of Na I D show a trend with OH velocity and broad [CII] FWHM. These observations suggest that the atomic and molecular gas phases of the outflow are connected.

We discuss the dense molecular gas in central regions of nearby Seyfert galaxies, and report new arcsec resolution observations of HCN(1-0) and HCO$^+$(1-0) for 3 objects. In NGC 3079 the lines show complex profiles as a result of self-absorption and saturated continuum absorption. H$^{13}$CN reveals the continuum absorption profile, with a peak close to the galaxy's systemic velocity that traces disk rotation, and a second feature with a blue wing extending to $-350$km s$^{-1}$ that most likely traces a nuclear outflow. The morphological and spectral properties of the emission lines allow us to constrain the dense gas dynamics. We combine our kinematic analysis for these 3 objects, as well as another with archival data, with a previous comparable analysis of 4 other objects, to create a sample of 8 Seyferts. In 7 of these, the emission line kinematics imply thick disk structures on radial scales of $\sim$100pc, suggesting such structures are a common occurrence. We find a relation between the circumnuclear LHCN and Mdyn that can be explained by a gas fraction of 10% and a conversion factor {\alpha}HCN $\sim$ 10 between gas mass and HCN luminosity. Finally, adopting a different perspective to probe the physical properties of the gas around AGN, we report on an analysis of molecular line ratios which indicates that the clouds in this region are not self-gravitating.

The optical classification of a Seyfert galaxy and whether it is considered X-ray absorbed are often used interchangeably. But there are many borderline cases and also numerous examples where the optical and X-ray classifications appear to be in conflict. In this article we re-visit the relation between optical obscuration and X-ray absorption in AGNs. We make use of our "dust color" method (Burtscher et al. 2015) to derive the optical obscuration A_V and consistently estimated X-ray absorbing columns using 0.3--150 keV spectral energy distributions. We also take into account the variable nature of the neutral gas column N_H and derive the Seyfert sub-classes of all our objects in a consistent way. We show in a sample of 25 local, hard-X-ray detected Seyfert galaxies (log L_X / (erg/s) ~ 41.5 - 43.5) that there can actually be a good agreement between optical and X-ray classification. If Seyfert types 1.8 and 1.9 are considered unobscured, the threshold between X-ray unabsorbed and absorbed should be chosen at a column N_H = 10^22.3 / cm^2 to be consistent with the optical classification. We find that N_H is related to A_V and that the N_H/A_V ratio is approximately Galactic or higher in all sources, as indicated previously. But in several objects we also see that deviations from the Galactic ratio are only due to a variable X-ray column, showing that (1) deviations from the Galactic N_H/A_V can simply be explained by dust-free neutral gas within the broad line region in some sources, that (2) the dust properties in AGNs can be similar to Galactic dust and that (3) the dust color method is a robust way to estimate the optical extinction towards the sublimation radius in all but the most obscured AGNs.

Emission from high-J CO lines in galaxies has long been proposed as a tracer of X-ray dominated regions (XDRs) produced by AGN. Of particular interest is the question of whether the obscuring torus, which is required by AGN unification models, can be observed via high-J CO cooling lines. Here we report on the analysis of a deep Herschel-PACS observation of an extremely high J CO transition (40-39) in the Seyfert 2 galaxy NGC 1068. The line was not detected, with a derived 3$\sigma$ upper limit of $2 \times 10^{-17}\,\text{W}\,\text{m}^{-2}$. We apply an XDR model in order to investigate whether the upper limit constrains the properties of a molecular torus in NGC 1068. The XDR model predicts the CO Spectral Line Energy Distributions for various gas densities and illuminating X-ray fluxes. In our model, the CO(40-39) upper limit is matched by gas with densities $\sim 10^{6}-10^{7}\,\text{cm}^{-3}$, located at $1.6-5\,\text{pc}$ from the AGN, with column densities of at least $10^{25}\,\text{cm}^{-2}$. At such high column densities, however, dust absorbs most of the CO(40-39) line emission at $\lambda = 65.69\, \mu$m. Therefore, even if NGC 1068 has a molecular torus which radiates in the CO(40-39) line, the dust can attenuate the line emission to below the PACS detection limit. The upper limit is thus consistent with the existence of a molecular torus in NGC 1068. In general, we expect that the CO(40-39) is observable in only a few AGN nuclei (if at all), because of the required high gas column density, and absorption by dust.

We present high resolution (0."4) IRAM PdBI and ALMA mm and submm observations of the (ultra) luminous infrared galaxies ((U)LIRGs) IRAS17208-0014, Arp220, IC860 and Zw049.057 that reveal intense line emission from vibrationally excited ($\nu_2$=1) J=3-2 and 4-3 HCN. The emission is emerging from buried, compact (r<17-70 pc) nuclei that have very high implied mid-infrared surface brightness $>$$5\times 10^{13}$ L$_{\odot}$ kpc$^{-2}$. These nuclei are likely powered by accreting supermassive black holes (SMBHs) and/or hot (>200 K) extreme starbursts. Vibrational, $\nu_2$=1, lines of HCN are excited by intense 14 micron mid-infrared emission and are excellent probes of the dynamics, masses, and physical conditions of (U)LIRG nuclei when H$_2$ column densities exceed $10^{24}$ cm$^{-2}$. It is clear that these lines open up a new interesting avenue to gain access to the most obscured AGNs and starbursts. Vibrationally excited HCN acts as a proxy for the absorbed mid-infrared emission from the embedded nuclei, which allows for reconstruction of the intrinsic, hotter dust SED. In contrast, we show strong evidence that the ground vibrational state ($\nu$=0), J=3-2 and 4-3 rotational lines of HCN and HCO$^+$ fail to probe the highly enshrouded, compact nuclear regions owing to strong self- and continuum absorption. The HCN and HCO$^+$ line profiles are double-peaked because of the absorption and show evidence of non-circular motions - possibly in the form of in- or outflows. Detections of vibrationally excited HCN in external galaxies are so far limited to ULIRGs and early-type spiral LIRGs, and we discuss possible causes for this. We tentatively suggest that the peak of vibrationally excited HCN emission is connected to a rapid stage of nuclear growth, before the phase of strong feedback.

We describe a complete volume limited sample of nearby active galaxies selected by their 14-195keV luminosity, and outline its rationale for studying the mechanisms regulating gas inflow and outflow. We describe also a complementary sample of inactive galaxies, selected to match the AGN host galaxy properties. The active sample appears to have no bias in terms of AGN type, the only difference being the neutral absorbing column which is two orders of magnitude greater for the Seyfert 2s. In the luminosity range spanned by the sample, log L_{14-195keV} [erg/s] = 42.4-43.7, the optically obscured and X-ray absorbed fractions are 50-65%. The similarity of these fractions to more distant spectroscopic AGN samples, although over a limited luminosity range, suggests that the torus does not strongly evolve with redshift. Our sample confirms that X-ray unabsorbed Seyfert 2s are rare, comprising not more than a few percent of the Seyfert 2 population. At higher luminosities, the optically obscured fraction decreases (as expected for the increasing dust sublimation radius), but the X-ray absorbed fraction changes little. We argue that the cold X-ray absorption in these Seyfert 1s can be accounted for by neutral gas in clouds that also contribute to the broad line region (BLR) emission; and suggest that a geometrically thick neutral gas torus co-exists with the BLR and bridges the gap to the dusty torus.

We combine two approaches to isolate the AGN luminosity at near-infrared wavelengths and relate the near-IR pure AGN luminosity to other tracers of the AGN. Using integral-field spectroscopic data of an archival sample of 51 local AGNs, we estimate the fraction of non-stellar light by comparing the nuclear equivalent width of the stellar 2.3 micron CO absorption feature with the intrinsic value for each galaxy. We compare this fraction to that derived from a spectral decomposition of the integrated light in the central arc second and find them to be consistent with each other. Using our estimates of the near-IR AGN light, we find a strong correlation with presumably isotropic AGN tracers. We show that a significant offset exists between type 1 and type 2 sources in the sense that type 1 sources are 7 (10) times brighter in the near-IR at log L_MIR = 42.5 (log L_X = 42.5). These offsets only becomes clear when treating infrared type 1 sources as type 1 AGNs. All AGNs have very red near-to-mid-IR dust colors. This, as well as the range of observed near-IR temperatures, can be explained with a simple model with only two free parameters: the obscuration to the hot dust and the ratio between the warm and hot dust areas. We find obscurations of A_V (hot) = 5 - 15 mag for infrared type 1 sources and A_V (hot) = 15 - 35 mag for type 2 sources. The ratio of hot dust to warm dust areas of about 1000 is nicely consistent with the ratio of radii of the respective regions as found by infrared interferometry.

GRAVITY is the second generation Very Large Telescope Interferometer instrument for precision narrow-angle astrometry and interferometric imaging in the Near Infra-Red (NIR). It shall provide precision astrometry of order 10 microarcseconds, and imaging capability at a few milliarcsecond resolution, and hence will revolutionise dynamical measurements of celestial objects. GRAVITY is currently in the last stages of its integration and tests in Garching at MPE, and will be delivered to the VLT Interferometer (VLTI) in 2015. We present here the instrument, with a particular focus on the components making use of fibres: integrated optics beam combiners, polarisation rotators, fibre differential delay lines, and the metrology.

GRAVITY is the four-beam, near- infrared, AO-assisted, fringe tracking, astrometric and imaging instrument for the Very Large Telescope Interferometer (VLTI). It is requiring the development of one of the most complex instrument software systems ever built for an ESO instrument. Apart from its many interfaces and interdependencies, one of the most challenging aspects is the overall performance and stability of this complex system. The three infrared detectors and the fast reflective memory network (RMN) recorder contribute a total data rate of up to 20 MiB/s accumulating to a maximum of 250 GiB of data per night. The detectors, the two instrument Local Control Units (LCUs) as well as the five LCUs running applications under TAC (Tools for Advanced Control) architecture, are interconnected with fast Ethernet, RMN fibers and dedicated fiber connections as well as signals for the time synchronization. Here we give a simplified overview of all subsystems of GRAVITY and their interfaces and discuss two examples of high-level applications during observations: the acquisition procedure and the gathering and merging of data to the final FITS file.

We report the detection of far-infrared (FIR) CO rotational emission from nearby active galactic nuclei (AGN) and starburst galaxies, as well as several merging systems and Ultra-Luminous Infrared Galaxies (ULIRGs). Using Herschel-PACS, we have detected transitions in the J$_{upp}$ = 14 - 20 range ($\lambda \sim$ 130 - 185 $\mu$m, $\nu \sim$ 1612 - 2300 GHz) with upper limits on (and in two cases, detections of) CO line fluxes up to J$_{upp}$ = 30. The PACS CO data obtained here provide the first well-sampled FIR extragalactic CO SLEDs for this range, and will be an essential reference for future high redshift studies. We find a large range in the overall SLED shape, even amongst galaxies of similar type, demonstrating the uncertainties in relying solely on high-J CO diagnostics to characterize the excitation source of a galaxy. Combining our data with low-J line intensities taken from the literature, we present a CO ratio-ratio diagram and discuss its potential diagnostic value in distinguishing excitation sources and physical properties of the molecular gas. The position of a galaxy on such a diagram is less a signature of its excitation mechanism, than an indicator of the presence (or absence) of warm, dense molecular gas. We then quantitatively analyze the CO emission from a subset of the detected sources with Large Velocity Gradient (LVG) radiative transfer models to fit the CO SLEDs. Using both single-component and two-component LVG models to fit the kinetic temperature, velocity gradient, number density and column density of the gas, we derive the molecular gas mass and the corresponding CO-to-H$_2$ conversion factor, $\alpha_{CO}$, for each respective source. For the ULIRGs we find $\alpha$ values in the canonical range 0.4 - 5 M$_\odot$/(K kms$^{-1}$pc$^2$), while for the other objects, $\alpha$ varies between 0.2 and 14.} Finally, we compare our best-fit LVG model ..

The VLTI instrument GRAVITY will provide very powerful astrometry by combining the light from four telescopes for two objects simultaneously. It will measure the angular separation between the two astronomical objects to a precision of 10 microarcseconds. This corresponds to a differential optical path difference (dOPD) between the targets of few nanometers and the paths within the interferometer have to be maintained stable to that level. For this purpose, the novel metrology system of GRAVITY will monitor the internal dOPDs by means of phase-shifting interferometry. We present the four-step phase-shifting concept of the metrology with emphasis on the method used for calibrating the phase shifts. The latter is based on a phase-step insensitive algorithm which unambiguously extracts phases in contrast to other methods that are strongly limited by non-linearities of the phase-shifting device. The main constraint of this algorithm is to introduce a robust ellipse fitting routine. Via this approach we are able to measure phase shifts in the laboratory with a typical accuracy of lambda/2000 or 1 nanometer of the metrology wavelength.

Herschel/PACS observations of 29 local (Ultra-)Luminous Infrared Galaxies, including both starburst and AGN-dominated sources as diagnosed in the mid-infrared/optical, show that the equivalent width of the absorbing OH 65 um Pi_{3/2} J=9/2-7/2 line (W_{eq}(OH65)) with lower level energy E_{low}~300 K, is anticorrelated with the [C ii]158 um line to far-infrared luminosity ratio, and correlated with the far-infrared luminosity per unit gas mass and with the 60-to-100 um far-infrared color. While all sources are in the active L_{IR}/M_{H2}>50 Lsun/Msun mode as derived from previous CO line studies, the OH65 absorption shows a bimodal distribution with a discontinuity at L_{FIR}/M_{H2}~100 Lsun/Msun. In the most buried sources, OH65 probes material partially responsible for the silicate 9.7 um absorption. Combined with observations of the OH 71 um Pi_{1/2} J=7/2-5/2 doublet (E_{low}~415 K), radiative transfer models characterized by the equivalent dust temperature, Tdust, and the continuum optical depth at 100 um, tau_{100}, indicate that strong [C ii]158 um deficits are associated with far-IR thick (tau_{100}>~0.7, N_{H}>~10^{24} cm^{-2}), warm (Tdust>~60 K) structures where the OH 65 um absorption is produced, most likely in circumnuclear disks/tori/cocoons. With their high L_{FIR}/M_{H2} ratios and columns, the presence of these structures is expected to give rise to strong [C ii] deficits. W_{eq}(OH65) probes the fraction of infrared luminosity arising from these compact/warm environments, which is >~30-50% in sources with high W_{eq}({OH65}). Sources with high W_{eq}({OH65}) have surface densities of both L_{IR} and M_{H2} higher than inferred from the half-light (CO or UV/optical) radius, tracing coherent structures that represent the most buried/active stage of (circum)nuclear starburst-AGN co-evolution.

We present the results of interferometric spectral line observations of Arp 220 at 3.5mm and 1.2mm from the Plateau de Bure Interferometer (PdBI), imaging the two nuclear disks in H$^{13}$CN$(1 - 0)$ and $(3 - 2)$, H$^{13}$CO$^+(1 - 0)$ and $(3 - 2)$, and HN$^{13}$C$(3 - 2)$ as well as SiO$(2 - 1)$ and $(6 - 5)$, HC$^{15}$N$(3 - 2)$, and SO$(6_6 - 5_5)$. The gas traced by SiO$(6 - 5)$ has a complex and extended kinematic signature including a prominent P Cygni profile, almost identical to previous observations of HCO$^+(3 - 2)$. Spatial offsets $0.1''$ north and south of the continuum centre in the emission and absorption of the SiO$(6 - 5)$ P Cygni profile in the western nucleus (WN) imply a bipolar outflow, delineating the northern and southern edges of its disk and suggesting a disk radius of $\sim40$ pc, consistent with that found by ALMA observations of Arp 220. We address the blending of SiO$(6 - 5)$ and H$^{13}$CO$^+(3 - 2)$ by considering two limiting cases with regards to the H$^{13}$CO$^+$ emission throughout our analysis. Large velocity gradient (LVG) modelling is used to constrain the physical conditions of the gas and to infer abundance ratios in the two nuclei. Our most conservative lower limit on the [H$^{13}$CN]/[H$^{13}$CO$^+$] abundance ratio is 11 in the WN, cf. 0.10 in the eastern nucleus (EN). Comparing these ratios to the literature we argue on chemical grounds for an energetically significant AGN in the WN driving either X-ray or shock chemistry, and a dominant starburst in the EN.

We present the installed and fully operational beam stabilization and fiber injection subsystem feeding the 2nd generation VLTI instrument GRAVITY. The interferometer GRAVITY requires an unprecedented stability of the VLTI optical train to achieve micro-arcsecond astrometry. For this purpose, GRAVITY contains four fiber coupler units, one per telescope. Each unit is equipped with actuators to stabilize the telescope beam in terms of tilt and lateral pupil displacement, to rotate the field, to adjust the polarization and to compensate atmospheric piston. A special roof-prism offers the possibility of on-axis as well as off-axis fringe tracking without changing the optical train. We describe the assembly, integration and alignment and the resulting optical quality and performance of the individual units. Finally, we present the closed-loop performance of the tip-tilt and pupil tracking achieved with the final systems in the lab.

Past observations of QSO host galaxies at z >6 have found cold gas and star formation on compact scales of a few kiloparsecs. We present new high sensitivity IRAM PdBI follow-up observations of the [CII] 158micron emission line and FIR continuum in the host galaxy of SDSS J1148+5152, a luminous QSO at redshift 6.4189. We find that a large fraction of the gas traced by [CII] is at high velocities, up to ~1400 km/s relative to the systemic velocity, confirming the presence of a major quasar-driven outflow indicated by previous observations. The outflow has a complex morphology and reaches a maximum projected radius of ~30 kpc. The extreme spatial extent of the outflow allows us, for the first time in an external galaxy, to estimate mass-loss rate, kinetic power and momentum rate of the outflow as a function of the projected distance from the nucleus and the dynamical time-scale. These trends reveal multiple outflow events during the past 100 Myr, although the bulk of the mass, energy and momentum appear to have been released more recently, within the past ~20 Myr. Surprisingly, we discover that also the quiescent gas at systemic velocity is extremely extended. More specifically, we find that, while 30% of the [CII] within v\in(-200, 200) km/s traces a compact component that is not resolved by our observations, 70% of the [CII] emission in this velocity range is extended, with a projected FWHM size of 17.4+-1.4 kpc. We detect FIR continuum emission associated with both the compact and the extended [CII] components, although the extended FIR emission has a FWHM of 11+-3 kpc, thus smaller than the extended [CII] source. Overall, our results indicate that the cold gas traced by [CII] is distributed up to r~30 kpc. A large fraction of extended [CII] is likely associated with star formation on large scales, but the [CII] source extends well beyond the FIR continuum.

GRAVITY is the second generation VLT Interferometer (VLTI) instrument for high-precision narrow-angle astrometry and phase-referenced interferometric imaging. The laser metrology system of GRAVITY is at the heart of its astrometric mode, which must measure the distance of 2 stars with a precision of 10 micro-arcseconds. This means the metrology has to measure the optical path difference between the two beam combiners of GRAVITY to a level of 5 nm. The metrology design presents some non-common paths that have consequently to be stable at a level of 1 nm. Otherwise they would impact the performance of GRAVITY. The various tests we made in the past on the prototype give us hints on the components responsible for this error, and on their respective contribution to the total error. It is however difficult to assess their exact origin from only OPD measurements, and therefore, to propose a solution to this problem. In this paper, we present the results of a semi-empirical modeling of the fibered metrology system, relying on theoretical basis, as well as on characterisations of key components. The modeling of the metrology system regarding various effects, e.g., temperature, waveguide heating or mechanical stress, will help us to understand how the metrology behave. The goals of this modeling are to 1) model the test set-ups and reproduce the measurements (as a validation of the modeling), 2) determine the origin of the non-common path errors, and 3) propose modifications to the current metrology design to reach the required 1nm stability.

We present in this paper the design and characterisation of a new sub-system of the VLTI 2nd generation instrument GRAVITY: the Calibration Unit. The Calibration Unit provides all functions to test and calibrate the beam combiner instrument: it creates two artificial stars on four beams, and dispose of four delay lines with an internal metrology. It also includes artificial stars for the tip-tilt and pupil guiding systems, as well as four metrology pick-up diodes, for tests and calibration of the corresponding sub-systems. The calibration unit also hosts the reference targets to align GRAVITY to the VLTI, and the safety shutters to avoid the metrology light to propagate in the VLTI-lab. We present the results of the characterisation and validtion of these differrent sub-units.

A deep 400-ksec ACIS-S observation of the nearest quasar known, Mrk 231, is combined with archival 120-ksec data obtained with the same instrument and setup to carry out the first ever spatially resolved spectral analysis of a hot X-ray emitting circumgalactic nebula around a quasar. The 65 x 50 kpc X-ray nebula shares no resemblance with the tidal debris seen at optical wavelengths. One notable exception is the small tidal arc 3.5 kpc south of the nucleus where excess soft X-ray continuum emission and Si XIII 1.8 keV line emission are detected, consistent with star formation and its associated alpha-element enhancement, respectively. An X-ray shadow is also detected at the location of the 15-kpc northern tidal tail. The hard X-ray continuum emission within 6 kpc of the center is consistent with being due entirely to the bright central AGN. The soft X-ray spectrum of the outer (>6 kpc) portion of the nebula is best described as the sum of two thermal components with T~3 and ~8 million K and spatially uniform super-solar alpha element abundances, relative to iron. This result implies enhanced star formation activity over ~10^8 yrs accompanied with redistribution of the metals on large scale. The low-temperature thermal component is not present within 6 kpc of the nucleus, suggesting extra heating in this region from the circumnuclear starburst, the central quasar, or the wide-angle quasar-driven outflow identified from optical IFU spectroscopy on a scale of >3 kpc. Significant azimuthal variations in the soft X-ray intensity are detected in the inner region where the outflow is present. The soft X-ray emission is weaker in the western quadrant, coincident with a deficit of Halpha and some of the largest columns of neutral gas outflowing from the nucleus. Shocks created by the interaction of the wind with the ambient ISM may heat the gas to high temperatures at this location. (abridged)

We analyze the applicability of far-infrared fine-structure lines [CII] 158 micron, [OI] 63 micron and [OIII] 88 micron to reliably trace the star formation rate (SFR) in a sample of low-metallicity dwarf galaxies from the Herschel Dwarf Galaxy Survey and compare with a broad sample of galaxies of various types and metallicities in the literature. We study the trends and scatter in the relation between the SFR (as traced by GALEX FUV and MIPS 24 micron) and far-infrared line emission, on spatially resolved and global galaxy scales, in dwarf galaxies. We assemble far-infrared line measurements from the literature and infer whether the far-infrared lines can probe the SFR (as traced by the total-infrared luminosity) in a variety of galaxy populations. In metal-poor dwarfs, the [OI] and [OIII] lines show the strongest correlation with the SFR with an uncertainty on the SFR estimates better than a factor of 2, while the link between [CII] emission and the SFR is more dispersed (uncertainty factor of 2.6). The increased scatter in the SFR-L([CII]) relation towards low metal abundances, warm dust temperatures, large filling factors of diffuse, highly ionized gas suggests that other cooling lines start to dominate depending on the density and ionization state of the gas. For the literature sample, we evaluate the correlations for a number of different galaxy populations. The [CII] and [OI] lines are considered to be reliable SFR tracers in starburst galaxies, recovering the star formation activity within an uncertainty of factor 2. [Abridged]

[Abridged] We study the evolution of the dust temperatures of galaxies in the SFR-M* plane up to z~2 using observations from the Herschel Space Observatory. Starting from a sample of galaxies with reliable star-formation rates (SFRs), stellar masses (M*) and redshift estimates, we grid the SFR-M* parameter space in several redshift ranges and estimate the mean Tdust of each SFR-M*-z bin. Dust temperatures are inferred using the stacked far-infrared flux densities of our SFR-M*-z bins. At all redshifts, Tdust increases with infrared luminosities (LIR), specific SFRs (SSFR; i.e., SFR/M*) and distances with respect to the main sequence (MS) of the SFR-M* plane (i.e., D_SSFR_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). The Tdust-SSFR and Tdust-D_SSFR_MS correlations are statistically more significant than the Tdust-LIR one. While the slopes of these three correlations are redshift-independent, their normalizations evolve from z=0 and z~2. We convert these results into a recipe to derive Tdust from SFR, M* and z. The existence of a strong Tdust-D_SSFR_MS correlation provides us with information on the dust and gas content of galaxies. (i) The slope of the Tdust-D__SSFR_MS correlation can be explained by the increase of the star-formation efficiency (SFE; SFR/Mgas) with D_SSFR_MS as found locally by molecular gas studies. (ii) At fixed D_SSFR_MS, the constant Tdust observed in galaxies probing large ranges in SFR and M* can be explained by an increase or decrease of the number of star-forming regions with comparable SFE enclosed in them. (iii) At high redshift, the normalization towards hotter temperature of the Tdust-D_SSFR_MS correlation can be explained by the decrease of the metallicities of galaxies or by the increase of the SFE of MS galaxies. All these results support the hypothesis that the conditions prevailing in the star-forming regions of MS and far-above-MS galaxies are different.

We study the properties of massive, galactic-scale outflows of molecular gas and investigate their impact on galaxy evolution. We present new IRAM PdBI CO(1-0) observations of local ULIRGs and QSO hosts: clear signature of massive and energetic molecular outflows, extending on kpc scales, is found in the CO(1-0) kinematics of four out of seven sources, with measured outflow rates of several 100 Msun/yr. We combine these new observations with data from the literature, and explore the nature and origin of massive molecular outflows within an extended sample of 19 local galaxies. We find that starburst-dominated galaxies have an outflow rate comparable to their SFR, or even higher by a factor of ~ 2-4, implying that starbursts can indeed be effective in removing cold gas from galaxies. Nevertheless, our results suggest that the presence of an AGN can boost the outflow rate by a large factor, which is found to increase with the L_AGN/L_bol ratio. The gas depletion time-scales due to molecular outflows are anti-correlated with the presence and luminosity of an AGN in these galaxies, and range from a few hundred million years in starburst galaxies, down to just a few million years in galaxies hosting powerful AGNs. In quasar hosts the depletion time-scales due to the outflow are much shorter than the depletion time-scales due to star formation. We estimate the outflow kinetic power and find that, for galaxies hosting powerful AGNs, it corresponds to about 5% of the AGN luminosity, as expected by models of AGN feedback. Moreover, we find that momentum rates of about 20 L_AGN/c are common among the AGN-dominated sources in our sample. For "pure" starburst galaxies our data tentatively support models in which outflows are mostly momentum-driven by the radiation pressure from young stars onto dusty clouds.

We present a study of the morphology and intensity of star formation in the host galaxies of eight Palomar-Green quasars using observations with the Hubble Space Telescope. Our observations are motivated by recent evidence for a close relationship between black hole growth and the stellar mass evolution in its host galaxy. We use narrow-band [O II] $\lambda$3727, H$\beta$, [O III] $\lambda$5007 and Pa$\alpha$ images, taken with the WFPC2 and NICMOS instruments, to map the morphology of line-emitting regions, and, after extinction corrections, diagnose the excitation mechanism and infer star-formation rates. Significant challenges in this type of work are the separation of the quasar light from the stellar continuum and the quasar-excited gas from the star-forming regions. To this end, we present a novel technique for image decomposition and subtraction of quasar light. Our primary result is the detection of extended line-emitting regions with sizes ranging from 0.5 to 5 kpc and distributed symmetrically around the nucleus, powered primarily by star formation. We determine star-formation rates of order a few tens of M$_\odot$/yr. The host galaxies of our target quasars have stellar masses of order $10^{11}$ M$_\odot$ and specific star formation rates on a par with those of M82 and luminous infrared galaxies. As such they fall at the upper envelope or just above the star-formation mass sequence in the specific star formation vs stellar mass diagram. We see a clear trend of increasing star formation rate with quasar luminosity, reinforcing the link between the growth of the stellar mass of the host and the black hole mass found by other authors.

We report on the Herschel/PACS observations of OH in Mrk 231, with detections in 9 doublets observed within the PACS range, and present radiative transfer models for the outflowing OH. Signatures of outflowing gas are found in up to 6 OH doublets with different excitation requirements. At least two outflowing components are identified, one with OH radiatively excited, and the other with low excitation, presumably spatially extended. Particularly prominent, the blue wing of the absorption detected in the in-ladder 2Pi_{3/2} J=9/2-7/2 OH doublet at 65 um, with E_lower=290 K, indicates that the excited outflowing gas is generated in a compact and warm (circum)nuclear region. Because the excited, outflowing OH gas in Mrk 231 is associated with the warm, far-IR continuum source, it is likely more compact (diameter of 200-300 pc) than that probed by CO and HCN. Nevertheless, its mass-outflow rate per unit of solid angle as inferred from OH is similar to that previously derived from CO, >~70x(2.5x10^{-6}/X_{OH}) Msun yr^{-1} sr^{-1}, where X_{OH} is the OH abundance relative to H nuclei. In spherical symmetry, this would correspond to >~850x(2.5x10^{-6}/X_{OH}) Msun yr^{-1}, though significant collimation is inferred from the line profiles. The momentum flux of the excited component attains ~15 L_{AGN}/c, with an OH column density of (1.5-3)x10^{17} cm^-2 and a mechanical luminosity of ~10^{11} Lsun. The detection of very excited OH peaking at central velocities indicates the presence of a nuclear reservoir of gas rich in OH, plausibly the 130-pc scale circumnuclear torus previously detected in OH megamaser emission, that may be feeding the outflow. An exceptional ^{18}OH enhancement, with OH/^{18}OH<~30 at both central and blueshifted velocities, is likely the result of interstellar-medium processing by recent starburst/SNe activity.

We present Herschel photometry and spectroscopy, carried out as part of the Herschel ULIRG survey (HERUS), and a model for the infrared to submillimetre emission of the ultraluminous infrared galaxy IRAS 08572+3915. This source shows one of the deepest known silicate absorption features and no polycyclic aromatic hydrocarbon (PAH) emission. The model suggests that this object is powered by an active galactic nucleus (AGN) with a fairly smooth torus viewed almost edge-on and a very young starburst. According to our model the AGN contributes about 90% of the total luminosity of 1.1 x 10^13 Lo, which is about a factor of five higher than previous estimates. The large correction of the luminosity is due to the anisotropy of the emission of the best fit torus. Similar corrections may be necessary for other local and high-z analogs. This correction implies that IRAS 08572+3915 at a redshift of 0.05835 may be the nearest hyperluminous infrared galaxy and probably the most luminous infrared galaxy in the local (z < 0.2) Universe. IRAS 08572+3915 shows a low ratio of [CII] to IR luminosity (log L_[CII]/L_{IR} < -3.8) and a [OI]63um to [CII]158um line ratio of about 1 that supports the model presented in this letter.

We combine IRAM Plateau de Bure Interferometer and Herschel PACS and SPIRE measurements to study the dust and gas contents of high-redshift star forming galaxies. We present new observations for a sample of 17 lensed galaxies at z=1.4-3.1, which allow us to directly probe the cold ISM of normal star-forming galaxies with stellar masses of ~10^10Msun, a regime otherwise not (yet) accessible by individual detections in Herschel and molecular gas studies. The lensed galaxies are combined with reference samples of sub-millimeter and normal z~1-2 star-forming galaxies with similar far-infrared photometry to study the gas and dust properties of galaxies in the SFR-M*-redshift parameter space. The mean gas depletion timescale of main sequence galaxies at z>2 is measured to be only ~450Myr, a factor of ~1.5 (~5) shorter than at z=1 (z=0), in agreement with a (1+z)^-1 scaling. The mean gas mass fraction at z=2.8 is 40+/-15% (44% after incompleteness correction), suggesting a flattening or even a reversal of the trend of increasing gas fractions with redshift recently observed up to z~2. The depletion timescale and gas fractions of the z>2 normal star-forming galaxies can be explained under the "equilibrium model" for galaxy evolution, in which the gas reservoir of galaxies is the primary driver of the redshift evolution of specific star formation rates. Due to their high star formation efficiencies and low metallicities, the z>2 lensed galaxies have warm dust despite being located on the star formation main sequence. At fixed metallicity, they also have a gas-to-dust ratio 1.7 times larger than observed locally when using the same standard techniques, suggesting that applying the local calibration of the relation between gas-to-dust ratio and metallicity to infer the molecular gas mass of high redshift galaxies may lead to systematic differences with CO-based estimates.

We present the first Herschel spectroscopic detections of the [OI]63 and [CII]158 micron fine-structure transitions, and a single para-H2O line from the 35 x 15 kpc^2 shocked intergalactic filament in Stephan's Quintet. The filament is believed to have been formed when a high-speed intruder to the group collided with clumpy intergroup gas. Observations with the PACS spectrometer provide evidence for broad (> 1000 km s^-1) luminous [CII] line profiles, as well as fainter [OI]63micron emission. SPIRE FTS observations reveal water emission from the p-H2O (111-000) transition at several positions in the filament, but no other molecular lines. The H2O line is narrow, and may be associated with denser intermediate-velocity gas experiencing the strongest shock-heating. The [CII]/PAH{tot) and [CII]/FIR ratios are too large to be explained by normal photo-electric heating in PDRs. HII region excitation or X-ray/Cosmic Ray heating can also be ruled out. The observations lead to the conclusion that a large fraction the molecular gas is diffuse and warm. We propose that the [CII], [OI] and warm H2 line emission is powered by a turbulent cascade in which kinetic energy from the galaxy collision with the IGM is dissipated to small scales and low-velocities, via shocks and turbulent eddies. Low-velocity magnetic shocks can help explain both the [CII]/[OI] ratio, and the relatively high [CII]/H2 ratios observed. The discovery that [CII] emission can be enhanced, in large-scale turbulent regions in collisional environments has implications for the interpretation of [CII] emission in high-z galaxies.

We report the results from a systematic search for molecular (OH-119 um) outflows with Herschel-PACS in a sample of 43 nearby (z < 0.3) galaxy mergers, mostly ultraluminous infrared galaxies (ULIRGs) and QSOs. We find that the character of the OH feature (strength of the absorption relative to the emission) correlates with that of the 9.7-um silicate feature, a measure of obscuration in ULIRGs. Unambiguous evidence for molecular outflows, based on the detection of OH absorption profiles with median velocities more blueshifted than -50 km/sec, is seen in 26 (70%) of the 37 OH-detected targets, suggesting a wide-angle (~145 degrees) outflow geometry. Conversely, unambiguous evidence for molecular inflows, based on the detection of OH absorption profiles with median velocities more redshifted than +50 km/sec, is seen in only 4 objects, suggesting a planar or filamentary geometry for the inflowing gas. Terminal outflow velocities of ~-1000 km/sec are measured in several objects, but median outflow velocities are typically ~-200 km s^{-1}. While the outflow velocities show no statistically significant dependence on the star formation rate, they are distinctly more blueshifted among systems with large AGN fractions and luminosities [log (L_AGN / L_sun) > 11.8 +/- 0.3]. The quasars in these systems play a dominant role in driving the molecular outflows. In contrast, the most AGN dominated systems, where OH is seen purely in emission, show relatively modest OH line widths, despite their large AGN luminosities, perhaps indicating that molecular outflows subside once the quasar has cleared a path through the obscuring material.

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

We report the PACS-100um/160um detections of a sample of 42 GALEX-selected and FIR-detected Lyman break galaxies (LBGs) at z ~ 1 located in the COSMOS field and analyze their ultra-violet (UV) to far-infrared (FIR) properties. The detection of these LBGs in the FIR indicates that they have a dust content high enough so that its emission can be directly detected. According to a spectral energy distribution (SED) fitting with stellar population templates to their UV-to-near-IR observed photometry, PACS-detected LBGs tend to be bigger, more massive, dustier, redder in the UV continuum, and UV-brighter than PACS-undetected LBGs. PACS-detected LBGs at z ~ 1 are mostly disk-like galaxies and are located over the green-valley and red sequence of the color-magnitude diagram of galaxies at their redshift. By using their UV and IR emission, we find that PACS-detected LBGs tend to be less dusty and have slightly higher total star-formation rates (SFRs) than other PACS-detected UV-selected galaxies within their same redshift range. As a consequence of the selection effect due to the depth of the FIR observations employed, all our PACS-detected LBGs are LIRGs. However, none of them are in the ULIRG regime, where the FIR observations are complete. The finding of ULIRGs-LBGs at higher redshifts suggests an evolution of the FIR emission of LBGs with cosmic time. In an IRX-$\beta$ diagram, PACS-detected LBGs at z ~ 1 tend to be located around the relation for local starburst similarly to other UV-selected PACS-detected galaxies at their same redshift. Consequently, the dust-correction factors obtained with their UV continuum slope allow to determine their total SFR, unlike at higher redshifts. However, the dust attenuation derived from UV to NIR SED fitting overestimates the total SFR for most of our PACS-detected LBGs in age-dependent way: the overestimation factor is higher in younger galaxies.

We present the first results of a survey of the [CII]158um emission line in 241 luminous infrared galaxies (LIRGs) comprising the Great Observatories All-sky Survey (GOALS) sample, obtained with the PACS instrument on board Herschel. The [CII] luminosities of the LIRGs in GOALS range from ~10^7 to 2x10^9 Lsun. We find that LIRGs show a tight correlation of [CII]/FIR with far-IR flux density ratios, with a strong negative trend spanning from ~10^-2 to 10^-4, as the average temperature of dust increases. We find correlations between the [CII]/FIR ratio and the strength of the 9.7um silicate absorption feature as well as with the luminosity surface density of the mid-IR emitting region (Sigma_MIR), suggesting that warmer, more compact starbursts have substantially smaller [CII]/FIR ratios. Pure star-forming (SF) LIRGs have a mean [CII]/FIR ~ 4x10^-3, while galaxies with low 6.2um PAH equivalent widths (EWs), indicative of the presence of active galactic nuclei (AGN), span the full range in [CII]/FIR. However, we show that even when only pure SF galaxies are considered, the [CII]/FIR ratio drops by an order of magnitude, from 10^-2 to 10^-3, with Sigma_MIR and Sigma_IR, implying that the [CII] luminosity is not a good indicator of the star formation rate (SFR) for most LIRGs, for it does not scale linearly with the warm dust emission. Moreover, even in LIRGs in which we detect an AGN in the mid-IR, the majority (2/3) of galaxies show [CII]/FIR >= 10^-3 typical of high 6.2um PAH EW sources, suggesting that most AGNs do not contribute significantly to the far-IR emission. We provide an empirical relation between the [CII]/FIR and the specific SFR (SSFR) for SF LIRGs. Finally, we present predictions for the starburst size based on the observed [CII] and far-IR luminosities which should be useful for comparing with results from future surveys of high-redshift galaxies with ALMA and CCAT.

We present high spatial resolution CO(1-0) mapping of the luminous infrared galaxy NGC6240 obtained with the IRAM - Plateau de Bure Interferometer. This source is a well-known early-stage merging system hosting two AGN. We find a broad CO(1-0) line profile, with maximum velocity 800 km/s and FWZI=1400 km/s, and displaying several kinematic components, revealing the complexity of the gas dynamics in this system. We detect a blueshifted CO emission with velocity in the range -200 and -500 km/s, peaked around the southern AGN, at the same position where the H2 outflow is located, and with a mass loss rate of ~500 Msun/yr. We interpret this blueshifted component as a outflow, originating from the southern nucleus. The spatial and spectral match strongly suggests that the CO outflow is connected to the H2 superwind located around the southern nucleus, and to the large scale CO outflow, with similar velocities, extended on scales of 10 kpc. The large mass loading factor (dM/dt /SFR ~ 10) of the molecular gas suggests that the outflow is likely driven by both SNa winds and the radiation of the southern AGN. We discovered a nuclear, redshifted CO emission peaking in the mid point of the two nuclei, as it is the case for the CO emission at the systemic velocity. The large velocity dispersion, which reaches its maximum (~500 Km/s) in the mid point between the two nuclei, suggests that the gas might be highly turbulent in this region, although the presence of an unresolved rotation component cannot be ruled out.

The Dwarf Galaxy Survey (DGS) program is studying low-metallicity galaxies using 230h of far-infrared (FIR) and submillimetre (submm) photometric and spectroscopic observations of the Herschel Space Observatory and draws to this a rich database of a wide range of wavelengths tracing the dust, gas and stars. This sample of 50 galaxies includes the largest metallicity range achievable in the local Universe including the lowest metallicity (Z) galaxies, 1/50 Zsun, and spans 4 orders of magnitude in star formation rates. The survey is designed to get a handle on the physics of the interstellar medium (ISM) of low metallicity dwarf galaxies, especially on their dust and gas properties and the ISM heating and cooling processes. The DGS produces PACS and SPIRE maps of low-metallicity galaxies observed at 70, 100, 160, 250, 350, and 500 mic with the highest sensitivity achievable to date in the FIR and submm. The FIR fine-structure lines, [CII] 158 mic, [OI] 63 mic, [OI] 145 mic, [OIII] 88 mic, [NIII] 57 mic and [NII] 122 and 205 mic have also been observed with the aim of studying the gas cooling in the neutral and ionized phases. The SPIRE FTS observations include many CO lines (J=4-3 to J=13-12), [NII] 205 mic and [CI] lines at 370 and 609 mic. This paper describes the sample selection and global properties of the galaxies, the observing strategy as well as the vast ancillary database available to complement the Herschel observations. The scientific potential of the full DGS survey is described with some example results included.

By making use of Herschel-PEP observations of the COSMOS and Extended Groth Strip fields, we have estimated the dependence of the clustering properties of FIR-selected sources on their 100um fluxes. Our analysis shows a tendency for the clustering strength to decrease with limiting fluxes: r0(S100um >8 mJy)~4.3 Mpc and r0(S100um >5 mJy)~5.8 Mpc. These values convert into minimum halo masses Mmin~10^{11.6} Msun for sources brighter than 8 mJy and Mmin~10^{12.4} Msun for S100um > 5 mJy galaxies. We show such an increase of the clustering strength to be due to an intervening population of z~2 sources, which are very strongly clustered and whose relative contribution, equal to about 10% of the total counts at S100um > 2 mJy, rapidly decreases for brighter flux cuts. By removing such a contribution, we find that z <~ 1 FIR galaxies have approximately the same clustering properties, irrespective of their flux level. The above results were then used to investigate the intrinsic dependence on cosmic epoch of the clustering strength of dusty star-forming galaxies between z~0 and z~2.5. In order to remove any bias in the selection process, the adopted sample only includes galaxies observed at the same rest-frame wavelength, lambda~60 um, which have comparable luminosities and therefore star-formation rates (SFR>~100 Msun/yr). Our analysis shows that the same amount of (intense) star forming activity takes place in extremely different environments at the different cosmological epochs. For z<~1 the hosts of such star forming systems are small, Mmin~10^{11} Msun, isolated galaxies. High (z~2) redshift star formation instead seems to uniquely take place in extremely massive/cluster-like halos, Mmin~10^{13.5} Msun, which are associated with the highest peaks of the density fluctuation field at those epochs. (abridged)

We present results from the deepest Herschel-PACS (Photodetector Array Camera and Spectrometer) far-infrared blank field extragalactic survey, obtained by combining observations of the GOODS (Great Observatories Origins Deep Survey) fields from the PACS Evolutionary Probe (PEP) and GOODS-Herschel key programmes. We describe data reduction and the construction of images and catalogues. In the deepest parts of the GOODS-S field, the catalogues reach 3-sigma depths of 0.9, 0.6 and 1.3 mJy at 70, 100 and 160 um, respectively, and resolve ~75% of the cosmic infrared background at 100um and 160um into individually detected sources. We use these data to estimate the PACS confusion noise, to derive the PACS number counts down to unprecedented depths and to determine the infrared luminosity function of galaxies down to LIR=10^11 Lsun at z~1 and LIR=10^12 Lsun at z~2, respectively. For the infrared luminosity function of galaxies, our deep Herschel far-infrared observations are fundamental because they provide more accurate infrared luminosity estimates than those previously obtained from mid-infrared observations. Maps and source catalogues (>3-sigma) are now publicly released. Combined with the large wealth of multi-wavelength data available for the GOODS fields, these data provide a powerful new tool for studying galaxy evolution over a broad range of redshifts.

The Great Observatories All-Sky LIRG Survey (GOALS) is a multiwavelength study of luminous infrared galaxies (LIRGs) in the local universe. Here we present low resolution Spitzer spectra covering 5-38um and provide a basic analysis of the mid-IR spectral properties for nearby LIRGs. In a companion paper, we discuss detailed fits to the spectra. The GOALS sample of 244 nuclei in 180 luminous and 22 ultraluminous IR galaxies represents a complete subset of the IRAS RBGS and covers a range of merger stages, morphologies and spectral types. The majority (>60%) of GOALS LIRGs have high 6.2um PAH equivalent widths (EQW > 0.4um) and low levels of silicate absorption (s_9.7um >-1.0). There is a general trend among the U/LIRGs for silicate depth and MIR slope to increase with LIR. U/LIRGs in the late stages of a merger also have on average steeper MIR slopes and higher levels of dust obscuration. Together these trends suggest that as gas & dust is funneled towards the center of a coalescing merger, the nuclei become more compact and obscured. The sources that depart from these correlations have very low PAH EQW (EQW < 0.1um) consistent with their MIR emission being dominated by an AGN. The most heavily dust obscured sources are the most compact in their MIR emission, suggesting that the obscuring (cool) dust is associated with the outer regions of the starburst. As the merger progresses a marked decline is seen for the fraction of high EQW (star formation dominated) sources while the fraction of composite sources increases but the fraction of AGN-dominated sources remains low. When compared to the MIR spectra of submillimeter galaxies (SMGs) at z~2, the average GOALS LIRG is more absorbed at 9.7um and has more PAH emission. However, when the AGN contributions to both the local LIRGs and the high-z SMGs are removed, the average local starbursting LIRG closely resembles the starbursting SMGs.

(Abridged) We study relationships between the SFR and the nuclear properties of X-ray selected AGNs out to z=2.5, using far-IR data in three extragalactic deep fields as part of the PACS Evolutionary Probe (PEP) program. Guided by studies of intrinsic infra-red AGN SEDs, we show that the majority of the FIR emission in AGNs is produced by cold dust heated by star-formation. We uncover characteristic redshift-dependent trends between the mean FIR luminosity (L_fir) and accretion luminosity (L_agn) of AGNs. At low AGN luminosities, accretion and SFR are uncorrelated at all redshifts, consistent with a scenario where most low-luminosity AGNs are primarily fueled by secular processes in their host galaxies. At high AGN luminosities, a significant correlation is observed between L_fir and L_agn, but only among AGNs at low and moderate redshifts (z<1). We interpret this as a signature of the increasing importance of major-mergers in driving both the growth of super-massive black holes (SMBHs) and global star-formation in their hosts at high AGN luminosities. However, we also find that the enhancement of SFR in luminous AGNs weakens or disappears at high redshifts (z>1). This suggests that the role of mergers in SMBH-galaxy co-evolution is less important at these epochs. At all redshifts, we find essentially no relationship between L_fir and nuclear obscuration across five orders of magnitude in obscuring column density, suggesting that various different mechanisms are likely to be responsible for obscuring X-rays in active galaxies. We explain our results within a scenario in which two different modes of SMBH fueling operate among low- and high-luminosity AGNs. We postulate, guided by emerging knowledge about the properties of high redshift galaxies, that the dominant mode of accretion among high-luminosity AGNs evolves with redshift.

Aims. The aim of this work is to further investigate the nature of PAH excitation and emission especially in the context of tracing star formation in a variety of extragalactic environments. Here we turn our attention to the energetic environment of the closest AGN in our sample, Centaurus A. Methods. Using ISAAC on the ESO VLT UT1 (Antu) we have made high spatial resolution 3.3 {\mu}m imaging observations of the central kiloparsec of CenA. These observations have been compared with star formation tracers in the near- and mid-infrared, as well as with mid-infrared tracers of nuclear activity. Results. The nucleus is not devoid of PAH emission, implying that the PAH particles are not destroyed in the nucleus as might be expected for such a harsh environment. However, we see the feature to continuum ratio decrease towards the AGN. As well, the 3.3 {\mu}m PAH feature emission generally traces the sites of star formation in Cen A, but in detail there are spatial offsets, consistent with an earlier study of the starburst galaxies NGC 253 and NGC 1808. However, the feature-to-continuum ratio does not drop at the positions of star formation as was previously seen in that earlier study. The cause for this difference remains uncertain. Finally, our data reveal possible evidence for a nearly face-on, circular or spiral, dust structure surrounding the nucleus.

We present design results of the 2nd generation VLTI instrument GRAVITY beam stabilization and light injection subsystems. Designed to deliver micro-arcsecond astrometry, GRAVITY requires an unprecedented stability of the VLTI optical train. To meet the astrometric requirements, we have developed a dedicated 'laser guiding system', correcting the longitudinal and lateral pupil position as well as the image jitter. The actuators for the correction are provided by four 'fiber coupler' units located in the GRAVITY cryostat. Each fiber coupler picks the light of one telescope and stabilizes the beam. Furthermore each unit provides field de-rotation, polarization analysis as well as atmospheric piston correction. Using a novel roof prism design offers the possibility of on-axis as well as off-axis fringe tracking without changing the optical path. Finally the stabilized beam is injected with minimized losses into single-mode fibers via parabolic mirrors. We present lab results of the first guiding- as well as the first fiber coupler prototype regarding the closed loop performance and the optical quality. Based on the lab results we discuss the on-sky performance of the system and the implications concerning the sensitivity of GRAVITY.

Mrk 231, the nearest (z = 0.0422) quasar, hosts both a galactic-scale wind and a nuclear-scale iron low-ionization broad absorption line (FeLoBAL) outflow. We recently obtained a far-ultraviolet (FUV) spectrum of this object covering ~1150 - 1470 A with the Cosmic Origins Spectrograph on board the Hubble Space Telescope. This spectrum is highly peculiar, highlighted by the presence of faint (~< 2% of predictions based on H-alpha), broad (>~ 10,000 km/s at the base), and highly blueshifted (centroid at ~ -3500 km/s) Ly-alpha emission. The FUV continuum emission is slightly declining at shorter wavelengths (consistent with F_lambda ~ lambda^1.7) and does not show the presence of any obvious photospheric or wind stellar features. Surprisingly, the FUV spectrum also does not show any unambiguous broad absorption features. It thus appears to be dominated by the AGN, rather than hot stars, and virtually unfiltered by the dusty FeLoBAL screen. The observed Ly-alpha emission is best explained if it is produced in the outflowing BAL cloud system, while the Balmer lines arise primarily from the standard broad emission line region seen through the dusty (A_V ~ 7 mag.) broad absorption line region. Two possible geometric models are discussed in the context of these new results.

We use the infrared excess (IRX) FIR/UV luminosity ratio to study the relation between the effective UV attenuation (A_IRX) and the UV spectral slope (beta) in a sample of 450 1<z<2.5 galaxies. The FIR data is from very deep Herschel observations in the GOODS fields that allow us to detect galaxies with SFRs typical of galaxies with log(M)>9.3. Thus, we are able to study galaxies on and even below the main SFR-stellar mass relation (main sequence). We find that main sequence galaxies form a tight sequence in the IRX--beta plane, which has a flatter slope than commonly used relations. This slope favors a SMC-like UV extinction curve, though the interpretation is model dependent. The scatter in the IRX-beta plane, correlates with the position of the galaxies in the SFR-M plane. Using a smaller sample of galaxies with CO gas masses, we study the relation between the UV attenuation and the molecular gas content. We find a very tight relation between the scatter in the IRX-beta plane and the specific attenuation (S_A), a quantity that represents the attenuation contributed by the molecular gas mass per young star. S_A is sensitive to both the geometrical arrangement of stars and dust, and to the compactness of the star forming regions. We use this empirical relation to derive a method for estimating molecular gas masses using only widely available integrated rest-frame UV and FIR photometry. The method produces gas masses with an accuracy between 0.12-0.16 dex in samples of normal galaxies between z~0 and z~1.5. Major mergers and sub-millimeter galaxies follow a different S_A relation.

We report on Herschel/PACS observations of absorption lines of OH+, H2O+ and H3O+ in NGC 4418 and Arp 220. Excited lines of OH+ and H2O+ with E_lower of at least 285 and \sim200 K, respectively, are detected in both sources, indicating radiative pumping and location in the high radiation density environment of the nuclear regions. Abundance ratios OH+/H2O+ of 1-2.5 are estimated in the nuclei of both sources. The inferred OH+ column and abundance relative to H nuclei are (0.5-1)x10^{16} cm-2 and \sim2x10^{-8}, respectively. Additionally, in Arp 220, an extended low excitation component around the nuclear region is found to have OH+/H2O+\sim5-10. H3O+ is detected in both sources with N(H3O+)\sim(0.5-2)x10^{16} cm-2, and in Arp 220 the pure inversion, metastable lines indicate a high rotational temperature of ~500 K, indicative of formation pumping and/or hot gas. Simple chemical models favor an ionization sequence dominated by H+ - O+ - OH+ - H2O+ - H3O+, and we also argue that the H+ production is most likely dominated by X-ray/cosmic ray ionization. The full set of observations and models leads us to propose that the molecular ions arise in a relatively low density (\gtrsim10^4 cm-3) interclump medium, in which case the ionization rate per H nucleus (including secondary ionizations) is zeta>10^{-13} s-1, a lower limit that is severalx10^2 times the highest rate estimates for Galactic regions. In Arp 220, our lower limit for zeta is compatible with estimates for the cosmic ray energy density inferred previously from the supernova rate and synchrotron radio emission, and also with the expected ionization rate produced by X-rays. In NGC 4418, we argue that X-ray ionization due to an AGN is responsible for the molecular ion production.

We investigate the effect of the "high-pass filter" data reduction technique on the Herschel PACS PSF and noise of the PACS maps at the 70, 100 and 160 um bands and in medium and fast scan speeds. This branch of the PACS Photometer pipeline is the most used for cosmological observations and for point-source observations.The calibration of the flux loss due to the median removal applied by the PACS pipeline (high-pass filter) is done via dedicated simulations obtained by "polluting" real PACS timelines with fake sources at different flux levels. The effect of the data reduction parameter settings on the final map noise is done by using selected observations of blank fields with high data redundancy. We show that the running median removal can cause significant flux losses at any flux level. We analyse the advantages and disadvantages of several masking strategies and suggest that a mask based on putting circular patches on prior positions is the best solution to reduce the amount of flux loss. We provide a calibration of the point-source flux loss for several masking strategies in a large range of data reduction parameters, and as a function of the source flux. We also show that, for stacking analysis, the impact of the high-pass filtering effect is to reduce significantly the clustering effect. The analysis of the global noise and noise components of the PACS maps shows that the dominant parameter in determining the final noise is the high-pass filter width. We also provide simple fitting functions to build the error map from the coverage map and to estimate the cross-correlation correction factor in a representative portion of the data reduction parameter space.

(Abridged) The Magellanic Clouds provide a nearby laboratory for metal-poor dwarf galaxies. The low dust abundance enhances the penetration of UV photons into the interstellar medium (ISM), resulting in a relatively larger filling factor of the ionized gas. Furthermore, there is likely a hidden molecular gas reservoir probed by the [CII]157um line. We present Herschel/PACS maps in several tracers, [CII], [OI]63um,145um, [NII]122um, [NIII]57um, and [OIII]88um in the HII region N11B in the Large Magellanic Cloud. Halpha and [OIII]5007A images were used as complementary data to investigate the effect of dust extinction. Observations were interpreted with photoionization models to infer the gas conditions and estimate the ionized gas contribution to the [CII] emission. Photodissociation regions (PDRs) are probed through polycyclic aromatic hydrocarbons (PAHs). We first study the distribution and properties of the ionized gas. We then constrain the origin of [CII]157um by comparing to tracers of the low-excitation ionized gas and of PDRs. [OIII] is dominated by extended emission from the high-excitation diffuse ionized gas; it is the brightest far-infrared line, ~4 times brighter than [CII]. The extent of the [OIII] emission suggests that the medium is rather fragmented, allowing far-UV photons to permeate into the ISM to scales of >30pc. Furthermore, by comparing [CII] with [NII], we find that 95% of [CII] arises in PDRs, except toward the stellar cluster for which as much as 15% could arise in the ionized gas. We find a remarkable correlation between [CII]+[OI] and PAH emission, with [CII] dominating the cooling in diffuse PDRs and [OI] dominating in the densest PDRs. The combination of [CII] and [OI] provides a proxy for the total gas cooling in PDRs. Our results suggest that PAH emission describes better the PDR gas heating as compared to the total infrared emission.

We present deep CO observations of NGC6240 performed with the IRAM Plateau de Bure Interferometer (PdBI). NGC6240 is the prototypical example of a major galaxy merger in progress, caught at an early stage, with an extended, strongly-disturbed butterfly-like morphology and the presence of a heavily obscured active nucleus in the core of each progenitor galaxy. The CO line shows a skewed profile with very broad and asymmetric wings detected out to velocities of -600 km/s and +800 km/s with respect to the systemic velocity. The PdBI maps reveal the existence of two prominent structures of blueshifted CO emission. One extends eastward, i.e. approximately perpendicular to the line connecting the galactic nuclei, over scales of ~7 kpc and shows velocities up to -400 km/s. The other extends southwestward out to ~7 kpc from the nuclear region, and has a velocity of -100 km/s with respect to the systemic one. Interestingly, redshifted emission with velocities 400 to 800 km/s is detected around the two nuclei, extending in the east-west direction, and partly overlapping with the eastern blue-shifted structure, although tracing a more compact region of size ~1.7 kpc. The overlap between the southwestern CO blob and the dust lanes seen in HST images, which are interpreted as tidal tails, indicates that the molecular gas is deeply affected by galaxy interactions. The eastern blueshifted CO emission is co-spatial with an Halpha filament that is associated with strong H2 and soft X-ray emission. The analysis of Chandra X-ray data provides strong evidence for shocked gas at the position of the Halpha emission. Its association with outflowing molecular gas supports a scenario where the molecular gas is compressed into a shock wave that propagates eastward from the nuclei. If this is an outflow, the AGN are likely the driving force.

[Abridged] We present maps of the main cooling lines of the neutral atomic gas ([OI] at 63 and 145 micron and [CII] at 158 micron) and in the [OIII] 88 micron line of the starburst galaxy M82, carried out with the PACS spectrometer on board the Herschel satellite. By applying PDR modeling we derive maps of the main ISM physical parameters, including the [CII] optical depth, at unprecedented spatial resolution (~300 pc). We can clearly kinematically separate the disk from the outflow in all lines. The [CII] and [OI] distributions are consistent with PDR emission both in the disk and in the outflow. Surprisingly, in the outflow, the atomic and the ionized gas traced by the [OIII] line both have a deprojected velocity of ~75 km/s, very similar to the average velocity of the outflowing cold molecular gas (~ 100 km/s) and several times smaller than the outflowing material detected in Halpha (~ 600 km/s). This suggests that the cold molecular and neutral atomic gas and the ionized gas traced by the [OIII] 88 micron line are dynamically coupled to each other but decoupled from the Halpha emitting gas. We propose a scenario where cold clouds from the disk are entrained into the outflow by the winds where they likely evaporate, surviving as small, fairly dense cloudlets (n_H\sim 500-1000 cm^-3, G_0\sim 500- 1000, T_gas\sim300 K). We show that the UV photons provided by the starburst are sufficient to excite the PDR shells around the molecular cores. The mass of the neutral atomic gas in the outflow is \gtrsim 5-12x 10^7 M_sun to be compared with that of the molecular gas (3.3 x 10^8 M_sun) and of the Halpha emitting gas (5.8 x 10^6 M_sun). The mass loading factor, (dM/dt)/SFR, of the molecular plus neutral atomic gas in the outflow is ~ 2. Energy and momentum driven outflow models can explain the data equally well, if all the outflowing gas components are taken into account.

GRAVITY is a second generation VLTI instrument, combining the light of four telescopes and two objects simultaneously. The main goal is to obtain astrometrically accurate information. Besides correctly measured stellar phases this requires the knowledge of the instrumental differential phase, which has to be measured optically during the astronomical observations. This is the purpose of a dedicated metrology system. The GRAVITY metrology covers the full optical path, from the beam combiners up to the reference points in the beam of the primary telescope mirror, minimizing the systematic uncertainties and providing a proper baseline in astrometric terms. Two laser beams with a fixed phase relation travel backward the whole optical chain, creating a fringe pattern in any plane close to a pupil. By temporal encoding the phase information can be extracted at any point by means of flux measurements with photo diodes. The reference points chosen sample the pupil at typical radii, eliminating potential systematics due differential focus. We present the final design and the performance estimate, which is in accordance with the overall requirements for GRAVITY.

Results of observations with the Spitzer, Hubble, GALEX, Chandra, and XMM-Newton space telescopes are presented for the Luminous Infrared Galaxy (LIRG) merger Mrk 266. The SW (Seyfert 2) and NE (LINER) nuclei reside in galaxies with Hubble types SBb (pec) and S0/a (pec), respectively. Both galaxies have L > L*, and they are inferred to each contain a ~2.5x10^8 M_sun black hole. Mrk 266 SW is likely the primary source of a bright Fe K-alpha line detected from the system, consistent with the reflection-dominated X-ray spectrum of a heavily obscured AGN. Optical knots embedded in an arc with aligned radio continuum radiation, combined with luminous H_2 line emission, provide evidence for a radiative bow shock in an AGN-driven outflow surrounding the NE nucleus. Soft X-ray emission modeled as shock-heated plasma is co-spatial with radio continuum emission between the galaxies. Mid-IR diagnostics indicate roughly equal contributions of AGN and starburst radiation powering the bolometric luminosity. Approximately 120 star clusters have been detected, with most having estimated ages < 50 Myr. Detection of 24 micron emission extending ~34 arcsec (20 kpc) north of the galaxies is interpreted as ~2x10^7 M_sun of dust entrained in an outflowing superwind; at optical wavelengths this region is resolved into a fragmented morphology indicative of Rayleigh-Taylor instabilities in an expanding shell of ionized gas. Mrk 266 demonstrates that the dust "blow-out" phase can begin in a LIRG well before the galaxies fully coalesce during a subsequent ULIRG phase, and rapid gas consumption in luminous dual AGNs (kpc scale separations) early in the merger process may explain the paucity of detected binary QSOs (sub-pc scale orbits) in large surveys. An evolutionary sequence is proposed representing a progression from dual to binary AGNs, accompanied by an increase in observed L_x/L_ir ratios by 10^4 or more.

Massive AGN-driven outflows are invoked by AGN-galaxy co-evolutionary models to suppress both star formation and black hole accretion. Massive molecular outflows have recently been revealed in some AGN hosts. However, the physical properties and structure of these AGN-driven molecular outflows are still poorly constrained. Here we present new IRAM PdBI observations of Mrk231, the closest quasar known, targeting both the CO(1-0) and CO(2-1) transitions. We detect broad wings in both transitions, tracing a massive molecular outflow with velocities up to 800 km/s. The wings are spatially resolved at high significance level (5-11 sigma), indicating that the molecular outflow extends on the kpc scale. The CO(2-1)/CO(1-0) ratio of the red broad wings is consistent with the ratio observed in the narrow core, while the blue broad wing is less excited than the core. The latter result suggests that quasar driven outflow models invoking shocks (which would predict higher gas excitation) are not appropriate to describe the bulk of the outflow in Mrk231. However, we note that within the central 700 pc the CO(2-1)/CO(1-0) ratio of the red wing is slightly, but significantly, higher than in the line core, suggesting that shocks may play a role in the central region. We also find that the average size of the outflow anticorrelates with the critical density of the transition used as a wind tracer. This indicates that, although diffuse and dense clumps coexist in the outflowing gas, dense outflowing clouds have shorter lifetime and that they evaporate into the diffuse component along the outflow or, more simply, that diffuse clouds are more efficiently accelerated to larger distances by radiation pressure.

Most theoretical models invoke quasar driven outflows to quench star formation in massive galaxies, this feedback mechanism is required to account for the population of old and passive galaxies observed in the local universe. The discovery of massive, old and passive galaxies at z=2, implies that such quasar feedback onto the host galaxy must have been at work very early on, close to the reionization epoch. We have observed the [CII]158um transition in SDSSJ114816.64+525150.3 that, at z=6.4189, is one of the most distant quasars known. We detect broad wings of the line tracing a quasar-driven massive outflow. This is the most distant massive outflow ever detected and is likely tracing the long sought quasar feedback, already at work in the early Universe. The outflow is marginally resolved on scales of about 16 kpc, implying that the outflow can really affect the whole galaxy, as required by quasar feedback models. The inferred outflow rate, dM/dt > 3500 Msun/yr, is the highest ever found. At this rate the outflow can clean the gas in the host galaxy, and therefore quench star formation, in a few million years.

We report the detection of far-IR CO rotational emission from the prototypical Seyfert 2 galaxy NGC 1068. Using Herschel-PACS, we have detected 11 transitions in the J_upper=14-30 (E_upper/k_B = 580-2565 K) range, all of which are consistent with arising from within the central 10" (700 pc). The detected transitions are modeled as arising from 2 different components: a moderate excitation (ME) component close to the galaxy systemic velocity, and a high excitation (HE) component that is blueshifted by ~80 km s^{-1}. We employ a large velocity gradient (LVG) model and derive n_H2~10^{5.6} cm^{-3}, T_kin~170 K, and M_H2~10^{6.7} M_sun for the ME component, and n_H2~10^{6.4} cm^{-3}, T_kin~570 K, and M_H2~10^{5.6} M_sun for the HE component, although for both components the uncertainties in the density and mass are plus/minus (0.6-0.9) dex. We compare the CO line profiles with those of other molecular tracers observed at higher spatial and spectral resolution, and find that the ME transitions are consistent with these lines arising in the ~200 pc diameter ring of material traced by H_2 1-0 S(1) observations. The blueshift of the HE lines may also be consistent with the bluest regions of this H_2 ring, but a better kinematic match is found with a clump of infalling gas ~40 pc north of the AGN. We discuss the prospects of placing the HE component near the AGN, and conclude that while the moderate thermal pressure precludes an association with the ~1 pc radius H_2O maser disk, the HE component could potentially be located only a few parsecs more distant from the AGN, and might then provide the N_H~10^{25} cm^{-2} column obscuring the nuclear hard X-rays. Finally, we also report sensitive upper limits extending up to J_upper=50, which place constraints on a previous model prediction for the CO emission from the X-ray obscuring torus. [Abridged]

The analysis of the physical properties of low-redshift Ly$\alpha$ emitters (LAEs) can provide clues in the study of their high-redshift analogues. At $z \sim 0.3$, LAEs are bright enough to be detected over almost the entire electromagnetic spectrum and it is possible to carry out a more precise and complete study than at higher redshifts. In this study, we examine the UV and IR emission, dust attenuation, SFR and morphology of a sample of 23 GALEX-discovered star-forming (SF) LAEs at $z \sim 0.3$ with direct UV (GALEX), optical (ACS) and FIR (PACS and MIPS) data. Using the same UV and IR limiting luminosities, we find that LAEs at $z\sim 0.3$ tend to be less dusty, have slightly higher total SFRs, have bluer UV continuum slopes, and are much smaller than other galaxies that do not exhibit Ly$\alpha$ emission in their spectrum (non-LAEs). These results suggest that at $z \sim 0.3$ Ly$\alpha$ photons tend to escape from small galaxies with low dust attenuation. Regarding their morphology, LAEs belong to Irr/merger classes, unlike non-LAEs. Size and morphology represent the most noticeable difference between LAEs and non-LAEs at $z \sim 0.3$. Furthermore, the comparison of our results with those obtained at higher redshifts indicates that either the Ly$\alpha$ technique picks up different kind of galaxies at different redshifts or that the physical properties of LAEs are evolving with redshift.

We compare the average star formation (SF) activity in X-ray selected AGN hosts with mass-matched control inactive galaxies,including star forming and quiescent sources, at 0.5<z<2.5. Recent observations carried out by PACS, the 60-210um Herschel photometric camera, in GOODS-S, GOODS-N and COSMOS allow us to unbiasedly estimate the far-IR luminosity, and hence the SF properties, of the two samples. Accurate AGN host stellar masses are measured by decomposing their total emission into the stellar and nuclear components. We find a higher average SF activity in AGN hosts with respect to non-AGNs. The level of SF enhancement is modest (~0.26dex at ~3sigma) at low X-ray luminosities (Lx<~10^43.5erg/s) and more pronounced (0.56dex at >10sigma) for bright AGNs. However, when comparing to star forming galaxies only, AGN hosts are broadly consistent with the locus of their `main sequence'. We investigate the relative far-IR luminosity distributions of active and inactive galaxies, and find a higher fraction of PACS detected, hence normal and highly star forming systems among AGN hosts. Although different interpretations are possible, we explain our findings as a consequence of a twofold AGN growth path: faint AGNs evolve through secular processes, with instantaneous AGN accretion not tightly linked to the current total SF in the host, while luminous AGNs co-evolve with their hosts through periods of enhanced AGN activity and SF, possibly through major mergers. While an increased SF with respect to non-AGNs of similar mass is expected in the latter, we interpret the modest SF offsets measured in low-Lx AGN hosts as either a) generated by non-synchronous accretion and SF histories in a merger scenario or b) due to possible connections between instantaneous SF and accretion that can be induced by smaller scale (non-major merger) mechanisms. Far-IR luminosity distributions favour the latter scenario.

Herschel/PACS spectroscopy of the luminous infrared galaxies NGC4418 and Arp220 reveals high excitation in H2O, OH, HCN, and NH3. In NGC4418, absorption lines were detected with E_low>800 K (H2O), 600 K (OH), 1075 K (HCN), and 600 K (NH3), while in Arp220 the excitation is somewhat lower. While outflow signatures in moderate excitation lines are seen in Arp220 as reported in previous studies, in NGC4418 the lines tracing its outer regions are redshifted relative to the nucleus, suggesting an inflow with Mdot<~12 Msun yr^{-1}. Both galaxies have warm (Tdust>~100 K) nuclear continuum components, together with a more extended component that is much more prominent and massive in Arp220. A chemical dichotomy is found in both sources: on the one hand, the nuclear regions have high H2O abundances, ~10^{-5}, and high HCN/H2O and HCN/NH3 column density ratios of 0.1-0.4 and 2-5, respectively, indicating a chemistry typical of evolved hot cores where grain mantle evaporation has occurred. On the other hand, the high OH abundance, with OH/H2O ratios of ~0.5, indicates the effects of X-rays and/or cosmic rays. The nuclear media have surface brightnesses >~10^{13} Lsun/kpc^2 and are estimated to be thick (N_H>~10^{25} cm^{-2}). While NGC4418 shows weak absorption in H2^{18}O and ^{18}OH, with a ^{16}O-to-^{18}O ratio of >~250-500, the strong absorption of the rare isotopologues in Arp220 indicates ^{16}O-to-^{18}O of 70-130. Further away from the nuclear regions, the H2O abundance decreases to <~10^{-7} and the OH/H2O ratio is 2.5-10. Despite the different scales of NGC4418, Arp220, and Mrk231, preliminary evidence is found for an evolutionary sequence from infall, hot-core like chemistry, and solar oxygen isotope ratio to high velocity outflow, disruption of the hot core chemistry and cumulative high mass stellar processing of 18O.

In this work we analyze the physical properties of a sample of 56 spectroscopically selected star-forming (SF) Ly$\alpha$ emitting galaxies at 2.0$\lesssim$z$\lesssim$3.5 using both a spectral energy distribution (SED) fitting procedure from rest-frame UV to mid-IR and direct 160$\mu$m observations taken with the Photodetector Array Camera & Spectrometer (PACS) instrument onboard \emph{Herschel Space Observatory}. We define LAEs as those Ly$\alpha$ emitting galaxies whose rest-frame Ly$\alpha$ equivalent widths (Ly$\alpha$ EW$_{rest-frame}$) are above 20\AA, the typical threshold in narrow-band searches. Ly$\alpha$ emitting galaxies with Ly$\alpha$ EW$_{rest-frame}$ are called non-LAEs. As a result of an individual SED fitting for each object, we find that the studied sample of LAEs contains galaxies with ages mostly below 100Myr and a wide variety of dust attenuations, SFRs, and stellar masses. The heterogeneity in the physical properties is also seen in the morphology, ranging from bulge-like galaxies to highly clumpy systems. In this way, we find that LAEs at 2.0$\lesssim$z$\lesssim$3.5 are very diverse, and do not have a bimodal nature, as suggested in previous works. Furthermore, the main difference between LAEs and non-LAEs is their dust attenuation, because LAEs are not as dusty as non-LAEs. On the FIR side, four galaxies of the sample (two LAEs and two non-LAEs) have PACS-FIR counterparts. Their total IR luminosity place all of them in the ULIRG regime and are all dusty objects, with A$_{1200}$$\gtrsim$4mag. This is an indication from direct FIR measurements that dust and Ly$\alpha$ emission are not mutually exclusive. This population of red and dusty LAEs is not seen at z$\sim$0.3, suggesting an evolution with redshift of the IR nature of galaxies selected via their Ly$\alpha$ emission.

We study a sample of 61 submillimetre galaxies (SMGs) selected from ground-based surveys, with known spectroscopic redshifts and observed with Herschel as part of the PACS Evolutionary Probe (PEP) and the Herschel Multi-tiered Extragalactic Survey (HerMES) key programmes. We use the broad far-infrared wavelength coverage (100-600um) provided by the combination of PACS and SPIRE observations. Using a power-law temperature distribution model to derive infrared luminosities and dust temperatures, we measure a dust emissivity spectral index for SMGs of beta=2.0+/-0.2. Our results unveil the diversity of the SMG population. Some SMGs exhibit extreme infrared luminosities of ~10^13 Lsun and relatively warm dust components, while others are fainter (~10^12 Lsun) and are biased towards cold dust temperatures. The extreme infrared luminosities of some SMGs (LIR>10^12.7 Lsun, 26/61 systems) imply SFRs of >500Msun yr^-1. Such high SFRs are difficult to reconcile with a secular mode of star formation, and may instead correspond to a merger-driven stage in the evolution of these galaxies. Another observational argument in favour of this scenario is the presence of dust temperatures warmer than that of SMGs of lower luminosities (~40K as opposed to ~25K), consistent with observations of local ULIRGs triggered by major mergers and with results from hydrodynamic simulations of major mergers combined with radiative transfer calculations. Luminous SMGs are also offset from normal star-forming galaxies in the stellar mass-SFR plane, suggesting that they are undergoing starburst events with short duty cycles, compatible with the major merger scenario. On the other hand, a significant fraction of the low infrared luminosity SMGs have cold dust temperatures, are located close to the main sequence of star formation, and thus might be evolving through a secular mode of star formation. [abridged]

We combine Herschel-PACS data from the PEP program with Spitzer 24 um and 16 um photometry and ultra deep IRS mid-infrared spectra, to measure the mid- to far-infrared spectral energy distribution (SED) of 0.7<z<2.5 normal star forming galaxies around the main sequence (the redshift-dependent relation of star formation rate and stellar mass). Our deep data confirm from individual far-infrared detections that z~2 star formation rates are overestimated if based on 24 um fluxes and SED templates that are calibrated via local trends with luminosity. Galaxies with similar ratios of rest-frame nuLnu(8) to 8-1000 um infrared luminosity (LIR) tend to lie along lines of constant offset from the main sequence. We explore the relation between SED shape and offset in specific star formation rate (SSFR) from the redshift-dependent main sequence. Main sequence galaxies tend to have a similar nuLnu(8)/LIR regardless of LIR and redshift, up to z~2.5, and nuLnu(8)/LIR decreases with increasing offset above the main sequence in a consistent way at the studied redshifts. We provide a redshift-independent calibration of SED templates in the range of 8--60 um as a function of log(SSFR) offset from the main sequence. Redshift dependency enters only through the evolution of the main sequence with time. Ultra deep IRS spectra match these SED trends well and verify that they are mostly due to a change in ratio of PAH to LIR rather than continua of hidden AGN. Alternatively, we discuss the dependence of nuLnu(8)/LIR on LIR. Same nuLnu(8)/LIR is reached at increasingly higher LIR at higher redshift, with shifts relative to local by 0.5 and 0.8 dex in log(LIR) at redshifts z~1 and z~2. Corresponding SED template calibrations are provided for use if no stellar masses are in hand. For most of those z~2 star forming galaxies that also host an AGN, the mid-infrared is dominated by the star forming component.

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.

Star formation is arguably the most important physical process in the cosmos. It is a fundamental driver of galaxy evolution and the ultimate source of most of the energy emitted by galaxies. A correct interpretation of star formation rate (SFR) measures is therefore essential to our understanding of galaxy formation and evolution. Unfortunately, however, no single SFR estimator is universally available or even applicable in all circumstances: the numerous galaxies found in deep surveys are often too faint (or too distant) to yield significant detections with most standard SFR measures, and until now there have been no global, multi-band observations of nearby galaxies that span all the conditions under which star-formation is taking place. To address this need in a systematic way, we have undertaken a multi-band survey of all types of star-forming galaxies in the local Universe. This project, the Star Formation Reference Survey (SFRS), is based on a statistically valid sample of 369 nearby galaxies that span all existing combinations of dust temperature, SFR, and specific SFR. Furthermore, because the SFRS is blind with respect to AGN fraction and environment it serves as a means to assess the influence of these factors on SFR. Our panchromatic global flux measurements (including GALEX FUV+NUV, SDSS ugriz, 2MASS JHKs, Spitzer 3-8{\mu}m, and others) furnish uniform SFR measures and the context in which their reliability can be assessed. This paper describes the SFRS survey strategy, defines the sample, and presents the multi-band photometry collected to date.

One remaining open question regarding the physical properties of Ly$\alpha$ emitters (LAEs) is their dust content and its evolution with redshift. The variety of results is large and with those reported by now is difficult to establish clear relations between dust, other fundamental parameters of galaxies (star-formation rate, metallicity or age) and redshift. In this Letter, we report \emph{Herschel} PACS-100$\mu$m, PACS-160$\mu$m and \emph{Spitzer} MIPS-24$\mu$m detections of a sample of spectroscopically GALEX-selected LAEs at z$\sim$0.3 and $\sim$1.0. Five out of ten and one out of two LAEs are detected in, at least, one PACS band at z$\sim$0.3 and $\sim$1.0, respectively. These measurements have a great importance given that they allow us to quantify, for the first time, the dust content of LAEs from direct FIR observations. MIPS-24$\mu$m detections allow us to determine IR properties of the PACS-undetected LAEs. We obtain that mid-IR/FIR detected star-forming (SF) LAEs at z$\sim$0.3 have dust content within 0.75$\lesssim$ $A_{1200\AA}$ $\lesssim$2.0, with a median value of A$_{1200\textrm{\AA}}$$\sim$1.1. This range broadens out to 0.75$\lesssim$ $A_{1200\AA}$ $\lesssim$2.5 when considering those LAEs at z$\sim$1.0. Only one SF LAE is undetected both in MIPS-24$\mu$m and PACS, with $A_{1200\AA}$ $\lesssim$0.75. These results seem to be larger than those reported for high-redshift LAEs and, therefore, although an evolutionary trend is not clearly seen, it could point out that low-redshift LAEs are dustier than high-redshift ones. However, the diverse methods used could introduce a systematic offset in the results.

We compare multi-wavelength SFR indicators out to z~3 in GOODS-South. Our analysis uniquely combines U-to-8um photometry from FIREWORKS, MIPS 24um and PACS 70, 100, and 160um photometry from the PEP survey, and Ha spectroscopy from the SINS survey. We describe a set of conversions that lead to a continuity across SFR indicators. A luminosity-independent conversion from 24um to total infrared luminosity yields estimates of LIR that are in the median consistent with the LIR derived from PACS photometry, albeit with significant scatter. Dust correction methods perform well at low to intermediate levels of star formation. They fail to recover the total amount of star formation in systems with large SFR_IR/SFR_UV ratios, typically occuring at the highest SFRs (SFR_UV+IR \gtrsim 100 Msun/yr) and redshifts (z \gtrsim 2.5) probed. Finally, we confirm that Ha-based SFRs at 1.5<z<2.6 are consistent with SFR_SED and SFR_UV+IR provided extra attenuation towards HII regions is taken into account (Av,neb = Av,continuum / 0.44). With the cross-calibrated SFR indicators in hand, we perform a consistency check on the star formation histories inferred from SED modeling. We compare the observed SFR-M relations and mass functions at a range of redshifts to equivalents that are computed by evolving lower redshift galaxies backwards in time. We find evidence for underestimated stellar ages when no stringent constraints on formation epoch are applied. We demonstrate how resolved SED modeling, or alternatively deep UV data, may help to overcome this bias. The age bias is most severe for galaxies with young stellar populations, and reduces towards older systems. Finally, our analysis suggests that SFHs typically vary on timescales that are long (at least several 100 Myr) compared to the galaxies' dynamical time.

Deep far-infrared photometric surveys studying galaxy evolution and the nature of the cosmic infrared background are a key strength of the Herschel mission. We describe the scientific motivation for the PACS Evolutionary Probe (PEP) guaranteed time key program and its role in the complement of Herschel surveys, and the field selection which includes popular multiwavelength fields such as GOODS, COSMOS, Lockman Hole, ECDFS, EGS. We provide an account of the observing strategies and data reduction methods used. An overview of first science results illustrates the potential of PEP in providing calorimetric star formation rates for high redshift galaxy populations, thus testing and superseeding previous extrapolations from other wavelengths, and enabling a wide range of galaxy evolution studies.

The cosmic infrared background (CIB) includes roughly half of the energy radiated by all galaxies at all wavelengths across cosmic time, as observed at the present epoch. The PACS Evolutionary Probe (PEP) survey is exploited here to study the CIB and its redshift differential, at 70, 100 and 160 micron, where the background peaks. Combining PACS observations of the GOODS-S, GOODS-N, Lockman Hole and COSMOS areas, we define number counts spanning over more than two orders of magnitude in flux: from ~1 mJy to few hundreds mJy. Stacking of 24 micron sources and P(D) statistics extend the analysis down to ~0.2 mJy. Taking advantage of the wealth of ancillary data in PEP fields, differential number counts and CIB are studied up to z=5. Based on these counts, we discuss the effects of confusion on PACS blank field observations and provide confusion limits for the three bands considered. The total CIB surface brightness emitted above PEP 3 sigma flux limits is 4.52 +/- 1.18, 8.35 +/- 0.95 and 9.49 +/- 0.59 [nW/m2/sr] at 70, 100, and 160 micron, respectively. These values correspond to 58 +/- 7% and 74 +/- 5% of the COBE/DIRBE CIB direct measurements at 100 and 160 micron. Employing the P(D) analysis, these fractions increase to ~65% and ~89%. More than half of the resolved CIB was emitted at redshift z<=1. The 50%-light redshifts lie at z=0.58, 0.67 and 0.73 at the three PACS wavelengths. The distribution moves towards earlier epochs at longer wavelengths: while the 70 micron CIB is mainly produced by z<=1.0 objects, the contribution of z>1.0 sources reaches 50% at 160 micron. Most of the CIB resolved in the three PACS bands was emitted by galaxies with infrared luminosities in the range 1e11-1e12 L(sun).

ABRIDGED-This paper presents the first direct estimate of the 3D clustering properties of far-infrared sources up to z~3. This has been possible thanks to the Pacs Evolutionary Probe (PEP) survey of the GOODS South field performed with the PACS instrument onboard the Herschel Satellite. An analysis of the two-point correlation function over the whole redshift range spanned by the data reports for the correlation length, r_0~6.3 Mpc and r_0~6.7 Mpc, respectively at 100um and 160um, corresponding to dark matter halo masses M>~10^{12.4} M_sun. Objects at z~2 instead seem to be more strongly clustered, with r_0~19 Mpc and r_0~17 Mpc in the two considered PACS channels. This dramatic increase of the correlation length between z~1 and z~2 is connected with the presence of a wide, M>~10^{14} M_sun, filamentary structure which includes more than 50% of the sources detected at z~2. An investigation of the properties of such sources indicates the possibility for boosted star-forming activity in those which reside within the overdense environment with respect of more isolated galaxies found in the same redshift range. Lastly, we also present our results on the evolution of the relationship between luminous and dark matter in star-forming galaxies between z~1 and z~2. We find that the increase of (average) stellar mass in galaxies <M*> between z~1 and z~2 is about a factor 10 lower than that of the dark matter haloes hosting such objects (<M*>[z~1]/<M*>[z~2] ~ 0.4 vs M_{halo}[z~1]/M_{halo}[z~2] ~ 0.04). Our findings agree with the evolutionary picture of downsizing whereby massive galaxies at z~2 were more actively forming stars than their z~1 counterparts, while at the same time contained a lower fraction of their mass in the form of luminous matter.

The Herschel Dwarf Galaxy Survey investigates the interplay of star formation activity and the the metal-poor gas and dust of dwarf galaxies using FIR and submillimetre imaging spectroscopic and photometric observations in the 50 to 550mu window of the Herschel Space Observatory. The dust SEDs are well constrained with the new Herschel and MIR Spitzer data. A submillimetre excess is often found in low metallicity galaxies, which,if tracing very cold dust, would highlight large dust masses not easily reconciled in some cases, given the low metallicities and expected gas-to-dust mass ratios. The galaxies are also mapped in the FIR fine-structure lines (63 and 145mu OI, 158mu CII, 122 and 205mu NII, 88mu OIII) probing the low density ionised gas, the HII regions and photodissociation regions. While still early in the Herschel mission we can already see, along with earlier studies, that line ratios in the metal-poor ISM differ remarkably from those in the metal-rich starburst environments. In dwarf galaxies, L[CII]/L(CO) (>10^4) is at least an order of magnitude greater than in the most metal-rich starburst galaxies. The enhanced [CII] arises from the larger photodissociation region where H2, not traced by the CO, can exist. The 88mu [OIII] line usually dominates the FIR line emission over galaxy-wide scales in dwarf galaxies, not the 158mu [CII] line which is the dominant FIR cooling line in metal-rich galaxies. All of the FIR lines together can contribute 1% to 2% of the L(TIR). The Herschel Dwarf Galaxy survey will provide statistical information on the nature of the dust and gas in low metallicity galaxies, elucidating the origin of the submm excess in dwarf galaxies, and help determine a ([CII] +CO) to H2 conversion factor, thus providing observational constraints on chemical evolution models of galaxies.

Mass outflows driven by stars and active galactic nuclei are a key element in many current models of galaxy evolution. They may produce the observed black hole-galaxy mass relation and regulate and quench both star formation in the host galaxy and black hole accretion. However, observational evidence of such feedback processes through outflows of the bulk of the star forming molecular gas is still scarce. Here we report the detection of massive molecular outflows, traced by the hydroxyl molecule (OH), in far-infrared spectra of ULIRGs obtained with Herschel-PACS as part of the SHINING key project. In some of these objects the (terminal) outflow velocities exceed 1000 km/s, and their outflow rates (up to ~1200 M_sol/yr) are several times larger than their star formation rates. We compare the outflow signatures in different types of ULIRGs and in starburst galaxies to address the issue of the energy source (AGN or starburst) of these outflows. We report preliminary evidence that ULIRGs with a higher AGN luminosity (and higher AGN contribution to L_IR) have higher terminal velocities and shorter gas depletion time scales. The outflows in the observed ULIRGs are able to expel the cold gas reservoirs from the centres of these objects within ~1E6-1E8 years.

We analyse new spatially resolved integral field spectroscopic H- and K-band data at a resolution of 0.3arcsec (100pc) and re-analyse interferometric CO(2-1) line observations of the prototypical merging system Arp220. We find that the majority of the K-band luminosity is due to a 10Myr old starburst, with a significant contribution from an underlying ca 1Gyr old stellar population, and a small contribution from stars less than 8Myr old. The Calzetti (2000) reddening law provides the best fit to photometric datapoints spanning 0.45um -- 2.12um. Furthermore, estimates of the bolometric luminosity from IRAS fluxes in conjunction with our stellar population analysis indicate that we observe less than 10% of the emitted K-band light. The stellar and CO(2-1) kinematic centre of the western nucleus coincides with the compact hot dust emission, indicating that the latter marks the centre of the gravitational potential. In the eastern nucleus, the CO(2-1) data are well matched by a model in which the gas orbits around the peak of the dust emission. This, and the similarity of the K-band tracer kinematics, shows that despite the irregular morphology, the eastern nucleus is also a kinematically coherent structure. Comparison of the extinction map and EW_CO and EW_Br-gamma maps indicates that the lower half of the eastern nucleus is significantly more extincted than the upper half, suggesting that the lower half is buried in the larger-scale gas disk.

We report initial results from the far-infrared fine structure line observations of a sample of 44 local starbursts, Seyfert galaxies and infrared luminous galaxies obtained with the PACS spectrometer on board Herschel. We show that the ratio between the far-infrared luminosity and the molecular gas mass, Lfir/MH2, is a much better proxy for the relative brightness of the far-infrared lines than Lfir alone. Galaxies with high Lfir/MH2 ratios tend to have weaker fine structure lines relative to their far-infrared continuum than galaxies with Lfir/MH2 < 80 Lsun/Msun. A deficit of the [C II] 158 micron line relative to Lfir was previously found with the ISO satellite, but now we show for the first time that this is a general aspect of all far-infrared fine structure lines, regardless of their origin in the ionized or neutral phase of the interstellar medium. The Lfir/MH2 value where these line deficits start to manifest is similar to the limit that separates between the two modes of star formation recently found in galaxies on the basis of studies of their gas-star formation relations. Our finding that the properties of the interstellar medium are also significantly different in these regimes provides independent support for the different star forming relations in normal disk galaxies and major merger systems. We use the spectral synthesis code Cloudy to model the emission of the lines. The expected increase of the ionization parameter with Lfir/MH2 can simultaneously explain the line deficits in the [C II], [N II] and [O I] lines.

We present a statistical analysis of the mid-infrared (MIR) spectra of 248 luminous infrared (IR) galaxies (LIRGs) which comprise the Great Observatories All-sky LIRG Survey (GOALS) observed with the Infrared Spectrograph (IRS) on-board the Spitzer Space Telescope. The GOALS sample enables a direct measurement of the relative contributions of star-formation and active galactic nuclei (AGN) to the total IR emission from a large sample of local LIRGs. The AGN contribution to the MIR emission (f-AGN) is estimated by employing several diagnostics based on the properties of the [NeV], [OIV] and [NeII] fine structure gas emission lines, the 6.2 microns PAH and the shape of the MIR continuum. We find that 18% of all LIRGs contain an AGN and that in 10% of all sources the AGN contributes more than 50% of the total IR luminosity. Summing up the total IR luminosity contributed by AGN in all our sources suggests that AGN supply ~12% of the total energy emitted by LIRGs. The average spectrum of sources with an AGN looks similar to the average spectrum of sources without an AGN, but it has lower PAH emission and a flatter MIR continuum. AGN dominated LIRGs have higher IR luminosities, warmer MIR colors and are found in interacting systems more often than pure starbursts LIRGs. However we find no linear correlations between these properties and f-AGN. We used the IRAC colors of LIRGs to confirm that finding AGN on the basis of their MIR colors may miss ~40% of AGN dominated (U)LIRGs

We present spatially resolved integral field spectroscopic K-band data at a resolution of 0.13" (60pc) and interferometric CO(2-1) line observations of the prototypical merging system NGC6240. Despite the clear rotational signature, the stellar kinematics in the two nuclei are dominated by dispersion. We use Jeans modelling to derive the masses and the mass-to-light ratios of the nuclei. Combining the luminosities with the spatially resolved Br-gamma equivalent width shows that only 1/3 of the K-band continuum from the nuclei is associated with the most recent star forming episode; and that less than 30% of the system's bolometric luminosity and only 9% of its stellar mass is due to this starburst. The star formation properties, calculated from typical merger star formation histories, demonstrate the impact of different assumptions about the star formation history. The properties of the nuclei, and the existence of a prominent old stellar population, indicate that the nuclei are remnants of the progenitor galaxies' bulges.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

Ly-alpha emitters (LAEs) are seen everywhere in the redshift domain from local to z~7. Far-infrared (FIR) counterparts of LAEs at different epochs could provide direct clues on dust content, extinction, and spectral energy distribution (SED) for these galaxies. We search for FIR counterparts of LAEs that are optically detected in the GOODS-North field at redshift z~2.2 using data from the Herschel Space Telescope with the Photodetector Array Camera and Spectrometer (PACS). The LAE candidates were isolated via color-magnitude diagram using the medium-band photometry from the ALHAMBRA Survey, ancillary data on GOODS-North, and stellar population models. According to the fitting of these spectral synthesis models and FIR/optical diagnostics, most of them seem to be obscured galaxies whose spectra are AGN-dominated. From the analysis of the optical data, we have observed a fraction of AGN or composite over source total number of ~0.75 in the LAE population at z~2.2, which is marginally consistent with the fraction previously observed at z=2.25 and even at low redshift (0.2<z<0.45), but significantly different from the one observed at redshift ~3, which could be compatible either with a scenario of rapid change in the AGN fraction between the epochs involved or with a non detection of obscured AGN in other z=2-3 LAE samples due to lack of deep FIR observations. We found three robust FIR (PACS) counterparts at z~2.2 in GOODS-North. This demonstrates the possibility of finding dust emission in LAEs even at higher redshifts.

We present first insights into the far-IR properties for a sample of IRAC and MIPS-24um detected Lyman Break Galaxies (LBGs) at z ~ 3, as derived from observations in the northern field of the Great Observatories Origins Survey (GOODS-N) carried out with the PACS instrument on board the Herschel Space Observatory. Although none of our galaxies are detected by Herschel, we employ a stacking technique to construct, for the first time, the average spectral energy distribution of infrared luminous LBGs from UV to radio wavelengths. We derive a median IR luminosity of L_{IR} = 1.6 x 10^12 Lo, placing the population in the class of ultra luminous infrared galaxies (ULIRGs). Complementing our study with existing multi-wavelength data, we put constraints on the dust temperature of the population and find that for their L_{IR}, MIPS-LBGs are warmer than submm-luminous galaxies while they fall in the locus of the L_{IR}-T_{d} relation of the local ULIRGs. This, along with estimates based on the average SED, explains the marginal detection of LBGs in current sub-mm surveys and suggests that these latter studies introduce a bias towards the detection of colder ULIRGs in the high-z universe, while missing high-z ULIRGS with warmer dust.

We present a search for Herschel-PACS counterparts of dust-obscured, high-redshift objects previously selected at submillimeter and millimeter wavelengths in the Great Observatories Origins Deep Survey North field. We detect 22 of 56 submillimeter galaxies (39%) with a SNR of >=3 at 100 micron down to 3.0 mJy, and/or at 160 micron down to 5.7 mJy. The fraction of SMGs seen at 160 micron is higher than that at 100 micron. About 50% of radio-identified SMGs are associated with PACS sources. We find a trend between the SCUBA/PACS flux ratio and redshift, suggesting that these flux ratios could be used as a coarse redshift indicator. PACS undetected submm/mm selected sources tend to lie at higher redshifts than the PACS detected ones. A total of 12 sources (21% of our SMG sample) remain unidentified and the fact that they are blank fields at Herschel-PACS and VLA 20 cm wavelength may imply higher redshifts for them than for the average SMG population (e.g., z>3-4). The Herschel-PACS imaging of these dust-obscured starbursts at high-redshifts suggests that their far-infrared spectral energy distributions have significantly different shapes than template libraries of local infrared galaxies.

We use deep observations taken with the Photodetector Array Camera and Spectrometer (PACS), on board the Herschel satellite as part of the PACS evolutionary probe (PEP) guaranteed project along with submm ground-based observations to measure the dust mass of a sample of high-z submillimeter galaxies (SMGs). We investigate their dust content relative to their stellar and gas masses, and compare them with local star-forming galaxies. High-z SMGs are dust rich, i.e. they have higher dust-to-stellar mass ratios compared to local spiral galaxies (by a factor of 30) and also compared to local ultraluminous infrared galaxies (ULIRGs, by a factor of 6). This indicates that the large masses of gas typically hosted in SMGs have already been highly enriched with metals and dust. Indeed, for those SMGs whose gas mass is measured, we infer dust-to-gas ratios similar or higher than local spirals and ULIRGs. However, similarly to other strongly star-forming galaxies in the local Universe and at high-z, SMGs are characterized by gas metalicities lower (by a factor of a few) than local spirals, as inferred from their optical nebular lines, which are generally ascribed to infall of metal-poor gas. This is in contrast with the large dust content inferred from the far-IR and submm data. In short, the metalicity inferred from the dust mass is much higher (by more than an order of magnitude) than that inferred from the optical nebular lines. We discuss the possible explanations of this discrepancy and the possible implications for the investigation of the metalicity evolution at high-z.

Gravitational lensing by massive galaxy clusters allows study of the population of intrinsically faint infrared galaxies that lie below the sensitivity and confusion limits of current infrared and submillimeter telescopes. We present ultra-deep PACS 100 and 160 microns observations toward the cluster lens Abell 2218, to penetrate the Herschel confusion limit. We derive source counts down to a flux density of 1 mJy at 100 microns and 2 mJy at 160 microns, aided by strong gravitational lensing. At these levels, source densities are 20 and 10 beams/source in the two bands, approaching source density confusion at 160 microns. The slope of the counts below the turnover of the Euclidean-normalized differential curve is constrained in both bands and is consistent with most of the recent backwards evolutionary models. By integrating number counts over the flux range accessed by Abell 2218 lensing (0.94-35 mJy at 100 microns and 1.47-35 mJy at 160 microns, we retrieve a cosmic infrared background (CIB) surface brightness of ~8.0 and ~9.9 nW m^-2 sr^-1, in the respective bands. These values correspond to 55% (+/- 24%) and 77% (+/- 31%) of DIRBE direct measurements. Combining Abell 2218 results with wider/shallower fields, these figures increase to 62% (+/- 25%) and 88% (+/- 32%) CIB total fractions, resolved at 100 and 160 microns, disregarding the high uncertainties of DIRBE absolute values.

The Photodetector Array Camera and Spectrometer (PACS) is one of the three science instruments on ESA's far infrared and submillimetre observatory. It employs two Ge:Ga photoconductor arrays (stressed and unstressed) with 16x25 pixels, each, and two filled silicon bolometer arrays with 16x32 and 32x64 pixels, respectively, to perform integral-field spectroscopy and imaging photometry in the 60-210\mu\ m wavelength regime. In photometry mode, it simultaneously images two bands, 60-85\mu\ m or 85-125\mu\m and 125-210\mu\ m, over a field of view of ~1.75'x3.5', with close to Nyquist beam sampling in each band. In spectroscopy mode, it images a field of 47"x47", resolved into 5x5 pixels, with an instantaneous spectral coverage of ~1500km/s and a spectral resolution of ~175km/s. We summarise the design of the instrument, describe observing modes, calibration, and data analysis methods, and present our current assessment of the in-orbit performance of the instrument based on the Performance Verification tests. PACS is fully operational, and the achieved performance is close to or better than the pre-launch predictions.

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

The star formation rate (SFR) is a key parameter in the study of galaxy evolution. The accuracy of SFR measurements at z~2 has been questioned following a disagreement between observations and theoretical models. The latter predict SFRs at this redshift that are typically a factor 4 or more lower than the measurements. We present star-formation rates based on calorimetric measurements of the far-infrared (FIR) luminosities for massive 1.5<z<2.5, normal star-forming galaxies (SFGs), which do not depend on extinction corrections and/or extrapolations of spectral energy distributions. The measurements are based on observations in GOODS-N with the Photodetector Array Camera & Spectrometer (PACS) onboard Herschel, as part of the PACS Evolutionary Probe (PEP) project, that resolve for the first time individual SFGs at these redshifts at FIR wavelengths. We compare FIR-based SFRs to the more commonly used 24 micron and UV SFRs. We find that SFRs from 24 micron alone are higher by a factor of ~4-7.5 than the true SFRs. This overestimation depends on luminosity: gradually increasing for log L(24um)>12.2 L_sun. The SFGs and AGNs tend to exhibit the same 24 micron excess. The UV SFRs are in closer agreement with the FIR-based SFRs. Using a Calzetti UV extinction correction results in a mean excess of up to 0.3 dex and a scatter of 0.35 dex from the FIR SFRs. The previous UV SFRs are thus confirmed and the mean excess, while narrowing the gap, is insufficient to explain the discrepancy between the observed SFRs and simulation predictions.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

We present Herschel-PACS observations of rest-frame mid-infrared and far-infrared spectral line emissions from two lensed, ultra-luminous infrared galaxies at high redshift: MIPS J142824.0+352619 (MIPS J1428), a starburst-dominated system at z = 1.3, and IRAS F10214+4724 (F10214), a source at z = 2.3 hosting both star-formation and a luminous AGN. We have detected [OI]63 micron and [OIII]52 micron in MIPS J1428, and tentatively [OIII]52 micron in F10214. Together with the recent ZEUS-CSO [CII]158 micron detection in MIPS J1428 we can for the first time combine [OI], [CII] and far-IR (FIR) continuum measurements for photo-dissociation (PDR) modeling of an ultra-luminous (L_IR > 10^12 L_sun) star forming galaxy at the peak epoch of cosmic star formation. We find that MIPS J1428, contrary to average local ULIRGs, does not show a deficit in [OI] relative to FIR. The combination of far-UV flux G_0 and gas density n (derived from the PDR models), as well as the star formation efficiency (derived from CO and FIR) is similar to normal or starburst galaxies, despite the high infrared luminosity of this system. In contrast, F10214 has stringent upper limits on [OIV] and [SIII], and an [OIII]/FIR ratio at least an order of magnitude lower than local starbursts or AGN, similar to local ULIRGs.

Sensitive Herschel far-infrared observations can break degeneracies that were inherent to previous studies of star formation in high-z AGN hosts. Combining PACS 100 and 160um observations of the GOODS-N field with 2Msec Chandra data, we detect ~20% of X-ray AGN individually at >3sig. The host far-infrared luminosity of AGN with L2-10~10^43erg/s increases with redshift by an order of magnitude from z=0 to z~1. In contrast, there is little dependence of far-infrared luminosity on AGN luminosity, for L2-10<~10^44erg/s AGN at z>~1. We do not find a dependence of far-infrared luminosity on X-ray obscuring column, for our sample which is dominated by L2-10<10^44erg/s AGN. In conjunction with properties of local and luminous high-z AGN, we interpret these results as reflecting the interplay between two paths of AGN/host coevolution. A correlation of AGN luminosity and host star formation is traced locally over a wide range of luminosities and also extends to luminous high z AGN. This correlation reflects an evolutionary connection, likely via merging. For lower AGN luminosities, star formation is similar to that in non-active massive galaxies and shows little dependence on AGN luminosity. The level of this secular, non-merger driven star formation increasingly dominates over the correlation at increasing redshift.

In this first paper on the results of our Herschel PACS survey of local Ultraluminous Infrared Galaxies (ULIRGs), as part of our SHINING survey of local galaxies, we present far-infrared spectroscopy of Mrk 231, the most luminous of the local ULIRGs, and a type 1 broad absorption line AGN. For the first time in a ULIRG, all observed far-infrared fine-structure lines in the PACS range were detected and all were found to be deficient relative to the far infrared luminosity by 1 - 2 orders of magnitude compared with lower luminosity galaxies. The deficits are similar to those for the mid-infrared lines, with the most deficient lines showing high ionization potentials. Aged starbursts may account for part of the deficits, but partial covering of the highest excitation AGN powered regions may explain the remaining line deficits. A massive molecular outflow, discovered in OH and 18OH, showing outflow velocities out to at least 1400 km/sec, is a unique signature of the clearing out of the molecular disk that formed by dissipative collapse during the merger. The outflow is characterized by extremely high ratios of 18O / 16O suggestive of interstellar medium processing by advanced starbursts.

We present Herschel/PACS spectroscopic maps of the dwarf galaxy NC4214 observed in 6 far infrared fine-structure lines: [C II] 158mu, [O III] 88mu, [O I] 63mu, [O I] 146mu, [N II] 122mu, and [N II] 205mu. The maps are sampled to the full telescope spatial resolution and reveal unprecedented detail on ~ 150 pc size scales. We detect [C II] emission over the whole mapped area, [O III] being the most luminous FIR line. The ratio of [O III]/[C II] peaks at about 2 toward the sites of massive star formation, higher than ratios seen in dusty starburst galaxies. The [C II]/CO ratios are 20 000 to 70 000 toward the 2 massive clusters, which are at least an order of magnitude larger than spiral or dusty starbursts, and cannot be reconciled with single-slab PDR models. Toward the 2 massive star-forming regions, we find that L[CII] is 0.5 to 0.8% of the LTIR . All of the lines together contribute up to 2% of LTIR . These extreme findings are a consequence of the lower metallicity and young, massive-star formation commonly found in dwarf galaxies. These conditions promote large-scale photodissociation into the molecular reservoir, which is evident in the FIR line ratios. This illustrates the necessity to move to multiphase models applicable to star-forming clusters or galaxies as a whole.

We exploit the deepest existing far-infrared (FIR) data obtained so far by Herschel at 100 and 160 um in the GOODS-N, as part of the PACS Evolutionary Probe (PEP) survey, to derive for the first time the evolution of the rest-frame 60-um, 90-um, and total IR luminosity functions (LFs) of galaxies and AGNs from z=0 to unprecedented high redshifts (z~2-3). The PEP LFs were computed using the 1/Vmax method. The FIR sources were classified by means of a detailed broad- band SED-fitting analysis and spectral characterisation. Based on the best-fit model results, k-correction and total IR (8-1000 um) luminosity were obtained for each source. LFs (monochromatic and total) were then derived for various IR populations separately in different redshift bins and compared to backward evolution model predictions. We detect strong evolution in the LF to at least z~2. Objects with SEDs similar to local spiral galaxies are the major contributors to the star formation density (SFD) at z< 0.3, then, as redshift increases, moderate SF galaxies - most likely containing a low-luminosity AGN - start dominating up to z ~= 1.5. At >1.5 the SFD is dominated by the contributions of starburst galaxies. In agreement with previous findings, the comoving IR LD derived from our data evolves approximately as (1 + z)^(3.8+/-0.3) up to z~1, there being some evidence of flattening up to z~2.

We use deep observations obtained with the Photodetector Array Camera and Spectrometer (PACS) onboard the Herschel space observatory to study the far-infrared (FIR) properties of submm and optically faint radio galaxies (SMGs and OFRGs). From literature we compiled a sample of 35 securely identified SMGs and nine OFRGs located in the GOODS-N and the A2218 fields. This sample is cross-matched with our PACS 100 um and 160 um multi-wavelength catalogs. About half of the galaxies in our sample are detected with PACS. The dust temperatures and the infrared luminosities of our galaxies are derived by fitting their PACS and SCUBA 850 um (only the upper limits for the OFRGs) flux densities with a single modified (beta=1.5) black body function. The median dust temperature of our SMG sample is T=36+/-8K while for our OFRG sample it is T=47+/-3K. For both samples, median dust temperatures derived from Herschel data agree well with previous estimates. In particular, Chapman et al. (2005) found a dust temperature of T=36+/-7K for a large sample of SMGs assuming the validity of the FIR/radio correlation. The agreement between our studies confirms that the local FIR/radio correlation effectively holds at high redshift even though we find <q>=2.17+/-0.19, a slightly lower value than that observed in local systems. The median IR luminosities of SMGs and OFRGs are 4.6*10^12 Lsun and 2.6*10^12 Lsun, respectively. We note that for both samples the IR luminosity estimates from the radio part of the spectral energy distribution are accurate, while estimates from the mid-IR are considerably (x3) more uncertain. Our observations confirm the remarkably high luminosities of SMGs and thus imply median star-formation rates of 960Msun yr^-1 for SMGs with S(850um)>5mJy and 460Msun yr^-1 for SMGs with S(850um)>2mJy, assuming a Chabrier IMF and no dominant AGN contribution to the far-infrared luminosity.

The constituents of the cosmic IR background (CIB) are studied at its peak wavelengths (100 and 160 um) by exploiting Herschel/PACS observations of the GOODS-N, Lockman Hole, and COSMOS fields in the PACS Evolutionary Probe (PEP) guaranteed-time survey. The GOODS-N data reach 3 sigma depths of ~3.0 mJy at 100 um and ~5.7 mJy at 160 um. At these levels, source densities are 40 and 18 beams/source, respectively, thus hitting the confusion limit at 160 um. Differential number counts extend from a few mJy up to 100-200 mJy, and are approximated as a double power law, with the break lying between 5 and 10 mJy. The available ancillary information allows us to split number counts into redshift bins. At z<=0.5 we isolate a class of luminous sources (L(IR)~1e11 Lsun), whose SEDs resemble late-spiral galaxies, peaking at ~130 um restframe and significantly colder than what is expected on the basis of pre-Herschel models. By integrating number counts over the whole covered flux range, we obtain a surface brightness of 6.36 +/- 1.67 and 6.58 +/-1.62 [nW m^-2 sr^-1] at 100 and 160 um, resolving ~45% and ~52% of the CIB, respectively. When stacking 24 um sources, the inferred CIB lies within 1.1 sigma and 0.5 sigma from direct measurements in the two bands, and fractions increase to 50% and 75%.Most of this resolved CIB fraction was radiated at z<=1.0, with 160 um sources found at higher redshift than 100 um ones.

We exploit deep observations of the GOODS-N field taken with PACS, on board of Herschel, as part of the PEP guaranteed time, to study the link between star formation and stellar mass in galaxies to z~2. Starting from a stellar mass-selected sample of ~4500 galaxies with mag[4.5mu]<23 (AB), we identify ~350 objects with a PACS detection at 100 or 160mu and ~1500 with only Spitzer 24 mu counterpart. Stellar masses and total IR luminosities (LIR) are estimated by fitting the SEDs. Consistently with other Herschel results, we find that LIR based only on 24 mu data is overestimated by a median factor ~1.8 at z~2, whereas it is underestimated (with our approach) up to a factor ~1.6 at 0.5<z<1.0. We then exploit this calibration to correct LIR based on the MIPS fluxes. These results clearly show how Herschel is fundamental to constrain LIR, and hence the SFR, of high redshift galaxies. Using the galaxies detected with PACS (and/or MIPS), we investigate the existence and evolution of the relations between the star formation rate (SFR), the specific star formation rate (SSFR=SFR/mass) and the stellar mass. Moreover, in order to avoid selection effects, we also repeat this study through a stacking analysis on the PACS images to fully exploit the far-IR information also for the Herschel and Spitzer undetected subsamples. We find that the SSFR-mass relation steepens with redshift, being almost flat at z<1.0 and reaching a slope of alpha=-0.50^(+0.13)_(-0.16) at z~2, at odds with recent works based on radio-stacking analysis at the same redshift. The mean SSFR of galaxies increases with redshift, by a factor ~15 for massive M>10^(11) (M_sun) galaxies from z=0 to z=2, and seems to flatten at z>1.5 in this mass range. Moreover, the most massive galaxies have the lowest SSFR at any z, implying that they have formed their stars earlier and more rapidly than their low mass counterparts.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

An analysis of data from the Spitzer Space Telescope, Hubble Space Telescope, Chandra X-ray Observatory, and AKARI Infrared Astronomy Satellite is presented for the z=0.036 merging galaxy system II Zw 096 (CGCG 448-020). Because II Zw 096 has an infrared luminosity of log(L_IR/L_sun) = 11.94, it is classified as a Luminous Infrared Galaxy (LIRG), and was observed as part of the Great Observatories All-sky LIRG Survey (GOALS). The Spitzer data suggest that 80% of the total infrared luminosity comes from an extremely compact, red source not associated with the nuclei of the merging galaxies. The Spitzer mid-infrared spectra indicate no high-ionization lines from a buried active galactic nucleus in this source. The strong detection of the 3.3 micron and 6.2 micron PAH emission features in the AKARI and Spitzer spectra also implies that the energy source of II Zw 096 is a starburst. Based on Spitzer infrared imaging and AKARI near-infrared spectroscopy, the star formation rate is estimated to be 120 M_sun/yr and > 45 M_sun/yr, respectively. Finally, the high-resolution B, I, and H-band images show many star clusters in the interacting system. The colors of these clusters suggest at least two populations - one with an age of 1-5 Myr and one with an age of 20-500 Myr, reddened by 0-2 magnitudes of visual extinction. The masses of these clusters span a range between 10^6-10^8 M_sun. This starburst source is reminiscent of the extra-nuclear starburst seen in NGC 4038/9 (the Antennae Galaxies) and Arp 299 but approximately an order of magnitude more luminous than the Antennae. The source is remarkable in that the off-nuclear infrared luminosity dominates the enitre system.

We present new low resolution Spitzer mid-infrared spectroscopy of a sample of 20 ROSAT selected local Narrow Line Seyfert 1 galaxies (NLS1s). We detect strong AGN continuum in all and clear PAH emission in 70% of the sources. The 6.2 micron PAH luminosity spans three orders of magnitudes, from ~10^(39) erg/s to ~10^(42) erg/s providing strong evidence for intense ongoing star formation in the circumnuclear regions of these sources. Using the IRS/Spitzer archive we gather a large number of additional NLS1s and their broad line counterparts (BLS1s) and constructed NLS1 and BLS1 sub-samples to compare them in various ways. The comparison shows a clear separation according to FWHM(H_beta) such that objects with narrower broad H_beta lines are the strongest PAH emitters. We test this division in various ways trying to remove biases due to luminosity and aperture size. Specifically, we find that star formation activity around NLS1 AGN is larger than around BLS1 of the same AGN luminosity. The above result seems to hold over the entire range of distance and luminosity. Moreover the star formation rate is higher in low black hole mass and high L/L_Edd systems indicating that black hole growth and star formation are occurring simultaneously.

Empirical "backward" galaxy evolution models for infrared bright galaxies are constrained using multi-band infrared surveys. We developed a new Monte-Carlo algorithm for this task, implementing luminosity dependent distribution functions for the galaxies' infrared spectral energy distributions (SEDs) and for the AGN contribution, allowing for evolution of these quantities. The adopted SEDs take into account the contributions of both starbursts and AGN to the infrared emission, for the first time in a coherent treatment rather than invoking separate AGN and star-forming populations. In the first part of the paper we consider the quantification of the AGN contribution for local universe galaxies, as a function of total infrared luminosity. It is made using a large sample of LIRGs and ULIRGs for which mid-infrared spectra are available in the Spitzer archive. In the second part we present the model. Our best-fit model adopts very strong luminosity evolution, $L=L_0(1+z)^{3.4}$, up to $z=2.3$, and density evolution, $\rho=\rho_0(1+z)^2$, up to $z=1$, for the population of infrared galaxies. At higher $z$, the evolution rates drop as $(1+z)^{-1}$ and $(1+z)^{-1.5}$ respectively. To reproduce mid-infrared to submillimeter number counts and redshift distributions, it is necessary to introduce both an evolution in the AGN contribution and an evolution in the luminosity-temperature relation. Our models are in plausible agreement with current photometry-based estimates of the typical AGN contribution as a function of mid-infrared flux, and well placed to be compared to upcoming Spitzer spectroscopic results. As an example of future applications, we use our best-fitting model to make predictions for surveys with Herschel.

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)

The Great Observatories All-sky LIRG Survey (GOALS) combines data from NASA's Spitzer, Chandra, Hubble and GALEX observatories, together with ground-based data into a comprehensive imaging and spectroscopic survey of over 200 low redshift Luminous Infrared Galaxies (LIRGs). The LIRGs are a complete subset of the IRAS Revised Bright Galaxy Sample (RBGS). The LIRGs targeted in GOALS span the full range of nuclear spectral types defined via traditional optical line-ratio diagrams as well as interaction stages. They provide an unbiased picture of the processes responsible for enhanced infrared emission in galaxies in the local Universe. As an example of the analytic power of the multi-wavelength GOALS dataset, we present data for the interacting system VV 340 (IRAS F14547+2449). Between 80-95% of the total far-infrared emission (or about 5E11 solar luminosities) originates in VV 340 North. While the IRAC colors of VV 340 North and South are consistent with star-forming galaxies, both the Spitzer IRS and Chandra ACIS data indicate the presence of a buried AGN in VV 340 North. The GALEX far and near-UV fluxes imply a extremely large infrared "excess" (IRX) for the system (IR/FUV = 81) which is well above the correlation seen in starburst galaxies. Most of this excess is driven by VV 340 N, which alone has an IR excess of nearly 400. The VV 340 system seems to be comprised of two very different galaxies - an infrared luminous edge-on galaxy (VV 340 North) that dominates the long-wavelength emission from the system and which hosts a buried AGN, and a face-on starburst (VV 340 South) that dominates the short-wavelength emission.

We report the results from a comprehensive study of 74 ultraluminous infrared galaxies (ULIRGs) and 34 Palomar-Green (PG) quasars within z ~ 0.3$ observed with the Spitzer Infrared Spectrograph (IRS). The contribution of nuclear activity to the bolometric luminosity in these systems is quantified using six independent methods that span a range in wavelength and give consistent results within ~ +/-10-15% on average. The average derived AGN contribution in ULIRGs is ~35-40%, ranging from ~15-35% among "cool" (f_25/f_60 =< 0.2) optically classified HII-like and LINER ULIRGs to ~50 and ~75% among warm Seyfert 2 and Seyfert 1 ULIRGs, respectively. This number exceeds ~80% in PG QSOs. ULIRGs fall in one of three distinct AGN classes: (1) objects with small extinctions and large PAH equivalent widths are highly starburst-dominated; (2) systems with large extinctions and modest PAH equivalent widths have larger AGN contributions, but still tend to be starburst-dominated; and (3) ULIRGs with both small extinctions and small PAH equivalent widths host AGN that are at least as powerful as the starbursts. The AGN contributions in class 2 ULIRGs are more uncertain than in the other objects, and we cannot formally rule out the possibility that these objects represent a physically distinct type of ULIRGs. A morphological trend is seen along the sequence (1)-(2)-(3), in general agreement with the standard ULIRG - QSO evolution scenario and suggestive of a broad peak in extinction during the intermediate stages of merger evolution. However, the scatter in this sequence, implies that black hole accretion, in addition to depending on the merger phase, also has a strong chaotic/random component, as in local AGN. (abridged)

We present Spitzer rest-frame mid-infrared spectroscopy of twelve z~2 mm-bright type 1 QSOs, selected from unlensed and lensed QSO samples and covering a range of AGN optical luminosities L_5100=10^45 to 10^47 erg/s. On top of the AGN continuum, we detect PAH emission from luminous star formation in nine objects individually as well as in the composite spectrum for the full sample. PAH luminosity and rest frame far-infrared luminosity correlate and extend the similar correlation for lower luminosity local QSOs. This provides strong evidence for intense star formation in the hosts of these mm-bright QSOs, sometimes exceeding 1000 Msun/yr and dominating their rest frame far-infrared emission. The PAH-based limit on star formation rates is lower for luminous z~2 QSOs that are not preselected for their mm emission. Partly dependent on systematic changes of the AGN dust covering factor and the dust spectral energy distribution of the AGN proper, the spectral energy distributions of mm-faint high-z QSOs may be AGN dominated out to rest frame far-infrared wavelengths. Towards the most luminous high-z QSOs, there is a flattening of the relation between star formation and AGN luminosity that is observed for lower redshift QSOs. No QSO in our sample has a PAH-measured star formation rate in excess of 3000 Msun/yr.

This paper addresses the origin of the silicate emission observed in PG QSOs, based on observations with the Spitzer Space Telescope. Scenarios based on the unified model suggest that silicate emission in AGN arises mainly from the illuminated faces of the clouds in the torus at temperatures near sublimation. However, detections of silicate emission in Type 2 QSOs, and the estimated cool dust temperatures, argue for a more extended emission region.To investigate this issue we present the mid-infrared spectra of 23 QSOs. These spectra, and especially the silicate emission features at ~10 and ~18 mu can be fitted using dusty narrow line region (NLR) models and a combination of black bodies. The bolometric luminosities of the QSOs allow us to derive the radial distances and covering factors for the silicate-emitting dust. The inferred radii are 100-200 times larger than the dust sublimation radius, much larger than the expected dimensions of the inner torus. Our QSO mid-IR spectra are consistent with the bulk of the silicate dust emission arising from the dust in the innermost parts of the NLR.

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.

The impact of active galactic nuclei on low-ionization nuclear emission-line regions (LINERs) remains a vigorous field of study. We present preliminary results from a study of the mid-infrared atomic emission lines of LINERs with the Spitzer Space Telescope. We assess the ubiquity and properties of AGN in LINERs using this data. We discuss what powers the mid-infrared emission lines and conclude that the answer depends unsurprisingly on the emission line ionization state and, more interestingly, on the infrared luminosity.

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 Spitzer IRS low resolution, mid-IR spectra of a sample of 25 high luminosity QSOs at 2<z<3.5. When combined with archival IRS observations of local, low luminosity type-I AGNs, the sample spans five orders of magnitude in luminosity. We find that the continuum dust thermal emission at lambda(rest)=6.7um is correlated with the optical luminosity, following the non-linear relation L(6.7um) propto L(5100A)^0.82. We also find an anti correlation between the ratio L(6.7um)/L(5100A) and the [OIII]5007A line luminosity. These effects are interpreted as a decreasing covering factor of the circumnuclear dust as a function of luminosity. Such a result is in agreement with the decreasing fraction of absorbed AGNs as a function of luminosity recently found in various surveys. We clearly detect the silicate emission feature in the average spectrum, but also in four individual objects. These are the Silicate emission in the most luminous objects obtained so far. When combined with the silicate emission observed in local, low luminosity type-I AGNs, we find that the silicate emission strength is correlated with luminosity. The silicate strength of all type-I AGNs also follows a positive correlation with the black hole mass and with the accretion rate. The Polycyclic Aromatic Hydrocarbon (PAH) emission features, expected from starburst activity, are not detected in the average spectrum of luminous, high-z QSOs. The upper limit inferred from the average spectrum points to a ratio between PAH luminosity and QSO optical luminosity significantly lower than observed in lower luminosity AGNs, implying that the correlation between star formation rate and AGN power saturates at high luminosities.

We report the first detection of the 6.2micron and 7.7micron infrared `PAH' emission features in the spectrum of a high redshift QSO, from the Spitzer-IRS spectrum of the Cloverleaf lensed QSO (H1413+117, z~2.56). The ratio of PAH features and rest frame far-infrared emission is the same as in lower luminosity star forming ultraluminous infrared galaxies and in local PG QSOs, supporting a predominantly starburst nature of the Cloverleaf's huge far-infrared luminosity (5.4E12 Lsun, corrected for lensing). The Cloverleaf's period of dominant QSO activity (Lbol ~ 7E13 Lsun) is coincident with an intense (star formation rate ~1000 Msun/yr) and short (gas exhaustion time ~3E7yr) star forming event.

We present rest frame mid-infrared spectroscopy of a sample of 13 submillimeter galaxies, obtained using the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope. The sample includes exclusively bright objects from blank fields and cluster lens assisted surveys that have accurate interferometric positions. We find that the majority of spectra are well fitted by a starburst template or by the superposition of PAH emission features and a weak mid-infrared continuum, the latter a tracer of Active Galactic Nuclei (including Compton-thick ones). We obtain mid-infrared spectroscopic redshifts for all nine sources detected with IRS. For three of them the redshifts were previously unknown. The median value of the redshift distribution is z~2.8 if we assume that the four IRS non-detections are at high redshift. The median for the IRS detections alone is z~2.7. Placing the IRS non-detections at similar redshift would require rest frame mid-IR obscuration larger than is seen in local ULIRGs. The rest frame mid-infrared spectra and mid- to far-infrared spectral energy distributions are consistent with those of local ultraluminous infrared galaxies, but scaled-up further in luminosity. The mid-infrared spectra support the scenario that submillimeter galaxies are sites of extreme star formation, rather than X-ray-obscured AGN, and represent a critical phase in the formation of massive galaxies.

We report results from the first mid-infrared spectroscopic study of a comprehensive sample of 33 LINERs, observed with the Spitzer Space Telescope. We compare the properties of two different LINER populations: infrared-faint LINERs, with LINER emission arising mostly in compact nuclear regions, and infrared-luminous LINERs, which often show spatially extended (non-AGN) LINER emission. We show that these two populations can be easily distinguished by their mid-infrared spectra in three different ways: (i) their mid-IR spectral energy distributions (SEDs), (ii) the emission features of polycyclic aromatic hydrocarbons (PAHs), and (iii) various combinations of IR fine-structure line ratios. IR-luminous LINERs show mid-IR SEDs typical of starburst galaxies, while the mid-IR SEDs of IR-faint LINERs are much bluer. PAH flux ratios are significantly different in the two groups. Fine structure emission lines from highly excited gas, such as [O IV], are detected in both populations, suggesting the presence of an additional AGN also in a large fraction of IR-bright LINERs, which contributes little to the combined mid-IR light. The two LINER groups occupy different regions of mid-infrared emission-line excitation diagrams. The positions of the various LINER types in our diagnostic diagrams provide important clues regarding the power source of each LINER type. Most of these mid-infrared diagnostics can be applied at low spectral resolution, making AGN- and starburst-excited LINERs distinguishable also at high redshifts.

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 have performed a spectroscopic study of 7 Type 2 QSOs using the mid-infrared spectrometer IRS on board the Spitzer Space Telescope. These are (to our knowledge) the first mid-IR spectra of X-ray selected QSO2s taken. The objects have been selected according to their high intrinsic luminosities and column densities in the X-rays. Their spectra strongly differ from template spectra of Type 2 AGN at lower luminosities. They do not exhibit strong PAH dust emission features from circum-nuclear star forming regions, typical for lower luminosity Type 2 Seyfert galaxies or other previously used QSO2 templates, such as the (Ultra)luminous Infrared Galaxy ((U)LIRG) NGC 6240. They also do not show the ice and silicate absorption features of highly luminous but deeply embedded compact nuclei seen in some ULIRGs. Instead they reveal a relatively featureless, rising continuum similar to luminous Type 1 AGN. We also find evidence for a 10 micron silicate feature in emission. Models of dusty tori in the AGN unification scenario predict this only for Type 1 AGN. The ratio of the AGN continuum luminosity at 6 micron to the absorption corrected 2-10keV X-ray AGN luminosity is very similar to that found in Seyfert galaxies. X-ray selected QSO2s are thus characterized by powerful AGN in hosts with a luminosity due to star formation < 1e11 L_sol. The dominance of the AGN light in the mid-IR spectra of QSO2s together with their flatter spectral energy distributions (SEDs) places important constraints on models of the cosmic infrared background and of the star formation history of the universe.

We present low resolution near-infrared spectroscopy of an unbiased sample of 24 ultraluminous infrared galaxies (ULIRGs), selected from samples previously observed spectroscopically in the mid-infrared with the Infrared Space Observatory (ISO). Qualitatively, the near-infrared spectra resemble those of starbursts. Only in one ULIRG, IRAS 04114-5117E, do we find spectroscopic evidence for AGN activity. The spectroscopic classification in the near-infrared is in very good agreement with the mid-infrared one. For a subset of our sample for which extinction corrections can be derived from Pa-alpha and Br-gamma, we find rather high Pa-alpha luminosities, in accordance with the powering source of these galaxies being star formation.[Fe] emission is strong in ULIRGs and may be linked to starburst and superwind activity. Additionally, our sample includes two unusual objects. The first, IRAS F00183-7111, exhibits extreme [Fe] emission and the second, IRAS F23578-5307, is according to our knowledge one of the most luminous infrared galaxies in H2 rotation-vibration emission.

We report the first detection of ~10 and ~18 micron silicate dust emissions in a low-luminosity active galactic nucleus (AGN), obtained in Spitzer-IRS 7-37 micron spectroscopy of the Type 1 LINER galaxy NGC3998. Silicate emissions in AGN have only recently been detected in several quasars. Our detection counters suggestions that silicate emissions are present only in the most luminous AGN. The silicate features may be signatures of a dusty ``obscuring torus'' viewed face-on as postulated for Type 1 AGN. However, the apparently cool (~200 K) dust is inconsistent with theoretical expectations of much hotter torus walls. Furthermore, not all Type 1 objects are silicate emission sources. Alternatively, the silicate emission may originate in dust not directly associated with a torus. We find that the long-wavelength (>20 micron) tail of the emission in NGC3998 is significantly weaker than in the sample of bright QSOs recently presented by Hao et al. The 10 micron profile in our NGC3998 spectrum is inconsistent with ``standard'' silicate ISM dust. This may indicate differences in the dust composition, grain size distribution, or degree of crystallization. The differences between NGC3998, QSOs, and Galactic templates suggest that there are significant environmental variations.

Using the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope, we have obtained rest frame mid-infrared spectroscopy of two bright submillimeter galaxies. SMMJ02399-0136 at z=2.81 shows a superposition of PAH emission features and a mid-infrared continuum, indicating significant and roughly equal contributions to its bolometric luminosity from star formation and from a Compton-thick AGN. We derive a new redshift of z=2.80 for MMJ154127+6616 from the IRS spectrum and find this object is dominated by starburst PAH emission. The rest frame mid- to far-infrared spectral energy distributions are consistent with these submillimeter galaxies being scaled up versions of local ultraluminous infrared galaxies. The mid-infrared spectra support the scenario that submillimeter galaxies are sites of extreme star formation and represent a key phase in the formation of massive galaxies.

(ABRIDGED) We have undertaken a programme of 3 micron narrowband imaging of a sample of nearby template galaxies covering a wide range of metallicity, star formation activity, and nuclear activity. In the present paper we present first results: high spatial resolution images of PAH emission and the adjacent continuum emission from the central regions of the nearby starburst galaxies NGC253 and NGC1808. Globally, the PAH emission is seen to peak on the central starburst regions of both sources. On smaller scales, however, we see no general spatial correlation or anti-correlation between the PAH emission and the location of sites of recent star formation, suggesting that the degree to which PAH emission traces starburst activity is more complicated than previously hypothesized based on results from data with lower spatial resolution. We do find spatial correlations, though, when we consider the PAH-to-continuum ratio, which is low at the positions of known super star clusters in NGC1808 as well as at the position of the IR peak in NGC253. We take this to imply a decrease in the efficiency of PAH emission induced by the star formation, caused by mechanical energy input into the ISM or photoionisation or photodissociation of the PAH molecules. All three hypotheses are discussed. In addition, for the first time we present observations of PAH emission in the superwind of a starbursting system (NGC253), providing strong support that winds are heavily mass-loaded and entrain substantial amounts of ambient ISM. We have also found a plausible connection between observed NaD absorption, H2, and PAHs above the plane of NGC253. This observation has important implications for enriching galaxy halos and possibly the intergalactic medium with small dust grains.

In order to constrain the properties of dense and warm gas around active galactic nuclei, we have searched Infrared Space Observatory spectra of local active galactic nuclei for the signature of the 4.7micron fundamental ro-vibrational band of carbon monoxide. Low resolution spectra of 31 AGN put upper limits on the presence of wide absorption bands corresponding to absorption by large columns of warm and dense gas against the nuclear dust continuum. High resolution (R~2500) spectra of NGC 1068 detect no significant absorption or emission in individual lines, to a 3sigma limit of 7% of the continuum. The limits set on CO absorption in local AGN are much lower than the recent Spitzer Space Telescope detection of strong CO absorption by dense and warm gas in the obscured ultraluminous infrared galaxy IRAS F00183-7111, despite evidence for dense material on parsec scales near an AGN in both types of objects. This suggests that such deep absorptions are not intimately related to the obscuring `torus' material invoked in local AGN, but rather are a signature of the peculiar conditions in the circumnuclear region of highly obscured infrared galaxies like IRAS F00183-7111. They may reflect fully covered rather than torus geometries.

We present Infrared Space Observatory mid- and far-infrared spectroscopy of the merging galaxy NGC 6240, an object presenting many aspects of importance for the role of star formation and AGN activity in [ultra]luminous infrared galaxies. The mid-infrared spectrum shows starburst indicators in the form of low excitation fine-structure line emission and aromatic `PAH' features. A strong high excitation [OIV] line is observed which most likely originates in the Narrow Line Region of an optically obscured AGN. NGC 6240 shows extremely powerful emission in the pure rotational lines of molecular hydrogen. We argue that this emission is mainly due to shocks in its turbulent central gas component and its starburst superwind. The total shock cooling in infrared emission lines accounts for ~0.6% of the bolometric luminosity, mainly through rotational H_2 emission and the [OI] 63micron line. We analyse several ways of estimating the luminosities of the starburst and the AGN in NGC 6240 and suggest that the contributions to its bolometric luminosity are most likely in the range 50-75% starburst and 25-50% AGN.

We present a study of the [CII] 157.74 um fine-structure line in a sample of 15 ultraluminous infrared (IR) galaxies (L_IR>10^12 Lsun; ULIRGs) using the Long Wavelength Spectrometer (LWS) on the Infrared Space Observatory (ISO). We confirm the observed order of magnitude deficit (compared to normal and starburst galaxies) in the strength of the [CII] line relative to the far-IR dust continuum emission found in our initial report (Luhman et al. 1998), but here with a sample that is twice as large. This result suggests that the deficit is a general phenomenon affecting 4/5 ULIRGs. We present an analysis using observations of generally acknowledged photodissociation region (PDR) tracers ([CII], [OI] 63 and 145 um, and FIR continuum emission), which suggests that a high UV flux G_o incident on a moderate density n PDR could explain the deficit. However, comparisons with other ULIRG observations, including CO (1-0), [CI] (1-0), and 6.2 um polycyclic aromatic hydrocarbon (PAH) emission, suggest that high G_o/n PDRs alone cannot produce a self-consistent solution that is compatible with all of the observations. We propose that non-PDR contributions to the FIR continuum can explain the apparent [CII] deficiency. Here, unusually high G_o and/or n physical conditions in ULIRGs as compared to those in normal and starburst galaxies are not required to explain the [CII] deficit. Dust-bounded photoionization regions, which generate much of the FIR emission but do not contribute significant [CII] emission, offer one possible physical origin for this additional non-PDR component. Such environments may also contribute to the observed suppression of FIR fine-structure emission from ionized gas and PAHs, as well as the warmer FIR colors found in ULIRGs. The implications for observations at higher redshifts are also revisited.

We present spectroscopy of MIR emission lines in twelve starburst regions, located in eleven starburst galaxies, for which a significant number of lines between 2.38 and 45micron were observed with the ISO Short Wavelength Spectrometer, with the intention of providing a reference resource for MIR spectra of starburst galaxies. The observation apertures were centred on actively star forming regions, including those which are inaccessible at optical wavelengths due to high levels of obscuration. We use this data set, which includes fine structure and hydrogen recombination lines, to investigate excitation and to derive gas phase abundances of neon, argon, and sulphur of the starburst galaxies. The derived Ne abundances span approximately an order of magnitude, up to values of ~3 times solar. The excitation ratios measured from the Ne and Ar lines correlate well with each other (positively) and with abundances (negatively). Both in excitation and abundance, a separation of objects with visible Wolf-Rayet features (high excitation, low abundance) is noted from those without (low excitation, high abundance). For a given abundance, the starbursts are of relatively lower excitation than a comparative sample of HII regions, possibly due to ageing stellar populations. By considering the abundance ratios of S with Ne and Ar we find that, in our higher metallicity systems, S is relatively underabundant by a factor of ~3. We discuss the origin of this deficit and favour depletion of S onto dust grains as a likely explanation. This weakness of the MIR fine structure lines of S has ramifications for future infrared missions such as SIRTF and Herschel since it indicates that the S lines are less favourable tracers of star formation than is suggested by nebular models which do not consider this effect.

The state of the matter that is obscuring a small circumnuclear region in AGN can be probed by observations of its broad emission lines. Infrared lines are particularly useful since they penetrate significant columns of obscuring matter, the properties of which can be constrained by comparing infrared and X-ray obscuration. We report on new 4um spectroscopy of a sample of 12 Seyfert 2 galaxies, probing for broad components to the Brackett alpha 4.05um line. Broad components are observed in 3 to 4 objects. All objects with a broad component exhibit relatively low X-ray obscuring columns, and our results are consistent with a Galactic ratio of 4um obscuration to the BLR and X-ray column. In combination with observations of a NON-Galactic ratio of VISUAL obscuration of BLRs and X-ray obscuring column in Seyferts, and interpreted in a unified AGN scheme, this result can be reconciled with two interpretations. Either the properties of dust near the AGN are modified towards larger grains, for example through coagulation, in a way that significantly flattens the optical/IR extinction curve, or the ratio of dust obscuration to X-ray column varies for different viewing angles with respect to the axis of symmetry of the putative torus. Our spectra also provide a survey of emission in the [SiIX] 3.94um coronal line, finding variation by an order of magnitude in its ratio to Br alpha. The first extragalactic detection of the [CaVII] 4.09um and [CaV] 4.16um coronal lines is reported in the spectrum of the Circinus galaxy.

In this paper we report on an XMM-Newton observation of the ultraluminous infrared QSO Mrk 1014. The X-ray observation reveals a power-law dominated (photon index of about 2.2) spectrum with a slight excess in the soft energy range. AGN and starburst emission models fit the soft excess emission equally well, however, the most plausible explanation is an AGN component as the starburst model parameter, temperature and luminosity, appear physically unrealistic. The mean luminosity of Mrk 1014 is about 2 times 10^44 erg s^-1. We have also observed excess emission at energies greater than 5 keV. This feature could be attributed to a broadened and redshifted iron complex, but deeper observations are required to constrain its origin. The light curve shows small scale variability over the 11 ks observation. There is no evidence of intrinsic absorption in Mrk 1014. The X-ray observations support the notion of an AGN dominated central engine. We establish the need for a longer observation to constrain more precisely the nature of the X-ray components.

We present medium resolution (R approx. 1500) ISO-SWS 2.4--45 micron spectra of a sample of 29 galaxies with active nuclei. This data set is rich in fine structure emission lines tracing the narrow line regions and (circum-)nuclear star formation regions, and it provides a coherent spectroscopic reference for future extragalactic studies in the mid-infrared. We use the data set to briefly discuss the physical conditions in the narrow line regions (density, temperature, excitation, line profiles) and to test for possible differences between AGN sub-types. Our main focus is on new tools for determining the propertibes of dusty galaxies and on the AGN-starburst connection. We present mid-IR line ratio diagrams which can be used to identify composite (starburst + AGN) sources and to distinguish between emission excited by active nuclei and emission from (circum-nuclear) star forming regions. For instance, line ratios of high to low excitation lines like [O IV]25.9um/[Ne II]12.8um, that have been used to probe for AGNs in dusty objects, can be examined in more detail and with better statistics now. In addition, we present two-dimensional diagnostic diagrams that are fully analogous to classical optical diagnostic diagrams, but better suited for objects with high extinction. Finally, we discuss correlations of mid-infrared line fluxes to the mid- and far-infrared continuum. We compare these relations to similar relations in starburst galaxies in order to examine the contribution of AGNs to the bolometric luminosities of their host galaxies. The spectra are available in electronic form from the authors.

We present low resolution mid infrared spectra of 16 ultraluminous infrared galaxies (ULIRGs) obtained with the CVF spectroscopy mode of ISOCAM on board the Infrared Space Observatory ISO. Our sample completes previous ISO spectroscopy of ultra- and hyperluminous infrared galaxies towards higher luminosities. The combined samples cover an infrared luminosity range of \~10^{12 -13.1} Lo. For about half of the high luminosity ULIRGs studied here, strong aromatic emission bands suggest starburst dominance. Other spectra are dominated by a strong AGN-related continuum with weak superposed emission features of uncertain nature. An improved method to quantitatively characterize the relative contribution of star formation and AGN activity to the mid-infrared emission of ULIRGs is presented. As dominant source of the bolometric luminosity, starbursts prevail at the lower end and AGNs at the higher end of this range. The transition between mostly starburst and mostly AGN powered occurs at ~10^{12.4}$ to $10^{12.5} Lo, and individual luminous starbursts are found up to ~10^{12.65} Lo.

We present ISOPHOT_S, ISOSWS and 8 to 13 micron ground based observations of Centaurus A that show prominent PAH and silicate features. These and other data are used to construct a model for the infrared continuum. We find that in a nuclear sized aperture (~4 arcsec, ~60 pc) the SED is characteristic of emission from a starburst and AGN torus; in larger apertures an additional component of cirrus emission is required. Based on our model, the torus diameter is estimated to be 3.6 pc and the best fitting inclination angle of the torus is 45 degrees. This result has implications for the detectability of tori in low power AGN and in particular for the use of the IRAS 60/25 micron flux ratio as an indicator of the torus inclination.

We present an inventory of mid-infrared spectral features detected in high resolution 2.4 to 45um ISO-SWS spectra of the galaxies M82, NGC253, Circinus, NGC1068, and the 30Doradus region of the LMC. We discuss their identifications and highlight possible relations between these features and the physical state of the interstellar medium in galaxies. In the 5--13 um region, emission from unidentified infrared bands (UIBs), usually ascribed to aromatic molecules, and apparent silicate absorption dominate the spectrum. We demonstrate the danger of overestimating the amount of silicate absorption by comparing the spectrum of M82 to the (absorption free) spectrum of the reflection nebula NGC7023. No strong silicate absorption is present in M82. The (very small grain) dust continuum under the UIB emission in our starburst templates can be modeled by a simple power law, starting at wavelengths between 8 and 9um. The spectra are valuable templates for future mid-infrared missions. We smooth our data to simulate low resolution spectra as obtained with ISOCAM-CVF, ISOPHOT-S, and in the future with the low resolution mode of SIRTF-IRS, and use our high spectral resolution information to highlight possible identification problems at low resolving power that are caused by coincidences of lines and features.

The ISO-SWS 2.5-45 um infrared spectroscopic observations of the nucleus of the Seyfert 2 galaxy NGC 1068 (see companion paper) are combined with a compilation of UV to IR narrow emission line data to determine the spectral energy distribution (SED) of the obscured extreme-UV continuum that photoionizes the narrow line emitting gas in the active galactic nucleus. We search a large grid of gas cloud models and SEDs for the combination that best reproduces the observed line fluxes and NLR geometry. Our best fit model reproduces the observed line fluxes to better than a factor of 2 on average and is in general agreement with the observed NLR geometry. It has two gas components that are consistent with a clumpy distribution of dense outflowing gas in the center and a more extended distribution of less dense and more clumpy gas farther out that has no net outflow. The best fit SED has a deep trough at ~4 Ryd, which is consistent with an intrinsic Big Blue Bump that is partially absorbed by ~6x10^19 cm^-2 of neutral hydrogen interior to the NLR.

We present ISO-SWS spectroscopy of NGC 1068 for the wavelength range 2.4 to 45um, detecting a total of 36 emission lines. Most of the observed transitions are fine structure and recombination lines originating in the narrow line region. We compare the line profiles of optical lines and reddening-insensitive infrared lines to constrain the dynamical structure and extinction properties of the NLR. The considerable differences found are most likely explained by two effects. (1) The spatial structure of the NLR is a combination of a highly ionized outflow cone and lower excitation extended emission. (2) Parts of the NLR, mainly in the receding part at velocities above systemic, are subject to extinction that is significantly suppressing optical emission. Line asymmetries and net blueshifts remain, however, even for infrared fine structure lines suffering very little obscuration. This may be either due to an intrinsic asymmetry of the NLR, or due to a very high column density obscuring component which is hiding part of the NLR even from infrared view. Mid-infrared emission of molecular hydrogen in NGC 1068 arises in a dense molecular medium at temperatures of a few hundred Kelvin that is most likely closely related to the warm and dense components seen in the near-infrared H2 transitions, and in millimeter wave tracers of molecular gas. Any emission of the putative pc-scale molecular torus is likely overwhelmed by this larger scale emission.

We report infrared spectroscopy of the prototypical Seyfert 2 galaxy NGC 1068, aiming at detection of broad components of hydrogen recombination lines that originate in the obscured broad-line region. Using the Short Wavelength Spectrometer on board the Infrared Space Observatory, we have observed for the first time the regions of Brackett beta 2.626um and Pfund alpha 7.460um, and present improved data for Brackett alpha 4.052um. No significant broad components are detected, implying an equivalent visual extinction to the broad-line region of at least 50 magnitudes and an obscuring column density of at least 10^23 cm^-2. While consistent with a highly obscured broad-line region, as required by the classical unified scenario, these limits are not yet significant enough to discriminate strongly between different torus models or to constrain properties of the gas causing the very large X-ray obscuration. We discuss the systematic limitations of infrared broad-line region searches and suggest that Brackett alpha may often be the most favorable transition for future searches.

ISO spectral observations of the supernova remnant RCW103 are presented. This object is the prototype of relatively young remnants with fast shocks interacting with dense interstellar medium. The spectrum is dominated by prominent lines of [NeII], [SiII], [FeII] and other low excitation species which provide, for the first time, a simple and reliable estimate of the gas abundances of refractory (Si, Fe, Ni) and non-refractory (Ne, P, S, Cl, Ar) species. Apart from nickel, all the derived abundances are close to solar, confirming that the shock has destroyed all dust grains. Like the optical nickel lines, [NiII] 6.64micron yields Ni abundances a factor ~10 solar which we propose results from a large underestimation of the computed Ni+ collision strengths. The observed intensities and velocity widths of ionic lines are compatible with emission from the post-shock region alone with only a very small (if any) contribution from the photoionized precursor. This result does not agree with shock models which predict that the precursor should emit powerful line emission, especially from highly ionized species. The possible consequence of this on the modelling of Seyfert spectra is briefly discussed.

We present ISO-SWS and ISOPHOT-S spectroscopy of the Seyfert galaxy NGC 4151. We detect a total of 17 fine structure emission lines emitted by a wide range of low- and high-excitation ions, two rotational lines of molecular hydrogen, and the Br beta HI line. We find that the mid-IR fine structure line profiles display blue asymmetries which are very similar to those observed in the optical lines produced in the narrow line region. Because the mid-infrared lines are much less sensitive to extinction than are the optical lines this similarity places strong constraints on scenarios which have been invoked to explain the optical line asymmetries. For example, we are able to rule out the simplest radial-motion-plus-dust scenarios for the production of the line asymmetries. Our preferred model is that of a central, geometrically thin but optically thick, obscuring screen of sub-arcsecond extension, enclosing a total hydrogen gas mass of >~5E6 Msun. This mass may be molecular. The weakness of `PAH' emission features in the low resolution spectrum is evidence that star formation plays a minor role in the circumnuclear region of NGC 4151. In a companion paper (Alexander et al., accepted by ApJ), we use the rich set of mid-infrared lines to determine the obscured photoionizing continuum produced by the active galactic nucleus.

The ISO-SWS infrared spectroscopic observations of the nucleus of Seyfert galaxy NGC4151, which are described in a companion paper, are used together with a compilation of UV to IR narrow emission line data to determine the spectral shape of the obscured extreme-UV continuum that photoionizes the narrow line emitting gas in the active galactic nucleus. We present a new method to determine the best fitting photoionizing continuum and emission line cloud model from a heterogeneous set of emission line data. For NGC4151, we find a best-fit model that reproduces the observed line fluxes to within a factor of 2 on average, and which is consistent with the observed geometry of the optical Narrow Line Region (NLR). Our model consists of a clumpy, optically thick (ionization bounded) gas distribution, with a hydrogen gas density of ~1000 cm^-3 and a volume filling factor of 6.5*10^-4. Our best fitting spectral energy distribution (SED) falls sharply beyond the Lyman limit and then rises sharply again towards 100 eV. In particular, it does not display a `Big Blue Bump' signature of a hot accretion disk. We find that this SED, which best reproduces the NLR line emission on the 100--500 pc scale, does not have enough UV photons to produce the observed BLR recombination emission from the inner 1 pc. This suggests that the BLR is photoionized by the intrinsic continuum source, which does have a strong UV component (perhaps a Big Blue Bump), but that this UV component is absorbed by material located between the NLR and BLR. Our analysis suggests that the absorber consists of \~5*10^19 cm^-2 of neutral hydrogen. Such an absorber was independently detected by UV absorption lines (Kriss et al. 1992, 1995).

We present an ISO SWS and ISOPHOT-S, mid-infrared spectroscopic survey of 15 ultra-luminous IRAS galaxies. We combine the survey results with a detailed case study, based on near-IR and mm imaging spectroscopy, of one of the sample galaxies (UGC 5101). We compare the near- and mid-IR characteristics of these ultra-luminous galaxies to ISO and literature data of thirty starburst and active galactic nuclei (AGN), template galaxies. We find that 1) 70-80% of the ultra-luminous IRAS galaxies in our sample are predominantly powered by recently formed massive stars. 20-30% are powered by a central AGN. These conclusions are based on a new infrared 'diagnostic diagram' involving the ratio of high to low excitation mid-IR emission lines on the one hand, and on the strength of the 7.7um PAH feature on the other hand. 2) at least half of the sources probably have simultaneously an active nucleus and starburst activity in a 1-2 kpc diameter circum-nuclear disk/ring. 3) the mid-infrared emitting regions are highly obscured. After correction for these extinctions, we estimate that the star forming regions in ULIRGs have ages between 10^7 and 10^8 years, similar to but somewhat larger than those found in lower luminosity starburst galaxies. 4) in the sample we have studied there is no obvious trend for the AGN component to dominate in the most compact, and thus most advanced mergers. Instead, at any given time during the merger evolution, the time dependent compression of the circum-nuclear interstellar gas, the accretion rate onto the central black hole and the associated radiation efficiency may determine whether star formation or AGN activity dominates the luminosity of the system.