Cold brown dwarfs are excellent analogs of widely separated, gas giant exoplanets, and provide insight into the potential atmospheric chemistry and physics we may encounter in objects discovered by future direct imaging surveys. We present a low resolution R $\sim$ 300 $M$-band spectroscopic sequence of seven brown dwarfs with effective temperatures between 750 K and 250 K along with Jupiter. These spectra reveal disequilibrium abundances of carbon monoxide (CO) produced by atmospheric quenching. We use the eddy diffusion coefficient (K$_{zz}$) to estimate the strength of vertical mixing in each object. The K$_{zz}$ values of cooler gaseous objects are close to their theoretical maximum and warmer objects show weaker mixing, likely due to less efficient convective mixing in primarily radiative layers. The CO-derived K$_{zz}$ values imply that disequilibrium phosphine (PH$_{3}$) should be easily observable in all of the brown dwarfs, but none as yet show any evidence for PH$_{3}$ absorption. We find that ammonia is relatively insensitive to atmospheric quenching at these effective temperatures. We are able to improve the fit to WISE 0855's $M$-band spectrum by including both CO and water clouds in the atmospheric model.

The origin of life on Earth seems to demand a highly reduced early atmosphere, rich in CH4, H2, and NH3, but geological evidence suggests that Earth's mantle has always been relatively oxidized and its emissions dominated by CO2 H2O, and N2. The paradox can be resolved by exploiting the reducing power inherent in the "late veneer," i.e., material accreted by Earth after the Moon-forming impact. Isotopic evidence indicates that the late veneer consisted of extremely dry, highly reduced inner solar system materials, suggesting that Earth's oceans were already present when the late veneer came. The major primary product of reaction between the late veneer's iron and Earth's water was H2. Ocean vaporizing impacts generate high pressures and long cooling times that favor CH4 and NH3. Impacts too small to vaporize the oceans are much less productive of CH4 and NH3, unless (i) catalysts were available to speed their formation, or (ii) additional reducing power was extracted from pre-existing crustal or mantle materials. The transient H2-CH4 atmospheres evolve photochemically to generate nitrogenated hydrocarbons at rates determined by solar radiation and hydrogen escape, on timescales ranging up to tens of millions of years and with cumulative organic production ranging up to half a kilometer. Roughly one ocean of hydrogen escapes. The atmosphere after the methane's gone is typically H2 and CO rich, with eventual oxidation to CO2 rate-limited by water photolysis and hydrogen escape.

The characterization of rocky, Earth-like planets is an important goal for future large ground- and space-based telescopes. In support of developing an efficient observational strategy, we have applied Bayesian statistical inference to interpret the albedo spectrum of cloudy true-Earth analogs that include a diverse spread in their atmospheric water vapor mixing ratios. We focus on detecting water-bearing worlds by characterizing their atmospheric water vapor content via the strong 0.94$\,\mu$m H$_2$O absorption feature, with several observational configurations. Water vapor is an essential signpost when assessing planetary habitability, and determining its presence is important in vetting whether planets are suitable for hosting life. We find that R=140 spectroscopy of the absorption feature combined with a same-phase green optical photometric point at $0.525-0.575\,\mu$m is capable of distinguishing worlds with less than $0.1\times$ Earth-like water vapor levels from worlds with $1\times$ Earth-like levels or greater at a signal-to-noise ratio of 5 or better with $2\sigma$ confidence. This configuration can differentiate between $0.01\times$ and $0.1\times$ Earth-like levels when the signal-to-noise ratio is 10 or better at the same confidence. However, strong constraints on the water vapor mixing ratio remained elusive with this configuration even at signal-to-noise of 15. We find that adding the same-phase optical photometric point does not significantly help characterize the H$_2$O mixing ratio, but does enable an upper limit on atmospheric ozone levels. Finally, we find that a 0.94$\,\mu$m photometric point, instead of spectroscopy, combined with the green-optical point, fails to produce meaningful information about atmospheric water content.

On hot Jupiter exoplanets, strong horizontal and vertical winds should homogenize the abundances of the important absorbers CH$_4$ and CO much faster than chemical reactions restore chemical equilibrium. This effect, typically neglected in general circulation models (GCMs), has been suggested as explanation for discrepancies between observed infrared lightcurves and those predicted by GCMs: On the nightsides of several hot Jupiters, GCMs predict outgoing fluxes that are too large, especially in the Spitzer 4.5 $\mu$m band. We modified the SPARC/MITgcm to include disequilibrium abundances of CH$_4$, CO and H$_2$O by assuming that the CH$_4$/CO ratio is constant throughout the simulation domain. We ran simulations of hot Jupiter HD 189733b with 8 CH$_4$/CO ratios. In the more likely CO-dominated regime, we find temperature changes $\geq$50-100 K compared to the equilibrium chemistry case across large regions. This effect is large enough to affect predicted emission spectra and should thus be included in GCMs of hot Jupiters with equilibrium temperatures between 600K and 1300K. We find that spectra in regions with strong methane absorption, including the Spitzer 3.6 and 8 $\mu$m bands, are strongly impacted by disequilibrium abundances. We expect chemical quenching to result in much larger nightside fluxes in the 3.6 $\mu$m band, in stark contrast to observations. Meanwhile, we find almost no effect on predicted observations in the 4.5 $\mu$m band, as the opacity changes due to CO and H$_2$O offset each other. We thus conclude that disequilibrium carbon chemistry cannot explain the observed low nightside fluxes in the 4.5 $\mu$m band.

We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as critical inputs to atmospheric physical and chemical tools. Here we provide an up-to-date and condensed description of areas where laboratory and/or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry, building off a larger 2016 white paper, and endorsed by the NAS Exoplanet Science Strategy report. Now is the ideal time for progress in these areas, but this progress requires better access to, understanding of, and training in the production of spectroscopic data as well as a better insight into chemical reaction kinetics both thermal and radiation-induced at a broad range of temperatures. Given that most published efforts have emphasized relatively Earth-like conditions, we can expect significant and enlightening discoveries as emphasis moves to the exotic atmospheres of exoplanets.

We highlight a physical effect that is often not considered that impacts the calculation of model spectra of planets at secondary eclipse, affecting both emission and reflection spectra. The radius of the emitting surface of the planet is not merely one value measured from a transit light curve, but is itself a function of wavelength, yet it is not directly measurable. At high precision, a similar effect is well-known in transit "transmission spectroscopy" but this related effect also impacts emission and reflection. As is well-appreciated, the photospheric radius can vary across $\sim$4-8 atmospheric scale heights, depending on atmospheric opacity and spectral resolution. This effect leads to a decreased weighting in model calculations at wavelengths where atmospheric opacity is low, and one sees more deeply into the atmosphere, to a smaller radius. The overall effect serves to mute emission spectra features for atmospheres with no thermal inversion but to enhance features for atmospheres with a thermal inversion. While this effect can be ignored for current \emph{Hubble} observations, it can lead to wavelength-dependent 10-20\% changes in planet-to-star flux ratios in the infrared at $R\sim~200-1000$ (readily achievable for JWST) for low-gravity hot Jupiters, although values of 5\% are more typical for the population. The effect is mostly controlled by the ratio of the atmospheric scale height to the planet radius, and can be important at any planetary temperature. Of known planets, the effect is largest for the cool "super-puffs" at very low surface gravity, where it can alter calculated flux ratios by over 100\%. We discuss complexities of including this photospheric radius effect in 1D and 3D atmosphere models.

Planned astronomical observatories of the 2020s will be capable of obtaining reflected light photometry and spectroscopy of cool extrasolar giant planets. Here we explain that such data are valuable both for understanding the origin and evolution of giant planets as a whole and for preparing for the interpretation of similar datasets from potentially habitable extrasolar terrestrial planets in the decades to follow.

The best-characterized exoplanets to date are planets on close-in transiting orbits around their host stars. The high level of irradiation and transiting geometry of these objects make them ideal targets for atmospheric investigations. However, the modest apertures of many current telescopes allow mostly low resolution spectra to be observed for transiting planets, failing to extract key physical and chemical properties of their atmospheres. Ground-based 30-meter class telescopes will set the stage for a substantial leap in our understanding of exoplanet atmospheres. We outline a two-pronged survey, recently submitted as a Key Science Project (KSP) for the US ELTs, which would yield unprecedented insight into the atmospheres of close-in exoplanets via combined observations with the GMT and TMT. (1) The first opportunity involves measuring the global-scale atmospheric circulation and planetary rotation for a sample of 40 hot Jupiters to glean insight into the unique radiative forcing regime governing highly-irradiated, tidally-locked giant planets. (2) The second opportunity involves measuring atmospheric mass-loss and extracting atmospheric composition and abundance ratio information for $\sim60$ sub-Neptunes and super-Earths (including candidate disintegrating planets) to constrain their formation and evolution histories. These efforts would be made possible by the unparalleled combination of high spectral resolution instrumentation and large aperture size of the ELTs. This survey would enable the first statistical study of atmospheric circulation in extrasolar giant planets, and would provide detections of trace gases and measurements of atmospheric escape in small-planet atmospheres, far exceeding the reach of \textit{JWST}.

Surveys reveal that terrestrial- to Neptune-sized planets (1 $< R <$ 4 R$_{\rm{Earth}}$) are the most common type of planets in our galaxy. Detecting and characterizing such small planets around nearby stars holds the key to understanding the diversity of exoplanets and will ultimately address the ubiquitousness of life in the universe. The following fundamental questions will drive research in the next decade and beyond: (1) how common are terrestrial to Neptune-sized planets within a few AU of their host star, as a function of stellar mass? (2) How does planet composition depend on planet mass, orbital radius, and host star properties? (3) What are the energy budgets, atmospheric dynamics, and climates of the nearest worlds? Addressing these questions requires: a) diffraction-limited spatial resolution; b) stability and achievable contrast delivered by adaptive optics; and c) the light-gathering power of extremely large telescopes (ELTs), as well as multi-wavelength observations and all-sky coverage enabled by a comprehensive US ELT Program. Here we provide an overview of the challenge, and promise of success, in detecting and comprehensively characterizing small worlds around the very nearest stars to the Sun with ELTs. This white paper extends and complements the material presented in the findings and recommendations published in the National Academy reports on Exoplanet Science Strategy and Astrobiology Strategy for the Search for Life in the Universe.

A new class of exoplanets has emerged: the ultra hot Jupiters, the hottest close-in gas giants. Most of them have weaker than expected spectral features in the $1.1-1.7\mu m$ bandpass probed by HST/WFC3 but stronger spectral features at longer wavelengths probed by Spitzer. This led previous authors to puzzling conclusions about the thermal structures and chemical abundances of these planets. Using the SPARC/MITgcm, we investigate how thermal dissociation, ionization, H$^-$ opacity and clouds shape the thermal structures and spectral properties of ultra hot Jupiters with a special focus on WASP-121b. We expand our findings to the whole population of ultra hot Jupiters through analytical quantification of the thermal dissociation and its influence on the strength of spectral features. We predict that most molecules are thermally dissociated and alkalies are ionized in the dayside photospheres of ultra hot Jupiters. This includes H$_{\rm 2}$O, TiO, VO, and H$_{\rm 2}$ but not CO, which has a stronger molecular bond. The vertical molecular gradient created by the dissociation significantly weakens the spectral features from water while the $4.5\mu m$ CO feature remain unchanged. The water band in the HST/WFC3 bandpass is further weakened by H$^-$ continuum opacity. Molecules are expected to recombine before reaching the limb, leading to order of magnitude variations of the chemical composition and cloud coverage between the limb and the dayside. Overall, molecular dissociation provides a qualitative understanding of the lack of strong spectral feature of water in the $1-2\mu m$ bandpass observed in most ultra hot Jupiters. Quantitatively, however, our model does not provide a satisfactory match to the WASP-121b emission spectrum. Together with WASP-33b and Kepler-33Ab, they seem the outliers among the population of ultra hot Jupiters in need of a more thorough understanding.

The coldest brown dwarf, WISE 0855, is the closest known planetary-mass, free-floating object and has a temperature nearly as cold as the solar system gas giants. Like Jupiter, it is predicted to have an atmosphere rich in methane, water, and ammonia, with clouds of volatile ices. WISE 0855 is faint at near-infrared wavelengths and emits almost all its energy in the mid-infrared. Skemer et al. 2016 presented a spectrum of WISE 0855 from 4.5-5.1 micron (M band), revealing water vapor features. Here, we present a spectrum of WISE 0855 in L band, from 3.4-4.14 micron. We present a set of atmosphere models that include a range of compositions (metallicities and C/O ratios) and water ice clouds. Methane absorption is clearly present in the spectrum. The mid-infrared color can be better matched with a methane abundance that is depleted relative to solar abundance. We find that there is evidence for water ice clouds in the M band spectrum, and we find a lack of phosphine spectral features in both the L and M band spectra. We suggest that a deep continuum opacity source may be obscuring the near-infrared flux, possibly a deep phosphorous-bearing cloud, ammonium dihyrogen phosphate. Observations of WISE 0855 provide critical constraints for cold planetary atmospheres, bridging the temperature range between the long-studied solar system planets and accessible exoplanets. JWST will soon revolutionize our understanding of cold brown dwarfs with high-precision spectroscopy across the infrared, allowing us to study their compositions and cloud properties, and to infer their atmospheric dynamics and formation processes.

Space-based high contrast imaging mission concepts for studying rocky exoplanets in reflected light are currently under community study. We develop an inverse modeling framework to estimate the science return of such missions given different instrument design considerations. By combining an exoplanet albedo model, an instrument noise model, and an ensemble Markov chain Monte Carlo sampler, we explore retrievals of atmospheric and planetary properties for Earth twins as a function of signal-to-noise ratio (SNR) and resolution ($R$). Our forward model includes Rayleigh scattering, single-layer water clouds with patchy coverage, and pressure-dependent absorption due to water vapor, oxygen, and ozone. We simulate data at $R = 70$ and $R = 140$ from 0.4-1.0 $\mu$m with SNR $ = 5, 10, 15, 20$ at 550 nm (i.e., for HabEx/LUVOIR-type instruments). At these same SNR, we simulate data for WFIRST paired with a starshade, which includes two photometric points between 0.48-0.6 $\mu$m and $R = 50$ spectroscopy from 0.6-0.97 $\mu$m. Given our noise model for WFIRST-type detectors, we find that weak detections of water vapor, ozone, and oxygen can be achieved with observations with at least $R = 70$ / SNR$\ = 15$, or $R = 140$ / SNR$\ = 10$ for improved detections. Meaningful constraints are only achieved with $R = 140$ / SNR$\ = 20$ data. The WFIRST data offer limited diagnostic information, needing at least SNR = 20 to weakly detect gases. Most scenarios place limits on planetary radius, but cannot constrain surface gravity and, thus, planetary mass.

Exoplanets orbiting M dwarf stars are a prime target in the search for life in the Universe. M dwarf stars are active, with powerful flares that could adversely impact prospects for life, though there are counter-arguments. Here, we turn flaring to advantage and describe ways in which it can be used to enhance the detectability of planets, in the absence of transits or a coronagraph, significantly expanding the accessible discovery and characterization space. Flares produce brief bursts of intense luminosity, after which the star dims. Due to the light travel time between the star and planet, the planet receives the high intensity pulse, which it re-emits through scattering (a light echo) or intrinsic emission when the star is much fainter, thereby increasing the planet's detectability. The planet's light echo emission can potentially be discriminated from that of the host star by means of a time delay, Doppler shift, spatial shift, and polarization, each of which can improve the contrast of the planet to the star. Scattered light can reveal the albedo spectrum of the planet to within a size scale factor, and is likely to be polarized. Intrinsic emission mechanisms include fluorescent pumping of multiple molecular hydrogen and neutral oxygen lines by intense LyAlpha and LyBeta flare emission, recombination radiation of ionized and photodissociated species, and atmospheric processes such as terrestrial upper atmosphere airglow and near infrared hydroxyl emission. We discuss the feasibility of detecting light echoes and find that under favorable circumstances, echo detection is possible.

Infrared radiation emitted from a planet contains information about the chemical composition and vertical temperature profile of its atmosphere. If upper layers are cooler than lower layers, molecular gases will produce absorption features in the planetary thermal spectrum. Conversely, if there is a stratosphere - where temperature increases with altitude - these molecular features will be observed in emission. It has been suggested that stratospheres could form in highly irradiated exoplanets, but the extent to which this occurs is unresolved both theoretically and observationally. A previous claim for the presence of a stratosphere remains open to question, owing to the challenges posed by the highly variable host star and the low spectral resolution of the measurements. Here we report a near-infrared thermal spectrum for the ultrahot gas giant WASP-121b, which has an equilibrium temperature of approximately 2,500 kelvin. Water is resolved in emission, providing a detection of an exoplanet stratosphere at 5-sigma confidence. These observations imply that a substantial fraction of incident stellar radiation is retained at high altitudes in the atmosphere, possibly by absorbing chemical species such as gaseous vanadium oxide and titanium oxide.

Brown dwarf spectra are rich in information revealing of the chemical and physical processes operating in their atmospheres. We apply a recently developed atmospheric retrieval tool to an ensemble of late T-dwarf (600-800K) near infrared spectra. With these spectra we are able to place direct constraints the molecular abundances of H$_2$O, CH$_4$, CO, CO$_2$, NH$_3$, H$_2$S, and Na+K, gravity, thermal structure (and effective temperature), photometric radius, and cloud optical depths. We find that ammonia, water, methane, and the alkali metals are present and well constrained in all 11 objects. From the abundance constraints we find no significant trend in the water, methane, or ammonia abundances with temperature, but find a very strong ($>$25$\sigma$) increasing trend in the alkali metal abundances with effective temperature, indicative of alkali rainout. We also find little evidence for optically thick clouds. With the methane and water abundances, we derive the intrinsic atmospheric metallicity and carbon-to-oxygen ratios. We find in our sample, that metallicities are typically sub solar and carbon-to-oxygen ratios are somewhat super solar, different than expectations from the local stellar population. We also find that the retrieved vertical thermal profiles are consistent with radiative equilibrium over the photospheric regions. Finally, we find that our retrieved effective temperatures are lower than previous inferences for some objects and that our radii are larger than expectations from evolutionary models, possibly indicative of un-resolved binaries. This investigation and methodology represents a paradigm in linking spectra to the determination of the fundamental chemical and physical processes governing cool brown dwarf atmospheres.

We present the first results from applying the spectral inversion technique in the cloudy L dwarf regime. Our new framework provides a flexible approach to modelling cloud opacity which can be built incrementally as the data requires, and improves upon previous retrieval experiments in the brown dwarf regime by allowing for scattering in two stream radiative transfer. Our first application of the tool to two mid-L dwarfs is able to reproduce their near-infrared spectra far more closely than grid models. Our retrieved thermal, chemical, and cloud profiles allow us to estimate $T_{\rm eff} = 1796^{+23}_{-25}$ K and $\log g = 5.21^{+0.05}_{-0.08}$ for 2MASS J05002100+0330501 and for 2MASSW J2224438-015852 we find $T_{\rm eff} = 1723^{+18}_{-19}$ K and $\log g = 5.31^{+0.04}_{-0.08}$, in close agreement with previous empirical estimates. Our best model for both objects includes an optically thick cloud deck which passes $\tau_{cloud} \geq 1$ (looking down) at a pressure of around 5 bar. The temperature at this pressure is too high for silicate species to condense, and we argue that corundum and/or iron clouds are responsible for this cloud opacity. Our retrieved profiles are cooler at depth, and warmer at altitude than the forward grid models that we compare, and we argue that some form of heating mechanism may be at work in the upper atmospheres of these L dwarfs. We also identify anomalously high CO abundance in both targets, which does not correlate with the warmth of our upper atmospheres or our choice of cloud model, and find similarly anomalous alkali abundance for one of our targets. These anomalies may reflect unrecognised shortcomings in our retrieval model, or inaccuracies in our gas phase opacities.

A number of exciting concepts are under development for Flagship, Probe class, Explorer class, and Suborbital class NASA missions in the ultraviolet/optical spectral ranges. These missions will depend on high performance silicon detector arrays being delivered affordably and in high numbers. In a focused effort we have advanced delta-doping technology to high throughput and high yield wafer-scale processing, encompassing a multitude of state-of-the-art silicon-based detector formats and designs. As part of this technology advancement and in preparation for upcoming missions, we have embarked on a number of field observations, instrument integrations, and independent evaluations of delta-doped arrays. In this paper, we present recent data and innovations from the Advanced Detectors and Systems program at JPL, including two-dimensional doping technology; our end-to-end post-fabrication processing of high performance UV/Optical/NIR arrays; and advanced coatings for detectors and optical elements. Additionally, we present examples of past, in-progress, and planned observations and deployments of delta-doped arrays.

Future space-based telescopes, such as the Wide-Field Infrared Survey Telescope (WFIRST), will observe the reflected-light spectra of directly imaged extrasolar planets. Interpretation of such data presents a number of novel challenges, including accounting for unknown planet radius and uncertain stellar illumination phase angle. Here we report on our continued development of Markov Chain Monte Carlo retrieval methods for addressing these issues in the interpretation of such data. Specifically we explore how the unknown planet radius and potentially poorly known observer-planet-star phase angle impacts retrievals of parameters of interest such as atmospheric methane abundance, cloud properties and surface gravity. As expected, the uncertainty in retrieved values is a strong function of signal-to-noise ratio (SNR) of the observed spectra, particularly for low metallicity atmospheres, which lack deep absorption signatures. Meaningful results may only be possible above certain SNR thresholds; for cases across a metallicity range of 1-50 times solar, we find that only an SNR of 20 systematically reproduces close to the correct methane abundance at all phase angles. However, even in cases where the phase angle is poorly known we find that the planet radius can be constrained to within a factor of two. We find that uncertainty in planet radius decreases at phase angles past quadrature, as the highly forward scattering nature of the atmosphere at these geometries limits the possible volume of phase space that relevant parameters can occupy. Finally, we present an estimation of possible improvement that can result from combining retrievals against observations at multiple phase angles.

We present Exo-Transmit, a software package to calculate exoplanet transmission spectra for planets of varied composition. The code is designed to generate spectra of planets with a wide range of atmospheric composition, temperature, surface gravity, and size, and is therefore applicable to exoplanets ranging in mass and size from hot Jupiters down to rocky super-Earths. Spectra can be generated with or without clouds or hazes with options to (1) include an optically thick cloud deck at a user-specified atmospheric pressure or (2) to augment the nominal Rayleigh scattering by a user-specified factor. The Exo-Transmit code is written in C and is extremely easy to use. Typically the user will only need to edit parameters in a single user input file in order to run the code for a planet of their choosing. Exo-Transmit is available publicly on Github with open-source licensing at https://github.com/elizakempton/Exo_Transmit .

The Neptune-mass GJ 436b is one of the most-studied transiting exoplanets with repeated measurements of both its thermal emission and transmission spectra. We build on previous studies to answer outstanding questions about this planet, including its potentially high metallicity and tidal heating of its interior. We present new observations of GJ 436b's thermal emission at 3.6 and 4.5 micron, which reduce uncertainties in estimates of GJ 436b's flux at those wavelengths and demonstrate consistency between Spitzer observations spanning more than 7 years. We analyze the Spitzer thermal emission photometry and Hubble WFC3 transmission spectrum in tandem. We use a powerful dual-pronged modeling approach, comparing these data to both self-consistent and retrieval models. We vary the metallicity, intrinsic luminosity from tidal heating, disequilibrium chemistry, and heat redistribution. We also study the effect of clouds and photochemical hazes on the spectra, but do not find strong evidence for either. The self-consistent and retrieval modeling combine to suggest that GJ 436b has a high atmospheric metallicity, with best fits at or above several hundred times solar metallicity, tidal heating warming its interior with best-fit intrinsic effective effective temperatures around 300--350 K, and disequilibrium chemistry. High metal-enrichments (>600x solar) can only occur from the accretion of rocky, rather than icy, material. Assuming Tint~300--350 K, we find that Q'~2x10^5--10^6, larger than Neptune's Q', and implying a long tidal circularization timescale for the planet's orbit. We suggest that Neptune-mass planets may be a more diverse class than previously imagined, with metal-enhancements potentially spanning several orders of magnitude, to perhaps over 1000x solar metallicity. High fidelity observations with instruments like JWST will be critical for characterizing this diversity.

Upcoming space-based coronagraphic instruments in the next decade will perform reflected light spectroscopy and photometry of cool, directly imaged extrasolar giant planets. We are developing a new atmospheric retrieval methodology to help assess the science return and inform the instrument design for such future missions, and ultimately interpret the resulting observations. Our retrieval technique employs a geometric albedo model coupled with both a Markov chain Monte Carlo Ensemble Sampler (emcee) and a multimodal nested sampling algorithm (MultiNest) to map the posterior distribution. This combination makes the global evidence calculation more robust for any given model, and highlights possible discrepancies in the likelihood maps. As a proof-of-concept, our current atmospheric model contains 1 or 2 cloud layers, methane as a major absorber, and a H$_2$-He background gas. This 6-to-9 parameter model is appropriate for Jupiter-like planets and can be easily expanded in the future. In addition to deriving the marginal likelihood distribution and confidence intervals for the model parameters, we perform model selection to determine the significance of methane and cloud detection as a function of expected signal-to-noise in the presence of spectral noise correlations. After internal validation, the method is applied to realistic spectra of Jupiter, Saturn, and HD 99492 c, a model observing target. We find that the presence or absence of clouds and methane can be determined with high confidence, while parameter uncertainties are model-dependent and correlated. Such general methods will also be applicable to the interpretation of direct imaging spectra of cloudy terrestrial planets.

We investigate the physical characteristics of the Solar System's proposed Planet Nine using modeling tools with a heritage in studying Uranus and Neptune. For a range of plausible masses and interior structures, we find upper limits on the intrinsic Teff, from ~35-50 K for masses of 5-20 M_Earth, and we also explore lower Teff values. Possible planetary radii could readily span from 3 to 6 R_Earth depending on the mass fraction of any H/He envelope. Given its cold temperature, the planet encounters significant methane condensation, which dramatically alters the atmosphere away from simple Neptune-like expectations. We find the atmosphere is strongly depleted in molecular absorption at visible wavelengths, suggesting a Rayleigh scattering atmosphere with a high geometric albedo approaching 0.75. We highlight two diagnostics for the atmosphere's temperature structure, the first being the value of the methane mixing ratio above the methane cloud. The second is the wavelength at which cloud scattering can be seen, which yields the cloud-top pressure. Surface reflection may be seen if the atmosphere is thin. Due to collision-induced opacity of H2 in the infrared, the planet would be extremely blue (instead of red) in the shortest wavelength WISE colors if methane is depleted, and would, in some cases, exist on the verge of detectability by WISE. For a range of models, thermal fluxes from ~3-5 microns are ~20 orders of magnitude larger than blackbody expectations. We report a search of the AllWISE Source Catalog for Planet Nine, but find no detection.

The recently discovered brown dwarf WISE 0855 presents our first opportunity to directly study an object outside the Solar System that is nearly as cold as our own gas giant planets. However the traditional methodology for characterizing brown dwarfs---near infrared spectroscopy---is not currently feasible as WISE 0855 is too cold and faint. To characterize this frozen extrasolar world we obtained a 4.5-5.2 $\mu$m spectrum, the same bandpass long used to study Jupiter's deep thermal emission. Our spectrum reveals the presence of atmospheric water vapor and clouds, with an absorption profile that is strikingly similar to Jupiter. The spectrum is high enough quality to allow the investigation of dynamical and chemical processes that have long been studied in Jupiter's atmosphere, but now on an extrasolar world.

As gas giant planets and brown dwarfs radiate away the residual heat from their formation, they cool through a spectral type transition from L to T, which encompasses the dissipation of cloud opacity and the appearance of strong methane absorption. While there are hundreds of known T-type brown dwarfs, the first generation of directly-imaged exoplanets were all L-type. Recently, Kuzuhara et al. (2013) announced the discovery of GJ 504 b, the first T dwarf exoplanet. GJ 504 b provides a unique opportunity to study the atmosphere of a new type of exoplanet with a ~500 K temperature that bridges the gap between the first directly imaged planets (~1000 K) and our own Solar System's Jupiter (~130 K). We observed GJ 504 b in three narrow L-band filters (3.71, 3.88, and 4.00 microns), spanning the red end of the broad methane fundamental absorption feature (3.3 microns) as part of the LEECH exoplanet imaging survey. By comparing our new photometry and literature photometry to a grid of custom model atmospheres, we were able to fit GJ 504 b's unusual spectral energy distribution for the first time. We find that GJ 504 b is well-fit by models with the following parameters: T_eff=544+/-10 K, g<600 m/s^2, [M/H]=0.60+/-0.12, cloud opacity parameter of f_sed=2-5, R=0.96+/-0.07 R_Jup, and log(L)=-6.13+/-0.03 L_Sun, implying a hot start mass of 3-30 M_jup for a conservative age range of 0.1-6.5 Gyr. Of particular interest, our model fits suggest that GJ 504 b has a super-stellar metallicity. Since planet formation can create objects with non-stellar metallicities, while binary star formation cannot, this result suggests that GJ 504 b formed like a planet, not like a binary companion.

We address the thermal history of the Earth after the Moon-forming impact, taking tidal heating and thermal blanketing by the atmosphere into account. The atmosphere sets an upper bound of ~100 W/m^2 on how quickly the Earth can cool. The liquid magma ocean cools over 2-10 Myrs, with longer times corresponding to high angular-momentum events. Tidal heating is focused mostly in mantle materials that are just beginning to freeze. The atmosphere's control over cooling sets up a negative feedback between viscosity-dependent tidal heating and temperature-dependent viscosity of the magma ocean. While the feedback holds, evolution of the Moon's orbit is limited by the modest radiative cooling rate of Earth's atmosphere. Orbital evolution is orders of magnitude slower than in conventional constant Q models, which promotes capture by resonances. The evection resonance is encountered early, when the Earth is molten. Capture by the evection resonance appears certain but unlikely to generate much eccentricity because it is encountered early when the Earth is molten and Q_Earth >> Q_Moon. Tidal dissipation in the Earth becomes more efficient (Q_Earth << Q_Moon) later when the Moon is between ~20 R_Earth and ~40 R_Earth. If lunar eccentricity grew great, this was when it did so, perhaps setting the table for some other process to leave its mark on the inclination of the Moon.

We explored two aspects of the problem of characterizing cool extrasolar giant planets in scattered optical light with a space based coronagraph. First, for a number of the known radial velocity (RV) giants we computed traditional forward models of their atmospheric structure and clouds, given various input assumptions, and computed model albedo spectra. Such models have been computed before, but mostly for generic planets. Our new models demonstrate that there is likely interesting spectral diversity among those planets that are most favorable for direct detection. Second, we applied a powerful Markov Chain Monte Carlo (MCMC) retrieval technique to synthetic noisy data of cool giants to better understand how well various atmospheric parameters--particularly molecular abundances and cloud properties--could be constrained. This is the first time such techniques have been applied to this problem. The process is time consuming, so only a dozen or so cases could be completed in the limited time available. Nevertheless the results clearly show that even at S/N ~ 5, scientifically interesting and valuable conclusions can be drawn about the properties of giant planet atmospheres from noisy spectra. We find that atmospheric abundances are best constrained when the planet gravity is bounded. Thus direct imaging observations of known radial velocity planets are extremely valuable as limits on a target planet's gravity can be obtained from astrometric imaging and reflectivity and mass-radius relationship arguments (for mass and radius, respectively). Further retrieval studies are clearly warranted and would be valuable to help guide instrument design decisions.

We present new calculations of Rosseland and Planck gaseous mean opacities relevant to the atmospheres of giant planets and ultracool dwarfs. Such calculations are used in modeling the atmospheres, interiors, formation, and evolution of these objects. Our calculations are an expansion of those presented in Freedman et al. (2008) to include lower pressures, finer temperature resolution, and also the higher metallicities most relevant for giant planet atmospheres. Calculations span 1 microbar to 300 bar, and 75 K to 4000 K, in a nearly square grid. Opacities at metallicities from solar to 50 times solar abundances are calculated. We also provide an analytic fit to the Rosseland mean opacities over the grid in pressure, temperature, and metallicity. In addition to computing mean opacities at these local temperatures, we also calculate them with weighting functions up to 7000 K, to simulate the mean opacities for incident stellar intensities, rather than locally thermally emitted intensities. The chemical equilibrium calculations account for the settling of condensates in a gravitational field and are applicable to cloud-free giant planet and ultracool dwarf atmospheres, but not circumstellar disks. We provide our extensive opacity tables for public use.

[Abridged] Aims: This work focuses on one lensed system, HATLAS J142935.3-002836 (H1429-0028), selected in the Herschel-ATLAS field. Gathering a rich, multi-wavelength dataset, we aim to confirm the lensing hypothesis and model the background source's morphology and dynamics, as well as to provide a full physical characterisation. Methods: Multi-wavelength high-resolution data is utilised to assess the nature of the system. A lensing-analysis algorithm which simultaneously fits different wavebands is adopted to characterise the lens. The background galaxy dynamical information is studied by reconstructing the 3-D source-plane of the ALMA CO(J:4-3) transition. Near-IR imaging from HST and Keck-AO allows to constrain rest-frame optical photometry independently for the foreground and background systems. Physical parameters (such as stellar and dust masses) are estimated via modelling of the spectral energy distribution taking into account source blending, foreground obscuration, and differential magnification. Results: The system comprises a foreground edge-on disk galaxy (at z_sp=0.218) with an almost complete Einstein ring around it. The background source (at z_sp=1.027) is magnified by a factor of ~8-10 depending on wavelength. It is comprised of two components and a tens of kpc long tidal tail resembling the Antennae merger. As a whole, the system is a massive stellar system (1.32[-0.41,+0.63] x1E11 Mo) forming stars at a rate of 394+-90 Mo/yr, and has a significant gas reservoir M_ISM = 4.6+-1.7 x1E10 Mo. Its depletion time due to star formation alone is thus expected to be tau_SF=M_ISM/SFR=117+-51 Myr. The dynamical mass of one of the components is estimated to be 5.8+-1.7 x1E10 Mo, and, together with the photometric total mass estimate, it implies that H1429-0028 is a major merger system (1:2.8[-1.5,+1.8]).

Brown dwarfs of a variety of spectral types have been observed to be photometrically variable. Previous studies have focused on objects at the L/T transition, where the iron and silicate clouds in L dwarfs break up or dissipate. However, objects outside of this transitional effective temperature regime also exhibit variability. Here, we present models for mid-late T dwarfs and Y dwarfs. We present models that include patchy salt and sulfide clouds as well as water clouds for the Y dwarfs. We find that for objects over 375 K, patchy cloud opacity would generate the largest amplitude variability within near-infrared spectral windows. For objects under 375 K, water clouds also become important and generate larger amplitude variability in the mid-infrared. We also present models in which we perturb the temperature structure at different pressure levels of the atmosphere to simulate hot spots. These models show the most variability in the absorption features between spectral windows. The variability is strongest at wavelengths that probe pressure levels at which the heating is the strongest. The most illustrative types of observations for understanding the physical processes underlying brown dwarf variability are simultaneous, multi-wavelength observations that probe both inside and outside of molecular absorption features.

Recent discoveries of variable brown dwarfs have provided us with a new window into their three-dimensional cloud structure. The highest variables are found at the L/T transition, where the cloud cover is thought to break up, but variability has been found to occur also for both cloudy L dwarfs and (mostly) cloud-free mid T dwarfs. We summarize results from recent HST programs measuring the spectral variability of brown dwarfs in the near-infrared and compare to results from ground-based programs. We discuss the patchy cloud structure of L/T transition objects, for which it is becoming increasingly certain that the variability does not arise from cloud holes into the deep hot regions but from varying cloud thickness. We present a new patchy cloud model to explain the spectral variability of 2MASSJ21392676+0220226. We also discuss the curious multi-wavelength variability behavior of the recently discovered very nearby early T dwarf WISE J104915.57-531906.1B (Luhman 16B) and the mid T dwarf 2MASS J22282889-431026.

The formation of clouds affects brown dwarf and planetary atmospheres of nearly all effective temperatures. Iron and silicate condense in L dwarf atmospheres and dissipate at the L/T transition. Minor species such as sulfides and salts condense in mid-late T dwarfs. For brown dwarfs below Teff=450 K, water condenses in the upper atmosphere to form ice clouds. Currently over a dozen objects in this temperature range have been discovered, and few previous theoretical studies have addressed the effect of water clouds on brown dwarf or exoplanetary spectra. Here we present a new grid of models that include the effect of water cloud opacity. We find that they become optically thick in objects below Teff=350-375 K. Unlike refractory cloud materials, water ice particles are significantly non-gray absorbers; they predominantly scatter at optical wavelengths through J band and absorb in the infrared with prominent features, the strongest of which is at 2.8 microns. H2O, NH3, CH4, and H2 CIA are dominant opacity sources; less abundant species such as may also be detectable, including the alkalis, H2S, and PH3. PH3, which has been detected in Jupiter, is expected to have a strong signature in the mid-infrared at 4.3 microns in Y dwarfs around Teff=450 K; if disequilibrium chemistry increases the abundance of PH3, it may be detectable over a wider effective temperature range than models predict. We show results incorporating disequilibrium nitrogen and carbon chemistry and predict signatures of low gravity in planetary- mass objects. Lastly, we make predictions for the observability of Y dwarfs and planets with existing and future instruments including the James Webb Space Telescope and Gemini Planet Imager.

It is now understood that the accretion of terrestrial planets naturally involves giant collisions, the moon-forming impact being a well known example. In the aftermath of such collisions the surface of the surviving planet is very hot and potentially detectable. Here we explore the atmospheric chemistry, photochemistry, and spectral signatures of post-giant-impact terrestrial planets enveloped by thick atmospheres consisting predominantly of CO2, and H2O. The atmospheric chemistry and structure are computed self-consistently for atmospheres in equilibrium with hot surfaces with composition reflecting either the bulk silicate Earth (which includes the crust, mantle, atmosphere and oceans) or Earth's continental crust. We account for all major molecular and atomic opacity sources including collision-induced absorption. We find that these atmospheres are dominated by H2O and CO2, while the formation of CH4, and NH3 is quenched due to short dynamical timescales. Other important constituents are HF, HCl, NaCl, and SO2. These are apparent in the emerging spectra, and can be indicative that an impact has occurred. The use of comprehensive opacities results in spectra that are a factor of 2 lower in surface brightness in the spectral windows than predicted by previous models. The estimated luminosities show that the hottest post-giant-impact planets will be detectable with near-infrared coronagraphs on the planned 30m-class telescopes. The 1-4um region will be most favorable for such detections, offering bright features and better contrast between the planet and a potential debris disk. We derive cooling timescales on the order of 10^5-10^6 Myrs, based on the modeled effective temperatures. This leads to the possibility of discovering tens of such planets in future surveys.

We present observations of SPT-S J053816-5030.8, a gravitationally-lensed dusty star forming galaxy (DSFG) at z = 2.7817, first discovered at millimeter wavelengths by the South Pole Telescope. SPT 0538-50 is typical of the brightest sources found by wide-field millimeter-wavelength surveys, being lensed by an intervening galaxy at moderate redshift (in this instance, at z = 0.441). We present a wide array of multi-wavelength spectroscopic and photometric data on SPT 0538-50, including data from ALMA, Herschel PACS and SPIRE, Hubble, Spitzer, VLT, ATCA, APEX, and the SMA. We use high resolution imaging from HST to de-blend SPT 0538-50, separating DSFG emission from that of the foreground lens. Combined with a source model derived from ALMA imaging (which suggests a magnification factor of 21 +/- 4), we derive the intrinsic properties of SPT 0538-50, including the stellar mass, far-IR luminosity, star formation rate, molecular gas mass, and - using molecular line fluxes - the excitation conditions within the ISM. The derived physical properties argue that we are witnessing compact, merger-driven star formation in SPT 0538-50, similar to local starburst galaxies, and unlike that seen in some other DSFGs at this epoch.

Upon its completion the Herschel ATLAS (H-ATLAS) will be the largest submillimetre survey to date, detecting close to half-a-million sources. It will only be possible to measure spectroscopic redshifts for a small fraction of these sources. However, if the rest-frame spectral energy distribution (SED) of a typical H-ATLAS source is known, this SED and the observed Herschel fluxes can be used to estimate the redshifts of the H-ATLAS sources without spectroscopic redshifts. In this paper, we use a subset of 40 H-ATLAS sources with previously measured redshifts in the range 0.5<z<4.2 to derive a suitable average template for high redshift H-ATLAS sources. We find that a template with two dust components T_c = 23.9 K, T_h = 46.9 K and ratio of mass of cold dust to mass of warm dust of 30.1) provides a good fit to the rest-frame fluxes of the sources in our calibration sample. We use a jackknife technique to estimate the accuracy of the redshifts estimated with this template, finding a root mean square of Delta z/(1+z) = 0.26. For sources for which there is prior information that they lie at z > 1 we estimate that the rms of Delta z/(1+z) = 0.12. We have used this template to estimate the redshift distribution for the sources detected in the H-ATLAS equatorial fields, finding a bimodal distribution with a mean redshift of 1.2, 1.9 and 2.5 for 250, 350 and 500 um selected sources respectively. \end{abstract}

The Sunyaev-Zel'dovich (SZ) effect spectrum crosses through a null where dT_CMB = 0 near nu_0 = 217 GHz. In a cluster of galaxies, nu_0 can be shifted from the canonical thermal SZ effect value by corrections to the SZ effect scattering due to the properties of the inter-cluster medium. We have measured the SZ effect in the hot galaxy cluster RX J1347.5-1145 with Z-Spec, an R ~ 300 grating spectrometer sensitive between 185 and 305 GHz. These data comprise a high spectral resolution measurement around the null of the SZ effect and clearly exhibit the transition from negative to positive dT_CMB over the Z-Spec band. The SZ null position is measured to be nu_0 = 225.8 \pm 2.5 (stat.) \pm 1.2 (sys.) GHz, which differs from the canonical null frequency by 3.0 sigma and is evidence for modifications to the canonical thermal SZ effect shape. Assuming the measured shift in nu_0 is due only to relativistic corrections to the SZ spectrum, we place the limit T_e = 17.1 \pm 5.3 keV from the zero-point measurement alone. By simulating the response of the instrument to the sky, we are able to generate likelihood functions in {y_0, T_e, v_pec} space. For v_pec = 0 km/s, we measure the best fitting SZ model to be y_0 = 4.6 (+0.6, -0.9) x 10^-4, T_e,0 = 15.2 (+12,-7.4) keV. When v_pec is allowed to vary, a most probable value of v_pec = +450 \pm 810 km/s is found.

Massive present-day early-type (elliptical and lenticular) galaxies probably gained the bulk of their stellar mass and heavy elements through intense, dust-enshrouded starbursts - that is, increased rates of star formation - in the most massive dark matter halos at early epochs. However, it remains unknown how soon after the Big Bang such massive starburst progenitors exist. The measured redshift distribution of dusty, massive starbursts has long been suspected to be biased low in redshift owing to selection effects, as confirmed by recent findings of systems out to redshift z~5. Here we report the identification of a massive starburst galaxy at redshift 6.34 through a submillimeter color-selection technique. We unambiguously determined the redshift from a suite of molecular and atomic fine structure cooling lines. These measurements reveal a hundred billion solar masses of highly excited, chemically evolved interstellar medium in this galaxy, which constitutes at least 40% of the baryonic mass. A "maximum starburst" converts the gas into stars at a rate more than 2,000 times that of the Milky Way, a rate among the highest observed at any epoch. Despite the overall downturn of cosmic star formation towards the highest redshifts, it seems that environments mature enough to form the most massive, intense starbursts existed at least as early as 880 million years after the Big Bang.

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 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 high-resolution maps of stars, dust, and molecular gas in a strongly lensed submillimeter galaxy (SMG) at z = 3.259. HATLAS12--00 is selected from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) as a strong lens candidate mainly based on its unusually high 500um flux density (~300 mJy). It is the only high-redshift Planck detection in the 130 deg^2 H-ATLAS Phase 1 area. Keck Adaptive Optics images reveal a quadruply imaged galaxy in the K-band while the Submillimeter Array and the Extended Very Large Array show doubly imaged 880um and CO(1-0) sources, indicating differentiated distributions of the various components in the galaxy. In the source plane, the stars reside in three major kpc-scale clumps extended over ~1.6 kpc, the dust in a compact (~1 kpc) region ~3 kpc north of the stars, and the cold molecular gas in an extended (~7 kpc) disk ~5 kpc northeast of the stars. The emission from the stars, dust, and gas are magnified by ~17, 8, and 7 times, respectively, by four lensing galaxies at z ~ 1. Intrinsically, the lensed galaxy is a warm (T_dust ~ 40-65 K), hyper-luminous (L_IR ~ 1.7e13 Lsun; SFR ~ 2000 Msun/yr), gas-rich (M_gas/M_baryon ~ 70%), young (M_stellar/SFR ~ 20 Myr), and short-lived (M_gas/SFR ~ 40 Myr) starburst, without a significant active galactic nucleus. With physical properties similar to unlensed z > 2 SMGs, HATLAS12--00 offers a detailed view of a typical SMG through a powerful cosmic microscope.

We report the detection of an extremely bright ($\sim$37 mJy at 1100 $\mu$m and $\sim$91 mJy at 880 $\mu$m) submillimeter galaxy (SMG), AzTEC-ASTE-SXDF1100.001 (hereafter referred to as SXDF1100.001 or Orochi), discovered in 1100 $\mu$m observations of the Subaru/XMM-Newton Deep Field using AzTEC on ASTE. Subsequent CARMA 1300 $\mu$m and SMA 880 $\mu$m observations successfully pinpoint the location of Orochi and suggest that it has two components, one extended (FWHM of $\sim$ 4$^{\prime\prime}$) and one compact (unresolved). Z-Spec on CSO has also been used to obtain a wide band spectrum from 190 to 308 GHz, although no significant emission/absorption lines are found. The derived upper limit to the line-to-continuum flux ratio is 0.1--0.3 (2 $\sigma$) across the Z-Spec band. Based on the analysis of the derived spectral energy distribution from optical to radio wavelengths of possible counterparts near the SMA/CARMA peak position, we suggest that Orochi is a lensed, optically dark SMG lying at $z \sim 3.4$ behind a foreground, optically visible (but red) galaxy at $z \sim 1.4$. The deduced apparent (i.e., no correction for magnification) infrared luminosity ($L_{\rm IR}$) and star formation rate (SFR) are $6 \times 10^{13}$ $L_{\odot}$ and 11000 $M_{\odot}$ yr$^{-1}$, respectively, assuming that the $L_{\rm IR}$ is dominated by star formation. These values suggest that Orochi will consume its gas reservoir within a short time scale ($3 \times 10^{7}$ yr), which is indeed comparable to those in extreme starbursts like the centres of local ULIRGs.

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 present the rest-frame 200--320 \mm\ spectrum of the z=3.91 quasar \apm, obtained with Z-Spec at the Caltech Submillimeter Observatory. In addition to the \jeight\ to \jthirteen\ CO rotational transitions which dominate the CO cooling, we find six transitions of water originating at energy levels ranging up to 643 K. Most are first detections at high redshift, and we have confirmed one transition with CARMA. The CO cooling is well-described by our XDR model, assuming L$_{\rm 1-100\,keV}\sim1\times10^{46}\rm\,erg\,s^{-1}$, and that the gas is distributed over a 550-pc sizescale, per the now-favored $\mu$=4 lensing model. The total observed cooling in water corresponds to 6.5$\times10^{9}$ \ls, comparable to that of CO. We compare the water spectrum with that of Mrk 231, finding that the intensity ratios among the high-lying lines are similar, but with a total luminosity scaled up by a factor of $\sim$50. Using this scaling, we estimate an average water abundance relative to \hh\ of 1.4$\times10^{-7}$, a good match to the prediction of the chemical network in the XDR model. As with Mrk 231, the high-lying water transitions are excited radiatively via absorption in the rest-frame far-infrared, and we show that the powerful dust continuum in \apm\ is more than sufficient to pump this massive reservoir of warm water vapor.

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $\mu \sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

We report on the redshift measurement and CO line excitation of HERMES J105751.1+573027 (HLSW-01), a strongly lensed submillimeter galaxy discovered in Herschel/SPIRE observations as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). HLSW-01 is an ultra-luminous galaxy with an intrinsic far-infrared luminosity of 1.4x10^(13) solar luminosities, and is lensed by a massive group of galaxies into at least four images with a total magnification of 10.9+/-0.7. With the 100 GHz instantaneous bandwidth of the Z-Spec instrument on the Caltech Submillimeter Observatory, we robustly identify a redshift of z=2.958+/-0.007 for this source, using the simultaneous detection of four CO emission lines (J = 7-6, J = 8-7, J = 9-8, and J = 10-9). Combining the measured line fluxes for these high-J transitions with the J = 1-0, J = 3-2 and J = 5-4 line fluxes measured with the Green Bank Telescope, the Combined Array for Research in Millimeter Astronomy, and the Plateau de Bure Interferometer, respectively, we model the physical properties of the molecular gas in this galaxy. We find that the full CO spectral line energy distribution is well described by warm, moderate-density gas with Tkin = 86-235 K and n(H2) = (1.1-3.5)x10^3 cm^(-3). However, it is possible that the highest-J transitions are tracing a small fraction of very dense gas in molecular cloud cores, and two-component models that include a warm/dense molecular gas phase with Tkin ~ 200 K, n(H2) ~ 10^5 cm^(-3) are also consistent with these data. Higher signal-to-noise measurements of the J(upper) > 7 transitions with high spectral resolution, combined with high spatial resolution CO maps, are needed to improve our understanding of the gas excitation, morphology, and dynamics of this interesting high-redshift galaxy.

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $\mu=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.

Partial frequency redistribution (PRD), describing the formation of the line profile, has negligible observational effects for optical depths smaller than ~10^3, at the resolving power of most current instruments. However, when the spectral resolution is sufficiently high, PRD modeling becomes essential in interpreting the line shapes and determining the total line fluxes. We demonstrate the effects of PRD on the H_2 line profiles observed at high spectral resolution by the Far-Ultraviolet Spectroscopic Explorer (FUSE) in the atmospheres of Jupiter and Saturn. In these spectra, the asymmetric shapes of the lines in the Lyman (v"- 6) progression pumped by the solar Ly-beta are explained by coherent scattering of the photons in the line wings. We introduce a simple computational approximation to mitigate the numerical difficulties of radiative transfer with PRD, and show that it reproduces the exact radiative transfer solution to better than 10%. The lines predicted by our radiative transfer model with PRD, including the H_2 density and temperature distribution as a function of height in the atmosphere, are in agreement with the line profiles observed by FUSE. We discuss the observational consequences of PRD, and show that this computational method also allows us to include PRD in modeling the continuum pumped H_2 fluorescence, treating about 4000 lines simultaneously.

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.

The mid-infrared (MIR) spectra of dense photodissociation regions (PDRs) are typically dominated by emission from polycyclic aromatic hydrocarbons (PAHs) and the lowest pure rotational states of molecular hydrogen (H2); two species which are probes of the physical properties of gas and dust in intense UV radiation fields. We utilize the high angular resolution of the Infrared Spectrograph on the Spitzer Space Telescope to construct spectral maps of the PAH and H2 features for three of the best studied PDRs in the galaxy, NGC 7023, NGC 2023 and IC 63. We present spatially resolved maps of the physical properties, including the H2 ortho-to-para ratio, temperature, and G_o/n_H. We also present evidence for PAH dehydrogenation, which may support theories of H2 formation on PAH surfaces, and a detection of preferential self-shielding of ortho-H2. All PDRs studied exhibit average temperatures of ~500 - 800K, warm H2 column densities of ~10^20 cm^-2, G_o/n_H ~ 0.1 - 0.8, and ortho-to-para ratios of ~ 1.8. We find that while the average of each of these properties is consistent with previous single value measurements of these PDRs, when available, the addition of spatial resolution yields a diversity of values with gas temperatures as high as 1500 K, column densities spanning ~ 2 orders of magnitude, and extreme ortho-to-para ratios of < 1 and > 3.

Observations of four comets made with the Far Ultraviolet Spectroscopic Explorer show the rotational envelope of the (0,0) band of the CO Hopfield-Birge system (C - X) at 1088 A to consist of both "cold" and "hot" components, the "cold" component accounting for ~75% of the flux and with a rotational temperature in the range 55-75 K. We identify the "hot" component as coming from the dissociation of CO2 into rotationally "hot" CO, with electron impact dissociation probably dominant over photodissociation near the nucleus. An additional weak, broad satellite band is seen centered near the position of the P(40) line that we attribute to CO fluorescence from a non-thermal high J rotational population produced by photodissociation of formaldehyde into CO and H2. This process also leaves the H2 preferentially populated in excited vibrational levels which are identified by fluorescent H2 lines in the spectrum excited by solar OVI 1031.9 and solar Lyman-alpha. The amount of H2 produced by H2CO dissociation is comparable to the amount produced by photodissociation of H2O. Electron impact excitation of CO, rather than resonance fluorescence, appears to be the primary source of the observed (B - X) (0,0) band at 1151 A.

We present the modeling of partial frequency redistribution (PRD) effects for the fluorescent emission lines of molecular hydrogen, the general computational approximations, and the applications to planetary atmospheres, as well as interstellar medium. Our model is applied to FUSE observations of Jupiter, Saturn, and reflection nebulae, allowing an independent confirmation of the H2 abundance and the structure of planetary atmospheres.

We explore the design of a space mission, Project Lyman, which has the goal of quantifying the ionization history of the universe from the present epoch to a redshift of z ~ 3. Observations from WMAP and SDSS show that before a redshift of z >~ 6 the first collapsed objects, possibly dwarf galaxies, emitted Lyman continuum (LyC) radiation shortward of 912 A, reionizing most of the universe. How LyC escapes from galactic environments, whether it induces positive or negative feedback on the local and global collapse of structures, and the role played by clumping, molecules, metallicity and dust are major unanswered theoretical questions, requiring observational constraint. Numerous intervening Lyman limit systems, which frustrate the detection of LyC from high z objects, thin below z ~ 3 where there are a few objects with apparently very high fesc. At low z there are only controversial detections and a handful of upper limits. A wide-field multi-object spectroscopic survey with moderate spectral and spatial resolution can quantify fesc within diverse spatially resolved galactic environments over redshifts with significant evolution in galaxy assemblage and quasar activity. It can also calibrate LyC escape against Ly-alpha escape, providing an essential tool to JWST for probing the beginnings of reionization. We present calculations showing the evolution of the characteristic apparent magnitude of star-forming galaxy luminosity functions at 900 A, as a function of redshift and assumed escape fraction to determine the required aperture for detecting LyC. We review our efforts to build a pathfinding dual order multi-object spectro/telescope with a (0.5deg)^2 field-of-view, using a GSFC microshutter array, and crossed delay-line micro-channel plate detector.

We present a preliminary calibration and flight performance of the Long-Slit Imaging Dual Order Spectrograph (LIDOS), a rocket-borne instrument with a large dynamic range in the 900 - 1700A bandpass. The instrument observes UV-bright objects with a CCD channel and fainter nebulosity with an MCP detector. The image quality and the detector quantum efficiencies were determined using the calibration and test equipment at the Johns Hopkins University, and further monitored using an on-board electron-impact calibration lamp. We review results from each of the three flights of the instrument.

The rich structure of the Fourth Positive System (A-X) of carbon monoxide accounts for many of the spectral features seen in long slit HST-STIS observations of comets 153P/Ikeya-Zhang, C/2001 Q4 (NEAT), and C/2000 WM1 (LINEAR), as well as in the HST-GHRS spectrum of comet C/1996 B2 Hyakutake. A detailed CO fluorescence model is developed to derive the CO abundances in these comets by simultaneously fitting all of the observed A-X bands. The model includes the latest values for the oscillator strengths and state parameters, and accounts for optical depth effects due to line overlap and self-absorption. The model fits yield radial profiles of CO column density that are consistent with a predominantly native source for all the comets observed by STIS. The derived CO abundances relative to water in these comets span a wide range, from 0.44% for C/2000 WM1 (LINEAR), 7.2% for 153P/Ikeya-Zhang, 8.8% for C/2001 Q4 (NEAT) to 20.9% for C/1996 B2 (Hyakutake). The subtraction of the CO spectral features using this model leads to the first identification of a molecular hydrogen line pumped by solar HI Lyman-beta longward of 1200A in the spectrum of comet 153P/Ikeya-Zhang. (Abridged)

We present a determination of the molecular gas mass in the AU Microscopii circumstellar disk. Direct detection of a gas component to the AU Mic disk has proven elusive, with upper limits derived from ultraviolet absorption line and submillimeter CO emission studies. Fluorescent emission lines of H2, pumped by the OVI 1032 resonance line through the C-X (1 -- 1) Q(3) 1031.87 \AA\ transition, are detected by the Far Ultraviolet Spectroscopic Explorer. These lines are used to derive the H2 column density associated with the AU Mic system. The derived column density is in the range N(H2) = 1.9 x 10^{17} - 2.8 x 10^{15} cm^{-2}, roughly two orders of magnitude lower than the upper limit inferred from absorption line studies. This range of column densities reflects the range of H2 excitation temperature consistent with the observations, T(H2) = 800 -- 2000 K, derived from the presence of emission lines excited by OVI in the absence of those excited by LyA. Within the observational uncertainties, the data are consistent with the H2 gas residing in the disk. The inferred N(H2) range corresponds to H2-to-dust ratios of < 1/30:1 and a total M(H2) = 4.0 x 10^{-4} - 5.8 x 10^{-6} Earth masses. We use these results to predict the intensity of the associated rovibrational emission lines of H2 at infrared wavelengths covered by ground-based instruments, HST-NICMOS, and the Spitzer-IRS.

High resolution spectroscopy with FUSE and STIS of atomic and molecular velocity stratification in the nebular outflow of M27 challenge models for the abundance kinematics in planetary nebulae. The simple picture of a very high speed (~ 1000 km/s), high ionization, radiation driven stellar wind surrounded by a slower (~ 10 km/s) mostly molecular outflow, with low ionization and neutral atomic species residing at the wind interaction interface, is not supported... We find ...there is a fast (33 -- 65 km/s) low ionization zone, surrounding a slower (<~ 33 km/s) high ionization zone and, at the transition velocity (33 km/s), vibrationally excited H_2 is intermixed with a predominately neutral atomic medium... Far-UV continuum fluorescence of H_2 is not detected, but Lyman alpha (Lya) fluorescence is present. The diffuse nebular medium is inhospitable to molecules and dust. Maintaining the modest equilibrium abundance of H_2 (N(H_2)/N(HI) << 1) in the diffuse nebular medium requires a source of H_2, mostly likely the clumpy nebular medium. The stellar SED shows no signs of reddening (E(B-V) < 0.01), but paradoxically measurements of Ha/Hb ... indicate E(B-V) ~ 0.1. ...the apparent enhancement of Ha/Hb in the absence of dust may result from a two step process of H_2 ionization by Lyman continuum (Lyc) photons followed by dissociative recombination (H_2 + gamma -> H_2^+ + e -> H(1s) + H (nl)), which ultimately produces fluorescence of Ha and Lya. In the optically thin limit at the inferred radius of the velocity transition we find dissociation of H_2 by stellar Lyc photons is an order of magnitude more efficient than spontaneous dissociation by far-UV photons. We suggest that the importance of this H_2 destruction process in HII regions has been overlooked.

We present new infrared observations of the emission/reflection nebula IC 405 obtained with the Spitzer Space Telescope. Infrared images in the four IRAC bands (3.6, 4.5, 5.8, and 8.0 um) and two MIPS bands (24 and 70 um) are complemented by IRS spectroscopy (5-30 um) of two nebular filaments. The IRAC (8.0 um) and MIPS imaging shows evidence of a bow shock associated with the runaway O9.5V star, HD 34078, created by the interaction between the star and nebular material. The ratio of emission at 24 to 70 um is higher in the immediate vicinity of HD 34078 than in the outer filaments, providing evidence for elevated dust temperatures (T_d > 90 K) in the shock region. The nebular imaging reveals that the morphology is band dependent, with varying contributions from aromatic emission features, H2, and dust emission. Nebular spectroscopy is used to quantify these contributions, showing several aromatic emission bands between 6-14 um, the S(5), S(3), S(2), and S(1) pure rotational emission lines of H2, and atomic fine structure lines of Ne, S, and Ar. The low-dispersion spectra provide constraints on the ionization state of the large molecules responsible for the aromatic infrared features. H2 rotational temperatures of the two bright nebular filaments are determined from the observed line strengths. An average T(H2) ~ 400 K is inferred, with evidence for additional non-uniform excitation by UV photons in the intense radiation field of HD 34078. The photoexcitation hypothesis is supported by direct measurement of the far-UV H2 fluorescence spectrum, obtained with FUSE.

One of the goals of the Hubble Space Telescope program to observe periodic comet 9P/Tempel 1 in conjunction with NASA's Deep Impact mission was to study the generation and evolution of the gaseous coma resulting from the impact. For this purpose, the Solar Blind Channel of the Advanced Camera for Surveys was used with the F140LP filter which is sensitive primarily to the ultraviolet emission (>1400 A) from the CO Fourth Positive system. Following the impact we detected an increase in brightness, which if all due to CO corresponds to 1.5 x 10^31 molecules or a mass of 6.6 x 10^5 kg, an amount that would normally be produced by 7-10 hours of quiescent outgassing from the comet. This number is less than or equal to 10% of the number of water molecules excavated, and suggests that the volatile content of the material excavated by the impact did not differ significantly from the surface or near sub-surface material responsible for the quiescent outgassing of the comet.

We report the first observation of Ly$\alpha$ pumped molecular hydrogen emission lines in planetary nebulae. The H$_{2}$ emission observed in the ultraviolet spectra of NGC 6853 and NGC 3132 can be explained by excitation of vibrationally hot H$_{2}$ by Ly$\alpha$ photons. Constraints are placed on the nebular Ly$\alpha$ emission profile, as well as the molecular hydrogen temperature, column density and turbulent motion. These parameters are similar for the two nebulae, pointing to similar physical conditions in these objects. The ro-vibrational cascade following Ly$\alpha$ pumping is predicted to have low surface brightness signatures in the visible and near infrared.