Search

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., Comment: Accepted in A&A

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., Comment: 38 pages, 34 figures, published in A&A on October 24, 2023

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 $\rho=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}$., Comment: Published in ApJL

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., Comment: 15 pages, 10 figures. Accepted to MNRAS. Repository with source code: https://github.com/ekaterinailin/flare-locations-ensembles-science Repository with data: https://zenodo.org/record/7996929

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, which have been the favored targets of eight JWST Cycle 1 programs. While 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 constrain stellar heterogeneities and gather up to 2$\times$ more transits per JWST hour spent. We conclude that in such systems planets cannot be studied in isolation by small programs, thus large-scale community-supported programs should be supported to enable the efficient and robust exploration of terrestrial exoplanets in the JWST era.

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., Comment: 31 pages, 24 figures, accepted for publication at MNRAS

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 $\mu$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)., Comment: Corrected typos; revised Section 3, 4, and 5 to reflect reanalysis, results unchanged. Submitted to AAS Journals Nov 11th, 2020; favorable referee report received Jan 3rd; final draft accepted for publication in the AJ Apr 19th

Study Analysis Group 21 (SAG21) of the 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, and plages. This SAG has 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 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")., Comment: 122 pages, 30 figures

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.

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.

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., Comment: 43 pages, 5 figures. To appear in ApJ. The author list has two parts: the authors that made significant contributions to the writing and/or coordination of the paper, followed by maintainers of and contributors to the Astropy Project. The position in the author list does not correspond to contributions to the Astropy Project as a whole

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., Comment: Consistent with published version which contains significantly improved detection and rejection of contaminant objects. Comments welcome

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., Comment: 16 pages, 17 figures; Accepted to AJ

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., Comment: 31 pages, 24 figures, accepted for publication at MNRAS

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 $\mu$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)., Comment: Corrected typos; revised Section 3, 4, and 5 to reflect reanalysis, results unchanged. Submitted to AAS Journals Nov 11th, 2020; favorable referee report received Jan 3rd; final draft accepted for publication in the AJ Apr 19th

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., Comment: Invited review in press at RASTI. Based on the ExoPAG SAG21 report (arXiv:2201.09905v1) and refined via feedback from three reviewers. 75 pages, 30 figures, 5 tables

General circulation models (GCMs) provide context for interpreting multiwavelength, multiphase 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 deg on the sphere and with 40 radial points) with multiwavelength radiative transfer (30 k-table bins) running for 3000 Earth days (864 000 time-steps) takes about 19–26 d to complete depending on the type of GPU

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., Comment: Published in Nature Astronomy on 27 January 2022, accepted on 1 December 2021 (32 pages, 21 figures, 3 tables)

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., Comment: Accepted in A&A

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., Comment: Proceedings of the Small Satellite Conference, Science/Mission Payloads, SSC21-VI-02 (2021)

We obtained high-resolution spectra of the ultra-cool M-dwarf TRAPPIST-1 during the transit of its planet `b' using two high dispersion near-infrared spectrographs, IRD instrument on the Subaru 8.2m telescope and HPF instrument on the 10m Hobby-Eberly Telescope. These spectroscopic observations are complemented by a photometric transit observation for planet `b' using the APO/ARCTIC, which assisted us to capture the correct transit times for our transit spectroscopy. Using the data obtained by the new IRD and HPF observations, as well as the prior transit observations of planets `b', `e' and `f' from IRD, we attempt to constrain the atmospheric escape of the planet using the He I triplet 10830 {\AA} absorption line. We do not detect evidence for any primordial extended H-He atmospheres in all three planets. To limit any planet related absorption, we place an upper limit on the equivalent widths of <7.754 m{\AA} for planet `b', <10.458 m{\AA} for planet `e', and <4.143 m{\AA} for planet `f' at 95% confidence from the IRD data, and <3.467 m{\AA} for planet `b' at 95% confidence from HPF data. Using these limits along with a solar-like composition isothermal Parker wind model, we attempt to constrain the mass-loss rates for the three planets. For TRAPPIST-1b, our models exclude the highest possible energy-limited rate for a wind temperature <5000 K. This non-detection of extended atmospheres having low mean-molecular weight in all three planets aids in further constraining their atmospheric composition by steering the focus towards the search of high molecular weight species in their atmospheres., Comment: 11 pages; 4 figures; Accepted for publication in AJ

"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., Comment: Published in the Journal of Open Source Software. Comments (still) welcome. Software available at https://docs.exoplanet.codes

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., Comment: 10 pages, accepted as a Letter to the Editor in A&A

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 CHaracterising ExOPlanet Satellite (CHEOPS). The targets are relatively faint for CHEOPS, 11 mag1000 km.

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., Comment: Published in Nature Astronomy on 27 January 2022, accepted on 1 December 2021 (32 pages, 21 figures, 3 tables)

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., Comment: Accepted in A&A

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., Comment: Proceedings of the Small Satellite Conference, Science/Mission Payloads, SSC21-VI-02 (2021)

We obtained high-resolution spectra of the ultra-cool M-dwarf TRAPPIST-1 during the transit of its planet `b' using two high dispersion near-infrared spectrographs, IRD instrument on the Subaru 8.2m telescope and HPF instrument on the 10m Hobby-Eberly Telescope. These spectroscopic observations are complemented by a photometric transit observation for planet `b' using the APO/ARCTIC, which assisted us to capture the correct transit times for our transit spectroscopy. Using the data obtained by the new IRD and HPF observations, as well as the prior transit observations of planets `b', `e' and `f' from IRD, we attempt to constrain the atmospheric escape of the planet using the He I triplet 10830 {\AA} absorption line. We do not detect evidence for any primordial extended H-He atmospheres in all three planets. To limit any planet related absorption, we place an upper limit on the equivalent widths of <7.754 m{\AA} for planet `b', <10.458 m{\AA} for planet `e', and <4.143 m{\AA} for planet `f' at 95% confidence from the IRD data, and <3.467 m{\AA} for planet `b' at 95% confidence from HPF data. Using these limits along with a solar-like composition isothermal Parker wind model, we attempt to constrain the mass-loss rates for the three planets. For TRAPPIST-1b, our models exclude the highest possible energy-limited rate for a wind temperature <5000 K. This non-detection of extended atmospheres having low mean-molecular weight in all three planets aids in further constraining their atmospheric composition by steering the focus towards the search of high molecular weight species in their atmospheres., Comment: 11 pages; 4 figures; Accepted for publication in AJ

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., Comment: 10 pages, accepted as a Letter to the Editor in A&A

"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., Comment: Published in the Journal of Open Source Software. Comments (still) welcome. Software available at https://docs.exoplanet.codes

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}$., Comment: 12 pages, 8 figures, 2 tables. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://www.nature.com/articles/s41550-021-01444-7. Equations (1) to (43) have not changed since v1. Differences between v1 and this version (v4) are in Figure 4 and the inferred parameter values associated with it

Astronomers have discovered thousands of planets outside the Solar System1, most of which orbit stars that will eventually evolve into red giants and then into white dwarfs. During the red giant phase, any close-orbiting planets will be engulfed by the star2, but more distant planets can survive this phase and remain in orbit around the white dwarf3,4. Some white dwarfs show evidence for rocky material floating in their atmospheres5, in warm debris disks6–9 or orbiting very closely10–12, which has been interpreted as the debris of rocky planets that were scattered inwards and tidally disrupted13. Recently, the discovery of a gaseous debris disk with a composition similar to that of ice giant planets14 demonstrated that massive planets might also find their way into tight orbits around white dwarfs, but it is unclear whether these planets can survive the journey. So far, no intact planets have been detected in close orbits around white dwarfs. Here we report the observation of a giant planet candidate transiting the white dwarf WD 1856+534 (TIC 267574918) every 1.4 days. We observed and modelled the periodic dimming of the white dwarf caused by the planet candidate passing in front of the star in its orbit. The planet candidate is roughly the same size as Jupiter and is no more than 14 times as massive (with 95 per cent confidence). Other cases of white dwarfs with close brown dwarf or stellar companions are explained as the consequence of common-envelope evolution, wherein the original orbit is enveloped during the red giant phase and shrinks owing to friction. In this case, however, the long orbital period (compared with other white dwarfs with close brown dwarf or stellar companions) and low mass of the planet candidate make common-envelope evolution less likely. Instead, our findings for the WD 1856+534 system indicate that giant planets can be scattered into tight orbits without being tidally disrupted, motivating the search for smaller transiting planets around

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., Comment: 12 pages, 8 figures, accepted to AJ

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., Comment: 11 pages, 8 figures, accepted in MNRAS

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., Comment: 21 pages, 11 figures, accepted in ApJ, software in review at JOSS https://github.com/openjournals/joss-reviews/issues/2103

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., Comment: 37 pages, 12 figures, 1 table. Accepted for publication in The Astrophysical Journal

Magnetic activity in stars manifests as dark spots on their surfaces that modulate the brightness observed by telescopes. These light curves contain important information on stellar rotation. However, the accurate estimation of rotation periods is computationally expensive due to scarce ground truth information, noisy data, and large parameter spaces that lead to degenerate solutions. We harness the power of deep learning and successfully apply Convolutional Neural Networks to regress stellar rotation periods from Kepler light curves. Geometry-preserving time-series to image transformations of the light curves serve as inputs to a ResNet-18 based architecture which is trained through transfer learning. The McQuillan catalog of published rotation periods is used as ansatz to groundtruth. We benchmark the performance of our method against a random forest regressor, a 1D CNN, and the Auto-Correlation Function (ACF) - the current standard to estimate rotation periods. Despite limiting our input to fewer data points (1k), our model yields more accurate results and runs 350 times faster than ACF runs on the same number of data points and 10,000 times faster than ACF runs on 65k data points. With only minimal feature engineering our approach has impressive accuracy, motivating the application of deep learning to regress stellar parameters on an even larger scale, Comment: Accepted for the Machine Learning and the Physical Sciences Workshop, NeurIPS 2020

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., Comment: Final version to be published in the Planetary Sciences Journal. 56 pages, 30 figures. Data from the paper and a complete table of forecast JWST times may be found at https://github.com/ericagol/TRAPPIST1_Spitzer

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

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 $\alpha$ Cygni variables in our sample, of which three are newly discovered., Comment: 29 pages, 13 figures, 8 tables. Accepted for publication in ApJ. Comments welcome

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., Comment: 12 pages, 8 figures, accepted to AJ

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., Comment: 11 pages, 8 figures, accepted in MNRAS

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., Comment: 21 pages, 11 figures, accepted in ApJ, software in review at JOSS https://github.com/openjournals/joss-reviews/issues/2103

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., Comment: 37 pages, 12 figures, 1 table. Accepted for publication in The Astrophysical Journal

Magnetic activity in stars manifests as dark spots on their surfaces that modulate the brightness observed by telescopes. These light curves contain important information on stellar rotation. However, the accurate estimation of rotation periods is computationally expensive due to scarce ground truth information, noisy data, and large parameter spaces that lead to degenerate solutions. We harness the power of deep learning and successfully apply Convolutional Neural Networks to regress stellar rotation periods from Kepler light curves. Geometry-preserving time-series to image transformations of the light curves serve as inputs to a ResNet-18 based architecture which is trained through transfer learning. The McQuillan catalog of published rotation periods is used as ansatz to groundtruth. We benchmark the performance of our method against a random forest regressor, a 1D CNN, and the Auto-Correlation Function (ACF) - the current standard to estimate rotation periods. Despite limiting our input to fewer data points (1k), our model yields more accurate results and runs 350 times faster than ACF runs on the same number of data points and 10,000 times faster than ACF runs on 65k data points. With only minimal feature engineering our approach has impressive accuracy, motivating the application of deep learning to regress stellar parameters on an even larger scale, Comment: Accepted for the Machine Learning and the Physical Sciences Workshop, NeurIPS 2020

Astronomers have discovered thousands of planets outside the Solar System1, most of which orbit stars that will eventually evolve into red giants and then into white dwarfs. During the red giant phase, any close-orbiting planets will be engulfed by the star2, but more distant planets can survive this phase and remain in orbit around the white dwarf3,4. Some white dwarfs show evidence for rocky material floating in their atmospheres5, in warm debris disks6–9 or orbiting very closely10–12, which has been interpreted as the debris of rocky planets that were scattered inwards and tidally disrupted13. Recently, the discovery of a gaseous debris disk with a composition similar to that of ice giant planets14 demonstrated that massive planets might also find their way into tight orbits around white dwarfs, but it is unclear whether these planets can survive the journey. So far, no intact planets have been detected in close orbits around white dwarfs. Here we report the observation of a giant planet candidate transiting the white dwarf WD 1856+534 (TIC 267574918) every 1.4 days. We observed and modelled the periodic dimming of the white dwarf caused by the planet candidate passing in front of the star in its orbit. The planet candidate is roughly the same size as Jupiter and is no more than 14 times as massive (with 95 per cent confidence). Other cases of white dwarfs with close brown dwarf or stellar companions are explained as the consequence of common-envelope evolution, wherein the original orbit is enveloped during the red giant phase and shrinks owing to friction. In this case, however, the long orbital period (compared with other white dwarfs with close brown dwarf or stellar companions) and low mass of the planet candidate make common-envelope evolution less likely. Instead, our findings for the WD 1856+534 system indicate that giant planets can be scattered into tight orbits without being tidally disrupted, motivating the search for smaller transiting planets around

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 $\mu$m transmission spectrum that is featureless at low resolution., Comment: Published in AJ, updated with proof corrections. 17 pages, 8 figures

We validate the discovery of a 2 Earth radii sub-Neptune-size planet around the nearby high proper motion M2.5-dwarf G 9-40 (EPIC 212048748), using high-precision near-infrared (NIR) radial velocity (RV) observations with the Habitable-zone Planet Finder (HPF), precision diffuser-assisted ground-based photometry with a custom narrow-band photometric filter, and adaptive optics imaging. At a distance of $d=27.9\mathrm{pc}$, G 9-40b is the second closest transiting planet discovered by K2 to date. The planet's large transit depth ($\sim$3500ppm), combined with the proximity and brightness of the host star at NIR wavelengths (J=10, K=9.2) makes G 9-40b one of the most favorable sub-Neptune-sized planet orbiting an M-dwarf for transmission spectroscopy with JWST, ARIEL, and the upcoming Extremely Large Telescopes. The star is relatively inactive with a rotation period of $\sim$29 days determined from the K2 photometry. To estimate spectroscopic stellar parameters, we describe our implementation of an empirical spectral matching algorithm using the high-resolution NIR HPF spectra. Using this algorithm, we obtain an effective temperature of $T_{\mathrm{eff}}=3404\pm73$K, and metallicity of $\mathrm{[Fe/H]}=-0.08\pm0.13$. Our RVs, when coupled with the orbital parameters derived from the transit photometry, exclude planet masses above $11.7 M_\oplus$ with 99.7% confidence assuming a circular orbit. From its radius, we predict a mass of $M=5.0^{+3.8}_{-1.9} M_\oplus$ and an RV semi-amplitude of $K=4.1^{+3.1}_{-1.6}\mathrm{m\:s^{-1}}$, making its mass measurable with current RV facilities. We urge further RV follow-up observations to precisely measure its mass, to enable precise transmission spectroscopic measurements in the future., Comment: Accepted for publication in AJ, 22 pages, 15 figures

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., Comment: Accepted for publication in AJ

The 2020s will be the most data-rich decade of astronomy in history. As the scale and complexity of our surveys increase, the problem of scheduling becomes more critical. We must develop high-quality scheduling approaches, implement them as open-source software, and begin linking the typically separate stages of observation and data analysis., Comment: 10 pages; Astro2020 APC White Paper