Search

Earth-like planets orbiting M-dwarf stars, M-Earths, are currently the best targets to search for signatures of life. Life as we know it requires water. The habitability of M-Earths is jeopardized by water loss to space: high flux from young M-dwarf stars can drive the loss of 3–20 Earth oceans from otherwise habitable planets. We develop a 0-D box model for Earth-mass terrestrial exoplanets, orbiting within the habitable zone, which tracks water loss to space and exchange between reservoirs during an early surface magma ocean phase and the longer deep-water cycling phase. A key feature is the duration of the surface magma ocean, assumed concurrent with the runaway greenhouse. This time-scale can discriminate between desiccated planets, planets with desiccated mantles but substantial surface water, and planets with significant water sequestered in the mantle. A longer-lived surface magma ocean helps M-Earths retain water: dissolution of water in the magma provides a barrier against significant loss to space during the earliest, most active stage of the host M-dwarf, depending on the water saturation limit of the magma. Although a short-lived basal magma ocean can be beneficial to surface habitability, a long-lived basal magma ocean may sequester significant water in the mantle at the detriment of surface habitability. We find that magma oceans and deep-water cycling can maintain or recover habitable surface conditions on Earth-like planets at the inner edge of the habitable zone around late M-dwarf stars – these planets would otherwise be desiccated if they form with less than ∼10 terrestrial oceans of water. [ABSTRACT FROM AUTHOR]

The James Webb Space Telescope (JWST) will provide an opportunity to investigate the atmospheres of potentially habitable planets. Aerosols significantly mute molecular features in transit spectra because they prevent light from probing the deeper layers of the atmosphere. Earth occasionally has stratospheric/high tropospheric clouds at 15–20 km that could substantially limit the observable depth of the underlying atmosphere. We use solar occultations of Earth's atmosphere to create synthetic JWST transit spectra of Earth analogues orbiting dwarf stars. Unlike previous investigations, we consider both clear and cloudy sightlines from the SCISAT satellite. We find that the maximum difference in effective thickness of the atmosphere between a clear and globally cloudy atmosphere is 8.5 km at 2.28   |$\mu$| m, with a resolution of 0.02  |$\mu$| m. After incorporating the effects of refraction and Pandexo's noise modelling, we find that JWST would not be able to detect Earth-like stratospheric clouds if an exo-Earth was present in the TRAPPIST-1 system, as the cloud spectrum differs from the clear spectrum by a maximum of 10 ppm. These stratospheric clouds are also not robustly detected by TauREx when performing spectral retrieval for a cloudy TRAPPIST-1 planet. However, if an Earth-sized planet were to orbit in a white dwarf's habitable zone, then we predict that JWST 's NIRSpec would be able to detect its stratospheric clouds after only four transits. We conclude that stratospheric clouds would not impede JWST transit spectroscopy or the detection of biosignatures for Earth-like atmospheres. [ABSTRACT FROM AUTHOR]

Astronomers have discovered a handful of exoplanets with rocky bulk compositions but orbiting that orbit so close to their host star that the surface of the planet must be at least partially molten. It is expected that the dayside of such “lava planets” harbors a rock-vapor atmosphere that flows quickly toward the airless nightsideâ€"this partial atmosphere is critical to the interpretation of lava planet observations, but transports negligible heat toward the nightside. As a result, the surface temperature of the magma ocean may range from 3000 K near the substellar point down to 1500 K near the dayâ€"night terminator. We use simple models incorporating the thermodynamics and geochemistry of partial melt to predict the physical and chemical properties of the magma ocean as a function of the distance from the substellar point. Our principal findings are that: (1) the dayside magma ocean is much deeper than previously thought, probably extending down to the coreâ€"mantle boundary below the substellar point of an Earth-sized planet; (2) much of the dayside is only partially molten, leading to gradients in the surface chemistry of the magma ocean; and (3) the temperature at the base of the silicate mantle is as important as the surface temperature. In the most extreme cases, lava planet interiors could be cold enough such that thermal stratification below the substellar point is gravitationally stable. These findings have important implications for the dynamics of the magma ocean, as well as the composition and dynamics of the atmosphere. [ABSTRACT FROM AUTHOR]

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

The field of exoplanet atmospheric characterization is trending towards comparative studies involving many planetary systems, and using Bayesian hierarchical modelling is a natural next step. Here we demonstrate two use cases. We first use hierarchical modelling to quantify variability in repeated observations by reanalysing a suite of 10 Spitzer secondary eclipse observations of the hot Jupiter XO-3 b. We compare three models: one where we fit 10 separate eclipse depths, one where we use a single eclipse depth for all 10 observations, and a hierarchical model. By comparing the widely applicable information criterion of each model, we show that the hierarchical model is preferred over the others. The hierarchical model yields less scatter across the suite of eclipse depths – and higher precision on the individual eclipse depths – than does fitting the observations separately. We find that the hierarchical eclipse depth uncertainty is larger than the uncertainties on the individual eclipse depths, which suggests either slight astrophysical variability or that single eclipse observations underestimate the true eclipse depth uncertainty. Finally, we fit a suite of published dayside brightness measurements for 37 planets using a hierarchical model of brightness temperature versus irradiation temperature. The hierarchical model gives tighter constraints on the individual brightness temperatures than the non-hierarchical model. Although we tested hierarchical modelling on Spitzer eclipse data of hot Jupiters, it is applicable to observations of smaller planets like hot Neptunes and super-Earths, as well as for photometric and spectroscopic transit or phase-curve observations. [ABSTRACT FROM AUTHOR]

Transit searches have uncovered Earth-size planets orbiting so close to their host star that their surface should be molten, so-called lava planets. We present idealized simulations of the atmosphere of lava planet K2-141b and calculate the return flow of material via circulation in the magma ocean. We then compare how pure Na, SiO, or SiO2 atmospheres would impact future observations. The more volatile Na atmosphere is thickest followed by SiO and SiO2, as expected. Despite its low vapour pressure, we find that a SiO2 atmosphere is easier to observe via transit spectroscopy due to its greater scale height near the day–night terminator and the planetary radial velocity and acceleration are very high, facilitating high dispersion spectroscopy. The special geometry that arises from very small orbits allows for a wide range of limb observations for K2-141b. After determining the magma ocean depth, we infer that the ocean circulation required for SiO steady-state flow is only 10−4 m s−1, while the equivalent return flow for Na is several orders of magnitude greater. This suggests that a steady-state Na atmosphere cannot be sustained and that the surface will evolve over time. [ABSTRACT FROM AUTHOR]

Water cycling between Earth's mantle and surface has previously been modelled and extrapolated to rocky exoplanets, but these studies neglected the host star. M-dwarf stars are more common than Sun-like stars and at least as likely to host temperate rocky planets (M-Earths). However, M dwarfs are active throughout their lifetimes; specifically, X-ray and extreme ultraviolet (XUV) radiation during their early evolution can cause rapid atmospheric loss on orbiting planets. The increased bolometric flux reaching M-Earths leads to warmer, moister upper atmospheres, while XUV radiation can photodissociate water molecules and drive hydrogen and oxygen escape to space. Here, we present a coupled model of deep-water cycling and water loss to space on M-Earths to explore whether these planets can remain habitable despite their volatile evolution. We use a cycling parametrization accounting for the dependence of mantle degassing on seafloor pressure, the dependence of regassing on mantle temperature, and the effect of water on mantle viscosity and thermal evolution. We assume the M dwarf's XUV radiation decreases exponentially with time, and energy-limited water loss with 30 per cent efficiency. We explore the effects of cycling and loss to space on planetary water inventories and water partitioning. Planet surfaces desiccated by loss can be rehydrated, provided there is sufficient water sequestered in the mantle to degas once loss rates diminish at later times. For a given water loss rate, the key parameter is the mantle overturn time-scale at early times: if the mantle overturn time-scale is longer than the loss time-scale, then the planet is likely to keep some of its water. [ABSTRACT FROM AUTHOR]

In this work, we outline a new method for estimating dayside effective temperatures of exoplanets and associated uncertainties using Gaussian process (GP) regression. By applying our method to simulated observations, we show that the GP method estimates uncertainty more robustly than other model-independent approaches. We find that unbiased estimates of effective temperatures can be made using as few as three broad-band measurements (white-light HST WFC3 and the two warm Spitzer IRAC channels), although we caution that estimates made using only IRAC can be significantly biased. We then apply our GP method to the 12 hot Jupiters in the literature whose secondary eclipse depths have been measured by WFC3 and IRAC channels 1 and 2: CoRoT-2 b; HAT-P-7 b; HD 189733 b; HD 209458 b; Kepler-13A b; TrES-3 b; WASP-4 b; WASP-12 b; WASP-18 b; WASP-33 b; WASP-43 b; and WASP-103 b. We present model-independent dayside effective temperatures for these planets, with uncertainty estimates that range from ±66 to ±136 K. [ABSTRACT FROM AUTHOR]

The Atmospheric Chemistry Experiment Fourier Transform Spectrometer on the SCISAT satellite has been measuring infrared (IR) transmission spectra of Earth during Solar occultations since 2004. We use these data to build an IR transit spectrum of Earth. Regions of low atmospheric opacity, known as windows, are of particular interest, as they permit observations of the planet's lower atmosphere. Even in the absence of clouds or refraction, imperfect transmittance leads to a minimum effective thickness of h min ≈ 4 km in the 10–12  |$\mu \mathrm{m}$| opacity window at a spectral resolution of R  = 103. None the less, at R  = 105, the maximum transmittance at the surface is around |${70}{{{\ \rm per\ cent}}}$|⁠. In principle, one can probe the troposphere of an Earth-like planet via high-dispersion transit spectroscopy in the mid-IR; in practice aerosols and/or refraction likely make this impossible. We simulate the transit spectrum of an Earth-like planet in the TRAPPIST-1 system. We find that a long-term near-IR (NIR) campaign with the James Webb Space Telescope (JWST) could readily detect CO2 and H2O, establishing the presence of an atmosphere. A mid-IR campaign or longer NIR campaign would be more challenging, but in principle could detect the biosignatures O3 and CH4. [ABSTRACT FROM AUTHOR]

The disc-integrated reflected brightness of an exoplanet changes as a function of time due to orbital and rotational motions coupled with an inhomogeneous albedo map. We have previously derived analytic reflected light curves for spherical harmonic albedo maps in the special case of a synchronously rotating planet on an edge-on orbit (Cowan, Fuentes & Haggard). In this paper, we present analytic reflected light curves for the general case of a planet on an inclined orbit, with arbitrary spin period and non-zero obliquity. We do so for two different albedo basis maps: bright points (δ-maps), and spherical harmonics ( $Y_\mathrm{ l}^m$ -maps). In particular, we use Wigner D-matrices to express an harmonic light curve for an arbitrary viewing geometry as a non-linear combination of harmonic light curves for the simpler edge-on, synchronously rotating geometry. These solutions will enable future exploration of the degeneracies and information content of reflected light curves, as well as fast calculation of light curves for mapping exoplanets based on time-resolved photometry. To these ends, we make available Exoplanet Analytic Reflected Lightcurves, a simple open-source code that allows rapid computation of reflected light curves. [ABSTRACT FROM AUTHOR]

In addition to the transits of a planet in front of its star and the eclipses of the planet by its star, researchers have reported flux variations at the orbital frequency and its harmonics: planetary reflection and/or emission and Doppler beaming of starlight produce one peak per orbit, while ellipsoidal variations of a tidally distorted star and/or planet produce two maxima per orbit. Researchers have also reported significant photometric variability at three times the orbital frequency, as yet unexplained. Reflected phase variations of homogeneous planets only contain power at the orbital frequency and its even harmonics. We show that odd harmonics can, however, be produced by an edge-on planet with a time-variable map, or an inclined planet with a north-south (N-S) asymmetric map. For tidally locked planets with thick atmospheres, either of these scenarios entail weather: planets with zero obliquity experience N-S symmetric stellar forcing. North-south asymmetry would therefore suggest stochastic localized features, i.e. weather. However, we find that previous claims of large-amplitude odd modes in Kepler photometry are artefacts of removing planetary transits rather than modelling them. The only reliable claims of odd harmonics remain HAT-P-7b and Kepler-13Ab, for which the third mode amplitude is 6-8 per cent of the planetary flux. Although time-variable albedo maps could in principle explain these odd harmonics, upper-limits on the infrared variability of other hot Jupiters make this scenario unlikely. We recommend further studying the tidal effects of close-in planets on their host stars, as this remains the only plausible hypothesis. [ABSTRACT FROM AUTHOR]

Kaib and Cowan (Kaib, N.A., Cowan, N.B. [2015]. Icarus 252, 161–174) recently used terrestrial planet formation simulations to conclude that the moon-forming impactor (Theia) had only a ∼5% or less chance of having the same oxygen isotope composition as Earth, while Mastrobuono-Battisti et al. (Mastrobuono-Battisti, A., Perets, H.B., Raymond, S.N. [2015]. Nature 520, 212–215) used seemingly similar simulations and methods to arrive at a higher value of ∼20% or more. Here we derive the results of both papers from a single set of simulations. Compared to Kaib and Cowan (Kaib, N.A., Cowan, N.B. [2015]. Icarus 252, 161–174), the analysis of Mastrobuono-Battisti et al. (Mastrobuono-Battisti, A., Perets, H.B., Raymond, S.N. [2015]. Nature 520, 212–215) systematically yields more massive Theia analogs and imposes flatter Δ 17 O gradients across the original protoplanetary disk. Both of these effects diminish isotopic differences between Earth and Theia analogs. While it is notoriously difficult to produce systems resembling our actual terrestrial planets, the analysis of Kaib and Cowan (Kaib, N.A., Cowan, N.B. [2015]. Icarus 252, 161–174) more often selects and analyzes Earth and Mars analogs at orbital locations near the real planets. Given this, we conclude that the greater isotopic differences between Earth and Theia found in Kaib and Cowan (Kaib, N.A., Cowan, N.B. [2015]. Icarus 252, 161–174) better reflect the predictions of terrestrial planet formation models. Finally, although simulation uncertainties and a terrestrial contribution to Moon formation enhance the fraction of Theia analogs consistent with the canonical giant impact hypothesis, this fraction still remains in the 5–8% range. [ABSTRACT FROM AUTHOR]

The final stage of terrestrial planet formation consists of several hundred approximately lunar mass bodies accreting into a few terrestrial planets. This final stage is stochastic, making it hard to predict which parts of the original planetesimal disk contributed to each of our terrestrial planets. Here we present an extensive suite of terrestrial planet formation simulations that allows quantitative analysis of this process. Although there is a general correlation between a planet’s location and the initial semi-major axes of its constituent planetesimals, we concur with previous studies that Venus, Earth, and Mars analogs have overlapping, stochastic feeding zones. We quantify the feeding zone width, Δ a , as the mass-weighted standard deviation of the initial semi-major axes of the planetary embryos and planetesimals that make up the final planet. The size of a planet’s feeding zone in our simulations does not correlate with its final mass or semi-major axis, suggesting there is no systematic trend between a planet’s mass and its volatile inventory. Instead, we find that the feeding zone of any planet more massive than 0.1 M ⊕ is roughly proportional to the radial extent of the initial disk from which it formed: Δ a ≈ 0.25 ( a max - a min ) , where a min and a max are the inner and outer edge of the initial planetesimal disk. These wide stochastic feeding zones have significant consequences for the origin of the Moon, since the canonical scenario predicts the Moon should be primarily composed of material from Earth’s last major impactor (Theia), yet its isotopic composition is indistinguishable from Earth. In particular, we find that the feeding zones of Theia analogs are significantly more stochastic than the planetary analogs. Depending on our assumed initial distribution of oxygen isotopes within the planetesimal disk, we find a ∼ 5% or less probability that the Earth and Theia will form with an isotopic difference equal to or smaller than the Earth and Moon’s. In fact we predict that every planetary mass body should be expected to have a unique isotopic signature. In addition, we find paucities of massive Theia analogs and high velocity Moon-forming collisions, two recently proposed explanations for the Moon’s isotopic composition. Our work suggests that there is still no scenario for the Moon’s origin that explains its isotopic composition with a high probability event. [ABSTRACT FROM AUTHOR]

Large terrestrial planets are expected to have muted topography and deep oceans, implying that most super-Earths should be entirely covered in water, so-called waterworlds. This is important because waterworlds lack a silicate weathering thermostat so their climate is predicted to be less stable than that of planets with exposed continents. In other words, the continuously habitable zone for waterworlds is much narrower than for Earth-like planets. A planet's water is partitioned, however, between a surface reservoir, the ocean, and an interior reservoir, the mantle. Plate tectonics transports water between these reservoirs on geological timescales. Degassing of melt at mid-ocean ridges and serpentinization of oceanic crust depend negatively and positively on seafloor pressure, respectively, providing a stabilizing feedback on long-term ocean volume. Motivated by Earth's approximately steady-state deep water cycle, we develop a two-box model of the hydrosphere and derive steady-state solutions to the water partitioning on terrestrial planets. Critically, hydrostatic seafloor pressure is proportional to surface gravity, so super-Earths with a deep water cycle will tend to store more water in the mantle. We conclude that a tectonically active terrestrial planet of any mass can maintain exposed continents if its water mass fraction is less than ∼0.2%, dramatically increasing the odds that super-Earths are habitable. The greatest source of uncertainty in our study is Earth's current mantle water inventory: the greater its value, the more robust planets are to inundation. Lastly, we discuss how future missions can test our hypothesis by mapping the oceans and continents of massive terrestrial planets. [ABSTRACT FROM AUTHOR]

We report Warm Spitzer full-orbit phase observations of WASP-12b at 3.6 and 4.5 µm. This extremely inflated hot Jupiter is thought to be overflowing its Roche lobe, undergoing mass loss and accretion onto its host star, and has been claimed to have a C/O ratio in excess of unity. We are able to measure the transit depths, eclipse depths, thermal and ellipsoidal phase variations at both wavelengths. The large-amplitude phase variations, combined with the planet's previously measured dayside spectral energy distribution, are indicative of non-zero Bond albedo and very poor day--night heat redistribution. The transit depths in the mid-infrared--(Rp/R*)² = 0.0123(3) and 0.0111(3) at 3.6 and 4.5 µm, respectively--indicate that the atmospheric opacity is greater at 3.6 than at 4.5 µm, in disagreement with model predictions, irrespective of C/O ratio. The secondary eclipse depths are consistent with previous studies: Fday/F* = 0.0038(4) and 0.0039(3) at 3.6 and 4.5 µm, respectively. We do not detect ellipsoidal variations at 3.6 µm, but our parameter uncertainties--estimated via prayer-bead Monte Carlo--keep this non-detection consistent with model predictions. At 4.5 µm, on the other hand, we detect ellipsoidal variations that are much stronger than predicted. If interpreted as a geometric effect due to the planet's elongated shape, these variations imply a 3:2 ratio for the planet's longest:shortest axes and a relatively bright day-night terminator. If we instead presume that the 4.5 µm ellipsoidal variations are due to uncorrected systematic noise and we fix the amplitude of the variations to zero, the best-fit 4.5 µm transit depth becomes commensurate with the 3.6 µm depth, within the uncertainties. The relative transit depths are then consistent with a solar composition and short scale height at the terminator. Assuming zero ellipsoidal variations also yields a much deeper 4.5 µm eclipse depth, consistent with a solar composition and modest temperature inversion. We suggest future observations that could distinguish between these two scenarios. [ABSTRACT FROM AUTHOR]

Thermal phase variations of short-period planets indicate that they are not spherical cows: day-to-night temperature contrasts range from hundreds to thousands of degrees, rivaling their vertical temperature contrasts. Nonetheless, the emergent spectra of short-period planets have typically been fit using one-dimensional (1D) spectral retrieval codes that only account for vertical temperature gradients. The popularity of 1D spectral retrieval codes is easy to understand: they are robust and have a rich legacy in solar system atmospheric studies. Exoplanet researchers have recently introduced multidimensional retrieval schemes to interpret the spectra of short-period planets, but these codes are necessarily more complex and computationally expensive than their 1D counterparts. In this paper we present an alternative: phase-dependent spectral observations are inverted to produce longitudinally resolved spectra that can then be fit using standard 1D spectral retrieval codes. We test this scheme on the iconic phase-resolved spectra of WASP-43b and on simulated observations for the James Webb Space Telescope (JWST) using the open-source Pyrat Bay 1D spectral retrieval framework. Notably, we take the model complexity of the simulations one step further from previous studies by allowing for longitudinal variations in composition in addition to temperature. We show that performing 1D spectral retrieval on longitudinally resolved spectra is more accurate than applying 1D spectral retrieval codes to disk-integrated emission spectra, even though this is identical in terms of computational load. We find that for the extant Hubble and Spitzer observations of WASP-43b, the difference between the two approaches is negligible, but JWST phase measurements should be treated with longitudinally re solved spect ral retrieval (ReSpect). [ABSTRACT FROM AUTHOR]

We have developed an open-source pipeline for the analysis of Spitzer /IRAC channel 1 and 2 time-series photometry, incorporating some of the most popular decorrelation methods. We applied this pipeline to new phase curve observations of ultra-hot Jupiters MASCARA-1b and KELT-16b, and we performed the first comprehensive reanalysis of 15 phase curves. We find that MASCARA-1b and KELT-16b have phase offsets of |$6^{+11}_{-11}~^{\circ }$| W and |$38^{+16}_{-15}~^{\circ }$| W, dayside temperatures of |$2952^{+100}_{-97}$|   and |$3070^{+160}_{-150}$|  K, and nightside temperatures of |$1300^{+340}_{-340}$| and |$1900^{+430}_{-440}$|  K, respectively. We confirm a strong correlation between dayside and irradiation temperatures with a shallower dependence for nightside temperature. We also find evidence that the normalized phase curve amplitude (peak-to-trough divided by eclipse depth) is correlated with stellar effective temperature. In addition, while our different models often retrieve similar parameters, significant differences occasionally arise between them, as well as between our preferred model and the literature values. Nevertheless, our preferred models are consistent with published phase offsets to within −8 ± 21 degrees (−1.6 ± 3.2 sigma), and normalized phase curve amplitudes are on average reproduced to within −0.01 ± 0.24 (−0.1 ± 1.6 sigma). Finally, we find that BLISS performs best in most cases, but not all; we therefore recommend future analyses consider numerous detector models to ensure an optimal fit and to assess model dependencies. [ABSTRACT FROM AUTHOR]

The JWST Early Release Observations (ERO) included a NIRISS/SOSS (0.6–2.8  μ m) transit of the ∼ 850 K Saturn-mass exoplanet HAT-P-18 b. Initial analysis of these data reported detections of water, escaping helium and haze. However, active K dwarfs like HAT-P-18 possess surface heterogeneities – star-spots and faculae – that can complicate the interpretation of transmission spectra, and indeed, a spot-crossing event is present in HAT-P-18 b's NIRISS/SOSS light curves. Here, we present an extensive reanalysis and interpretation of the JWST ERO transmission spectrum of HAT-P-18 b, as well as HST /WFC3 and Spitzer /IRAC transit observations. We detect H2O (12.5σ), CO2 (7.3σ), a cloud deck (7.4σ), and unocculted star-spots (5.8σ), alongside hints of Na (2.7σ). We do not detect the previously reported CH4 (log  CH4 < −6 to 2σ). We obtain excellent agreement between three independent retrieval codes, which find a sub-solar H2O abundance (log  H2O ≈ −4.4 ± 0.3). However, the inferred CO2 abundance (log  CO2 ≈ −4.8 ± 0.4) is significantly super-solar and requires further investigation into its origin. We also introduce new stellar heterogeneity considerations by fitting for the active regions' surface gravities – a proxy for the effects of magnetic pressure. Finally, we compare our JWST inferences to those from HST /WFC3 and Spitzer /IRAC. Our results highlight the exceptional promise of simultaneous planetary atmosphere and stellar heterogeneity constraints in the era of JWST and demonstrate that JWST transmission spectra may warrant more complex treatments of the transit light source effect. [ABSTRACT FROM AUTHOR]

NASA is engaged in planning for a Habitable Worlds Observatory (HabWorlds ), a coronagraphic space mission to detect rocky planets in habitable zones and establish their habitability. Surface liquid water is central to the definition of planetary habitability. Photometric and polarimetric phase curves of starlight reflected by an exoplanet can reveal ocean glint, rainbows, and other phenomena caused by scattering by clouds or atmospheric gas. Direct imaging missions are optimized for planets near quadrature, but HabWorlds ' coronagraph may obscure the phase angles where such optical features are strongest. The range of accessible phase angles for a given exoplanet will depend on the planet's orbital inclination and/or the coronagraph's inner working angle (IWA). We use a recently created catalog relevant to HabWorlds of 164 stars to estimate the number of exo-Earths that could be searched for ocean glint, rainbows, and polarization effects due to Rayleigh scattering. We find that the polarimetric Rayleigh scattering peak is accessible in most of the exo-Earth planetary systems. The rainbow due to water clouds at phase angles of ∼20○ − 60○ would be accessible with HabWorlds  for a planet with an Earth equivalent instellation in ∼46 systems, while the ocean glint signature at phase angles of ∼130○ − 170○ would be accessible in ∼16 systems, assuming an IWA = 62 mas (3λ/ D). Improving the IWA = 41 mas (2λ/ D) increases accessibility to rainbows and glints by factors of approximately 2 and 3, respectively. By observing these scattering features, HabWorlds could detect a surface ocean and water cycle, key indicators of habitability. [ABSTRACT FROM AUTHOR]

The future is now – after its long-awaited launch in 2021 December, JWST began science operations in 2022 July and is already revolutionizing exoplanet astronomy. The Early Release Observations (ERO) programme was designed to provide the first images and spectra from JWST , covering a multitude of science cases and using multiple modes of each on-board instrument. Here, we present transmission spectroscopy observations of the hot-Saturn WASP-96 b with the Single Object Slitless Spectroscopy (SOSS) mode of the Near Infrared Imager and Slitless Spectrograph, observed as part of the ERO programme. As the SOSS mode presents some unique data reduction challenges, we provide an in-depth walk-through of the major steps necessary for the reduction of SOSS data: including background subtraction, correction of 1/ f noise, and treatment of the trace order overlap. We furthermore offer potential routes to correct for field star contamination, which can occur due to the SOSS mode's slitless nature. By comparing our extracted transmission spectrum with grids of atmosphere models, we find an atmosphere metallicity between 1× and 5× solar, and a solar carbon-to-oxygen ratio. Moreover, our models indicate that no grey cloud deck is required to fit WASP-96 b's transmission spectrum, but find evidence for a slope shortward of 0.9 μm, which could either be caused by enhanced Rayleigh scattering or the red wing of a pressure-broadened Na feature. Our work demonstrates the unique capabilities of the SOSS mode for exoplanet transmission spectroscopy and presents a step-by-step reduction guide for this new and exciting instrument. [ABSTRACT FROM AUTHOR]

Precise measurements of chemical abundances in planetary atmospheres are necessary to constrain the formation histories of exoplanets. A recent study of WASP-127 b, a close-in puffy sub-Saturn orbiting its solar-type host star in 4.2 d, using Hubble Space Telescope (HST) and Spitzer revealed a feature-rich transmission spectrum with strong excess absorption at 4.5 μm. However, the limited spectral resolution and coverage of these instruments could not distinguish between CO and/or CO2 absorption causing this signal, with both low and high C/O ratio scenarios being possible. Here we present near-infrared (0.9–2.5 μm) transit observations of WASP-127 b using the high-resolution SPIRou spectrograph, with the goal to disentangle CO from CO2 through the 2.3 μm CO band. With SPIRou, we detect H2O at a t -test significance of 5.3 σ and observe a tentative (3σ) signal consistent with OH absorption. From a joint SPIRou + HST + Spitzer retrieval analysis, we rule out a CO-rich scenario by placing an upper limit on the CO abundance of log10[CO] <−4.0, and estimate a log10[CO2] of −3.7 |$^{+0.8}_{-0.6}$| , which is the level needed to match the excess absorption seen at 4.5 μm. We also set abundance constraints on other major C-, O-, and N-bearing molecules, with our results favouring low C/O (0.10 |$^{+0.10}_{-0.06}$|), disequilibrium chemistry scenarios. We further discuss the implications of our results in the context of planet formation. Additional observations at high and low resolution will be needed to confirm these results and better our understanding of this unusual world. [ABSTRACT FROM AUTHOR]

In order to understand the climate on terrestrial planets orbiting nearby Sun-like stars, one would like to know their thermal inertia. We use a global climate model to simulate the thermal phase variations of Earth analogs and test whether these data could distinguish between planets with different heat storage and heat transport characteristics. In particular, we consider a temperate climate with polar ice caps (like the modern Earth) and a snowball state where the oceans are globally covered in ice. We first quantitatively study the periodic radiative forcing from, and climatic response to, rotation, obliquity, and eccentricity. Orbital eccentricity and seasonal changes in albedo cause variations in the global-mean absorbed flux. The responses of the two climates to these global seasons indicate that the temperate planet has 3× the bulk heat capacity of the snowball planet due to the presence of liquid water oceans. The obliquity seasons in the temperate simulation are weaker than one would expect based on thermal inertia alone; this is due to cross-equatorial oceanic and atmospheric energy transport. Thermal inertia and cross-equatorial heat transport have qualitatively different effects on obliquity seasons, insofar as heat transport tends to reduce seasonal amplitude without inducing a phase lag. For an Earth-like planet, however, this effect is masked by the mixing of signals from low thermal inertia regions (sea ice and land) with that from high thermal inertia regions (oceans), which also produces a damped response with small phase lag. We then simulate thermal light curves as they would appear to a high-contrast imaging mission (TPF-I/Darwin). In order of importance to the present simulations, which use modern-Earth orbital parameters, the three drivers of thermal phase variations are (1) obliquity seasons, (2) diurnal cycle, and (3) global seasons. Obliquity seasons are the dominant source of phase variations for most viewing angles. A pole-on observer would measure peak-to-trough amplitudes of 13% and 47% for the temperate and snowball climates, respectively. Diurnal heating is important for equatorial observers (∼5% phase variations), because the obliquity effects cancel to first order from that vantage. Finally, we compare the prospects of optical versus thermal direct imaging missions for constraining the climate on exoplanets and conclude that while zero- and one-dimensional models are best served by thermal measurements, second-order models accounting for seasons and planetary thermal inertia would require both optical and thermal observations. [ABSTRACT FROM AUTHOR]

It is likely that unambiguous habitable zone terrestrial planets of unknown water content will soon be discovered. Water content helps determine surface land fraction, which influences planetary weathering behavior. This is important because the silicate-weathering feedback determines the width of the habitable zone in space and time. Here a low-order model of weathering and climate, useful for gaining qualitative understanding, is developed to examine climate evolution for planets of various land-ocean fractions. It is pointed out that, if seafloor weathering does not depend directly on surface temperature, there can be no weathering-climate feedback on a waterworld. This would dramatically narrow the habitable zone of a waterworld. Results from our model indicate that weathering behavior does not depend strongly on land fraction for partially ocean-covered planets. This is powerful because it suggests that previous habitable zone theory is robust to changes in land fraction, as long as there is some land. Finally, a mechanism is proposed for a waterworld to prevent complete water loss during a moist greenhouse through rapid weathering of exposed continents. This process is named a “waterworld self-arrest,” and it implies that waterworlds can go through a moist greenhouse stage and end up as planets like Earth with partial ocean coverage. This work stresses the importance of surface and geologic effects, in addition to the usual incident stellar flux, for habitability. [ABSTRACT FROM AUTHOR]

A planet's emission spectrum contains information about atmospheric composition and structure. We compare the Bayesian Information Criterion (BIC) of blackbody fits and idealized spectral retrieval fits for the 44 planets with published eclipse measurements in multiple thermal wavebands, mostly obtained with the Spitzer Space Telescope. The evidence for spectral features depends on eclipse depth uncertainties. Spitzer has proven capable of eclipse precisions better than 10−4 when multiple eclipses are analysed simultaneously, but this feat has only been performed four times. It is harder to self-calibrate photometry when a single occultation is reduced and analysed in isolation; we find that such measurements have not passed the test of repeatability. Single-eclipse measurements either have an uncertainty floor of 5 × 10−4, or their uncertainties have been underestimated by a factor of 3. If one adopts these empirical uncertainties for single-eclipse measurements, then the evidence for molecular features all but disappears: blackbodies have better BIC than spectral retrieval for all planets, save HD 189733b, and the few planets poorly fit by blackbodies are also poorly fit by self-consistent radiative transfer models. This suggests that the features in extant broad-band emission spectra are due to astrophysical and instrumental noise rather than molecular bands. Claims of stratospheric inversions, disequilibrium chemistry, and high C/O ratios based solely on photometry are premature. We recommend that observers be cautious of error estimates from self-calibration of small data sets, and that modellers compare the evidence for spectral models to that of simpler models such as blackbodies. [ABSTRACT FROM AUTHOR]

The ESA-Ariel mission will include a tier dedicated to exoplanet phase curves corresponding to ∼ 10 % of the science time. We present here the current observing strategy for studying exoplanet phase curves with Ariel. We define science questions, requirements and a list of potential targets. We also estimate the precision of phase curve reconstruction and atmospheric retrieval using simulated phase curves. Based on this work, we found that full-orbit phase variations for 35-40 exoplanets could be observed during the 3.5-yr mission. This statistical sample would provide key constraints on atmospheric dynamics, composition, thermal structure and clouds of warm exoplanets, complementary to the scientific yield from spectroscopic transits/eclipses measurements. [ABSTRACT FROM AUTHOR]

Reflected light photometry of terrestrial exoplanets could reveal the presence of oceans and continents, hence placing direct constraints on the current and long-term habitability of these worlds. Inferring the albedo map of a planet from its observed light curve is challenging because different maps may yield indistinguishable light curves. This degeneracy is aggravated by changing clouds. It has previously been suggested that disc-integrated photometry spanning multiple days could be combined to obtain a cloud-free surface map of an exoplanet. We demonstrate this technique as part of a Bayesian retrieval by simultaneously fitting for the fixed surface map of a planet and the time-variable overlying clouds. We test this approach on synthetic data and then apply it to real disc-integrated observations of the Earth. We find that 8 d of continuous synthetic observations are sufficient to reconstruct a faithful low-resolution surface albedo map, without needing to make assumptions about cloud physics. For light curves with negligible photometric uncertainties, the minimal top-of-atmosphere albedo at a location is a good estimate of its surface albedo. When applied to observations from the Earth Polychromating Imaging Camera aboard the Deep Space Climate Observatory spacecraft, our approach removes only a small fraction of clouds. We attribute this difficulty to the full-phase geometry of observations combined with the short correlation length for Earth clouds. For exoplanets with Earth-like climatology, it may be hard to do much better than a cloud-averaged map. We surmise that cloud removal will be most successful for exoplanets imaged near quarter phase that harbour large cloud systems. [ABSTRACT FROM AUTHOR]

In this paper, we describe a uniform analysis of eight transits and eleven secondary eclipses of the extrasolar planet GJ 436b obtained in the 3.6, 4.5, and 8.0 μm bands using the IRAC instrument on the Spitzer Space Telescope between UT 2007 June 29 and UT 2009 February 4. We find that the best-fit transit depths for visits in the same bandpass can vary by as much as 8% of the total (4.7σ significance) from one epoch to the next. Although we cannot entirely rule out residual detector effects or a time-varying, high-altitude cloud layer in the planet's atmosphere as the cause of these variations, we consider the occultation of active regions on the star in a subset of the transit observations to be the most likely explanation. We find that for the deepest 3.6 μm transit the in-transit data have a higher standard deviation than the out-of-transit data, as would be expected if the planet occulted a star spot. We also compare all published transit observations for this object and find that transits observed in the infrared typically have smaller timing offsets than those observed in visible light. In this case, the three deepest Spitzer transits are all measured within a period of five days, consistent with a single epoch of increased stellar activity. We reconcile the presence of magnetically active regions with the lack of significant visible or infrared flux variations from the star by proposing that the star's spin axis is tilted with respect to our line of sight and that the planet's orbit is therefore likely to be misaligned. In contrast to the results reported by Beaulieu et al., we find no convincing evidence for methane absorption in the planet's transmission spectrum. If we exclude the transits that we believe to be most affected by stellar activity, we find that we prefer models with enhanced CO and reduced methane, consistent with GJ 436b's dayside composition from Stevenson et al. It is also possible that all transits are significantly affected by this activity, in which case it may not be feasible to characterize the planet's transmission spectrum using broadband photometry obtained over multiple epochs. These observations serve to illustrate the challenges associated with transmission spectroscopy of planets orbiting late-type stars; we expect that other systems, such as GJ 1214, may display comparably variable transit depths. We compare the limb-darkening coefficients predicted by PHOENIX and ATLAS stellar atmosphere models and discuss the effect that these coefficients have on the measured planet-star radius ratios given GJ 436b's near-grazing transit geometry. Our measured 8 μm secondary eclipse depths are consistent with a constant value, and we place a 1σ upper limit of 17% on changes in the planet's dayside flux in this band. These results are consistent with predictions from general circulation models for this planet, which find that the planet's dayside flux varies by a few percent or less in the 8 μm band. Averaging over the eleven visits gives us an improved estimate of 0.0452% ± 0.0027% for the secondary eclipse depth; we also examine residuals from the eclipse ingress and egress and place an upper limit on deviations caused by a non-uniform surface brightness for GJ 436b. We combine timing information from our observations with previously published data to produce a refined orbital ephemeris and determine that the best-fit transit and eclipse times are consistent with a constant orbital period. We find that the secondary eclipse occurs at a phase of 0.58672 ± 0.00017, corresponding to e cos(ω) = 0.13754 ± 0.00027, where e is the planet's orbital eccentricity and ω is the longitude of pericenter. We also present improved estimates for other system parameters, including the orbital inclination, a/R⋆, and the planet-star radius ratio. [ABSTRACT FROM AUTHOR]

Using a large (14 857), homogenously selected sample of cluster galaxies identified in the Sloan Digital Sky Survey Data Release 4, we investigate the impact of cluster membership and local density on the stellar mass–gas phase metallicity relation (MZR). We show that stellar metallicities are not suitable for this work, being relatively insensitive to subtle changes in the MZR. Accurate nebular abundances can be obtained for 1318 cluster galaxies in our sample and we show that these galaxies are drawn from clusters that are fully representative of the parent sample in terms of mass, size, velocity dispersion and richness. By comparing the MZR of the cluster galaxies with a sample of control galaxies matched in mass, redshift, fibre covering fraction and rest-frame colour cluster galaxies are found to have, on average, higher metallicities by up to 0.04 dex. The magnitude of this offset does not depend strongly on galactic half-light radius or cluster properties such as velocity dispersion or cluster mass. The effect of local density on the MZR is investigated, using the presence of a near neighbour and both two- and three-dimensional density estimators. For all three metrics, it is found that the cluster galaxies in locally rich environments have higher median metallicities by up to ∼0.05 dex than those in locally poor environments (or without a near neighbour). Control (non-cluster) galaxies at locally high densities exhibit similar metal enhancements. Taken together, these results show that galaxies in clusters are, on average, slightly more metal rich than the field, but that this effect is driven by local overdensity and not simply cluster membership. [ABSTRACT FROM AUTHOR]

Rotational modulations of emission spectra in brown dwarf and exoplanet atmospheres show that clouds are often distributed non-uniformly in these ultracool atmospheres. The spatial heterogeneity in cloud distribution demonstrates the impact of atmospheric dynamics on cloud formation and evolution. In this study, we update the Hubble Space Telescope (HST) time-series data analysis of the previously reported rotational modulations of WISEP J004701+680352—an unusually red late-L brown dwarf with a spectrum similar to that of the directly imaged planet HR 8799e. We construct a self-consistent spatially heterogeneous cloud model to explain the HST and the Spitzer time-series observations, as well as the time-averaged spectra of WISE 0047. In the heterogeneous cloud model, a cloud thickness variation of around one pressure scale height explains the wavelength dependence in the HST near-IR spectral variability. By including disequilibrium CO/CH4 chemistry, our models also reproduce the redder color of WISE 0047 compared to that of field brown dwarfs. We discuss the impact of vertical cloud structure on atmospheric profile and estimate the minimum eddy diffusivity coefficient for other objects with redder colors. Our data analysis and forward modeling results demonstrate that time-series spectrophotometry with a broad wavelength coverage is a powerful tool for constraining heterogeneous atmospheric structure. [ABSTRACT FROM AUTHOR]

The exoplanet WASP-12b is the prototype for the emerging class of ultrahot, Jupiter-mass exoplanets. Past models have predicted – and near-ultraviolet observations have shown – that this planet is losing mass. We present an analysis of two sets of 3.6 and 4.5  |$\mu \rm{m}$| Spitzer phase curve observations of the system which show clear evidence of infrared radiation from gas stripped from the planet, and the gas appears to be flowing directly toward or away from the host star. This accretion signature is only seen at 4.5  |$\mu \rm{m}$|⁠ , not at 3.6  |$\mu \rm{m}$|⁠ , which is indicative either of CO emission at the longer wavelength or blackbody emission from cool, ≲600 K gas. It is unclear why WASP-12b is the only ultrahot Jupiter to exhibit this mass-loss signature, but perhaps WASP-12b's orbit is decaying as some have claimed, while the orbits of other exoplanets may be more stable; alternatively, the high-energy irradiation from WASP-12A may be stronger than the other host stars. We also find evidence for phase offset variability at the level of 6.4σ (46.2°) at 3.6  |$\mu \rm{m}$|⁠. [ABSTRACT FROM AUTHOR]

We report secondary eclipse photometry of the hot Jupiter XO-3b in the 4.5 μm band taken with the Infrared Array Camera on the Spitzer Space Telescope. We measure individual eclipse depths and center of eclipse times for a total of 12 secondary eclipses. We fit these data simultaneously with two transits observed in the same band in order to obtain a global best-fit secondary eclipse depth of 0.1580% ± 0.0036% and a center of eclipse phase of 0.67004 ± 0.00013. We assess the relative magnitude of variations in the dayside brightness of the planet by measuring the size of the residuals during ingress and egress from fitting the combined eclipse light curve with a uniform disk model and place an upper limit of 0.05%. The new secondary eclipse observations extend the total baseline from one and a half years to nearly three years, allowing us to place an upper limit on the periastron precession rate of 2.9 × 10–3 deg day–1— the tightest constraint to date on the periastron precession rate of a hot Jupiter. We use the new transit observations to calculate improved estimates for the system properties, including an updated orbital ephemeris. We also use the large number of secondary eclipses to obtain the most stringent limits to date on the orbit-to-orbit variability of an eccentric hot Jupiter and demonstrate the consistency of multiple-epoch Spitzer observations. [ABSTRACT FROM AUTHOR]

We present secondary eclipse photometry of the extrasolar planet WASP-5b taken in the 3.6 and 4.5μm bands with the Spitzer Space Telescope's Infrared Array Camera as part of the extended warm mission. By estimating the depth of the secondary eclipse in these two bands we can place constraints on the planet's atmospheric pressure–temperature profile and chemistry. We measure secondary eclipse depths of 0.197% ± 0.028% and 0.237% ± 0.024% in the 3.6μm and 4.5μm bands, respectively. For the case of a solar-composition atmosphere and chemistry in local thermal equilibrium, our observations are bestmatched by models showing a hot dayside and, depending on our choice of model, a weak thermal inversion or no inversion at all.We measure a mean offset from the predicted center of eclipse of 3.7 ± 1.8 minutes, corresponding to e cos ω = 0.0025 ± 0.0012 and consistent with a circular orbit. We conclude that the planet's orbit is unlikely to have been perturbed by interactions with another body in the system as claimed by Fukui et al. [ABSTRACT FROM AUTHOR]

We present Warm Spitzer/IRAC secondary eclipse time series photometry of three short-period transiting exoplanets, HAT-P-3b, HAT-P-4b and HAT-P-12b, in both the available 3.6 and 4.5μm bands. HAT-P-3b and HAT-P-4b are Jupiter-mass objects orbiting an early K and an early G dwarf star, respectively. For HAT-P-3b we find eclipse depths of 0.112%-0.030%+0.015%(3.6μm) and 0.094%-0.009%+0.016%(4.5μm). The HAT-P-4b values are 0.142%-0.016%+0.014%(3.6μm) and 0.122%-0.014%+0.012%(4.5μm). The two planets' photometry is consistent with inefficient heat redistribution from their day to night sides (and low albedos), but it is inconclusive about possible temperature inversions in their atmospheres. HAT-P-12b is a Saturn-mass planet and is one of the coolest planets ever observed during secondary eclipse, along with the hot Neptune GJ 436b and the hot Saturn WASP-29b. We are able to place 3σ upper limits on the secondary eclipse depth of HAT-P-12b in both wavelengths: <0.042% (3.6μm) and <0.085% (4.5μm).We discuss these results in the context of the Spitzer secondary eclipse measurements of GJ 436b and WASP-29b. It is possible that we do not detect the eclipses of HAT-P-12b due to high eccentricity, but find that weak planetary emission in these wavelengths is a more likely explanation. We place 3σ upper limits on the ∣e cos ω∣ quantity (where e is eccentricity and ω is the argument of periapsis) for HAT-P-3b (<0.0081) and HAT-P-4b (<0.0042), based on the secondary eclipse timings. [ABSTRACT FROM AUTHOR]

We present the first secondary eclipse and phase curve observations for the highly eccentric hot Jupiter HAT-P-2b in the 3.6, 4.5, 5.8, and 8.0μm bands of the Spitzer Space Telescope. The 3.6 and 4.5μm data sets span an entire orbital period of HAT-P-2b (P = 5.6334729 d), making them the longest continuous phase curve observations obtained to date and the first full-orbit observations of a planet with an eccentricity exceeding 0.2. We present an improved non-parametric method for removing the intrapixel sensitivity variations in Spitzer data at 3.6 and 4.5μm that robustly maps position-dependent flux variations. We find that the peak in planetary flux occurs at 4.39 ± 0.28, 5.84 ± 0.39, and 4.68 ± 0.37 hr after periapse passage with corresponding maxima in the planet/star flux ratio of 0.1138% ± 0.0089%, 0.1162% ± 0.0080%, and 0.1888% ± 0.0072% in the 3.6, 4.5, and 8.0μm bands, respectively. Our measured secondary eclipse depths of 0.0996%±0.0072%, 0.1031%±0.0061%, 0.071%+0.029% -0.013%, and 0.1392%±0.0095% in the 3.6, 4.5, 5.8, and 8.0μm bands, respectively, indicate that the planet cools significantly from its peak temperature before we measure the dayside flux during secondary eclipse. We compare our measured secondary eclipse depths to the predictions from a one-dimensional radiative transfer model, which suggests the possible presence of a transient day side inversion in HAT-P-2b's atmosphere near periapse.We also derive improved estimates for the system parameters, including its mass, radius, and orbital ephemeris. Our simultaneous fit to the transit, secondary eclipse, and radial velocity data allows us to determine the eccentricity (e = 0.50910 ± 0.00048) and argument of periapse (ω = 188. ° 09 ± 0. ° 39) of HAT-P-2b's orbit with a greater precision than has been achieved for any other eccentric extrasolar planet. We also find evidence for a long-term linear trend in the radial velocity data. This trend suggests the presence of another substellar companion in the HAT-P-2 system, which could have caused HAT-P-2b to migrate inward to its present-day orbit via the Kozai mechanism. [ABSTRACT FROM AUTHOR]

We analyze Warm Spitzer/Infrared Array Camera observations of the secondary eclipses of three planets, XO-4b, HAT-P-6b, and HAT-P-8b. We measure secondary eclipse amplitudes at 3.6 µm and 4.5 µm for each target. XO-4b exhibits a stronger eclipse depth at 4.5 µm than at 3.6 µm, which is consistent with the presence of a temperature inversion. HAT-P-8b shows a stronger eclipse amplitude at 3.6 µm and is best described by models without a temperature inversion. The eclipse depths of HAT-P-6b can be fitted with models with a small or no temperature inversion. We consider our results in the context of a postulated relationship between stellar activity and temperature inversion and a relationship between irradiation level and planet dayside temperature, as discussed by Knutson et al. and Cowan & Agol, respectively. Our results are consistent with these hypotheses, but do not significantly strengthen them. To measure accurate secondary eclipse central phases, we require accurate ephemerides. We obtain primary transit observations and supplement them with publicly available observations to update the orbital ephemerides of the three planets. Based on the secondary eclipse timing, we set upper boundaries for e cos(ω) for HAT-P-6b, HAT-P-8b, and XO-4b and find that the values are consistent with circular orbits. [ABSTRACT FROM AUTHOR]

ABSTRACT The influence of a galaxy's environment on its evolution has been studied and compared extensively in the literature, although differing techniques are often used to define environment. Most methods fall into two broad groups: those that use nearest neighbours to probe the underlying density field and those that use fixed apertures. The differences between the two inhibit a clean comparison between analyses and leave open the possibility that, even with the same data, different properties are actually being measured. In this work, we apply 20 published environment definitions to a common mock galaxy catalogue constrained to look like the local Universe. We find that nearest-neighbour-based measures best probe the internal densities of high-mass haloes, while at low masses the interhalo separation dominates and acts to smooth out local density variations. The resulting correlation also shows that nearest-neighbour galaxy environment is largely independent of dark matter halo mass. Conversely, aperture-based methods that probe superhalo scales accurately identify high-density regions corresponding to high-mass haloes. Both methods show how galaxies in dense environments tend to be redder, with the exception of the largest apertures, but these are the strongest at recovering the background dark matter environment. We also warn against using photometric redshifts to define environment in all but the densest regions. When considering environment, there are two regimes: the 'local environment' internal to a halo best measured with nearest neighbour and 'large-scale environment' external to a halo best measured with apertures. This leads to the conclusion that there is no universal environment measure and the most suitable method depends on the scale being probed. [ABSTRACT FROM AUTHOR]

‘Hot Jupiter’ extrasolar planets are expected to be tidally locked because they are close (<0.05 astronomical units, where 1 au is the average Sun–Earth distance) to their parent stars, resulting in permanent daysides and nightsides. By observing systems where the planet and star periodically eclipse each other, several groups have been able to estimate the temperatures of the daysides of these planets. A key question is whether the atmosphere is able to transport the energy incident upon the dayside to the nightside, which will determine the temperature at different points on the planet’s surface. Here we report observations of HD 189733, the closest of these eclipsing planetary systems, over half an orbital period, from which we can construct a ‘map’ of the distribution of temperatures. We detected the increase in brightness as the dayside of the planet rotated into view. We estimate a minimum brightness temperature of 973 ± 33 K and a maximum brightness temperature of 1,212 ± 11 K at a wavelength of 8 μm, indicating that energy from the irradiated dayside is efficiently redistributed throughout the atmosphere, in contrast to a recent claim for another hot Jupiter. Our data indicate that the peak hemisphere-integrated brightness occurs 16 ± 6° before opposition, corresponding to a hotspot shifted east of the substellar point. The secondary eclipse (when the planet moves behind the star) occurs 120 ± 24 s later than predicted, which may indicate a slightly eccentric orbit. [ABSTRACT FROM AUTHOR]

The hot Jupiter HD 209458b is particularly amenable to detailed study as it is among the brightest transiting exoplanet systems currently known (V-mag = 7.65; K-mag = 6.308) and has a large planet-to-star contrast ratio. HD 209458b is predicted to be in synchronous rotation about its host star with a hot spot that is shifted eastward of the substellar point by superrotating equatorial winds. Here we present the first full-orbit observations of HD 209458b, in which its 4.5 μm emission was recorded with Spitzer/IRAC. Our study revises the previous 4.5 μm measurement of HD 209458b's secondary eclipse emission downward by ∼35% to , changing our interpretation of the properties of its dayside atmosphere. We find that the hot spot on the planet's dayside is shifted eastward of the substellar point by 40.°9 ± 6.°0, in agreement with circulation models predicting equatorial superrotation. HD 209458b's dayside (Tbright = 1499 ± 15 K) and nightside (Tbright = 972 ± 44 K) emission indicate a day-to-night brightness temperature contrast smaller than that observed for more highly irradiated exoplanets, suggesting that the day-to-night temperature contrast may be partially a function of the incident stellar radiation. The observed phase curve shape deviates modestly from global circulation model predictions potentially due to disequilibrium chemistry or deficiencies in the current hot CH4 line lists used in these models. Observations of the phase curve at additional wavelengths are needed in order to determine the possible presence and spatial extent of a dayside temperature inversion, as well as to improve our overall understanding of this planet's atmospheric circulation. [ABSTRACT FROM AUTHOR]

Kepler-13Ab (= KOI-13.01) is a unique transiting hot Jupiter. It is one of very few known short-period planets orbiting a hot A-type star, making it one of the hottest planets currently known. The availability of Kepler data allows us to measure the planet's occultation (secondary eclipse) and phase curve in the optical, which we combine with occultations observed by warm Spitzer at 4.5 μm and 3.6 μm and a ground-based occultation observation in the Ks band (2.1 μm). We derive a day-side hemisphere temperature of 2750 ± 160 K as the effective temperature of a black body showing the same occultation depths. Comparing the occultation depths with one-dimensional planetary atmosphere models suggests the presence of an atmospheric temperature inversion. Our analysis shows evidence for a relatively high geometric albedo, Ag = . While measured with a simplistic method, a high Ag is supported also by the fact that the one-dimensional atmosphere models underestimate the occultation depth in the optical. We use stellar spectra to determine the dilution, in the four wide bands where occultation was measured, due to the visual stellar binary companion 1.″15 ± 0.″05 away. The revised stellar parameters measured using these spectra are combined with other measurements, leading to revised planetary mass and radius estimates of Mp = 4.94-8.09 MJ and Rp = 1.406 ± 0.038 RJ. Finally, we measure a Kepler midoccultation time that is 34.0 ± 6.9 s earlier than expected based on the midtransit time and the delay due to light-travel time and discuss possible scenarios. [ABSTRACT FROM AUTHOR]

We present new, full-orbit observations of the infrared phase variations of the canonical hot Jupiter HD 189733b obtained in the 3.6 and 4.5 μm bands using the Spitzer Space Telescope. When combined with previous phase curve observations at 8.0 and 24 μm, these data allow us to characterize the exoplanet's emission spectrum as a function of planetary longitude and to search for local variations in its vertical thermal profile and atmospheric composition. We utilize an improved method for removing the effects of intrapixel sensitivity variations and robustly extracting phase curve signals from these data, and we calculate our best-fit parameters and uncertainties using a wavelet-based Markov Chain Monte Carlo analysis that accounts for the presence of time-correlated noise in our data. We measure a phase curve amplitude of 0.1242% ± 0.0061% in the 3.6 μm band and 0.0982% ± 0.0089% in the 4.5 μm band, corresponding to brightness temperature contrasts of 503 ± 21 K and 264 ± 24 K, respectively. We find that the times of minimum and maximum flux occur several hours earlier than predicted for an atmosphere in radiative equilibrium, consistent with the eastward advection of gas by an equatorial super-rotating jet. The locations of the flux minima in our new data differ from our previous observations at 8 μm, and we present new evidence indicating that the flux minimum observed in the 8 μm is likely caused by an overshooting effect in the 8 μm array. We obtain improved estimates for HD 189733b's dayside planet-star flux ratio of 0.1466% ± 0.0040% in the 3.6 μm band and 0.1787% ± 0.0038% in the 4.5 μm band, corresponding to brightness temperatures of 1328 ± 11 K and 1192 ± 9 K, respectively; these are the most accurate secondary eclipse depths obtained to date for an extrasolar planet. We compare our new dayside and nightside spectra for HD 189733b to the predictions of one-dimensional radiative transfer models from Burrows et al. and conclude that fits to this planet's dayside spectrum provide a reasonably accurate estimate of the amount of energy transported to the night side. Our 3.6 and 4.5 μm phase curves are generally in good agreement with the predictions of general circulation models for this planet from Showman et al., although we require either excess drag or slower rotation rates in order to match the locations of the measured maxima and minima in the 4.5, 8.0, and 24 μm bands. We find that HD 189733b's 4.5 μm nightside flux is 3.3σ smaller than predicted by these models, which assume that the chemistry is in local thermal equilibrium. We conclude that this discrepancy is best explained by vertical mixing, which should lead to an excess of CO and correspondingly enhanced 4.5 μm absorption in this region. This result is consistent with our constraints on the planet's transmission spectrum, which also suggest excess absorption in the 4.5 μm band at the day-night terminator. [ABSTRACT FROM AUTHOR]

We report the results of an analysis of all Spitzer/MIPS 24 μm observations of HD 209458b, one of the touchstone objects in the study of irradiated giant planet atmospheres. Altogether, we analyze two and a half transits, three eclipses, and a 58 hr near-continuous observation designed to detect the planet's thermal phase curve. The results of our analysis are: (1) a mean transit depth of 1.484% ± 0.033%, consistent with previous measurements and showing no evidence of variability in transit depth at the 3% level. (2) A mean eclipse depth of 0.338% ± 0.026%, somewhat higher than that previously reported for this system; this new value brings observations into better agreement with models. From this eclipse depth we estimate an average dayside brightness temperature of 1320 ± 80 K; the dayside flux shows no evidence of variability at the 12% level. (3) Eclipses in the system occur 32 ± 129 s earlier than would be expected from a circular orbit, which constrains the orbital quantity ecos ω to be 0.00004 ± 0.00033. This result is fully consistent with a circular orbit and sets an upper limit of 140 m s–1 (3σ) on any eccentricity-induced velocity offset during transit. The phase curve observations (including one of the transits) exhibit an anomalous trend similar to the detector ramp seen in previous Spitzer/IRAC observations; by modeling this ramp we recover the system parameters for this transit. The long-duration photometry which follows the ramp and transit exhibits a gradual ∼0.2% decrease in flux over ∼30 hr. This effect is similar to that seen in pre-launch calibration data taken with the 24 μm array and is better fit by an instrumental model than a model invoking planetary emission. The large uncertainties associated with this poorly understood, likely instrumental effect prevent us from usefully constraining the planet's thermal phase curve. Our observations highlight the need for a thorough understanding of detector-related instrumental effects on long timescales when making the high-precision mid-infrared measurements planned for future missions such as EChO, SPICA, and the James Webb Space Telescope. [ABSTRACT FROM AUTHOR]