Using the EMIR instrument on the IRAM 30m telescope, we conducted a spectroscopic redshift search of seven z$_{\rm phot}$ $\sim$ 4 sub-millimetre bright galaxies selected from the Herschel Bright Sources (HerBS) sample with fluxes at 500 $\mu$m greater than 80 mJy. For four sources, we obtained spectroscopic redshifts between 3.4 < z < 4.1 through the detection of multiple CO-spectral lines with J $\leq$ 3. Later, we detected low-J transitions for two of these sources with the GBT including the CO(1-0) transition. For the remaining three sources, more data are needed to determine the spectroscopic redshift unambiguously. The measured CO luminosities and line widths suggest that all these sources are gravitationally lensed. These observations demonstrate that the 2 mm window is indispensable to confirm robust spectroscopic redshifts for z < 4 sources. Finally, we present an efficient graphical method to correctly identify spectroscopic redshifts.

We measure the 850-$\mu$m source densities of 46 candidate protoclusters selected from the Planck High-z catalogue (PHz) and the Planck Catalogue of Compact Sources (PCCS) that were followed up with Herschel-SPIRE and SCUBA-2. This paper aims to search for overdensities of 850-$\mu$m sources in order to select the fields that are most likely to be genuine protoclusters. Of the 46 candidate protoclusters, 25 have significant overdensities ($>$5 times the field counts), 11 have intermediate overdensities (3--5 times the field counts) and 10 have no overdensity ($<$3 times the field counts) of 850-$\mu$m sources. We find that the enhanced number densities are unlikely to be the result of sample variance. Compared with the number counts of another sample selected from Planck's compact source catalogues, this [PHz+PCCS]-selected sample has a higher fraction of candidate protoclusters with significant overdensities, though both samples show overdensities of 850-$\mu$m sources above intermediate level. Based on the estimated star-formation rate densities (SFRDs), we suggest that both samples can efficiently select protoclusters with starbursting galaxies near the redshift at which the global field SFRD peaks ($2 < z < 3$). Based on the confirmation of overdensities found here, future follow-up observations on other PHz targets may greatly increase the number of genuine DSFG-rich clusters/protoclusters.

The fate of cooling gas in the centers of galaxy clusters and groups is still not well understood, as is also the case for the complex processes of triggering star formation in central dominant galaxies (CDGs), re-heating of cooled gas by AGN, and the triggering/feeding of supermassive black hole outbursts. We present CO observations of the early type galaxy NGC 5044, which resides at the center of an X-ray bright group with a moderate cooling flow. For our analysis we combine CO(2-1) data from the 7m antennae of the Atacama Compact Array (ACA), and the ACA total power array (TP). We demonstrate, using the 7m array data, that we can recover the total flux inferred from IRAM 30m single dish observations, which corresponds to a total molecular mass of about 4x10^7 Msun. Most of the recovered flux is blueshifted with respect to the galaxy rest frame and is extended on kpc-scales, suggesting low filling factor dispersed clouds. We find 8 concentrations of molecular gas out to a radius of 10 arcsec (1.5 kpc), which we identify with giant molecular clouds. The total molecular gas mass is more centrally concentrated than the X-ray emitting gas, but extended in the north-east/south-west direction beyond the IRAM 30m beam. We also compare the spatial extent of the molecular gas to the Halpha emission: The CO emission coincides with the very bright Halpha region in the center. We do not detect CO emission in the fainter Halpha regions. Furthermore, we find two CO absorption features spatially located at the center of the galaxy, within 5 pc projected distance of the AGN, infalling at 255 and 265 km/s relative to the AGN. This indicates that the two giant molecular clouds seen in absorption are most likely within the sphere of influence of the supermassive black hole.

We carry out a comprehensive Bayesian correlation analysis between hot halos and direct masses of supermassive black holes (SMBHs), by retrieving the X-ray plasma properties (temperature, luminosity, density, pressure, masses) over galactic to cluster scales for 85 diverse systems. We find new key scalings, with the tightest relation being the $M_\bullet-T_{\rm x}$, followed by $M_\bullet-L_{\rm x}$. The tighter scatter (down to 0.2 dex) and stronger correlation coefficient of all the X-ray halo scalings compared with the optical counterparts (as the $M_\bullet-\sigma_{\rm e}$) suggest that plasma halos play a more central role than stars in tracing and growing SMBHs (especially those that are ultramassive). Moreover, $M_\bullet$ correlates better with the gas mass than dark matter mass. We show the important role of the environment, morphology, and relic galaxies/coronae, as well as the main departures from virialization/self-similarity via the optical/X-ray fundamental planes. We test the three major channels for SMBH growth: hot/Bondi-like models have inconsistent anti-correlation with X-ray halos and too low feeding; cosmological simulations find SMBH mergers as sub-dominant over most of the cosmic time and too rare to induce a central-limit-theorem effect; the scalings are consistent with chaotic cold accretion (CCA), the rain of matter condensing out of the turbulent X-ray halos that sustains a long-term self-regulated feedback loop. The new correlations are major observational constraints for models of SMBH feeding/feedback in galaxies, groups, and clusters (e.g., to test cosmological hydrodynamical simulations), and enable the study of SMBHs not only through X-rays, but also via the Sunyaev-Zel'dovich effect (Compton parameter), lensing (total masses), and cosmology (gas fractions).

We present a combined radio/X-ray analysis of the poorly studied galaxy cluster Abell 2495 (z=0.07923) based on new EVLA and Chandra data. We also analyze and discuss Halpha emission and optical continuum data retrieved from the literature. We find an offset of 6 kpc between the cluster BCG (MCG+02-58-021) and the peak of the X-ray emission, suggesting that the cooling process is not taking place on the central galaxy nucleus. We propose that sloshing of the ICM could be responsible for this separation. Furthermore, we detect a second, 4 kpc offset between the peak of the Halpha emission and that of the X-ray emission. Optical images highlight the presence of a dust filament extending up to 6 kpc in the cluster BCG, and allow us to estimate a dust mass within the central 7 kpc of 1.7e+5 Msun. Exploiting the dust to gas ratio and the L_Halpha-M_mol relation, we argue that a significant amount (up to 10^9 Msun) of molecular gas should be present in the BCG of this cluster. We also investigate the presence of ICM depressions, finding two putative systems of cavities; the inner pair is characterized by t age = 18 Myr and P cav = 1.2e+43 erg/s, the outer one by t age = 53 Myr and P cav = 5.6e+42 erg/s. Their age difference appears to be consistent with the free-fall time of the central cooling gas and with the offset timescale estimated with the Halpha kinematic data, suggesting that sloshing is likely playing a key role in this environment. Furthermore, the cavities' power analysis shows that the AGN energy injection is able to sustain the feedback cycle, despite cooling being offset from the BCG nucleus.

We present SCUBA-2 850-$\mu$m observations of 13 candidate starbursting protoclusters selected using Planck and Herschel data. The cumulative number counts of the 850-$\mu$m sources in 9/13 of these candidate protoclusters show significant overdensities compared to the field, with the probability $<$10$^{-2}$ assuming the sources are randomly distributed in the sky. Using the 250-, 350-, 500- and 850-$\mu$m flux densities, we estimate the photometric redshifts of individual SCUBA-2 sources by fitting spectral energy distribution (SED) templates with an MCMC method. The photometric redshift distribution, peaking at $2<z<3$, is consistent with that of known $z>2$ protoclusters and the peak of the cosmic star-formation rate density (SFRD). We find that the 850-$\mu$m sources in our candidate protoclusters have infrared luminosities of $L_{\mathrm{IR}}\gtrsim$10$^{12}L_{\odot}$ and star-formation rates of SFR=(500-1,500)$M_{\odot}$yr$^{-1}$. By comparing with results in the literature considering only Herschel photometry, we conclude that our 13 candidate protoclusters can be categorised into four groups: six of them being high-redshift starbursting protoclusters, one being a lower-redshift cluster/protocluster, three being protoclusters that contain lensed DSFG(s) or are rich in 850-$\mu$m sources, and three regions without significant Herschel or SCUBA-2 source overdensities. The total SFRs of the candidate protoclusters are found to be comparable or higher than those of known protoclusters, suggesting our sample contains some of the most extreme protocluster population. We infer that cross-matching Planck and Herschel data is a robust method for selecting candidate protoclusters with overdensities of 850-$\mu$m sources.

A large fraction of the baryons and most of the metals in the Universe are unaccounted for. They likely lie in extended galaxy halos, galaxy groups, and the cosmic web, and measuring their nature is essential to understanding galaxy formation. These environments have virial temperatures >10^5.5 K, so the gas should be visible in X-rays. Here we show the breakthrough capabilities of grating spectrometers to 1) detect these reservoirs of hidden metals and mass, and 2) quantify hot gas flows, turbulence, and rotation around the Milky Way and external galaxies. Grating spectrometers are essential instruments for future X-ray missions, and existing technologies provide 50-1500-fold higher throughput compared to current orbiting instruments.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

The HIgh-Resolution Mid-infrarEd Spectrometer (HIRMES) is the 3rd Generation Instrument for the Stratospheric Observatory For Infrared Astronomy (SOFIA), currently in development at the NASA Goddard Space Flight Center (GSFC), and due for commissioning in 2019. By combining direct-detection Transition Edge Sensor (TES) bolometer arrays, grating-dispersive spectroscopy, and a host of Fabry-Perot tunable filters, HIRMES will provide the ability for High Resolution (R~100,000), Mid-Resolution (R~10,000), and Low-Resolution (R~600) slit-spectroscopy, and 2D Spectral Imaging (R~2000 at selected wavelengths) over the 25 - 122 {\mu}m mid-far infrared waveband. The driving science application is the evolution of proto-planetary systems via measurements of water-vapor, water-ice, deuterated hydrogen (HD), and neutral oxygen lines. However, HIRMES has been designed to be as flexible as possible to cover a wide range of science cases that fall within its phase-space, all whilst reaching sensitivities and observing powers not yet seen thus far on SOFIA, providing unique observing capabilities which will remain unmatched for decades.

We present the detection at 89 $\mu$m (observed frame) of the {\it Herschel}-selected gravitationally lensed starburst galaxy HATLASJ1429-0028 (also known as G15v2.19) in 15 minutes with the High-resolution Airborne Wideband Camera-plus (HAWC+) onboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). The spectacular lensing system consists of an edge-on foreground disk galaxy at $z$ = 0.22 and a nearly complete Einstein ring of an intrinsic ultra-luminous infrared galaxy at $z$ = 1.03. Is this high IR luminosity powered by pure star formation (SF) or also an active galactic nucleus (AGN)? Previous nebular line diagnostics indicate that it is star-formation dominated. We perform a 27-band multi-wavelength spectral energy distribution modeling (SED) including the new SOFIA/HAWC+ data to constrain the fractional AGN contribution to the total IR luminosity. The AGN fraction in the IR turns out to be negligible. In addition, J1429-0028 serves as a testbed for comparing SED results from different models/templates and SED codes (MAGPHYS, SED3FIT, and CIGALE). We stress that star formation history is the dominant source of uncertainty in the derived stellar mass (as high as a factor of $\sim$ 10) even in the case of extensive photometric coverage. Furthermore, the detection of a source at $z$ $\sim$ 1 with SOFIA/HAWC+ demonstrates the potential of utilizing this facility for distant galaxy studies including the decomposition of SF/AGN components, which cannot be accomplished with other current facilities.

We present new ALMA CO(2--1) observations of two well studied group-centered elliptical galaxies: NGC~4636 and NGC~5846. In addition, we include a revised analysis of Cycle 0 ALMA observations of the central galaxy in the NGC~5044 group that has been previously published. We find evidence that molecular gas, in the form of off-center orbiting clouds, is a common presence in bright group-centered galaxies (BGG). CO line widths are $\gtrsim 10$ times broader than Galactic molecular clouds, and using the reference Milky Way $X_{CO}$, the total molecular mass ranges from as low as $2.6\times 10^5 M_\odot$ in NGC~4636 to $6.1\times 10^7 M_\odot$ in NGC~5044. With these parameters the virial parameters of the molecular structures is $\gg 1$. Complementary observations of NGC~5846 and NGC~4636 using the ALMA Compact Array (ACA) do not exhibit any detection of a CO diffuse component at the sensitivity level achieved by current exposures. The origin of the detected molecular features is still uncertain, but these ALMA observations suggest that they are the end product of the hot gas cooling process and not the result of merger events. Some of the molecular clouds are associated with dust features as revealed by HST dust extinction maps suggesting that these clouds formed from dust-enhanced cooling. The global nonlinear condensation may be triggered via the chaotic turbulent field or buoyant uplift. The large virial parameter of the molecular structures and correlation with the warm ($10^3 - 10^5 K$)/hot ($\ge10^6$) phase velocity dispersion provide evidence that they are unbound giant molecular associations drifting in the turbulent field, consistently with numerical predictions of the chaotic cold accretion process. Alternatively, the observed large CO line widths may be generated by molecular gas flowing out from cloud surfaces due to heating by the local hot gas atmosphere.

We propose a novel method to constrain turbulence and bulk motions in massive galaxies, groups and clusters, exploring both simulations and observations. As emerged in the recent picture of the top-down multiphase condensation, the hot gaseous halos are tightly linked to all other phases in terms of cospatiality and thermodynamics. While hot halos (10^7 K) are perturbed by subsonic turbulence, warm (10^4 K) ionized and neutral filaments condense out of the turbulent eddies. The peaks condense into cold molecular clouds (< 100 K) raining in the core via chaotic cold accretion (CCA). We show all phases are tightly linked via the ensemble (wide-aperture) velocity dispersion along the line of sight. The correlation arises in complementary long-term AGN feedback simulations and high-resolution CCA runs, and is corroborated by the combined Hitomi and new IFU measurements in Perseus cluster. The ensemble multiphase gas distributions are characterized by substantial spectral line broadening (100-200 km/s) with mild line shift. On the other hand, pencil-beam detections sample the small-scale clouds displaying smaller broadening and significant line shift up to several 100 km/s, with increased scatter due to the turbulence intermittency. We present new ensemble sigma_v of the warm Halpha+[NII] gas in 72 observed cluster/group cores: the constraints are consistent with the simulations and can be used as robust proxies for the turbulent velocities, in particular for the challenging hot plasma (otherwise requiring extremely long X-ray exposures). We show the physically motivated criterion C = t_cool/t_eddy ~ 1 best traces the condensation extent region and presence of multiphase gas in observed clusters/groups. The ensemble method can be applied to many available datasets and can substantially advance our understanding of multiphase halos in light of the next-generation multiwavelength missions.

We present high-spatial-resolution ($\sim 0.12''$ or $\approx 800 \, {\rm pc}$ at $z = 4.5$) ALMA $870\,\mu$m dust continuum observations of a sample of 44 ultrared dusty star-forming galaxies (DSFGs) selected from the H-ATLAS and HerMES far-infrared surveys because of their red colors from 250 to 500 $\mu$m: $S_{500} / S_{250} > 1.5$ and $S_{500} / S_{350} > 1.0$. With photometric redshifts in the range $z \sim 4$-6, our sample includes the most luminous starbursting systems in the early Universe known so far, with total obscured star-formation rates (SFRs) of up to $\sim 4,500 \, M_\odot \, {\rm yr}^{-1}$, as well as a population of lensed, less intrinsically luminous sources. The lower limit on the number of ultrared DSFGs at 870 $\mu$m (with flux densities measured from the ALMA maps and thus not affected by source confusion) derived in this work is in reasonable agreement with models of galaxy evolution, whereas there have been reports of conflicts at 500 $\mu$m (where flux densities are derived from SPIRE). Ultrared DSFGs have a variety of morphologies (from relatively extended disks with smooth radial profiles, to compact sources, both isolated and interacting) and an average size, $\theta_{\rm FWHM}$, of $1.46 \pm 0.41\, {\rm kpc}$, considerably smaller than the values reported in previous work for less-luminous DSFGs at lower redshifts. The size and the estimated gas-depletion times of our sources are compatible with their being the progenitors of the most massive, compact, red-and-dead galaxies at $z \sim 2$-3, and ultimately of local ultra-massive elliptical galaxies or massive galaxy clusters. We are witnessing the birth of the high-mass tail of the red sequence of galaxies.

The atmospheres filling massive galaxies, groups, and clusters display remarkable similarities with rainfalls. Such plasma halos are shaped by AGN heating and subsonic turbulence (~150 km/s), as probed by Hitomi. The new 3D high-resolution simulations show the soft X-ray (< 1 keV) plasma cools rapidly via radiative emission at the high-density interface of the turbulent eddies, stimulating a top-down condensation cascade of warm, $10^4$ K filaments. The ionized (optical/UV) filaments extend up to several kpc and form a skin enveloping the neutral filaments (optical/IR/21-cm). The peaks of the warm filaments further condense into cold molecular clouds (<50 K; radio) with total mass up to several $10^7$ M$_\odot$, i.e., 5/50$\times$ the neutral/ionized masses. The multiphase structures inherit the chaotic kinematics and are dynamically supported. In the inner 500 pc, the clouds collide in inelastic way, mixing angular momentum and leading to chaotic cold accretion (CCA). The BHAR can be modeled via quasi-spherical viscous accretion with collisional mean free path ~100 pc. Beyond the inner kpc region pressure torques drive the angular momentum transport. In CCA, the BHAR is recurrently boosted up to 2 dex compared with the disc evolution, which arises as turbulence is subdominant. The CCA BHAR distribution is lognormal with pink noise power spectrum characteristic of fractal phenomena. The rapid self-similar CCA variability can explain the light curve variability of AGN and HMXBs. An improved criterium to trace thermal instability is proposed. The 3-phase CCA reproduces crucial observations of cospatial multiphase gas in massive galaxies, as Chandra X-ray images, SOAR H$\alpha$ warm filaments and kinematics, Herschel [C$^+$] emission, and ALMA giant molecular associations. CCA plays key role in AGN feedback, AGN unification/obscuration, the evolution of BHs, galaxies, and clusters.

We use high-resolution Herschel-PACS data of 2 nearby elliptical galaxies, IC1459 & NGC2768 to characterize their dust and stellar content. IC1459 & NGC2768 have an unusually large amount of dust for elliptical galaxies (1-3 x 10^5 Msun), this dust is also not distributed along the stellar content. Using data from GALEX (ultraviolet) to PACS (far-infrared), we analyze the spectral energy distribution (SED) of these galaxies with CIGALEMC as a function of the projected position, binning images in 7.2" pixels. From this analysis, we derive maps of SED parameters, such as the metallicity, the stellar mass, the fraction of young star and the dust mass. The larger amount of dust in FIR maps seems related in our model to a larger fraction of young stars which can reach up to 4% in the dustier area. The young stellar population is fitted as a recent (~ 0.5 Gyr) short burst of star formation for both galaxies. The metallicities, which are fairly large at the center of both galaxies, decrease with the radial distance with fairly steep gradient for elliptical galaxies.

In this work we present IRAM-30m telescope observations of a sample of bulge-dominated galaxies with large dust lanes, which have had a recent minor merger. We find these galaxies are very gas rich, with H2 masses between 4x10^8 and 2x10^10 Msun. We use these molecular gas masses, combined with atomic gas masses from an accompanying paper, to calculate gas-to-dust and gas-to-stellar mass ratios. The gas-to-dust ratios of our sample objects vary widely (between ~50 and 750), suggesting many objects have low gas-phase metallicities, and thus that the gas has been accreted through a recent merger with a lower mass companion. We calculate the implied minor companion masses and gas fractions, finding a median predicted stellar mass ratio of ~40:1. The minor companion likely had masses between ~10^7 - 10^10 Msun. The implied merger mass ratios are consistent with the expectation for low redshift gas-rich mergers from simulations. We then go on to present evidence that (no matter which star-formation rate indicator is used) our sample objects have very low star-formation efficiencies (star-formation rate per unit gas mass), lower even than the early-type galaxies from ATLAS3D which already show a suppression. This suggests that minor mergers can actually suppress star-formation activity. We discuss mechanisms that could cause such a suppression, include dynamical effects induced by the minor merger.

We explore the properties of early-type galaxies (ETGs), including ellipticals (E) and lenticulars (S0), in rich environments such as clusters of galaxies (Virgo and Coma). The L_24/L_K distribution of ETGs in both Virgo and Coma clusters shows that some S0s have a much larger L_24/L_K ratio (0.5 to ~2 dex) than the bulk of the ETG population. This could be interpreted as an enhanced star formation rate in these lenticulars. We compare the optical colors of galaxies in these two clusters and investigate the nature of these sources with a large L24/L_K ratio by looking at their spatial distribution within the cluster, by analyzing their optical spectra and by looking at their optical colors compared to late-types. We obtain 10 Coma and 3 Virgo early-type sources with larger L24/L_K ratios than the bulk of their population. We call these sources Mid-Infrared Enhanced Galaxies (MIEGs). In Coma, they are mostly located in the South-West part of the cluster where a substructure is falling onto the main cluster. MIEGs present lower g-r color than the rest of the ETG sample, because of a blue continuum. We interpret the excess L24/L_K ratio as evidence for an enhanced star-formation induced as a consequence of their infall into the main cluster.

Chaotic cold accretion (CCA) profoundly differs from classic black hole accretion models. Using 3D high-resolution simulations, we probe the impact of rotation on the hot and cold accretion flow in a typical massive galaxy. In the hot mode, with or without turbulence, the pressure-dominated flow forms a geometrically thick rotational barrier, suppressing the accretion rate to ~1/3 of the Bondi rate. When radiative cooling is dominant, the gas loses pressure support and quickly circularizes in a cold thin disk. In the more common state of a turbulent and heated atmosphere, CCA drives the dynamics if the gas velocity dispersion exceeds the rotational velocity, i.e., turbulent Taylor number < 1. Extended multiphase filaments condense out of the hot phase via thermal instability and rain toward the black hole, boosting the accretion rate up to 100 times the Bondi rate. Initially, turbulence broadens the angular momentum distribution of the hot gas, allowing the cold phase to condense with prograde or retrograde motion. Subsequent chaotic collisions between the cold filaments, clouds, and a clumpy variable torus promote the cancellation of angular momentum, leading to high accretion rates. The simulated sub-Eddington accretion rates cover the range inferred from AGN cavity observations. CCA predicts inner flat X-ray temperature and $r^{-1}$ density profiles, as recently discovered in M 87 and NGC 3115. The synthetic H{\alpha} images reproduce the main features of cold gas observations in massive ellipticals, as the line fluxes and the filaments versus disk morphology. Such dichotomy is key for the long-term AGN feedback cycle. As gas cools, filamentary CCA develops and boosts AGN heating; the cold mode is thus reduced and the rotating disk remains the sole cold structure. Its consumption leaves the atmosphere in hot mode with suppressed accretion and feedback, reloading the cycle.

Expanding X-ray cavities observed in hot gas atmospheres of many galaxy groups and clusters generate shock waves and turbulence that are primary heating mechanisms required to avoid uninhibited radiatively cooling flows which are not observed. However, we show here that the evolution of buoyant cavities also stimulates radiative cooling of observable masses of low-temperature gas. During their early evolution, radiative cooling occurs in the wakes of buoyant cavities in two locations: in thin radial filaments parallel to the buoyant velocity and more broadly in gas compressed beneath rising cavities. Radiation from these sustained compressions removes entropy from the hot gas. Gas experiencing the largest entropy loss cools first, followed by gas with progressively less entropy loss. Most cooling occurs at late times, $\sim 10^8-10^9$ yrs, long after the X-ray cavities have disrupted and are impossible to detect. During these late times, slightly denser low entropy gas sinks slowly toward the centers of the hot atmospheres where it cools intermittently, forming clouds near the cluster center. Single cavities of energy $10^{57}-10^{58}$ ergs in the atmosphere of the NGC 5044 group create $10^8 - 10^9$ $M_{\odot}$ of cooled gas, exceeding the mass of extended molecular gas currently observed in that group. The cooled gas clouds we compute share many attributes with molecular clouds recently observed in NGC 5044 with ALMA: self-gravitationally unbound, dust-free, quasi-randomly distributed within a few kpc around the group center.

A short 30 minute ALMA observation of the early-type galaxy NGC 5044, which resides at the center of an X-ray bright group with a moderate cooling flow, has detected 24 molecular structures within the central 2.5 kpc. The masses of the molecular structures vary from 3e5 to 1e7 Mo3 and the CO(2-1) linewidths vary from 15 to 65 km/s. Given the large CO(2-1) linewidths, the observed structures are likely giant molecular associations (GMAs) and not individual molecular clouds (GMCs). Only a few of the GMAs are spatially resolved with the cycle 0 ALMA beam and the average density of these GMAs yields a GMC volume filling factor of about 15%. The observed masses of the resolved GMAs are insufficient for them to be gravitationally bound, however, the most massive GMA does contain a less massive component with a linewidth of 5.5 km/s (typical of an individual virialized GMC). We also show that the GMAs cannot be pressure confined by the hot gas. Given the observed CO(2-1) linewidths of the GMAs (i.e., the velocity dispersion of the embedded GMCs) they will likely disperse on a timescale of about 12 Myr, which is less than the central cooling time of the hot gas, so the embedded GMCs within a GMA must condense out of the hot gas at the same time and arise from local concentrations of thermally unstable parcels of hot gas. There are no indications of any disk-like molecular structures and all indications suggest that the molecular gas follows ballistic trajectories after condensing out of the thermally unstable hot gas. The 230 GHz luminosity of the central continuum source is 500 times greater than its low frequency radio luminosity and probably reflects a recent accretion event by the central supermassive black hole. The spectrum of the central continuum source also exhibits an absorption feature.

The Stratospheric Observatory for Infrared Astronomy (SOFIA) has recently concluded a set of engineering flights for Observatory performance evaluation. These in-flight opportunities are viewed as a first comprehensive assessment of the Observatory's performance and are used to guide future development activities, as well as to identify additional Observatory upgrades. Pointing stability was evaluated, including the image motion due to rigid-body and flexible-body telescope modes as well as possible aero-optical image motion. We report on recent improvements in pointing stability by using an active mass damper system installed on the telescope. Measurements and characterization of the shear layer and cavity seeing, as well as image quality evaluation as a function of wavelength have also been performed. Additional tests targeted basic Observatory capabilities and requirements, including pointing accuracy, chopper evaluation and imager sensitivity. This paper reports on the data collected during these flights and presents current SOFIA Observatory performance and characterization.

We compute the properties of a sample of 221 local early-type galaxies with a spectral energy distribution (SED) modeling software, CIGALEMC. Concentrating on the star forming activity and dust contents, we derive parameters such as the specific star formation rate, the dust luminosity, dust mass and temperature. 52 % of our sample is composed of elliptical (E) galaxies and 48 % of lenticular (S0) galaxies. We find a larger proportion of S0 galaxies among galaxies with a large specific star formation rate (sSFR) and large specific dust emission. The stronger activity of S0 galaxies is confirmed by larger dust masses. We investigate the relative proportion of active galactic nuclei (AGN) and star-forming (SF) galaxies in our sample using spectroscopic SDSS data and near-infrared selection techniques, and find a larger proportion of AGN dominated galaxy in the S0 sample than the E one. This could corroborate a scenario where blue galaxies evolve into red ellipticals by passing through a S0 AGN active period while quenching its star formation. Finally, we find a good agreement comparing our estimates with color indicators.

It is commonly thought that AGN feedback can break the self-similar scaling relations of galaxies, groups, and clusters. Using high-resolution 3D hydrodynamic simulations, we isolate the impact of AGN feedback on the $L_{\rm x}-T_{\rm x} $ relation, testing the two archetypal and common regimes, self-regulated mechanical feedback and a quasar thermal blast. We find that AGN feedback has severe difficulty in breaking the relation in a consistent way. The similarity breaking is directly linked to the gas evacuation within $R_{500}$, while the central cooling times are inversely proportional to the core density. Breaking self-similarity implies thus breaking the cool core, morphing all systems to non-cool-core objects, which is in clear contradiction with the observed data populated by several cool-core systems. Self-regulated feedback, which quenches cooling flows and preserves cool cores, prevents the dramatic evacuation and similarity breaking at any scale; the relation scatter is also limited. The impulsive thermal blast can break the core-included $L_{\rm x}-T_{\rm x}$ at $T_{500} < 1$ keV, but substantially empties and overheats the halo, generating a perennial non-cool-core group, as experienced by cosmological simulations. Even with partial evacuation, massive systems remain overheated. We show the action of purely AGN feedback is to lower the luminosity and heating the gas, perpendicular to the fit.

We report on deep near-infrared observations obtained with the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST) of the first five confirmed gravitational lensing events discovered by the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We succeed in disentangling the background galaxy from the lens to gain separate photometry of the two components. The HST data allow us to significantly improve on previous constraints of the mass in stars of the lensed galaxy and to perform accurate lens modelling of these systems, as described in the accompanying paper by Dye et al. We fit the spectral energy distributions of the background sources from near-IR to millimetre wavelengths and use the magnification factors estimated by Dye et al. to derive the intrinsic properties of the lensed galaxies. We find these galaxies to have star-formation rates of approximately 400 to 2000 M_sol/yr, with approximately (6-25)x10^10 M_sol of their baryonic mass already turned into stars. At these rates of star formation, all remaining molecular gas will be exhausted in less than 100 Myr, reaching a final mass in stars of a few 10^11 M_sol. These galaxies are thus proto-ellipticals caught during their major episode of star formation, and observed at the peak epoch z=1.5-3 of the cosmic star formation history of the Universe.

The Hubble morphological sequence from early to late galaxies corresponds to an increasing rate of specific star formation. The Hubble sequence also follows a banana-shaped correlation between 24 and 70 micron luminosities, both normalized with the K-band luminosity. We show that this correlation is significantly tightened if galaxies with central AGN emission are removed, but the cosmic scatter of elliptical galaxies in both 24 and 70 micron luminosities remains significant along the correlation. We find that the 24 micron variation among ellipticals correlates with stellar metallicity, reflecting emission from hot dust in winds from asymptotic giant branch stars of varying metallicity. Infrared surface brightness variations in elliptical galaxies indicate that the K - 24 color profile is U-shaped for reasons that are unclear. In some elliptical galaxies cold interstellar dust emitting at 70 and 160 microns may arise from recent gas-rich mergers. However, we argue that most of the large range of 70 micron luminosity in elliptical galaxies is due to dust transported from galactic cores by feedback events in (currently IR-quiet) active galactic nuclei. Cooler dusty gas naturally accumulates in the cores of elliptical galaxies due to dust-cooled local stellar mass loss and may accrete onto the central black hole, releasing energy. AGN-heated gas can transport dust in cores 5-10 kpc out into the hot gas atmospheres where it radiates extended 70 micron emission but is eventually destroyed by sputtering. This, and some modest star formation, defines a cycle of dust creation and destruction. Elliptical galaxies evidently undergo large transient excursions in the banana plot in times comparable to the sputtering time or AGN duty cycle, 10 Myrs. Normally regarded as passive, elliptical galaxies are the most active galaxies in the IR color-color correlation.

We present two large, nearby (0.013$\le$z$\le$0.06) samples of Early-Type Galaxies (ETGs): a visually classified sample of 220 ETGs, created using source-matched data from the Galaxy and Mass Assembly (GAMA) database with FIR/sub-mm detections from $Herschel$-ATLAS; and a visually classified sample of 551 ETGs which are undetected with $Herschel$-ATLAS. Active galactic nuclei (AGN) are removed from our samples using optical emission line diagnostics. These samples are scrutinised to determine characteristics of sub-mm detected versus undetected ETGs. We find similarities in the stellar mass distributions of the two ETG samples but testing other properties uncovers significant differences. The sub-mm detected sample is shown to have lower concentration and S\'ersic indices than those in the undetected sample - a result which may be linked to the presence of dust in the former. Optical and UV-optical colours are also shown to be much bluer, indicating that the dust is linked with recent star formation. The intrinsic effective radii are on average 1.5 times larger for the sub-mm detected ETGs. Surface densities and groups data from the GAMA database are examined for the two samples, leading to the conclusion that dusty ETGs inhabit sparser environments than non-dusty ETGs in the nearby universe, although environments of the brightest ETGs are shown to differ the least. Modified Planck functions are fit to the H-ATLAS detected PACS and SPIRE fluxes for ETGs with sub-mm flux densities of at least 3$\sigma$ in the 350$\mu$m SPIRE band, giving a resultant mean cold dust temperature of T$_{d}$=22.1K, with a range of 9-30K. The corresponding mean dust mass is 1.8$\times10^{7}$M$_{\odot}$, with a range of (0.08-35.0)$\times10^{7}$M$_{\odot}$. The dust masses calculated from these fits, normalised by stellar mass, are shown to increase with decreasing stellar mass and bluer colours.

Using IRAM PdBI we report the detection of H2O in six new lensed ultra-luminous starburst galaxies at high redshift, discovered in the Herschel H-ATLAS survey. The sources are detected either in the 2_{02}-1_{11} or 2_{11}-2_{02} H_2O emission lines with integrated line fluxes ranging from 1.8 to 14 Jy.km/s. The corresponding apparent luminosities are mu x L_H2O ~ 3-12 x 10^8 Lo, where mu is the lensing magnification factor (3 < mu < 12). These results confirm that H2O lines are among the strongest molecular lines in such galaxies, with intensities almost comparable to those of the high-J CO lines, and same profiles and line widths (200-900 km/s) as the latter. With the current sensitivity of PdBI, H2O can therefore easily be detected in high-z lensed galaxies (with F(500um) > 100 mJy) discovered in the Herschel surveys. Correcting the luminosities for lensing amplification, L_H2O is found to have a strong dependence on the IR luminosity, varying as ~L_IR^{1.2}. This relation which needs to be confirmed with better statistics, may indicate a role of radiative (IR) excitation of the H2O lines, and implies that high-z galaxies with L_IR >~ 10^13 Lo tend to be very strong emitters in H2O, that have no equivalent in the local universe.

We have measured the clustering properties of low-redshift (z < 0.3) sub-mm galaxies detected at 250 micron in the Herschel-ATLAS Science Demonstration Phase (SDP) field. We selected a sample for which we have high-quality spectroscopic redshifts, obtained from reliably matching the 250-micron sources to a complete (for r < 19.4) sample of galaxies from the GAMA database. Both the angular and spatial clustering strength are measured for all z < 0.3 sources as well as for five redshift slices with thickness delta z=0.05 in the range 0.05 < z < 0.3. Our measured spatial clustering length r_0 is comparable to that of optically-selected, moderately star-forming (blue) galaxies: we find values around 5 Mpc. One of the redshift bins contains an interesting structure, at z = 0.164.

We present new observations from Z-Spec, a broadband 185-305 GHz spectrometer, of five sub-millimeter bright lensed sources selected from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration phase (SDP) catalog. We construct a redshift finding algorithm using combinations of the signal-to-noise of all the lines falling in the Z-Spec bandpass to determine redshifts with high confidence, even in cases where the signal-to-noise in individual lines is low. We measure the dust continuum in all sources and secure CO redshifts for four out of five (z~1.5-3). In one source, SDP.17, we tentatively identify two independent redshifts and a water line, confirmed at z=2.308. Our sources have properties characteristic of dusty starburst galaxies, with magnification-corrected star formation rates of 10^(2-3) M_sun/yr. Lower limits for the dust masses (~a few 10^8 M_sun) and spatial extents (~1 kpc equivalent radius) are derived from the continuum spectral energy distributions, corresponding to dust temperatures between 54 and 69K. In the LTE approximation, we derive relatively low CO excitation temperatures (< 100 K) and optical depths (tau<1). Performing a non-LTE excitation analysis using RADEX, we find that the CO lines measured by Z-Spec (from J=4->3 to 10->9, depending on the galaxy) localize the best solutions to either a high-temperature / low-density region, or a low-temperature / high-density region near the LTE solution, with the optical depth varying accordingly. Observations of additional CO lines, CO(1-0) in particular, are needed to constrain the non-LTE models.

We present a pan-chromatic analysis of an unprecedented sample of 1402 250 micron-selected galaxies at z < 0.5 (mean z = 0.24) from the Herschel-ATLAS survey. We complement our Herschel 100-500 micron data with UV-K-band photometry from the Galaxy And Mass Assembly (GAMA) survey and apply the MAGPHYS energy-balance technique to produce pan-chromatic SEDs for a representative sample of 250 micron selected galaxies spanning the most recent 5 Gyr of cosmic history. We derive estimates of physical parameters, including star formation rates, stellar masses, dust masses and infrared luminosities. The typical H-ATLAS galaxy at z < 0.5 has a far-infrared luminosity in the range 10^10 - 10^12 Lsolar (SFR: 1-50 Msolar/yr) thus is broadly representative of normal star forming galaxies over this redshift range. We show that 250 micron-selected galaxies contain a larger mass of dust at a given infra-red luminosity or star formation rate than previous samples selected at 60 micron from IRAS. We derive typical SEDs for H-ATLAS galaxies, and show that the emergent SED shape is most sensitive to specific star formation rate. The optical-UV SEDs also become more reddened due to dust at higher redshifts. Our template SEDs are significantly cooler than existing infra-red templates. They may therefore be most appropriate for inferring total IR luminosities from moderate redshift submillimetre selected samples and for inclusion in models of the lower redshift submillimetre galaxy populations.

We present high-spatial resolution imaging obtained with the Submillimeter Array (SMA) at 880um and the Keck Adaptive Optics (AO) system at Ks-band of a gravitationally lensed sub-millimeter galaxy (SMG) at z=4.243 discovered in the Herschel-Astrophysical Terahertz Large Area Survey. The SMA data (angular resolution ~0.6") resolve the dust emission into multiple lensed images, while the Keck AO Ks-band data (angular resolution ~0.1") resolve the lens into a pair of galaxies separated by 0.3". We present an optical spectrum of the foreground lens obtained with the Gemini-South telescope that provides a lens redshift of z_lens = 0.595 +/- 0.005. We develop and apply a new lens modeling technique in the visibility plane that shows that the SMG is magnified by a factor of mu = 4.1 +/- 0.2 and has an intrinsic infrared (IR) luminosity of L_IR = (2.1 +/- 0.2) x 10^13 Lsun. We measure a half-light radius of the background source of r_s = 4.4 +/- 0.5 kpc which implies an IR luminosity surface density of Sigma_IR = (3.4 +/- 0.9) x 10^11 Lsun kpc^-2, a value that is typical of z > 2 SMGs but significantly lower than IR luminous galaxies at z~0. The two lens galaxies are compact (r_lens ~ 0.9 kpc) early-types with Einstein radii of theta_E1 = 0.57 +/- 0.01 and theta_E2 = 0.40 +/- 0.01 that imply masses of M_lens1 = (7.4 +/- 0.5) x 10^10 Msun and M_lens2 = (3.7 +/- 0.3) x 10^10 Msun. The two lensing galaxies are likely about to undergo a dissipationless merger, and the mass and size of the resultant system should be similar to other early-type galaxies at z~0.6. This work highlights the importance of high spatial resolution imaging in developing models of strongly lensed galaxies discovered by Herschel.

We use spitzer-IRAC data to identify near-infrared counterparts to submillimeter galaxies detected with Herschel-SPIRE at 250um in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). Using a likelihood ratio analysis we identify 146 reliable IRAC counterparts to 123 SPIRE sources out of the 159. We find that, compared to the field population, the SPIRE counterparts occupy a distinct region of 3.6 and 4.5um color-magnitude space, and we use this property to identify a further 23 counterparts to 13 SPIRE sources. The IRAC identification rate of 86% is significantly higher than those that have been demonstrated with wide-field ground-based optical and near-IR imaging of Herschel fields. We estimate a false identification rate of 3.6%, corresponding to 4 to 5 sources. Among the 73 counterparts that are undetected in SDSS, 57 have both 3.6 and 4.5um coverage. Of these 43 have [3.6] - [4.5]> 0 indicating that they are likely to be at z > 1.4. Thus, ~ 40% of identified SPIRE galaxies are likely to be high redshift (z > 1.4) sources. We discuss the statistical properties of the IRAC-identified SPIRE galaxy sample including far-IR luminosities, dust temperatures, star-formation rates, and stellar masses. The majority of our detected galaxies have 10^10 to 10^11 L_sun total IR luminosities and are not intense starbursting galaxies as those found at z ~ 2, but they have a factor of 2 to 3 above average specific star-formation rates compared to near-IR selected galaxy samples.

We identify near-infrared Ks band counterparts to Herschel-ATLAS sub-mm sources, using a preliminary object catalogue from the VISTA VIKING survey. The sub-mm sources are selected from the H-ATLAS Phase 1 catalogue of the GAMA 9h field, which includes all objects detected at 250, 350 or 500 um with the SPIRE instrument. We apply and discuss a likelihood ratio (LR) method for VIKING candidates within a search radius of 10" of the 22,000 SPIRE sources with a 5 sigma detection at 250 um. We find that 11,294(51%) of the SPIRE sources have a best VIKING counterpart with a reliability $R\ge 0.8$, and the false identification rate of these is estimated to be 4.2%. We expect to miss ~5% of true VIKING counterparts. There is evidence from Z-J and J-Ks colours that the reliable counterparts to SPIRE galaxies are marginally redder than the field population. We obtain photometric redshifts for ~68% of all (non-stellar) VIKING candidates with a median redshift of 0.405. Comparing to the results of the optical identifications supplied with the Phase I catalogue, we find that the use of medium-deep near-infrared data improves the identification rate of reliable counterparts from 36% to 51%.

A fundamental gap in the current understanding of galaxies concerns the thermodynamical evolution of the ordinary, baryonic matter. On one hand, radiative emission drastically decreases the thermal energy content of the interstellar plasma (ISM), inducing a slow cooling flow toward the centre. On the other hand, the active galactic nucleus (AGN) struggles to prevent the runaway cooling catastrophe, injecting huge amount of energy in the ISM. The present study intends to deeply investigate the role of mechanical AGN feedback in (isolated or massive) elliptical galaxies, extending and completing the mass range of tested cosmic environments. Our previously successful feedback models, in galaxy clusters and groups, demonstrated that AGN outflows, self-regulated by cold gas accretion, are able to properly quench the cooling flow, without destroying the cool core. Via 3D hydrodynamic simulations (FLASH 3.3), including also stellar evolution, we show that massive mechanical AGN outflows can indeed solve the cooling flow problem for the entire life of the galaxy, at the same time reproducing typical observational features and constraints, such as buoyant underdense bubbles, elliptical shock cocoons, sonic ripples, dredge-up of metals, subsonic turbulence, and extended filamentary or nuclear cold gas. In order to avoid overheating and totally emptying the isolated galaxy, the frequent mechanical AGN feedback should be less powerful and efficient (~1.e-4), compared to the heating required for more massive and bound ellipticals surrounded by the intragroup medium (efficiency ~1.e-3).

The Stratospheric Observatory for Infrared Astronomy (SOFIA) is an airborne observatory consisting of a specially modified Boeing 747SP with a 2.7-m telescope, flying at altitudes as high as 13.7 km (45,000 ft). Designed to observe at wavelengths from 0.3 micron to 1.6 mm, SOFIA operates above 99.8 % of the water vapor that obscures much of the infrared and submillimeter. SOFIA has seven science instruments under development, including an occultation photometer, near-, mid-, and far-infrared cameras, infrared spectrometers, and heterodyne receivers. SOFIA, a joint project between NASA and the German Aerospace Center DLR, began initial science flights in 2010 December, and has conducted 30 science flights in the subsequent year. During this early science period three instruments have flown: the mid-infrared camera FORCAST, the heterodyne spectrometer GREAT, and the occultation photometer HIPO. This article provides an overview of the observatory and its early performance.

We report measurements of the carbon monoxide ground state rotational transition (12C16O J = 1--0) with the Zpectrometer ultra-wideband spectrometer on the 100-m diameter Green Bank Telescope. The sample comprises 11 galaxies with redshifts between z = 2.1 and 3.5 from a total sample of 24 targets identified by Herschel-ATLAS photometric colors from the SPIRE instrument. Nine of the CO measurements are new redshift determinations, substantially adding to the number of detections of galaxies with rest-frame peak submillimeter emission near 100um. The CO detections confirm the existence of massive gas reservoirs within these luminous dusty star-forming galaxies (DSFGs). The CO redshift distribution of the 350um-selected galaxies is strikingly similar to the optical redshifts of 850um-selected submillimeter galaxies (SMGs) in 2.1 < z < 3.5. Spectroscopic redshifts break a temperature-redshift degeneracy; optically thin dust models fit to the far-infrared photometry indicate characteristic dust temperatures near 34 K for most of the galaxies we detect in CO. Detections of two warmer galaxies and statistically significant nondetections hint at warmer or molecule-poor DSFGs with redshifts difficult determine from from Herschel-SPIRE photometric colors alone. Many of the galaxies identified by H-ATLAS photometry are expected to be amplified by foreground gravitational lenses. Analysis of CO linewidths and luminosities provides a method for finding approximate gravitational lens magnifications mu from spectroscopic data alone, yielding mu ~ 3--20. Corrected for magnification, most galaxy luminosities are consistent with an ultra-luminous infrared galaxy (ULIRG) classification, but three are candidate hyper-LIRGs with luminosities greater than 10^13 L_sun.

We describe the infrared properties of sources detected over ~36 deg^2 of sky in the GAMA 15-hr equatorial field, using data from both the Herschel Astrophysical Terahertz Large-Area Survey (H-ATLAS) and Wide-field Infrared Survey (WISE). With 5-sigma point-source depths of 34 and 0.048 mJy at 250 micron and 3.4 micron, respectively, we are able to identify 50.6% of the H-ATLAS sources in the WISE survey, corresponding to a surface density of ~630 deg^{-2}. Approximately two-thirds of these sources have measured spectroscopic or optical/near-IR photometric redshifts of z<1. For sources with spectroscopic redshifts at z<0.3, we find a linear correlation between the infrared luminosity at 3.4 micron and that at 250 micron, with +-50% scatter over ~1.5 orders of magnitude in luminosity, ~10^9 - 10^{10.5} L_sun. By contrast, the matched sources without previously measured redshifts (r>~20.5) have 250-350 micron flux density ratios that suggest either high-redshift galaxies (z>~1.5) or optically faint low-redshift galaxies with unusually low temperatures (T<~20). Their small 3.4-250 micron flux ratios favor a high-redshift galaxy population, as only the most actively star-forming galaxies at low redshift (e.g., Arp 220) exhibit comparable flux density ratios. Furthermore, we find a relatively large AGN fraction (~30%) in a 12 micron flux-limited subsample of H-ATLAS sources, also consistent with there being a significant population of high-redshift sources in the no-redshift sample.

We use the Herschel-ATLAS survey to conduct the first large-scale statistical study of the submm properties of optically selected galaxies. Using ~80,000 r-band selected galaxies from 126 deg^2 of the GAMA survey, we stack into submm imaging at 250, 350 and 500{\mu}m to gain unprecedented statistics on the dust emission from galaxies at z < 0.35. We find that low redshift galaxies account for 5% of the cosmic 250{\mu}m background (4% at 350{\mu}m; 3% at 500{\mu}m), of which approximately 60% comes from 'blue' and 20% from 'red' galaxies (rest-frame g - r). We compare the dust properties of different galaxy populations by dividing the sample into bins of optical luminosity, stellar mass, colour and redshift. In blue galaxies we find that dust temperature and luminosity correlate strongly with stellar mass at a fixed redshift, but red galaxies do not follow these correlations and overall have lower luminosities and temperatures. We make reasonable assumptions to account for the contaminating flux from lensing by red sequence galaxies and conclude that galaxies with different optical colours have fundamentally different dust emission properties. Results indicate that while blue galaxies are more luminous than red galaxies due to higher temperatures, the dust masses of the two samples are relatively similar. Dust mass is shown to correlate with stellar mass, although the dust/stellar mass ratio is much higher for low stellar mass galaxies, consistent with the lowest mass galaxies having the highest specific star formation rates. We stack the 250{\mu}m/NUV luminosity ratio, finding results consistent with greater obscuration of star formation at lower stellar mass and higher redshift. Submm luminosities and dust masses of all galaxies are shown to evolve strongly with redshift, indicating a fall in the amount of obscured star formation in ordinary galaxies over the last four billion years.

While the selection of strongly lensed galaxies with 500{\mu}m flux density S_500>100 mJy has proven to be rather straightforward (Negrello et al. 2010), for many applications it is important to analyze samples larger than the ones obtained when confining ourselves to such a bright limit. Moreover, only by probing to fainter flux densities is possible to exploit strong lensing to investigate the bulk of the high-z star-forming galaxy population. We describe HALOS (the Herschel-ATLAS Lensed Objects Selection), a method for efficiently selecting fainter candidate strongly lensed galaxies, reaching a surface density of ~1.5-2 deg^-2, i.e. a factor of about 4 to 6 higher than that at the 100 mJy flux limit. HALOS will allow the selection of up to ~1000 candidate strongly lensed galaxies (with amplifications \mu>2) over the full H-ATLAS survey area. Applying HALOS to the H-ATLAS Science Demonstration Phase field (~14.4 deg^2) we find 31 candidate strongly lensed galaxies, whose candidate lenses are identified in the VIKING near-infrared catalog. Using the available information on candidate sources and candidate lenses we tentatively estimate a ~72% purity of the sample. The redshift distribution of the candidate lensed sources is close to that reported for most previous surveys for lensed galaxies, while that of candidate lenses extends to substantially higher redshifts than found in the other surveys. The counts of candidate strongly lensed galaxies are also in good agreement with model predictions (Lapi et al. 2011). Even though a key ingredient of the method is the deep near-infrared VIKING photometry, we show that H-ATLAS data alone allow the selection of a similarly deep sample of candidate strongly lensed galaxies with an efficiency close to 50%; a slightly lower surface density (~1.45 deg^-2) can be reached with a ~70% efficiency.

We present the dust properties and star-formation histories of local submillimetre-selected galaxies in Herschel-ATLAS, classified by optical morphology. The early-type galaxies (ETGs) that are detected contain as much dust as typical spirals, and form a unique sample that has been blindly selected at submillimetre wavelengths. Comparing H-ATLAS galaxies to a control sample of optically selected galaxies, we find 5.5% of luminous ETGs are detected in H-ATLAS. The H-ATLAS ETGs contain a significant mass of cold dust: the mean dust mass is 5.5x10^7 Msun, with individual galaxies ranging from 9x10^5-4x10^8 Msun. This is comparable to that of spirals in our sample, and is an order of magnitude more dust than that found for the control ETGs, which have a median dust mass inferred from stacking of (0.8-4.0)x10^6 Msun. The ETGs detected in H-ATLAS have bluer NUV-r colours, higher specific star-formation rates and younger stellar populations than ETGs which are optically selected, and may be transitioning from the blue cloud to the red sequence. We also find that H-ATLAS and control ETGs inhabit similar low-density environments. We conclude that the dust in H-ATLAS and control ETGs cannot be solely from stellar sources, and a large contribution from dust formed in the ISM or external sources is required. Alternatively, dust destruction may not be as efficient as predicted. We also explore the properties of the most passive spiral galaxies in our sample with SSFR<10^-11/yr. We find these passive spirals have lower dust-to-stellar mass ratios, higher stellar masses and older stellar population ages than normal spirals. The passive spirals inhabit low density environments similar to those of the normal spiral galaxies in our sample. This shows that the processes which turn spirals passive do not occur solely in the intermediate density environments of group and cluster outskirts. (Abridged)

Exploiting the Herschel-ATLAS Science Demonstration Phase (SDP) survey data, we have determined the luminosity functions (LFs) at rest-frame wavelengths of 100 and 250 micron and at several redshifts z>1, for bright sub-mm galaxies with star formation rates (SFR) >100 M_sun/yr. We find that the evolution of the comoving LF is strong up to z~2.5, and slows down at higher redshifts. From the LFs and the information on halo masses inferred from clustering analysis, we derived an average relation between SFR and halo mass (and its scatter). We also infer that the timescale of the main episode of dust-enshrouded star formation in massive halos (M_H>3*10^12 M_sun) amounts to ~7*10^8 yr. Given the SFRs, which are in the range 10^2-10^3 M_sun/yr, this timescale implies final stellar masses of order of 10^11-10^12 M_sun. The corresponding stellar mass function matches the observed mass function of passively evolving galaxies at z>1. The comparison of the statistics for sub-mm and UV selected galaxies suggests that the dust-free, UV bright phase, is >10^2 times shorter than the sub-mm bright phase, implying that the dust must form soon after the onset of star formation. Using a single reference Spectral Energy Distribution (SED; the one of the z~2.3 galaxy SMM J2135-0102), our simple physical model is able to reproduce not only the LFs at different redshifts > 1 but also the counts at wavelengths ranging from 250 micron to ~1 mm. Owing to the steepness of the counts and their relatively broad frequency range, this result suggests that the dispersion of sub-mm SEDs of z>1 galaxies around the reference one is rather small.

The Herschel survey, H-ATLAS, with its large areal coverage, has recently discovered a number of bright, strongly lensed high-z submillimeter galaxies. The strong magnification makes it possible to study molecular species other than CO, which are otherwise difficult to observe in high-z galaxies. Among the lensed galaxies already identified by H-ATLAS, the source J090302.9-014127B (SDP.17b) at z = 2.305 is remarkable due to its excitation conditions and a tentative detection of the H2O 202-111 emission line (Lupu et al. 2010). We report observations of this line in SDP.17b using the IRAM interferometer equipped with its new 277- 371GHz receivers. The H2O line is detected at a redshift of z = 2.3049+/-0.0006, with a flux of 7.8+/-0.5 Jy km s-1 and a FWHM of 250+/-60 km s-1. The new flux is 2.4 times weaker than the previous tentative detection, although both remain marginally consistent within 1.6-sigma. The intrinsic line luminosity and ratio of H2O(202-111)/CO8-7 seem comparable with those of the nearby starburst/enshrouded-AGN Mrk 231, suggesting that SDP.17b could also host a luminous AGN. The detection of a strong H2O 202-111 line in SDP.17b implies an efficient excitation mechanism of the water levels that must occur in very dense and warm interstellar gas.

We present Herschel-SPIRE Fourier Transform Spectrometer (FTS) and radio follow-up observations of two Herschel-ATLAS (H-ATLAS) detected strongly lensed distant galaxies. In one of the targeted galaxies H-ATLAS J090311.6+003906 (SDP.81) we detect [OIII] 88\mum and [CII] 158\mum lines at a signal-to-noise ratio of ~5. We do not have any positive line identification in the other fainter target H-ATLAS J091305.0-005343 (SDP.130). Currently SDP.81 is the faintest sub-mm galaxy with positive line detections with the FTS, with continuum flux just below 200 mJy in the 200-600 \mum wavelength range. The derived redshift of SDP.81 from the two detections is z=3.043 +/-0.012, in agreement with ground-based CO measurements. This is the first detection by Herschel of the [OIII] 88\mum line in a galaxy at redshift higher than 0.05. Comparing the observed lines and line ratios with a grid of photo-dissociation region (PDR) models with different physical conditions, we derive the PDR cloud density n ~ 2000 cm^{-3} and the far-UV ionizing radiation field G_0 ~ 200 (in units of the Habing field -- the local Galactic interstellar radiation field of 1.6x10^{-6} W/m^2). Using the CO derived molecular mass and the PDR properties we estimate the effective radius of the emitting region to be 500-700 pc. These characteristics are typical for star-forming, high redshift galaxies. The radio observations indicate that SDP.81 deviates significantly from the local FIR/radio correlation, which hints that some fraction of the radio emission is coming from an AGN. The constraints on the source size from millimiter-wave observations put a very conservative upper limit of the possible AGN contribution to less than 33%. These indications, together with the high [OIII]/FIR ratio and the upper limit of [OI] 63\mum/[CII] 158\mum suggest that some fraction of the ionizing radiation is likely to originate from an AGN.

We introduce a fast Markov Chain Monte Carlo (MCMC) exploration of the astrophysical parameter space using a modified version of the publicly available code CIGALE (Code Investigating GALaxy emission). The original CIGALE builds a grid of theoretical Spectral Energy Distribution (SED) models and fits to photometric fluxes from Ultraviolet (UV) to Infrared (IR) to put contraints on parameters related to both formation and evolution of galaxies. Such a grid-based method can lead to a long and challenging parameter extraction since the computation time increases exponentially with the number of parameters considered and results can be dependent on the density of sampling points, which must be chosen in advance for each parameter. Markov Chain Monte Carlo methods, on the other hand, scale approximately linearly with the number of parameters, allowing a faster and more accurate exploration of the parameter space by using a smaller number of efficiently chosen samples. We test our MCMC version of the code CIGALE (called CIGALEMC) with simulated data. After checking the ability of the code to retrieve the input parameters used to build the mock sample, we fit theoretical SEDs to real data from the well known and studied SINGS sample. We discuss constraints on the parameters and show the advantages of our MCMC sampling method in terms of accuracy of the results and optimization of CPU time.

We present the first direct and unbiased measurement of the evolution of the dust mass function of galaxies over the past 5 billion years of cosmic history using data from the Science Demonstration Phase of the Herschel-ATLAS. The sample consists of galaxies selected at 250{\mu}m which have reliable counterparts from SDSS at z < 0.5, and contains 1867 sources. Dust masses are calculated using both a single temperature grey-body model for the spectral energy distribution and also using a model with multiple temperature components. The dust temperature for either model shows no trend with redshift. Splitting the sample into bins of redshift reveals a strong evolution in the dust properties of the most massive galaxies. At z = 0.4 - 0.5, massive galaxies had dust masses about five times larger than in the local Universe. At the same time, the dust-to-stellar mass ratio was about 3-4 times larger, and the optical depth derived from fitting the UV-sub-mm data with an energy balance model was also higher. This increase in the dust content of massive galaxies at high redshift is difficult to explain using standard dust evolution models and requires a rapid gas consumption timescale together with either a more top-heavy IMF, efficient mantle growth, less dust destruction or combinations of all three. This evolution in dust mass is likely to be associated with a change in overall ISM mass, and points to an enhanced supply of fuel for star formation at earlier cosmic epochs.

We report ground-based follow-up observations of the exceptional source, ID141, one the brightest sources detected so far in the H-ATLAS cosmological survey. ID141 was observed using the IRAM 30-meter telescope and Plateau de Bure interferometer (PdBI), the Submillimeter Array (SMA) and the Atacama Pathfinder Experiment (APEX) submillimeter telescope to measure the dust continuum and emission lines of the main isotope of carbon monoxide and carbon ([C I] and [C II]). The detection of strong CO emission lines with the PdBI confirms that ID141 is at high redshift (z=4.243 +/- 0.001). The strength of the continuum and emission lines suggests that ID141 is gravitationally lensed. The width (Delta V (FWHM) ~ 800 km/s}) and asymmetric profiles of the CO and carbon lines indicate orbital motion in a disc or a merger. The properties derived for ID141 are compatible with a ultraluminous (L_FIR ~ 8.5 +/- 0.3 x 10^13/mu_L Lsun, where mu_L is the amplification factor, dense (n ~ 10^4 cm^-3) and warm (T_kin ~ 40K) starburst galaxy, with an estimated star-formation rate of (0.7 to 1.7) x 10^4/mu_L Msun/yr. The carbon emission lines indicate a dense (n ~ 10^4 cm^-3) Photo-Dominated Region, illuminated by a far-UV radiation field a few thousand times more intense than that in our Galaxy. In conclusion, the physical properties of the high-z galaxy, ID141, are remarkably similar to those of local ultraluminous infrared galaxies.

We present a technique to identify optical counterparts of 250 um-selected sources from the Herschel-ATLAS survey. Of the 6621 250 um > 32 mJy sources in our science demonstration catalogue we find that ~60 percent have counterparts brighter than r=22.4 mag in the Sloan Digital Sky Survey. Applying a likelihood ratio technique we are able to identify 2423 of the counterparts with a reliability R > 0.8. This is approximately 37 percent of the full 250 micron catalogue. We have estimated photometric redshifts for each of these 2423 reliable counterparts, while 1099 also have spectroscopic redshifts collated from several different sources, including the GAMA survey. We estimate the completeness of identifying counterparts as a function of redshift, and present evidence that 250 um-selected Herschel-ATLAS galaxies have a bimodal redshift distribution. Those with reliable optical identifications have a redshift distribution peaking at z ~ 0.25 +/- 0.05, while sub-mm colours suggest that a significant fraction with no counterpart above the r-band limit have z > 1. We also suggest a method for selecting populations of strongly-lensed high redshift galaxies. Our identifications are matched to UV--NIR photometry from the GAMA survey, and these data are available as part of the Herschel-ATLAS public data release.

We investigate the ultraviolet and optical properties and environment of low redshift galaxies detected in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration data. We use the Sloan Digital Sky Survey seventh release and the Galaxy And Mass Assembly database to select galaxies with r_Petro < 19.0 mag in the redshift range 0.02 < z < 0.2 and look for their submillimeter counterparts in H-ATLAS. Our results show that at low redshift, H-ATLAS detects mainly blue/star-forming galaxies with a minor contribution from red systems which are highly obscured by dust. In addition we find that the colour of a galaxy rather than the local density of its environment determines whether it is detectable by H-ATLAS. The average dust temperature of galaxies that are simultaneously detected by both PACS and SPIRE is 25K \pm 4K, independent of environment. This analysis provides a glimpse of the potential of the H-ATLAS data to investigate the submillimeter properties of galaxies in the local universe.

We report the detection of a significant excess in the surface density of far-infrared sources from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) within ~1 Mpc of the centres of 66 optically-selected clusters of galaxies in the SDSS with <z>~0.25. From the analysis of the multiwavelength properties of their counterparts we conclude that the far-infrared emission is associated with dust-obscured star formation and/or active galactic nuclei within galaxies in the clusters themselves. The excess reaches a maximum at a radius of ~0.8 Mpc, where we find 1.0\pm0.3 S_250um>34 mJy sources on average per cluster above what would be expected for random field locations. If the far-infrared emission is dominated by star formation (as opposed to AGN) then this corresponds to an average star formation rate of ~7 M_sun/yr per cluster in sources with L_IR>5d10 L_sun. Although lensed sources make a negligible contribution to the excess signal, a fraction of the sources around the clusters could be gravitationally lensed, and we have identified a sample of potential cases of cluster-lensed Herschel sources that could be targeted in follow-up studies.

The Herschel-ATLAS is a survey of 550 square degrees with the Herschel Space Observatory in five far--infrared and submillimetre bands. The first data for the survey, observations of a field 4x4 sq. degrees in size, were taken during the Science Demonstration Phase, and reach a 5 sigma noise level of 33 mJy/beam at 250 microns. This paper describes the source extraction methods used to create the corresponding Science Demonstration Phase catalogue, which contains 6876 sources, selected at 250 microns, within ~14 sq. degrees. SPIRE sources are extracted using a new method specifically developed for Herschel data; PACS counterparts of these sources are identified using circular apertures placed at the SPIRE positions. Aperture flux densities are measured for sources identified as extended after matching to optical wavelengths. The reliability of this catalogue is also discussed, using full simulated maps at the three SPIRE bands. These show that a significant number of sources at 350 and 500 microns have undergone flux density enhancements of up to a factor of ~2, due mainly to source confusion. Correction factors are determined for these effects. The SDP dataset and corresponding catalogue will be available from http://www.h-atlas.org/.

We use the science demonstration field data of the Herschel-ATLAS to study how star formation, traced by the far-infrared Herschel data, is related to both the accretion luminosity and redshift of quasars selected from the Sloan Digital Sky Survey and the 2SLAQ survey. By developing a maximum likelihood estimator to investigate the presence of correlations between the far-infrared and optical luminosities we find evidence that the star-formation in quasar hosts is correlated with both redshift and quasar accretion luminosity. Assuming a relationship of the form L_IR \propto L_QSO^{\theta} (1 + z)^{\zeta}, we find {\theta} = 0.22 +/- 0.08 and {\zeta} = 1.6 +/- 0.4, although there is substantial additional uncertainty in {\zeta} of order +/- 1, due to uncertainties in the host galaxy dust temperature. We find evidence for a large intrinsic dispersion in the redshift dependence, but no evidence for intrinsic dispersion in the correlation between L_QSO and L_IR, suggesting that the latter may be due to a direct physical connection between star formation and black hole accretion. This is consistent with the idea that both the quasar activity and star formation are dependent on the same reservoir of cold gas, so that they are both affected by the influx of cold gas during mergers or heating of gas via feedback processes.

We use multiwavelength data from the Galaxy And Mass Assembly (GAMA) and Herschel ATLAS (H-ATLAS) surveys to compare the relationship between various dust obscuration measures in galaxies. We explore the connections between the ultraviolet (UV) spectral slope, $\beta$, the Balmer decrement, and the far infrared (IR) to $150\,$nm far ultraviolet (FUV) luminosity ratio. We explore trends with galaxy mass, star formation rate (SFR) and redshift in order to identify possible systematics in these various measures. We reiterate the finding of other authors that there is a large scatter between the Balmer decrement and the $\beta$ parameter, and that $\beta$ may be poorly constrained when derived from only two broad passbands in the UV. We also emphasise that FUV derived SFRs, corrected for dust obscuration using $\beta$, will be overestimated unless a modified relation between $\beta$ and the attenuation factor is used. Even in the optimum case, the resulting SFRs have a significant scatter, well over an order of magnitude. While there is a stronger correlation between the IR to FUV luminosity ratio and $\beta$ parameter than with the Balmer decrement, neither of these correlations are particularly tight, and dust corrections based on $\beta$ for high redshift galaxy SFRs must be treated with caution. We conclude with a description of the extent to which the different obscuration measures are consistent with each other as well as the effects of including other galactic properties on these correlations.

The Stratospheric Observatory for Infrared Astronomy (SOFIA), a joint U.S./German project, is a 2.5-meter infrared airborne telescope carried by a Boeing 747-SP that flies in the stratosphere at altitudes as high as 45,000 feet (13.72 km). This facility is capable of observing from 0.3 {\mu}m to 1.6 mm with an average transmission greater than 80 percent. SOFIA will be staged out of the NASA Dryden Flight Research Center aircraft operations facility at Palmdale, CA. The SOFIA Science Mission Operations (SMO) will be located at NASA Ames Research Center, Moffett Field, CA. First science flights began in 2010 and a full operations schedule of up to one hundred 8 to 10 hour flights per year will be reached by 2014. The observatory is expected to operate until the mid 2030's. SOFIAs initial complement of seven focal plane instruments includes broadband imagers, moderate-resolution spectrographs that will resolve broad features due to dust and large molecules, and high-resolution spectrometers capable of studying the kinematics of atomic and molecular gas at sub-km/s resolution. We describe the SOFIA facility and outline the opportunities for observations by the general scientific community and for future instrumentation development. The operational characteristics of the SOFIA first-generation instruments are summarized. The status of the flight test program is discussed and we show First Light images obtained at wavelengths from 5.4 to 37 \"im with the FORCAST imaging camera. Additional information about SOFIA is available at http://www.sofia.usra.edu and http://www.sofia.usra.edu/Science/docs/SofiaScienceVision051809-1.pdf

We present physical properties of two submillimeter selected gravitationally lensed sources, identified in the Herschel Astrophysical Terahertz Large Area Survey. These submillimeter galaxies (SMGs) have flux densities > 100 mJy at 500 um, but are not visible in existing optical imaging. We fit light profiles to each component of the lensing systems in Spitzer IRAC 3.6 and 4.5 um data and successfully disentangle the foreground lens from the background source in each case, providing important constraints on the spectral energy distributions (SEDs) of the background SMG at rest-frame optical-near-infrared wavelengths. The SED fits show that these two SMGs have high dust obscuration with Av ~4 to 5 and star formation rates of ~100 M_sun/yr. They have low gas fractions and low dynamical masses compared to 850 um selected galaxies.

The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) has uncovered a population of strongly-lensed submillimeter galaxies (SMGs). The Zpectrometer instrument on the Green Bank Telescope (GBT) was used to measure the redshifts and constrain the masses of the cold molecular gas reservoirs for two candidate high-redshift lensed sources. We derive CO(1-0) redshifts of z=3.042+/-0.001 and z=2.625+/-0.001, and measure molecular gas masses of (1--3)x10^{10}Msun, corrected for lens amplification and assuming a conversion factor of alhpa=0.8 Msun(K km/s pc^2)^{-1}. We find typical L(IR)/L'(CO) ratios of 120+/-40 and 140+/-50 Lsun (K km/s pc^2)^{-1}, which are consistent with those found for local ULIRGs and other high-redshift SMGs. From analysis of published data, we find no evidence for enhanced L(IR)/L'(CO(1-0)) ratios for the SMG population in comparison to local ULIRGs. The GBT results highlight the power of using the CO lines to derive blind redshifts, which is challenging for the SMG population at optical wavelengths given their high obscuration.

We measure the projected cross-correlation between low redshift (z < 0.5) far-IR selected galaxies in the SDP field of the Herschel-ATLAS (H-ATLAS) survey and optically selected galaxies from the Galaxy and Mass Assembly (GAMA) redshift survey. In order to obtain robust correlation functions, we restrict the analysis to a subset of 969 out of 6900 H-ATLAS galaxies, which have reliable optical counterparts with r<19.4 mag and well-determined spectroscopic redshifts. The overlap region between the two surveys is 12.6 sq. deg; the matched sample has a median redshift of z ~ 0.2. The cross-correlation of GAMA and H-ATLAS galaxies within this region can be fitted by a power law, with correlation length r_0 ~ 4.63 +/- 0.51 Mpc. Comparing with the corresponding auto-correlation function of GAMA galaxies within the SDP field yields a relative bias (averaged over 2-8 Mpc) of H-ATLAS and GAMA galaxies of b_H/b_G ~ 0.6. Combined with clustering measurements from previous optical studies, this indicates that most of the low redshift H-ATLAS sources are hosted by halos with masses comparable to that of the Milky Way. The correlation function appears to depend on the 250 um luminosity, L_250, with bright (median luminosity \nu L_250 ~ 1.6 x 10^10 L_sun) objects being somewhat more strongly clustered than faint (\nu L_250 ~ 4.0 x 10^9 L_sun) objects. This implies that galaxies with higher dust-obscured star formation rates are hosted by more massive halos.

Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty starforming galaxies. However the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

We study the Spectral Energy Distribution (SED) and the power spectrum of Galactic cirrus emission observed in the 14 deg^2 Science Demonstration Phase field of the Herschel-ATLAS using Herschel and IRAS data from 100 to 500 um. We compare the SPIRE 250, 350 and 500um maps with IRAS 100um emission, binned in 6' pixels. We assume a modified black-body SED with dust emissivity parameter beta (F ~ lambda^(-beta)) and a single dust temperature T_d, and find that the dust temperature and emissivity index varies over the science demonstration field as 10< T_rm < 25 K and 1 < beta< 4. The latter values are somewhat higher than the range of beta often quoted in the literature (1< beta< 2). We estimate the mean values of these parameters to be T_d=19.0 +/- 2.4 K and beta = 1.4 +/- 0.4. In regions of bright cirrus emission, we find that the dust has similar temperatures with T_d = 18.0 +/- 2.5 K, and similar values of beta, ranging from 1.4 +- 0.5 to 1.9+/- 0.5. We show that T_d and beta associated with diffuse cirrus emission are anti-correlated and can be described by the relationship: beta(T_d) = NT_d^alpha with [N=116+/-38, alpha=-1.4+/1 0.1]. The strong correlation found in this analysis is not just limited to high density clumps of cirrus emission as seen in previous studies, but is also seen in diffuse cirrus in low density regions. To provide an independent measure of $T_{\rm d}$ and $\beta$, we obtain the angular power spectrum of the cirrus emission in the {\it IRAS} and SPIRE maps, which is consistent with a power spectrum of the form P(k)=P_0(k/k_0)^gamma where gamma = \^a H R2.6+/-m 0.2 for scales of 50-200' in the SPIRE maps. The cirrus rms fluctuation amplitude at angular scales of 100' is consistent with a modified blackbody SED with T_d = 20.1+/- 0.9 K and beta = 1.3+/- 0.2, in agreement with the values obtained above.

We use the Herschel-ATLAS science demonstration data to investigate the star-formation properties of radio-selected galaxies in the GAMA-9h field as a function of radio luminosity and redshift. Radio selection at the lowest radio luminosities, as expected, selects mostly starburst galaxies. At higher radio luminosities, where the population is dominated by AGN, we find that some individual objects are associated with high far-infrared luminosities. However, the far-infrared properties of the radio-loud population are statistically indistinguishable from those of a comparison population of radio-quiet galaxies matched in redshift and K-band absolute magnitude. There is thus no evidence that the host galaxies of these largely low-luminosity (Fanaroff-Riley class I), and presumably low-excitation, AGN, as a population, have particularly unusual star-formation histories. Models in which the AGN activity in higher-luminosity, high-excitation radio galaxies is triggered by major mergers would predict a luminosity-dependent effect that is not seen in our data (which only span a limited range in radio luminosity) but which may well be detectable with the full Herschel-ATLAS dataset.

We use data from the Herschel-ATLAS to investigate the evolution of the far-infrared--radio correlation over the redshift range 0<z<0.5. Using the total far-infrared luminosity of all >5sigma sources in the Herschel-ATLAS Science Demonstration Field and cross-matching these data with radio data from the Faint Images of the Radio Sky at Twenty-Centimetres (FIRST) survey and the NRAO VLA Northern Sky Survey (NVSS), we obtain 104 radio counterparts to the Herschel sources. With these data we find no evidence for evolution in the far-infrared--radio correlation over the redshift range 0<z<0.5, where the median value for the ratio between far-infrared and radio luminosity, $q_{\rm IR}$, over this range is $q_{\rm IR} = 2.40\pm 0.12$ (and a mean of $q_{\rm IR}=2.52 \pm 0.03$ accounting for the lower limits), consistent with both the local value determined from {\em IRAS} and values derived from surveys targeting the high-redshift Universe. By comparing the radio fluxes of our sample measured from both FIRST and NVSS we show that previous results suggesting an increase in the value of $q_{\rm IR}$ from high to low redshift may be the result of resolving out extended emission of the low-redshift sources with relatively high-resolution interferometric data, although AGN contamination could still play a significant role. We also find tentative evidence that the longer wavelength, cooler dust is heated by an evolved stellar population which does not trace the star-formation rate as closely as the shorter wavelength $\ltsim 250~\mu$m emission or the radio emission, supporting suggestions based on detailed models of individual galaxies.

We describe the reduction of data taken with the PACS instrument on board the Herschel Space Observatory in the Science Demonstration Phase of the Herschel-ATLAS (H-ATLAS) survey, specifically data obtained for a 4x4-deg^2 region using Herschel's fast-scan (60 arcsec/s) parallel mode. We describe in detail a pipeline for data reduction using customised procedures within HIPE from data retrieval to the production of science-quality images. We found that the standard procedure for removing Cosmic-Ray glitches also removed parts of bright sources and so implemented an effective two-stage process to minimise these problems. The pronounced 1/f noise is removed from the timelines using 3.4- and 2.5-arcmin boxcar high-pass filters at 100 and 160-um. Empirical measurements of the point-spread function (PSF) are used to determine the encircled energy fraction as a function of aperture size. For the 100- and 160-um bands, the effective PSFs are ~9 and ~13 arcsec (FWHM), and the 90-per-cent encircled energy radii are 13 and 18 arcsec. Astrometric accuracy is good to ~<2 arcsec. The noise in the final maps is correlated between neighbouring pixels and rather higher than advertised prior to launch. For a pair of cross-scans, the 5-sigma point-source sensitivities are 125-165 mJy for 9-13-arcsec radius apertures at 100-um and 150-240 mJy for 13-18-arcsec radius apertures at 160-um.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We have determined the luminosity function of 250um-selected galaxies detected in the ~14 sq.deg science demonstration region of the Herschel-ATLAS project out to a redshift of z=0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250um luminosity density out to z=0.2 where it is 3.6+1.4-0.9 times higher than the local value.

Aims: We aim to demonstrate that the Herschel ATLAS (H-ATLAS) is suitable for a blind and unbiased survey for debris disks by identifying candidate debris disks associated with main sequence stars in the initial science demonstration field of the survey. We show that H-ATLAS reveals a population of far-infrared/sub-mm sources that are associated with stars or star-like objects on the SDSS main-sequence locus. We validate our approach by comparing the properties of the most likely candidate disks to those of the known population. Methods: We use a photometric selection technique to identify main sequence stars in the SDSS DR7 catalogue and a Bayesian Likelihood Ratio method to identify H-ATLAS catalogue sources associated with these main sequence stars. Following this photometric selection we apply distance cuts to identify the most likely candidate debris disks and rule out the presence of contaminating galaxies using UKIDSS LAS K-band images. Results: We identify 78 H-ATLAS sources associated with SDSS point sources on the main-sequence locus, of which two are the most likely debris disk candidates: H-ATLAS J090315.8 and H-ATLAS J090240.2. We show that they are plausible candidates by comparing their properties to the known population of debris disks. Our initial results indicate that bright debris disks are rare, with only 2 candidates identified in a search sample of 851 stars. We also show that H-ATLAS can derive useful upper limits for debris disks associated with Hipparcos stars in the field and outline the future prospects for our debris disk search programme.

We present colour-colour diagrams of detected sources in the Herschel-ATLAS Science Demonstration Field from 100 to 500 microns using both PACS and SPIRE. We fit isothermal modified black bodies to the spectral energy distribution (SED) to extract the dust temperature of sources with counterparts in Galaxy And Mass Assembly (GAMA) or SDSS surveys with either a spectroscopic or a photometric redshift. For a subsample of 330 sources detected in at least three FIR bands with a significance greater than 3 $\sigma$, we find an average dust temperature of $(28 \pm 8)$K. For sources with no known redshift, we populate the colour-colour diagram with a large number of SEDs generated with a broad range of dust temperatures and emissivity parameters, and compare to colours of observed sources to establish the redshift distribution of this sample. For another subsample of 1686 sources with fluxes above 35 mJy at 350 microns and detected at 250 and 500 microns with a significance greater than 3$\sigma$, we find an average redshift of $2.2 \pm 0.6$.

Aims. The Herschel-ATLAS survey (H-ATLAS) will be the largest area survey to be undertaken by the Herschel Space Observatory. It will cover 550 sq. deg. of extragalactic sky at wavelengths of 100, 160, 250, 350 and 500 microns when completed, reaching flux limits (5 sigma) from 32 to 145mJy. We here present galaxy number counts obtained for SPIRE observations of the first ~14 sq. deg. observed at 250, 350 and 500 microns. Methods. Number counts are a fundamental tool in constraining models of galaxy evolution. We use source catalogs extracted from the H-ATLAS maps as the basis for such an analysis. Correction factors for completeness and flux boosting are derived by applying our extraction method to model catalogs and then applied to the raw observational counts. Results. We find a steep rise in the number counts at flux levels of 100-200mJy in all three SPIRE bands, consistent with results from BLAST. The counts are compared to a range of galaxy evolution models. None of the current models is an ideal fit to the data but all ascribe the steep rise to a population of luminous, rapidly evolving dusty galaxies at moderate to high redshift.

We present a derivation of the star formation rate per comoving volume of quasar host galaxies, derived from stacking analyses of far-infrared to mm-wave photometry of quasars with redshifts 0<z<6 and absolute I-band magnitudes -22>I_AB>-32. We use the science demonstration observations of the first ~16 deg^2 from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) in which there are 240 quasars from the Sloan Digital Sky Survey (SDSS) and a further 171 from the 2dF-SDSS LRG and QSO (2SLAQ) survey. We supplement this data with a compilation of data from IRAS, ISO, Spitzer, SCUBA and MAMBO. H-ATLAS alone statistically detects the quasars in its survey area at >5sigma at 250, 350 and 500um. From the compilation as a whole we find striking evidence of downsizing in quasar host galaxy formation: low-luminosity quasars with absolute magnitudes in the range -22>I_AB>-24 have a comoving star formation rate (derived from 100um rest-frame luminosities) peaking between redshifts of 1 and 2, while high-luminosity quasars with I_AB<-26 have a maximum contribution to the star formation density at z~3. The volume-averaged star formation rate of -22>I_AB>-24 quasars evolves as (1+z)^{2.3 +/- 0.7} at z<2, but the evolution at higher luminosities is much faster reaching (1+z)^{10 +/- 1} at -26>I_AB>-28. We tentatively interpret this as a combination of a declining major merger rate with time and gas consumption reducing fuel for both black hole accretion and star formation.

To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$\mu$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $\gamma$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.

We present measurements of the angular correlation function of galaxies selected from the first field of the H-ATLAS survey. Careful removal of the background from galactic cirrus is essential, and currently dominates the uncertainty in our measurements. For our 250 micron-selected sample we detect no significant clustering, consistent with the expectation that the 250 micron-selected sources are mostly normal galaxies at z<~ 1. For our 350 micron and 500 micron-selected samples we detect relatively strong clustering with correlation amplitudes A of 0.2 and 1.2 at 1', but with relatively large uncertainties. For samples which preferentially select high redshift galaxies at z~2-3 we detect significant strong clustering, leading to an estimate of r_0 ~ 7-11 h^{-1} Mpc. The slope of our clustering measurements is very steep, delta~2. The measurements are consistent with the idea that sub-mm sources consist of a low redshift population of normal galaxies and a high redshift population of highly clustered star-bursting galaxies.

We use Herschel PACS and SPIRE observations of the edge-on spiral galaxy UGC 4754, taken as part of the H-ATLAS SDP observations, to investigate the dust energy balance in this galaxy. We build detailed SKIRT radiative models based on SDSS and UKIDSS maps and use these models to predict the far-infrared emission. We find that our radiative transfer model underestimates the observed FIR emission by a factor of two to three. Similar discrepancies have been found for other edge-on spiral galaxies based on IRAS, ISO, and SCUBA data. Thanks to the good sampling of the SED at FIR wavelengths, we can rule out an underestimation of the FIR emissivity as the cause for this discrepancy. Instead we support highly obscured star formation that contributes little to the optical extinction as a more probable explanation.

We describe the infrared properties of a large sample of early type galaxies, comparing data from the Spitzer archive with Ks-band emission from 2MASS. While most representations of this data result in correlations with large scatter, we find a remarkably tight relation among colors formed by ratios of luminosities in Spitzer-MIPS (24, 70 and 160 um) bands and the Ks-band. Remarkably, this correlation among E and S0 galaxies follows that of nearby normal galaxies of all morphological types. In particular, the tight infrared color-color correlation for S0 galaxies alone follows that of the entire Hubble sequence of normal galaxies, roughly in order of galaxy type from ellipticals to spirals to irregulars. The specific star formation rate of S0 galaxies estimated from the 24um luminosity increases with decreasing Ks-band luminosity (or stellar mass) from essentially zero, as with most massive ellipticals, to rates typical of irregular galaxies. Moreover, the luminosities of the many infrared-luminous S0 galaxies can significantly exceed those of the most luminous (presumably post-merger) E galaxies. Star formation rates in the most infrared-luminous S0 galaxies approach 1-10 solar masses per year. Consistently with this picture we find that while most early-type galaxies populate an infrared red sequence, about 24% of the objects (mostly S0s) are in an infrared blue cloud together with late type galaxies. For those early-type galaxies also observed at radio frequencies we find that the far-infrared luminosities correlate with the mass of neutral and molecular hydrogen, but the scatter is large. This scatter suggests that the star formation may be intermittent or that similar S0 galaxies with cold gaseous disks of nearly equal mass can have varying radial column density distributions that alter the local and global SF rates.

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands.

We present AGN feedback in the interesting cases of two groups: AWM 4 and NGC 5044. AWM 4 is characterized by a combination of properties which seems to defy the paradigm for AGN heating in cluster cores: a flat inner temperature profile indicative of a past, major heating episode which completely erased the cool core, as testified by the high central cooling time (> 3 Gyrs) and by the high central entropy level (~ 50 keV cm^2), and yet an active central radio galaxy with extended radio lobes out to 100 kpc, revealing recent feeding of the central massive black hole. A recent Chandra observation has revealed the presence of a compact cool corona associated with the BCG, solving the puzzle of the apparent lack of low entropy gas surrounding a bright radio source, but opening the question of its origin. NGC 5044 shows in the inner 10 kpc a pair of cavities together with a set of bright filaments. The cavities are consistent with a recent AGN outburst as also indicated by the extent of dust and H_alpha emission even though the absence of extended 1.4 GHz emission remains to be explained. The soft X-ray filaments coincident with H_alpha and dust emission are cooler than those which do not correlate with optical and infrared emission, suggesting that dust-aided cooling can contribute to the overall cooling. For the first time sloshing cold fronts at the scale of a galaxy group have been observed in this object.

We discuss infrared Spitzer observations of early type galaxies in the SAURON sample at 24, 60 and 170 microns. When compared with 2MASS Ks band luminosities, lenticular (S0) galaxies exhibit a much wider range of mid to far-infrared luminosities then elliptical (E) galaxies. Mid and far-infrared emission from E galaxies is a combination of circumstellar or interstellar emission from local mass-losing red giant stars, dust buoyantly transported from the galactic cores into distant hot interstellar gas and dust accreted from the environment. The source of mid and far-IR emission in S0 galaxies is quite different and is consistent with low levels of star formation, 0.02 - 0.2 Msol/yr, in cold, dusty gaseous disks. The infrared 24micron-70micron color is systematically lower for (mostly S0) galaxies with known molecular disks. Our observations support the conjecture that cold dusty gas in some S0 galaxies is created by stellar mass loss at approximately the same rate that it is consumed by star formation, so the mass depletion of these disks by star formation will be slow. Unlike E galaxies, the infrared luminosities of S0 galaxies correlate with both the mass of molecular gas and the stellar Hbeta spectral index, and all are related to the recent star formation rate. However, star formation rates estimated from the Hbeta emission line luminosities L_{Hbeta} in SAURON S0 galaxies are generally much smaller. Since L_{Hbeta} does not correlate with 24 microns emission from dust heated by young stars, optical emission lines appear to be a poor indicator of star formation rates in SAURON S0 galaxies. The absence of Hbeta emission may be due to a relative absence of OB stars in the initial mass function or to dust absorption of Hbeta emission lines.

We present a two-dimensional analysis of the bright nearby galaxy group NGC 5044 using the currently available Chandra and XMM data. In the inner 10 kpc a pair of cavities are evident together with a set of bright X-ray filaments. If the cavities are interpreted as gas displaced by relativistic plasma inflated by an AGN, even in the absence of extended 1.4 GHz emission, this would be consistent with a recent outburst as also indicated by the extent of dust and H_alpha emission. The soft X-ray filaments coincident with H_alpha and dust emission are cooler than the ones which do not correlate with optical and infrared emission. We suggest that dust-aided cooling contributes to form warm (T =10^4 K) gas, emitting H_alpha radiation. At 31 kpc and 67 kpc a pair of cold fronts are present, indicative of sloshing due to a dynamical perturbation caused by accretion of a less massive group, also suggested by the peculiar velocity of the brightest galaxy NGC 5044 with respect to the mean group velocity.

We describe photometry at mid-infrared passbands (1.2 - 24 microns) for a sample of 18 elliptical galaxies. All surface brightness distributions resemble de Vaucouleurs profiles, indicating that most of the emission arises from the photospheres or circumstellar regions of red giant stars. The spectral energy distribution peaks near 1.6 microns, but the half-light or effective radius has a pronounced minimum near the K band (2.15 microns). Apart from the 24 micron passband, all sample-averaged radial color profiles have measurable slopes within about twice the (K band) effective radius. Evidently this variation arises because of an increase in stellar metallicity toward the galactic cores. For example, the sampled-averaged color profile (K - 5.8 microns) has a positive slope although no obvious absorption feature is observed in spectra of elliptical galaxies near 5.8 microns. This, and the minimum in the effective radius, suggests that the K band may be anomalously luminous in metal-rich stars in galaxy cores. Unusual radial color profiles involving the 24 micron passband may suggest that some 24 micron emission comes from interstellar not circumstellar dust grains.

Spitzer observations of extended dust in two optically normal elliptical galaxies provide a new confirmation of buoyant feedback outflow in the hot gas atmospheres around these galaxies. AGN feedback energy is required to prevent wholesale cooling and star formation in these group-centered galaxies. In NGC 5044 we observe interstellar (presumably PAH) emission at 8 microns out to about 5 kpc. Both NGC 5044 and 4636 have extended 70 microns emission from cold dust exceeding that expected from stellar mass loss. The sputtering lifetime of this extended dust in the ~1keV interstellar gas, ~10^7 yrs, establishes the time when the dust first entered the hot gas. Evidently the extended dust originated in dusty disks or clouds, commonly observed in elliptical galaxy cores, that were disrupted, heated and buoyantly transported outward. The surviving central dust in NGC 5044 and 4636 has been disrupted into many small filaments. It is remarkable that the asymmetrically extended 8 micron emission in NGC 5044 is spatially coincident with Halpha+[NII] emission from warm gas. A calculation shows that dust-assisted cooling in buoyant hot gas moving out from the galactic core can cool within a few kpc in about ~10^7 yrs, explaining the optical line emission observed. The X-ray images of both galaxies are disturbed. All timescales for transient activity - restoration of equilibrium and buoyant transport in the hot gas, dynamics of surviving dust fragments, and dust sputtering - are consistent with a central release of feedback energy in both galaxies about 10^7 yrs ago.

Far-infrared Spitzer observations of elliptical galaxies are inconsistent with simple steady state models of dust creation in red giant stars and destruction by grain sputtering in the hot interstellar gas at T ~ 10^7 K. The flux at 24 microns correlates with optical fluxes, suggesting that this relatively hot dust is largely circumstellar. But fluxes at 70 and 160 microns do not correlate with optical fluxes. Elliptical galaxies with similar L_B have luminosities at 70 and 160 microns (L_70 and L_160) that vary over a factor ~ 100, implying an additional source of dust unrelated to that produced by ongoing local stellar mass loss. Neither L_70/L_B nor L_160/L_B correlate with the stellar age or metallicity. Optical line fluxes from warm gas at T ~ 10^4 K correlate weakly with L_70 and L_160, suggesting that the dust may be responsible for cooling this gas. Many normal elliptical galaxies have emission at 70 microns that is extended to 5-10 kpc. Extended far-infrared emission with sputtering lifetimes of ~10^8 yrs is difficult to maintain by mergers with gas-rich galaxies. Instead, we propose that this cold dust is buoyantly transported from reservoirs of dust in the galactic cores which is supplied by mass loss from stars in the core. Intermittent energy outbursts from AGNs can drive the buoyant outflow.

In a sample of thirty normal elliptical galaxies observed with the Spitzer IRS, one galaxy, NGC4697, shows strong PAH emission, but with an apparently weak 7.7 micron feature. We find that the PAH emission is confined to the central regions of the galaxy and that once a quiescent elliptical galaxy spectrum is subtracted, the PAH feature ratios are normal. We show that the PAH spectrum resembles the diffuse ISM of our galaxy rather than an HII region, and is not indicative of a starburst galaxy. We suggest that the PAHs in NGC4697 are consistent with a recent but now past burst of star formation.

The mean ages of early-type galaxies obtained from the analysis of optical spectra, give a mean age of 8 Gyr at z = 0, with 40% being younger than 6 Gyr. Independent age determinations are possible by using infrared spectra (5-21 microns), which we have obtained with the Infrared Spectrograph on the Spitzer Observatory. This age indicator is based on the collective mass loss rate of stars, where mass loss from AGB stars produces a silicate emission feature at 9-12 microns. This feature decreases more rapidly than the shorter wavelength continuum as a stellar population ages, providing an age indicator. From observations of 30 nearby early-type galaxies, 29 show a spectral energy distribution dominated by stars and one has significant emission from the ISM and is excluded. The infrared age indicators for the 29 galaxies show them all to be old, with a mean age of about 10 Gyr and a standard deviation of only a few Gyr. This is consistent with the ages inferred from the values of M/L_B, but is inconsistent with the ages derived from the optical line indices, which can be much younger. All of these age indicators are luminosity-weighted and should be correlated, even if multiple-age components are considered. The inconsistency indicates that there is a significant problem with either the infrared and the M/L_B ages, which agree, or with the ages inferred from the optical absorption lines.

Mid-infrared observations (3.6 - 24 microns) of normal giant elliptical galaxies with the Spitzer space telescope are consistent with pure populations of very old stars with no evidence of younger stars. Most of the stars in giant elliptical galaxies are old but the mean stellar age determined from Balmer absorption in optical spectra can appear much younger due to a small admixture of younger stars. The mean stellar age can also be determined from the spectral energy distribution in the mid-infrared which decreases with time relative to the optical emission and shifts to shorter wavelengths. The observed flux ratios F_8um/F_3.6um and F_24um/F_3.6um for elliptical galaxies with the oldest Balmer line ages are lower than predicted by recent models of single stellar populations. For ellipticals with the youngest Balmer line ages in our sample, 3-5 Gyrs, the flux ratios F_24um/F_3.6um are identical to those of the oldest stars. When theoretical mid-IR spectra of old (12 Gyr) and young stellar populations are combined, errors in the F_24um/F_3.6um observations are formally inconsistent with a mass fraction of young stars that exceeds ~1%. This is less than the fraction of young stars expected in discussions of recent surveys of elliptical galaxies at higher redshifts. However, this inconsistancy between Balmer line ages and those inferred from mid-IR observations must be regarded as provisional until more accurate observations and theoretical spectra become available. Finally, there is no evidence to date that central disks or patches of dust commonly visible in optical images of elliptical galaxies contribute sensibly to the mid-IR spectrum.

We present infrared (IR) spectrophotometry ($R \simeq 180$) of three Herbig Ae/Be stars surrounded by possible protoplanetary disks: HD 150193, HD100546 and HD 179218. We construct a mid-IR spectral energy distributions (SED) for each object by using $7.6 - 13.2$ \micron HIFOGS spectra, 2.4 -- 45 \micron\ spectrophotometry from the {\it ISO} SWS, the 12, 25, 60, and 100 \micron\ photometric points from IRAS, and for HD 179218, photometric bolometric data points from the Mt. Lemmon Observing Facility. The SEDs are modeled by using an expanded version of the \citet{chigol97} two-layer, radiative and hydrostatic equilibrium, passive disk. This expanded version includes the emission from Mg-pure crystalline olivine (forsterite) grains in the disk surface layer. HD 150193 contains no crystals while HD 100546 and HD 179218 respectively show evidence of having crystalline silicates in the surface layers of their disks. We find that the inner region of HD100546 has a 37% higher crystalline-to-amorphous silicate ratio in its inner disk region ($\leq 5$ AU) compared to the outer disk region, while the inner disk region of HD 179218 has a 84% higher crystalline-to-amorphous silicate ratio in its inner disk region ($\leq 5$ AU) compared to the outer region. HD 150193 is best-fit by a small disk ($\sim 5$ AU in radius) while HD 100546 and HD 179218 are best fit by larger disks ($\sim 150$ AU in radius).

The mid-infrared (10-20 um) luminosity of elliptical galaxies is dominated by the integrated emission from circumstellar dust in red giant stars. As a single stellar population evolves, the rate of dusty mass loss from red giant stars decreases with time, so the mid-infrared luminosity should also decline with stellar age. To seek such a correlation, we have used archival ISO observations to determine surface brightness profiles and central fluxes at 15 um in 17 early-type galaxies for which stellar ages have been determined from optical spectral indices. The radial surface brightness distributions at 15 um generally follow the stellar de Vaucouleurs profile as expected. We find that the surface brightness ratio mu_{15um}/mu_{I-band} is systematically higher in elliptical galaxies with ages < 5 Gyrs and in galaxies that exhibit evidence of recent mergers. Within the accuracy of our observations, mu_{15um}/mu_{I-band} shows no age dependence for ages > 5 Gyrs. The corresponding flux ratios F_{15um}/F_{I-band} within apertures scaled to the effective radius (R_e/8) are proportional to the mu_{15um}/mu_{I-band} ratios at larger galactic radii, indicating that no 15 um emission is detected from central dust clouds visible in optical images in some of our sample galaxies. Emission at 15 um is observed in non-central massive clouds of dust and cold gas in NGC1316, an elliptical galaxy that is thought to have had a recent merger. Recent {\it Spitzer Space Telescope} data also indicate the presence of PAH emission at 8 um. Several ellipticals have extended regions of 15 um emission that have no obvious counterparts at other frequencies.

Using spatial-spectral data cubes of reflection nebulae obtained by ISOCAM, we have observed a shift in the central wavelength of the 7.7 um band within several reflection nebulae. This band, composed of components at 7.85 um and 7.65 um, shows a centroid shift from 7.75 um near the edge of the nebulae to 7.65 um towards the center of the nebulae as the shorter wavelength component becomes relatively stronger. The behavior of the 7.7 um band center can be explained either by assuming that anions are the origin of the 7.85 um band and cations the 7.65 um band, or that the band center wavelength depends on the chemical nature of the PAHs. The ratio of the 11.3/7.7 um bands also changes with distance from the central star, first rising from the center towards the edge of the nebula, then falling at the largest distances from the star, consistent with the 11.3/7.7 um band ratio being controlled by the PAH ionization state.

We describe far infrared observations of early-type galaxies selected from the ISO archive. This rather inhomogeneous sample includes 39 giant elliptical galaxies and 14 S0 (or later) galaxies. These galaxies were observed using a variety of different observing modes, each of which requires different and often rather elaborate photometric reduction procedures. The ISO background data agrees well with COBE-DIRBE results to which we have renormalized our calibrations. As a further check, the ISO fluxes from galaxies at 60 and 100 microns agree very well with those previously observed with IRAS at these wavelengths. The spatial resolution of ISO is several times greater than that of IRAS and the ISO observations extend out to 200 microns which views a significantly greater mass of colder dust not assessable to IRAS. Most of the galaxies are essentially point sources at ISO resolution, but a few are clearly extended at FIR wavelengths with image sizes that increase with FIR wavelength. The integrated far infrared luminosities do not correlate with optical luminosities, suggesting that the dust may have an external, merger-related origin. In general the far infrared spectral energy distributions can be modeled with dust at two temperatures, ~ 43 K and ~ 20 K, which probably represent limits of a continuous range of temperatures. The colder dust component dominates the total mass of dust, which is typically more than ten times larger than the dust masses previously estimated for the same galaxies using IRAS observations. For S0 galaxies we find that the optically normalized far-infrared luminosity $L_{FIR} / L_B$ correlates strongly with the mid-infrared luminosity $L_{15\mu\rm{m}} / L_B$, but that correlation is weaker for elliptical galaxies.

We present far-IR ISO observations of three early-type galaxies in the Virgo cluster. The data were recorded using the ISOPHOT instrument in both the P32 oversampled maps and the P39/37 sparse maps. The maps reach the limiting sensitivity of the ISOPHOT instrument at 60, 90, and 180$\mu$m. Two of the most prominent elliptical galaxies in Virgo -- NGC 4472 and NGC 4649 -- clearly show no emission at far-IR wavelengths at a level of few tens of mJy, but NGC 4636 is detected at all three wavelengths. We have computed the far-IR emission from dust for NGC4472 and NGC4636 under the assumption that dusty outflows from evolving red giant stars are continuously supplying dust to the interstellar medium and that the grains, once diffused into the interstellar medium, are sputtered away by collisions with ions. While the calculated fluxes are consistent with the observed upper limits for NGC4472, the dust emission detected in NGC4636 support the hypothesis that additional dust has been accreted in a very recent ($\lta$ few $10^8$ yrs) merger with a dusty, gas-rich galaxy.

We present ISOPHOT P32 oversampled maps and P37/39 sparse maps, of three bright elliptical galaxies in the Virgo Cluster. The maps reach the limiting sensitivity of the ISOPHOT instrument at 60, 100, 170 and 200 microns. Two elliptical galaxies show no emission at all far-IR ISOPHOT wavelengths at a level of few tens of mJy. The null detection provides a test of the evolution of dust in elliptical galaxies and its size distribution and composition. As previous studies have shown, in many elliptical galaxies both IRAS and ISO have detected mid-IR excess 6-15 micron emission relative to the stellar continuum indicating emission from circumstellar dust. Under the assumption that these dusty outflows from evolving red giant stars and planetary nebulae are continuously supplying dust to the interstellar medium, we have computed the infrared luminosity at the ISOPHOT bands appropriate for NGC4472. The null far-IR ISOPHOT observations exceed the far-IR flux expected from dust expelled from a normal old stellar population.

Early-type galaxies exhibit thermal and molecular resonance emission from dust that is shed and heated through stellar mass loss as a subset of the population moves through the AGB phase of evolution. Because this emission can give direct insight into stellar evolution in addition to galactic stellar mass loss and ISM injection rates, we conducted a program to search for this signature emission with CAM on ISO. We obtained 6-15 micron imaging observations in six narrow bands for nine elliptical galaxies; every galaxy is detected in every band. For wavelengths shorter than 9 microns, the spectra are well matched by a blackbody, originating from the K and M stars that dominate the integrated light of elliptical galaxies. However, at wavelengths between 9 and 15 microns, the galaxies display excess emission relative to the stellar photospheric radiation. Additional data taken with the fine resolution circular variable filter on one source clearly shows broad emission from 9 to 15 microns, peaking around 10 microns. This result is consistent with the known, broad silicate feature at 9.7 microns, originating in the circumstellar envelopes of AGB stars. This emission is compared with studies of Galactic and LMC AGB stars to derive cumulative mass loss rates. In general, these mass loss rates agree with the expected ~0.8 solar masses per year value predicted by stellar evolutionary models. Both the photospheric and circumstellar envelope emission follow a de Vaucouleurs' R^{1/4} law, supporting the conclusion that the mid-infrared excess emission originates in the stellar component of the galaxies and acts as a tracer of AGB mass loss and mass injection into the ISM.