BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Kevin Heng

Kevin Heng appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.GAphysics.ao-phastro-ph.SRastro-ph.IMastro-ph.HEphysics.chem-phphysics.class-phastro-phastro-ph.COphysics.atom-phphysics.flu-dyn

Catalog footprint

What is connected

73works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

73 published item(s)

preprint2024arXiv

Small but mighty: High-resolution spectroscopy of ultra-hot Jupiter atmospheres with compact telescopes. KELT-9 b's transmission spectrum with Wendelstein's FOCES Spectrograph

When observing transmission spectra produced by atmospheres of ultra-hot Jupiters, large telescopes are typically the instrument of choice due to the very weak signal of the planet's atmosphere. This study aims to alleviate the desire for large telescopes by illustrating that the same science is possible with smaller telescope classes. We use the cross-correlation technique to showcase the potential of the high-resolution spectrograph FOCES at Wendelstein Observatory and demonstrate its potential to resolve the atmosphere of the ultra-hot Jupiter, KELT-9 b. A performance comparison is conducted between FOCES and HARPS-N spectrographs, considering both single transit and combined observations over three nights. With FOCES, we have detected seven species in KELT-9 b's atmosphere: Ti II, Fe I, Fe II, Na I, Mg I, Na II, Cr II, Sc II. Although HARPS-N surpasses FOCES in performance, our results reveal that smaller telescope classes are capable of resolving ultra-hot Jupiter atmospheres. This broadens the scope of potential studies, allowing for investigations into phenomena like temporal variations in atmospheric signals and the atmospheric loss characteristics of these close-in planets.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Diverse Carbonates in Exoplanet Oceans Promote the Carbon Cycle

Carbonate precipitation in oceans is essential for the carbonate-silicate cycle (inorganic carbon cycle) to maintain temperate climates. By considering the thermodynamics of carbonate chemistry, we demonstrate that the ocean pH decreases by approximately 0.5 for a factor of 10 increase in the atmospheric carbon dioxide content. The upper and lower limits of ocean pH are within 1-4 of each other, where the upper limit is buffered by carbonate precipitation and defines the ocean pH when the carbon cycle operates. If the carbonate compensation depth (CCD) resides above the ocean floor, then carbonate precipitation and the carbon cycle cease to operate. The CCD is deep (>40 km) for high ocean temperature and high atmospheric carbon dioxide content. Key divalent carbonates of magnesium, calcium and iron produce an increasingly wider parameter space of deep CCDs, suggesting that chemical diversity promotes the carbon cycle. The search for life from exoplanets will benefit by including chemically more diverse targets than Earth twins.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A pair of Sub-Neptunes transiting the bright K-dwarf TOI-1064 characterised with CHEOPS

We report the discovery and characterisation of a pair of sub-Neptunes transiting the bright K-dwarf TOI-1064 (TIC 79748331), initially detected in TESS photometry. To characterise the system, we performed and retrieved CHEOPS, TESS, and ground-based photometry, HARPS high-resolution spectroscopy, and Gemini speckle imaging. We characterise the host star and determine $T_{\rm eff, \star}=4734\pm67$ K, $R_{\star}=0.726\pm0.007$ $R_{\odot}$, and $M_{\star}=0.748\pm0.032$ $M_{\odot}$. We present a novel detrending method based on PSF shape-change modelling and demonstrate its suitability to correct flux variations in CHEOPS data. We confirm the planetary nature of both bodies and find that TOI-1064 b has an orbital period of $P_{\rm b}=6.44387\pm0.00003$ d, a radius of $R_{\rm b}=2.59\pm0.04$ $R_{\oplus}$, and a mass of $M_{\rm b}=13.5_{-1.8}^{+1.7}$ $M_{\oplus}$, whilst TOI-1064 c has an orbital period of $P_{\rm c}=12.22657^{+0.00005}_{-0.00004}$ d, a radius of $R_{\rm c}=2.65\pm0.04$ $R_{\oplus}$, and a 3$σ$ upper mass limit of 8.5 ${\rm M_{\oplus}}$. From the high-precision photometry we obtain radius uncertainties of $\sim$1.6%, allowing us to conduct internal structure and atmospheric escape modelling. TOI-1064 b is one of the densest, well-characterised sub-Neptunes, with a tenuous atmosphere that can be explained by the loss of a primordial envelope following migration through the protoplanetary disc. It is likely that TOI-1064 c has an extended atmosphere due to the tentative low density, however further RVs are needed to confirm this scenario and the similar radii, different masses nature of this system. The high-precision data and modelling of TOI-1064 b are important for planets in this region of mass-radius space, and it allows us to identify a trend in bulk density-stellar metallicity for massive sub-Neptunes that may hint at the formation of this population of planets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Astrometric Accelerations as Dynamical Beacons: Discovery and Characterization of HIP 21152 B, the First T-Dwarf Companion in the Hyades

Benchmark brown dwarf companions with well-determined ages and model-independent masses are powerful tools to test substellar evolutionary models and probe the formation of giant planets and brown dwarfs. Here, we report the independent discovery of HIP~21152~B, the first imaged brown dwarf companion in the Hyades, and conduct a comprehensive orbital and atmospheric characterization of the system. HIP~21152 was targeted in an ongoing high-contrast imaging campaign of stars exhibiting proper motion changes between Hipparcos and Gaia, and was also recently identified by Bonavita et al. (2022) and Kuzuhara et al. (2022). Our Keck/NIRC2 and SCExAO/CHARIS imaging of HIP~21152 revealed a comoving companion at a separation of $0.37^{\prime\prime}$ (16 au). We perform a joint orbit fit of all available relative astrometry and radial velocities together with the Hipparcos-Gaia proper motions, yielding a dynamical mass of $24^{+6}_{-4}\,\mathrm{M_{Jup}}$, which is $1{-}2σ$ lower than evolutionary model predictions. Hybrid grids that include the evolution of cloud properties best reproduce the dynamical mass. We also identify a comoving wide-separation ($1837^{\prime\prime}$ or $7.9 \times 10^4 \, \mathrm{au}$) early-L dwarf with an inferred mass near the hydrogen-burning limit. Finally, we analyze the spectra and photometry of HIP~21152~B using the Saumon & Marley (2008) atmospheric models and a suite of retrievals. The best-fit grid-based models have $f_{\mathrm{sed}}=2$, indicating the presence of clouds, $T_{\mathrm{eff}}=1400 \, \mathrm{K}$, and $\log{g}=4.5 \, \mathrm{dex}$. These results are consistent with the object's spectral type of $\mathrm{T0\pm1}$. As the first benchmark brown dwarf companion in the Hyades, HIP~21152~B joins the small but growing number of substellar companions with well-determined ages and dynamical masses.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Exoplanet atmospheres at high resolution through a modest-size telescope. Fe II in MASCARA-2b and KELT-9b with FIES on the Nordic Optical Telescope

Ground-based, high-resolution spectrographs are providing us with an unprecedented view of the dynamics and chemistry of the atmospheres of planets outside the Solar System. While there is a large number of stable and precise high-resolution spectrographs on modest-size telescopes, it is the spectrographs at observatories with apertures larger than 3.5 metres that dominate the atmospheric follow-up of exoplanets. In this work, we explore the potential of characterising exoplanetary atmospheres with FIES, a high-resolution spectrograph at the 2.56 metre Nordic Optical Telescope. We observed two transits of MASCARA-2 b (also known as KELT-20 b) and one transit of KELT-9 b to search for atomic iron, a species that has been recently discovered in both neutral and ionised forms in the atmospheres of these ultra-hot Jupiters using large telescopes. Using a cross-correlation method, we detect a signal of Fe II at the $4.5σ$ and $4.0σ$ level in the transits of MASCARA-2 b. We also detect Fe II in the transit of KELT-9 b at the $8.5σ$ level. Although we do not find any significant Doppler shift in the signal of MASCARA-2 b, we do measure a moderate blueshift (3-6 km/s) of the feature in KELT-9 b, which might be a manifestation of high-velocity winds transporting Fe II from the planetary dayside to the nightside. Our work demonstrates the feasibility of investigating exoplanet atmospheres with FIES, potentially unlocking a wealth of additional atmosphere detections with this and other high-resolution spectrographs mounted on similar-size telescopes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

How do we optimally sample model grids of exoplanet spectra?

The construction and implementation of atmospheric model grids is a popular tool in exoplanet characterisation. These typically vary a number of parameters linearly, containing one model for every combination of parameter values. Here we investigate alternative methods of sampling parameters, including random sampling and Latin hypercube (LH) sampling, and how these compare to linearly sampled grids. We use a random forest to analyse the performance of these grids for two different models, as well as investigate the information content of the particular model grid from Goyal et al. 2019. We also use nested-sampling to implement mock atmospheric retrievals on simulated JWST transmission spectra by interpolating on linearly sampled model grids. Our results show that random or LH sampling out-performs linear sampling in parameter predictability for our higher dimensional models, requiring fewer models in the grid, and thus allowing for more computationally intensive forward models to be used. We also find that using a traditional retrieval with interpolation on a linear grid can produce biased posterior distributions, especially for parameters with non-linear effects on the spectrum. In particular, we advise caution when performing linear interpolation on the C/O ratio, cloud properties, and metallicity. Finally, we find that the information content analysis of the grid from Goyal et al. 2019 is able to highlight key areas of the spectra where the presence or absence of certain molecules can be detected, providing good indicators for parameters such as temperature and C/O ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Meta-modelling the climate of dry tide locked rocky planets

Rocky planets hosted by close-in extrasolar systems are likely to be tidally locked in 1:1 spin-orbit resonance, a configuration where they exhibit permanent dayside and nightside. Because of the resulting day-night temperature gradient, the climate and large-scale circulation of these planets are strongly determined by their atmospheric stability against collapse, which designates the runaway condensation of greenhouse gases. To better constrain the surface conditions of rocky planets located in the habitable zone of their host star, it is therefore crucial to elucidate the mechanisms that govern the day-night heat redistribution. As a first attempt to bridge the gap between multiple modelling approaches ranging from idealised models to 3-D General Circulation Models (GCM), we developed a General Circulation Meta-Model (GCMM) able to reproduce both the closed-form solutions provided by analytical models and the numerical solutions obtained from GCM simulations. We used this approach to characterise the atmospheric stability of Earth-sized rocky planets with dry atmospheres containing CO2, and we benchmarked it against 3-D GCM simulations using THOR GCM. We observe that the collapse pressure below which collapse occurs can vary by ~40% around the value predicted by analytical scaling laws depending on the mechanisms taken into account among radiative transfer, atmospheric dynamics, and turbulent diffusion. Particularly, we find (i) that the turbulent diffusion taking place in the dayside planetary boundary layer (PBL) globally tends to warm up the nightside surface hemisphere except in the transition zone between optically thin and optically thick regimes, (ii) that the PBL also significantly affects the day-night advection timescale, and (iii) that the slow rotator approximation holds from the moment that the normalised equatorial Rossby deformation radius is greater than 2.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Retrieval Study of Brown Dwarfs Across the L-T Sequence

A large suite of 228 atmospheric retrievals is performed on a curated sample of 19 brown dwarfs spanning the L0 to T8 spectral types using the open-source Helios-r2 retrieval code, which implements the method of short characteristics for radiative transfer and a finite-element description of the temperature-pressure profile. Surprisingly, we find that cloud-free and cloudy (both gray and non-gray) models are equally consistent with the archival SpeX data from the perspective of Bayesian model comparison. Only upper limits for cloud properties are inferred if log-uniform priors are assumed, but the cloud optical depth becomes constrained if a uniform prior is used. Water is detected in all 19 objects and methane is detected in all of the T dwarfs, but no obvious trend exists across effective temperature. As carbon monoxide is only detected in a handful of objects, the inferred carbon-to-oxygen ratios are unreliable. The retrieved radius generally decreases with effective temperature, but the values inferred for some T dwarfs are implausibly low and may indicate missing physics or chemistry in the models. For the early L dwarfs, the retrieved surface gravity depends on whether the gray or non-gray cloud model is preferred. Future data are necessary for constraining cloud properties and the vertical variation of chemical abundances, the latter of which is needed for distinguishing between the chemical instability versus traditional cloud interpretation of the L-T transition.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The THOR+HELIOS general circulation model: multi-wavelength radiative transfer with accurate scattering by clouds/hazes

General circulation models (GCMs) provide context for interpreting multi-wavelength, multi-phase data of the atmospheres of tidally locked exoplanets. In the current study, the non-hydrostatic THOR GCM is coupled with the HELIOS radiative transfer solver for the first time, supported by an equilibrium chemistry solver (FastChem), opacity calculator (HELIOS-K) and Mie scattering code (LX-MIE). To accurately treat the scattering of radiation by medium-sized to large aerosols/condensates, improved two-stream radiative transfer is implemented within a GCM for the first time. Multiple scattering is implemented using a Thomas algorithm formulation of the two-stream flux solutions, which decreases the computational time by about 2 orders of magnitude compared to the iterative method used in past versions of HELIOS. As a case study, we present four GCMs of the hot Jupiter WASP-43b, where we compare the temperature, velocity, entropy, and streamfunction, as well as the synthetic spectra and phase curves, of runs using regular versus improved two-stream radiative transfer and isothermal versus non-isothermal layers. While the global climate is qualitatively robust, the synthetic spectra and phase curves are sensitive to these details. A THOR+HELIOS WASP-43b GCM (horizontal resolution of about 4 degrees on the sphere and with 40 radial points) with multi-wavelength radiative transfer (30 k-table bins) running for 3000 Earth days (864,000 time steps) takes about 19-26 days to complete depending on the type of GPU.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Titanium oxide and chemical inhomogeneity in the atmosphere of the exoplanet WASP-189b

The temperature of an atmosphere decreases with increasing altitude, unless a shortwave absorber exists that causes a temperature inversion. Ozone plays this role in the Earth`s atmosphere. In the atmospheres of highly irradiated exoplanets, shortwave absorbers are predicted to be titanium oxide (TiO) and vanadium oxide (VO). Detections of TiO and VO have been claimed using both low and high spectral resolution observations, but later observations have failed to confirm these claims or overturned them. Here we report the unambiguous detection of TiO in the ultra-hot Jupiter WASP-189b using high-resolution transmission spectroscopy. This detection is based on applying the cross-correlation technique to many spectral lines of TiO from 460 to 690 nm. Moreover, we report detections of metals, including neutral and singly ionised iron and titanium, as well as chromium, magnesium, vanadium and manganese (Fe, Fe+, Ti, Ti+, Cr, Mg, V, Mn). The line positions of the detected species differ, which we interpret as a consequence of spatial gradients in their chemical abundances, such that they exist in different regions or dynamical regimes. This is direct observational evidence for the three-dimensional thermo-chemical stratification of an exoplanet atmosphere derived from high-resolution ground-based spectroscopy.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Uncovering the true periods of the young sub-Neptunes orbiting TOI-2076

Context: TOI-2076 is a transiting three-planet system of sub-Neptunes orbiting a bright (G = 8.9 mag), young ($340\pm80$ Myr) K-type star. Although a validated planetary system, the orbits of the two outer planets were unconstrained as only two non-consecutive transits were seen in TESS photometry. This left 11 and 7 possible period aliases for each. Aims: To reveal the true orbits of these two long-period planets, precise photometry targeted on the highest-probability period aliases is required. Long-term monitoring of transits in multi-planet systems can also help constrain planetary masses through TTV measurements. Methods: We used the MonoTools package to determine which aliases to follow, and then performed space-based and ground-based photometric follow-up of TOI-2076 c and d with CHEOPS, SAINT-EX, and LCO telescopes. Results: CHEOPS observations revealed a clear detection for TOI-2076 c at $P=21.01538^{+0.00084}_{-0.00074}$ d, and allowed us to rule out three of the most likely period aliases for TOI-2076 d. Ground-based photometry further enabled us to rule out remaining aliases and confirm the $P=35.12537\pm0.00067$ d alias. These observations also improved the radius precision of all three sub-Neptunes to $2.518\pm0.036$, $3.497\pm0.043$, and $3.232\pm0.063$ $R_\oplus$. Our observations also revealed a clear anti-correlated TTV signal between planets b and c likely caused by their proximity to the 2:1 resonance, while planets c and d appear close to a 5:3 period commensurability, although model degeneracy meant we were unable to retrieve robust TTV masses. Their inflated radii, likely due to extended H-He atmospheres, combined with low insolation makes all three planets excellent candidates for future comparative transmission spectroscopy with JWST.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Comparative Study of Atmospheric Chemistry with VULCAN

We present an update of the open-source photochemical kinetics code VULCAN (Tsai et al. 2017; https://github.com/exoclime/VULCAN) to include C-H-N-O-S networks and photochemistry. Additional new features are advection transport, condensation, various boundary conditions, and temperature-dependent UV cross-sections. First, we validate our photochemical model for hot Jupiter atmospheres by performing an intercomparison of HD 189733b models between Moses et al. (2011), Venot et al. (2012), and VULCAN, to diagnose possible sources of discrepancy. Second, we set up a model of Jupiter extending from the deep troposphere to upper stratosphere to verify the kinetics for low temperature. Our model reproduces hydrocarbons consistent with observations, and the condensation scheme successfully predicts the locations of water and ammonia ice clouds. We show that vertical advection can regulate the local ammonia distribution in the deep atmosphere. Third, we validate the model for oxidizing atmospheres by simulating Earth and find agreement with observations. Last, VULCAN is applied to four representative cases of extrasolar giant planets: WASP-33b, HD 189733b, GJ 436b, and 51 Eridani b. We look into the effects of the C/O ratio and chemistry of titanium/vanadium species for WASP-33b; we revisit HD 189733b for the effects of sulfur and carbon condensation; the effects of internal heating and vertical mixing ($K_{\textrm{zz}}$) are explored for GJ 436b; we test updated planetary properties for 51 Eridani b with S$_8$ condensates. We find sulfur can couple to carbon or nitrogen and impact other species such as hydrogen, methane, and ammonia. The observable features of the synthetic spectra and trends in the photochemical haze precursors are discussed for each case.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Physically-motivated basis functions for temperature maps of exoplanets

Thermal phase curves of exoplanet atmospheres have revealed temperature maps as a function of planetary longitude, often by sinusoidal decomposition of the phase curve. We construct a framework for describing two-dimensional temperature maps of exoplanets with mathematical basis functions derived for a fluid layer on a rotating, heated sphere with drag/friction, which are generalizations of spherical harmonics. These basis functions naturally produce physically-motivated temperature maps for exoplanets with few free parameters. We investigate best practices for applying this framework to temperature maps of hot Jupiters by splitting the problem into two parts: (1) we constrain the temperature map as a function of latitude by tuning the basis functions to reproduce general circulation model (GCM) outputs, since disk-integrated phase curve observations do not constrain this dimension; and (2) we infer the temperature maps of real hot Jupiters using original reductions of several Spitzer phase curves, which directly constrain the temperature variations with longitude. The resulting phase curves can be described with only three free parameters per bandpass -- an efficiency improvement over the usual five or so used to describe sinusoidal decompositions of phase curves. Upon obtaining the hemispherically averaged dayside and nightside temperatures, the standard approach would be to use zero-dimensional box models to infer the Bond albedo and redistribution efficiency. We elucidate the limitation of these box models by demonstrating that negative Bond albedos may be obtained due to a choice of boundary condition on the nightside temperature. We propose generalized definitions for the Bond albedo and heat redistribution efficiency for use with two-dimensional (2D) temperature maps. Open-source software called kelp is provided to efficiently compute these phase curves.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Giant Planet Candidate Transiting a White Dwarf

Astronomers have discovered thousands of planets outside the solar system, most of which orbit stars that will eventually evolve into red giants and then into white dwarfs. During the red giant phase, any close-orbiting planets will be engulfed by the star, but more distant planets can survive this phase and remain in orbit around the white dwarf. Some white dwarfs show evidence for rocky material floating in their atmospheres, in warm debris disks, or orbiting very closely, which has been interpreted as the debris of rocky planets that were scattered inward and tidally disrupted. Recently, the discovery of a gaseous debris disk with a composition similar to ice giant planets demonstrated that massive planets might also find their way into tight orbits around white dwarfs, but it is unclear whether the planets can survive the journey. So far, the detection of intact planets in close orbits around white dwarfs has remained elusive. Here, we report the discovery of a giant planet candidate transiting the white dwarf WD 1856+534 (TIC 267574918) every 1.4 days. The planet candidate is roughly the same size as Jupiter and is no more than 14 times as massive (with 95% confidence). Other cases of white dwarfs with close brown dwarf or stellar companions are explained as the consequence of common-envelope evolution, wherein the original orbit is enveloped during the red-giant phase and shrinks due to friction. In this case, though, the low mass and relatively long orbital period of the planet candidate make common-envelope evolution less likely. Instead, the WD 1856+534 system seems to demonstrate that giant planets can be scattered into tight orbits without being tidally disrupted, and motivates searches for smaller transiting planets around white dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Atmospheric stability and collapse on tidally locked rocky planets

Over large timescales, a terrestrial planet may be driven towards spin-orbit synchronous rotation by tidal forces. In this particular configuration, the planet exhibits permanent dayside and nightside, which may induce strong day-night temperature gradients. The nightside temperature depends on the efficiency of the day-night heat redistribution and determines the stability of the atmosphere against collapse. To better constrain the atmospheric stability, climate, and surface conditions of rocky planets located in the habitable zone of their host star, it is thus crucial to understand the complex mechanism of heat redistribution. Building on early works and assuming dry thermodynamics, we developed a hierarchy of analytic models taking into account the coupling between radiative transfer, dayside convection, and large-scale atmospheric circulation in the case of slowly rotating planets. There are two types of these models: a zero-dimensional two-layer approach and a two-column radiative-convective-subsiding-upwelling (RCSU) model. They yield analytical solutions and scaling laws characterising the dependence of the collapse pressure on physical features, which are compared to the results obtained by early works using 3D global climate models (GCMs). The analytical theory captures (i) the dependence of temperatures on atmospheric opacities and scattering in the shortwave and in the longwave, (ii) the behaviour of the collapse pressure observed in GCM simulations at low stellar fluxes that are due to the non-linear dependence of the atmospheric opacity on the longwave optical depth at the planet's surface, (iii) the increase of stability generated by dayside sensible heating, and (iv) the decrease of stability induced by the increase of the planet size.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Exploring the atmospheric dynamics of the extreme ultra-hot Jupiter KELT-9b using TESS photometry

We carry out a phase-curve analysis of the KELT-9 system using photometric observations from NASA's Transiting Exoplanet Survey Satellite (TESS). The measured secondary eclipse depth and peak-to-peak atmospheric brightness modulation are $650^{+14}_{-15}$ ppm and $566\pm16$ ppm, respectively. The planet's brightness variation reaches maximum $31\pm5$ minutes before the midpoint of the secondary eclipse, indicating a $5\overset{\circ}{.}2\pm0\overset{\circ}{.}9$ eastward shift in the dayside hot spot from the substellar point. We also detect stellar pulsations on KELT-9 with a period of $7.58695\pm0.00091$ hours. The dayside emission of KELT-9b in the TESS bandpass is consistent with a blackbody brightness temperature of $4600\pm100$ K. The corresponding nightside brightness temperature is $3040\pm100$ K, comparable to the dayside temperatures of the hottest known exoplanets. In addition, we detect a significant phase-curve signal at the first harmonic of the orbital frequency and a marginal signal at the second harmonic. While the amplitude of the first harmonic component is consistent with the predicted ellipsoidal distortion modulation assuming equilibrium tides, the phase of this photometric variation is shifted relative to the expectation. Placing KELT-9b in the context of other exoplanets with phase-curve observations, we find that the elevated nightside temperature and relatively low day-night temperature contrast agree with the predictions of atmospheric models that include H$_{2}$ dissociation and recombination. The nightside temperature of KELT-9b implies an atmospheric composition containing about 50% molecular and 50% atomic hydrogen at 0.1 bar, a nightside emission spectrum that deviates significantly from a blackbody, and a 0.5-2.0 $μ$m transmission spectrum that is featureless at low resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Helios-r2 -- A new Bayesian, open-source retrieval model for brown dwarfs and exoplanet atmospheres

We present an improved, hybrid CPU-GPU atmospheric retrieval code, Helios-r2, which is applicable to medium-resolution emission spectra of brown dwarfs, in preparation for precision atmospheric spectroscopy in the era of the James Webb Space Telescope. The model is available as open-source code on the Exoclimes Simulation Platform. We subject Helios-r2 to a battery of tests of varying difficulty. The simplest test involves a mock retrieval on a forward model generated using the same radiative transfer technique, the same implementation of opacities, and the same chemistry model. The least trivial test involves a mock retrieval on synthetic spectra from the Sonora model grid, which uses a different radiative transfer technique, a different implementation of opacities, and a different chemistry model. A calibration factor, which is included to capture uncertainties in the brown dwarf radius, distance to the brown dwarf and flux calibration of the spectrum, may compensate, sometimes erroneously, for discrepancies in modeling choices and implementation. We analyze spectra of the benchmark brown dwarf GJ 570 D and the binary brown dwarf companions in the Epsilon Indi system. The retrieved surface gravities are consistent with previous studies and/or values inferred from dynamical masses (for Epsilon Indi Ba and Bb only). There remains no clear criterion on how to reject unphysical values of the retrieved brown dwarf radii. The inferred radii and corresponding masses should be taken with great caution. The retrieved carbon-to-oxygen ratios and metallicity depend on whether chemical equilibrium is assumed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Information content of JWST-NIRSPEC transmission spectra of warm Neptunes

Warm Neptunes offer a rich opportunity for understanding exo-atmospheric chemistry. With the upcoming James Webb Space Telescope (JWST), there is a need to elucidate the balance between investments in telescope time versus scientific yield. We use the supervised machine learning method of the random forest to perform an information content analysis on a 11-parameter model of transmission spectra from the various NIRSpec modes. The three bluest medium-resolution NIRSpec modes (0.7 - 1.27 microns, 0.97 - 1.84 microns, 1.66 - 3.07 microns) are insensitive to the presence of CO. The reddest medium-resolution mode (2.87 - 5.10 microns) is sensitive to all of the molecules assumed in our model: CO, CO2, CH4, C2H2, H2O, HCN and NH3. It competes effectively with the three bluest modes on the information encoded on cloud abundance and particle size. It is also competitive with the low-resolution prism mode (0.6 - 5.3 microns) on the inference of every parameter except for the temperature and ammonia abundance. We recommend astronomers to use the reddest medium-resolution NIRSpec mode for studying the atmospheric chemistry of 800-1200 K warm Neptunes; its corresponding high-resolution counterpart offers diminishing returns. We compare our findings to previous JWST information content analyses that favor the blue orders, and suggest that the reliance on chemical equilibrium could lead to biased outcomes if this assumption does not apply. A simple, pressure-independent diagnostic for identifying chemical disequilibrium is proposed based on measuring the abundances of H2O, CO and CO2.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Interpreting High-Resolution Spectroscopy of Exoplanets Using Cross-Correlations and Supervised Machine Learning

We present a new method for performing atmospheric retrieval on ground-based, high-resolution data of exoplanets. Our method combines cross-correlation functions with a random forest, a supervised machine learning technique, to overcome challenges associated with high-resolution data. A series of cross-correlation functions are concatenated to give a "CCF-sequence" for each model atmosphere, which reduces the dimensionality by a factor of ~100. The random forest, trained on our grid of ~65,000 models, provides a likelihood-free method of retrieval. The pre-computed grid spans 31 values of both temperature and metallicity, and incorporates a realistic noise model. We apply our method to HARPS-N observations of the ultra-hot Jupiter KELT-9b, and obtain a metallicity consistent with solar (logM = $-0.2\pm0.2$). Our retrieved transit chord temperature (T = $6000^{+0}_{-200}$K) is unreliable as the ion cross-correlations lie outside of the training set, which we interpret as being indicative of missing physics in our atmospheric model. We compare our method to traditional nested-sampling, as well as other machine learning techniques, such as Bayesian neural networks. We demonstrate that the likelihood-free aspect of the random forest makes it more robust than nested-sampling to different error distributions, and that the Bayesian neural network we tested is unable to reproduce complex posteriors. We also address the claim in Cobb et al. (2019) that our random forest retrieval technique can be over-confident but incorrect. We show that this is an artefact of the training set, rather than the machine learning method, and that the posteriors agree with those obtained using nested-sampling.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Outstanding Challenges of Exoplanet Atmospheric Retrievals

Spectral retrieval has long been a powerful tool for interpreting planetary remote sensing observations. Flexible, parameterised, agnostic models are coupled with inversion algorithms in order to infer atmospheric properties directly from observations, with minimal reliance on physical assumptions. This approach, originally developed for application to Earth satellite data and subsequently observations of other Solar System planets, has been recently and successfully applied to transit, eclipse and phase curve spectra of transiting exoplanets. In this review, we present the current state-of-the-art in terms of our ability to accurately retrieve information about atmospheric chemistry, temperature, clouds and spatial variability; we discuss the limitations of this, both in the available data and modelling strategies used; and we recommend approaches for future improvement.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The lifetimes of planetary debris discs around white dwarfs

The lifetime of a planetary disc which orbits a white dwarf represents a crucial input parameter into evolutionary models of that system. Here we apply a purely analytical formalism to estimate lifetimes of the debris phase of these discs, before they are ground down into dust or are subject to sublimation from the white dwarf. We compute maximum lifetimes for three different types of white dwarf discs, formed from (i) radiative YORP breakup of exo-asteroids along the giant branch phases at 2-100 au, (ii) radiation-less spin-up disruption of these minor planets at $\sim 1.5-4.5R_{\odot}$, and (iii) tidal disruption of minor or major planets within about $1.3R_{\odot}$. We display these maximum lifetimes as a function of disc mass and extent, constituent planetesimal properties, and representative orbital excitations of eccentricity and inclination. We find that YORP discs with masses up to $10^{24}$ kg live long enough to provide a reservoir of surviving cm-sized pebbles and m- to km-sized boulders that can be perturbed intact to white dwarfs with cooling ages of up to 10 Gyr. Debris discs formed from the spin or tidal disruption of these minor planets or major planets can survive in a steady state for up to respectively 1 Myr or 0.01 Myr, although most tidal discs would leave a steady state within about 1 yr. Our results illustrate that dust-less planetesimal transit detections are plausible, and would provide particularly robust evolutionary constraints. Our formalism can easily be adapted to individual systems and future discoveries.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

THOR 2.0: Major Improvements to the Open-Source General Circulation Model

THOR is the first open-source general circulation model (GCM) developed from scratch to study the atmospheres and climates of exoplanets, free from Earth- or Solar System-centric tunings. It solves the general non-hydrostatic Euler equations (instead of the primitive equations) on a sphere using the icosahedral grid. In the current study, we report major upgrades to THOR, building upon the work of Mendonça et al. (2016). First, while the Horizontally Explicit Vertically Implicit (HEVI) integration scheme is the same as that described in Mendonça et al. (2016), we provide a clearer description of the scheme and improved its implementation in the code. The differences in implementation between the hydrostatic shallow (HSS), quasi-hydrostatic deep (QHD) and non-hydrostatic deep (NHD) treatments are fully detailed. Second, standard physics modules are added: two-stream, double-gray radiative transfer and dry convective adjustment. Third, THOR is tested on additional benchmarks: tidally-locked Earth, deep hot Jupiter, acoustic wave, and gravity wave. Fourth, we report that differences between the hydrostatic and non-hydrostatic simulations are negligible in the Earth case, but pronounced in the hot Jupiter case. Finally, the effects of the so-called "sponge layer", a form of drag implemented in most GCMs to provide numerical stability, are examined. Overall, these upgrades have improved the flexibility, user-friendliness, and stability of THOR.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Supervised Machine Learning for Intercomparison of Model Grids of Brown Dwarfs: Application to GJ 570D and the Epsilon Indi B Binary System

Self-consistent model grids of brown dwarfs involve complex physics and chemistry, and are often computed using proprietary computer codes, making it challenging to identify the reasons for discrepancies between model and data as well as between the models produced by different research groups. In the current study, we demonstrate a novel method for analyzing brown dwarf spectra, which combines the use of the Sonora, AMES-Cond and HELIOS model grids with the supervised machine learning method of the random forest. Besides performing atmospheric retrieval, the random forest enables information content analysis of the three model grids as a natural outcome of the method, both individually on each grid and by comparing the grids against one another, via computing large suites of mock retrievals. Our analysis reveals that the different choices made in modelling the alkali line shapes hinder the use of the alkali lines as gravity indicators. Nevertheless, the spectrum longward of 1.2 micron encodes enough information on the surface gravity to allow its inference from retrieval. Temperature may be accurately and precisely inferred independent of the choice of model grid, but not the surface gravity. We apply random forest retrieval to three objects: the benchmark T7.5 brown dwarf GJ 570D; and Epsilon Indi Ba (T1.5 brown dwarf) and Bb (T6 brown dwarf), which are part of a binary system and have measured dynamical masses. For GJ 570D, the inferred effective temperature and surface gravity are consistent with previous studies. For Epsilon Indi Ba and Bb, the inferred surface gravities are broadly consistent with the values informed by the dynamical masses.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Cloudiness Index for Transiting Exoplanets Based on the Sodium and Potassium Lines: Tentative Evidence for Hotter Atmospheres Being Less Cloudy at Visible Wavelengths

We present a dimensionless index that quantifies the degree of cloudiness of the atmosphere of a transiting exoplanet. Our cloudiness index is based on measuring the transit radii associated with the line center and wing of the sodium or potassium line. In deriving this index, we revisited the algebraic formulae for inferring the isothermal pressure scale height from transit measurements. We demonstrate that the formulae of Lecavelier et al. and Benneke & Seager are identical: the former is inferring the temperature while assuming a value for the mean molecular mass and the latter is inferring the mean molecular mass while assuming a value for the temperature. More importantly, these formulae cannot be used to distinguish between cloudy and cloudfree atmospheres. We derive values of our cloudiness index for a small sample of 7 hot Saturns/Jupiters taken from Sing et al. We show that WASP-17b, WASP-31b and HAT-P-1b are nearly cloudfree at visible wavelengths. We find the tentative trend that more irradiated atmospheres tend to have less clouds consisting of sub-micron-sized particles. We also derive absolute sodium and/or potassium abundances $\sim 10^2$ cm$^{-3}$ for WASP-17b, WASP-31b and HAT-P-1b (and upper limits for the other objects). Higher-resolution measurements of both the sodium and potassium lines, for a larger sample of exoplanetary atmospheres, are needed to confirm or refute this trend.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A map of the large day-night temperature gradient of a super-Earth exoplanet

Over the past decade, observations of giant exoplanets (Jupiter-size) have provided key insights into their atmospheres, but the properties of lower-mass exoplanets (sub-Neptune) remain largely unconstrained because of the challenges of observing small planets. Numerous efforts to observe the spectra of super-Earths (exoplanets with masses of one to ten times that of Earth) have so far revealed only featureless spectra. Here we report a longitudinal thermal brightness map of the nearby transiting super-Earth 55 Cancri e revealing highly asymmetric dayside thermal emission and a strong day-night temperature contrast. Dedicated space-based monitoring of the planet in the infrared revealed a modulation of the thermal flux as 55 Cancri e revolves around its star in a tidally locked configuration. These observations reveal a hot spot that is located 41 +- 12 degrees east of the substellar point (the point at which incident light from the star is perpendicular to the surface of the planet). From the orbital phase curve, we also constrain the nightside brightness temperature of the planet to 1380 +- 400 kelvin and the temperature of the warmest hemisphere (centred on the hot spot) to be about 1300 kelvin hotter (2700 +- 270 kelvin) at a wavelength of 4.5 microns, which indicates inefficient heat redistribution from the dayside to the nightside. Our observations are consistent with either an optically thick atmosphere with heat recirculation confined to the planetary dayside, or a planet devoid of atmosphere with low-viscosity magma flows at the surface.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optical phase curves as diagnostics for aerosol composition in exoplanetary atmospheres

Optical phase curves have become one of the common probes of exoplanetary atmospheres, but the information they encode has not been fully elucidated. Building on a diverse body of work, we upgrade the Flexible Modeling System (FMS) to include scattering in the two-stream, dual-band approximation and generate plausible, three-dimensional structures of irradiated atmospheres to study the radiative effects of aerosols or condensates. In the optical, we treat the scattering of starlight using a generalisation of Beer's law that allows for a finite Bond albedo to be prescribed. In the infrared, we implement the two-stream solutions and include scattering via an infrared scattering parameter. We present a suite of four-parameter general circulation models for Kepler-7b and demonstrate that its climatology is expected to be robust to variations in optical and infrared scattering. The westward and eastward shifts of the optical and infrared phase curves, respectively, are shown to be robust outcomes of the simulations. Assuming micron-sized particles and a simplified treatment of local brightness, we further show that the peak offset of the optical phase curve is sensitive to the composition of the aerosols or condensates. However, to within the measurement uncertainties, we cannot distinguish between aerosols made of silicates (enstatite or forsterite), iron, corundum or titanium oxide, based on a comparison to the measured peak offset ($41^\circ \pm 12^\circ$) of the optical phase curve of Kepler-7b. Measuring high-precision optical phase curves will provide important constraints on the atmospheres of cloudy exoplanets and reduce degeneracies in interpreting their infrared spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Raman Scattering by Molecular Hydrogen and Nitrogen in Exoplanetary Atmospheres

An important source of opacity in exoplanet atmospheres at short visible and near-UV wavelengths is Rayleigh scattering of light on molecules. It is accompanied by a related, albeit weaker process -- Raman scattering. We analyze the signatures of Raman scattering imprinted in the reflected light and the geometric albedo of exoplanets, which could provide information about atmospheric properties. Raman scattering affects the geometric albedo spectra of planets in following ways. Firstly, it causes filling-in of strong absorption lines in the incident radiation, thus producing sharp peaks in the albedo. Secondly, it shifts the wavelengths of spectral features in the reflected light causing the so-called Raman ghost lines. Raman scattering can also cause a broadband reduction of the albedo due to wavelength shifting of a stellar spectrum with red spectral index. Observing the Raman peaks in the albedo could be used to measure the column density of gas, thus providing constrains on the presence of clouds in the atmosphere. Observing the Raman ghost lines could be used to spectroscopically identify the main scatterer in the atmosphere, even molecules like H$_2$ or N$_2$ that do not have prominent spectral signatures in the optical wavelength range. If detected, ghost lines could also provide information about the temperature of the atmosphere. In this paper we investigate the effects of Raman scattering in hydrogen- and nitrogen-dominated atmospheres. We analyze the feasibility of detecting the signatures of Raman scattering with the existing and future observational facilities, and of using these signatures as probes of exoplanetary atmospheres.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Shear-driven instabilities and shocks in the atmospheres of hot Jupiters

General circulation models of the atmosphere of hot Jupiters have shown the existence of a supersonic eastward equatorial jet. In this paper, we investigate the effects of compressibility on the atmospheric dynamics by solving the standard Euler equations. This is done by means of a series of simulations performed in the framework of the equatorial beta-plane approximation using the finite volume shock-capturing code RAMSES. At low resolution, we recover the classical results described in the literature: we find a strong and steady supersonic equatorial jet of a few km/s that displays no signature of shocks. We next show that the jet zonal velocity depends significantly on the grid meridional resolution. When that resolution is fine enough to properly resolve the jet, the latter is subject to a Kelvin-Helmholtz instability. The jet zonal mean velocity displays regular oscillations with a typical timescale of few days and a significant amplitude of about 15% of the jet velocity. We also find compelling evidence for the development of a vertical shear instability at pressure levels of a few bars. It seems to be responsible for an increased downward kinetic energy flux, significantly affecting the temperature of the deep atmosphere, and appears to act as a form of drag on the equatorial jet. This instability also creates velocity fluctuations that propagate upward and steepen into weak shocks at pressure levels of a few mbars. We conclude that hot Jupiter equatorial jets are potentially unstable to both a barotropic Kelvin-Helmholtz instability and a vertical shear instability. Upon confirmation using more realistic models, both instabilities could result in significant time variability of the atmospheric winds, may provide a small scale dissipation mechanism in the flow, and might have consequences for the internal evolution of hot Jupiters.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Need for Laboratory Work to Aid in The Understanding of Exoplanetary Atmospheres

Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements stem from significant investments in new missions and facilities, such as JWST and the several planned ground-based extremely large telescopes. However, while exoplanet science currently has a wide range of sophisticated models that can be applied to the tide of forthcoming observations, the trajectory for preparing these models for the upcoming observational challenges is unclear. Thus, our ability to maximize the insights gained from the next generation of observatories is not certain. In many cases, uncertainties in a path towards model advancement stems from insufficiencies in the laboratory data that serve as critical inputs to atmospheric physical and chemical tools. We outline a number of areas where laboratory or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry. Specifically highlighted are needs for: (1) molecular opacity linelists with parameters for a diversity of broadening gases, (2) extended databases for collision-induced absorption and dimer opacities, (3) high spectral resolution opacity data for relevant molecular species, (4) laboratory studies of haze and condensate formation and optical properties, (5) significantly expanded databases of chemical reaction rates, and (6) measurements of gas photo-absorption cross sections at high temperatures. We hope that by meeting these needs, we can make the next two decades of exoplanet science as productive and insightful as the previous two decades. (abr)

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Transiting Exoplanet Studies and Community Targets for JWST's Early Release Science Program

The James Webb Space Telescope will revolutionize transiting exoplanet atmospheric science due to its capability for continuous, long-duration observations and its larger collecting area, spectral coverage, and spectral resolution compared to existing space-based facilities. However, it is unclear precisely how well JWST will perform and which of its myriad instruments and observing modes will be best suited for transiting exoplanet studies. In this article, we describe a prefatory JWST Early Release Science (ERS) program that focuses on testing specific observing modes to quickly give the community the data and experience it needs to plan more efficient and successful future transiting exoplanet characterization programs. We propose a multi-pronged approach wherein one aspect of the program focuses on observing transits of a single target with all of the recommended observing modes to identify and understand potential systematics, compare transmission spectra at overlapping and neighboring wavelength regions, confirm throughputs, and determine overall performances. In our search for transiting exoplanets that are well suited to achieving these goals, we identify 12 objects (dubbed "community targets") that meet our defined criteria. Currently, the most favorable target is WASP-62b because of its large predicted signal size, relatively bright host star, and location in JWST's continuous viewing zone. Since most of the community targets do not have well-characterized atmospheres, we recommend initiating preparatory observing programs to determine the presence of obscuring clouds/hazes within their atmospheres. Measurable spectroscopic features are needed to establish the optimal resolution and wavelength regions for exoplanet characterization. Other initiatives from our proposed ERS program include testing the instrument brightness limits and performing phase-curve observations.(Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Non-isothermal Theory for Interpreting Sodium Lines in Transmission Spectra of Exoplanets

We present a theory for interpreting the sodium lines detected in transmission spectra of exoplanetary atmospheres. Previous analyses employed the isothermal approximation and dealt only with the transit radius. By recognising the absorption depth and the transit radius as being independent observables, we develop a theory for jointly interpreting both quantities, which allows us to infer the temperatures and number densities associated with the sodium lines. We are able to treat a non-isothermal situation with a constant temperature gradient. Our novel diagnostics take the form of simple-to-use algebraic formulae and require measurements of the transit radii (and their corresponding absorption depths) at line center and in the line wing for both sodium lines. We apply our diagnostics to the HARPS data of HD 189733b, confirm the upper atmospheric heating reported by Huitson et al. (2012), derive a temperature gradient of $0.4376 \pm 0.0154$ K km$^{-1}$ and find densities $\sim 1$ to $10^4$ cm$^{-3}$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An Integral View of Fast Shocks around Supernova 1006

Supernova remnants are among the most spectacular examples of astrophysical pistons in our cosmic neighborhood. The gas expelled by the supernova explosion is launched with velocities ~ 1000 km/s into the ambient, tenuous interstellar medium, producing shocks that excite hydrogen lines. We have used an optical integral-field spectrograph to obtain high-resolution spatial-spectral maps that allow us to study in detail the shocks in the northwestern rim of supernova 1006. The two-component Hα line is detected at 133 sky locations. Variations in the broad line widths and the broad-to-narrow line intensity ratios across tens of atomic mean free paths suggest the presence of suprathermal protons, the potential seed particles for generating high-energy cosmic-rays.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Atmospheric Chemistry for Astrophysicists: A Self-consistent Formalism and Analytical Solutions for Arbitrary C/O

We present a self-consistent formalism for computing and understanding the atmospheric chemistry of exoplanets from the viewpoint of an astrophysicist. Starting from the first law of thermodynamics, we demonstrate that the van't Hoff equation (which describes the equilibrium constant), Arrhenius equation (which describes the rate coefficients) and procedures associated with the Gibbs free energy (minimisation, rescaling) have a common physical and mathematical origin. We address an ambiguity associated with the equilibrium constant, which is used to relate the forward and reverse rate coefficients, and restate its two definitions. By necessity, one of the equilibrium constants must be dimensionless and equate to an exponential function involving the Gibbs free energy, while the other is a ratio of rate coefficients and must therefore possess physical units. We demonstrate that the Arrhenius equation takes on a functional form that is more general than previously stated without recourse to tagging on ad hoc functional forms. Finally, we derive analytical models of chemical systems, in equilibrium, with carbon, hydrogen and oxygen. We include acetylene and are able to reproduce several key trends, versus temperature and carbon-to-oxygen ratio, published in the literature. The rich variety of behavior that mixing ratios exhibit as a function of the carbon-to-oxygen ratio is merely the outcome of stoichiometric book-keeping and not the direct consequence of temperature or pressure variations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Can we constrain interior structure of rocky exoplanets from mass and radius measurements?

We present an inversion method based on Bayesian analysis to constrain the interior structure of terrestrial exoplanets, in the form of chemical composition of the mantle and core size. Specifically, we identify what parts of the interior structure of terrestrial exoplanets can be determined from observations of mass, radius, and stellar elemental abundances. We perform a full probabilistic inverse analysis to formally account for observational and model uncertainties and obtain confidence regions of interior structure models. This enables us to characterize how model variability depends on data and associated uncertainties. We test our method on terrestrial solar system planets and find that our model predictions are consistent with independent estimates. Furthermore, we apply our method to synthetic exoplanets up to 10 Earth masses and up to 1.7 Earth radii as well as to exoplanet Kepler-36b. Importantly, the inversion strategy proposed here provides a framework for understanding the level of precision required to characterize the interior of exoplanets. Our main conclusions are: (1) observations of mass and radius are sufficient to constrain core size; (2) stellar elemental abundances (Fe, Si, Mg) are key constraints to reduce degeneracy in interior structure models and to constrain mantle composition; (3) the inherent degeneracy in determining interior structure from mass and radius observations does not only depend on measurement accuracies but also on the actual size and density of the exoplanet. We argue that precise observations of stellar elemental abundances are central in order to place constraints on planetary bulk composition and to reduce model degeneracy. [...]

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Carbon Dioxide in Exoplanetary Atmospheres: Rarely Dominant Compared to Carbon Monoxide and Water in Hot, Hydrogen-dominated Atmospheres

We present a comprehensive study of the abundance of carbon dioxide in exoplanetary atmospheres in hot, hydrogen-dominated atmospheres. We construct novel analytical models of systems in chemical equilibrium that include carbon monoxide, carbon dioxide, water, methane and acetylene and relate the equilibrium constants of the chemical reactions to temperature and pressure via the tabulated Gibbs free energies. We prove that such chemical systems may be described by a quintic equation for the mixing ratio of methane. By examining the abundances of these molecules across a broad range of temperatures (spanning equilibrium temperatures from 600 to 2500 K), pressures (via temperature-pressure profiles that explore albedo and opacity variations) and carbon-to-oxygen ratios, we conclude that carbon dioxide is subdominant compared to carbon monoxide and water. Atmospheric mixing does not alter this conclusion if carbon dioxide is subdominant everywhere in the atmosphere. Carbon dioxide and carbon monoxide may attain comparable abundances if the metallicity is greatly enhanced, but this property is negated by temperatures above 1000 K. For hydrogen-dominated atmospheres, our generic result has the implication that retrieval studies may wish to set the subdominance of carbon dioxide as a prior of the calculation and not let its abundance completely roam free as a fitting parameter, because it directly affects the inferred value of the carbon-to-oxygen ratio and may produce unphysical conclusions. We discuss the relevance of these implications for the hot Jupiter WASP-12b and suggest that some of the previous results are chemically impossible. The relative abundance of carbon dioxide to acetylene is potentially a sensitive diagnostic of the carbon-to-oxygen ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

HELIOS-K: An Ultrafast, Open-source Opacity Calculator for Radiative Transfer

We present an ultrafast opacity calculator that we name HELIOS-K. It takes a line list as an input, computes the shape of each spectral line and provides an option for grouping an enormous number of lines into a manageable number of bins. We implement a combination of Algorithm 916 and Gauss-Hermite quadrature to compute the Voigt profile, write the code in CUDA and optimise the computation for graphics processing units (GPUs). We restate the theory of the k-distribution method and use it to reduce $\sim 10^5$ to $10^8$ lines to $\sim 10$ to $10^4$ wavenumber bins, which may then be used for radiative transfer, atmospheric retrieval and general circulation models. The choice of line-wing cutoff for the Voigt profile is a significant source of error and affects the value of the computed flux by $\sim 10\%$. This is an outstanding physical (rather than computational) problem, due to our incomplete knowledge of pressure broadening of spectral lines in the far line wings. We emphasize that this problem remains regardless of whether one performs line-by-line calculations or uses the k-distribution method and affects all calculations of exoplanetary atmospheres requiring the use of wavelength-dependent opacities. We elucidate the correlated-k approximation and demonstrate that it applies equally to inhomogeneous atmospheres with a single atomic/molecular species or homogeneous atmospheres with multiple species. Using a NVIDIA K20 GPU, HELIOS-K is capable of computing an opacity function with $\sim 10^5$ spectral lines in $\sim 1$ second and is publicly available as part of the Exoclimes Simulation Platform (ESP; www.exoclime.org).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mapping High-velocity H-alpha and Lyman-alpha Emission from Supernova 1987A

We present new {\it Hubble Space Telescope} images of high-velocity H-$α$ and Lyman-$α$ emission in the outer debris of SN~1987A. The H-$α$ images are dominated by emission from hydrogen atoms crossing the reverse shock. For the first time we observe emission from the reverse shock surface well above and below the equatorial ring, suggesting a bipolar or conical structure perpendicular to the ring plane. Using the H$α$ imaging, we measure the mass flux of hydrogen atoms crossing the reverse shock front, in the velocity intervals ($-$7,500~$<$~$V_{obs}$~$<$~$-$2,800 km s$^{-1}$) and (1,000~$<$~$V_{obs}$~$<$~7,500 km s$^{-1}$), $\dot{M_{H}}$ = 1.2~$\times$~10$^{-3}$ M$_{\odot}$ yr$^{-1}$. We also present the first Lyman-$α$ imaging of the whole remnant and new $Chandra$ X-ray observations. Comparing the spatial distribution of the Lyman-$α$ and X-ray emission, we observe that the majority of the high-velocity Lyman-$α$ emission originates interior to the equatorial ring. The observed Lyman-$α$/H-$α$ photon ratio, $\langle$$R(Lα/ Hα)$$\rangle$ $\approx$~17, is significantly higher than the theoretically predicted ratio of $\approx$ 5 for neutral atoms crossing the reverse shock front. We attribute this excess to Lyman-$α$ emission produced by X-ray heating of the outer debris. The spatial orientation of the Lyman-$α$ and X-ray emission suggests that X-ray heating of the outer debris is the dominant Lyman-$α$ production mechanism in SN 1987A at this phase in its evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Destruction of the Circumstellar Ring of SN 1987A

We present imaging and spectroscopic observations with HST and VLT of the ring of SN 1987A from 1994 to 2014. After an almost exponential increase of the shocked emission from the hotspots up to day ~8,000 (~2009), both this and the unshocked emission are now fading. From the radial positions of the hotspots we see an acceleration of these up to 500-1000 km/s, consistent with the highest spectroscopic shock velocities from the radiative shocks. In the most recent observations (2013 and 2014), we find several new hotspots outside the inner ring, excited by either X-rays from the shocks or by direct shock interaction. All of these observations indicate that the interaction with the supernova ejecta is now gradually dissolving the hotspots. We predict, based on the observed decay, that the inner ring will be destroyed by ~2025.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Analytical Models of Exoplanetary Atmospheres. I. Atmospheric Dynamics via the Shallow Water System

Within the context of exoplanetary atmospheres, we present a comprehensive linear analysis of forced, damped, magnetized shallow water systems, exploring the effects of dimensionality, geometry (Cartesian, pseudo-spherical and spherical), rotation, magnetic tension and hydrodynamic and magnetic sources of friction. Across a broad range of conditions, we find that the key governing equation for atmospheres and quantum harmonic oscillators are identical, even when forcing (stellar irradiation), sources of friction (molecular viscosity, Rayleigh drag and magnetic drag) and magnetic tension are included. The global atmospheric structure is largely controlled by a single, key parameter that involves the Rossby and Prandtl numbers. This near-universality breaks down when either molecular viscosity or magnetic drag acts non-uniformly across latitude or a poloidal magnetic field is present, suggesting that these effects will introduce qualitative changes to the familiar chevron-shaped feature witnessed in simulations of atmospheric circulation. We also find that hydrodynamic and magnetic sources of friction have dissimilar phase signatures and affect the flow in fundamentally different ways, implying that using Rayleigh drag to mimic magnetic drag is inaccurate. We exhaustively lay down the theoretical formalism (dispersion relations, governing equations and time-dependent wave solutions) for a broad suite of models. In all situations, we derive the steady state of an atmosphere, which is relevant to interpreting infrared phase and eclipse maps of exoplanetary atmospheres. We elucidate a pinching effect that confines the atmospheric structure to be near the equator. Our suite of analytical models may be used to decisively develop physical intuition and as a reference point for three-dimensional, magnetohydrodynamic (MHD) simulations of atmospheric circulation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Analytical Models of Exoplanetary Atmospheres. II. Radiative Transfer via the Two-stream Approximation

We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion (which describes radiative transfer in the deep interior) and solutions for the temperature-pressure profiles. Generally, the problem is mathematically under-determined unless a set of closures (Eddington coefficients) is specified. We demonstrate that the hemispheric (or hemi-isotropic) closure naturally derives from the radiative transfer equation if energy conservation is obeyed, while the Eddington closure produces spurious enhancements of both reflected light and thermal emission. We concoct recipes for implementing two-stream radiative transfer in stand-alone numerical calculations and general circulation models. We use our two-stream solutions to construct toy models of the runaway greenhouse effect. We present a new solution for temperature-pressure profiles with a non-constant optical opacity and elucidate the effects of non-isotropic scattering in the optical and infrared. We derive generalized expressions for the spherical and Bond albedos and the photon deposition depth. We demonstrate that the value of the optical depth corresponding to the photosphere is not always 2/3 (Milne's solution) and depends on a combination of stellar irradiation, internal heat and the properties of scattering both in optical and infrared. Finally, we derive generalized expressions for the total, net, outgoing and incoming fluxes in the convective regime.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Atmospheric Dynamics of Hot Exoplanets

The characterization of exoplanetary atmospheres has come of age in the last decade, as astronomical techniques now allow for albedos, chemical abundances, temperature profiles and maps, rotation periods and even wind speeds to be measured. Atmospheric dynamics sets the background state of density, temperature and velocity that determines or influences the spectral and temporal appearance of an exoplanetary atmosphere. Hot exoplanets are most amenable to these characterization techniques; in the present review, we focus on highly-irradiated, large exoplanets (the "hot Jupiters"), as astronomical data begin to confront theoretical questions. We summarize the basic atmospheric quantities inferred from the astronomical observations. We review the state of the art by addressing a series of current questions and look towards the future by considering a separate set of exploratory questions. Attaining the next level of understanding will require a concerted effort of constructing multi-faceted, multi-wavelength datasets for benchmark objects. Understanding clouds presents a formidable obstacle, as they introduce degeneracies into the interpretation of spectra, yet their properties and existence are directly influenced by atmospheric dynamics. Confronting general circulation models with these multi-faceted, multi-wavelength datasets will help us understand these and other degeneracies. The coming decade will witness a decisive confrontation of theory and simulation by the next generation of astronomical data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constraining the Atmospheric Composition of the Day-Night Terminators of HD 189733b : Atmospheric Retrieval with Aerosols

A number of observations have shown that Rayleigh scattering by aerosols dominates the transmission spectrum of HD 189733b at wavelengths shortward of 1 $μ$m. In this study, we retrieve a range of aerosol distributions consistent with transmission spectroscopy between 0.3-24 $μ$m that were recently re-analyzed by Pont et al. (2013). To constrain the particle size and the optical depth of the aerosol layer, we investigate the degeneracies between aerosol composition, temperature, planetary radius, and molecular abundances that prevent unique solutions for transit spectroscopy. Assuming that the aerosol is composed of MgSiO$_3$, we suggest that a vertically uniform aerosol layer over all pressures with a monodisperse particle size smaller than about 0.1 $μ$m and an optical depth in the range 0.002-0.02 at 1 $μ$m provides statistically meaningful solutions for the day/night terminator regions of HD 189733b. Generally, we find that a uniform aerosol layer provide adequate fits to the data if the optical depth is less than 0.1 and the particle size is smaller than 0.1 $μ$m, irrespective of the atmospheric temperature, planetary radius, aerosol composition, and gaseous molecules. Strong constraints on the aerosol properties are provided by spectra at wavelengths shortward of 1 $μ$m as well as longward of 8 $μ$m, if the aerosol material has absorption features in this region. We show that these are the optimal wavelengths for quantifying the effects of aerosols, which may guide the design of future space observations. The present investigation indicates that the current data offer sufficient information to constrain some of the aerosol properties of HD189733b, but the chemistry in the terminator regions remains uncertain.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High Resolution Transmission Spectroscopy as a Diagnostic for Jovian Exoplanet Atmospheres: Constraints from Theoretical Models

We present high resolution transmission spectra of giant planet atmospheres from a coupled 3-D atmospheric dynamics and transmission spectrum model that includes Doppler shifts which arise from winds and planetary motion. We model jovian planets covering more than two orders of magnitude in incident flux, corresponding to planets with 0.9 to 55 day orbital periods around solar-type stars. The results of our 3-D dynamical models reveal certain aspects of high resolution transmission spectra that are not present in simple 1-D models. We find that the hottest planets experience strong substellar to anti-stellar (SSAS) winds, resulting in transmission spectra with net blue shifts of up to 3 km s$^{-1}$, whereas less irradiated planets show almost no net Doppler shifts. Compared to 1-D models, peak line strengths are significantly reduced for the hottest atmospheres owing to Doppler broadening from a combination of rotation (which is faster for close-in planets under the assumption of tidal locking) and atmospheric winds. Finally, high resolution transmission spectra may be useful in studying the atmospheres of exoplanets with optically thick clouds since line cores for very strong transitions should remain optically thick to very high altitude. High resolution transmission spectra are an excellent observational test for the validity of 3-D atmospheric dynamics models, because they provide a direct probe of wind structures and heat circulation. Ground-based exoplanet spectroscopy is currently on the verge of being able to verify some of our modeling predictions, most notably the dependence of SSAS winds on insolation. We caution that interpretation of high resolution transmission spectra based on 1-D atmospheric models may be inadequate, as 3-D atmospheric motions can produce a noticeable effect on the absorption signatures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An Integral View of Balmer-dominated Shocks in Supernova Remnants

We present integral-field spectroscopic observations with the VIMOS-IFU at the VLT of fast (2000-3000 km/s) Balmer-dominated shocks surrounding the northwestern rim of the remnant of supernova 1006. The high spatial and spectral resolution of the instrument enable us to show that the physical characteristics of the shocks exhibit a strong spatial variation over few atomic scale lengths across 133 sky locations. Our results point to the presence of a population of non-thermal protons (10-100 keV) which might well be the seed particles for generating high-energy cosmic rays. We also present observations of Tycho's supernova remnant taken with the narrow-band tunable filter imager OSIRIS at the GTC and the Fabry-Perot interferometer GHaFaS at the WHT to resolve respectively the broad and narrow Hα lines across a large part of the remnant.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Atmospheric Retrieval Analysis of the Directly Imaged Exoplanet HR 8799b

Directly-imaged exoplanets are unexplored laboratories for the application of the spectral and temperature retrieval method, where the chemistry and composition of their atmospheres are inferred from inverse modeling of the available data. As a pilot study, we focus on the extrasolar gas giant HR 8799b for which more than 50 data points are available. We upgrade our non-linear optimal estimation retrieval method to include a phenomenological model of clouds that requires the cloud optical depth and monodisperse particle size to be specified. Previous studies have focused on forward models with assumed values of the exoplanetary properties; there is no consensus on the best-fit values of the radius, mass, surface gravity and effective temperature of HR 8799b. We show that cloudfree models produce reasonable fits to the data if the atmosphere is of super-solar metallicity and non-solar elemental abundances. Intermediately cloudy models with moderate values of the cloud optical depth and micron-sized particles provide an equally reasonable fit to the data and require a lower mean molecular weight. We report our best-fit values for the radius, mass, surface gravity and effective temperature of HR 8799b. The mean molecular weight is about 3.8, while the carbon-to-oxygen ratio is about unity due to the prevalence of carbon monoxide. Our study emphasizes the need for robust claims about the nature of an exoplanetary atmosphere to be based on analyses involving both photometry and spectroscopy and inferred from beyond a few photometric data points, such as are typically reported for hot Jupiters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Debris discs around M stars: non-existence versus non-detection

Motivated by the reported dearth of debris discs around M stars, we use survival models to study the occurrence of planetesimal discs around them. These survival models describe a planetesimal disc with a small number of parameters, determine if it may survive a series of dynamical processes and compute the associated infrared excess. For the WISE satellite, we demonstrate that the dearth of debris discs around M stars may be attributed to the small semi-major axes generally probed if either: 1. the dust grains behave like blackbodies emitting at a peak wavelength coincident with the observed one; 2. or the grains are hotter than predicted by their blackbody temperatures and emit at peak wavelengths that are shorter than the observed one. At these small distances from the M star, planetesimals are unlikely to survive or persist for time scales of 300 Myr or longer if the disc is too massive. Conversely, our survival models allow for the existence of a large population of low-mass debris discs that are too faint to be detected with current instruments. However, our interpretation becomes less clear and large infrared excesses are allowed if only one of these scenarios holds: 3. the dust grains are hotter than blackbody and predominantly emit at the observed wavelength; 4. or are blackbody in nature and emit at peak wavelengths longer than the observed one. Both scenarios imply that the parent planetesimals reside at larger distances from the star than inferred if the dust grains behaved like blackbodies. In all scenarios, we show that the infrared excesses detected at 22 and 70 microns from AU Mic are easily reconciled with its young age. We elucidate the conditions under which stellar wind drag may be neglected when considering dust populations around M stars. The WISE satellite should be capable of detecting debris discs around young M stars with ages on the order of 10 Myr.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Inference of Inhomogeneous Clouds in an Exoplanet Atmosphere

We present new visible and infrared observations of the hot Jupiter Kepler-7b to determine its atmospheric properties. Our analysis allows us to 1) refine Kepler-7b's relatively large geometric albedo of Ag=0.35+-0.02, 2) place upper limits on Kepler-7b thermal emission that remains undetected in both Spitzer bandpasses and 3) report a westward shift in the Kepler optical phase curve. We argue that Kepler-7b's visible flux cannot be due to thermal emission or Rayleigh scattering from H2 molecules. We therefore conclude that high altitude, optically reflective clouds located west from the substellar point are present in its atmosphere. We find that a silicate-based cloud composition is a possible candidate. Kepler-7b exhibits several properties that may make it particularly amenable to cloud formation in its upper atmosphere. These include a hot deep atmosphere that avoids a cloud cold trap, very low surface gravity to suppress cloud sedimentation, and a planetary equilibrium temperature in a range that allows for silicate clouds to potentially form in the visible atmosphere probed by Kepler. Our analysis does not only present evidence of optically thick clouds on Kepler-7b but also yields the first map of clouds in an exoplanet atmosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Deep Blue Color of HD189733b: Albedo Measurements with Hubble Space Telescope/Space Telescope Imaging Spectrograph at Visible Wavelengths

We present a secondary eclipse observation for the hot Jupiter HD189733b across the wavelength range 290-570nm made using the Space Telescope Imaging Spectrograph on the Hubble Space Telescope. We measure geometric albedos of Ag = 0.40 \pm 0.12 across 290-450nm and Ag < 0.12 across 450-570nm at 1-sigma confidence. The albedo decrease toward longer wavelengths is also apparent when using six wavelength bins over the same wavelength range. This can be interpreted as evidence for optically thick reflective clouds on the dayside hemisphere with sodium absorption suppressing the scattered light signal beyond ~450nm. Our best-fit albedo values imply that HD189733b would appear a deep blue color at visible wavelengths.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Understanding Trends Associated with Clouds in Irradiated Exoplanets

Unlike previously explored relationships between the properties of hot Jovian atmospheres, the geometric albedo and the incident stellar flux do not exhibit a clear correlation, as revealed by our re-analysis of Q0 to Q14 Kepler data. If the albedo is primarily associated with the presence of clouds in these irradiated atmospheres, a holistic modeling approach needs to relate the following properties: the strength of stellar irradiation (and hence the strength and depth of atmospheric circulation), the geometric albedo (which controls both the fraction of starlight absorbed and the pressure level at which it is predominantly absorbed) and the properties of the embedded cloud particles (which determine the albedo). The anticipated diversity in cloud properties renders any correlation between the geometric albedo and the stellar flux to be weak and characterized by considerable scatter. In the limit of vertically uniform populations of scatterers and absorbers, we use an analytical model and scaling relations to relate the temperature-pressure profile of an irradiated atmosphere and the photon deposition layer and to estimate if a cloud particle will be lofted by atmospheric circulation. We derive an analytical formula for computing the albedo spectrum in terms of the cloud properties, which we compare to the measured albedo spectrum of HD 189733b by Evans et al. (2013). Furthermore, we show that whether an optical phase curve is flat or sinusoidal depends on whether the particles are small or large as defined by the Knudsen number. This may be an explanation for why Kepler-7b exhibits evidence for the longitudinal variation in abundance of condensates, while Kepler-12b shows no evidence for the presence of condensates, despite the incident stellar flux being similar for both exoplanets.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Excitation and charge transfer in hydrogen-proton collisions at 5--80 keV and application to astrophysical shocks

In astrophysical regimes where the collisional excitation of hydrogen atoms is relevant, the cross sections for the interactions of hydrogen atoms with electrons and protons are necessary for calculating line profiles and intensities. In particular, at relative velocities exceeding ~1000 km/s, collisional excitation by protons dominates over that by electrons. Surprisingly, the hydrogen-proton cross sections at these velocities do not exist for atomic levels of n >= 4, forcing researchers to utilize extrapolation via inaccurate scaling laws. In this study, we present a faster and improved algorithm for computing cross sections for the hydrogen-proton collisional system, including excitation and charge transfer to the n >= 2 levels of the hydrogen atom. We develop a code named BDSCx which directly solves the Schrodinger equation with variable (but non-adaptive) resolution and utilizes a hybrid spatial-Fourier grid. Our novel hybrid grid reduces the number of grid points needed from ~4000 n^6 (for a "brute force", Cartesian grid) to ~2000 n^4 and speeds up the computation by a factor ~50 for calculations going up to n = 4 . We present (l,m)-resolved results for charge-transfer and excitation final states for n = 2--4 and for projectile energies of 5--80 keV, as well as fitting functions for the cross sections. The ability to accurately compute proton-hydrogen cross sections to n = 4 allows us to calculate the Balmer decrement, the ratio of Balmer alpha to Balmer beta line intensities. We find that the Balmer decrement starts to increase beyond its largely constant value of 2--3 below 10 keV, reaching values of 4--5 at 5 keV, thus complicating its use as a diagnostic of dust extinction when fast (~1000$ km/s) shocks are impinging upon the ambient interstellar medium.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the Existence of Shocks in Irradiated Exoplanetary Atmospheres

Supersonic flows are expected to exist in the atmospheres of irradiated exoplanets, but the question of whether shocks develop lingers. Specifically, it reduces to whether continuous flow in a closed loop may become supersonic and if some portions of the supersonic flow steepen into shocks. We first demonstrate that continuous, supersonic flow may exist in two flavors: isentropic and non-isentropic, with shocks being included in the latter class of solutions. Supersonic flow is a necessary but insufficient condition for shocks to develop. The development of a shock requires the characteristics of neighboring points in a flow to intersect. We demonstrate that the intersection of characteristics may be quantified via knowledge of the Mach number. Finally, we examine 3D simulations of hot Jovian atmospheres and demonstrate that shock formation is expected to occur mostly on the dayside hemisphere, upstream of the substellar point, because the enhanced temperatures near the substellar point provide a natural pressure barrier for the returning flow. Understanding the role of shocks in irradiated exoplanetary atmospheres is relevant to correctly modeling observables such as the peak offsets of infrared phase curves.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the Stability of Super-Earth Atmospheres

We investigate the stability of super Earth atmospheres around M stars using a 7-parameter, analytical framework. We construct stability diagrams in the parameter space of exoplanetary radius versus semi-major axis and elucidate the regions in which the atmospheres are stable against the condensation of their major constituents, out of the gas phase, on their permanent nightside hemispheres. We find that super Earth atmospheres which are nitrogen-dominated ("Earth-like") occupy a smaller region of allowed parameter space, compared to hydrogen-dominated atmospheres, because of the dual effects of diminished advection and enhanced radiative cooling. Furthermore, some super Earths which reside within the habitable zones of M stars may not possess stable atmospheres, depending on the mean molecular weight and infrared photospheric pressure of their atmospheres. We apply our stability diagrams to GJ 436b and GJ 1214b, and demonstrate that atmospheric compositions with high mean molecular weights are disfavoured if these exoplanets possess solid surfaces and shallow atmospheres. Finally, we construct stability diagrams tailored to the Kepler dataset, for G and K stars, and predict that about half of the exoplanet candidates are expected to habour stable atmospheres if Earth-like conditions are assumed. We include 55 Cancri e and CoRoT-7b in our stability diagram for G stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Effects of Irradiation on Hot Jovian Atmospheres: Heat Redistribution and Energy Dissipation

Hot Jupiters, due to the proximity to their parent stars, are subjected to a strong irradiating flux which governs their radiative and dynamical properties. We compute a suite of 3D circulation models with dual-band radiative transfer, exploring a relevant range of irradiation temperatures, both with and without temperature inversions. We find that, for irradiation temperatures T \lesssim 2000 K, heat redistribution is very efficient, producing comparable day- and night-side fluxes. For Tirr \approx 2200-2400 K, the redistribution starts to break down, resulting in a high day-night flux contrast. Our simulations indicate that the efficiency of redistribution is primarily governed by the ratio of advective to radiative timescales. Models with temperature inversions display a higher day-night contrast due to the deposition of starlight at higher altitudes, but we find this opacity-driven effect to be secondary compared to the effects of irradiation. The hotspot offset from the substellar point is large when insolation is weak and redistribution is efficient, and decreases as redistribution breaks down. The atmospheric flow can be potentially subjected to the Kelvin-Helmholtz instability (as indicated by the Richardson number) only in the uppermost layers, with a depth that penetrates down to pressures of a few millibars at most. Shocks penetrate deeper, down to several bars in the hottest model. Ohmic dissipation generally occurs down to deeper levels than shock dissipation (to tens of bars), but the penetration depth varies with the atmospheric opacity. The total dissipated Ohmic power increases steeply with the strength of the irradiating flux and the dissipation depth recedes into the atmosphere, favoring radius inflation in the most irradiated objects.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Influence of Atmospheric Scattering and Absorption on Ohmic Dissipation in Hot Jupiters

Using semi-analytical, one-dimensional models, we elucidate the influence of scattering and absorption on the degree of Ohmic dissipation in hot Jovian atmospheres. With the assumption of Saha equilibrium, the variation in temperature is the main driver of the variations in the electrical conductivity, induced current and Ohmic power dissipated. Atmospheres possessing temperature inversions tend to dissipate most of the Ohmic power superficially, at high altitudes, whereas those without temperature inversions are capable of greater dissipation deeper down. Scattering in the optical range of wavelengths tends to cool the lower atmosphere, thus reducing the degree of dissipation at depth. Purely absorbing cloud decks (in the infrared), of a finite extent in height, allow for localized reductions in dissipation and may reverse a temperature inversion if they are dense and thick enough, thus greatly enhancing the dissipation at depth. If Ohmic dissipation is the mechanism for inflating hot Jupiters, then variations in the atmospheric opacity (which may be interpreted as arising from variations in metallicity and cloud/haze properties) and magnetic field strength naturally produce a scatter in the measured radii at a given strength of irradiation. Future work will determine if these effects are dominant over evolutionary effects, which also contribute a scatter to the measured radii.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Atmospheric circulation of tidally locked exoplanets: a suite of benchmark tests for dynamical solvers

The complexity of atmospheric modelling and its inherent non-linearity, together with the limited amount of data of exoplanets available, motivate model intercomparisons and benchmark tests. In the geophysical community, the Held-Suarez test is a standard benchmark for comparing dynamical core simulations of the Earth's atmosphere with different solvers, based on statistically-averaged flow quantities. In the present study, we perform analogues of the Held-Suarez test for tidally-locked exoplanets with the GFDL-Princeton Flexible Modeling System (FMS) by subjecting both the spectral and finite difference dynamical cores to a suite of tests, including the standard benchmark for Earth, a hypothetical tidally-locked Earth, a "shallow" hot Jupiter model and a "deep" model of HD 209458b. We find qualitative and quantitative agreement between the solvers for the Earth, tidally-locked Earth and shallow hot Jupiter benchmarks, but the agreement is less than satisfactory for the deep model of HD 209458b. Further investigation reveals that closer agreement may be attained by arbitrarily adjusting the values of the horizontal dissipation parameters in the two solvers, but it remains the case that the magnitude of the horizontal dissipation is not easily specified from first principles. Irrespective of radiative transfer or chemical composition considerations, our study points to limitations in our ability to accurately model hot Jupiter atmospheres with meteorological solvers at the level of ten percent for the temperature field and several tens of percent for the velocity field. Direct wind measurements should thus be particularly constraining for the models. Our suite of benchmark tests also provides a reference point for researchers wishing to adapt their codes to study the atmospheric circulation regimes of tidally-locked Earths/Neptunes/Jupiters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Atmospheric circulation of tidally locked exoplanets: II. Dual-band radiative transfer and convective adjustment

Improving upon our purely dynamical work, we present three-dimensional simulations of the atmospheric circulation on Earth-like (exo)planets and hot Jupiters using the GFDL-Princeton Flexible Modeling System (FMS). As the first steps away from the dynamical benchmarks of Heng, Menou & Phillipps (2011), we add dual-band radiative transfer and dry convective adjustment schemes to our computational setup. Our treatment of radiative transfer assumes stellar irradiation to peak at a wavelength shorter than and distinct from that at which the exoplanet re-emits radiation ("shortwave" versus "longwave"), and also uses a two-stream approximation. Convection is mimicked by adjusting unstable lapse rates to the dry adiabat. The bottom of the atmosphere is bounded by a uniform slab with a finite thermal inertia. For our models of hot Jupiters, we include an analytical formalism for calculating temperature-pressure profiles, in radiative equilibrium, which accounts for the effect of collision-induced absorption via a single parameter. We discuss our results within the context of: the predicted temperature-pressure profiles and the absence/presence of a temperature inversion; the possible maintenance, via atmospheric circulation, of the putative high-altitude, shortwave absorber expected to produce these inversions; the angular/temporal offset of the hot spot from the substellar point, its robustness to our ignorance of hyperviscosity and hence its utility in distinguishing between different hot Jovian atmospheres; and various zonal-mean flow quantities. Our work bridges the gap between three-dimensional simulations which are purely dynamical and those which incorporate multi-band radiative transfer, thus contributing to the construction of a required hierarchy of three-dimensional theoretical models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Estimating the mass of the debris disc in HD 69830

We present a method to estimate the mass of the debris disc in the HD 69830 system, which also hosts three exoplanets with Neptune-like minimum masses. By considering the range of published stellar ages, we interpret the infrared emission from the debris disc as originating from a steady state, collisional cascade of dust grains. Using dynamical survival models subjected to observational constraints, we estimate the allowed range of disc masses. If the disc has an age of about 1 Gyr, then its mass is several times that of our asteroid belt. The maximum allowed age for the disc and the number of planetesimals it contains are determined by the assumed value for the binding energy of the planetesimals. If one insists on interpreting the disc as being transient, then this mass estimate becomes an upper limit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gliese 581g as a scaled-up version of Earth: atmospheric circulation simulations

We use three-dimensional simulations to study the atmospheric circulation on the first Earth-sized exoplanet discovered in the habitable zone of an M star. We treat Gliese 581g as a scaled-up version of Earth by considering increased values for the exoplanetary radius and surface gravity, while retaining terrestrial values for parameters which are unconstrained by current observations. We examine the long-term, global temperature and wind maps near the surface of the exoplanet --- the climate. The specific locations for habitability on Gliese 581g depend on whether the exoplanet is tidally-locked and how fast radiative cooling occurs on a global scale. Independent of whether the existence of Gliese 581g is confirmed, our study highlights the use of general circulation models to quantify the atmospheric circulation on potentially habitable, Earth-sized exoplanets, which will be the prime targets of exoplanet discovery and characterization campaigns in the next decade.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

HST-COS Observations of Hydrogen, Helium, Carbon and Nitrogen Emission from the SN 1987A Reverse Shock

We present the most sensitive ultraviolet observations of Supernova 1987A to date. Imaging spectroscopy from the Hubble Space Telescope-Cosmic Origins Spectrograph shows many narrow (dv \sim 300 km/s) emission lines from the circumstellar ring, broad (dv \sim 10 -- 20 x 10^3 km/s) emission lines from the reverse shock, and ultraviolet continuum emission. The high signal-to-noise (> 40 per resolution element) broad LyA emission is excited by soft X-ray and EUV heating of mostly neutral gas in the circumstellar ring and outer supernova debris. The ultraviolet continuum at λ> 1350A can be explained by HI 2-photon emission from the same region. We confirm our earlier, tentative detection of NV λ1240 emission from the reverse shock and we present the first detections of broad HeII \lambda1640, CIV \lambda1550, and NIV] \lambda1486 emission lines from the reverse shock. The helium abundance in the high-velocity material is He/H = 0.14 +/- 0.06. The NV/H-alpha line ratio requires partial ion-electron equilibration (T_{e}/T_{p} \approx 0.14 - 0.35). We find that the N/C abundance ratio in the gas crossing the reverse shock is significantly higher than that in the circumstellar ring, a result that may be attributed to chemical stratification in the outer envelope of the supernova progenitor. The N/C abundance ratio may have been stratified prior to the ring expulsion, or this result may indicate continued CNO processing in the progenitor subsequent to the expulsion of the circumstellar ring.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the effects of clouds and hazes in the atmospheres of hot Jupiters: semi-analytical temperature-pressure profiles

Motivated by the work of Guillot (2010), we present a semi-analytical formalism for calculating the temperature-pressure profiles in hot Jovian atmospheres which includes the effects of clouds/hazes and collision-induced absorption. Using the dual-band approximation, we assume that stellar irradiation and thermal emission from the hot Jupiter occur at distinct wavelengths ("shortwave" versus "longwave"). For a purely absorbing cloud/haze, we demonstrate its dual effect of cooling and warming the upper and lower atmosphere, respectively, which modifies, in a non-trivial manner, the condition for whether a temperature inversion is present in the upper atmosphere. The warming effect becomes more pronounced as the cloud/haze deck resides at greater depths. If it sits below the shortwave photosphere, the warming effect becomes either more subdued or ceases altogether. If shortwave scattering is present, its dual effect is to warm and cool the upper and lower atmosphere, respectively, thus counteracting the effects of enhanced longwave absorption by the cloud/haze. We make a tentative comparison of a 4-parameter model to the temperature-pressure data points inferred from the observations of HD 189733b and estimate that its Bond albedo is approximately 10%. Besides their utility in developing physical intuition, our semi-analytical models are a guide for the parameter space exploration of hot Jovian atmospheres via three-dimensional simulations of atmospheric circulation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Dependence of Brown Dwarf Radii on Atmospheric Metallicity and Clouds: Theory and Comparison with Observations

Employing realistic and consistent atmosphere boundary conditions, we have generated evolutionary models for brown dwarfs and very-low-mass stars (VLMs) for different metallicities ([Fe/H]), with and without clouds. We find that the spread in radius at a given mass and age can be as large as $\sim$10% to $\sim$25%, with higher-metallicity, higher-cloud-thickness atmospheres resulting quite naturally in larger radii. For each 0.1 dex increase in [Fe/H], radii increase by $\sim$1% to $\sim$2.5%, depending upon age and mass. We also find that, while for smaller masses and older ages brown dwarf radii decrease with increasing helium fraction ($Y$) (as expected), for more massive brown dwarfs and a wide range of ages they increase with helium fraction. The increase in radius in going from $Y=0.25$ to $Y=0.28$ can be as large as $\sim$0.025 \rj\ ($\sim$2.5%). Furthermore, we find that for VLMs an increase in atmospheric metallicity from 0.0 to 0.5 dex, increases radii by $\sim$4%, and from -0.5 to 0.5 dex by $\sim$10%. Therefore, we suggest that opacity due to higher metallicity might naturally account for the apparent radius anomalies in some eclipsing VLM systems. Ten to twenty-five percent variations in radius exceed errors stemming from uncertainities in the equation of state alone. This serves to emphasize that transit and eclipse measurements of brown dwarf radii constrain numerous effects collectively, importantly including the atmosphere and condensate cloud models, and not just the equation of state. At all times, one is testing a multi-parameter theory, and not a universal radius$-$mass relation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Long-lived planetesimal discs

We investigate the survival of planetesimal discs over Gyr timescales, using a unified approach that is applicable to all Keplerian discs of solid bodies -- dust grains, asteroids, planets, etc. Planetesimal discs can be characterized locally by four parameters: surface density, semi-major axis, planetesimal size and planetesimal radial velocity dispersion. Any planetesimal disc must have survived all dynamical processes, including gravitational instability, dynamical chaos, gravitational scattering, physical collisions, and radiation forces, that would lead to significant evolution over its lifetime. These processes lead to a rich set of constraints that strongly restrict the possible properties of long-lived discs. Within this framework, we also discuss the detection of planetesimal discs using radial velocity measurements, transits, microlensing, and the infrared emission from the planetesimals themselves or from dust generated by planetesimal collisions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Observing Supernova 1987A with the Refurbished Hubble Space Telescope

Observations with the Hubble Space Telescope (HST), conducted since 1990, now offer an unprecedented glimpse into fast astrophysical shocks in the young remnant of supernova 1987A. Comparing observations taken in 2010 using the refurbished instruments on HST with data taken in 2004, just before the Space Telescope Imaging Spectrograph failed, we find that the Ly-a and H-a lines from shock emission continue to brighten, while their maximum velocities continue to decrease. We observe broad blueshifted Ly-a, which we attribute to resonant scattering of photons emitted from hotspots on the equatorial ring. We also detect NV~λλ1239,1243 A line emission, but only to the red of Ly-A. The profiles of the NV lines differ markedly from that of H-a, suggesting that the N^{4+} ions are scattered and accelerated by turbulent electromagnetic fields that isotropize the ions in the collisionless shock.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Vortices as nurseries for planetesimal formation in protoplanetary discs

Turbulent, two-dimensional, hydrodynamic flows are characterized by the emergence of coherent, long-lived vortices without a need to invoke special initial conditions. Vortices have the ability to sequester particles, with typical radii from about 1 mm to 10 cm, that are slightly decoupled from the gas. A generic feature of discs with surface density and effective temperature profiles that are decreasing, power-law functions of radial distance is that four vortex zones exist for a fixed particle size. In particular, two of the zones form an annulus at intermediate radial distances within which small particles reside. Particle capture by vortices occurs on a dynamical time scale near and at the boundaries of this annulus. As the disc ages and the particles grow via coagulation, the size of the annulus shrinks. Older discs prefer to capture smaller particles because the gas surface density decreases with time, a phenomenon we term "vortex aging". More viscous, more dust-opaque and/or less massive discs can have vortices that age faster and trap a broader range of particle sizes throughout the lifetime of the disc. Thus, how efficiently a disc retains its mass in solids depends on the relative time scales between coagulation and vortex aging. If vortices form in protoplanetary discs, they are important in discs with typical masses and for particles that are likely to condense out of the protostellar nebula. Particle capture also occurs at distances relevant to planet formation. Future infrared, submillimetre and centimetre observations of grain opacity as a function of radial distance will test the hypothesis that vortices serve as nurseries for particle growth in protoplanetary discs.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Balmer-Dominated Shocks: A Concise Review

A concise and critical review of Balmer-dominated shocks (BDSs) is presented, summarizing the state of theory and observations, including models with/without shock precursors and their synergy with atomic physics. Observations of BDSs in supernova remnants are reviewed on an object-by-object basis. The relevance of BDSs towards understanding the acceleration of cosmic rays in shocks is emphasized. Probable and possible detections of BDSs in astrophysical objects other than supernova remnants, including pulsar wind nebulae and high-redshift galaxies, are described. The case for the continued future of studying BDSs in astrophysics is made, including their relevance towards understanding electron-ion temperature equilibration in collisionless shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Magnetohydrodynamic Shallow Water Waves: Linear Analysis

We present a linear analysis of inviscid, incompressible, magnetohydrodynamic (MHD) shallow water systems. In spherical geometry, a generic property of such systems is the existence of five wave modes. Three of them (two magneto-Poincare modes and one magneto-Rossby mode) are previously known. The other two wave modes are strongly influenced by the magnetic field and rotation, and have substantially lower angular frequencies; as such, we term them "magnetostrophic modes". We obtain analytical functions for the velocity, height and magnetic field perturbations in the limit that the magnitude of the MHD analogue of Lamb's parameter is large. On a sphere, the magnetostrophic modes reside near the poles, while the other modes are equatorially confined. Magnetostrophic modes may be an ingredient in explaining the frequency drifts observed in Type I X-ray bursts from neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Planetesimal Disk Microlensing

Motivated by debris disk studies, we investigate the gravitational microlensing of background starlight by a planetesimal disk around a foreground star. We use dynamical survival models to construct a plausible example of a planetesimal disk and study its microlensing properties using established ideas of microlensing by small bodies. When a solar-type source star passes behind a planetesimal disk, the microlensing light curve may exhibit short-term, low-amplitude residuals caused by planetesimals several orders of magnitude below Earth mass. The minimum planetesimal mass probed depends on the photometric sensitivity and the size of the source star, and is lower when the planetesimal lens is located closer to us. Planetesimal lenses may be found more nearby than stellar lenses because the steepness of the planetesimal mass distribution changes how the microlensing signal depends on the lens/source distance ratio. Microlensing searches for planetesimals require essentially continuous monitoring programs that are already feasible and can potentially set constraints on models of debris disks, the supposed extrasolar analogues of Kuiper belts.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Dust Echoes from the Ambient Medium of Gamma-Ray Bursts

Long gamma-ray bursts (GRBs) are likely associated with the collapse of massive stars, which produce dust and are born in dusty environments. Absorption and scattering of ultraviolet/X-ray photons from the prompt, optical flash and afterglow emission of the GRB produce dust echoes. We perform time-dependent calculations of these echoes, accounting for the evolution of the dust grain distribution due to selective grain destruction by the GRB radiation, and for off-axis beaming. We explore cloud configurations of differing density and size -- the echo light curve and spectrum depend on the cloud radius, with larger clouds peaking at longer wavelengths. For a region ~ 3 pc in size with hydrogen density ~ 1000 per cubic centimeter, the echo spectrum peaks at ~ 3.6 microns and ~ 8.8 eV for thermal and scattered components, respectively. Dust echoes should be detectable with the Very Large Telescope up to z ~ 0.1, IRAC onboard the Spitzer Space Telescope up to z ~ 0.2, and NICMOS onboard the Hubble Space Telescope up to z ~ 0.3. Furthermore, the shape of the echo light curve allows one to infer: the jet opening angle; the inclination of the jet axis with respect to the line of sight; the size of the dust-emitting region. For sources with symmetric, bipolar jets, dust echoes exhibit two bumps in the light curve, making them easily distinguishable from the rebrightening due to an underlying supernova.

0 citations0 reviewsNo code linkedOpen access