The decline of star formation in massive low-redshift galaxies, often referred to as quenching, has been attributed to a variety of factors. Some proposals suggest that erupting active galactic nuclei may strip galaxies of their interstellar medium, and thus the ability to form stars. Here, we note that, whereas star formation is universal in small, low-redshift galaxies, fractional duty cycles of star formation steadily decline in galaxies of increasing mass, although star formation may not cease entirely. We show that, when infall of gas from extragalactic space ceases, galaxies of high stellar mass appear to sustain star formation on gas liberated in mass loss from evolved low- and intermediate-mass stars admixed with occasional Type II supernova ejecta. This model quantitatively accounts for the universal limiting metallicity plateau at a ratio of oxygen to hydrogen atoms, Z(O) = n(O)/n(H) = 0.0013, characterizing high-mass intermittently star-forming galaxies. We show that, when fractional duty cycles are specifically taken into account, the star formation rates of galaxies on this plateau correspond to mass loss rates from evolving stars in rough agreement with observed estimates. Far-infrared continuum and fine-structure line observations, as well as molecular data, may soon be able to resolve whether or not low levels of sporadic star formation can be sustained indefinitely in massive galaxies.

We describe an equilibrium model that links the metallicity of low-redshift galaxies to stellar evolution models. It enables the testing of different stellar initial mass functions and metal yields against observed galaxy metallicities. We show that the metallicities of more than 80,000 Sloan Digital Sky Survey (SDSS) galaxies in the low-redshift range $0.07\leq z\leq 0.3$ considerably constrain stellar evolution models that simultaneously relate galaxy stellar mass, metallicity, and star formation rates (SFRs) to the infall rate of low-metallicity extragalactic gas and outflow of enriched matter. A feature of our model is that it encompasses both the active star forming phases of a galaxy and epochs during which the same galaxy may lie fallow. We show that the galaxy-mass-metallicity-star-formation relation can be traced to infall of extragalactic gas mixing with native gas from host galaxies to form stars of observed metallicities, the most massive of which eject oxygen into extragalactic space. Most consequential among our findings is that, on average, extragalactic infall accounts for one half of the gas required for star formation, a ratio that is remarkably constant across galaxies with stellar masses ranging at least from $M* = 2 \times 10^9$ to $6\times 10^{10} M_{\odot}$. This leads us to propose that star formation is initiated when extragalactic infall roughly doubles the mass of marginally stable interstellar clouds. The processes described may also account quantitatively for the metallicity of extragalactic space, though to check this the fraction of extragalactic baryons will need to be more firmly established.

Context. The high intensities and narrow bandwidths exhibited by some astronomical masers make them ideal tools for studying star-forming giant molecular clouds. The water maser transition $J_{K_{-1}K_{1}}=5_{32}-4_{41}$ at 620.701 GHz can only be observed from above Earth's strongly absorbing atmosphere; its emission has recently been detected from space. Aims. We sought to further characterize the star-forming environment of Orion KL by investigating the linear polarisation of a source emitting a narrow 620.701 GHz maser feature with the heterodyne spectrometer HIFI on board the Herschel Space Observatory. Methods. High-resolution spectral datasets were collected over a thirteen month period beginning in 2011 March, to establish not only the linear polarisation but also the temporal variability of the source. Results. Within a $3\sigma$ uncertainty, no polarisation was detected to an upper limit of approximately 2%. These results are compared with coeval linear polarisation measurements of the 22.235 GHz $J_{K_{-1}K_{1}}=6_{16}-5_{23}$ maser line from the Effelsberg 100-m radio telescope, typically a much stronger maser transition. Although strongly polarised emission is observed for one component of the 22.235 GHz maser at 7.2 km s$^{-1}$, a weaker component at the same velocity as the 620.701 GHz maser at 11.7 km s$^{-1}$ is much less polarised.

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

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

Working with 108,786 Sloan Digital Sky Survey low redshift galaxies we have examined the relation between galaxy mass, metallicity, radius, and star formation rates primarily in the central portions of galaxies. We subdivided the redshift range covered in our sample, 0.07<z<0.3, into three narrower redshift bins, and three sets of radial size. We show that for 72% of the galaxies the observed gas metallicities, Zx, are consistent with (i) a quantitative physical relation for star formation through episodic infall of gas of metallicity Zi = 0.125x10^-3 +/- 1.25x10^-3; (ii) thorough mixing of infalling and native gas before onset of star formation; (iii) a star formation rate (SFR) proportional to the 3/2 power of the infalling mass rate, Mi; and (iv) intermittent quiescent phases devoid of star formation during which the native gas in a galaxy exhibits a characteristic elevated gas metallicity, Z0, dependent on galaxy mass, M*, and a characteristic ratio of stellar mass to native mass of gas, Mg. Most if not all our star-forming galaxies with M* < 2.0x10^10 Msun, and many with M* > 2.0x10^10 Msun and large radii appear fed by infall. Smaller massive galaxies with high Zx and high star formation rates show more complex behavior. A mean-field-theory toy model for the physics of infall accounts for the (SFR) \propto Mi^(3/2) relation and permits us to estimate the mean densities and velocities of clumps of baryonic matter traversing the dark matter halos in which the SDSS galaxies may be embedded.

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

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out deep extragalactic surveys, one of whose aims is to establish spectral energy distributions (SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength region. We report observations of the (IR/SMM) emission from the Lockman North field (LN) and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy images in the wavelength range covered by Herschel generally represent a blend with contributions from neighboring galaxies, we present sets of galaxies in each field especially free of blending at 250, 350, and 500 microns. We identify the cumulative emission of these galaxies and the fraction of the far infrared cosmic background radiation they contribute. Our surveys reveal a number of highly luminous galaxies at redshift z ~< 3 and a novel relationship between infrared and visible emission that shows a dependence on luminosity and redshift.

CONTEXT: Water vapour maser emission from evolved oxygen-rich stars remains poorly understood. Additional observations, including polarisation studies and simultaneous observation of different maser transitions may ultimately lead to greater insight. AIMS: We have aimed to elucidate the nature and structure of the VY CMa water vapour masers in part by observationally testing a theoretical prediction of the relative strengths of the 620.701 GHz and the 22.235 GHz maser components of ortho water vapour. METHODS: In its high-resolution mode (HRS) the Herschel Heterodyne Instrument for the Infrared (HIFI) offers a frequency resolution of 0.125 MHz, corresponding to a line-of-sight velocity of 0.06 km/s, which we employed to obtain the strength and linear polarisation of maser spikes in the spectrum of VY CMa at 620.701 GHz. Simultaneous ground based observations of the 22.235 GHz maser with the Max-Planck-Institut f\"ur Radioastronomie 100-meter telescope at Effelsberg, provided a ratio of 620.701 GHz to 22.235 GHz emission. RESULTS:We report the first astronomical detection to date of water vapour maser emission at 620.701 GHz. In VY CMa both the 620.701 and the 22.235 GHz polarisation are weak. At 620.701 GHz the maser peaks are superposed on what appears to be a broad emission component, jointly ejected asymmetrically from the star. We observed the 620.701 GHz emission at two epochs 21 days apart, both to measure the potential direction of linearly polarised maser components and to obtain a measure of the longevity of these components. Although we do not detect significant polarisation levels in the core of the line, they rise up to approximately 6% in its wings.

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.

We carried out an unbiased, spectroscopic survey using the low-resolution module of the infrared spectrograph (IRS) on board Spitzer targeting two 2.6 square arcminute regions in the GOODS-North field. IRS was used in spectral mapping mode with 5 hours of effective integration time per pixel. One region was covered between 14 and 21 microns and the other between 20 and 35 microns. We extracted spectra for 45 sources. About 84% of the sources have reported detections by GOODS at 24 microns, with a median F_nu(24um) ~ 100 uJy. All but one source are detected in all four IRAC bands, 3.6 to 8 microns. We use a new cross-correlation technique to measure redshifts and estimate IRS spectral types; this was successful for ~60% of the spectra. Fourteen sources show significant PAH emission, four mostly SiO absorption, eight present mixed spectral signatures (low PAH and/or SiO) and two show a single line in emission. For the remaining 17, no spectral features were detected. Redshifts range from z ~ 0.2 to z ~ 2.2, with a median of 1. IR Luminosities are roughly estimated from 24 microns flux densities, and have median values of 2.2 x 10^{11} L_{\odot} and 7.5 x 10^{11} L_{\odot} at z ~ 1 and z ~ 2 respectively. This sample has fewer AGN than previous faint samples observed with IRS, which we attribute to the fainter luminosities reached here.

We present mid-infrared spectral maps of the NGC 1333 star forming region, obtained with the the Infrared Spectrometer on board the Spitzer Space Telescope. Eight pure H2 rotational lines, from S (0) to S (7), are detected and mapped. The H2 emission appears to be associated with the warm gas shocked by the multiple outflows present in the region. A comparison between the observed intensities and the predictions of detailed shock models indicates that the emission arises in both slow (12 - 24 km/s) and fast (36 - 53 km/s) C-type shocks with an initial ortho-to-para ratio of ~ 1. The present H2 ortho-to-para ratio exhibits a large degree of spatial variations. In the post-shocked gas, it is usually about 2, i.e. close to the equilibrium value (~ 3). However, around at least two outflows, we observe a region with a much lower (~ 0.5) ortho-to-para ratio. This region probably corresponds to gas which has been heated-up recently by the passage of a shock front, but whose ortho-to-para has not reached equilibrium yet. This, together with the low initial ortho-to-para ratio needed to reproduce the observed emission, provide strong evidence that H2 is mostly in para form in cold molecular clouds. The H2 lines are found to contribute to 25 - 50% of the total outflow luminosity, and thus can be used to ascertain the importance of star formation feedback on the natal cloud. From these lines, we determine the outflow mass loss rate and, indirectly, the stellar infall rate, the outflow momentum and the kinetic energy injected into the cloud over the embedded phase. The latter is found to exceed the binding energy of individual cores, suggesting that outflows could be the main mechanism for core disruption.

This white paper is submitted to the Astronomy and Astrophysics 2010 Decadal Survey (Astro2010)1 Committee on the State of the Profession to emphasize the potential of the Stratospheric Observatory for Infrared Astronomy (SOFIA) to contribute to the training of instrumentalists and observers, and to related technology developments. This potential goes beyond the primary mission of SOFIA, which is to carry out unique, high priority astronomical research. SOFIA is a Boeing 747SP aircraft with a 2.5 meter telescope. It will enable astronomical observations anywhere, any time, and at most wavelengths between 0.3 microns and 1.6 mm not accessible from ground-based observatories. These attributes, accruing from the mobility and flight altitude of SOFIA, guarantee a wealth of scientific return. Its instrument teams (nine in the first generation) and guest investigators will do suborbital astronomy in a shirt-sleeve environment. The project will invest $10M per year in science instrument development over a lifetime of 20 years. This, frequent flight opportunities, and operation that enables rapid changes of science instruments and hands-on in-flight access to the instruments, assure a unique and extensive potential - both for training young instrumentalists and for encouraging and deploying nascent technologies. Novel instruments covering optical, infrared, and submillimeter bands can be developed for and tested on SOFIA by their developers (including apprentices) for their own observations and for those of guest observers, to validate technologies and maximize observational effectiveness.

When and how did galaxies form and their metals accumulate? Over the last decade, this has moved from an archeological question to a live investigation: there is now a broad picture of the evolution of galaxies in dark matter halos: their masses, stars, metals and supermassive blackholes. Galaxies have been found and studied in which these formation processes are taking place most vigorously, all the way back in cosmic time to when the intergalactic medium (IGM) was still largely neutral. However, the details of how and why the interstellar medium (ISM) in distant galaxies cools, is processed, recycled and enriched in metals by stars, and fuels active galactic nuclei (AGNs) remain uncertain. In particular, the cooling of gas to fuel star formation, and the chemistry and physics of the most intensely active regions is hidden from view at optical wavelengths, but can be seen and diagnosed at mid- & far-infrared (IR) wavelengths. Rest-frame IR observations are important first to identify the most luminous, interesting and important galaxies, secondly to quantify accurately their total luminosity, and finally to use spectroscopy to trace the conditions in the molecular and atomic gas out of which stars form. In order to map out these processes over the full range of environments and large-scale structures found in the universe - from the densest clusters of galaxies to the emptiest voids - we require tools for deep, large area surveys, of millions of galaxies out to z~5, and for detailed follow-up spectroscopy. The necessary tools can be realized technically. Here, we outline the requirements for gathering the crucial information to build, validate and challenge models of galaxy evolution.

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

By providing sensitive sub-arcsecond images and integral field spectroscopy in the 25 - 400 micron wavelength range, a far-IR interferometer will revolutionize our understanding of planetary system formation, reveal otherwise-undetectable planets through the disk perturbations they induce, and spectroscopically probe the atmospheres of extrasolar giant planets in orbits typical of most of the planets in our solar system. The technical challenges associated with interferometry in the far-IR are greatly relaxed relative to those encountered at shorter wavelengths or when starlight nulling is required. A structurally connected far-IR interferometer with a maximum baseline length of 36 m can resolve the interesting spatial structures in nascent and developed exoplanetary systems and measure exozodiacal emission at a sensitivity level critical to TPF-I mission planning. The Space Infrared Interferometric Telescope was recommended in the Community Plan for Far-IR/Submillimeter Space Astronomy, studied as a Probe-class mission, and estimated to cost 800M dollars. The scientific communities in Europe, Japan, and Canada have also demonstrated a keen interest in far-IR interferometry through mission planning workshops and technology research, suggesting the possibility of an international collaborative effort.

On 4 July 2005 at 5:52 UT the Deep Impact mission successfully completed its goal to hit the nucleus of 9P/Tempel 1 with an impactor, forming a crater on the nucleus and ejecting material into the coma of the comet. NASA's Submillimeter Wave Astronomy Satellite (SWAS) observed the 1(10)-1(01) ortho-water ground-state rotational transition in comet 9P/Tempel 1 before, during, and after the impact. No excess emission from the impact was detected by SWAS and we derive an upper limit of 1.8e7 kg on the water ice evaporated by the impact. However, the water production rate of the comet showed large natural variations of more than a factor of three during the weeks before and after the impact. Episodes of increased activity with Q(H2O)~1e28 molecule/s alternated with periods with low outgassing (Q(H2O)<~5e27 molecule/s). We estimate that 9P/Tempel 1 vaporized a total of N~4.5e34 water molecules (~1.3e9 kg) during June-September 2005. Our observations indicate that only a small fraction of the nucleus of Tempel 1 appears to be covered with active areas. Water vapor is expected to emanate predominantly from topographic features periodically facing the Sun as the comet rotates. We calculate that appreciable asymmetries of these features could lead to a spin-down or spin-up of the nucleus at observable rates.

On 4 July 2005 at 1:52 UT the Deep Impact mission successfully completed its goal to hit the nucleus of 9P/Tempel 1 with an impactor, forming a crater on the nucleus and ejecting material into the coma of the comet. The 370 kg impactor collided with the sunlit side of the nucleus with a relative velocity of 10.2 km/s. NASA's Submillimeter Wave Astronomy Satellite (SWAS) observed the 1(10)-1(01) ortho-water ground-state rotational transition in comet 9P/Tempel 1 before, during, and after the impact. No excess emission from the impact was detected by SWAS. However, the water production rate of the comet showed large natural variations of more than a factor of three during the weeks before the impact.

We report new observations of the far infrared and submillimeter water vapor emission of NML Cygnus based on data gathered with the Infrared Space Observatory and the Submillimeter Wave Astronomy Satellite. We compare the emission from NML Cyg to that previously published for VY CMa and W Hya in an attempt to establish the validity of recently proposed models for the outflow from evolved stars. The data obtained support the contention by Ivezi\'{c} & Elitzur (1995, 1997) that the atmospheres of evolved stars obey a set of scaling laws in which the optical depth of the outflow is the single most significant scaling parameter, affecting both the radiative transfer and the dynamics of the outflow. Specifically, we provide observations comparing the water vapor emission from NML Cyg, VY CMa and W Hya, and find, to the extent permitted by the quality of our data, that the results are in reasonable agreement with a model developed by Zubko & Elitzur (2000). Using this model we derive a mass loss based on the dust opacities, spectral line fluxes, and outflow velocities of water vapor observed in the atmospheres of these oxygen-rich giants. For VY CMa and NML Cyg we also obtain an estimate of the stellar mass.

Astronomical observations of the orbital angular momentum of photons, a property of electromagnetic radiation that has come to the fore in recent years, have apparently never been attempted. Here, I show that measurements of this property of photons have a number of astrophysical applications.

A prediction of standard inflationary cosmology is that the elemental composition of the medium out of which the earliest stars and galaxies condensed consisted primarily of hydrogen and helium 4He with small admixtures of deuterium, lithium 7Li, and 3He. The most red-shifted quasars, galaxies, and Ly-alpha absorbers currently observed, however, all exhibit at least some admixture of heavier elements, as do the most ancient stars in the Galaxy. Here we examine ways in which the abundance of these same elements, if present before the epoch of population III formation, might be observationally established or ruled out.

ESPI has been proposed for direct imaging and spectral analysis of giant planets orbiting solar-type stars. ESPI extends the concept suggested by Nisenson and Papaliolios (2001) for a square aperture apodized telescope that has sufficient dynamic range to directly detect exo-planets. With a 1.5 M square mirror, ESPI can deliver high dynamic range imagery as close as 0.3 arcseconds to bright sources, permitting a sensitive search for exoplanets around nearby stars and a study of their characteristics in reflected light.

The question "How did we get here and what will the future bring?" captures the human imagination and the attention of the National Academy of Science's Astronomy and Astrophysics Survey Commitee (AASC). Fulfillment of this "fundamental goal" requires astronomers to have sensitive, high angular and spectral resolution observations in the far-infrared/submillimeter (far-IR/sub-mm) spectral region. With half the luminosity of the universe and vital information about galaxy, star and planet formation, observations in this spectral region require capabilities similar to those currently available or planned at shorter wavelengths. In this paper we summarize the scientific motivation, some mission concepts and technology requirements for far-IR/sub-mm space interferometers that can be developed in the 2010-2020 timeframe.

We report the first detection of thermal water vapor emission in the 557 GHz, $1_{10} - 1_{01}$ ground state transition of ortho-H$_2$O toward VY Canis Majoris. In observations obtained with the Submillimeter Wave Astronomy Satellite (SWAS), we measured a flux of $\sim 450 $Jy, in a spectrally resolved line centered on a velocity $v_{LSR} = 25 $km s$^{-1}$ with a full width half maximum of $\sim 35$ km s$^{-1}$, somewhat dependent on the assumed line shape. We analyze the line shape in the context of three different radial outflow models for which we provide analytical expressions. We also detected a weaker 557 GHz emission line from W Hydrae. We find that these and other H$_2$O emission line strengths scale as suggested by Zubko and Elitzur (2000).

We combine spectra of VY CMa obtained with the short- and long-wavelength spectrometers, SWS and LWS, on the Infrared Space Observatory to provide a first detailed continuum spectrum of this highly luminous star. The circumstellar dust cloud through which the star is observed is partially self-absorbing, which makes for complex computational modeling. We review previous work and comment on the range of uncertainties about the physical traits and mineralogical composition of the modeled disk. We show that these uncertainties significantly affect the modeling of the outflow and the estimated mass loss. In particular, we demonstrate that a variety of quite diverse models can produce good fits to the observed spectrum. If the outflow is steady, and the radiative repulsion on the dust cloud dominates the star's gravitational attraction, we show that the total dust mass-loss rate is $\sim 4\times 10^{-6}M_{\odot}$ yr$^{-1}$, assuming that the star is at a distance of 1.5 kpc. Several indications, however, suggest that the outflow from the star may be spasmodic. We discuss this and other problems facing the construction of a physically coherent model of the dust cloud and a realistic mass-loss analysis.

Harwit, Protheroe, and Biermann (1999) recently proposed that Bose-Einstein photon bunching might significantly affect the interpretation of Cerenkov counts of TeV gamma photons. Here, we show that a combination of two recent results of Aharonian et al. (2000) and Aharonian et al. (2001) permits us to set new, more stringent upper limits of $\lesssim 10%$ on the fractional amount of photon bunching in the 7-10 TeV radiation from Markarian 501. Potential bunching at even higher energies should nevertheless continue to be investigated for this and other TeV sources, since a clear understanding of TeV energy spectra is required to unambiguously determine the spectral energy density of the mid-infrared extragalactic background.

The FIRBACK (Far Infrared BACKground) survey is one of the deepest imaging surveys carried out at 170 microns with ISOPHOT onboard ISO, and is aimed at the study of the structure of the Cosmic Far Infrared Background. This paper provides the analysis of resolved sources. After a validated process of data reduction and calibration, we perform intensive simulations to optimize the source extraction, measure the confusion noise (sigma_c = 45 mJy), and give the photometric and astrometric accuracies. 196 galaxies with flux S > 3 sigma_c are detected in the area of 3.89 square degrees. Counts of sources with flux S > 4 sigma_c present a steep slope of 3.3 +/- 0.6 on a differential "logN-logS" plot between 180 and 500 mJy. As a consequence, the confusion level is high and will impact dramatically on future IR deep surveys. This strong evolution, compared with a slope of 2.5 from Euclidian geometry, is in line with models implying a strongly evolving Luminous Infrared Galaxy population. The resolved sources account for less than 10% of the Cosmic Infrared Background at 170 microns, which is expected to be resolved into sources in the 1 to 10 mJy range.

We mapped a 13 by 22 pc region in emission from 492 GHz [C I] and, for the first time, 551 GHz ^{13}CO(5-4) in the giant molecular cloud M17SW, using the Submillimeter Wave Astronomy Satellite. The morphologies of the [C I] and ^{13}CO emission are strikingly similar. The extent and intensity of the [C I] and ^{13}CO(5-4) emission is explained as arising from photodissociation regions on the surfaces of embedded molecular clumps. Modeling of the ^{13}CO(5-4) emission in comparison to ^{13}CO(1-0) indicates a temperature gradient across the cloud, peaking to at least 63 K near the M17 ionization front and decreasing to at least 20 K at the western edge of the cloud. We see no correlation between gas density and column density. The beam-averaged column density of C I in the core is 1x10^{18} cm^-2, and the mean column density ratio N(C I)/N(CO) is about 0.4. The variations of N(C I)/N(CO) with position in M17SW indicate a similar clump size distribution throughout the cloud.

We have used the Submillimeter Wave Astronomy Satellite to map the ground-state 1_{10}-1_{01} transition of ortho-water at 557 GHz in the Orion molecular cloud. Water emission was detected in Orion over an angular extent of about 20 arcmin, or nearly 3 pc. The water emission is relatively weak, with line widths (3-6 km s^{-1}) and V_{LSR} velocities (9-11 km s^{-1}) consistent with an origin in the cold gas of the molecular ridge. We find that the ortho-water abundance relative to H_2 in the extended gas in Orion varies between 1 and 8x10^{-8}, with an average of 3x10^{-8}. The absence of detectable narrow-line ortho-H_2^{18}O emission is used to set a 3-sigma upper limit on the relative ortho-water abundance of 7x10^{-8}.

We present a 17-point map of the M17SW cloud core in the 1_{10}-1_{01} transition of ortho-water at 557 GHz obtained with the Submillimeter Wave Astronomy Satellite. Water emission was detected in 11 of the 17 observed positions. The line widths of the water emission vary between 4 and 9 km s^{-1}, and are similar to other emission lines that arise in the M17SW core. A direct comparison is made between the spatial extent of the water emission and the ^{13}CO J = 5\to4 emission; the good agreement suggests that the water emission arises in the same warm, dense gas as the ^{13}CO emission. A spectrum of the H_2^{18}O line was also obtained at the center position of the cloud core, but no emission was detected. We estimate that the average abundance of ortho-water relative to H_2 within the M17 dense core is approximately 1x10^{-9}, 30 times smaller than the average for the Orion core. Toward the H II region/molecular cloud interface in M17SW the ortho-water abundance may be about 5 times larger than in the dense core.

We present Submillimeter Wave Astronomy Satellite (SWAS) observations of the 1_{10}-1_{01} transition of ortho-water at 557 GHz toward 12 molecular cloud cores. The water emission was detected in NGC 7538, Rho Oph A, NGC 2024, CRL 2591, W3, W3(OH), Mon R2, and W33, and was not detected in TMC-1, L134N, and B335. We also present a small map of the water emission in S140. Observations of the H_2^{18}O line were obtained toward S140 and NGC 7538, but no emission was detected. The abundance of ortho-water relative to H_2 in the giant molecular cloud cores was found to vary between 6x10^{-10} and 1x10^{-8}. Five of the cloud cores in our sample have previous water detections; however, in all cases the emission is thought to arise from hot cores with small angular extents. The water abundance estimated for the hot core gas is at least 100 times larger than in the gas probed by SWAS. The most stringent upper limit on the ortho-water abundance in dark clouds is provided in TMC-1, where the 3-sigma upper limit on the ortho-water fractional abundance is 7x10^{-8}.

FIRBACK is a one of the deepest surveys performed at 170 microns with ISOPHOT onboard ISO, and is aimed at the study of cosmic far infrared background sources. About 300 galaxies are detected in an area of four square degrees, and source counts present a strong slope of 2.2 on an integral "logN-logS" plot, which cannot be due to cosmological evolution if no K-correction is present. The resolved sources account for less than 10% of the Cosmic Infrared Background at 170 microns. In order to understand the nature of the sources contributing to the CIB, and to explain deep source counts at other wavelengths, we have developed a phenomenological model, which constrains in a simple way the luminosity function evolution with redshift, and fits all the existing deep source counts from the mid-infrared to the submillimetre range. Images, materials and papers available on the FIRBACK web: http://wwwfirback.ias.u-psud.fr wwwfirback.ias.u-psud.fr

A Cosmic Far-InfraRed Background (CFIRB) has long been predicted that would traces the intial phases of galaxy formation. It has been first detected by Puget et al.(1996) using COBE data and has been later confirmed by several recent studies (Fixsen et al. 1998, Hauser et al. 1998, Lagache et al. 1999). We will present a new determination of the CFIRB that uses for the first time, in addition to COBE data, two independent gas tracers: the HI survey of Leiden/Dwingeloo (hartmann, 1998) and the WHAM H$_{\alpha}$ survey (Reynolds et al 1998). We will see that the CFIRB above 100 micron is now very well constrained. The next step is to see if we can detect its fluctuations. To search for the CFIRB fluctuations, we have used the FIRBACK observations. FIRBACK is a deep cosmological survey conducted at 170 micron with ISOPHOT (Dole et al., 2000). We show that the emission of unresolved extra-galactic sources clearly dominates, at arcminute scales, the background fluctuations in the lowest galactic emission regions. This is the first detection of the CFIRB fluctuations.

We present the results of the five mid-IR 15 microns (12-18 microns LW3 band) ISOCAM Guaranteed Time Extragalactic Surveys performed in the regions of the Lockman Hole and Marano Field. The roughly 1000 sources detected, 600 of which have a flux above the 80 % completeness limit, guarantee a very high statistical significance for the integral and differential source counts from 0.1 mJy up to 5 mJy. By adding the ISOCAM surveys of the HDF-North and South (plus flanking fields) and the lensing cluster A2390 at low fluxes and IRAS at high fluxes, we cover four decades in flux from 50 microJy to 0.3 Jy. The slope of the differential counts is very steep (alpha =-3.0) in the flux range 0.4-4 mJy, hence much above the Euclidean expectation of alpha =-2.5. When compared with no-evolution models based on IRAS, our counts show a factor of 10 excess at 400 microJy, and a fast convergence, with alpha =-1.6 at lower fluxes.

The recent detection of gamma radiation from Mkn 501 at energies as high as 25 TeV suggests stringent upper bounds on the diffuse, far infrared, extragalactic radiation density. The production of electron-positron pairs through photon-photon collisions would prevent gamma photons of substantially higher energies from reaching us across distances of order 100 Mpc. However, coherently arriving TeV or sub-TeV gammas - Bose-Einstein condensations of photons at these energies - could mimic the Cherenkov shower signatures of extremely energetic gammas. To better understand such events, we describe their observational traits and discuss how they might be generated.

We report the detection of numerous far-infrared emission lines of water vapor toward the supergiant star VY Canis Majoris. A 29.5 - 45 micron grating scan of VY CMa, obtained using the Short Wavelength Spectrometer (SWS) of the Infrared Space Observatory (ISO) at a spectral resolving power of approximately 2000, reveals at least 41 spectral features due to water vapor that together radiate a total luminosity ~ 25 solar luminosities. In addition to pure rotational transitions within the ground vibrational state, these features include rotational transitions within the (010) excited vibrational state. The spectrum also shows the doublet Pi 1/2 (J=5/2) <-- doublet Pi 3/2 (J=3/2) OH feature near 34.6 micron in absorption. Additional SWS observations of VY CMa were carried out in the instrument's Fabry-Perot mode for three water transitions: the 7(25)-6(16) line at 29.8367 micron, the 4(41)-3(12) line 31.7721 micron, and the 4(32)-3(03) line at 40.6909 micron. The higher spectral resolving power of approximately 30,000 thereby obtained permits the line profiles to be resolved spectrally for the first time and reveals the "P Cygni" profiles that are characteristic of emission from an outflowing envelope.

FIRBACK is one of the deepest cosmological surveys performed in the far infrared, using ISOPHOT. We describe this survey, its data reduction and analysis. We present the maps of fields at 175 microns. We point out some first results: source identifications with radio and mid infrared, and source counts at 175 microns. These two results suggest that half of the FIRBACK sources are probably at redshifts greater than 1. We also present briefly the large follow-up program.

ISOCAM extragalactic mid-infrared deep surveys have detected a population of strong IR emitters ten times more numerous than expected if there were no evolution from z<0.2 (IRAS) up to the maximum redshift of these galaxies (z=1.5). The mid-IR cosmic background produced by these galaxies (2.35+/-0.8 nW m^-2 sr^-1, at 15 microns above 50 microJy) is larger than 30 per cent of the energy radiated in the I band by the optical galaxies detected in the HDF, which are two orders of magnitude more numerous. This fraction is much higher than in the local universe (z<0.2) as probed by IRAS, where all integrated IR emission from 8 to 1000 microns makes 30 per cent of the optical starlight. Even assuming a conservative SED, they produce a major contribution to the 140 mu m DIRBE background measured by Hauser et al. (1998, 25.1+/-7 nW m^-2 sr^-1) and Lagache et al. (1999, 15.3+/-9.5 nW m^ -2 sr^-1). This is both a confirmation of the strong IR cosmic background that was first detected by Puget et al. (1996) and the first identification of the galaxies responsible for a large fraction of this background at 140 microns. We find that these bright (several 10^11 Lsol) and massive (<M>= 1.5x10^11 Msol) IR galaxies have a typical redshift of z=0.7 and optical colors similar to field galaxies. The fraction of galaxies with morphological signs of interactions is larger at higher redshift. This is consistent with the nature of local luminous IR galaxies found by IRAS, which often exhibit evidence of galaxy interaction or merging and indicates that interactions should play a crucial role in the evolution of galaxies as expected from bottom-up scenarios.

A major goal of modern astrophysics is to understand the processes by which the universe evolved from its initial simplicity, as seen in measurements of the Cosmic Microwave Background, to the universe we see today, with complexity on all scales. While the diffuse background measurements of COBE reveal the importance of the far infrared and submillimeter in early galaxy and star formation, the understanding of the development of complex structure requires high resolution imaging and spectroscopy. We present a concept for a space mission called SPECS, the Submillimeter Probe of the Evolution of Cosmic Structure, which provides these capabilities. SPECS is a cold, spatial and spectral Michelson interferometer with adjustable baselines ranging up to 1 km. It has Hubble sensitivity and angular resolution in the far IR and submillimeter, spectral resolution up to 10,000, and a 14' field of view. SPECS will be able to image thermal dust continuum and infrared cooling and diagnostic line emission over a wide range of redshifts, providing extinction-free astrophysical probes of young galaxies and early cosmic structures, and measures of the luminosity and heavy element formation history of the universe. SPECS would also have the potential to improve vastly our knowledge of protostars, protoplanetary systems, Active Galactic Nuclei and other objects in the local universe. We recommend that a concerted effort be made during the next decade to develop and test certain critical technologies (photon-counting far IR detectors, formation flying spacecraft, cold, lightweight mirrors, and active coolers), so that SPECS can be deployed in about 2015.

The density of the far infrared / submillimeter (FIR/SMM) diffuse extragalactic radiation field has recently been determined from COBE data. Nearly simultaneously, deep FIR/SMM surveys have detected substantial numbers of optically unidentified sources, which have led to the proposal that galaxies and protogalaxies at red shifts z = 2 to 4 may account for an appreciable fraction of the background. Here, I show that, if the reported radiation levels are generated through nucleosynthesis, most of this energy must have been produced at epochs $z \lesssim 2$. Hubble Deep Field data cited by Madau et al. (1998) indicate that the bulk of the integrated extragalactic background must have been generated even more recently at $z < 1$.

We have detected the S(1), S(2), S(3), S(4) and S(5) pure rotational lines of molecular hydrogen toward the outflow source HH 54, using the Short Wavelength Spectrometer (SWS) on board the Infrared Space Observatory (ISO). The observed H2 line ratios indicate the presence of warm molecular gas with an H2 density of at least 100,000 per cubic centimeter and a temperature ~ 650 K in which the ortho- to para-H2 ratio is only 1.2 +/- 0.4, significantly smaller than the equilibrium ratio of 3 expected in gas at that temperature. These observations imply that the measured ortho- to para-H2 ratio is the legacy of an earlier stage in the thermal history of the gas when the gas had reached equilibrium at a temperature < 90 K. Based upon the expected timescale for equilibration, we argue that the non-equilibrium ortho- to para-H2 ratio observed in HH 54 serves as a chronometer that places a conservative upper limit ~ 5000 yr on the period for which the emitting gas has been warm. The S(2)/S(1) and S(3)/S(1) H2 line ratios measured toward HH 54 are consistent with recent theoretical models of Timmermann for the conversion of para- to ortho-H2 behind slow, `C'-type shocks, but only if the preshock ortho- to para-H2 ratio was < 0.2.

Using the Long Wavelength Spectrometer (LWS) onboard the Infrared Space Observatory (ISO), we have observed thermal water vapor emission from a roughly circular field of view approximately 75 arc seconds in diameter centered on the Orion BN-KL region. The Fabry-Perot line strengths, line widths, and spectral line shifts observed in eight transitions between 71 and 125 microns show good agreement with models of thermal emission arising from a molecular cloud subjected to a magnetohydrodynamic C-type shock. Both the breadth and the relative strengths of the observed lines argue for emission from a shock rather than from warm quiescent gas in the Orion core. Though one of the eight transitions appears anomalously strong, and may be subject to the effects of radiative pumping, the other seven indicate an H2O/H2 abundance ratio of order 5E-4, and a corresponding gas-phase oxygen-to-hydrogen abundance ratio of order 4E-4. Given current estimates of the interstellar, gas-phase, oxygen and carbon abundances in the solar vicinity, this value is consistent with theoretical shock models that predict the conversion into water of all the gas-phase oxygen that is not bound as CO. The overall cooling provided by rotational transitions of H2O in this region appears to be comparable to the cooling through rotational lines of CO, but is an order of magnitude lower than cooling through H2 emission. However, the model that best fits our observations shows cooling by H2O and CO dominant in that portion of the post-shock region where temperatures are below ~ 800 K and neither vibrational nor rotational radiative cooling by H2 is appreciable.