Showing 1–35 of 35 results for author: Lupu, R E

The impact of ultraviolet heating and cooling on the dynamics and observability of lava planet atmospheres

Authors: T. Giang Nguyen, Nicolas B. Cowan, Raymond T. Pierrehumbert, Roxana E. Lupu, John E. Moores

Abstract: Lava planets have non-global, condensible atmospheres similar to icy bodies within the solar system. Because they depend on interior dynamics, studying the atmospheres of lava planets can lead to understanding unique geological processes driven by their extreme environment. Models of lava planet atmospheres have thus far focused on either radiative transfer or hydrodynamics. In this study, we coup… ▽ More Lava planets have non-global, condensible atmospheres similar to icy bodies within the solar system. Because they depend on interior dynamics, studying the atmospheres of lava planets can lead to understanding unique geological processes driven by their extreme environment. Models of lava planet atmospheres have thus far focused on either radiative transfer or hydrodynamics. In this study, we couple the two processes by introducing ultraviolet and infrared radiation to a turbulent boundary layer model. We also test the effect of different vertical temperature profiles on atmospheric dynamics. Results from the model show that UV radiation affects the atmosphere much more than IR. UV heating and cooling work together to produce a horizontally isothermal atmosphere away from the sub-stellar point regardless of the vertical temperature profile. We also find that stronger temperature inversions induce stronger winds and hence cool the atmosphere. Our simulated transmission spectra of the bound atmosphere show a strong SiO feature in the UV that would be challenging to observe in the planet's transit spectrum due to the precision required. Our simulated emission spectra are more promising, with significant SiO spectral features at 4.5 and 9 $μ$m that can be observed with the James Webb Space Telescope. Different vertical temperature profiles produce discernible dayside emission spectra, but not in the way one would expect. △ Less

Submitted 9 May, 2022; v1 submitted 28 April, 2022; originally announced April 2022.

Comments: Accepted to MNRAS

Following the Lithium: Tracing Li-bearing Molecules Across Age, Mass, and Gravity in Brown Dwarfs

Authors: Ehsan Gharib-Nezhad, Mark S. Marley, Natasha E. Batalha, Channon Visscher, Richard S. Freedman, Roxana E. Lupu

Abstract: Lithium is an important element for the understanding of ultracool dwarfs because it is lost to fusion at masses above $\sim 68\, M_{\rm J}$. Hence, the presence or absence of atomic Li has served as an indicator of the nearby H-burning boundary at about $75\,M_{\rm J}$ between brown-dwarfs and very low-mass stars. Historically the "Lithium test", a search for the presence and strength of the Li l… ▽ More Lithium is an important element for the understanding of ultracool dwarfs because it is lost to fusion at masses above $\sim 68\, M_{\rm J}$. Hence, the presence or absence of atomic Li has served as an indicator of the nearby H-burning boundary at about $75\,M_{\rm J}$ between brown-dwarfs and very low-mass stars. Historically the "Lithium test", a search for the presence and strength of the Li line at 670.8 nm, has been a marker if an object has a substellar mass with stellar-like spectral energy distribution (e.g., a late-type M dwarf). While the Li test could in principle also be used to distinguish masses of later-type L-T dwarfs, Li is predominantly no longer found as an atomic gas, but rather a molecular species such as LiH, LiF, LiOH, and LiCl in their cooler atmospheres. L- and T-type brown dwarfs are also quite faint at 670 nm and thus challenging targets for high resolution spectroscopy. But only recently have experimental molecular line lists become available for the molecular Li species, allowing molecular Li mass discrimination. In this study, we generated the latest opacity of each of these Li-bearing molecules and performed thermochemical equilibrium atmospheric composition calculation of the abundance of these molecules. Finally, we computed thermal emission spectra for a series of radiative-convective equilibrium models of cloudy and cloudless brown dwarf atmospheres (with $T_{\rm eff}=$ 500--2400~K, and $\log g$=4.0, 4.5, 5.0) to understand where the presence or absence of atmospheric lithium-bearing species is most easily detected as a function of brown dwarf mass and age. After atomic Li, the best spectral signatures were found to be LiF at $10.5-12.5$~\micron and LiCl at $14.5-18.5$ $\micron$. LiH also shows a narrow feature at $\sim 9.38$ $\micron$. △ Less

Submitted 1 June, 2021; originally announced June 2021.

Comments: 13 pages, 7 figures, 1 table, In final revision at ApJ. Comments welcome

The Nancy Grace Roman Space Telescope Coronagraph Instrument (CGI) Technology Demonstration

Authors: N. Jeremy Kasdin, Vanessa P. Bailey, Bertrand Mennesson, Robert T. Zellem, Marie Ygouf, Jason Rhodes, Thomas Luchik, Feng Zhao, A J Eldorado Riggs, Young-Joon Seo, John Krist, Brian Kern, Hong Tang, Bijan Nemati, Tyler D. Groff, Neil Zimmerman, Bruce Macintosh, Margaret Turnbull, John Debes, Ewan S. Douglas, Roxana E. Lupu

Abstract: The Coronagraph Instrument (CGI) on the Nancy Grace Roman Space Telescope will demonstrate the high-contrast technology necessary for visible-light exoplanet imaging and spectroscopy from space via direct imaging of Jupiter-size planets and debris disks. This in-space experience is a critical step toward future, larger missions targeted at direct imaging of Earth-like planets in the habitable zone… ▽ More The Coronagraph Instrument (CGI) on the Nancy Grace Roman Space Telescope will demonstrate the high-contrast technology necessary for visible-light exoplanet imaging and spectroscopy from space via direct imaging of Jupiter-size planets and debris disks. This in-space experience is a critical step toward future, larger missions targeted at direct imaging of Earth-like planets in the habitable zones of nearby stars. This paper presents an overview of the current instrument design and requirements, highlighting the critical hardware, algorithms, and operations being demonstrated. We also describe several exoplanet and circumstellar disk science cases enabled by these capabilities. A competitively selected Community Participation Program team will be an integral part of the technology demonstration and could perform additional CGI observations beyond the initial tech demo if the instrument performance warrants it. △ Less

Submitted 2 March, 2021; originally announced March 2021.

Comments: SPIE Space Telescopes and Instrumentation 2020: Optical, Infrared, and Millimeter Wave, 114431U (15 December 2020)

Report number: Proc. SPIE 11443

Prospects for Characterizing the Haziest Sub-Neptune Exoplanets with High Resolution Spectroscopy

Authors: Callie E. Hood, Jonathan J. Fortney, Michael R. Line, Emily C. Martin, Caroline V. Morley, Jayne L. Birkby, Zafar Rustamkulov, Roxana E. Lupu, Richard S. Freedman

Abstract: Observations to characterize planets larger than Earth but smaller than Neptune have led to largely inconclusive interpretations at low spectral resolution due to hazes or clouds that obscure molecular features in their spectra. However, here we show that high-resolution spectroscopy (R $\sim$ 25,000 to 100,000) enables one to probe the regions in these atmospheres above the clouds where the cores… ▽ More Observations to characterize planets larger than Earth but smaller than Neptune have led to largely inconclusive interpretations at low spectral resolution due to hazes or clouds that obscure molecular features in their spectra. However, here we show that high-resolution spectroscopy (R $\sim$ 25,000 to 100,000) enables one to probe the regions in these atmospheres above the clouds where the cores of the strongest spectral lines are formed. We present models of transmission spectra for a suite of GJ1214b-like planets with thick photochemical hazes covering 1 - 5 $μ$m at a range of resolutions relevant to current and future ground-based spectrographs. Furthermore, we compare the utility of the cross-correlation function that is typically used with a more formal likelihood-based approach, finding that only the likelihood based method is sensitive to the presence of haze opacity. We calculate the signal-to-noise of these spectra, including telluric contamination, required to robustly detect a host of molecules such as CO, CO$_{2}$, H$_{2}$O, and CH$_{4}$, and photochemical products like HCN, as a function of wavelength range and spectral resolution. Spectra in M band require the lowest S/N$_{res}$ to detect multiple molecules simultaneously. CH$_{4}$ is only observable for the coolest models ($T_{\rm{eff}} =$ 412 K) and only in the L band. We quantitatively assess how these requirements compare to what is achievable with current and future instruments, demonstrating that characterization of small cool worlds with ground-based high resolution spectroscopy is well within reach. △ Less

Submitted 25 August, 2020; originally announced August 2020.

Comments: Submitted to AAS Journals, revised to reflect referee comments. Posting of this manuscript on the arXiv was coordinated with S. Ghandi et al

Detecting and Characterizing Water Vapor in the Atmospheres of Earth Analogs through Observation of the 0.94 Micron Feature in Reflected Light

Authors: Adam J. R. W. Smith, Y. Katherina Feng, Jonathan J. Fortney, Tyler D. Robinson, Mark S. Marley, Roxana E. Lupu, Nikole K. Lewis

Abstract: 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… ▽ More 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$\,μ$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\,μ$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σ$ 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$\,μ$m photometric point, instead of spectroscopy, combined with the green-optical point, fails to produce meaningful information about atmospheric water content. △ Less

Submitted 20 November, 2019; originally announced November 2019.

Comments: Accepted to AJ

Exploring A Photospheric Radius Correction to Model Secondary Eclipse Spectra for Transiting Exoplanets

Authors: Jonathan J. Fortney, Roxana E. Lupu, Caroline V. Morley, Richard S. Freedman, Callie Hood

Abstract: 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… ▽ More 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. △ Less

Submitted 29 March, 2019; originally announced April 2019.

Comments: Submitted to ApJ Letters, revised to reflect referee comments

The Effect of 3D Transport-Induced Disequilibrium Carbon Chemistry on the Atmospheric Structure and Phase Curves and Emission Spectra of Hot Jupiter HD 189733b

Authors: Maria E Steinrueck, Vivien Parmentier, Adam P Showman, Joshua D Lothringer, Roxana E Lupu

Abstract: 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… ▽ More 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 $μ$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 $μ$m bands, are strongly impacted by disequilibrium abundances. We expect chemical quenching to result in much larger nightside fluxes in the 3.6 $μ$m band, in stark contrast to observations. Meanwhile, we find almost no effect on predicted observations in the 4.5 $μ$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 $μ$m band. △ Less

Submitted 21 May, 2019; v1 submitted 6 August, 2018; originally announced August 2018.

Comments: 28 pages, 12 figures; accepted for publication in ApJ. Key changes since v1: updated introduction and discussion to include recent literature more thoroughly; added paragraph and figure to clarify change in water abundance; ran simulations for additional 600 days simulation time (leading to small quantitative changes in some figures but no qualitative changes)

Characterizing Earth Analogs in Reflected Light: Atmospheric Retrieval Studies for Future Space Telescopes

Authors: Y. Katherina Feng, Tyler D. Robinson, Jonathan J. Fortney, Roxana E. Lupu, Mark S. Marley, Nikole K. Lewis, Bruce Macintosh, Michael R. Line

Abstract: 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 expl… ▽ More 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 $μ$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 $μ$m and $R = 50$ spectroscopy from 0.6-0.97 $μ$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. △ Less

Submitted 16 March, 2018; originally announced March 2018.

Comments: Resubmitted to AAS Journals after incorporating reviewer feedback. 26 pages, 18 figure, 9 tables

Direct detection and characterization of M-dwarf planets using light echoes

Authors: William B. Sparks, Richard L. White, Roxana E. Lupu, Holland C. Ford

Abstract: 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… ▽ More 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. △ Less

Submitted 3 January, 2018; originally announced January 2018.

High Efficiency UV/Optical/NIR Detectors for Large Aperture Telescopes and UV Explorer Missions: Development of and Field Observations with Delta-doped Arrays

Authors: Shouleh Nikzad, April D. Jewell, Michael E. Hoenk, Todd Jones, John Hennessy, Tim Goodsall, Alexander Carver, Charles Shapiro, Samuel R. Cheng, Erika Hamden, Gillian Kyne, D. Christopher Martin, David Schiminovich, Paul Scowen, Kevin France, Stephan McCandliss, Roxana E. Lupu

Abstract: 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… ▽ More 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. △ Less

Submitted 12 January, 2017; v1 submitted 14 December, 2016; originally announced December 2016.

Comments: submitted to JATIS

Exo-Transmit: An Open-Source Code for Calculating Transmission Spectra for Exoplanet Atmospheres of Varied Composition

Authors: Eliza M. -R. Kempton, Roxana E. Lupu, Albert Owusu-Asare, Patrick Slough, Bryson Cale

Abstract: 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 w… ▽ More 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 . △ Less

Submitted 11 November, 2016; originally announced November 2016.

Comments: 8 pages, 7 figures, 1 table, submitted to PASP

The Hunt for Planet Nine: Atmosphere, Spectra, Evolution, and Detectability

Authors: Jonathan J. Fortney, Mark S. Marley, Gregory Laughlin, Nadine Nettelmann, Caroline V. Morley, Roxana E. Lupu, Channon Visscher, Pavle Jeremic, Wade G. Khadder, Mason Hargrave

Abstract: 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_Eart… ▽ More 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. △ Less

Submitted 6 June, 2016; v1 submitted 25 April, 2016; originally announced April 2016.

Comments: Accepted to ApJ Letters, with minor changes from the submitted version

Developing Atmospheric Retrieval Methods for Direct Imaging Spectroscopy of Gas Giants in Reflected Light I: Methane Abundances and Basic Cloud Properties

Authors: Roxana E. Lupu, Mark S. Marley, Nikole Lewis, Michael Line, Wesley A. Traub, Kevin Zahnle

Abstract: 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 techniqu… ▽ More 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. △ Less

Submitted 16 September, 2016; v1 submitted 18 April, 2016; originally announced April 2016.

Comments: 37 pages, 30 figures, accepted for publication in AJ

Gaseous Mean Opacities for Giant Planet and Ultracool Dwarf Atmospheres over a Range of Metallicities and Temperatures

Authors: Richard S. Freedman, Jacob Lustig-Yaeger, Jonathan J. Fortney, Roxana E. Lupu, Mark S. Marley, Katharina Lodders

Abstract: 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 th… ▽ More 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. △ Less

Submitted 24 June, 2021; v1 submitted 29 August, 2014; originally announced September 2014.

Comments: This version fixes a units description error that was only in the arxiv version, but is correct in the ApJS version. In equations 4 and 5 pressure is in dyne/cm2. Online tables available at the ApJS website: https://iopscience.iop.org/article/10.1088/0067-0049/214/2/25

Herschel-ATLAS and ALMA: HATLAS J142935.3-002836, a lensed major merger at redshift 1.027

Authors: Hugo Messias, Simon Dye, Neil Nagar, Gustavo Orellana, R. Shane Bussmann, Jae Calanog, Helmut Dannerbauer, Hai Fu, Edo Ibar, Andrew Inohara, R. J. Ivison, Mattia Negrello, Dominik A. Riechers, Yun-Kyeong Sheen, Simon Amber, Mark Birkinshaw, Nathan Bourne, Dave L. Clements, Asantha Cooray, Gianfranco De Zotti, Ricardo Demarco, Loretta Dunne, Stephen Eales, Simone Fleuren, Roxana E. Lupu , et al. (7 additional authors not shown)

Abstract: [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 asses… ▽ More [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]). △ Less

Submitted 18 June, 2014; originally announced June 2014.

Comments: submitted to A&A, 20 pages, 12 figures

Journal ref: A&A 568, A92 (2014)

Cloud structure of brown dwarfs from spectroscopic variability observations

Authors: Esther Buenzli, Mark S. Marley, Daniel Apai, Roxana E. Lupu

Abstract: 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 spectr… ▽ More 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. △ Less

Submitted 1 June, 2014; originally announced June 2014.

Comments: 4 pages, 1 figure, submitted to the proceedings of "Gaia and the unseen - the brown dwarf question", Torino, 24-26 March 2014, to be published in Memorie della Societa' Astronomica Italiana (SAIt), eds Ricky Smart, David Barrado, Jackie Faherty

The Atmospheres of Earth-like Planets after Giant Impact Events

Authors: R. E. Lupu, Kevin Zahnle, Mark S. Marley, Laura Schaefer, Bruce Fegley, Caroline Morley, Kerri Cahoy, Richard Freedman, Jonathan J. Fortney

Abstract: 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… ▽ More 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. △ Less

Submitted 6 January, 2014; originally announced January 2014.

Comments: 51 pages, 16 figures, preprint format; ApJ submitted

SPT 0538-50: Physical conditions in the ISM of a strongly lensed dusty star-forming galaxy at z=2.8

Authors: M. S. Bothwell, J. E. Aguirre, S. C. Chapman, D. P. Marrone, J. D. Vieira, M. L. N. Ashby, M. Aravena, B. A. Benson, J. J. Bock, C. M. Bradford, M. Brodwin, J. Carlstrom, T. M. Crawford, C. de Breuck, T. P. Downes, C. D. Fassnacht, A. H. Gonzalez, T. R. Greve, B. Gullberg, Y. Hezaveh, G. P. Holder, W. L. Holzapfel, E. Ibar, R. Ivison, J. Kamenetzky , et al. (16 additional authors not shown)

Abstract: 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 w… ▽ More 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. △ Less

Submitted 12 September, 2013; originally announced September 2013.

Comments: 16 pages, 11 figures. Accepted for publication in ApJ

H-ATLAS: estimating redshifts of Herschel sources from sub-mm fluxes

Authors: E. A. Pearson, S. Eales, L. Dunne, J. Gonzalez Nuevo, S. Maddox, J. E. Aguirre, M. Baes, A. J. Baker, N. Bourne, C. M. Bradford, C. J. R. Clark, A. Cooray, A. Dariush, G. De Zotti, S. Dye, D. Frayer, H. L. Gomez, A. I. Harris, R. Hopwood, E. Ibar, R. J. Ivison, M. Jarvis, M. Krips, A. Lapi, R. E. Lupu , et al. (7 additional authors not shown)

Abstract: 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 es… ▽ More 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} △ Less

Submitted 26 August, 2013; originally announced August 2013.

Comments: 13 pages, 11 figures

A Dust-Obscured Massive Maximum-Starburst Galaxy at a Redshift of 6.34

Authors: Dominik A. Riechers, C. M. Bradford, D. L. Clements, C. D. Dowell, I. Perez-Fournon, R. J. Ivison, C. Bridge, A. Conley, Hai Fu, J. D. Vieira, J. Wardlow, J. Calanog, A. Cooray, P. Hurley, R. Neri, J. Kamenetzky, J. E. Aguirre, B. Altieri, V. Arumugam, D. J. Benford, M. Bethermin, J. Bock, D. Burgarella, A. Cabrera-Lavers, S. C. Chapman , et al. (39 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 15 April, 2013; originally announced April 2013.

Comments: 35 pages, 19 figures, 5 tables, including Supplementary Information (version formatted by the authors). To appear in Nature. Under press embargo until 18:00 London time/13:00 US Eastern Time on 17 April 2013

Journal ref: Nature 496:329-333,2013

Redshift Determination and CO Line Excitation Modeling for the Multiply-Lensed Galaxy HLSW-01

Authors: K. S. Scott, R. E. Lupu, J. E. Aguirre, R. Auld, H. Aussel, A. J. Baker, A. Beelen, J. Bock, C. M. Bradford, D. Brisbin, D. Burgarella, J. M. Carpenter, P. Chanial, S. C. Chapman, D. L. Clements, A. Conley, A. Cooray, P. Cox, C. D. Dowell, S. Eales, D. Farrah, A. Franceschini, D. T. Frayer, R. Gavazzi, J. Glenn , et al. (38 additional authors not shown)

Abstract: 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 gal… ▽ More 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. △ Less

Submitted 20 April, 2011; originally announced April 2011.

Comments: 18 pages, 5 figures, 2 table. Accepted for publication in ApJ

Observations and modeling of H_2 fluorescence with partial frequency redistribution in giant planet atmospheres

Authors: Roxana E. Lupu, Paul D. Feldman, Stephan R. McCandliss, Darrell F. Strobel

Abstract: 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 effec… ▽ More 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. △ Less

Submitted 25 February, 2011; originally announced February 2011.

Comments: 17 pages, accepted for publication in ApJ

The Detection of a Population of Submillimeter-Bright, Strongly-Lensed Galaxies

Authors: Mattia Negrello, R. Hopwood, G. De Zotti, A. Cooray, A. Verma, J. Bock, D. T. Frayer, M. A. Gurwell, A. Omont, R. Neri, H. Dannerbauer, L. L. Leeuw, E. Barton, J. Cooke, S. Kim, E. da Cunha, G. Rodighiero, P. Cox, D. G. Bonfield, M. J. Jarvis, S. Serjeant, R. J. Ivison, S. Dye, I. Aretxaga, D. H. Hughes , et al. (64 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 4 November, 2010; originally announced November 2010.

Comments: This is the author's version of the work. It is posted here by permission of the AAAS for personal use, not for redistribution. The definitive version was published in Science VOL 330, (Nov 5th, 2010), doi:10.1126/science.1193420

Spitzer Mapping of PAHs and H2 in Photodissociation Regions

Authors: Brian T. Fleming, Kevin France, Roxana E. Lupu, Stephan R. McCandliss

Abstract: 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 S… ▽ More 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. △ Less

Submitted 30 September, 2010; originally announced October 2010.

Comments: 14 figures

Journal ref: 2010ApJ...725..159F

Measurements of CO redshifts with Z-Spec for lensed submillimeter galaxies discovered in the H-ATLAS survey

Authors: R. E. Lupu, K. S. Scott, J. E. Aguirre, I. Aretxaga, R. Auld, E. Barton, A. Beelen, F. Bertoldi, J. J. Bock, D. Bonfield, C. M. Bradford, S. Buttiglione, A. Cava, D. L. Clements, J. Cooke, A. Cooray, H. Dannerbauer, A. Dariush, G. De Zotti, L. Dunne, S. Dye, S. Eales, D. Frayer, J. Fritz, J. Glenn , et al. (28 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 17 August, 2012; v1 submitted 29 September, 2010; originally announced September 2010.

Comments: 28 pages, 15 figures, accepted for publication in ApJ, October 1, 2012, Issue 757 - 2

The Far Ultraviolet Spectral Signatures of Formaldehyde and Carbon Dioxide in Comets

Authors: Paul D. Feldman, Roxana E. Lupu, Stephan R. McCandliss, Harold A. Weaver

Abstract: 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… ▽ More 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. △ Less

Submitted 11 May, 2009; originally announced May 2009.

Comments: 21 pages, 7 figures, accepted for publication in the Astrophysical Journal

Journal ref: Astrophys.J.699:1104-1112,2009

Modeling H2 Fluorescence in Planetary Atmospheres with Partial Frequency Redistribution

Authors: R. E. Lupu, P. D. Feldman, S. R. McCandliss, K. France

Abstract: 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 an… ▽ More 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. △ Less

Submitted 10 December, 2008; originally announced December 2008.

Comments: To appear in Proc. Future Directions in Ultraviolet Spectroscopy, Annapolis, Oct. 20-22 2008

Project Lyman

Authors: Stephan R. McCandliss, Jeffrey W. Kruk, William P. Blair, Mary Elizabeth Kaiser, Paul D. Feldman, Gerhardt R. Meurer, William V. Dixon, David J. Sahnow, David A. Neufeld, Roxana E. Lupu, Brian Fleming, Stephen A. Smee, B. G. Andersson, Samuel H. Moseley, Alexander S. Kutyrev, Mary J. Li, George Sonneborn, Oswald H. W. Siegmund, John V. Vallerga, Barry Y. Welsh, Massimo Stiavelli, Rogier A. Windhorst, Alice E. Shapley

Abstract: 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… ▽ More 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. △ Less

Submitted 15 July, 2008; originally announced July 2008.

Comments: SPIE oral paper 7011-76 presented at Astronomical Telescopes 2008 - 23 -- 28 June Marseille, France, 12 page, 5 figures

Calibration and flight performance of the long-slit imaging dual order spectrograph

Authors: Roxana E. Lupu, Stephan R. McCandliss, Brian Fleming, Kevin France, Paul D. Feldman, Shouleh Nikzad

Abstract: 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 calibratio… ▽ More 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. △ Less

Submitted 10 June, 2008; originally announced June 2008.

Comments: 12 pages. to appear in Proc. SPIE 7011

The Fourth Positive System of Carbon Monoxide in the Hubble Space Telescope Spectra of Comets

Authors: Roxana E. Lupu, Paul D. Feldman, Harold A. Weaver, Gian-Paolo Tozzi

Abstract: 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 simultaneou… ▽ More 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) △ Less

Submitted 22 August, 2007; originally announced August 2007.

Comments: 12 pages, 11 figures, ApJ accepted

A Low-Mass H2 Component to the AU Microscopii Circumstellar Disk

Authors: Kevin France, Aki Roberge, Roxana E. Lupu, Seth Redfield, Paul D. Feldman

Abstract: 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 Å transition, are de… ▽ More 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 Å 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. △ Less

Submitted 6 July, 2007; originally announced July 2007.

Comments: ApJ, accepted. 9 pages, 4 figures, emulateapj

Molecular and Atomic Excitation Stratification in the Outflow of the Planetary Nebula M27

Authors: Stephan R. McCandliss, Kevin France, Roxana E. Lupu, Eric B. Burgh, Kenneth Sembach, Jeffrey Kruk, B. -G. Andersson, Paul D. Feldman

Abstract: 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 sp… ▽ More 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. △ Less

Submitted 15 January, 2007; originally announced January 2007.

Comments: In emulateapj 2 column, 28 pages total, 18 figures. Accepted for publication in the Astrophysical Journal on 5 January 2007. Abstract abridged

Journal ref: Astrophys.J.659:1291-1316,2007

A Cometary Bow Shock and Mid-Infrared Emission Variations Revealed in Spitzer Observations of HD 34078 and IC 405

Authors: Kevin France, Stephan R. McCandliss, Roxana E. Lupu

Abstract: 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… ▽ More 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. △ Less

Submitted 31 October, 2006; originally announced October 2006.

Comments: ApJ accepted, uses emulateapj, images compressed for astro-ph

Journal ref: Astrophys.J.655:920-939,2007

Carbon Monoxide in Comet 9P/Tempel 1 before and after the Deep Impact Encounter

Authors: Paul D. Feldman, Roxana E. Lupu, Stephan R. McCandliss, Harold A. Weaver, Michael F. A'Hearn, Michael J. S. Belton, Karen J. Meech

Abstract: 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 (>1… ▽ More 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. △ Less

Submitted 10 July, 2006; originally announced July 2006.

Comments: 10 pages, 4 postscript figures. Accepted for publication in ApJ (Letters)

Journal ref: Astrophys.J.647:L61-L64,2006

Discovery of Lyman-alpha Pumped Molecular Hydrogen Emission in the Planetary Nebulae NGC 6853 and NGC 3132

Authors: R. E. Lupu, K. France, S. R. McCandliss

Abstract: We report the first observation of Ly$α$ 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$α$ photons. Constraints are placed on the nebular Ly$α$ emission profile, as well as the molecular hydrogen temperature, column density and turbul… ▽ More We report the first observation of Ly$α$ 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$α$ photons. Constraints are placed on the nebular Ly$α$ 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$α$ pumping is predicted to have low surface brightness signatures in the visible and near infrared. △ Less

Submitted 24 February, 2006; originally announced February 2006.

Comments: emulate ApJ, 9 pages, 8 figures, 9 tables, ApJ accepted

Journal ref: Astrophys.J. 644 (2006) 981-989