Submillimeter rotational lines of H2O are a powerful probe in warm gas regions of the ISM, tracing scales and structures ranging from kpc disks to the most compact and dust-obscured regions of galactic nuclei. The ortho-H2O(423-330) line at 448 GHz, which was recently detected in a local luminous infrared galaxy (Pereira-Santaella et al. 2017), offers a unique constraint on the excitation conditions and ISM properties in deeply buried galaxy nuclei since the line requires high far-IR optical depths to be excited. In this letter, we report the first high-redshift detection of the 448 GHz H2O(423-330) line using ALMA, in a strongly lensed submillimeter galaxy (SMG) at z=3.63. After correcting for magnification, the luminosity of the 448 GHz H2O line is ~10^6 L_sun. In combination with three other previously detected H2O lines, we build a model that "resolves" the dusty ISM structure of the SMG, and find that it is composed of a ~1 kpc optically thin (optical depth at 100{\mu}m {\tau}_{100}~0.3) disk component with dust temperature T_{dust} \approx 50 K emitting a total infrared power of 5e12 L_sun with surface density \Sigma_{IR}=4e11 L_sun kpc^{-2}, and a very compact (0.1 kpc) heavily dust-obscured ({\tau}_{100} \gtrsim 1) nuclear core with very warm dust (100 K) and \Sigma_{IR}=8e12 L_sun kpc^{-2}. The H2O abundance in the core component, X_{H2O}~(0.3-5)e{-5}, is at least one order of magnitude higher than in the disk component. The optically thick core has the characteristic properties of an Eddington-limited starburst, providing evidence that radiation pressure on dust is capable of supporting the ISM in buried nuclei at high redshifts. The multi-component ISM structure revealed by our models illustrates that dust and molecules such as H2O are present in regions characterized by highly differing conditions and scales, extending from the nucleus to more extended regions of SMGs.

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.

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)

Deep Learning (DL) algorithms for morphological classification of galaxies have proven very successful, mimicking (or even improving) visual classifications. However, these algorithms rely on large training samples of labelled galaxies (typically thousands of them). A key question for using DL classifications in future Big Data surveys is how much of the knowledge acquired from an existing survey can be exported to a new dataset, i.e. if the features learned by the machines are meaningful for different data. We test the performance of DL models, trained with Sloan Digital Sky Survey (SDSS) data, on Dark Energy survey (DES) using images for a sample of $\sim$5000 galaxies with a similar redshift distribution to SDSS. Applying the models directly to DES data provides a reasonable global accuracy ($\sim$ 90%), but small completeness and purity values. A fast domain adaptation step, consisting in a further training with a small DES sample of galaxies ($\sim$500-300), is enough for obtaining an accuracy > 95% and a significant improvement in the completeness and purity values. This demonstrates that, once trained with a particular dataset, machines can quickly adapt to new instrument characteristics (e.g., PSF, seeing, depth), reducing by almost one order of magnitude the necessary training sample for morphological classification. Redshift evolution effects or significant depth differences are not taken into account in this study.

We describe the SDSS-IV MaNGA PyMorph Photometric (MPP-VAC) and MaNGA Deep Learning Morphology (MDLM-VAC) Value Added Catalogs. The MPP-VAC provides photometric parameters from S\'ersic and S\'ersic+Exponential fits to the 2D surface brightness profiles of the MaNGA DR15 galaxy sample. Compared to previous PyMorph analyses of SDSS imaging, our analysis of the MaNGA DR15 incorporates three improvements: the most recent SDSS images; modified criteria for determining bulge-to-disk decompositions; and the fits in MPP-VAC have been eye-balled, and re-fit if necessary, for additional reliability. A companion catalog, the MDLM-VAC, provides Deep Learning-based morphological classifications for the same galaxies. The MDLM-VAC includes a number of morphological properties (e.g., a TType, and a finer separation between elliptical and S0 galaxies). Combining the MPP- and MDLM-VACs allows to show that the MDLM morphological classifications are more reliable than previous work. It also shows that single-S\'ersic fits to late- and early-type galaxies are likely to return S\'ersic indices of $n \le 2$ and $\ge 4$, respectively, and this correlation between $n$ and morphology extends to the bulge component as well. While the former is well-known, the latter contradicts some recent work suggesting little correlation between $n$-bulge and morphology. Combining both VACs with MaNGA's spatially resolved spectroscopy allows us to study how the stellar angular momentum depends on morphological type. We find correlations between stellar kinematics, photometric properties, and morphological type even though the spectroscopic data played no role in the construction of the MPP- and MDLM-VACs.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

Aiming to characterise the properties of the molecular gas in the ultraluminous infrared galaxy Mrk273 and its outflow, we used the NOEMA interferometer to image the dense gas molecular tracers HCN, HCO+, HNC, HOC+ and HC3N at 86GHz and 256GHz with angular resolutions of 4.9x4.5 arcsec (3.7x3.4 kpc) and 0.61x0.55 arcsec (460x420 pc). We also modelled the flux of several H2O lines observed with Herschel using a radiative transfer code that includes excitation by collisions as well as by far-infrared photons. The disk of the Mrk273 north nucleus has two components with decoupled kinematics. The gas in the outer parts (1.5 kpc) rotates with a south-east to north-west direction, while in the inner disk (300 pc) follows a north-east to south-west rotation. The central 300 pc, which hosts a compact starburst region, is filled with dense and warm gas, contains a dynamical mass of (4-5)x10^9M_sun, a luminosity of L'_HCN=(3-4)x10^8 K km/s pc^2, and a dust temperature of 55 K. At the very centre, a compact core with R~50 pc has a luminosity of L_IR=4x10^11L_sun (30% of the total infrared luminosity), and a dust temperature of 95 K. The core is expanding at low velocities ~50-100 km/s, probably affected by the outflowing gas. We detect the blue-shifted component of the outflow, while the red-shifted counterpart remains undetected in our data. Its cold and dense phase reaches fast velocities up to ~1000 km/s, while the warm outflowing gas has more moderate maximum velocities of ~600 km/s. The outflow is detected as far as 460 pc from the centre in the northern direction, and has a mass of dense gas <8x10^8M_sun. The difference between the position angles of the inner disk (~70 degree) and the outflow (~10 degree) indicates that the outflow is likely powered by the AGN, and not by the starburst. Regarding the chemistry, we measure an extremely low HCO+/HOC+ ratio of 10+-5 in the inner disk of Mrk273.

We report molecular gas observations of IRAS 20100-4156 and IRAS 03158+4227, two local ultraluminous infrared galaxies (ULIRGs) hosting some of the fastest and most massive molecular outflows known. Using ALMA and PdBI observations, we spatially resolve the CO(1-0) emission from the outflowing molecular gas in both and find maximum outflow velocities of $ v_{\rm max} \sim 1600$ and $\sim 1700$ km/s for IRAS 20100-4156 and IRAS 03158+4227, respectively. We find total gas mass outflow rates of $\dot M_{\rm OF} \sim 670$ and $\sim 350$ Msun/yr, respectively, corresponding to molecular gas depletion timescales $\tau^{\rm dep}_{\rm OF} \sim 11$ and $\sim 16$ Myr. This is nearly 3 times shorter than the depletion timescales implied by star formation, $\tau^{\rm dep}_{\rm SFR} \sim 33$ and $\sim 46$ Myr, respectively. To determine the outflow driving mechanism, we compare the starburst ($L_{*}$) and AGN ($L_{\rm AGN}$) luminosities to the outflowing energy and momentum fluxes, using mid-infrared spectral decomposition to discern $L_{\rm AGN}$. Comparison to other molecular outflows in ULIRGs reveals that outflow properties correlate similarly with $L_{*}$ and $L_{\rm IR}$ as with $L_{\rm AGN}$, indicating that AGN luminosity alone may not be a good tracer of feedback strength and that a combination of AGN and starburst activity may be driving the most powerful molecular outflows. We also detect the OH 1.667 GHz maser line from both sources and demonstrate its utility in detecting molecular outflows.

Within a galaxy the stellar mass-to-light ratio $\Upsilon_*$ is not constant. Spatially resolved kinematics of nearby early-type galaxies suggest that allowing for a variable initial mass function (IMF) returns significantly larger $\Upsilon_*$ gradients than if the IMF is held fixed. If $\Upsilon_*$ is greater in the central regions, then ignoring the IMF-driven gradient can overestimate $M_*^{\rm dyn}$ by as much as a factor of two for the most massive galaxies, though stellar population estimates $M_*^{\rm SP}$ are also affected. Large $\Upsilon_*$-gradients have four main consequences: First, $M_*^{\rm dyn}$ cannot be estimated independently of stellar population synthesis models. Second, if there is a lower limit to $\Upsilon_*$ and gradients are unknown, then requiring $M_*^{\rm dyn}=M_*^{\rm SP}$ constrains them. Third, if gradients are stronger in more massive galaxies, then $M_*^{\rm dyn}$ and $M_*^{\rm SP}$ can be brought into agreement, not by shifting $M_*^{\rm SP}$ upwards by invoking constant bottom-heavy IMFs, as advocated by a number of recent studies, but by revising $M_*^{\rm dyn}$ estimates in the literature downwards. Fourth, accounting for $\Upsilon_*$ gradients changes the high-mass slope of the stellar mass function $\phi(M_*^{\rm dyn})$, and reduces the associated stellar mass density. These conclusions potentially impact estimates of the need for feedback and adiabatic contraction, so our results highlight the importance of measuring $\Upsilon_*$ gradients in larger samples.

Inflation is the leading theory of the first instant of the universe. Inflation, which postulates that the universe underwent a period of rapid expansion an instant after its birth, provides convincing explanation for cosmological observations. Recent advancements in detector technology have opened opportunities to explore primordial gravitational waves generated by the inflation through B-mode (divergent-free) polarization pattern embedded in the Cosmic Microwave Background anisotropies. If detected, these signals would provide strong evidence for inflation, point to the correct model for inflation, and open a window to physics at ultra-high energies. LiteBIRD is a satellite mission with a goal of detecting degree-and-larger-angular-scale B-mode polarization. LiteBIRD will observe at the second Lagrange point with a 400 mm diameter telescope and 2,622 detectors. It will survey the entire sky with 15 frequency bands from 40 to 400 GHz to measure and subtract foregrounds. The U.S. LiteBIRD team is proposing to deliver sub-Kelvin instruments that include detectors and readout electronics. A lenslet-coupled sinuous antenna array will cover low-frequency bands (40 GHz to 235 GHz) with four frequency arrangements of trichroic pixels. An orthomode-transducer-coupled corrugated horn array will cover high-frequency bands (280 GHz to 402 GHz) with three types of single frequency detectors. The detectors will be made with Transition Edge Sensor (TES) bolometers cooled to a 100 milli-Kelvin base temperature by an adiabatic demagnetization refrigerator.The TES bolometers will be read out using digital frequency multiplexing with Superconducting QUantum Interference Device (SQUID) amplifiers. Up to 78 bolometers will be multiplexed with a single SQUID amplidier. We report on the sub-Kelvin instrument design and ongoing developments for the LiteBIRD mission.

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.

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

We present a morphological catalogue for $\sim$ 670,000 galaxies in the Sloan Digital Sky Survey in two flavours: T-Type, related to the Hubble sequence, and Galaxy Zoo 2 (GZ2 hereafter) classification scheme. By combining accurate existing visual classification catalogues with machine learning, we provide the largest and most accurate morphological catalogue up to date. The classifications are obtained with Deep Learning algorithms using Convolutional Neural Networks (CNNs). We use two visual classification catalogues, GZ2 and Nair & Abraham (2010), for training CNNs with colour images in order to obtain T-Types and a series of GZ2 type questions (disk/features, edge-on galaxies, bar signature, bulge prominence, roundness and mergers). We also provide an additional probability enabling a separation between pure elliptical (E) from S0, where the T-Type model is not so efficient. For the T-Type, our results show smaller offset and scatter than previous models trained with support vector machines. For the GZ2 type questions, our models have large accuracy (> 97\%), precision and recall values (> 90\%) when applied to a test sample with the same characteristics as the one used for training. The catalogue is publicly released with the paper.

Spatially resolved kinematics of nearby galaxies has shown that the ratio of dynamical- to stellar population-based estimates of the mass of a galaxy ($M_*^{\rm JAM}/M_*$) correlates with $\sigma_e$, if $M_*$ is estimated using the same IMF for all galaxies and the stellar M/L ratio within each galaxy is constant. This correlation may indicate that, in fact, the IMF is more dwarf-rich for galaxies with large $\sigma$. We use this correlation to estimate a dynamical or IMF-corrected stellar mass, $M_*^{\rm \alpha_{JAM}}$, from $M_{*}$ and $\sigma_e$ for a sample of $6 \times 10^5$ SDSS galaxies for which spatially resolved kinematics is not available. We also compute the `virial' mass estimate $k(n,R)\,R_e\,\sigma_R^2/G$, where $n$ is the Sersic index, in the SDSS and ATLAS$^{\rm 3D}$ samples. We show that an $n$-dependent correction must be applied to the $k(n,R)$ values provided by Prugniel & Simien (1997). Our analysis also shows that the shape of the velocity dispersion profile in the ATLAS$^{\rm 3D}$ sample varies weakly with $n$: $(\sigma_R/\sigma_e) = (R/R_e)^{-\gamma(n)}$. The resulting stellar mass functions, based on $M_*^{\rm \alpha_{JAM}}$ and the recalibrated virial mass, are in good agreement. If the $M_*^{\rm \alpha_{JAM}}/M_* - \sigma_e$ correlation is indeed due to the IMF, and stellar M/L gradients can be ignored, then our $\phi(M_*^{\rm \alpha_{JAM}})$ is an estimate of the stellar mass function in which $\sigma_e$-dependent variations in the IMF across the population have been accounted for. Using a Fundamental Plane based observational proxy for $\sigma_e$ produces comparable results. By demonstrating that cheaper proxies are sufficiently accurate, our analysis should enable a more reliable census of the mass in stars for large galaxy samples, at a fraction of the cost. Our results are provided in tabular form.

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.

The physical processes driving the chemical evolution of galaxies in the last $\sim 11\, \rm{Gyr}$ cannot be understood without directly probing the dust-obscured phase of star-forming galaxies and active galactic nuclei. This phase, hidden to optical tracers, represents the bulk of star formation and black hole accretion activity in galaxies at $1 < z < 3$. Spectroscopic observations with a cryogenic infrared (IR) observatory like SPICA will be sensitive enough to peer through the dust-obscured regions of galaxies and access the rest-frame mid- to far-IR range in galaxies at high-$z$. This wavelength range contains a unique suite of spectral lines and dust features that serve as proxies for the abundances of heavy elements and the dust composition, providing tracers with a feeble response to both extinction and temperature. In this work, we investigate how SPICA observations could be exploited to understand key aspects in the chemical evolution of galaxies: the assembly of nearby galaxies based on the spatial distribution of heavy element abundances, the global content of metals in galaxies reaching the knee of the luminosity function up to $z \sim 3$, and the dust composition of galaxies at high-$z$. Possible synergies with facilities available in the late 2020s are also discussed.

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.

The Sloan Digital Sky Survey pipeline photometry underestimates the brightnesses of the most luminous galaxies. This is mainly because (i) the SDSS overestimates the sky background and (ii) single or two-component Sersic-based models better fit the surface brightness profile of galaxies, especially at high luminosities, than does the de Vaucouleurs model used by the SDSS pipeline. We use the PyMorph photometric reductions to isolate effect (ii) and show that it is the same in the full sample as in small group environments, and for satellites in the most massive clusters as well. None of these are expected to be significantly affected by intracluster light (ICL). We only see an additional effect for centrals in the most massive halos, but we argue that even this is not dominated by ICL. Hence, for the vast majority of galaxies, the differences between PyMorph and SDSS pipeline photometry cannot be ascribed to the semantics of whether or not one includes the ICL when describing the stellar mass of massive galaxies. Rather, they likely reflect differences in star formation or assembly histories. Failure to account for the SDSS underestimate has significantly biased most previous estimates of the SDSS luminosity and stellar mass functions, and therefore Halo Model estimates of the z ~ 0.1 relation between the mass of a halo and that of the galaxy at its center. We also show that when one studies correlations, at fixed group mass, with a quantity which was not used to define the groups, then selection effects appear. We show why such effects arise, and should not be mistaken for physical effects.

We quantify the systematic effects on the stellar mass function which arise from assumptions about the stellar population, as well as how one fits the light profiles of the most luminous galaxies at z ~ 0.1. When comparing results from the literature, we are careful to separate out these effects. Our analysis shows that while systematics in the estimated comoving number density which arise from different treatments of the stellar population remain of order < 0.5 dex, systematics in photometry are now about 0.1 dex, despite recent claims in the literature. Compared to these more recent analyses, previous work based on Sloan Digital Sky Survey (SDSS) pipeline photometry leads to underestimates of rho_*(> M_*) by factors of 3-10 in the mass range 10^11 - 10^11.6 M_Sun, but up to a factor of 100 at higher stellar masses. This impacts studies which match massive galaxies to dark matter halos. Although systematics which arise from different treatments of the stellar population remain of order < 0.5 dex, our finding that systematics in photometry now amount to only about 0.1 dex in the stellar mass density is a significant improvement with respect to a decade ago. Our results highlight the importance of using the same stellar population and photometric models whenever low and high redshift samples are compared.

A number of recent estimates of the total luminosities of galaxies in the SDSS are significantly larger than those reported by the SDSS pipeline. This is because of a combination of three effects: one is simply a matter of defining the scale out to which one integrates the fit when defining the total luminosity, and amounts on average to < 0.1 mags even for the most luminous galaxies. The other two are less trivial and tend to be larger; they are due to differences in how the background sky is estimated and what model is fit to the surface brightness profile. We show that PyMorph sky estimates are fainter than those of the SDSS DR7 or DR9 pipelines, but are in excellent agreement with the estimates of Blanton et al. (2011). Using the SDSS sky biases luminosities by more than a few tenths of a magnitude for objects with half-light radii > 7 arcseconds. In the SDSS main galaxy sample these are typically luminous galaxies, so they are not necessarily nearby. This bias becomes worse when allowing the model more freedom to fit the surface brightness profile. When PyMorph sky values are used, then two component Sersic-Exponential fits to E+S0s return more light than single component deVaucouleurs fits (up to ~0.2 mag), but less light than single Sersic fits (0.1 mag). Finally, we show that PyMorph fits of Meert et al. (2015) to DR7 data remain valid for DR9 images. Our findings show that, especially at large luminosities, these PyMorph estimates should be preferred to the SDSS pipeline values.

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 a combined morphological and X-ray analysis of Was 49, an isolated, dual AGN system notable for the presence of a dominant AGN Was 49b in the disk of the primary galaxy Was 49a, at a projected radial distance of 8 kpc from the nucleus. Using X-ray data from Chandra, NuSTAR, and Swift, we find that this AGN has a bolometric luminosity of L_bol ~ 2 x 10^45 erg/s, with a black hole mass of M_BH=1.3^{+2.9}_{-0.9} x 10^8 M_Sol. Despite its large mass, our analysis of optical data from the Discovery Channel Telescope shows that the supermassive black hole is hosted by a stellar counterpart with a mass of only 5.6^{+4.9}_{-2.6} x 10^9 M_Sol, making the SMBH potentially larger than expected from SMBH-galaxy scaling relations, and the stellar counterpart exhibits a morphology that is consistent with dwarf elliptical galaxies. Our analysis of the system in the r and K bands indicates that Was 49 is a minor merger, with a mass ratio of Was 49a to Was 49b between 1:7 and 1:15. This is in contrast with findings that the most luminous merger-triggered AGNs are found in major mergers, and that minor mergers predominantly enhance AGN activity in the primary galaxy.

Aims. We probe the physical conditions in the core of Arp 299A and try to put constraints to the nature of its nuclear power source. Methods. We used Herschel Space Observatory far-infrared and submillimeter observations of H2O and OH rotational lines in Arp 299A to create a multi-component model of the galaxy. In doing this, we employed a spherically symmetric radiative transfer code. Results. Nine H2O lines in absorption and eight in emission as well as four OH doublets in absorption and one in emission, are detected in Arp 299A. No lines of the 18O isotopologues, which have been seen in compact obscured nuclei of other galaxies, are detected. The absorption in the ground state OH doublet at 119 {\mu}m is found redshifted by ~175 km/s compared to other OH and H2O lines, suggesting a low excitation inflow. We find that at least two components are required in order to account for the excited molecular line spectrum. The inner component has a radius of 20-25 pc, a very high infrared surface brightness (> 3e13 Lsun/kpc^2), warm dust (Td > 90 K), and a large H2 column density (NH2 > 1e24 cm^-2). The outer component is larger (50-100 pc) with slightly cooler dust (70-90 K). In addition, a much more extended inflowing component is required to also account for the OH doublet at 119 {\mu}m. Conclusions. The Compton-thick nature of the core makes it difficult to determine the nature of the buried power source, but the high surface brightness indicates that it is either an active galactic nucleus and/or a dense nuclear starburst. The high OH/H2O ratio in the nucleus indicates that ion-neutral chemistry induced by X-rays or cosmic-rays is important. Finally we find a lower limit to the 16O/18O ratio of 400 in the nuclear region, possibly indicating that the nuclear starburst is in an early evolutionary stage, or that it is fed through a molecular inflow of, at most, solar metallicity.

Host galaxy identification is a crucial step for modern supernova (SN) surveys such as the Dark Energy Survey (DES) and the Large Synoptic Survey Telescope (LSST), which will discover SNe by the thousands. Spectroscopic resources are limited, so in the absence of real-time SN spectra these surveys must rely on host galaxy spectra to obtain accurate redshifts for the Hubble diagram and to improve photometric classification of SNe. In addition, SN luminosities are known to correlate with host-galaxy properties. Therefore, reliable identification of host galaxies is essential for cosmology and SN science. We simulate SN events and their locations within their host galaxies to develop and test methods for matching SNe to their hosts. We use both real and simulated galaxy catalog data from the Advanced Camera for Surveys General Catalog and MICECATv2.0, respectively. We also incorporate "hostless" SNe residing in undetected faint hosts into our analysis, with an assumed hostless rate of 5%. Our fully automated algorithm is run on catalog data and matches SNe to their hosts with 91% accuracy. We find that including a machine learning component, run after the initial matching algorithm, improves the accuracy (purity) of the matching to 97% with a 2% cost in efficiency (true positive rate). Although the exact results are dependent on the details of the survey and the galaxy catalogs used, the method of identifying host galaxies we outline here can be applied to any transient survey.

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.

The collapse of the core of a star is expected to produce gravitational radiation. While this process will usually produce a luminous supernova, the optical signatue could be subluminous and a direct collapse to a black hole, with the star just disappearing, is possible. The gravitational wave event GW150914 reported by the LIGO Virgo Collaboration (LVC) on 2015 September 16, was detected by a burst analysis and whose high probability spatial localization included the Large Magellanic Cloud. Shortly after the announcement of the event, we used the Dark Energy Camera to observe 102 deg$^2$ of the localization area, including a 38 deg$^2$ area centered on the LMC. Using a catalog of 152 LMC luminous red supergiants, candidates to undergo a core collapse without a visible supernova, we find that the positions of 144 of these are inside our images, and that all are detected - none have disappeared. There are other classes of candidates: we searched existing catalogs of red supergiants, yellow supergiants, Wolf-Rayet stars, and luminous blue variable stars, recovering all that were inside the imaging area. Based on our observations, we conclude that it is unlikely that GW150914 was caused by the core collapse of a supergiant in the LMC, consistent with the LIGO Collaboration analyses of the gravitational wave form as best described by a binary black hole merger. We discuss how to generalize this search for future very nearby core collapse candidates.

We report initial results of a deep search for an optical counterpart to the gravitational wave event GW150914, the first trigger from the Advanced LIGO gravitational wave detectors. We used the Dark Energy Camera (DECam) to image a 102 deg$^2$ area, corresponding to 38% of the initial trigger high-probability sky region and to 11% of the revised high-probability region. We observed in i and z bands at 4-5, 7, and 24 days after the trigger. The median $5\sigma$ point-source limiting magnitudes of our search images are i=22.5 and z=21.8 mag. We processed the images through a difference-imaging pipeline using templates from pre-existing Dark Energy Survey data and publicly available DECam data. Due to missing template observations and other losses, our effective search area subtends 40 deg$^{2}$, corresponding to 12% total probability in the initial map and 3% of the final map. In this area, we search for objects that decline significantly between days 4-5 and day 7, and are undetectable by day 24, finding none to typical magnitude limits of i= 21.5,21.1,20.1 for object colors (i-z)=1,0,-1, respectively. Our search demonstrates the feasibility of a dedicated search program with DECam and bodes well for future research in this emerging field.

Using the largest single-survey sample of Type Ia supernovae (SNe Ia) to date, we study the relationship between properties of SNe Ia and those of their host galaxies, focusing primarily on correlations with Hubble residuals (HR). Our sample consists of 345 photometrically-classified or spectroscopically-confirmed SNeIa discovered as part of the SDSS-II Supernova Survey (SDSS-SNS). This analysis utilizes host-galaxy spectroscopy obtained during the SDSS-I/II spectroscopic survey and from an ancillary program on the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) that obtained spectra for nearly all host galaxies of SDSS-II SN candidates. In addition, we use photometric host-galaxy properties from the SDSS-SNS data release (Sako et al. 2014) such as host stellar mass and star-formation rate. We confirm the well-known relation between HR and host-galaxy mass and find a 3.6{\sigma} significance of a non-zero linear slope. We also recover correlations between HR and host-galaxy gas-phase metallicity and specific star-formation rate as they are reported in the literature. With our large dataset, we examine correlations between HR and multiple host-galaxy properties simultaneously and find no evidence of a significant correlation. We also independently analyze our spectroscopically-confirmed and photometrically-classified SNe Ia and comment on the significance of similar combined datasets for future surveys.

Using the largest single-survey sample of Type Ia supernovae (SNe Ia) to date, we study the relationship between properties of SNe Ia and those of their host galaxies, focusing primarily on correlations with Hubble residuals (HR). Our sample consists of 345 photometrically-classified or spectroscopically-confirmed SNeIa discovered as part of the SDSS-II Supernova Survey (SDSS-SNS). This analysis utilizes host-galaxy spectroscopy obtained during the SDSS-I/II spectroscopic survey and from an ancillary program on the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) that obtained spectra for nearly all host galaxies of SDSS-II SN candidates. In addition, we use photometric host-galaxy properties from the SDSS-SNS data release (Sako et al. 2014) such as host stellar mass and star-formation rate. We confirm the well-known relation between HR and host-galaxy mass and find a 3.6{\sigma} significance of a non-zero linear slope. We also recover correlations between HR and host-galaxy gas-phase metallicity and specific star-formation rate as they are reported in the literature. With our large dataset, we examine correlations between HR and multiple host-galaxy properties simultaneously and find no evidence of a significant correlation. We also independently analyze our spectroscopically-confirmed and photometrically-classified SNe Ia and comment on the significance of similar combined datasets for future surveys.

We describe an algorithm for identifying point-source transients and moving objects on reference-subtracted optical images containing artifacts of processing and instrumentation. The algorithm makes use of the supervised machine learning technique known as Random Forest. We present results from its use in the Dark Energy Survey Supernova program (DES-SN), where it was trained using a sample of 898,963 signal and background events generated by the transient detection pipeline. After reprocessing the data collected during the first DES-SN observing season (Sep. 2013 through Feb. 2014) using the algorithm, the number of transient candidates eligible for human scanning decreased by a factor of 13.4, while only 1 percent of the artificial Type Ia supernovae (SNe) injected into search images to monitor survey efficiency were lost, most of which were very faint events. Here we characterize the algorithm's performance in detail, and we discuss how it can inform pipeline design decisions for future time-domain imaging surveys, such as the Large Synoptic Survey Telescope and the Zwicky Transient Facility. An implementation of the algorithm and the training data used in this paper are available at http://portal.nersc.gov/project/dessn/autoscan.

We describe the difference imaging pipeline (DiffImg) used to detect transients in deep images from the Dark Energy Survey Supernova program (DES-SN) in its first observing season from Aug 2013 through Feb 2014. DES-SN is a search for transients in which ten 3-deg^2 fields are repeatedly observed in the g,r,i,z passbands with a cadence of about 1 week. The observing strategy has been optimized to measure high-quality light curves and redshifts for thousands of Type Ia supernova (SN Ia) with the goal of measuring dark energy parameters. The essential DiffImg functions are to align each search image to a deep reference image, do a pixel-by-pixel subtraction, and then examine the subtracted image for significant positive detections of point-source objects. The vast majority of detections are subtraction artifacts, but after selection requirements and image filtering with an automated scanning program, there are 130 detections per deg^2 per observation in each band, of which only 25% are artifacts. Of the 7500 transients discovered by DES-SN in its first observing season, each requiring a detection on at least 2 separate nights, Monte Carlo simulations predict that 27% are expected to be supernova. Another 30% of the transients are artifacts, and most of the remaining transients are AGN and variable stars. Fake SNe Ia are overlaid onto the images to rigorously evaluate detection efficiencies, and to understand the DiffImg performance. The DiffImg efficiency measured with fake SNe agrees well with expectations from a Monte Carlo simulation that uses analytical calculations of the fluxes and their uncertainties. In our 8 "shallow" fields with single-epoch 50% completeness depth 23.5, the SN Ia efficiency falls to 1/2 at redshift z 0.7, in our 2 "deep" fields with mag-depth 24.5, the efficiency falls to 1/2 at z 1.1.

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.

Context. The LIRG Zw 049.057 contains a compact obscured nucleus where a considerable amount of the galaxy's luminosity is generated. This nucleus contains a dusty environment that is rich in molecular gas. One approach to probing this kind of environment and to revealing what is hidden behind the dust is to study the rotational lines of molecules that couple well with the IR radiation emitted by the dust. Methods. We observed Zw 049.057 with PACS and SPIRE onboard the Herschel Space Observatory in rotational lines of H2O, H218O, OH, 18OH, and [O I]. We modeled the unresolved core of the galaxy using a spherically symmetric radiative transfer code. Results. We present the full SPIRE FTS spectrum of Zw 049.057, along with relevant spectral scans in the PACS range. We find that a minimum of two different components (nuclear and extended) are required in order to account for the rich molecular line spectrum. The nuclear component has a radius of 10-30 pc, a very high infrared surface brightness (1e14 Lsun/kpc2), warm dust (Td > 100 K), and a very large H2 column density (NH2 = 1e24-1e25 cm-2). The modeling also indicates high nuclear H2O (5e-6) and OH (4e-6) abundances relative to H2 as well as a low 16O/18O-ratio of 50-100. We also find a prominent infall signature in the [O I] line. We tentatively detect a 500 km/s outflow in the H2O 313->202 line. Conclusions. The high surface brightness of the core indicates the presence of either a buried active galactic nucleus or a very dense nuclear starburst.nThe H2O abundance is comparable to that of other compact (ultra-)luminous infrared galaxies such as NGC 4418 and Arp 220 - and also to hot cores in the Milky Way. The enhancement of 18O is a possible indicator that the nucleus of Zw 049.057 is in a similar evolutionary stage as the nuclei of Arp 220 - and more advanced than NGC 4418.

We present DES13S2cmm, the first spectroscopically-confirmed superluminous supernova (SLSN) from the Dark Energy Survey (DES). We briefly discuss the data and search algorithm used to find this event in the first year of DES operations, and outline the spectroscopic data obtained from the European Southern Observatory (ESO) Very Large Telescope to confirm its redshift (z = 0.663 +/- 0.001 based on the host-galaxy emission lines) and likely spectral type (type I). Using this redshift, we find M_U_peak = -21.05 +0.10 -0.09 for the peak, rest-frame U-band absolute magnitude, and find DES13S2cmm to be located in a faint, low metallicity (sub-solar), low stellar-mass host galaxy (log(M/M_sun) = 9.3 +/- 0.3); consistent with what is seen for other SLSNe-I. We compare the bolometric light curve of DES13S2cmm to fourteen similarly well-observed SLSNe-I in the literature and find it possesses one of the slowest declining tails (beyond +30 days rest frame past peak), and is the faintest at peak. Moreover, we find the bolometric light curves of all SLSNe-I studied herein possess a dispersion of only 0.2-0.3 magnitudes between +25 and +30 days after peak (rest frame) depending on redshift range studied; this could be important for 'standardising' such supernovae, as is done with the more common type Ia. We fit the bolometric light curve of DES13S2cmm with two competing models for SLSNe-I - the radioactive decay of 56Ni, and a magnetar - and find that while the magnetar is formally a better fit, neither model provides a compelling match to the data. Although we are unable to conclusively differentiate between these two physical models for this particular SLSN-I, further DES observations of more SLSNe-I should break this degeneracy, especially if the light curves of SLSNe-I can be observed beyond 100 days in the rest frame of the supernova.

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

(Ultra) Luminous Infrared Galaxies ((U)LIRGs) are objects characterized by their extreme infrared (8-1000 $\mu$m) luminosities ($L_{LIRG}>10^{11} $L$_\odot$ and $L_{ULIRG}>10^{12}$ L$_\odot$). The Herschel Comprehensive ULIRG Emission Survey (HerCULES; PI van der Werf) presents a representative flux-limited sample of 29 (U)LIRGs that spans the full luminosity range of these objects (10$^{11}\leq L_\odot \geq10^{13}$). With the \emph{Herschel Space Observatory}, we observe [CII] 157 $\mu$m, [OI] 63 $\mu$m, and [OI] 145 $\mu$m line emission with PACS, CO J=4-3 through J=13-12, [CI] 370 $\mu$m, and [CI] 609 $\mu$m with SPIRE, and low-J CO transitions with ground-based telescopes. The CO ladders of the sample are separated into three classes based on their excitation level. In 13 of the galaxies, the [OI] 63 $\mu$m emission line is self absorbed. Comparing the CO excitation to the IRAS 60/100 $\mu$m ratio and to far infrared luminosity, we find that the CO excitation is more correlated to the far infrared colors. We present cooling budgets for the galaxies and find fine-structure line flux deficits in the [CII], [SiII], [OI], and [CI] lines in the objects with the highest far IR fluxes, but do not observe this for CO $4\leq J_{upp}\leq13$. In order to study the heating of the molecular gas, we present a combination of three diagnostic quantities to help determine the dominant heating source. Using the CO excitation, the CO J=1-0 linewidth, and the AGN contribution, we conclude that galaxies with large CO linewidths always have high-excitation CO ladders, and often low AGN contributions, suggesting that mechanical heating is important.

Supermassive black holes (SMBHs) are found ubiquitously in large, bulge-dominated galaxies throughout the local universe, yet little is known about their presence and properties in bulgeless and low-mass galaxies. This is a significant deficiency, since the mass distribution and occupation fraction of nonstellar black holes provide important observational constraints on SMBH seed formation theories and many dwarf galaxies have not undergone major mergers that would erase information on their original black hole population. Using data from the Wide-field Infrared Survey Explorer, we discovered hundreds of bulgeless and dwarf galaxies that display mid-infrared signatures of extremely hot dust highly suggestive of powerful accreting massive black holes, despite having no signatures of black hole activity at optical wavelengths. Here we report, in our first follow-up X-ray investigation of this population, that the irregular dwarf galaxy J132932.41+323417.0 (z = 0.0156) contains a hard, unresolved X-ray source detected by XMM-Newton with luminosity L_2-10 keV = 2.4 x 10^40 erg/s, over two orders of magnitude greater than that expected from star formation, strongly suggestive of the presence of an accreting massive black hole. While enhanced X-ray emission and hot dust can be produced in extremely low metallicity environments, J132932.41+323417.0 is not extremely metal poor (~40% solar). With a stellar mass of 2.0 x 10^8 M_Sol, this galaxy is similar in mass to the Small Magellanic Cloud, and is one of the lowest mass galaxies with evidence for a massive nuclear black hole currently known.

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 results from a Chandra study of ultraluminous X-ray sources (ULXs) in a sample of 17 nearby (D_L<60 Mpc) luminous infrared galaxies (LIRGs), selected to have star formation rates (SFRs) in excess of 7 M_sun yr^-1 and low foreground Galactic column densities (N_H < 5*10^20 cm^-2). A total of 53 ULXs were detected and we confirm that this is a complete catalogue of ULXs for the LIRG sample. We examine the evolution of ULX spectra with luminosity by stacking the spectra of individual objects in three luminosity bins, finding a distinct change in spectral index at luminosity ~2 *10^39 erg s^-1. This may be a change in spectrum as 10 M_sun black holes transit from a ~Eddington to a super-Eddington accretion regime, and is supported by a plausible detection of partially-ionised absorption imprinted on the spectrum of the luminous ULX (L_X ~5*10^39 erg s^-1) CXOU J024238.9-000055 in NGC 1068, consistent with the highly ionised massive wind that we would expect to see driven by a super-Eddington accretion flow. This sample shows a large deficit in the number of ULXs detected per unit SFR (0.2 ULXs M_sun^-1 yr^-1). This deficit also manifests itself as a lower differential X-ray luminosity function normalisation for the LIRG sample than for samples of other star forming galaxies. We show that it is unlikely that this deficit is a purely observational effect. Part of this deficit might be attributable to the high metallicity of the LIRGs impeding the production efficiency of ULXs and/or a lag between the star formation starting and the production of ULXs; however, we argue that the evidence -- including very low N_ULX/L_FIR, and an even lower ULX incidence in the central regions of the LIRGs -- shows that the main culprit for this deficit is likely to be the high column of gas and dust in these galaxies, that fuels the high SFR but also acts to obscure many ULXs from our view.

We present a framework for the interpretation of the far-infrared spectra of galaxies in which we have expanded the model parameters compared with previous work by varying the ionization parameter $U$, column density $N$(H), and gas density at the cloud face $n$(H$^{+}$) for a central starburst or AGN. We compare these models carried out with the $Cloudy$ spectral synthesis code to trends in line-to-total far-infrared luminosity ratios, far-infrared fine-structure line ratios, IRAS colors, and OH and H$_{2}$O column densities with trends found in the well-studied sample of ten nearby galaxies from the IRAS Bright Galaxy Sample with infrared luminosities greater than 10$^{10}$ L$_{\odot}$ and IRAS 60 micron fluxes equal to or greater than that of the nearby ULIRG Arp 220. We find that the spectral sequence extending from normal starburst-type emission line spectra to ULIRG-type absorption-dominated spectra with significant absorption from excited levels, can be best explained by simultaneously increasing the hydrogen column density, from as low as 10$^{21}$ cm$^{-2}$ to as high as 10$^{24.8}$ cm$^{-2}$ or greater, and the ionization parameter, from as low as 10$^{-4}$ to as high as 1. The starburst models best reproduce most of the sequence, while AGN models are somewhat better able to produce the high OH and H$_{2}$O column densities in Arp 220. Our results suggest that the molecular ISM in ULIRG-like, absorption-dominated systems is located close to and at least partially obscures the source of power throughout much of the far-infrared, which must be taken into account in order to properly interpret diagnostics of both their sources of power and of feedback.

We present FIR-CO luminosity relations ($\log L_{\rm FIR} = \alpha \log L'_{\rm CO} + \beta$) for the full CO rotational ladder from J=1-0 to J=13-12 for 62 local (z < 0.1) (Ultra) Luminous Infrared Galaxies (LIRGs) using data from Herschel SPIRE-FTS and ground-based telescopes. We extend our sample to high redshifts (z > 1) by including 35 (sub)-millimeter selected dusty star forming galaxies from the literature with robust CO observations. The addition of luminous starbursts at high redshifts enlarge the range of the FIR-CO luminosity relations towards the high-IR-luminosity end while also significantly increasing the small amount of mid-/high-J CO line data available prior to Herschel. This new data-set (both in terms of IR luminosity and J-ladder) reveals linear FIR-CO luminosity relations ($\alpha \sim 1$) for J=1-0 up to J=5-4, with a nearly constant normalisation ($\beta \sim 2$). This is expected from the (also) linear FIR-(molecular line) relations found for the dense gas tracer lines (HCN and CS), as long as the dense gas mass fraction does not vary strongly within our (merger/starburst)-dominated sample. However from J=6-5 and up to J=13-12 we find an increasingly sub-linear slope and higher normalization constant with increasing J. We argue that these are caused by a warm (~100K) and dense ($>10^4{\rm cm^{-3}}$) gas component whose thermal state is unlikely to be maintained by star formation powered far-UV radiation fields (and thus is no longer directly tied to the star formation rate). We suggest that mechanical heating (e.g., supernova driven turbulence and shocks), and not cosmic rays, is the more likely source of energy for this component. The global CO spectral line energy distributions (SLEDs), which remain highly excited from J=6-5 up to J=13-12, are found to be a generic feature of the (U)LIRGs in our sample, and further support the presence of this gas component.

Recent observational studies have shown that H2O emission at (rest) submillimeter wavelengths is ubiquitous in infrared galaxies, both in the local and in the early Universe, suggestive of far-infrared pumping of H2O by dust in warm regions. In this work, models are presented that show that (i) the highest-lying H2O lines (E_{upper}>400 K) are formed in very warm (T_{dust}>~90 K) regions and require high H2O columns (N_{H2O}>~3x10^{17} cm^{-2}), while lower lying lines can be efficiently excited with T_{dust}~45-75 K and N_{H2O}~(0.5-2)x10^{17} cm^{-2}; (ii) significant collisional excitation of the lowest lying (E_{upper}<200 K) levels, which enhances the overall L_{H2O}-L_{IR} ratios, is identified in sources where the ground-state para-H2O 1_{11}-0_{00} line is detected in emission; (iii) the H2O-to-infrared (8-1000 um) luminosity ratio is expected to decrease with increasing T_{dust} for all lines with E_{upper}<~300 K, as has recently been reported in a sample of LIRGs, but increases with T_{dust} for the highest lying H2O lines (E_{upper}>400 K); (iv) we find theoretical upper limits for L_{H2O}/L_{IR} in warm environments, owing to H2O line saturation; (v) individual models are presented for two very different prototypical galaxies, the Seyfert 2 galaxy NGC 1068 and the nearest ultraluminous infrared galaxy Arp 220, showing that the excited submillimeter H2O emission is dominated by far-infrared pumping in both cases; (vi) the L_{H2O}-L_{IR} correlation previously reported in observational studies indicates depletion or exhaustion time scales, t_{dep}=Sigma_{gas}/Sigma_{SFR}, of <~12 Myr for star-forming sources where lines up to E_{upper}=300 K are detected, in agreement with the values previously found for (U)LIRGs from HCN millimeter emission...

In contrast to massive, bulge hosting galaxies, very few supermassive black holes (SMBHs) are known in either low mass, or bulgeless galaxies. Such a population could provide clues to the origins of SMBHs and to secular pathways for their growth. Using the all-sky Wide-Field Infrared Survey Explorer (WISE) survey, and bulge-to-disk decompositions from the Sloan Digital Sky Survey (SDSS) Data Release 7, we report the discovery of a population of local (z<0.3) bulgeless disk galaxies with extremely red mid-infrared colors highly suggestive of a dominant active galactic nucleus (AGN), despite having no optical AGN signatures in their SDSS spectra. Using various mid-infrared selection criteria from the literature, there are between 30 to over 300 bulgeless galaxies with possible AGNs. Other known scenarios that can heat the dust to high temperatures do not appear to explain the observed colors of this sample. If these galaxies are confirmed to host AGNs, this study will provide a breakthrough in characterizing the properties of SMBHs in the low bulge mass regime and in understanding their relation with their host galaxies. Mid-infrared selection identifies AGNs that dominate their host galaxy's emission and therefore reveal a different AGN population than is uncovered by optical studies. We find that the fraction of all galaxies identified as candidate AGNs by WISE is highest at lower stellar masses and drops dramatically in higher mass galaxies, in striking contrast to the findings from optical studies.

This paper describes the data release of the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey conducted between 2005 and 2007. Light curves, spectra, classifications, and ancillary data are presented for 10,258 variable and transient sources discovered through repeat ugriz imaging of SDSS Stripe 82, a 300 deg2 area along the celestial equator. This data release is comprised of all transient sources brighter than r~22.5 mag with no history of variability prior to 2004. Dedicated spectroscopic observations were performed on a subset of 889 transients, as well as spectra for thousands of transient host galaxies using the SDSS-III BOSS spectrographs. Photometric classifications are provided for the candidates with good multi-color light curves that were not observed spectroscopically. From these observations, 4607 transients are either spectroscopically confirmed, or likely to be, supernovae, making this the largest sample of supernova candidates ever compiled. We present a new method for SN host-galaxy identification and derive host-galaxy properties including stellar masses, star-formation rates, and the average stellar population ages from our SDSS multi-band photometry. We derive SALT2 distance moduli for a total of 1443 SN Ia with spectroscopic redshifts as well as photometric redshifts for a further 677 purely-photometric SN Ia candidates. Using the spectroscopically confirmed subset of the three-year SDSS-II SN Ia sample and assuming a flat Lambda-CDM cosmology, we determine Omega_M = 0.315 +/- 0.093 (statistical error only) and detect a non-zero cosmological constant at 5.7 sigmas.

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 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 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 present the first central velocity dispersions (sigma_o) measured from the 0.85 micron Calcium II Triplet (CaT) for 8 advanced (i.e. single nuclei) local (z < 0.15) Ultraluminous Infrared Galaxies (ULIRGs). First, these measurements are used to test the prediction that the "sigma-Discrepancy," in which the CaT sigma_o is systematically larger than the sigma_o obtained from the 1.6 or 2.3 micron stellar CO band-heads, extends to ULIRG luminosities. Next, we combine the CaT data with rest-frame I-band photometry obtained from archival Hubble Space Telescope data and the Sloan Digital Sky Survey (SDSS) to derive dynamical properties for the 8 ULIRGs. These are then compared to the dynamical properties of 9,255 elliptical galaxies from the SDSS within the same redshift volume and of a relatively nearby (z < 0.4) sample of 53 QSO host galaxies. A comparison is also made between the I-band and H-band dynamical properties of the ULIRGs. We find four key results: 1) the sigma-Discrepancy extends to ULIRG luminosities; 2) at I-band ULIRGs lie on the Fundamental Plane (FP) in a region consistent with the most massive elliptical galaxies and not low-intermediate mass ellipticals as previously reported in the near-infrared; 3) the I-band M/L of ULIRGs are consistent with an old stellar population, while at H-band ULIRGs appear significantly younger and less massive; and 4) we derive an I-band Kormendy Relation from the SDSS ellipticals and demonstrate that ULIRGs and QSO host galaxies are dynamically similar.

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

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.

[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.

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]

We study SN 2006oz, a newly-recognized member of the class of H-poor, super-luminous supernovae. We present multi-color light curves from the SDSS-II SN Survey, that cover the rise time, as well as an optical spectrum that shows that the explosion occurred at z~0.376. We fitted black body functions to estimate the temperature and radius evolution of the photosphere and used the parametrized code SYNOW to model the spectrum. We constructed a bolometric light curve and compared it with explosion models. The very early light curves show a dip in the g- and r-bands and a possible initial cooling phase in the u-band before rising to maximum light. The bolometric light curve shows a precursor plateau with a duration of 6-10 days in the rest-frame. A lower limit of M_u < -21.5 can be placed on the absolute peak luminosity of the SN, while the rise time is constrained to be at least 29 days. During our observations, the emitting sphere doubled its radius to 2x10^15 cm, while the temperature remained hot at 15000 K. As for other similar SNe, the spectrum is best modeled with elements including O II and Mg II, while we tentatively suggest that Fe III might be present. We suggest that the precursor plateau might be related to a recombination wave in a circumstellar medium (CSM) and discuss whether this is a common property of all similar explosions. The subsequent rise can be equally well described by input from a magnetar or by ejecta-CSM interaction, but the models are not well constrained owing to the lack of post-maximum observations, and CSM interaction has difficulties accounting for the precursor plateau self-consistently. Radioactive decay is less likely to be the mechanism that powers the luminosity. The host galaxy, detected in deep imaging with the 10 m GTC, is a moderately young and star-forming, but not a starburst, galaxy. It has an absolute magnitude of M_g = -16.9.

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.

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

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

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

We present 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.

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

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.

Mergers in the local universe present a unique opportunity for studying the transformations of galaxies in detail. Presented here are recent results, based on multi-wavelength, high-resolution imaging and medium resolution spectroscopy, which demonstrate how star-formation and the presence of Red Supergiants and/or Asymptotic Giant Branch stars has lead to a serious underestimation of the dynamical masses of infrared-bright galaxies. The dominance of a nuclear disk of young stars in the near-infrared bands, where dust obscuration does not block their signatures, can severely bias the global properties measured in a galaxy, including mass. This explains why past studies of gas-rich Luminous & Ultraluminous Infrared Galaxies, which have measured dynamical masses using the 1.62 or 2.29 micron CO band-heads, have found that these galaxies are forming m < m* ellipticals. On the other hand, precisely because of dust obscuration, I-band photometry and central velocity dispersions obtained with the Calcium II triplet at 0.85 microns reflect the global properties of the mergers and suggest that all types of merger remnants, including infrared-bright ones, will form m > m* ellipticals. Moreover, merger remnants, including LIRGs, are placed on the I-band Fundamental Plane for the first time and appear to be virtually indistinguishable from elliptical galaxies.

Many studies and simulations suggest gas-rich mergers do not contribute significantly to the overall star-formation rate and total mass function of galaxies. The velocity dispersions (sigma) of Luminous & Ultraluminous Infrared Galaxies measured using the 1.62 or 2.29 micron CO bandheads imply they will form m < m* ellipticals. Yet, sigma's obtained with the Calcium II triplet (CaT) at 0.85 micron suggest all types of mergers will form m > m* ellipticals. Presented here are recent results, based on high-resolution imaging and multi-wavelength spectroscopy, which demonstrate the dominance of a nuclear disk of Red Supergiants (RSG) or Asymptotic Giant Branch (AGB) stars in the near-infrared bands, where dust obscuration does not sufficiently block their signatures. The presence of these stars severely biases the dynamical mass. At I-band, where dust can sufficiently block RSG or AGB stars, LIRGs populate the Fundamental Plane over a large dynamic range and are virtually indistinguishable from elliptical galaxies.

The observed faintness of infrared fine-structure line emission along with the warm far-infrared (FIR) colors of ultraluminous infrared galaxies (ULIRGs) is a long-standing problem. In this work, we calculate the line and continuum properties of a cloud exposed to an Active Galactic Nucleus (AGN) and starburst spectral energy distribution (SED). We use an integrated modeling approach, predicting the spectrum of ionized, atomic, and molecular environments in pressure equilibrium. We find that the effects of high ratios of impinging ionizing radiation density to particle density (i.e. high ionization parameters, or U) can reproduce many ULIRG observational characteristics. Physically, as U increases, the fraction of UV photons absorbed by dust increases, corresponding to fewer photons available to photoionize and heat the gas, producing what is known as a "dust-bounded" nebula. We show that high U effects can explain the "[C II] deficit", the ~1 dex drop in the [C II] 158 micron /FIR ratio seen in ULIRGs when compared to starburst or normal galaxies. Additionally, by increasing U through increasing the ionizing photon flux, warmer dust and thus higher IRAS F(60)/F(100) ratios result. High U effects also predict an increase in [O I]63 micron /[C II] 158 micron and a gradual decline in [O III] 88 micron /FIR, similar to the magnitude of the trends observed, and yield a reasonable fit to [Ne V]14 micron /FIR ratio AGN observations.

The South Pole Acoustic Test Setup (SPATS) has been built to evaluate the acoustic characteristics of the South Pole ice in the 10 to 100 kHz frequency range so that the feasibility and specific design of an acoustic neutrino detection array at South Pole can be evaluated. SPATS consists of three vertical strings that were deployed in the upper 400 meters of the South Pole ice cap in January 2007, using the upper part of IceCube holes. The strings form a triangular array with the longest baseline 421 meters. Each of them has 7 stages with one transmitter and one sensor module. Both are equipped with piezoelectric ceramic elements in order to produce or detect sound. Analog signals are brought to the surface on electric cables where they are digitized by a PC-based data acquisition system. The data from all three strings are collected on a master-PC in a central facility, from which they are sent to the northern hemisphere via a satellite link or locally stored on tape. A technical overview of the SPATS detector and its performance is presented.

We present source code for the computer algebra system Mathematica that analyzes the distribution of nearby Galaxies using SDSS data. Download instructions are given, thus within 10 minutes, the reader can verify that galaxies are distributed in an essentially non-homogeneous manner and cluster on 2-dimensional structures. The short code uses a simple method inspired by Minkowski functionals: the distances to the next neighbors are calculated and compared to random distributions in three and two dimensions. The observed distance distribution corresponds clearly to the latter case. The paper may also be helpful for nonexpert scientists to get started with SDSS data analysis.

The ISO/LWS far-infrared spectrum of the ultraluminous galaxy Mkn 231 shows OH and H_2O lines in absorption from energy levels up to 300 K above the ground state, and emission in the [O I] 63 micron and [C II] 158 micron lines. Our analysis shows that OH and H_2O are radiatively pumped by the far-infrared continuum emission of the galaxy. The absorptions in the high-excitation lines require high far-infrared radiation densities, allowing us to constrain the properties of the underlying continuum source. The bulk of the far-infrared continuum arises from a warm (T_dust=70-100 K), optically thick (tau_100micron=1-2) medium of effective diameter 200-400 pc. In our best-fit model of total luminosity L_IR, the observed OH and H2O high-lying lines arise from a luminous (L/L_IR~0.56) region with radius ~100 pc. The high surface brightness of this component suggests that its infrared emission is dominated by the AGN. The derived column densities N(OH)>~10^{17} cm^{-2} and N(H_2O)>~6x10^{16} cm^{-2} may indicate XDR chemistry, although significant starburst chemistry cannot be ruled out. The lower-lying OH, [C II] 158 micron, and [O I] 63 micron lines arise from a more extended (~350 pc) starburst region. We show that the [C II] deficit in Mkn 231 is compatible with a high average abundance of C+ because of an extreme overall luminosity to gas mass ratio. Therefore, a [C II] deficit may indicate a significant contribution to the luminosity by an AGN, and/or by extremely efficient star formation.

Astrophysical neutrinos in the EeV range (particularly those generated by the interaction of cosmic rays with the cosmic microwave background) promise to be a valuable tool to study astrophysics and particle physics at the highest energies. Much could be learned from temporal, spectral, and angular distributions of ~100 events, which could be collected by a detector with ~100 km^3 effective volume in a few years. Scaling the optical Cherenkov technique to this scale is prohibitive. However, using the thick ice sheet available at the South Pole, the radio and acoustic techniques promise to provide sufficient sensitivity with sparse instrumentation. The best strategy may be a hybrid approach combining all three techniques. A new array of acoustic transmitters and sensors, the South Pole Acoustic Test Setup, was installed in three IceCube holes in January 2007. The purpose of SPATS is to measure the attenuation length, background noise, and sound speed for 10-100 kHz acoustic waves. Favorable results would pave the way for a large hybrid array. SPATS is the first array to study the possibility of acoustic neutrino detection in ice, the medium expected to be best for the purpose. First results from SPATS are presented.

We report results of a recently-completed pre-Formulation Phase study of SPIRIT, a candidate NASA Origins Probe mission. SPIRIT is a spatial and spectral interferometer with an operating wavelength range 25 - 400 microns. SPIRIT will provide sub-arcsecond resolution images and spectra with resolution R = 3000 in a 1 arcmin field of view to accomplish three primary scientific objectives: (1) Learn how planetary systems form from protostellar disks, and how they acquire their inhomogeneous composition; (2) characterize the family of extrasolar planetary systems by imaging the structure in debris disks to understand how and where planets of different types form; and (3) learn how high-redshift galaxies formed and merged to form the present-day population of galaxies. Observations with SPIRIT will be complementary to those of the James Webb Space Telescope and the ground-based Atacama Large Millimeter Array. All three observatories could be operational contemporaneously.

We present the first systematic investigation of the [NeV] (14um/24um) and [SIII] (18um/33um) infrared line flux ratios, traditionally used to estimate the density of the ionized gas, in a sample of 41 Type 1 and Type 2 active galactic nuclei (AGNs) observed with the Infrared Spectrograph on board Spitzer. The majority of galaxies with both [NeV] lines detected have observed [NeV] line flux ratios consistent with or below the theoretical low density limit, based on calculations using currently available collision strengths and ignoring absorption and stimulated emission. We find that Type 2 AGNs have lower line flux ratios than Type 1 AGNs and that all of the galaxies with line flux ratios below the low density limit are Type 2 AGNs. We argue that differential infrared extinction to the [NeV] emitting region due to dust in the obscuring torus is responsible for the ratios below the low density limit and we suggest that the ratio may be a tracer of the inclination angle of the torus to our line of sight. Because the temperature of the gas, the amount of extinction, and the effect of absorption and stimulated emission on the line ratios are all unknown, we are not able to determine the electron densities associated with the [NeV] line flux ratios for the objects in our sample. We also find that the [SIII] emission from the galaxies in our sample is extended and originates primarily in star forming regions. Since the emission from low-ionization species is extended, any analysis using line flux ratios from such species obtained from slits of different sizes is invalid for most nearby galaxies.

The SPEAR (or 'FIMS') instrumentation has been used to conduct the first large-scale spectral mapping of diffuse cosmic far ultraviolet (FUV, 900-1750 AA) emission, including important diagnostics of interstellar hot (10^4 K - 10^6 K) and photoionized plasmas, H_2, and dust scattered starlight. The instrumentation's performance has allowed for the unprecedented detection of astrophysical diffuse far UV emission lines. A spectral resolution of 550 and an imaging resolution of 5' is achieved on-orbit in the Short (900 - 1175 AA) and Long (1335 - 1750 AA) bandpass channels within their respective 7.4 deg x 4.3' and 4.0 deg x 4.6' fields of view. We describe the SPEAR imaging spectrographs, their performance, and the nature and handling of their data.

The first complete far-infrared spectrum (43-197um) of the Seyfert 2 galaxy NGC 1068 as observed with the LWS onboard ISO is reported. In addition to the expected ionic fine structure emission lines, the OH rotational lines at 79, 119 and 163um were all detected in emission. The line intensities were modelled together with ISO SWS and optical and UV line intensities from the literature, considering two independent emission components: the AGN component and the starburst component in the circumnuclear ring of 3 kpc in size. Using the UV to mid-IR emission line spectrum to constrain the nuclear ionizing continuum, we have confirmed previous results: a canonical power-law ionizing spectrum is a poorer fit than one with a deep absorption trough, while the presence of a big blue bump is ruled out. Based on starburst synthesis models, we were able to fit the mid- and far-infrared fine-structure line emission, with a low ionization parameter (logU=-3.5) and low densities (n=100 cm^-3). Combining AGN and starburst components, we succeeded in modeling the overall UV to far-IR atomic spectrum of NGC1068, reproducing the line fluxes to within a factor 2.0 on average with a standard deviation of 1.3 and the overall continuum as the sum of the contribution of the thermal dust emission in the ionized and neutral components. TheOH 119 um emission indicates that the line is collisionally excited and arises in a warm and dense region. The OH emission has been modeled using spherically symmetric, non-local, non-LTE radiative transfer models. The models indicate that the bulk of the emission arises from the nuclear region, although some extended contribution from the starburst is not ruled out. The OH abundance in the nuclear region is expected to be $\sim10^{-5}$, characteristic of X-ray dominated regions.

ISO/LWS grating observations of the ultraluminous infrared galaxy Arp 220 shows absorption in molecular lines of OH, H$_2$O, CH, NH, and NH$_3$, as well as in the [O I] 63 $\mu$m line and emission in the [C II] 158 $\mu$m line. We have modeled the continuum and the emission/absorption of all observed features by means of a non-local radiative transfer code. The continuum from 25 to 1300 $\mu$m is modeled as a warm (106 K) nuclear region that is optically thick in the far-infrared, attenuated by an extended region (size 2$''$) that is heated mainly through absorption of nuclear infrared radiation. The molecular absorption in the nuclear region is characterized by high excitation due to the high infrared radiation density. The OH column densities are high toward the nucleus and the extended region. The H$_2$O column density is also high toward the nucleus and lower in the extended region. The column densities in a halo are similar to what are found in the diffuse clouds toward Sgr B2 near the Galactic Center. Most notable are the high column densities found for NH and NH$_3$ toward the nucleus, with values of $\sim1.5\times10^{16}$ cm$^{-2}$ and $\sim3\times10^{16}$ cm$^{-2}$, respectively. A combination of PDRs in the extended region and hot cores with enhanced \hdo photodissociation and a possible shock contribution in the nuclei may explain the relative column densities of OH and \hdo, whereas the nitrogen chemistry may be strongly affected by cosmic ray ionization. The [C II] 158 $\mu$m line is well reproduced by our models and its ``deficit'' relative to the CII/FIR ratio in normal and starburst galaxies is suggested to be mainly a consequence of the dominant non-PDR component of far-infrared radiation, although our models alone cannot rule out extinction effects in the nuclei.

We report 60 and 90 micron observations of 7 millisecond pulsars with ISO. The pulsar PSR B1257+12 is orbited by three planets, and other millisecond pulsars may be orbited by dust disks that represent planets that failed to form or their residue. We do not detect any infrared emission from the 7 pulsars in our sample, and typical upper limits are 100 mJy. Using a simple model, we constrain the typical dust disk mass to be less than of order 100 M_Earth, assuming that the heating of any putative dust disk would be coupled only weakly to the pulsar's emission. If the planets around PSR B1257+12 are composed largely of metals, our limits are probably an order of magnitude above plausible values for the disk mass in metals. Future observations with the Spitzer Space Telescope should be capable of probing into the range of disk masses that could plausibly give rise to planets.

The Infrared Space Observatory Long wavelength Spectrometerhas been used to map distribution of the emission from a sample of 22 atomic, molecular and ionised lines toward the Circumnuclear Disk at the Galactic Centre. The CND disc is clearly seen in the maps of molecular lines such as CO and OH, whilst the central region dominates in other atomic and ionised lines. The spectrum toward Sgr A star is best represented by the sum of a 58 K blackbody, superposed with 22 identifiable emission or absorption features, including four lines each attributed to CO and OH, two broad features that may be indicative of a complex of solid state features, two H2O lines, and the rest being various atomic or ionised atomic lines.

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 report 60 and 90 micron observations of 7 millisecond pulsars with the ISOPHOT instrument and describe our analysis procedures. No pulsars were detected, and typical (3\sigma) upper limits are 150 mJy. We combine our results with others in the literature and use them to place constraints on the existence of protoplanetary or dust disks around millisecond pulsars.

Far-infrared spectroscopy obtained with the ISO LWS has shown that there is strong variation (more than 2 orders of magnitude) in the [C II]/FIR ratios in galaxies extending from blue compact dwarfs, to normal and starburst galaxies, down to elliptical and ultraluminous galaxies (ULIGs). The variation in the relative [C II] 158 micron line strength has been attributed to low metallicity in blue compact dwarfs, high <G_o>/n for normal galaxies and ULIGs, soft radiation fields in ellipticals, and extinction or enhanced abundance of dust in ionized regions in ULIGs. Full ISO/LWS far-infrared (43 - 197 micron) spectra of six nearby IR-bright galaxies reveal a dramatic progression of decreasing fine-structure line emission from ionized species to molecular (OH and H2O) absorption line dominated. The archetypical ULIG, Arp 220, lies at the absorption line dominated end of this sequence. For Arp 220, radiative transfer models indicate that it is optically thick in the FIR and that the water molecules observed in absorption are radiatively excited. If extinction plays a role in the sequence it appears from this analysis that the affected regions are heavily obscured even in the far-infrared, while the detected line emission is not more obscured in ULIGs than in starbursts. Linear correlation between polycyclic aromatic hydrocarbon (PAH) 6.2 micron feature strength and the [C II] 158 micron line strength in starbursts and ULIGs suggests a similar effect for these emitting species, and that the detected PAH emission is not more obscured in ULIGs than in lower luminosity starbursts.

We present the first full FIR spectrum of Centaurus A (NGC 5128) from 43 - 196.7 um. The data was obtained with the ISO Long Wavelength Spectrometer (LWS). We conclude that the FIR emission in a 70 arcsec beam centred on the nucleus is dominated by star formation rather than AGN activity. The flux in the far-infrared lines is ~ 1 % of the total FIR: the [CII] line flux is ~ 0.4 % FIR and the [OI] line is ~ 0.2 %, with the remainder arising from [OIII], [NII] and [NIII] lines. These are typical values for starburst galaxies. The ratio of the [NIII]/[NII] line intensities from the HII regions in the dust lane can be modelled as a ~ 6 million year old starburst. This suggests that the galaxy underwent either a recent merger or a merger which triggered a series of bursts. We estimate that < 5 % of the observed [CII] arises in the cold neutral medium (CNM) and that ~ 10 % arises in the warm ionized medium (WIM). The main contributors to the [CII] emission are the PDRs, which are located throughout the dust lane and in regions beyond where the bulk of the molecular material lies. On scales of ~ 1 kpc the average physical properties of the PDRs are modelled with a gas density, n ~ 1000 cm^-3, an incident far-UV field, G ~ 100 times the local Galactic field, and a gas temperature of ~ 250 K.

Based on far-infrared spectroscopy of a small sample of nearby infrared-bright and ultraluminous infrared galaxies (ULIRGs) with the ISO Long Wavelength Spectrometer, we find a dramatic progression in ionic/atomic fine-structure emission line and molecular/atomic absorption line characteristics in these galaxies extending from strong [O III]52,88 and [N III]57 micron line emission to detection of only faint [C II]158 micron line emission from gas in photodissociation regions in the ULIRGs. The molecular absorption spectra show varying excitation as well, extending from galaxies in which the molecular population mainly occupies the ground state to galaxies in which there is significant population in higher levels. In the case of the prototypical ULIRG, the merger galaxy Arp 220, the spectrum is dominated by absorption lines of OH, H2O, CH, and [O I]. Low [O III]88 micron line flux relative to the integrated far-infrared flux correlates with low excitation and does not appear to be due to far-infrared extinction or to density effects. A progression toward soft radiation fields or very dusty HII regions may explain these effects.

The ROSAT Bright Survey (RBS) aims to completely optically identify the more than 2000 brightest sources detected in the ROSAT all-sky survey at galactic latitudes |b| > 30 degr (excluding LMC, SMC, Virgo cluster). This paper presents a subsample of 66 bright point-like ROSAT survey sources with almost hard PSPC spectra, the hardness ratio HR1 is > 0.5 for most of the sources. The subsample could be nearly completely identified by low-resolution optical spectroscopy with the following breakdown into object classes: 31 Seyfert galaxies, 22 BL Lac candidates, 5 clusters of galaxies 1 cataclysmic variable, and 5 bright stars. Only one object remained unidentified and one X-ray source was a spurious detect ion. The redshift distribution peaks around 0.06 for the Seyferts and around 0.13 for the BL Lac candidates. Observations with medium spectral resolution were obtained for most of the new S eyfert galaxies. A large fraction (20 objects) are type 1 Seyfert galaxies, the other fraction includes Seyfert galaxies of type 1.5 - 1.8 (5 objects), two LINERs, and 4 possible narrow-line Seyfert 1 galaxies (NLS1). About one third of the new Seyfert's have nearby companion galaxies displaying either emission or absorption lines at the same redshift. Among them are a couple of systems showing direct morphological evidence for interaction. The large fraction of interacting galaxies among our sample suggests a scenario where interaction is the main trigger of AGN activity.

We present the first complete far-infrared spectrum (43 to 197 um) of M82, the brightest infrared galaxy in the sky, taken with the Long Wavelength Spectrometer of the Infrared Space Observatory (ISO). We detected seven fine structure emission lines, [OI] 63 and 145 um, [OIII] 52 and 88 um, [NII] 122 um, [NIII] 57 um and [CII] 158 um, and fit their ratios to a combination starburst and photo-dissociation region (PDR) model. The best fit is obtained with HII regions with n = 250 cm^{-3} and an ionization parameter of 10^{-3.5} and PDRs with n = 10^{3.3} cm^{-3} and a far-ultraviolet flux of G_o = 10^{2.8}. We applied both continuous and instantaneous starburst models, with our best fit being a 3-5 Myr old instantaneous burst model with a 100 M_o cut-off. We also detected the ground state rotational line of OH in absorption at 119.4 um. No excited level OH transitions are apparent, indicating that the OH is almost entirely in its ground state with a column density ~ 4x10^{14} cm^{-2}. The spectral energy distribution over the LWS wavelength range is well fit with a 48 K dust temperature and an optical depth, tau_{Dust} proportional to lambda^{-1}.