-
Differentiable Modeling of Planet and Substellar Atmosphere: High-Resolution Emission, Transmission, and Reflection Spectroscopy with ExoJAX2
Authors:
Hajime Kawahara,
Yui Kawashima,
Shotaro Tada,
Hiroyuki Tako Ishikawa,
Ko Hosokawa,
Yui Kasagi,
Takayuki Kotani,
Kento Masuda,
Stevanus Nuguroho,
Motohide Tamura,
Hibiki Yama,
Daniel Kitzmann,
Nicolas Minesi,
Brett M. Morris
Abstract:
Modeling based on differentiable programming holds great promise for astronomy, as it can employ techniques such as Hamiltonian Monte Carlo, gradient-based optimization, and other machine learning techniques. This new programming paradigm has motivated us to develop the first auto-differentiable spectrum model of exoplanets and brown dwarfs, ExoJAX (Kawahara et al. 2022). ExoJAX is designed to dir…
▽ More
Modeling based on differentiable programming holds great promise for astronomy, as it can employ techniques such as Hamiltonian Monte Carlo, gradient-based optimization, and other machine learning techniques. This new programming paradigm has motivated us to develop the first auto-differentiable spectrum model of exoplanets and brown dwarfs, ExoJAX (Kawahara et al. 2022). ExoJAX is designed to directly calculate cross-sections as functions of temperature and pressure, rather than interpolating tabulated data, to minimize errors in high-dispersion spectra modeling. However, its application was primarily proof-of-concept and limited to narrowband high-dispersion emission spectroscopy. In this paper, we have enhanced the differentiable opacity calculation using a new fast and memory-efficient algorithm, and have developed differentiable radiative transfer schemes, including emission, transmission, and reflection spectroscopy. These enhancements significantly expand the range of applications, as demonstrated through actual atmospheric retrievals: high-dispersion emission spectra of the brown dwarf GL229 B, medium-dispersion transmission spectra of the hot Saturn WASP-39 b from JWST, and high-dispersion reflection spectra of Jupiter. We obtained a C/O ratio for GL229 B consistent with its host star, constrained WASP-39 b's radial velocity from molecular fine structures at original resolution ($R \sim 3,000$), and estimated Jupiter's metallicity consistent with previous studies.
△ Less
Submitted 9 October, 2024;
originally announced October 2024.
-
An atlas of resolved spectral features in the transmission spectrum of WASP-189 b with MAROON-X
Authors:
B. Prinoth,
H. J. Hoeijmakers,
B. M. Morris,
M. Lam,
D. Kitzmann,
E. Sedaghati,
J. V. Seidel,
E. K. H. Lee,
B. Thorsbro,
N. W. Borsato,
Y. C. Damasceno,
S. Pelletier,
A. Seifahrt
Abstract:
Exoplanets in the ultra-hot Jupiter regime provide an excellent laboratory for testing the impact of stellar irradiation on the dynamics and chemical composition of gas giant atmospheres. In this study, we observed two transits of the ultra-hot Jupiter WASP-189 b with MAROON-X/Gemini-North to probe its high-altitude atmospheric layers, using strong absorption lines. We derived posterior probabilit…
▽ More
Exoplanets in the ultra-hot Jupiter regime provide an excellent laboratory for testing the impact of stellar irradiation on the dynamics and chemical composition of gas giant atmospheres. In this study, we observed two transits of the ultra-hot Jupiter WASP-189 b with MAROON-X/Gemini-North to probe its high-altitude atmospheric layers, using strong absorption lines. We derived posterior probability distributions for the planetary and stellar parameters by calculating the stellar spectrum behind the planet at every orbital phase during the transit. This was used to correct the Rossiter-McLaughlin imprint on the transmission spectra. Using differential transmission spectroscopy, we detect strong absorption lines of Ca+, Ba+, Na, H$α$, Mg, Fe, and Fe+, providing an unprecedented and detailed view of the atmospheric chemical composition. Ca+ absorption is particularly well suited for analysis through time-resolved narrow-band spectroscopy, owing to its transition lines formed in high-altitude layers. The spectral absorption lines show no significant blueshifts that would indicate high-altitude day-to-night winds, and further analysis is needed to investigate the implications for atmospheric dynamics. These high signal-to-noise observations provide a benchmark data set for testing high-resolution retrievals and the assumptions of atmospheric models. We also simulate observations of WASP-189 b with ANDES/ELT, and show that ANDES will be highly sensitive to the individual absorption lines of a myriad of elements and molecules, including TiO and CO.
△ Less
Submitted 13 March, 2024;
originally announced March 2024.
-
Observations of scattered light from exoplanet atmospheres
Authors:
Brett M. Morris,
Kevin Heng,
Daniel Kitzmann
Abstract:
Optical phase curves of hot Jupiters can reveal global scattering properties. We implement a Bayesian inference framework for optical phase curves with flux contributions from: reflected light from a potentially inhomogeneous atmosphere, thermal emission, ellipsoidal variations, Doppler beaming, and stellar rotation via a Gaussian process in the time domain. We probe for atmospheric homogeneity an…
▽ More
Optical phase curves of hot Jupiters can reveal global scattering properties. We implement a Bayesian inference framework for optical phase curves with flux contributions from: reflected light from a potentially inhomogeneous atmosphere, thermal emission, ellipsoidal variations, Doppler beaming, and stellar rotation via a Gaussian process in the time domain. We probe for atmospheric homogeneity and time-variability using the reflected light inferences for highly precise Kepler light curves of five hot Jupiters. We also investigate the scattering properties which constrain the most likely condensates in the inhomogeneous atmospheres. Cross validation prefers inhomogeneous albedo distributions for Kepler-7 b and Kepler-41 b, and a weak preference for inhomogeneity for KOI-13 b. None of the five planets exhibit significant variations in geometric albedo on one-year timescales, in agreement with theoretical expectations. We show that analytic reflected light phase curves with isotropic multiple scattering are in excellent agreement with full Rayleigh multiple scattering calculations, allowing for accelerated and analytic inference. In a case study of Kepler-41 b, we identify perovskite, forsterite, and enstatite as possible scattering species consistent with the reflected light phase curves, with condensate particle radii in the range 0.01-0.1 micron.
△ Less
Submitted 24 January, 2024;
originally announced January 2024.
-
A roadmap for the atmospheric characterization of terrestrial exoplanets with JWST
Authors:
TRAPPIST-1 JWST Community Initiative,
:,
Julien de Wit,
Ren� Doyon,
Benjamin V. Rackham,
Olivia Lim,
Elsa Ducrot,
Laura Kreidberg,
Bj�rn Benneke,
Ignasi Ribas,
David Berardo,
Prajwal Niraula,
Aishwarya Iyer,
Alexander Shapiro,
Nadiia Kostogryz,
Veronika Witzke,
Micha�l Gillon,
Eric Agol,
Victoria Meadows,
Adam J. Burgasser,
James E. Owen,
Jonathan J. Fortney,
Franck Selsis,
Aaron Bello-Arufe,
Zo� de Beurs
, et al. (58 additional authors not shown)
Abstract:
Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets. While JWST Cycle 1 observations have started to yield preliminary insights into the planets, they have also revealed that their atmospheric exploration requires a bet…
▽ More
Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets. While JWST Cycle 1 observations have started to yield preliminary insights into the planets, they have also revealed that their atmospheric exploration requires a better understanding of their host star. Here, we propose a roadmap to characterize the TRAPPIST-1 system -- and others like it -- in an efficient and robust manner. We notably recommend that -- although more challenging to schedule -- multi-transit windows be prioritized to mitigate the effects of stellar activity and gather up to twice more transits per JWST hour spent. We conclude that, for such systems, planets cannot be studied in isolation by small programs, but rather need large-scale, jointly space- and ground-based initiatives to fully exploit the capabilities of JWST for the exploration of terrestrial planets.
△ Less
Submitted 22 July, 2024; v1 submitted 24 October, 2023;
originally announced October 2023.
-
Time-resolved transmission spectroscopy of the ultra-hot Jupiter WASP-189 b
Authors:
Bibiana Prinoth,
H. Jens Hoeijmakers,
Stefan Pelletier,
Daniel Kitzmann,
Brett M. Morris,
Andreas Seifahrt,
David Kasper,
Heidi H. Korhonen,
Madeleine Burheim,
Jacob L. Bean,
Bj�rn Benneke,
Nicholas W. Borsato,
Madison Brady,
Simon L. Grimm,
Rafael Luque,
Julian St�rmer,
Brian Thorsbro
Abstract:
Ultra-hot Jupiters are tidally locked with their host stars dividing their atmospheres into a hot dayside and a colder nightside. As the planet moves through transit, different regions of the atmosphere rotate into view revealing different chemical regimes. High-resolution spectrographs can observe asymmetries and velocity shifts, and offer the possibility for time-resolved spectroscopy. In this s…
▽ More
Ultra-hot Jupiters are tidally locked with their host stars dividing their atmospheres into a hot dayside and a colder nightside. As the planet moves through transit, different regions of the atmosphere rotate into view revealing different chemical regimes. High-resolution spectrographs can observe asymmetries and velocity shifts, and offer the possibility for time-resolved spectroscopy. In this study, we search for other atoms and molecules in the planet`s transmission spectrum and investigate asymmetric signals. We analyse and combine eight transits of the ultra-hot Jupiter WASP-189 b taken with the HARPS, HARPS-N, ESPRESSO and MAROON-X high-resolution spectrographs. Using the cross-correlation technique, we search for neutral and ionised atoms, and oxides and compare the obtained signals to model predictions. We report significant detections for H, Na, Mg, Ca, Ca+, Ti, Ti+, TiO, V, Cr, Mn, Fe, Fe+, Ni, Sr, Sr+, and Ba+. Of these, Sr, Sr+, and Ba+ are detected for the first time in the transmission spectrum of WASP-189 b. In addition, we robustly confirm the detection of titanium oxide based on observations with HARPS and HARPS-N using the follow-up observations performed with MAROON-X and ESPRESSO. By fitting the orbital traces of the detected species by means of time-resolved spectroscopy using a Bayesian framework, we infer posterior distributions for orbital parameters as well as lineshapes. Our results indicate that different species must originate from different regions of the atmosphere to be able to explain the observed time dependence of the signals. Throughout the course of the transit, most signal strengths are expected to increase due to the larger atmospheric scale height at the hotter trailing terminator. For some species, however, the signals are instead observed to weaken due to ionisation for atoms and their ions, or the dissociation of molecules on the dayside.
△ Less
Submitted 3 November, 2023; v1 submitted 8 August, 2023;
originally announced August 2023.
-
Investigating the visible phase-curve variability of 55 Cnc e
Authors:
E. A. Meier Vald�s,
B. M. Morris,
B. -O. Demory,
A. Brandeker,
D. Kitzmann,
W. Benz,
A. Deline,
H. -G. Flor�n,
S. G. Sousa,
V. Bourrier,
V. Singh,
K. Heng,
A. Strugarek,
D. J. Bower,
N. J�ggi,
L. Carone,
M. Lendl,
K. Jones,
A. V. Oza,
O. D. S. Demangeon,
Y. Alibert,
R. Alonso,
G. Anglada,
J. Asquier,
T. B�rczy
, et al. (65 additional authors not shown)
Abstract:
55 Cnc e is an ultra-short period super-Earth transiting a Sun-like star. Previous observations in the optical range detected a time-variable flux modulation that is phased with the planetary orbital period, whose amplitude is too large to be explained by reflected light and thermal emission alone. The goal of the study is to investigate the origin of the variability and timescale of the phase-cur…
▽ More
55 Cnc e is an ultra-short period super-Earth transiting a Sun-like star. Previous observations in the optical range detected a time-variable flux modulation that is phased with the planetary orbital period, whose amplitude is too large to be explained by reflected light and thermal emission alone. The goal of the study is to investigate the origin of the variability and timescale of the phase-curve modulation in 55 Cnc e. To this end, we used the CHaracterising ExOPlanet Satellite (CHEOPS), whose exquisite photometric precision provides an opportunity to characterise minute changes in the phase curve from one orbit to the next. CHEOPS observed 29 individual visits of 55 Cnc e between March 2020 and February 2022. Based on these observations, we investigated the different processes that could be at the origin of the observed modulation. In particular, we built a toy model to assess whether a circumstellar torus of dust driven by radiation pressure and gravity might match the observed flux variability timescale. We find that the phase-curve amplitude and peak offset of 55 Cnc e do vary between visits. The sublimation timescales of selected dust species reveal that silicates expected in an Earth-like mantle would not survive long enough to explain the observed phase-curve modulation. We find that silicon carbide, quartz, and graphite are plausible candidates for the circumstellar torus composition because their sublimation timescales are long. The extensive CHEOPS observations confirm that the phase-curve amplitude and offset vary in time.We find that dust could provide the grey opacity source required to match the observations. However, the data at hand do not provide evidence that circumstellar material with a variable grain mass per unit area causes the observed variability. Future observations with the James Webb Space Telescope promise exciting insights into this iconic super-Earth.
△ Less
Submitted 27 July, 2023; v1 submitted 12 July, 2023;
originally announced July 2023.
-
Early Spectroscopy and Dense Circumstellar Medium Interaction in SN 2023ixf
Authors:
K. Azalee Bostroem,
Jeniveve Pearson,
Manisha Shrestha,
David J. Sand,
Stefano Valenti,
Saurabh W. Jha,
Jennifer E. Andrews,
Nathan Smith,
Giacomo Terreran,
Elizabeth Green,
Yize Dong,
Michael Lundquist,
Joshua Haislip,
Emily T. Hoang,
Griffin Hosseinzadeh,
Daryl Janzen,
Jacob E. Jencson,
Vladimir Kouprianov,
Emmy Paraskeva,
Nicolas E. Meza Retamal,
Daniel E. Reichart,
Iair Arcavi,
Alceste Z. Bonanos,
Michael W. Coughlin,
Ross Dobson
, et al. (31 additional authors not shown)
Abstract:
We present the optical spectroscopic evolution of SN~2023ixf seen in sub-night cadence spectra from 1.18 to 14 days after explosion. We identify high-ionization emission features, signatures of interaction with material surrounding the progenitor star, that fade over the first 7 days, with rapid evolution between spectra observed within the same night. We compare the emission lines present and the…
▽ More
We present the optical spectroscopic evolution of SN~2023ixf seen in sub-night cadence spectra from 1.18 to 14 days after explosion. We identify high-ionization emission features, signatures of interaction with material surrounding the progenitor star, that fade over the first 7 days, with rapid evolution between spectra observed within the same night. We compare the emission lines present and their relative strength to those of other supernovae with early interaction, finding a close match to SN~2020pni and SN~2017ahn in the first spectrum and SN~2014G at later epochs. To physically interpret our observations we compare them to CMFGEN models with confined, dense circumstellar material around a red supergiant progenitor from the literature. We find that very few models reproduce the blended \NC{} emission lines observed in the first few spectra and their rapid disappearance thereafter, making this a unique diagnostic. From the best models, we find a mass-loss rate of $10^{-3}-10^{-2}$ \mlunit{}, which far exceeds the mass-loss rate for any steady wind, especially for a red supergiant in the initial mass range of the detected progenitor. These mass-loss rates are, however, similar to rates inferred for other supernovae with early circumstellar interaction. Using the phase when the narrow emission features disappear, we calculate an outer dense radius of circumstellar material $R_\mathrm{CSM, out}\sim5\times10^{14}~\mathrm{cm}$ and a mean circumstellar material density of $ρ=5.6\times10^{-14}~\mathrm{g\,cm^{-3}}$. This is consistent with the lower limit on the outer radius of the circumstellar material we calculate from the peak \Halpha{} emission flux, $R_\text{CSM, out}\gtrsim9\times10^{13}~\mathrm{cm}$.
△ Less
Submitted 12 December, 2023; v1 submitted 16 June, 2023;
originally announced June 2023.
-
Flaring Latitudes in Ensembles of Low Mass Stars
Authors:
Ekaterina Ilin,
Ruth Angus,
Rodrigo Luger,
Brett M. Morris,
Florian U. Jehn
Abstract:
The distribution of small-scale magnetic fields in stellar photospheres is an important ingredient in our understanding of the magnetism of low mass stars. Their spatial distribution connects the field generated in the stellar interior with the outer corona and the large scale field, and thereby affects the space weather of planets. Unfortunately, we lack techniques that can locate them on most lo…
▽ More
The distribution of small-scale magnetic fields in stellar photospheres is an important ingredient in our understanding of the magnetism of low mass stars. Their spatial distribution connects the field generated in the stellar interior with the outer corona and the large scale field, and thereby affects the space weather of planets. Unfortunately, we lack techniques that can locate them on most low-mass stars. One strategy is to localize field concentrations using the flares that occur in their vicinity.
We explore a new method that adapts the spot simulation software fleck to study the modulation of flaring times as a function of active latitude. We use empirical relations to construct flare light curves similar to those available from Kepler and the Transiting Exoplanet Survey Satellite (TESS), search them for flares, and use the waiting times between flares to determine the location of active latitudes.
We find that the mean and standard deviation of the waiting time distribution provide a unique diagnostic of flaring latitudes as a function of the number of active regions. Latitudes are best recovered when stars have three or less active regions that flare repeatedly, and active latitude widths below 20 deg; when either increases, the information about the active latitude location is gradually lost. We demonstrate our technique on a sample of flaring G dwarfs observed with the Kepler satellite, and furthermore suggest that combining ensemble methods for spots and flares could overcome the limitations of each individual technique for the localization of surface magnetic fields.
△ Less
Submitted 3 June, 2023;
originally announced June 2023.
-
Interior heating of rocky exoplanets from stellar flares with application to TRAPPIST-1
Authors:
Alexander Grayver,
Dan J. Bower,
Joachim Saur,
Caroline Dorn,
Brett M. Morris
Abstract:
Many stars of different spectral types with planets in the habitable zone are known to emit flares. Until now, studies that address the long-term impact of stellar flares and associated Coronal Mass Ejections (CMEs) assumed that the planet's interior remains unaffected by interplanetary CMEs, only considering the effect of plasma/UV interactions on the atmosphere of planets. Here, we show that the…
▽ More
Many stars of different spectral types with planets in the habitable zone are known to emit flares. Until now, studies that address the long-term impact of stellar flares and associated Coronal Mass Ejections (CMEs) assumed that the planet's interior remains unaffected by interplanetary CMEs, only considering the effect of plasma/UV interactions on the atmosphere of planets. Here, we show that the magnetic flux carried by flare-associated CMEs results in planetary interior heating by Ohmic dissipation and leads to a variety of interior--exterior interactions. We construct a physical model to study this effect and apply it to the TRAPPIST-1 star whose flaring activity has been constrained by Kepler observations. Our model is posed in a stochastic manner to account for uncertainty and variability in input parameters. Particularly for the innermost planets, our results suggest that the heat dissipated in the silicate mantle is both of sufficient magnitude and longevity to drive geological processes and hence facilitate volcanism and outgassing of the TRAPPIST-1 planets. Furthermore, our model predicts that Joule heating can further be enhanced for planets with an intrinsic magnetic field compared to those without. The associated volcanism and outgassing may continuously replenish the atmosphere and thereby mitigate the erosion of the atmosphere caused by the direct impact of flares and CMEs. To maintain consistency of atmospheric and geophysical models, the impact of stellar flares and CMEs on atmospheres of close-in exoplanetary systems needs to be studied in conjunction with the effect on planetary interiors.
△ Less
Submitted 8 December, 2022; v1 submitted 11 November, 2022;
originally announced November 2022.
-
55 Cancri e's occultation captured with CHEOPS
Authors:
B. -O. Demory,
S. Sulis,
E. Meier Valdes,
L. Delrez,
A. Brandeker,
N. Billot,
A. Fortier,
S. Hoyer,
S. G. Sousa,
K. Heng,
M. Lendl,
A. Krenn,
B. M. Morris,
J. A. Patel,
Y. Alibert,
R. Alonso,
G. Anglada,
T. Barczy,
D. Barrado,
S. C. C. Barros,
W. Baumjohann,
M. Beck,
T. Beck,
W. Benz,
X. Bonfils
, et al. (51 additional authors not shown)
Abstract:
Past occultation and phase-curve observations of the ultra-short period super-Earth 55 Cnc e obtained at visible and infrared wavelengths have been challenging to reconcile with a planetary reflection and emission model. In this study, we analyse a set of 41 occultations obtained over a two-year timespan with the CHEOPS satellite. We report the detection of 55 Cnc e's occultation with an average d…
▽ More
Past occultation and phase-curve observations of the ultra-short period super-Earth 55 Cnc e obtained at visible and infrared wavelengths have been challenging to reconcile with a planetary reflection and emission model. In this study, we analyse a set of 41 occultations obtained over a two-year timespan with the CHEOPS satellite. We report the detection of 55 Cnc e's occultation with an average depth of $12\pm3$ ppm. We derive a corresponding 2-$σ$ upper limit on the geometric albedo of $A_g < 0.55$ once decontaminated from the thermal emission measured by Spitzer at 4.5$μ$m. CHEOPS's photometric performance enables, for the first time, the detection of individual occultations of this super-Earth in the visible and identifies short-timescale photometric corrugations likely induced by stellar granulation. We also find a clear 47.3-day sinusoidal pattern in the time-dependent occultation depths that we are unable to relate to stellar noise, nor instrumental systematics, but whose planetary origin could be tested with upcoming JWST occultation observations of this iconic super-Earth.
△ Less
Submitted 7 November, 2022;
originally announced November 2022.
-
The Mantis Network IV: A titanium cold-trap on the ultra-hot Jupiter WASP-121 b
Authors:
H. J. Hoeijmakers,
D. Kitzmann,
B. M. Morris,
B. Prinoth,
N. Borsato,
B. Thorsbro,
L. Pino,
E. K. H. Lee,
C. Akın,
J. V. Seidel,
J. L. Birkby,
R. Allart,
K. Heng
Abstract:
Observations of WASP-121 b have suggested an under-abundance of titanium and titanium-oxide from its terminator region. In this study, we aim to determine whether this depletion is global by investigating the day-side emission spectrum. We analyse 8 epochs of high-resolution spectra obtained with ESPRESSO, targeting orbital phases when the day-side is in view. We use a cross-correlation method to…
▽ More
Observations of WASP-121 b have suggested an under-abundance of titanium and titanium-oxide from its terminator region. In this study, we aim to determine whether this depletion is global by investigating the day-side emission spectrum. We analyse 8 epochs of high-resolution spectra obtained with ESPRESSO, targeting orbital phases when the day-side is in view. We use a cross-correlation method to search for various atoms, TiO and VO and compare to models. We constrain the velocities and phase-function of the emission signal using a Bayesian framework. We report significant detections of Ca I, V I, Cr I, Mn I, Fe I, Co I and Ni I, but not T i or TiO. Models containing Ti are unable to reproduce the data. The detected signals are consistent with the known orbital and systemic velocities and with peak emission originating from the sub-stellar point. We find that Ti is depleted from regions of the atmosphere where transmission and emission spectroscopy are sensitive. We interpret this as evidence for the night-side condensation of titanium, preventing it from being mixed back into the upper layers of the atmosphere elsewhere on the planet. Species with lower condensation temperatures are unaffected, implying sharp chemical transitions exist between ultra-hot Jupiters that have slight differences in temperature or dynamical properties. As TiO can act as a strong source of stratospheric heating, cold-trapping creates a coupling between the thermal structures on the day-side and night-side, and thus condensation chemistry needs to be included in global circulation models. Observed elemental abundances in hot Jupiters are not reliably representative of bulk abundances unless night-side condensation is accounted for or the planet is hot enough to avoid night-side cold-traps entirely. Planetary rotation may significantly lower the apparent orbital velocity of emission signals.
△ Less
Submitted 16 January, 2024; v1 submitted 23 October, 2022;
originally announced October 2022.
-
The stable climate of KELT-9b
Authors:
K. D. Jones,
B. M. Morris,
B. -O. Demory,
K. Heng,
M. J. Hooton,
N. Billot,
D. Ehrenreich,
S. Hoyer,
A. E. Simon,
M. Lendl,
O. D. S. Demangeon,
S. G. Sousa,
A. Bonfanti,
T. G. Wilson,
S. Salmon,
Sz. Csizmadia,
H. Parviainen,
G. Bruno,
Y. Alibert,
R. Alonso,
G. Anglada,
T. Bárczy,
D. Barrado y Navascues,
S. C. C. Barros,
W. Baumjohann
, et al. (54 additional authors not shown)
Abstract:
Even among the most irradiated gas giants, so-called ultra-hot Jupiters, KELT-9b stands out as the hottest planet thus far discovered with a dayside temperature of over 4500K. At these extreme irradiation levels, we expect an increase in heat redistribution efficiency and a low Bond albedo owed to an extended atmosphere with molecular hydrogen dissociation occurring on the planetary dayside. We pr…
▽ More
Even among the most irradiated gas giants, so-called ultra-hot Jupiters, KELT-9b stands out as the hottest planet thus far discovered with a dayside temperature of over 4500K. At these extreme irradiation levels, we expect an increase in heat redistribution efficiency and a low Bond albedo owed to an extended atmosphere with molecular hydrogen dissociation occurring on the planetary dayside. We present new photometric observations of the KELT-9 system throughout 4 full orbits and 9 separate occultations obtained by the 30cm space telescope CHEOPS. The CHEOPS bandpass, located at optical wavelengths, captures the peak of the thermal emission spectrum of KELT-9b. In this work we simultaneously analyse CHEOPS phase curves along with public phase curves from TESS and Spitzer to infer joint constraints on the phase curve variation, gravity-darkened transits, and occultation depth in three bandpasses, as well as derive 2D temperature maps of the atmosphere at three different depths. We find a day-night heat redistribution efficiency of $\sim$0.3 which confirms expectations of enhanced energy transfer to the planetary nightside due to dissociation and recombination of molecular hydrogen. We also calculate a Bond albedo consistent with zero. We find no evidence of variability of the brightness temperature of the planet, excluding variability greater than 1% (1$σ$).
△ Less
Submitted 9 August, 2022;
originally announced August 2022.
-
The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
Authors:
The Astropy Collaboration,
Adrian M. Price-Whelan,
Pey Lian Lim,
Nicholas Earl,
Nathaniel Starkman,
Larry Bradley,
David L. Shupe,
Aarya A. Patil,
Lia Corrales,
C. E. Brasseur,
Maximilian Nöthe,
Axel Donath,
Erik Tollerud,
Brett M. Morris,
Adam Ginsburg,
Eero Vaher,
Benjamin A. Weaver,
James Tocknell,
William Jamieson,
Marten H. van Kerkwijk,
Thomas P. Robitaille,
Bruce Merry,
Matteo Bachetti,
H. Moritz Günther,
Thomas L. Aldcroft
, et al. (111 additional authors not shown)
Abstract:
The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as…
▽ More
The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as of the recent major release, version 5.0, and provide major updates for the Project. We then discuss supporting a broader ecosystem of interoperable packages, including connections with several astronomical observatories and missions. We also revisit the future outlook of the Astropy Project and the current status of Learn Astropy. We conclude by raising and discussing the current and future challenges facing the Project.
△ Less
Submitted 28 June, 2022;
originally announced June 2022.
-
The Properties of Fast Yellow Pulsating Supergiants: FYPS Point the Way to Missing Red Supergiants
Authors:
Trevor Z. Dorn-Wallenstein,
Emily M. Levesque,
James R. A. Davenport,
Kathryn F. Neugent,
Brett M. Morris,
K. Azalee Bostroem
Abstract:
Fast yellow pulsating supergiants (FYPS) are a recently-discovered class of evolved massive pulsator. As candidate post-red supergiant objects, and one of the few classes of pulsating evolved massive stars, these objects have incredible potential to change our understanding of the structure and evolution of massive stars. Here we examine the lightcurves of a sample of 126 cool supergiants in the M…
▽ More
Fast yellow pulsating supergiants (FYPS) are a recently-discovered class of evolved massive pulsator. As candidate post-red supergiant objects, and one of the few classes of pulsating evolved massive stars, these objects have incredible potential to change our understanding of the structure and evolution of massive stars. Here we examine the lightcurves of a sample of 126 cool supergiants in the Magellanic Clouds observed by the Transiting Exoplanet Survey Satellite (\tess~) in order to identify pulsating stars. After making quality cuts and filtering out contaminant objects, we examine the distribution of pulsating stars in the Hertzprung-Russel (HR) diagram, and find that FYPS occupy a region above $\log L/L_\odot \gtrsim 5.0$. This luminosity boundary corresponds to stars with initial masses of $\sim$18-20 $M_\odot$, consistent with the most massive red supergiant progenitors of supernovae (SNe) II-P, as well as the observed properties of SNe IIb progenitors. This threshold is in agreement with the picture that FYPS are post-RSG stars. Finally, we characterize the behavior of FYPS pulsations as a function of their location in the HR diagram. We find low frequency pulsations at higher effective temperatures, higher frequency pulsations at lower temperatures, with a transition between the two behaviors at intermediate temperatures. The observed properties of FYPS make them fascinating objects for future theoretical study.
△ Less
Submitted 16 November, 2022; v1 submitted 23 June, 2022;
originally announced June 2022.
-
Modeling Stellar Surface Features on Subgiant Star with an M-dwarf Companion
Authors:
Maria C. Schutte,
Leslie Hebb,
Simon Lowry,
John Wisniewski,
Suzanne L. Hawley,
Suvrath Mahadevan,
Brett M. Morris,
Paul Robertson,
Graeme Rohn,
Gudmundur Stefansson
Abstract:
Understanding magnetic activity on the surface of stars other than the Sun is important for exoplanet analyses to properly characterize an exoplanet's atmosphere and to further characterize stellar activity on a wide range of stars. Modeling stellar surface features of a variety of spectral types and rotation rates are key to understanding of the magnetic activity of these stars. Using data from K…
▽ More
Understanding magnetic activity on the surface of stars other than the Sun is important for exoplanet analyses to properly characterize an exoplanet's atmosphere and to further characterize stellar activity on a wide range of stars. Modeling stellar surface features of a variety of spectral types and rotation rates are key to understanding of the magnetic activity of these stars. Using data from Kepler, we use the starspot modeling program STarSPot (STSP) to measure the position and size of spots for KOI-340 which is an eclipsing binary consisting of a subgiant star ($\rm T_{\rm eff} = 5593 \pm 27 K; R_{*} = 1.98 \pm 0.05 R_{\odot}$) with an M-dwarf companion ($M_{*} = 0.214 \pm 0.006 M_{\odot}$). STSP uses a novel technique to measure the spot positions and radii by using the transiting secondary to study and model individual active regions on the stellar surface using high-precision photometry. We find the average size of spot features on KOI-340's primary is $\sim$10% the radius of the star, i.e. two times larger than the mean size of Solar-maximum sunspots. The spots on KOI-340 are present at every longitude and show possible signs of differential rotation. The minimum fractional spotted area of KOI-340's primary is $2\substack{+12\\ -2} \%$ while the spotted area of the Sun is at most 0.2%. One transit of KOI-340 shows a signal in the transit consistent with a plage; this plage occurs right before a dark spot indicating the plage and spot might be co-located on the surface of the star.
△ Less
Submitted 17 May, 2022;
originally announced May 2022.
-
Weak evidence for variable occultation depth of 55 Cnc e with TESS
Authors:
E. A. Meier Valdés,
B. M. Morris,
R. D. Wells,
N. Schanche,
B. -O. Demory
Abstract:
55 Cnc e is in a 0.73 day orbit transiting a Sun-like star. It has been observed that the occultation depth of this Super-Earth, with a mass of 8$M_{\bigoplus}$ and radius of 2$R_{\bigoplus}$, changes significantly over time at mid-infrared wavelengths. Observations with Spitzer measured a change in its day-side brightness temperature of 1200 K, possibly driven by volcanic activity, magnetic star-…
▽ More
55 Cnc e is in a 0.73 day orbit transiting a Sun-like star. It has been observed that the occultation depth of this Super-Earth, with a mass of 8$M_{\bigoplus}$ and radius of 2$R_{\bigoplus}$, changes significantly over time at mid-infrared wavelengths. Observations with Spitzer measured a change in its day-side brightness temperature of 1200 K, possibly driven by volcanic activity, magnetic star-planet interaction, or the presence of a circumstellar torus of dust. Previous evidence for the variability in occultation was in the infrared range. Here we aim to explore if the variability exists also in the optical. TESS observed 55 Cnc during sectors 21, 44 and 46. We carefully detrend the data and fit a transit and occultation model for each sector in a Markov Chain Monte Carlo routine. In a later stage we use the Leave-One-Out Cross-Validation statistic to compare with a model of constant occultation for the complete set and a model with no occultation. We report an occultation depth of 8$\pm$2.5 ppm for the complete set of TESS observations. In particular, we measured a depth of 15$\pm$4 ppm for sector 21, while for sector 44 we detect no occultation. In sector 46 we measure a weak occultation of 8$\pm$5 ppm. The occultation depth varies from one sector to the next between 1.6 and 3.4 $σ$ significance. We derive the possible contribution on reflected light and thermal emission, setting an upper limit on the geometric albedo. Based on our model comparison the presence of an occultation is favoured considerably over no occultation, where the model with varying occultation across sectors takes most of the statistical weight. Our analysis confirms a detection of the occultation in TESS. Moreover, our results weakly lean towards a varying occultation depth between each sector, while the transit depth is constant across visits.
△ Less
Submitted 17 May, 2022;
originally announced May 2022.
-
The THOR+HELIOS general circulation model: multi-wavelength radiative transfer with accurate scattering by clouds/hazes
Authors:
Russell Deitrick,
Kevin Heng,
Urs Schroffenegger,
Daniel Kitzman,
Simon L. Grimm,
Matej Malik,
João M. Mendonça,
Brett M. Morris
Abstract:
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 a…
▽ More
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.
△ Less
Submitted 7 March, 2022; v1 submitted 4 March, 2022;
originally announced March 2022.
-
CHEOPS geometric albedo of the hot Jupiter HD 209458b
Authors:
A. Brandeker,
K. Heng,
M. Lendl,
J. A. Patel,
B. M. Morris,
C. Broeg,
P. Guterman,
M. Beck,
P. F. L. Maxted,
O. Demangeon,
L. Delrez,
B. -O. Demory,
D. Kitzmann,
N. C. Santos,
V. Singh,
Y. Alibert,
R. Alonso,
G. Anglada,
T. Bárczy,
D. Barrado y Navascues,
S. C. C. Barros,
W. Baumjohann,
T. Beck,
W. Benz,
N. Billot
, et al. (52 additional authors not shown)
Abstract:
We report the detection of the secondary eclipse of the hot Jupiter HD 209458b in optical/visible light using the CHEOPS space telescope. Our measurement of 20.4 +/- 3.3 ppm translates into a geometric albedo of A_g = 0.096 +/- 0.016. The previously estimated dayside temperature of about 1500 K implies that our geometric albedo measurement consists predominantly of reflected starlight and is large…
▽ More
We report the detection of the secondary eclipse of the hot Jupiter HD 209458b in optical/visible light using the CHEOPS space telescope. Our measurement of 20.4 +/- 3.3 ppm translates into a geometric albedo of A_g = 0.096 +/- 0.016. The previously estimated dayside temperature of about 1500 K implies that our geometric albedo measurement consists predominantly of reflected starlight and is largely uncontaminated by thermal emission. This makes the present result one of the most robust measurements of A_g for any exoplanet. Our calculations of the bandpass-integrated geometric albedo demonstrate that the measured value of A_g is consistent with a cloud-free atmosphere, where starlight is reflected via Rayleigh scattering by hydrogen molecules, and the water and sodium abundances are consistent with stellar metallicity. We predict that the bandpass-integrated TESS geometric albedo is too faint to detect and that a phase curve of HD 209458b observed by CHEOPS would have a distinct shape associated with Rayleigh scattering if the atmosphere is indeed cloud free.
△ Less
Submitted 23 February, 2022;
originally announced February 2022.
-
Transit Timing Variations for AU Microscopii b & c
Authors:
Justin M. Wittrock,
Stefan Dreizler,
Michael A. Reefe,
Brett M. Morris,
Peter P. Plavchan,
Patrick J. Lowrance,
Brice-Olivier Demory,
James G. Ingalls,
Emily A. Gilbert,
Thomas Barclay,
Bryson L. Cale,
Karen A. Collins,
Kevin I. Collins,
Ian J. M. Crossfield,
Diana Dragomir,
Jason D. Eastman,
Mohammed El Mufti,
Dax Feliz,
Jonathan Gagne,
Eric Gaidos,
Peter Gao,
Claire S. Geneser,
Leslie Hebb,
Christopher E. Henze,
Keith D. Horne
, et al. (20 additional authors not shown)
Abstract:
We explore the transit timing variations (TTVs) of the young (22 Myr) nearby AU Mic planetary system. For AU Mic b, we introduce three Spitzer (4.5 $μ$m) transits, five TESS transits, 11 LCO transits, one PEST transit, one Brierfield transit, and two transit timing measurements from Rossiter-McLaughlin observations; for AU Mic c, we introduce three TESS transits. We present two independent TTV ana…
▽ More
We explore the transit timing variations (TTVs) of the young (22 Myr) nearby AU Mic planetary system. For AU Mic b, we introduce three Spitzer (4.5 $μ$m) transits, five TESS transits, 11 LCO transits, one PEST transit, one Brierfield transit, and two transit timing measurements from Rossiter-McLaughlin observations; for AU Mic c, we introduce three TESS transits. We present two independent TTV analyses. First, we use EXOFASTv2 to jointly model the Spitzer and ground-based transits and to obtain the midpoint transit times. We then construct an O--C diagram and model the TTVs with Exo-Striker. Second, we reproduce our results with an independent photodynamical analysis. We recover a TTV mass for AU Mic c of 10.8$^{+2.3}_{-2.2}$ M$_{E}$. We compare the TTV-derived constraints to a recent radial-velocity (RV) mass determination. We also observe excess TTVs that do not appear to be consistent with the dynamical interactions of b and c alone, and do not appear to be due to spots or flares. Thus, we present a hypothetical non-transiting "middle-d" candidate exoplanet that is consistent with the observed TTVs, the candidate RV signal, and would establish the AU Mic system as a compact resonant multi-planet chain in a 4:6:9 period commensurability. These results demonstrate that the AU Mic planetary system is dynamically interacting producing detectable TTVs, and the implied orbital dynamics may inform the formation mechanisms for this young system. We recommend future RV and TTV observations of AU Mic b and c to further constrain the masses and to confirm the existence of possible additional planet(s).
△ Less
Submitted 20 April, 2022; v1 submitted 11 February, 2022;
originally announced February 2022.
-
The Effect of Stellar Contamination on Low-resolution Transmission Spectroscopy: Needs Identified by NASA's Exoplanet Exploration Program Study Analysis Group 21
Authors:
Benjamin V. Rackham,
N�stor Espinoza,
Svetlana V. Berdyugina,
Heidi Korhonen,
Ryan J. MacDonald,
Benjamin T. Montet,
Brett M. Morris,
Mahmoudreza Oshagh,
Alexander I. Shapiro,
Yvonne C. Unruh,
Elisa V. Quintana,
Robert T. Zellem,
D�niel Apai,
Thomas Barclay,
Joanna K. Barstow,
Giovanni Bruno,
Ludmila Carone,
Sarah L. Casewell,
Heather M. Cegla,
Serena Criscuoli,
Catherine Fischer,
Damien Fournier,
Mark S. Giampapa,
Helen Giles,
Aishwarya Iyer
, et al. (36 additional authors not shown)
Abstract:
Study Analysis Group 21 (SAG21) of NASA's Exoplanet Exploration Program Analysis Group (ExoPAG) was organized to study the effect of stellar contamination on space-based transmission spectroscopy, a method for studying exoplanetary atmospheres by measuring the wavelength-dependent radius of a planet as it transits its star. Transmission spectroscopy relies on a precise understanding of the spectru…
▽ More
Study Analysis Group 21 (SAG21) of NASA's Exoplanet Exploration Program Analysis Group (ExoPAG) was organized to study the effect of stellar contamination on space-based transmission spectroscopy, a method for studying exoplanetary atmospheres by measuring the wavelength-dependent radius of a planet as it transits its star. Transmission spectroscopy relies on a precise understanding of the spectrum of the star being occulted. However, stars are not homogeneous, constant light sources but have temporally evolving photospheres and chromospheres with inhomogeneities like spots, faculae, plages, granules, and flares. This SAG brought together an interdisciplinary team of more than 100 scientists, with observers and theorists from the heliophysics, stellar astrophysics, planetary science, and exoplanetary atmosphere research communities, to study the current research needs that can be addressed in this context to make the most of transit studies from current NASA facilities like HST and JWST. The analysis produced 14 findings, which fall into three Science Themes encompassing (1) how the Sun is used as our best laboratory to calibrate our understanding of stellar heterogeneities ("The Sun as the Stellar Benchmark"), (2) how stars other than the Sun extend our knowledge of heterogeneities ("Surface Heterogeneities of Other Stars") and (3) how to incorporate information gathered for the Sun and other stars into transit studies ("Mapping Stellar Knowledge to Transit Studies"). In this invited review, we largely reproduce the final report of SAG21 as a contribution to the peer-reviewed literature.
△ Less
Submitted 17 March, 2023; v1 submitted 24 January, 2022;
originally announced January 2022.
-
Titanium oxide and chemical inhomogeneity in the atmosphere of the exoplanet WASP-189b
Authors:
Bibiana Prinoth,
H. Jens Hoeijmakers,
Daniel Kitzmann,
Elin Sandvik,
Julia V. Seidel,
Monika Lendl,
Nicholas W. Borsato,
Brian Thorsbro,
David R. Anderson,
David Barrado,
Kateryna Kravchenko,
Romain Allart,
Vincent Bourrier,
Heather M. Cegla,
David Ehrenreich,
Chloe Fisher,
Christophe Lovis,
Andrea Guzm�n-Mesa,
Simon Grimm,
Matthew Hooton,
Brett M. Morris,
Maria Oreshenko,
Lorenzo Pino,
Kevin Heng
Abstract:
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 spect…
▽ More
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.
△ Less
Submitted 30 January, 2022; v1 submitted 24 November, 2021;
originally announced November 2021.
-
Physically-motivated basis functions for temperature maps of exoplanets
Authors:
Brett M. Morris,
Kevin Heng,
Kathryn Jones,
Caroline Piaulet,
Brice-Olivier Demory,
Daniel Kitzmann,
H. Jens Hoeijmakers
Abstract:
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 spheri…
▽ More
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.
△ Less
Submitted 22 October, 2021;
originally announced October 2021.
-
Spi-OPS: Spitzer and CHEOPS confirm the near-polar orbit of MASCARA-1 b and reveal a hint of dayside reflection
Authors:
M. J. Hooton,
S. Hoyer,
D. Kitzmann,
B. M. Morris,
A. M. S. Smith,
A. Collier Cameron,
D. Futyan,
P. F. L. Maxted,
D. Queloz,
B. -O. Demory,
K. Heng,
M. Lendl,
J. Cabrera,
Sz. Csizmadia,
A. Deline,
H. Parviainen,
S. Salmon,
S. Sulis,
T. G. Wilson,
A. Bonfanti,
A. Brandeker,
O. D. S. Demangeon,
M. Oshagh,
C. M. Persson,
G. Scandariato
, et al. (60 additional authors not shown)
Abstract:
The light curves of tidally locked hot Jupiters transiting fast-rotating, early-type stars are a rich source of information about both the planet and star, with full-phase coverage enabling a detailed atmospheric characterisation of the planet. Although it is possible to determine the true spin-orbit angle $Ψ$, a notoriously difficult parameter to measure, from any transit asymmetry resulting from…
▽ More
The light curves of tidally locked hot Jupiters transiting fast-rotating, early-type stars are a rich source of information about both the planet and star, with full-phase coverage enabling a detailed atmospheric characterisation of the planet. Although it is possible to determine the true spin-orbit angle $Ψ$, a notoriously difficult parameter to measure, from any transit asymmetry resulting from gravity darkening induced by the stellar rotation, the correlations that exist between the transit parameters have led to large disagreements in published values of $Ψ$ for some systems. We aimed to study these phenomena in the light curves of the ultra-hot Jupiter MASCARA-1 b. We obtained optical CHEOPS transit and occultation light curves of MASCARA-1 b, and analysed them jointly with a Spitzer/IRAC 4.5 $μ$m full-phase curve. When fitting the CHEOPS and Spitzer transits together, the degeneracies are greatly diminished and return results consistent with previously published Doppler tomography. Placing priors informed by the tomography achieves even better precision, allowing a determination of $Ψ=72.1^{+2.5}_{-2.4}$ deg. From the occultations and phase variations, we derived dayside and nightside temperatures of $3062^{+66}_{-68}$ K and $1720\pm330$ K, respectively. In addition, we could separately derive geometric albedo $A_g=0.171^{+0.066}_{-0.068}$ and spherical albedo $A_s=0.266^{+0.097}_{-0.100}$ from the CHEOPS data, and Bond albedo $A_B=0.057^{+0.083}_{-0.101}$ from the Spitzer phase curve. Where possible, priors informed by Doppler tomography should be used when fitting transits of fast-rotating stars, though multi-colour photometry may also unlock an accurate measurement of $Ψ$. Our approach to modelling the phase variations at different wavelengths provides a template for how to separate thermal emission from reflected light in spectrally resolved JWST phase curves.
△ Less
Submitted 4 February, 2022; v1 submitted 10 September, 2021;
originally announced September 2021.
-
The Pandora SmallSat: Multiwavelength Characterization of Exoplanets and their Host Stars
Authors:
Elisa V. Quintana,
Knicole D. Colón,
Gregory Mosby,
Joshua E. Schlieder,
Pete Supsinskas,
Jordan Karburn,
Jessie L. Dotson,
Thomas P. Greene,
Christina Hedges,
Dániel Apai,
Thomas Barclay,
Jessie L. Christiansen,
Néstor Espinoza,
Susan E. Mullally,
Emily A. Gilbert,
Kelsey Hoffman,
Veselin B. Kostov,
Nikole K. Lewis,
Trevor O. Foote,
James Mason,
Allison Youngblood,
Brett M. Morris,
Elisabeth R. Newton,
Joshua Pepper,
Benjamin V. Rackham
, et al. (2 additional authors not shown)
Abstract:
Pandora is a SmallSat mission designed to study the atmospheres of exoplanets, and was selected as part of NASA's Astrophysics Pioneers Program. Transmission spectroscopy of transiting exoplanets provides our best opportunity to identify the makeup of planetary atmospheres in the coming decade. Stellar brightness variations due to star spots, however, can impact these measurements and contaminate…
▽ More
Pandora is a SmallSat mission designed to study the atmospheres of exoplanets, and was selected as part of NASA's Astrophysics Pioneers Program. Transmission spectroscopy of transiting exoplanets provides our best opportunity to identify the makeup of planetary atmospheres in the coming decade. Stellar brightness variations due to star spots, however, can impact these measurements and contaminate the observed spectra. Pandora's goal is to disentangle star and planet signals in transmission spectra to reliably determine exoplanet atmosphere compositions. Pandora will collect long-duration photometric observations with a visible-light channel and simultaneous spectra with a near-IR channel. The broad-wavelength coverage will provide constraints on the spot and faculae covering fractions of low-mass exoplanet host stars and the impact of these active regions on exoplanetary transmission spectra. Pandora will subsequently identify exoplanets with hydrogen- or water-dominated atmospheres, and robustly determine which planets are covered by clouds and hazes. Pandora will observe at least 20 exoplanets with sizes ranging from Earth-size to Jupiter-size and host stars spanning mid-K to late-M spectral types. The project is made possible by leveraging investments in other projects, including an all-aluminum 0.45-meter Cassegrain telescope design, and a NIR sensor chip assembly from the James Webb Space Telescope. The mission will last five years from initial formulation to closeout, with one-year of science operations. Launch is planned for the mid-2020s as a secondary payload in Sun-synchronous low-Earth orbit. By design, Pandora has a diverse team, with over half of the mission leadership roles filled by early career scientists and engineers, demonstrating the high value of SmallSats for developing the next generation of space mission leaders.
△ Less
Submitted 19 August, 2021; v1 submitted 13 August, 2021;
originally announced August 2021.
-
CHEOPS Precision Phase Curve of the Super-Earth 55 Cnc e
Authors:
B. M. Morris,
L. Delrez,
A. Brandeker,
A. C. Cameron,
A. E. Simon,
D. Futyan,
G. Olofsson,
S. Hoyer,
A. Fortier,
B. -O. Demory,
M. Lendl,
T. G. Wilson,
M. Oshagh,
K. Heng,
D. Ehrenreich,
S. Sulis,
Y. Alibert,
R. Alonso,
G. Anglada Escudé,
D. Barrado,
S. C. C. Barros,
W. Baumjohann,
M. Beck,
T. Beck,
A. Bekkelien
, et al. (57 additional authors not shown)
Abstract:
55 Cnc e is a transiting super-Earth (radius $1.88\rm\,R_\oplus$ and mass $8\rm\, M_\oplus$) orbiting a G8V host star on a 17-hour orbit. Spitzer observations of the planet's phase curve at 4.5 $μ$m revealed a time-varying occultation depth, and MOST optical observations are consistent with a time-varying phase curve amplitude and phase offset of maximum light. Both broadband and high-resolution s…
▽ More
55 Cnc e is a transiting super-Earth (radius $1.88\rm\,R_\oplus$ and mass $8\rm\, M_\oplus$) orbiting a G8V host star on a 17-hour orbit. Spitzer observations of the planet's phase curve at 4.5 $μ$m revealed a time-varying occultation depth, and MOST optical observations are consistent with a time-varying phase curve amplitude and phase offset of maximum light. Both broadband and high-resolution spectroscopic analyses are consistent with either a high mean molecular weight atmosphere or no atmosphere for planet e. A long term photometric monitoring campaign on an independent optical telescope is needed to probe the variability in this system. We seek to measure the phase variations of 55 Cnc e with a broadband optical filter with the 30 cm effective aperture space telescope CHEOPS and explore how the precision photometry narrows down the range of possible scenarios. We observed 55 Cnc for 1.6 orbital phases in March of 2020. We designed a phase curve detrending toolkit for CHEOPS photometry which allows us to study the underlying flux variations of the 55 Cnc system. We detected a phase variation with a full-amplitude of $72 \pm 7$ ppm but do not detect a significant secondary eclipse of the planet. The shape of the phase variation resembles that of a piecewise-Lambertian, however the non-detection of the planetary secondary eclipse, and the large amplitude of the variations exclude reflection from the planetary surface as a possible origin of the observed phase variations. They are also likely incompatible with magnetospheric interactions between the star and planet but may imply that circumplanetary or circumstellar material modulate the flux of the system. Further precision photometry of 55 Cnc from CHEOPS will measure variations in the phase curve amplitude and shape over time this year.
△ Less
Submitted 14 June, 2021;
originally announced June 2021.
-
A CHEOPS White Dwarf Transit Search
Authors:
Brett M. Morris,
Kevin Heng,
Alexis Brandeker,
Andrew Swan,
Monika Lendl
Abstract:
White dwarf spectroscopy shows that nearly half of white dwarf atmospheres contain metals that must have been accreted from planetary material that survived the red giant phases of stellar evolution. We can use metal pollution in white dwarf atmospheres as flags, signalling recent accretion, in order to prioritize an efficient sample of white dwarfs to search for transiting material. We present a…
▽ More
White dwarf spectroscopy shows that nearly half of white dwarf atmospheres contain metals that must have been accreted from planetary material that survived the red giant phases of stellar evolution. We can use metal pollution in white dwarf atmospheres as flags, signalling recent accretion, in order to prioritize an efficient sample of white dwarfs to search for transiting material. We present a search for planetesimals orbiting six nearby white dwarfs with the CHEOPS spacecraft. The targets are relatively faint for CHEOPS, $11$ mag $< G < 12.8$ mag. We use aperture photometry data products from the CHEOPS mission as well as custom PSF photometry to search for periodic variations in flux due to transiting planetesimals. We detect no significant variations in flux that cannot be attributed to spacecraft systematics, despite reaching a photometric precision of $<2$ ppt in 60 s exposures on each target. We simulate observations to show that the small survey is sensitive primarily to Moon-sized transiting objects with periods $3$ hr $< P < 10$ hr, with radii $R \gtrsim 1000$ km.
△ Less
Submitted 17 May, 2021;
originally announced May 2021.
-
exoplanet: Gradient-based probabilistic inference for exoplanet data & other astronomical time series
Authors:
Daniel Foreman-Mackey,
Rodrigo Luger,
Eric Agol,
Thomas Barclay,
Luke G. Bouma,
Timothy D. Brandt,
Ian Czekala,
Trevor J. David,
Jiayin Dong,
Emily A. Gilbert,
Tyler A. Gordon,
Christina Hedges,
Daniel R. Hey,
Brett M. Morris,
Adrian M. Price-Whelan,
Arjun B. Savel
Abstract:
"exoplanet" is a toolkit for probabilistic modeling of astronomical time series data, with a focus on observations of exoplanets, using PyMC3 (Salvatier et al., 2016). PyMC3 is a flexible and high-performance model-building language and inference engine that scales well to problems with a large number of parameters. "exoplanet" extends PyMC3's modeling language to support many of the custom functi…
▽ More
"exoplanet" is a toolkit for probabilistic modeling of astronomical time series data, with a focus on observations of exoplanets, using PyMC3 (Salvatier et al., 2016). PyMC3 is a flexible and high-performance model-building language and inference engine that scales well to problems with a large number of parameters. "exoplanet" extends PyMC3's modeling language to support many of the custom functions and probability distributions required when fitting exoplanet datasets or other astronomical time series. While it has been used for other applications, such as the study of stellar variability, the primary purpose of "exoplanet" is the characterization of exoplanets or multiple star systems using time-series photometry, astrometry, and/or radial velocity. In particular, the typical use case would be to use one or more of these datasets to place constraints on the physical and orbital parameters of the system, such as planet mass or orbital period, while simultaneously taking into account the effects of stellar variability.
△ Less
Submitted 24 June, 2021; v1 submitted 5 May, 2021;
originally announced May 2021.
-
Near real-time precipitable water vapour monitoring for correcting near-infrared observations using satellite remote sensing
Authors:
E. A. Meier Vald�s,
B. M. Morris,
B. -O. Demory
Abstract:
In the search for small exoplanets orbiting cool stars whose spectral energy distributions peak in the near infrared, the strong absorption of radiation in this region due to water vapour in the atmosphere is a particularly adverse effect for the ground-based observations of cool stars. To achieve the photometric precision required to detect exoplanets in the near infrared, it is necessary to miti…
▽ More
In the search for small exoplanets orbiting cool stars whose spectral energy distributions peak in the near infrared, the strong absorption of radiation in this region due to water vapour in the atmosphere is a particularly adverse effect for the ground-based observations of cool stars. To achieve the photometric precision required to detect exoplanets in the near infrared, it is necessary to mitigate the impact of variable precipitable water vapour (PWV) on radial-velocity and photometric measurements. The aim is to enable global PWV correction by monitoring the amount of precipitable water vapour at zenith and along the line of sight of any visible target. We developed an open source Python package that uses Geostationary Operational Environmental Satellites (GOES) imagery data, which provides temperature and relative humidity at different pressure levels to compute near real-time PWV above any ground-based observatory covered by GOES every 5 minutes or 10 minutes depending on the location. We computed PWV values on selected days above Cerro Paranal (Chile) and San Pedro M�rtir (Mexico) to benchmark the procedure. We also simulated different pointing at test targets as observed from the sites to compute the PWV along the line of sight. To asses the accuracy of our method, we compared our results with the on-site radiometer measurements obtained from Cerro Paranal. Our results show that our publicly-available code proves to be a good supporting tool for measuring the local PWV for any ground-based facility within the GOES coverage, which will help in reducing correlated noise contributions in near-infrared ground-based observations that do not benefit from on-site PWV measurements.
△ Less
Submitted 9 March, 2021;
originally announced March 2021.
-
Closed-formed ab initio solutions of geometric albedos and reflected light phase curves of exoplanets
Authors:
Kevin Heng,
Brett M. Morris,
Daniel Kitzmann
Abstract:
Studying the albedos of the planets and moons of the Solar System dates back at least a century. Of particular interest is the relationship between the albedo measured at superior conjunction, known as the ``geometric albedo", and the albedo considered over all orbital phase angles, known as the ``spherical albedo". Determining the relationship between the geometric and spherical albedos usually i…
▽ More
Studying the albedos of the planets and moons of the Solar System dates back at least a century. Of particular interest is the relationship between the albedo measured at superior conjunction, known as the ``geometric albedo", and the albedo considered over all orbital phase angles, known as the ``spherical albedo". Determining the relationship between the geometric and spherical albedos usually involves complex numerical calculations and closed-form solutions are restricted to simple reflection laws. Here we report the discovery of closed-form solutions for the geometric albedo and integral phase function, which apply to any law of reflection that only depends on the scattering angle. The shape of a reflected light phase curve, quantified by the integral phase function, and the secondary eclipse depth, quantified by the geometric albedo, may now be self-consistently inverted to retrieve globally averaged physical parameters. Fully Bayesian phase curve inversions for reflectance maps and simultaneous light curve detrending may now be performed due to the efficiency of computation. Demonstrating these innovations for the hot Jupiter Kepler-7b, we infer a geometric albedo of $0.25^{+0.01}_{-0.02}$, a phase integral of $1.77 \pm0.07$, a spherical albedo of $0.44^{+0.02}_{-0.03}$ and a scattering asymmetry factor of $0.07^{+0.12}_{-0.11}$.
△ Less
Submitted 26 August, 2021; v1 submitted 3 March, 2021;
originally announced March 2021.
-
Refining the transit timing and photometric analysis of TRAPPIST-1: Masses, radii, densities, dynamics, and ephemerides
Authors:
Eric Agol,
Caroline Dorn,
Simon L. Grimm,
Martin Turbet,
Elsa Ducrot,
Laetitia Delrez,
Michael Gillon,
Brice-Olivier Demory,
Artem Burdanov,
Khalid Barkaoui,
Zouhair Benkhaldoun,
Emeline Bolmont,
Adam Burgasser,
Sean Carey,
Julien de Wit,
Daniel Fabrycky,
Daniel Foreman-Mackey,
Jonas Haldemann,
David M. Hernandez,
James Ingalls,
Emmanuel Jehin,
Zachary Langford,
Jeremy Leconte,
Susan M. Lederer,
Rodrigo Luger
, et al. (10 additional authors not shown)
Abstract:
We have collected transit times for the TRAPPIST-1 system with the Spitzer Space Telescope over four years. We add to these ground-based, HST and K2 transit time measurements, and revisit an N-body dynamical analysis of the seven-planet system using our complete set of times from which we refine the mass ratios of the planets to the star. We next carry out a photodynamical analysis of the Spitzer…
▽ More
We have collected transit times for the TRAPPIST-1 system with the Spitzer Space Telescope over four years. We add to these ground-based, HST and K2 transit time measurements, and revisit an N-body dynamical analysis of the seven-planet system using our complete set of times from which we refine the mass ratios of the planets to the star. We next carry out a photodynamical analysis of the Spitzer light curves to derive the density of the host star and the planet densities. We find that all seven planets' densities may be described with a single rocky mass-radius relation which is depleted in iron relative to Earth, with Fe 21 wt% versus 32 wt% for Earth, and otherwise Earth-like in composition. Alternatively, the planets may have an Earth-like composition, but enhanced in light elements, such as a surface water layer or a core-free structure with oxidized iron in the mantle. We measure planet masses to a precision of 3-5%, equivalent to a radial-velocity (RV) precision of 2.5 cm/sec, or two orders of magnitude more precise than current RV capabilities. We find the eccentricities of the planets are very small; the orbits are extremely coplanar; and the system is stable on 10 Myr timescales. We find evidence of infrequent timing outliers which we cannot explain with an eighth planet; we instead account for the outliers using a robust likelihood function. We forecast JWST timing observations, and speculate on possible implications of the planet densities for the formation, migration and evolution of the planet system.
△ Less
Submitted 14 January, 2021; v1 submitted 2 October, 2020;
originally announced October 2020.
-
The hot dayside and asymmetric transit of WASP-189b seen by CHEOPS
Authors:
M. Lendl,
Sz. Csizmadia,
A. Deline,
L. Fossati,
D. Kitzmann,
K. Heng,
S. Hoyer,
S. Salmon,
W. Benz,
C. Broeg,
D. Ehrenreich,
A. Fortier,
D. Queloz,
A. Bonfanti,
A. Brandeker,
A. Collier Cameron,
L. Delrez,
A. Garcia Muñoz,
M. J. Hooton,
P. F. L. Maxted,
B. M. Morris,
V. Van Grootel,
T. G. Wilson,
Y. Alibert,
R. Alonso
, et al. (80 additional authors not shown)
Abstract:
The CHEOPS space mission dedicated to exoplanet follow-up was launched in December 2019, equipped with the capacity to perform photometric measurements at the 20 ppm level. As CHEOPS carries out its observations in a broad optical passband, it can provide insights into the reflected light from exoplanets and constrain the short-wavelength thermal emission for the hottest of planets by observing oc…
▽ More
The CHEOPS space mission dedicated to exoplanet follow-up was launched in December 2019, equipped with the capacity to perform photometric measurements at the 20 ppm level. As CHEOPS carries out its observations in a broad optical passband, it can provide insights into the reflected light from exoplanets and constrain the short-wavelength thermal emission for the hottest of planets by observing occultations and phase curves. Here, we report the first CHEOPS observation of an occultation, namely, that of the hot Jupiter WASP-189b, a $M_P \approx 2 M_J$ planet orbiting an A-type star. We detected the occultation of WASP-189 b at high significance in individual measurements and derived an occultation depth of $dF = 87.9 \pm 4.3$ppm based on four occultations. We compared these measurements to model predictions and we find that they are consistent with an unreflective atmosphere heated to a temperature of $3435 \pm 27$K, when assuming inefficient heat redistribution. Furthermore, we present two transits of WASP-189b observed by CHEOPS. These transits have an asymmetric shape that we attribute to gravity darkening of the host star caused by its high rotation rate. We used these measurements to refine the planetary parameters, finding a $\sim25\%$ deeper transit compared to the discovery paper and updating the radius of WASP-189b to $1.619\pm0.021 R_J$. We further measured the projected orbital obliquity to be $λ= 86.4^{+2.9}_{-4.4}$deg, a value that is in good agreement with a previous measurement from spectroscopic observations, and derived a true obliquity of $Ψ= 85.4\pm4.3$deg. Finally, we provide reference values for the photometric precision attained by the CHEOPS satellite: for the V=6.6 mag star, and using a one-hour binning, we obtain a residual RMS between 10 and 17ppm on the individual light curves, and 5.7ppm when combining the four visits.
△ Less
Submitted 28 September, 2020;
originally announced September 2020.
-
A super-Earth and a sub-Neptune orbiting the bright, quiet M3 dwarf TOI-1266
Authors:
B. -O. Demory,
F. J. Pozuelos,
Y. Gomez Maqueo Chew,
L. Sabin,
R. Petrucci,
U. Schroffenegger,
S. L. Grimm,
M. Sestovic,
M. Gillon,
J. McCormac,
K. Barkaoui,
W. Benz,
A. Bieryla,
F. Bouchy,
A. Burdanov,
K. A. Collins,
J. de Wit,
C. D. Dressing,
L. J. Garcia,
S. Giacalone,
P. Guerra,
J. Haldemann,
K. Heng,
E. Jehin,
E. Jofre
, et al. (48 additional authors not shown)
Abstract:
We report the discovery and characterisation of a super-Earth and a sub-Neptune transiting the bright ($K=8.8$), quiet, and nearby (37 pc) M3V dwarf TOI-1266. We validate the planetary nature of TOI-1266 b and c using four sectors of TESS photometry and data from the newly-commissioned 1-m SAINT-EX telescope located in San Pedro Mártir (Mexico). We also include additional ground-based follow-up ph…
▽ More
We report the discovery and characterisation of a super-Earth and a sub-Neptune transiting the bright ($K=8.8$), quiet, and nearby (37 pc) M3V dwarf TOI-1266. We validate the planetary nature of TOI-1266 b and c using four sectors of TESS photometry and data from the newly-commissioned 1-m SAINT-EX telescope located in San Pedro Mártir (Mexico). We also include additional ground-based follow-up photometry as well as high-resolution spectroscopy and high-angular imaging observations. The inner, larger planet has a radius of $R=2.37_{-0.12}^{+0.16}$ R$_{\oplus}$ and an orbital period of 10.9 days. The outer, smaller planet has a radius of $R=1.56_{-0.13}^{+0.15}$ R$_{\oplus}$ on an 18.8-day orbit. The data are found to be consistent with circular, co-planar and stable orbits that are weakly influenced by the 2:1 mean motion resonance. Our TTV analysis of the combined dataset enables model-independent constraints on the masses and eccentricities of the planets. We find planetary masses of $M_\mathrm{p}$ = $13.5_{-9.0}^{+11.0}$ $\mathrm{M_{\oplus}}$ ($<36.8$ $\mathrm{M_{\oplus}}$ at 2-$σ$) for TOI-1266 b and $2.2_{-1.5}^{+2.0}$ $\mathrm{M_{\oplus}}$ ($<5.7$ $\mathrm{M_{\oplus}}$ at 2-$σ$) for TOI-1266 c. We find small but non-zero orbital eccentricities of $0.09_{-0.05}^{+0.06}$ ($<0.21$ at 2-$σ$) for TOI-1266 b and $0.04\pm0.03$ ($<0.10$ at 2-$σ$) for TOI-1266 c. The equilibrium temperatures of both planets are of $413\pm20$ K and $344\pm16$ K, respectively, assuming a null Bond albedo and uniform heat redistribution from the day-side to the night-side hemisphere. The host brightness and negligible activity combined with the planetary system architecture and favourable planet-to-star radii ratios makes TOI-1266 an exquisite system for a detailed characterisation.
△ Less
Submitted 9 September, 2020;
originally announced September 2020.
-
Short Term Variability of Evolved Massive Stars with TESS II: A New Class of Cool, Pulsating Supergiants
Authors:
Trevor Z. Dorn-Wallenstein,
Emily M. Levesque,
Kathryn F. Neugent,
James R. A. Davenport,
Brett M. Morris,
Keyan Gootkin
Abstract:
Massive stars briefly pass through the yellow supergiant (YSG) phase as they evolve redward across the HR diagram and expand into red supergiants (RSGs). Higher-mass stars pass through the YSG phase again as they evolve blueward after experiencing significant RSG mass loss. These post-RSG objects offer us a tantalizing glimpse into which stars end their lives as RSGs, and why. One telltale sign of…
▽ More
Massive stars briefly pass through the yellow supergiant (YSG) phase as they evolve redward across the HR diagram and expand into red supergiants (RSGs). Higher-mass stars pass through the YSG phase again as they evolve blueward after experiencing significant RSG mass loss. These post-RSG objects offer us a tantalizing glimpse into which stars end their lives as RSGs, and why. One telltale sign of a post-RSG object may be an instability to pulsations, depending on the star's interior structure. Here we report the discovery of five YSGs with pulsation periods faster than 1 day, found in a sample of 76 cool supergiants observed by \tess at two-minute cadence. These pulsating YSGs are concentrated in a HR diagram region not previously associated with pulsations; we conclude that this is a genuine new class of pulsating star, Fast Yellow Pulsating Supergiants (FYPS). For each FYPS, we extract frequencies via iterative prewhitening and conduct a time-frequency analysis. One FYPS has an extracted frequency that is split into a triplet, and the amplitude of that peak is modulated on the same timescale as the frequency spacing of the triplet; neither rotation nor binary effects are likely culprits. We discuss the evolutionary status of FYPS and conclude that they are candidate post-RSGs. All stars in our sample also show the same stochastic low-frequency variability (SLFV) found in hot OB stars and attributed to internal gravity waves. Finally, we find four $α$ Cygni variables in our sample, of which three are newly discovered.
△ Less
Submitted 26 August, 2020;
originally announced August 2020.
-
Hunt for Starspots in HARPS Spectra of G and K Stars
Authors:
Brett M. Morris,
H. Jens Hoeijmakers,
Daniel Kitzmann,
Brice-Olivier Demory
Abstract:
We present a method for detecting starspots on cool stars using the cross-correlation function (CCF) of high resolution molecular spectral templates applied to archival high-resolution spectra of G and K stars observed with HARPS/HARPS-N. We report non-detections of starspots on the Sun even when the Sun was spotted, the solar twin 18 Scorpii, and the very spotted Sun-like star HAT-P-11, suggestin…
▽ More
We present a method for detecting starspots on cool stars using the cross-correlation function (CCF) of high resolution molecular spectral templates applied to archival high-resolution spectra of G and K stars observed with HARPS/HARPS-N. We report non-detections of starspots on the Sun even when the Sun was spotted, the solar twin 18 Scorpii, and the very spotted Sun-like star HAT-P-11, suggesting that Sun-like starspot distributions will be invisible to the CCF technique, and should not produce molecular absorption signals which might be confused for signatures of exoplanet atmospheres. We detect strong TiO absorption in the T Tauri K-dwarfs LkCa 4 and AA Tau, consistent with significant coverage by cool regions. We show that despite the non-detections, the technique is sensitive to relatively small spot coverages on M dwarfs and large starspot areas on Sun-like stars.
△ Less
Submitted 14 May, 2020;
originally announced May 2020.
-
A Relationship Between Stellar Age and Spot Coverage
Authors:
Brett M. Morris
Abstract:
We investigate starspot distributions consistent with space-based photometry of F, G, and K stars in six stellar associations ranging in age from 10 Myr to 4 Gyr. We show that a simple light curve statistic called the "smoothed amplitude" is proportional to stellar age as $t^{-1/2}$, following a Skumanich-like spin-down relation. We marginalize over the unknown stellar inclinations by forward mode…
▽ More
We investigate starspot distributions consistent with space-based photometry of F, G, and K stars in six stellar associations ranging in age from 10 Myr to 4 Gyr. We show that a simple light curve statistic called the "smoothed amplitude" is proportional to stellar age as $t^{-1/2}$, following a Skumanich-like spin-down relation. We marginalize over the unknown stellar inclinations by forward modeling the ensemble of light curves for direct comparison with the Kepler, K2 and TESS photometry. We sample the posterior distributions for spot coverage with Approximate Bayesian Computation. We find typical spot coverages in the range 1-10% which decrease with increasing stellar age. The spot coverage is proportional to $t^n$ where $n =-0.37 \pm 0.16$, also statistically consistent with a Skumanich-like $t^{-1/2}$ decay of starspot coverage with age. We apply two techniques to estimate the spot coverage of young exoplanet-hosting stars likely to be targeted for transmission spectroscopy with the James Webb Space Telescope, and estimate the bias in exoplanet radius measurements due to varying starspot coverage.
△ Less
Submitted 21 February, 2020;
originally announced February 2020.
-
The Stellar Variability Noise Floor for Transiting Exoplanet Photometry with PLATO
Authors:
Brett M. Morris,
Monica G. Bobra,
Eric Agol,
Yu Jin Lee,
Suzanne L. Hawley
Abstract:
One of the main science motivations for the ESA PLAnetary Transit and Oscillations (PLATO) mission is to measure exoplanet transit radii with 3% precision. In addition to flares and starspots, stellar oscillations and granulation will enforce fundamental noise floors for transiting exoplanet radius measurements. We simulate light curves of Earth-sized exoplanets transiting continuum intensity imag…
▽ More
One of the main science motivations for the ESA PLAnetary Transit and Oscillations (PLATO) mission is to measure exoplanet transit radii with 3% precision. In addition to flares and starspots, stellar oscillations and granulation will enforce fundamental noise floors for transiting exoplanet radius measurements. We simulate light curves of Earth-sized exoplanets transiting continuum intensity images of the Sun taken by the HMI instrument aboard SDO to investigate the uncertainties introduced on the exoplanet radius measurements by stellar granulation and oscillations. After modeling the solar variability with a Gaussian process, we find that the amplitude of solar oscillations and granulation is of order 100 ppm -- similar to the depth of an Earth transit -- and introduces a fractional uncertainty on the depth of transit of 0.73% assuming four transits are observed over the mission duration. However, when we translate the depth measurement into a radius measurement of the planet, we find a much larger radius uncertainty of 3.6%. This is due to a degeneracy between the transit radius ratio, the limb-darkening, and the impact parameter caused by the inability to constrain the transit impact parameter in the presence of stellar variability. We find that surface brightness inhomogeneity due to photospheric granulation contributes a lower limit of only 2 ppm to the photometry in-transit. The radius uncertainty due to granulation and oscillations, combined with the degeneracy with the transit impact parameter, accounts for a significant fraction of the error budget of the PLATO mission, before detector or observational noise is introduced to the light curve. If it is possible to constrain the impact parameter or to obtain follow-up observations at longer wavelengths where limb-darkening is less significant, this may enable higher precision radius measurements.
△ Less
Submitted 19 February, 2020;
originally announced February 2020.
-
Temporal Evolution of Spatially-Resolved Individual Star Spots on a Planet-Hosting Solar-type Star: Kepler 17
Authors:
Kosuke Namekata,
James R. A. Davenport,
Brett M. Morris,
Suzanne L. Hawley,
Hiroyuki Maehara,
Yuta Notsu,
Shin Toriumi,
Kai Ikuta,
Shota Notsu,
Satoshi Honda,
Daisaku Nogami,
Kazunari Shibata
Abstract:
Star spot evolution is visible evidence of the emergence/decay of the magnetic field on stellar surface, and it is therefore important for the understanding of the underlying stellar dynamo and consequential stellar flares. In this paper, we report the temporal evolution of individual star spot area on the hot-Jupiter-hosting active solar-type star Kepler 17 whose transits occur every 1.5 days. Th…
▽ More
Star spot evolution is visible evidence of the emergence/decay of the magnetic field on stellar surface, and it is therefore important for the understanding of the underlying stellar dynamo and consequential stellar flares. In this paper, we report the temporal evolution of individual star spot area on the hot-Jupiter-hosting active solar-type star Kepler 17 whose transits occur every 1.5 days. The spot longitude and area evolution are estimated (1) from the stellar rotational modulations of Kepler data and (2) from the brightness enhancements during the exoplanet transits caused by existence of large star spots. As a result of the comparison, number of spots, spot locations, and the temporal evolution derived from the rotational modulations is largely different from those of in-transit spots. We confirm that although only two light curve minima appear per rotation, there are clearly many spots present on the star. We find that the observed differential intensity changes are sometimes consistent with the spot pattern detected by transits, but they sometimes do not match with each other. Although the temporal evolution derived from the rotational modulation differs from those of in-transit spots to a certain degree, the emergence/decay rates of in-transit spots are within an order of magnitude of those derived for sunspots as well as our previous research based only on rotational modulations. This supports a hypothesis that the emergence/decay of sunspots and extremely-large star spots on solar-type stars occur through the same underlying processes.
△ Less
Submitted 3 February, 2020;
originally announced February 2020.
-
Exploring the atmospheric dynamics of the extreme ultra-hot Jupiter KELT-9b using TESS photometry
Authors:
Ian Wong,
Avi Shporer,
Daniel Kitzmann,
Brett M. Morris,
Kevin Heng,
H. Jens Hoeijmakers,
Brice-Olivier Demory,
John P. Ahlers,
Megan Mansfield,
Jacob L. Bean,
Tansu Daylan,
Tara Fetherolf,
Joseph E. Rodriguez,
Björn Benneke,
George R. Ricker,
David W. Latham,
Roland Vanderspek,
Sara Seager,
Joshua N. Winn,
Jon M. Jenkins,
Christopher J. Burke,
Jessie Christiansen,
Zahra Essack,
Mark E. Rose,
Jeffrey C. Smith
, et al. (2 additional authors not shown)
Abstract:
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 ec…
▽ More
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.
△ Less
Submitted 30 July, 2020; v1 submitted 3 October, 2019;
originally announced October 2019.
-
Kepler Object of Interest Network III. Kepler-82f: A new non-transiting $21 M_\bigoplus$ planet from photodynamical modelling
Authors:
J. Freudenthal,
C. von Essen,
A. Ofir,
S. ~Dreizler,
E. Agol,
S. Wedemeyer,
B. M. Morris,
A. C. Becker,
H. J. Deeg,
S. Hoyer,
M. Mallonn,
K. Poppenhaeger,
E. Herrero,
I. Ribas,
P. Boumis,
A. Liakos
Abstract:
Context. The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised for follow-up observations of transiting planet candidate Kepler objects of interest (KOIs) with large transit timing variations (TTVs). The main goal of KOINet is the completion of their TTV curves as the Kepler telescope stopped observing the original Kepler field in 2013.
A…
▽ More
Context. The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised for follow-up observations of transiting planet candidate Kepler objects of interest (KOIs) with large transit timing variations (TTVs). The main goal of KOINet is the completion of their TTV curves as the Kepler telescope stopped observing the original Kepler field in 2013.
Aims. We ensure a comprehensive characterisation of the investigated systems by analysing Kepler data combined with new ground-based transit data using a photodynamical model. This method is applied to the Kepler-82 system leading to its first dynamic analysis.
Methods. In order to provide a coherent description of all observations simultaneously, we combine the numerical integration of the gravitational dynamics of a system over the time span of observations with a transit light curve model. To explore the model parameter space, this photodynamical model is coupled with a Markov chain Monte Carlo algorithm.
Results. The Kepler-82b/c system shows sinusoidal TTVs due to their near 2:1 resonance dynamical interaction. An additional chopping effect in the TTVs of Kepler-82c hints to a further planet near the 3:2 or 3:1 resonance. We photodynamically analysed Kepler long- and short-cadence data and three new transit observations obtained by KOINet between 2014 and 2018. Our result reveals a non-transiting outer planet with a mass of $m_f=20.9\pm1.0\;M_\bigoplus$ near the 3:2 resonance to the outermost known planet, Kepler-82c. Furthermore, we determined the densities of planets b and c to the significantly more precise values $ρ_b=0.98_{-0.14}^{+0.10}\;\text{g cm}^{-3}$ and $ρ_c=0.494_{-0.077}^{+0.066}\;\text{g cm}^{-3}$.
△ Less
Submitted 15 July, 2019;
originally announced July 2019.
-
Stellar Properties of Active G and K Stars: Exploring the Connection Between Starspots and Chromospheric Activity
Authors:
Brett M. Morris,
Jason L. Curtis,
Charli Sakari,
Suzanne L. Hawley,
Eric Agol
Abstract:
We gathered high resolution spectra for an ensemble of 55 bright active and inactive stars using the ARC 3.5 m Telescope Echelle Spectrograph at Apache Point Observatory ($R\approx$31,500). We measured spectroscopic effective temperatures, surface gravities and metallicities for most stars in the sample with SME and MOOG. Our stellar property results are consistent with the photometric effective t…
▽ More
We gathered high resolution spectra for an ensemble of 55 bright active and inactive stars using the ARC 3.5 m Telescope Echelle Spectrograph at Apache Point Observatory ($R\approx$31,500). We measured spectroscopic effective temperatures, surface gravities and metallicities for most stars in the sample with SME and MOOG. Our stellar property results are consistent with the photometric effective temperatures from the Gaia DR2 pipeline. We also measured their chromospheric $S$ and $\log R^\prime_\mathrm{HK}$ indices to classify the stars as active or inactive and study the connection between chromospheric activity and starspots. We then attempted to infer the starspot covering fractions on the active stars by modeling their spectra as a linear combination of hot and cool inactive stellar spectral templates. We find that it is critical to use precise colors of the stars to place stringent priors on the plausible spot covering fractions. The inferred spot covering fractions generally increase with the chromospheric activity indicator $\log R^\prime_\mathrm{HK}$, though we are largely insensitive to spot coverages smaller than $f_S \lesssim 20$%. We find a dearth of stars with small $\log R^\prime_\mathrm{HK}$ and significant spot coverages.
△ Less
Submitted 30 June, 2019;
originally announced July 2019.
-
Constraining Stellar Photospheres as an Essential Step for Transmission Spectroscopy of Small Exoplanets
Authors:
Benjamin V. Rackham,
Arazi Pinhas,
Dániel Apai,
Raphaëlle Haywood,
Heather Cegla,
Néstor Espinoza,
Johanna K. Teske,
Michael Gully-Santiago,
Gioia Rau,
Brett M. Morris,
Daniel Angerhausen,
Thomas Barclay,
Ludmila Carone,
P. Wilson Cauley,
Julien de Wit,
Shawn Domagal-Goldman,
Chuanfei Dong,
Diana Dragomir,
Mark S. Giampapa,
Yasuhiro Hasegawa,
Natalie R. Hinkel,
Renyu Hu,
Andrés Jordán,
Irina Kitiashvili,
Laura Kreidberg
, et al. (7 additional authors not shown)
Abstract:
Transmission spectra probe the atmospheres of transiting exoplanets, but these observations are also subject to signals introduced by magnetic active regions on host stars. Here we outline scientific opportunities in the next decade for providing useful constraints on stellar photospheres and inform interpretations of transmission spectra of the smallest ($R<4\,R_{\odot}$) exoplanets. We identify…
▽ More
Transmission spectra probe the atmospheres of transiting exoplanets, but these observations are also subject to signals introduced by magnetic active regions on host stars. Here we outline scientific opportunities in the next decade for providing useful constraints on stellar photospheres and inform interpretations of transmission spectra of the smallest ($R<4\,R_{\odot}$) exoplanets. We identify and discuss four primary opportunities: (1) refining stellar magnetic active region properties through exoplanet crossing events; (2) spectral decomposition of active exoplanet host stars; (3) joint retrievals of stellar photospheric and planetary atmospheric properties with studies of transmission spectra; and (4) continued visual transmission spectroscopy studies to complement longer-wavelength studies from $\textit{JWST}$. We make five recommendations to the Astro2020 Decadal Survey Committee: (1) identify the transit light source (TLS) effect as a challenge to precise exoplanet transmission spectroscopy and an opportunity ripe for scientific advancement in the coming decade; (2) include characterization of host star photospheric heterogeneity as part of a comprehensive research strategy for studying transiting exoplanets; (3) support the construction of ground-based extremely large telescopes (ELTs); (4) support multi-disciplinary research teams that bring together the heliophysics, stellar physics, and exoplanet communities to further exploit transiting exoplanets as spatial probes of stellar photospheres; and (5) support visual transmission spectroscopy efforts as complements to longer-wavelength observational campaigns with $\textit{JWST}$.
△ Less
Submitted 14 March, 2019;
originally announced March 2019.
-
The Solar Benchmark: Rotational Modulation of the Sun Reconstructed from Archival Sunspot Records
Authors:
Brett M. Morris,
James R. A. Davenport,
Helen A. C. Giles,
Leslie Hebb,
Suzanne L. Hawley,
Ruth Angus,
Peter A. Gilman,
Eric Agol
Abstract:
We use archival daily spot coverage measurements from Howard et al. (1984) to study the rotational modulation of the Sun as though it were a distant star. A quasi-periodic Gaussian process measures the solar rotation period $P_\mathrm{rot} = 26.3 \pm 0.1$ days, and activity cycle period $P_\mathrm{cyc} = 10.7 \pm 0.3$ years. We attempt to search for evidence of differential rotation in variations…
▽ More
We use archival daily spot coverage measurements from Howard et al. (1984) to study the rotational modulation of the Sun as though it were a distant star. A quasi-periodic Gaussian process measures the solar rotation period $P_\mathrm{rot} = 26.3 \pm 0.1$ days, and activity cycle period $P_\mathrm{cyc} = 10.7 \pm 0.3$ years. We attempt to search for evidence of differential rotation in variations of the apparent rotation period throughout the activity cycle and do not detect a clear signal of differential rotation, consistent with the null results of the hare-and-hounds exercise of Aigrain et al. (2015). The full reconstructed solar light curve is available online.
△ Less
Submitted 14 January, 2019;
originally announced January 2019.
-
astroquery: An Astronomical Web-Querying Package in Python
Authors:
Adam Ginsburg,
Brigitta M. Sipőcz,
C. E. Brasseur,
Philip S. Cowperthwaite,
Matthew W. Craig,
Christoph Deil,
James Guillochon,
Giannina Guzman,
Simon Liedtke,
Pey Lian Lim,
Kelly E. Lockhart,
Michael Mommert,
Brett M. Morris,
Henrik Norman,
Madhura Parikh,
Magnus V. Persson,
Thomas P. Robitaille,
Juan-Carlos Segovia,
Leo P. Singer,
Erik J. Tollerud,
Miguel de Val-Borro,
Ivan Valtchanov,
Julien Woillez,
the Astroquery collaboration
Abstract:
astroquery is a collection of tools for requesting data from databases hosted on remote servers with interfaces exposed on the internet, including those with web pages but without formal application program interfaces (APIs). These tools are built on the Python requests package, which is used to make HTTP requests, and astropy, which provides most of the data parsing functionality. astroquery modu…
▽ More
astroquery is a collection of tools for requesting data from databases hosted on remote servers with interfaces exposed on the internet, including those with web pages but without formal application program interfaces (APIs). These tools are built on the Python requests package, which is used to make HTTP requests, and astropy, which provides most of the data parsing functionality. astroquery modules generally attempt to replicate the web page interface provided by a given service as closely as possible, making the transition from browser-based to command-line interaction easy. astroquery has received significant contributions from throughout the astronomical community, including several significant contributions from telescope archives. astroquery enables the creation of fully reproducible workflows from data acquisition through publication. This paper describes the philosophy, basic structure, and development model of the astroquery package. The complete documentation for astroquery can be found at http://astroquery.readthedocs.io/.
△ Less
Submitted 14 January, 2019;
originally announced January 2019.
-
Pushing Automated Abundance Derivations Into the Cool Star Regime: A Test Using Three G and Two K Stars in Praesepe
Authors:
Marwan Gebran,
Marcel A. Ag�eros,
Keith Hawkins,
Simon C. Schuler,
Brett M. Morris
Abstract:
We present the results of an abundance analysis of three G and two K dwarfs in the Praesepe open cluster based on high-resolution, moderate signal-to-noise-ratio spectra obtained with the ARC 3.5-m Telescope at Apache Point Observatory. Using a Principle Component Analysis and the BACCHUS automated spectral analysis code, we determined stellar parameters and abundances of up to 24 elements for eac…
▽ More
We present the results of an abundance analysis of three G and two K dwarfs in the Praesepe open cluster based on high-resolution, moderate signal-to-noise-ratio spectra obtained with the ARC 3.5-m Telescope at Apache Point Observatory. Using a Principle Component Analysis and the BACCHUS automated spectral analysis code, we determined stellar parameters and abundances of up to 24 elements for each of our targets, which range in temperature from 6000 to 4600 K. The average derived iron abundance for the three G stars is 0.17+/-0.07 dex, consistent with the 0.12+/-0.04 dex derived by Boesgaard et al. (2013) for their sample of 11 solar-type Praesepe members, which included these G stars. To investigate the efficacy of using automated routines to derive the abundances of cooler main-sequence stars, we compared the abundances of the K dwarfs to those of the G dwarfs. Our abundances are consistent to <=0.1~dex for 13 of the 18 elements we report for all five of the stars, providing more evidence that G and K stars in a given open cluster are chemically homogeneous. The median difference between the mean G and K stars abundances is 0.08+/-0.05 dex, despite serious challenges with the noisier data for the fainter K dwarfs. Our results are encouraging for chemical tagging, as they indicate that it may be possible to use automated abundance determination techniques to identify chemically related main-sequence stars across larger temperature ranges than are usually considered in these experiments.
△ Less
Submitted 7 December, 2018;
originally announced December 2018.
-
Are Starspots and Plages Co-Located on Active G and K Stars?
Authors:
Brett M. Morris,
Jason L. Curtis,
Stephanie T. Douglas,
Suzanne L. Hawley,
Marcel A. Ag�eros,
Monica G. Bobra,
Eric Agol
Abstract:
We explore the connection between starspots and plages of three main-sequence stars by studying the chromospheric and photospheric activity over several rotation periods. We present simultaneous photometry and high-resolution ($R\sim 31,500$) spectroscopy of KIC 9652680, a young, superflare-producing G1 star with a rotation period of 1.4 days. Its Kepler light curve shows rotational modulation con…
▽ More
We explore the connection between starspots and plages of three main-sequence stars by studying the chromospheric and photospheric activity over several rotation periods. We present simultaneous photometry and high-resolution ($R\sim 31,500$) spectroscopy of KIC 9652680, a young, superflare-producing G1 star with a rotation period of 1.4 days. Its Kepler light curve shows rotational modulation consistent with a bright hemisphere followed by a relatively dark hemisphere, generating photometric variability with a semi-amplitude of 4%. We find that KIC 9652680 is darkest when its $S$-index of Ca II H & K emission is at its maximum. We interpret this anti-correlation between flux and $S$ to indicate that dark starspots in the photosphere are co-located with the bright plages in the chromosphere, as they are on the Sun. Moving to lower masses and slower rotators, we present K2 observations with simultaneous spectroscopy of EPIC 211928486 (K5V) and EPIC 211966629 (K4V), two active stars in the 650 Myr-old open cluster Praesepe. The K2 photometry reveals that both stars have rotation periods of 11.7 days; while their flux varies by 1 and 2% respectively, their Ca II H & K $S$-indices seem to hold relatively constant as a function of rotational phase. This suggests that extended chromospheric networks of plages are not concentrated into regions of emission centered on the starspots that drive rotational modulation, unlike KIC 9652680. We also note that the Ca II emission of EPIC 211928486 dipped and recovered suddenly over the duration of one rotation, suggesting that the evolution timescale of plages may be of order the rotation period.
△ Less
Submitted 12 September, 2018;
originally announced September 2018.
-
Non-detection of Contamination by Stellar Activity in the Spitzer Transit Light Curves of TRAPPIST-1
Authors:
Brett M. Morris,
Eric Agol,
Leslie Hebb,
Suzanne L. Hawley,
Michaël Gillon,
Elsa Ducrot,
Laetitia Delrez,
James Ingalls,
Brice-Olivier Demory
Abstract:
We apply the transit light curve self-contamination technique of Morris et al. (2018) to search for the effect of stellar activity on the transits of the ultracool dwarf TRAPPIST-1 with 2018 Spitzer photometry. The self-contamination method fits the transit light curves of planets orbiting spotted stars, allowing the host star to be a source of contaminating positive or negative flux which influen…
▽ More
We apply the transit light curve self-contamination technique of Morris et al. (2018) to search for the effect of stellar activity on the transits of the ultracool dwarf TRAPPIST-1 with 2018 Spitzer photometry. The self-contamination method fits the transit light curves of planets orbiting spotted stars, allowing the host star to be a source of contaminating positive or negative flux which influences the transit depths but not the ingress/egress durations. We find that none of the planets show statistically significant evidence for self-contamination by bright or dark regions of the stellar photosphere. However, we show that small-scale magnetic activity, analogous in size to the smallest sunspots, could still be lurking in the transit photometry undetected.
△ Less
Submitted 8 August, 2018;
originally announced August 2018.
-
Robust Transiting Exoplanet Radii in the Presence of Starspots from Ingress and Egress Durations
Authors:
Brett M. Morris,
Eric Agol,
Leslie Hebb,
Suzanne L. Hawley
Abstract:
We typically measure the radii of transiting exoplanets from the transit depth, which is given by the ratio of cross-sectional areas of the planet and star. However, if a star has dark starspots (or bright regions) distributed throughout the transit chord, the transit depth will be biased towards smaller (larger) values, and thus the inferred planet radius will be smaller (larger) if unaccounted f…
▽ More
We typically measure the radii of transiting exoplanets from the transit depth, which is given by the ratio of cross-sectional areas of the planet and star. However, if a star has dark starspots (or bright regions) distributed throughout the transit chord, the transit depth will be biased towards smaller (larger) values, and thus the inferred planet radius will be smaller (larger) if unaccounted for. We reparameterize the transit light curve to account for "self-contamination" by photospheric inhomogeneities by splitting the parameter $R_p/R_\star$ into two parameters: one for the radius ratio -- which controls the duration of ingress and egress -- and another which measures the possibly contaminated transit depth. We show that this is equivalent to the formulation for contamination by a second star (with positive or negative flux), and that it is sensitive to time-steady inhomogeneity of the stellar photosphere. We use synthetic light curves of spotted stars at high signal-to-noise to show that the radius recovered from measurement of the ingress/egress duration can recover the true radii of planets transiting spotted stars with axisymmetric spot distributions if the limb-darkening parameters are precisely known. We fit time-averaged high signal-to-noise transit light curves from Kepler and Spitzer of ten planets to measure the planet radii and search for evidence of spot distributions. We find that this sample has a range of measured depths and ingress durations which are self-consistent, providing no strong evidence for contamination by spots. However, there is suggestive evidence for occultation of starspots on Kepler-17, and that relatively bright regions are occulted by the planets of Kepler-412 and HD 80606. Future observations with the James Webb Space Telescope may enable this technique to yield accurate planetary radii in the presence of stellar inhomogeneities.
△ Less
Submitted 12 July, 2018;
originally announced July 2018.
-
The 0.8-4.5$μ$m broadband transmission spectra of TRAPPIST-1 planets
Authors:
E. Ducrot,
M. Sestovic,
B. M. Morris,
M. Gillon,
A. H. M. J. Triaud,
J. de Wit,
D. Thimmarayappa,
E. Agol,
Y. Almleaky,
A. Burdanov,
A. J. Burgasser,
L. Delrez,
B-O. Demory,
E. Jehin,
J. Leconte,
J. McCormac,
C. Murray,
D. Queloz,
F. Selsis,
S. Thompson,
V. Van Grootel
Abstract:
The TRAPPIST-1 planetary system represents an exceptional opportunity for the atmospheric characterization of temperate terrestrial exoplanets with the upcoming James Webb Space Telescope (JWST). Assessing the potential impact of stellar contamination on the planets' transit transmission spectra is an essential precursor step to this characterization. Planetary transits themselves can be used to s…
▽ More
The TRAPPIST-1 planetary system represents an exceptional opportunity for the atmospheric characterization of temperate terrestrial exoplanets with the upcoming James Webb Space Telescope (JWST). Assessing the potential impact of stellar contamination on the planets' transit transmission spectra is an essential precursor step to this characterization. Planetary transits themselves can be used to scan the stellar photosphere and to constrain its heterogeneity through transit depth variations in time and wavelength. In this context, we present our analysis of 169 transits observed in the optical from space with K2 and from the ground with the SPECULOOS and Liverpool telescopes. Combining our measured transit depths with literature results gathered in the mid/near-IR with Spitzer/IRAC and HST/WFC3, we construct the broadband transmission spectra of the TRAPPIST-1 planets over the 0.8-4.5 $μ$m spectral range. While planets b, d, and f spectra show some structures at the 200-300ppm level, the four others are globally flat. Even if we cannot discard their instrumental origins, two scenarios seem to be favored by the data: a stellar photosphere dominated by a few high-latitude giant (cold) spots, or, alternatively, by a few small and hot (3500-4000K) faculae. In both cases, the stellar contamination of the transit transmission spectra is expected to be less dramatic than predicted in recent papers. Nevertheless, based on our results, stellar contamination can still be of comparable or greater order than planetary atmospheric signals at certain wavelengths. Understanding and correcting the effects of stellar heterogeneity therefore appears essential to prepare the exploration of TRAPPIST-1's with JWST.
△ Less
Submitted 22 October, 2018; v1 submitted 3 July, 2018;
originally announced July 2018.
-
Kepler Object of Interest Network II. Photodynamical modelling of Kepler-9 over 8 years of transit observations
Authors:
J. Freudenthal,
C. von Essen,
S. Dreizler,
S. Wedemeyer,
E. Agol,
B. M. Morris,
A. C. Becker,
M. Mallonn,
S. Hoyer,
A. Ofir,
L. Tal Or,
H. J. Deeg,
E. Herrero,
I. Ribas,
S. Khalafinejad,
J. Hernández,
M. M. Rodríguez S
Abstract:
The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised to follow up transiting planet candidate KOIs with large transit timing variations (TTVs). Its main goal is to complete their TTV curves, as the Kepler telescope no longer observes the original Kepler field. Combining Kepler and new ground-based transit data we improve the modelling of t…
▽ More
The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised to follow up transiting planet candidate KOIs with large transit timing variations (TTVs). Its main goal is to complete their TTV curves, as the Kepler telescope no longer observes the original Kepler field. Combining Kepler and new ground-based transit data we improve the modelling of these systems. To this end, we have developed a photodynamical model, and we demonstrate its performance using the Kepler-9 system as an example. Our comprehensive analysis combines the numerical integration of the system's dynamics over the time span of the observations along with the transit light curve model. This model is coupled with a Markov chain Monte Carlo algorithm, allowing the exploration of the model parameter space. Applied to the Kepler-9 long cadence data, short cadence data and 13 new transit observations collected by KOINet between the years 2014 to 2017, our modelling provides well constrained predictions for the next transits and the system's parameters. We have determined the densities of the planets Kepler-9b and 9c to the very precise values of rho_b = 0.439 +/-0.023 g/cm3 and rho_c = 0.322 +/- 0.017 g/cm3. Our analysis reveals that Kepler-9c will stop transiting in about 30 years. This results from strong dynamical interactions between Kepler-9b and 9c, near 2:1 resonance, that leads to a periodic change in inclination. Over the next 30 years the inclination of Kepler-9c (-9b) will decrease (increase) slowly. This should be measurable by a substantial decrease (increase) in the transit duration, in as soon as a few years' time. Observations that contradict this prediction might indicate the presence of additional objects. If this prediction proves true, this behaviour opens up a unique chance to scan the different latitudes of a star.
△ Less
Submitted 29 June, 2018;
originally announced July 2018.
-
Possible Bright Starspots on TRAPPIST-1
Authors:
Brett M. Morris,
Eric Agol,
James R. A. Davenport,
Suzanne L. Hawley
Abstract:
The M8V star TRAPPIST-1 hosts seven roughly Earth-sized planets and is a promising target for exoplanet characterization. Kepler/K2 Campaign 12 observations of TRAPPIST-1 in the optical show an apparent rotational modulation with a 3.3 day period, though that rotational signal is not readily detected in the Spitzer light curve at 4.5 $μ$m. If the rotational modulation is due to starspots, persiste…
▽ More
The M8V star TRAPPIST-1 hosts seven roughly Earth-sized planets and is a promising target for exoplanet characterization. Kepler/K2 Campaign 12 observations of TRAPPIST-1 in the optical show an apparent rotational modulation with a 3.3 day period, though that rotational signal is not readily detected in the Spitzer light curve at 4.5 $μ$m. If the rotational modulation is due to starspots, persistent dark spots can be excluded from the lack of photometric variability in the Spitzer light curve. We construct a photometric model for rotational modulation due to photospheric bright spots on TRAPPIST-1 which is consistent with both the Kepler and Spitzer light curves. The maximum-likelihood model with three spots has typical spot sizes of $R_\mathrm{spot}/R_\star \approx 0.004$ at temperature $T_\mathrm{spot} \gtrsim 5300 \pm 200$ K. We also find that large flares are observed more often when the brightest spot is facing the observer, suggesting a correlation between the position of the bright spots and flare events. In addition, these flares may occur preferentially when the spots are increasing in brightness, which suggests that the 3.3 d periodicity may not be a rotational signal, but rather a characteristic timescale of active regions.
△ Less
Submitted 12 March, 2018;
originally announced March 2018.