Showing 1–50 of 54 results for author: Luger, R

Dissecting the Gaia HR diagram II. The vertical structure of the star formation history across the Solar Cylinder

Authors: Alessandro Mazzi, Léo Girardi, Michele Trabucchi, Julianne J. Dalcanton, Rodrigo Luger, Paola Marigo, Andrea Miglio, Guglielmo Costa, Yang Chen, Giada Pastorelli, Morgan Fouesneau, Simone Zaggia, Alessandro Bressan, Piero Dal Tio

Abstract: Starting from the Gaia DR3 HR diagram, we derive the star formation history (SFH) as a function of distance from the Galactic Plane within a cylinder centred on the Sun with a 200~pc radius and spanning 1.3~kpc above and below the Galaxy's midplane. We quantify both the concentration of the more recent star formation in the Galactic Plane, and the age-related increase in the scale height of the Ga… ▽ More Starting from the Gaia DR3 HR diagram, we derive the star formation history (SFH) as a function of distance from the Galactic Plane within a cylinder centred on the Sun with a 200~pc radius and spanning 1.3~kpc above and below the Galaxy's midplane. We quantify both the concentration of the more recent star formation in the Galactic Plane, and the age-related increase in the scale height of the Galactic Disc stellar component, which is well-described by power-laws with indices ranging from $1/2$ to $2/3$. The vertically-integrated star formation rate falls from $(1.147 \pm 0.039)\times10^{-8}\, \text{M}_{\odot} \text{yr}^{-1} \text{pc}^{-2}$ at earlier times down to $(6.2 \pm 3.0) \times10^{-9}\, \text{M}_{\odot} \text{yr}^{-1} \text{pc}^{-2}$ at present times, but we find a significant peak of star formation in the 2 to 3 Gyr age bin. The total mass of stars formed per unit area over time is $118.7 \pm 6.2\, \text{M}_{\odot} \text{pc}^{-2}$, which is nearly twice the present stellar mass derived from kinematics within 1~kpc from the Galactic Plane, implying a high degree of matter recycling in successive generations of stars. The method is then modified by adopting an age-dependent correlation between the SFH across the different slices, which results in less noisy and more symmetrical results without significantly changing the previously mentioned quantities. This appears to be a promising way to improve SFH recovery in external galaxies. △ Less

Submitted 23 September, 2023; originally announced September 2023.

Comments: Accepted for publication in MNRAS

Warm Jupiters Beyond the Tidal Synchronization Limit May Exhibit a Wide Range of Secondary Eclipse Depths

Authors: Emily Rauscher, Nicolas B. Cowan, Rodrigo Luger

Abstract: With JWST we can now characterize the atmospheres of planets on longer orbital planets, but this moves us into a regime where we cannot assume that tidal forces from the star have eroded planets' obliquities and synchronized their rotation rates. These rotation vectors may be tracers of formation and evolution histories and also enable a range of atmospheric circulation states. Here we delineate t… ▽ More With JWST we can now characterize the atmospheres of planets on longer orbital planets, but this moves us into a regime where we cannot assume that tidal forces from the star have eroded planets' obliquities and synchronized their rotation rates. These rotation vectors may be tracers of formation and evolution histories and also enable a range of atmospheric circulation states. Here we delineate the orbital space over which tidal synchronization and alignment assumptions may no longer apply and present three-dimensional atmospheric models of a hypothetical warm Jupiter over a range of rotation rates and obliquities. We simulate the secondary eclipses of this planet for different possible viewing orientations and times during its orbital, seasonal cycle. We find that the eclipse depth can be strongly influenced by rotation rate and obliquity through the timing of the eclipse relative to the planet's seasonal cycle, and advise caution in attempting to derive properties such as albedo or day-night transport from this measurement. We predict that if warm Jupiters beyond the tidal limit have intrinsic diversity in their rotation vectors, then it will manifest itself as dispersion in their secondary eclipse depths. We explore eclipse mapping as a way to uniquely constrain the rotation vector of warm Jupiters but find that the associated signals are likely at the edge of JWST performance. Nevertheless, as JWST begins to measure the secondary eclipses of longer orbital period planets, we should expect to observe the consequences of a wider range of rotation states and circulation patterns. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 14 pages, 6 figures, animated figures available in the journal's online version

Journal ref: Astronomical Journal (2023), 165 261

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.

Comments: 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

Joint Modeling of Radial Velocities and Photometry with a Gaussian Process Framework

Authors: Quang H. Tran, Megan Bedell, Daniel Foreman-Mackey, Rodrigo Luger

Abstract: Developments in the stability of modern spectrographs have led to extremely precise instrumental radial velocity (RV) measurements. For most stars, the detection limit of planetary companions with these instruments is expected to be dominated by astrophysical noise sources such as starspots. Correlated signals caused by rotationally-modulated starspots can obscure or mimic the Doppler shifts induc… ▽ More Developments in the stability of modern spectrographs have led to extremely precise instrumental radial velocity (RV) measurements. For most stars, the detection limit of planetary companions with these instruments is expected to be dominated by astrophysical noise sources such as starspots. Correlated signals caused by rotationally-modulated starspots can obscure or mimic the Doppler shifts induced by even the closest, most massive planets. This is especially true for young, magnetically active stars where stellar activity can cause fluctuation amplitudes of $\gtrsim$0.1 mag in brightness and $\gtrsim$100 m s$^{-1}$ in RV semi-amplitudes. Techniques that can mitigate these effects and increase our sensitivity to young planets are critical to improving our understanding of the evolution of planetary systems. Gaussian processes (GPs) have been successfully employed to model and constrain activity signals in individual cases. However, a principled approach of this technique, specifically for the joint modeling of photometry and RVs, has not yet been developed. In this work, we present a GP framework to simultaneously model stellar activity signals in photometry and RVs that can be used to investigate the relationship between both time series. Our method, inspired by the $\textit{FF}^\prime$ framework of (Aigrain et al. 2012), models spot-driven activity signals as the linear combinations of two independent latent GPs and their time derivatives. We also simulate time series affected by starspots by extending the $\texttt{starry}$ software (Luger et al. 2019) to incorporate time evolution of stellar features. Using these synthetic datasets, we show that our method can predict spot-driven RV variations with greater accuracy than other GP approaches. △ Less

Submitted 1 May, 2023; originally announced May 2023.

Comments: 19 pages, 10 figures

Backward Population Synthesis: Mapping the Evolutionary History of Gravitational-Wave Progenitors

Authors: Kaze W. K. Wong, Katelyn Breivik, Will M. Farr, Rodrigo Luger

Abstract: One promising way to extract information about stellar astrophysics from gravitational wave catalogs is to compare the catalog to the outputs of stellar population synthesis modeling with varying physical assumptions. The parameter space of physical assumptions in population synthesis is high-dimensional and the choice of parameters that best represents the evolution of a binary system may depend… ▽ More One promising way to extract information about stellar astrophysics from gravitational wave catalogs is to compare the catalog to the outputs of stellar population synthesis modeling with varying physical assumptions. The parameter space of physical assumptions in population synthesis is high-dimensional and the choice of parameters that best represents the evolution of a binary system may depend in an as-yet-to-be-determined way on the system's properties. Here we propose a pipeline to simultaneously infer zero-age main sequence properties and population synthesis parameter settings controlling modeled binary evolution from individual gravitational wave observations of merging compact binaries. Our pipeline can efficiently explore the high-dimensional space of population synthesis settings and progenitor system properties for each system in a catalog of gravitational wave observations. We apply our pipeline to observations in the third third LIGO-Virgo Gravitational-Wave Transient Catalog. We showcase the effectiveness of this pipeline with a detailed study of the progenitor properties and population synthesis settings that produce mergers like the observed GW150914. Our pipeline permits a measurement of the variation of population synthesis parameter settings with binary properties, if any; we present inferences for the recent GWTC-3 transient catalog that suggest that the stable mass transfer efficiency parameter may vary with primary black hole mass. △ Less

Submitted 8 June, 2022; originally announced June 2022.

Comments: Source code for compiling the paper is publicly hosted on github: https://github.com/kazewong/BackPop. Data is available on zenodo: 10.5281/zenodo.6624911

The Contribution of M-Dwarf Flares to the Thermal Escape of Potentially Habitable Planet Atmospheres

Authors: Laura Neves Ribeiro do Amaral, Rory Barnes, Antígona Segura, Rodrigo Luger

Abstract: The habitability of planets around M dwarfs ($\lesssim 0.5 M_\odot$) can be affected by the XUV (X rays + extreme UV) emission of these stars, with flares occasionally increasing the XUV flux by more than 2 orders of magnitude above quiescent levels. This wavelength range can warm and ionize terrestrial planets' upper atmospheres, which expands the planetary radius and promotes atmospheric loss. I… ▽ More The habitability of planets around M dwarfs ($\lesssim 0.5 M_\odot$) can be affected by the XUV (X rays + extreme UV) emission of these stars, with flares occasionally increasing the XUV flux by more than 2 orders of magnitude above quiescent levels. This wavelength range can warm and ionize terrestrial planets' upper atmospheres, which expands the planetary radius and promotes atmospheric loss. In this work, we study the contribution of the XUV flux due to flares on the atmospheric escape of Earth-like planets orbiting M dwarfs through numerical simulations. We considered the first Gyr of planets with initial surface water abundances between 1 and 10 terrestrial oceans (TO), a small primordial hydrogen envelope ($\le$ $10^{-3}$~$M_{\oplus}$), and with host star masses between 0.2 and 0.6 $M_{\odot}$. In this parameter range, we find that flares can remove up to two TO more than nonflaring stars, which, in some cases, translates to a doubling of the total water loss. We also find that flaring can increase atmospheric oxygen partial pressures by hundreds of bars in some cases. These results were obtained by adding a new module for flares to the \vplanet software package and upgrading its atmospheric escape module to account for Roche lobe overflow and radiation/recombination-limited escape. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: 24 pages, 13 figures

Hierarchical Bayesian Atmospheric Retrieval Modeling for Population Studies of Exoplanet Atmospheres: A Case Study on the Habitable Zone

Authors: Jacob Lustig-Yaeger, Kristin S. Sotzen, Kevin B. Stevenson, Rodrigo Luger, Erin M. May, L. C. Mayorga, Kathleen Mandt, Noam R. Izenberg

Abstract: With the growing number of spectroscopic observations and observational platforms capable of exoplanet atmospheric characterization, there is a growing need for analysis techniques that can distill information about a large population of exoplanets into a coherent picture of atmospheric trends expressed within the statistical sample. In this work, we develop a Hierarchical Bayesian Atmospheric Ret… ▽ More With the growing number of spectroscopic observations and observational platforms capable of exoplanet atmospheric characterization, there is a growing need for analysis techniques that can distill information about a large population of exoplanets into a coherent picture of atmospheric trends expressed within the statistical sample. In this work, we develop a Hierarchical Bayesian Atmospheric Retrieval (HBAR) model to infer population-level trends in exoplanet atmospheric characteristics. We demonstrate HBAR on the case of inferring a trend in atmospheric CO2 with incident stellar flux, predicted by the presence of a functioning carbonate-silicate weathering negative feedback cycle, an assumption upon which all calculations of the habitable zone (HZ) rest. Using simulated transmission spectra and JWST-quality observations of rocky planets with H2O, CO2, and N2 bearing atmospheres, we find that the predicted trend in CO2 causes subtle differences in the spectra of order 10 ppm in the 1-5 um range, underscoring the challenge inherent to testing this hypothesis. In the limit of highly precise data (100 stacked transits per planet), we show that our HBAR model is capable of inferring the population-level parameters that characterize the trend in CO2, and we demonstrate that the null hypothesis and other simpler trends can be rejected at high confidence. Although we find that this specific empirical test of the HZ may be prohibitively challenging in the JWST era, the HBAR framework developed in this work may find a more immediate usage for the analysis of gas giant spectra observed with JWST, Ariel, and other upcoming missions. △ Less

Submitted 1 February, 2022; originally announced February 2022.

Comments: 15 pages, 5 figures, accepted for publication in AJ

Applying the metallicity-dependent binary fraction to double white dwarf formation: Implications for LISA

Authors: Sarah Thiele, Katelyn Breivik, Robyn E. Sanderson, Rodrigo Luger

Abstract: Short-period double white dwarf (DWD) binaries will be the most prolific source of gravitational waves (GWs) for the Laser Interferometer Space Antenna (LISA). DWDs with GW frequencies below $\sim1$ mHz will be the dominant contributor to a stochastic foreground caused by overlapping GW signals. Population modeling of Galactic DWDs typically assumes a binary fraction of 50% and a log-uniform Zero… ▽ More Short-period double white dwarf (DWD) binaries will be the most prolific source of gravitational waves (GWs) for the Laser Interferometer Space Antenna (LISA). DWDs with GW frequencies below $\sim1$ mHz will be the dominant contributor to a stochastic foreground caused by overlapping GW signals. Population modeling of Galactic DWDs typically assumes a binary fraction of 50% and a log-uniform Zero Age Main Sequence (ZAMS) orbital period distribution. However, recent observations have shown that the binary fraction of close, solar-type stars exhibits a strong anti-correlation with metallicity which modulates the ZAMS orbital period distribution below $10^4$ days. In this study we perform the first simulation of the Galactic DWD population observable by LISA which incorporates an empirically-derived metallicity-dependent binary fraction, using the binary population synthesis suite COSMIC and a metallicity-dependent star formation history. We compare two models: one which assumes a metallicity-dependent binary fraction, and one with a binary fraction of 50%. We repeat our analysis for three different assumptions for Roche-lobe overflow interactions. We find that while metallicity impacts the evolution and intrinsic properties of our simulated DWD progenitor binaries, the LISA-resolvable populations of the two models remain roughly indistinguishable. However, the size of the total Galactic DWD population orbiting in the LISA frequency band is reduced by more than half when accounting for a metallicity-dependent binary fraction for two of our four variations, which also lowers the effective foreground. The LISA population remains unchanged in number for two variations, highlighting the sensitivity of the population to binary evolution prescriptions. △ Less

Submitted 16 December, 2022; v1 submitted 26 November, 2021; originally announced November 2021.

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

An Aligned Orbit for the Young Planet V1298 Tau b

Authors: Marshall C. Johnson, Trevor J. David, Erik A. Petigura, Howard T. Isaacson, Judah Van Zandt, Ilya Ilyin, Klaus Strassmeier, Matthias Mallonn, George Zhou, Andrew W. Mann, John H. Livingston, Rodrigo Luger, Fei Dai, Lauren M. Weiss, Teo Močnik, Steven Giacalone, Michelle L. Hill, Malena Rice, Sarah Blunt, Ryan Rubenzahl, Paul A. Dalba, Gilbert A. Esquerdo, Perry Berlind, Michael L. Calkins, Daniel Foreman-Mackey

Abstract: The alignment of planetary orbits with respect to the stellar rotation preserves information on their dynamical histories. Measuring this angle for young planets help illuminate the mechanisms that create misaligned orbits for older planets, as different processes could operate over timescales ranging from a few Myr to a Gyr. We present spectroscopic transit observations of the young exoplanet V12… ▽ More The alignment of planetary orbits with respect to the stellar rotation preserves information on their dynamical histories. Measuring this angle for young planets help illuminate the mechanisms that create misaligned orbits for older planets, as different processes could operate over timescales ranging from a few Myr to a Gyr. We present spectroscopic transit observations of the young exoplanet V1298 Tau b; we update the age of V1298 Tau to be $28\pm4$ Myr based on Gaia EDR3 measurements. We observed a partial transit with Keck/HIRES and LBT/PEPSI, and detected the radial velocity anomaly due to the Rossiter-McLaughlin effect. V1298 Tau~b has a prograde, well-aligned orbit, with $λ= 4_{-10}^{+7 \circ}$. By combining the spectroscopically-measured $v\sin i_{\star}$ and the phtometrically-measured rotation period of the host star we also find that the orbit is aligned in 3D, $ψ= 8_{-7}^{+4 \circ}$ deg. Finally, we combine our obliquity constraints with a previous measurement for the interior planet V1298 Tau c to constrain the mutual inclination between the two planets to be $i_{\mathrm{mut}}=0^{\circ} \pm 19^{\circ}$. This measurements adds to the growing number of well-aligned planets at young ages, hinting that misalignments may be generated over timescales of longer than tens of Myr. The number of measurements, however, is still small, and this population may not be representative of the older planets that have been observed to date. We also present the derivation of the relationship between $i_{\mathrm{mut}}$, $λ$, and $i$ for two planets. △ Less

Submitted 17 March, 2022; v1 submitted 20 October, 2021; originally announced October 2021.

Comments: Accepted for publication in AJ. 17 pages, 10 figures

Mapping stellar surfaces III: An Efficient, Scalable, and Open-Source Doppler Imaging Model

Authors: Rodrigo Luger, Megan Bedell, Daniel Foreman-Mackey, Ian J. M. Crossfield, Lily L. Zhao, David W. Hogg

Abstract: The study of stellar surfaces can reveal information about the chemical composition, interior structure, and magnetic properties of stars. It is also critical to the detection and characterization of extrasolar planets, in particular those targeted in extreme precision radial velocity (EPRV) searches, which must contend with stellar variability that is often orders of magnitude stronger than the p… ▽ More The study of stellar surfaces can reveal information about the chemical composition, interior structure, and magnetic properties of stars. It is also critical to the detection and characterization of extrasolar planets, in particular those targeted in extreme precision radial velocity (EPRV) searches, which must contend with stellar variability that is often orders of magnitude stronger than the planetary signal. One of the most successful methods to map the surfaces of stars is Doppler imaging, in which the presence of inhomogeneities is inferred from subtle line shape changes in high resolution stellar spectra. In this paper, we present a novel, efficient, and closed-form solution to the problem of Doppler imaging of stellar surfaces. Our model explicitly allows for incomplete knowledge of the local (rest frame) stellar spectrum, allowing one to learn differences from spectral templates while simultaneously mapping the stellar surface. It therefore works on blended lines, regions of the spectrum where line formation mechanisms are not well understood, or stars whose spots have intrinsically different spectra from the rest of the photosphere. We implement the model within the open source starry framework, making it fast, differentiable, and easy to use in both optimization and posterior inference settings. As a proof-of-concept, we use our model to infer the surface map of the brown dwarf WISE 1049-5319B, finding close agreement with the solution of Crossfield et al. (2014). We also discuss Doppler imaging in the context of EPRV studies and describe an interpretable spectral-temporal Gaussian process for stellar spectral variability that we expect will be important for EPRV exoplanet searches. △ Less

Submitted 12 October, 2021; originally announced October 2021.

Comments: 59 pages, 18 figures. Prepared using the showyourwork (github.com/rodluger/showyourwork) open source scientific article workflow

Efficient and precise transit light curves for rapidly-rotating, oblate stars

Authors: Shashank Dholakia, Rodrigo Luger, Shishir Dholakia

Abstract: We derive solutions to transit light curves of exoplanets orbiting rapidly-rotating stars. These stars exhibit significant oblateness and gravity darkening, a phenomenon where the poles of the star have a higher temperature and luminosity than the equator. Light curves for exoplanets transiting these stars can exhibit deviations from those of slowly-rotating stars, even displaying significantly as… ▽ More We derive solutions to transit light curves of exoplanets orbiting rapidly-rotating stars. These stars exhibit significant oblateness and gravity darkening, a phenomenon where the poles of the star have a higher temperature and luminosity than the equator. Light curves for exoplanets transiting these stars can exhibit deviations from those of slowly-rotating stars, even displaying significantly asymmetric transits depending on the system's spin-orbit angle. As such, these phenomena can be used as a protractor to measure the spin-orbit alignment of the system. In this paper, we introduce a novel semi-analytic method for generating model light curves for gravity-darkened and oblate stars with transiting exoplanets. We implement the model within the code package starry and demonstrate several orders of magnitude improvement in speed and precision over existing methods. We test the model on a TESS light curve of WASP-33, whose host star displays rapid rotation ($v \sin i_* = 86.4$ km/s). We subtract the host's $δ$-Scuti pulsations from the light curve, finding an asymmetric transit characteristic of gravity darkening. We find the projected spin orbit angle is consistent with Doppler tomography and constrain the true spin-orbit angle of the system as $\varphi=108.3^{+19.0}_{-15.4}$~$^{\circ}$. We demonstrate the method's uses in constraining spin-orbit inclinations of such systems photometrically with posterior inference. Lastly, we note the use of such a method for inferring the dynamical history of thousands of such systems discovered by TESS. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: 33 pages, 9 figures. Submitted to AAS Journals

The K2 Galactic Archaeology Program Data Release 3: Age-abundance patterns in C1-C8 and C10-C18

Authors: Joel C. Zinn, Dennis Stello, Yvonne Elsworth, Rafael A. García, Thomas Kallinger, Savita Mathur, Benoît Mosser, Marc Hon, Lisa Bugnet, Caitlin Jones, Claudia Reyes, Sanjib Sharma, Ralph Schönrich, Jack T. Warfield, Rodrigo Luger, Andrew Vanderburg, Chiaki Kobayashi, Marc H. Pinsonneault, Jennifer A. Johnson, Daniel Huber, Sven Buder, Meridith Joyce, Joss Bland-Hawthorn, Luca Casagrande, Geraint F. Lewis , et al. (6 additional authors not shown)

Abstract: We present the third and final data release of the K2 Galactic Archaeology Program (K2 GAP) for Campaigns C1-C8 and C10-C18. We provide asteroseismic radius and mass coefficients, $κ_R$ and $κ_M$, for $\sim 19,000$ red giant stars, which translate directly to radius and mass given a temperature. As such, K2 GAP DR3 represents the largest asteroseismic sample in the literature to date. K2 GAP DR3 s… ▽ More We present the third and final data release of the K2 Galactic Archaeology Program (K2 GAP) for Campaigns C1-C8 and C10-C18. We provide asteroseismic radius and mass coefficients, $κ_R$ and $κ_M$, for $\sim 19,000$ red giant stars, which translate directly to radius and mass given a temperature. As such, K2 GAP DR3 represents the largest asteroseismic sample in the literature to date. K2 GAP DR3 stellar parameters are calibrated to be on an absolute parallactic scale based on Gaia DR2, with red giant branch and red clump evolutionary state classifications provided via a machine-learning approach. Combining these stellar parameters with GALAH DR3 spectroscopy, we determine asteroseismic ages with precisions of $\sim 20-30\%$ and compare age-abundance relations to Galactic chemical evolution models among both low- and high-$α$ populations for $α$, light, iron-peak, and neutron-capture elements. We confirm recent indications in the literature of both increased Ba production at late Galactic times, as well as significant contribution to r-process enrichment from prompt sources associated with, e.g., core-collapse supernovae. With an eye toward other Galactic archaeology applications, we characterize K2 GAP DR3 uncertainties and completeness using injection tests, suggesting K2 GAP DR3 is largely unbiased in mass/age and with uncertainties of $2.9\%\,(\rm{stat.})\,\pm0.1\%\,(\rm{syst.})$ & $6.7\%\,(\rm{stat.})\,\pm0.3\%\,(\rm{syst.})$ in $κ_R$ & $κ_M$ for red giant branch stars and $4.7\%\,(\rm{stat.})\,\pm0.3\%\,(\rm{syst.})$ & $11\%\,(\rm{stat.})\,\pm0.9\%\,(\rm{syst.})$ for red clump stars. We also identify percent-level asteroseismic systematics, which are likely related to the time baseline of the underlying data, and which therefore should be considered in TESS asteroseismic analysis. △ Less

Submitted 10 March, 2022; v1 submitted 11 August, 2021; originally announced August 2021.

Comments: Published in ApJ; machine-readable tables available through the open access online journal article

Journal ref: ApJ 926 191 (2022)

H-alpha and Ca II Infrared Triplet Variations During a Transit of the 23 Myr Planet V1298 Tau c

Authors: Adina D. Feinstein, Benjamin T. Montet, Marshall C. Johnson, Jacob L. Bean, Trevor J. David, Michael A. Gully-Santiago, John H. Livingston, Rodrigo Luger

Abstract: Young transiting exoplanets (< 100 Myr) provide crucial insight into atmospheric evolution via photoevaporation. However, transmission spectroscopy measurements to determine atmospheric composition and mass loss are challenging due to the activity and prominent stellar disk inhomogeneities present on young stars. We observed a full transit of V1298 Tau c, a 23 Myr, 5.59$R_\oplus$ planet orbiting a… ▽ More Young transiting exoplanets (< 100 Myr) provide crucial insight into atmospheric evolution via photoevaporation. However, transmission spectroscopy measurements to determine atmospheric composition and mass loss are challenging due to the activity and prominent stellar disk inhomogeneities present on young stars. We observed a full transit of V1298 Tau c, a 23 Myr, 5.59$R_\oplus$ planet orbiting a young K0-K1.5 solar analogue with GRACES on Gemini-North. We were able to measure the Doppler tomographic signal of V1298 Tau c using the Ca II infrared triplet (IRT) and find a projected obliquity of $λ= 5^\circ \pm 15^\circ$. The tomographic signal is only seen in the chromospherically driven core of the Ca II IRT, which may be the result of star-planet interactions. Additionally, we find that excess absorption of the H-alpha line decreases smoothly during the transit. While this could be a tentative detection of hot gas escaping the planet, we find this variation is consistent with similar timescale observations of other young stars that lack transiting planets over similar timescales. We show this variation can also be explained by the presence of starspots with surrounding facular regions. More observations both in- and out-of the transits of V1298 Tau c are required to determine the nature of the Ca II IRT and H-alpha line variations. △ Less

Submitted 2 November, 2021; v1 submitted 2 July, 2021; originally announced July 2021.

Comments: 2021, AJ, 162, 213; 17 pages, 11 figures

Linearized Field Deblending: PSF Photometry for Impatient Astronomers

Authors: Christina Hedges, Rodrigo Luger, Jorge Martinez Palomera, Jessie Dotson, Geert Barentsen

Abstract: NASA's Kepler, K2 and TESS missions employ Simple Aperture Photometry (SAP) to derive time-series photometry, where an aperture is estimated for each star, and pixels containing each star are summed to create a single light curve. This method is simple, but in crowded fields the derived time-series can be highly contaminated. The alternate method of fitting a Point Spread Function (PSF) to the dat… ▽ More NASA's Kepler, K2 and TESS missions employ Simple Aperture Photometry (SAP) to derive time-series photometry, where an aperture is estimated for each star, and pixels containing each star are summed to create a single light curve. This method is simple, but in crowded fields the derived time-series can be highly contaminated. The alternate method of fitting a Point Spread Function (PSF) to the data is able to account for crowding, but is computationally expensive. In this paper, we present a new approach to extracting photometry from these time-series missions, which fits the PSF directly, but makes simplifying assumptions in order to greatly reduce the computation expense. Our method fixes the scene of the field in each image, estimates the PSF shape of the instrument with a linear model, and allows only source flux and position to vary. We demonstrate that our method is able to separate the photometry from blended targets in the Kepler dataset that are separated by less than a pixel. Our method is fast to compute, and fully accounts for uncertainties from degeneracies due to crowded fields. We name the method described in this work Linearized Field Deblending (LFD). We demonstrate our method on the false positive Kepler target \koi. We are able to separate the photometry of the two sources in the data, and demonstrate the contaminating transiting signal is consistent with a small, sub-stellar companion with a radius of $2.67R_{jup}$ ($0.27R_{sol}$). Our method is equally applicable to extracting photometry from NASA's TESS mission. △ Less

Submitted 15 June, 2021; originally announced June 2021.

Comments: 22 Pages, 9 Figures

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.

Comments: Published in the Journal of Open Source Software. Comments (still) welcome. Software available at https://docs.exoplanet.codes

Analytic Light Curves in Reflected Light: Phase Curves, Occultations, and Non-Lambertian Scattering for Spherical Planets and Moons

Authors: Rodrigo Luger, Eric Agol, Fran Bartolić, Daniel Foreman-Mackey

Abstract: We derive efficient, closed form, differentiable, and numerically stable solutions for the flux measured from a spherical planet or moon seen in reflected light, either in or out of occultation. Our expressions apply to the computation of scattered light phase curves of exoplanets, secondary eclipse light curves in the optical, or future measurements of planet-moon and planet-planet occultations,… ▽ More We derive efficient, closed form, differentiable, and numerically stable solutions for the flux measured from a spherical planet or moon seen in reflected light, either in or out of occultation. Our expressions apply to the computation of scattered light phase curves of exoplanets, secondary eclipse light curves in the optical, or future measurements of planet-moon and planet-planet occultations, as well as to photometry of solar system bodies. We derive our solutions for Lambertian bodies illuminated by a point source, but extend them to model illumination sources of finite angular size and rough surfaces with phase-dependent scattering. Our algorithm is implemented in Python within the open-source starry mapping framework and is designed with efficient gradient-based inference in mind. The algorithm is 4-5 orders of magnitude faster than direct numerical evaluation methods and about 10 orders of magnitude more precise. We show how the techniques developed here may one day lead to the construction of two-dimensional maps of terrestrial planet surfaces, potentially enabling the detection of continents and oceans on exoplanets in the habitable zone. △ Less

Submitted 10 March, 2021; originally announced March 2021.

Comments: 61 pages, 17 figures. To be submitted to AAS journals

Occultation mapping of Io's surface in the near-infrared I: Inferring static maps

Authors: Fran Bartolić, Rodrigo Luger, Daniel Foreman-Mackey, Robert R. Howell, Julie A. Rathbun

Abstract: Jupiter's moon Io is the most volcanically active body in the Solar System with hundreds of active volcanoes varying in intensity on different timescales. Io has been observed during occultations by other Galilean moons and Jupiter since the 1980s, using high-cadence near infrared photometry. These observations encode a wealth of information about the volcanic features on its surface. We built a g… ▽ More Jupiter's moon Io is the most volcanically active body in the Solar System with hundreds of active volcanoes varying in intensity on different timescales. Io has been observed during occultations by other Galilean moons and Jupiter since the 1980s, using high-cadence near infrared photometry. These observations encode a wealth of information about the volcanic features on its surface. We built a generative model for the observed occultations using the code starry which enables fast, analytic, and differentiable computation of occultation light curves in emitted and reflected light. Our probabilistic Bayesian model is able to recover known hotspots on the surface of Io using only two light curves and without any assumptions on the locations, shapes or the number of spots. The methods we have developed are also directly applicable to the problem of mapping the surfaces of stars and exoplanets. △ Less

Submitted 5 March, 2021; originally announced March 2021.

Comments: 38 pages, 14 figures. To be submitted to AAS Journals

starry_process: Interpretable Gaussian processes for stellar light curves

Authors: Rodrigo Luger, Daniel Foreman-Mackey, Christina Hedges

Abstract: In this note we present the starry_process code, which implements an interpretable Gaussian process (GP) for modeling variability in stellar light curves. As dark starspots rotate in and out of view, the total flux received from a distant star will change over time. Unresolved flux time series therefore encode information about the spatial structure of features on the stellar surface. The starry_p… ▽ More In this note we present the starry_process code, which implements an interpretable Gaussian process (GP) for modeling variability in stellar light curves. As dark starspots rotate in and out of view, the total flux received from a distant star will change over time. Unresolved flux time series therefore encode information about the spatial structure of features on the stellar surface. The starry_process software package allows one to easily model the flux variability due to starspots, whether one is interested in understanding the properties of these spots or marginalizing over the stellar variability when it is treated as a nuisance signal. The main difference between the GP implemented here and typical GPs used to model stellar variability is the explicit dependence of our GP on physical properties of the star, such as its period, inclination, and limb darkening coefficients, and on properties of the spots, such as their radius and latitude distributions. This code is the Python implementation of the interpretable GP algorithm developed in Luger, Foreman-Mackey, and Hedges (2021). △ Less

Submitted 9 February, 2021; v1 submitted 2 February, 2021; originally announced February 2021.

Comments: 4 pages, 3 figures. To be submitted to JOSS. Updated references

Mapping stellar surfaces II: An interpretable Gaussian process model for light curves

Authors: Rodrigo Luger, Daniel Foreman-Mackey, Christina Hedges

Abstract: The use of Gaussian processes (GPs) as models for astronomical time series datasets has recently become almost ubiquitous, given their ease of use and flexibility. GPs excel in particular at marginalization over the stellar signal in cases where the variability due to starspots rotating in and out of view is treated as a nuisance, such as in exoplanet transit modeling. However, these effective mod… ▽ More The use of Gaussian processes (GPs) as models for astronomical time series datasets has recently become almost ubiquitous, given their ease of use and flexibility. GPs excel in particular at marginalization over the stellar signal in cases where the variability due to starspots rotating in and out of view is treated as a nuisance, such as in exoplanet transit modeling. However, these effective models are less useful in cases where the starspot signal is of primary interest since it is not obvious how the parameters of the GP model are related to the physical properties of interest, such as the size, contrast, and latitudinal distribution of the spots. Instead, it is common practice to explicitly model the effect of individual starspots on the light curve and attempt to infer their properties via optimization or posterior inference. Unfortunately, this process is degenerate, ill-posed, and often computationally intractable when applied to stars with more than a few spots and/or to ensembles of many light curves. In this paper, we derive a closed-form expression for the mean and covariance of a Gaussian process model that describes the light curve of a rotating, evolving stellar surface conditioned on a given distribution of starspot sizes, contrasts, and latitudes. We demonstrate that this model is correctly calibrated, allowing one to robustly infer physical parameters of interest from one or more stellar light curves, including the typical radii and the mean and variance of the latitude distribution of starspots. Our GP has far-ranging implications for understanding the variability and magnetic activity of stars from both light curves and radial velocity (RV) measurements, as well as for robustly modeling correlated noise in both transiting and RV exoplanet searches. Our implementation is efficient, user-friendly, and open source, available as the Python package starry-process. △ Less

Submitted 9 February, 2021; v1 submitted 2 February, 2021; originally announced February 2021.

Comments: 91 pages, 30 figures. Updated references. Interactive visualization at http://starry-process.flatironinstitute.org

Multi-Wavelength Photometry Derived from Monochromatic Kepler Data

Authors: Christina Hedges, Rodrigo Luger, Jessie Dotson, Daniel Foreman-Mackey, Geert Barentsen

Abstract: The Kepler mission has provided a wealth of data, revealing new insights in time-domain astronomy. However, Kepler's single band-pass has limited studies to a single wavelength. In this work we build a data-driven, pixel-level model for the Pixel Response Function (PRF) of Kepler targets, modeling the image data from the spacecraft. Our model is sufficiently flexible to capture known detector effe… ▽ More The Kepler mission has provided a wealth of data, revealing new insights in time-domain astronomy. However, Kepler's single band-pass has limited studies to a single wavelength. In this work we build a data-driven, pixel-level model for the Pixel Response Function (PRF) of Kepler targets, modeling the image data from the spacecraft. Our model is sufficiently flexible to capture known detector effects, such as non-linearity, intra-pixel sensitivity variations, and focus change. In theory, the shape of the Kepler PRF should also be weakly wavelength dependent, due to optical chromatic aberration and wavelength dependent detector response functions. We are able to identify these predicted shape changes to the PRF using the residuals between Kepler data and our model. In this work, we show that these PRF changes correspond to wavelength variability in Kepler targets using a small sample of eclipsing binaries. Using our model, we demonstrate that pixel-level light curves of eclipsing binaries show variable eclipse depths, ellipsoidal modulation and limb darkening. These changes at the pixel level are consistent with multi-wavelength photometry. Our work suggests each pixel in the Kepler data of a single target has a different effective wavelength, ranging from $\approx$ 550-750 $nm$. In this proof of concept, we demonstrate our model, and discuss possible use cases for the wavelength dependent Pixel Response Function of Kepler. These use cases include characterizing variable systems, and vetting exoplanet discoveries at the pixel level. The chromatic PRF of Kepler is due to weak wavelength dependence in the optical systems and detector of the telescope, and similar chromatic PRFs are expected in other similar telescopes, notably the NASA TESS telescope. △ Less

Submitted 29 January, 2021; originally announced February 2021.

Comments: 20 pages, 17 figures, 1 table

Journal ref: Published 2021 January 29, The Astronomical Journal, Volume 161, Number 2

Mapping stellar surfaces I: Degeneracies in the rotational light curve problem

Authors: Rodrigo Luger, Daniel Foreman-Mackey, Christina Hedges, David W. Hogg

Abstract: Thanks to missions like Kepler and TESS, we now have access to tens of thousands of high precision, fast cadence, and long baseline stellar photometric observations. In principle, these light curves encode a vast amount of information about stellar variability and, in particular, about the distribution of starspots and other features on their surfaces. Unfortunately, the problem of inferring stell… ▽ More Thanks to missions like Kepler and TESS, we now have access to tens of thousands of high precision, fast cadence, and long baseline stellar photometric observations. In principle, these light curves encode a vast amount of information about stellar variability and, in particular, about the distribution of starspots and other features on their surfaces. Unfortunately, the problem of inferring stellar surface properties from a rotational light curve is famously ill-posed, as it often does not admit a unique solution. Inference about the number, size, contrast, and location of spots can therefore depend very strongly on the assumptions of the model, the regularization scheme, or the prior. The goal of this paper is twofold: (1) to explore the various degeneracies affecting the stellar light curve "inversion" problem and their effect on what can and cannot be learned from a stellar surface given unresolved photometric measurements; and (2) to motivate ensemble analyses of the light curves of many stars at once as a powerful data-driven alternative to common priors adopted in the literature. We further derive novel results on the dependence of the null space on stellar inclination and limb darkening and show that single-band photometric measurements cannot uniquely constrain quantities like the total spot coverage without the use of strong priors. This is the first in a series of papers devoted to the development of novel algorithms and tools for the analysis of stellar light curves and spectral time series, with the explicit goal of enabling statistically robust inference about their surface properties. △ Less

Submitted 9 February, 2021; v1 submitted 29 January, 2021; originally announced February 2021.

Comments: 33 pages, 12 figures. To be submitted to AAS Journals. Updated references

The K2 Galactic Archaeology Program Data Release 2: Asteroseismic results from campaigns 4, 6, & 7

Authors: Joel C. Zinn, Dennis Stello, Yvonne Elsworth, Rafael A. García, Thomas Kallinger, Savita Mathur, Benoît Mosser, Lisa Bugnet, Caitlin Jones, Marc Hon, Sanjib Sharma, Ralph Schönrich, Jack T. Warfield, Rodrigo Luger, Marc H. Pinsonneault, Jennifer A. Johnson, Daniel Huber, Victor Silva Aguirre, William J. Chaplin, Guy R. Davies, Andrea Miglio

Abstract: Studies of Galactic structure and evolution have benefitted enormously from Gaia kinematic information, though additional, intrinsic stellar parameters like age are required to best constrain Galactic models. Asteroseismology is the most precise method of providing such information for field star populations $\textit{en masse}$, but existing samples for the most part have been limited to a few nar… ▽ More Studies of Galactic structure and evolution have benefitted enormously from Gaia kinematic information, though additional, intrinsic stellar parameters like age are required to best constrain Galactic models. Asteroseismology is the most precise method of providing such information for field star populations $\textit{en masse}$, but existing samples for the most part have been limited to a few narrow fields of view by the CoRoT and Kepler missions. In an effort to provide well-characterized stellar parameters across a wide range in Galactic position, we present the second data release of red giant asteroseismic parameters for the K2 Galactic Archaeology Program (GAP). We provide $ν_{\mathrm{max}}$ and $Δν$ based on six independent pipeline analyses; first-ascent red giant branch (RGB) and red clump (RC) evolutionary state classifications from machine learning; and ready-to-use radius & mass coefficients, $κ_R$ & $κ_M$, which, when appropriately multiplied by a solar-scaled effective temperature factor, yield physical stellar radii and masses. In total, we report 4395 radius and mass coefficients, with typical uncertainties of $3.3\% \mathrm{\ (stat.)} \pm 1\% \mathrm{\ (syst.)}$ for $κ_R$ and $7.7\% \mathrm{\ (stat.)} \pm 2\% \mathrm{\ (syst.)}$ for $κ_M$ among RGB stars, and $5.0\% \mathrm{\ (stat.)} \pm 1\% \mathrm{\ (syst.)}$ for $κ_R$ and $10.5\% \mathrm{\ (stat.)} \pm 2\% \mathrm{\ (syst.)}$ for $κ_M$ among RC stars. We verify that the sample is nearly complete -- except for a dearth of stars with $ν_{\mathrm{max}} \lesssim 10-20μ$Hz -- by comparing to Galactic models and visual inspection. Our asteroseismic radii agree with radii derived from Gaia Data Release 2 parallaxes to within $2.2 \pm 0.3\%$ for RGB stars and $2.0 \pm 0.6\%$ for RC stars. △ Less

Submitted 7 December, 2020; originally announced December 2020.

Comments: Published in ApJS

Journal ref: The Astrophysical Journal Supplement Series, Volume 251, Issue 2, 23 pp. (2020)

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.

Comments: 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/

Quantifying the Stellar Halo's Response to the LMC's Infall with Spherical Harmonics

Authors: Emily C. Cunningham, Nicolas Garavito-Camargo, Alis J. Deason, Kathryn V. Johnston, Denis Erkal, Chervin F. P. Laporte, Gurtina Besla, Rodrigo Luger, Robyn E. Sanderson

Abstract: The vast majority of the mass in the Milky Way (MW) is in dark matter (DM); we therefore cannot directly observe the MW mass distribution, and have to use tracer populations in order to infer properties of the MW DM halo. However, MW halo tracers do not only feel the gravitational influence of the MW itself. Tracers can also be affected by MW satellites; Garavito-Camargo et al. (2019) (hereafter G… ▽ More The vast majority of the mass in the Milky Way (MW) is in dark matter (DM); we therefore cannot directly observe the MW mass distribution, and have to use tracer populations in order to infer properties of the MW DM halo. However, MW halo tracers do not only feel the gravitational influence of the MW itself. Tracers can also be affected by MW satellites; Garavito-Camargo et al. (2019) (hereafter GC19) demonstrate that the Large Magellanic Cloud (LMC) induces a density wake in the MW DM, resulting in large scale kinematic patterns in the MW stellar halo. In this work, we use spherical harmonic expansion (SHE) of the velocity fields of simulated stellar halos in an effort to disentangle perturbations on large scales (e.g., due to the LMC itself as well as the LMC-induced DM wake) and small scales (due to substructure). Using the GC19 simulations, we demonstrate how the different terms in the SHE of the stellar velocity field reflect the different wake components, and show that these signatures are a strong function of the LMC mass. An exploration of model halos built from accreted dwarfs Bullock & Johnston (2005) suggests that stellar debris from massive, recent accretion events can produce much more power in the velocity angular power spectra than the perturbation from the LMC-induced wake. We therefore consider two models for the Sagittarius (Sgr) stream -- the most recent, massive accretion event in the MW apart from the LMC -- and find that the angular power on large scales is generally dominated by the LMC-induced wake, even when Sgr is included. We conclude that SHE of the MW stellar halo velocity field may therefore be a useful tool in quantifying the response of the MW DM halo to the LMC's infall. △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: 25 pages, 13 figures. Accepted for publication in ApJ

The Young Planet DS Tuc Ab has a Low Obliquity

Authors: Benjamin T. Montet, Adina D. Feinstein, Rodrigo Luger, Megan E. Bedell, Michael A. Gully-Santiago, Johanna K. Teske, Sharon Xuesong Wang, R. Paul Butler, Erin Flowers, Stephen A. Shectman, Jeffrey D. Crane, Ian B. Thompson

Abstract: The abundance of short-period planetary systems with high orbital obliquities relative to the spin of their host stars is often taken as evidence that scattering processes play important roles in the formation and evolution of these systems. More recent studies have suggested that wide binary companions can tilt protoplanetary disks, inducing a high stellar obliquity that form through smooth proce… ▽ More The abundance of short-period planetary systems with high orbital obliquities relative to the spin of their host stars is often taken as evidence that scattering processes play important roles in the formation and evolution of these systems. More recent studies have suggested that wide binary companions can tilt protoplanetary disks, inducing a high stellar obliquity that form through smooth processes like disk migration. DS Tuc Ab, a transiting planet with an 8.138 day period in the 40 Myr Tucana-Horologium association, likely orbits in the same plane as its now-dissipated protoplanetary disk, enabling us to test these theories of disk physics. Here, we report on Rossiter-McLaughlin observations of one transit of DS Tuc Ab with the Planet Finder Spectrograph on the Magellan Clay Telescope at Las Campanas Observatory. We confirm the previously detected planet by modeling the planet transit and stellar activity signals simultaneously. We test multiple models to describe the stellar activity-induced radial velocity variations over the night of the transit, finding the obliquity to be low: $λ= 12 \pm 13$ degrees, suggesting that this planet likely formed through smooth disk processes and its protoplanetary disk was not significantly torqued by DS Tuc B. The specific stellar activity model chosen affects the results at the $\approx 5$ degree level. This is the youngest planet to be observed using this technique; we provide a discussion on best practices to accurately measure the observed signal of similar young planets. △ Less

Submitted 17 January, 2020; v1 submitted 8 December, 2019; originally announced December 2019.

Comments: 18 pages, four figures, AJ accepted

Four newborn planets transiting the young solar analog V1298 Tau

Authors: Trevor J. David, Erik A. Petigura, Rodrigo Luger, Daniel Foreman-Mackey, John H. Livingston, Eric E. Mamajek, Lynne A. Hillenbrand

Abstract: Exoplanets orbiting pre-main sequence stars are laboratories for studying planet evolution processes, including atmospheric loss, orbital migration, and radiative cooling. V1298 Tau, a young solar analog with an age of 23 $\pm$ 4 Myr, is one such laboratory. The star is already known to host a Jupiter-sized planet on a 24 day orbit. Here, we report the discovery of three additional planets --- all… ▽ More Exoplanets orbiting pre-main sequence stars are laboratories for studying planet evolution processes, including atmospheric loss, orbital migration, and radiative cooling. V1298 Tau, a young solar analog with an age of 23 $\pm$ 4 Myr, is one such laboratory. The star is already known to host a Jupiter-sized planet on a 24 day orbit. Here, we report the discovery of three additional planets --- all between the size of Neptune and Saturn --- based on our analysis of K2 Campaign 4 photometry. Planets c and d have sizes of 5.6 and 6.4 $R_\oplus$, respectively and with orbital periods of 8.25 and 12.40 days reside 0.25% outside of the nominal 3:2 mean-motion resonance. Planet e is 8.7 $R_\oplus$ in size but only transited once in the K2 time series and thus has a period longer than 36 days, but likely shorter than 223 days. The V1298 Tau system may be a precursor to the compact multiplanet systems found to be common by the Kepler mission. However, the large planet sizes stand in sharp contrast to the vast majority of Kepler multis which have planets smaller than 3 $R_\oplus$. Simple dynamical arguments suggest total masses of $<$28 $M_\oplus$ and $<$120 $M_\oplus$ for the c-d and d-b planet pairs, respectively. The implied low masses suggest that the planets may still be radiatively cooling and contracting, and perhaps losing atmosphere. The V1298 Tau system offers rich prospects for further follow-up including atmospheric characterization by transmission or eclipse spectroscopy, dynamical characterization through transit-timing variations, and measurements of planet masses and obliquities by radial velocities. △ Less

Submitted 10 October, 2019; originally announced October 2019.

Comments: accepted to ApJL. 18 pages, 4 figures, 1 table

Analytic Planetary Transit Light Curves and Derivatives for Stars with Polynomial Limb Darkening

Authors: Eric Agol, Rodrigo Luger, Daniel Foreman-Mackey

Abstract: We derive analytic, closed-form solutions for the light curve of a planet transiting a star with a limb darkening profile which is a polynomial function of the stellar elevation, up to arbitrary integer order. We provide improved analytic expressions for the uniform, linear, and quadratic limb-darkened cases, as well as novel expressions for higher order integer powers of limb darkening. The formu… ▽ More We derive analytic, closed-form solutions for the light curve of a planet transiting a star with a limb darkening profile which is a polynomial function of the stellar elevation, up to arbitrary integer order. We provide improved analytic expressions for the uniform, linear, and quadratic limb-darkened cases, as well as novel expressions for higher order integer powers of limb darkening. The formulae are crafted to be numerically stable over the expected range of usage. We additionally present analytic formulae for the partial derivatives of instantaneous flux with respect to the radius ratio, impact parameter, and limb darkening coefficients. These expressions are rapid to evaluate, and compare quite favorably in speed and accuracy to existing transit light curve codes. We also use these expressions to numerically compute the first partial derivatives of exposure-time averaged transit light curves with respect to all model parameters. An additional application is modeling eclipsing binary or eclipsing multiple star systems in cases where the stars may be treated as spherically symmetric. We provide code which implements these formulae in C++, Python, IDL, and Julia, with tests and examples of usage. △ Less

Submitted 11 October, 2019; v1 submitted 8 August, 2019; originally announced August 2019.

Comments: 57 pages, 18 figures. Accepted to AJ

Detection of Hundreds of New Planet Candidates and Eclipsing Binaries in K2 Campaigns 0-8

Authors: Ethan Kruse, Eric Agol, Rodrigo Luger, Daniel Foreman-Mackey

Abstract: We implement a search for exoplanets in campaigns zero through eight (C0-8) of the K2 extension of the Kepler spacecraft. We apply a modified version of the QATS planet search algorithm to K2 light curves produced by the EVEREST pipeline, carrying out the C0-8 search on $1.5 \times 10^5$ target stars with magnitudes in the range of Kp = 9-15. We detect 818 transiting planet candidates, of which 37… ▽ More We implement a search for exoplanets in campaigns zero through eight (C0-8) of the K2 extension of the Kepler spacecraft. We apply a modified version of the QATS planet search algorithm to K2 light curves produced by the EVEREST pipeline, carrying out the C0-8 search on $1.5 \times 10^5$ target stars with magnitudes in the range of Kp = 9-15. We detect 818 transiting planet candidates, of which 374 were undiscovered by prior searches, with {64,15,5,2,1} in {2,3,4,5,6}-planet multi-planet candidate systems, respectively. Of the new planets detected, 100 orbit M dwarfs, including one which is potentially rocky and in the habitable zone. 154 of our candidates reciprocally transit with our Solar System: they are geometrically aligned to see at least one Solar System planet transit. We find candidates which display transit timing variations and dozens of candidates on both period extremes with single transits or ultra-short periods. We point to evidence that our candidates display similar patterns in frequency and size-period relation as confirmed planets, such as tentative evidence for the radius gap. Confirmation of these planet candidates with follow-up studies will increase the number of K2 planets by up to 50%, and characterization of their host stars will improve statistical studies of planet properties. Our sample includes many planets orbiting bright stars amenable for radial velocity follow-up and future characterization with JWST. We also list the 579 eclipsing binary systems detected as part of this search. △ Less

Submitted 13 September, 2019; v1 submitted 24 July, 2019; originally announced July 2019.

Comments: Updated to match published version in ApJS. 37 pages, 32 figures, two very long tables. Machine readable tables available in source

EXOFASTv2: A public, generalized, publication-quality exoplanet modeling code

Authors: Jason D. Eastman, Joseph E. Rodriguez, Eric Agol, Keivan G. Stassun, Thomas G. Beatty, Andrew Vanderburg, B. Scott Gaudi, Karen A. Collins, Rodrigo Luger

Abstract: We present the next generation public exoplanet fitting software, EXOFASTv2. It is capable of fitting an arbitrary number of planets, radial velocity data sets, astrometric data sets, and/or transits observed with any combination of wavelengths. We model the star simultaneously in the fit and provide several state-of-the-art ways to constrain its properties, including taking advantage of the now-u… ▽ More We present the next generation public exoplanet fitting software, EXOFASTv2. It is capable of fitting an arbitrary number of planets, radial velocity data sets, astrometric data sets, and/or transits observed with any combination of wavelengths. We model the star simultaneously in the fit and provide several state-of-the-art ways to constrain its properties, including taking advantage of the now-ubiquitous all-sky catalog photometry and Gaia parallaxes. EXOFASTv2 can model the star by itself, too. Multi-planet systems are modeled self-consistently with the same underlying stellar mass that defines their semi-major axes through Kepler's law and the planetary period. Transit timing, duration, and depth variations can be modeled with a simple command line option. We explain our methodology and rationale as well as provide an improved version of the core transit model that is both 25\% faster and more accurate. We highlight several potential pitfalls in exoplanet modeling, including the handling of eccentricity in transit-only fits, that the standard exoplanet convention for $ω$ uses a left-handed coordinate system, contrary to most modern textbooks, how to avoid an important degeneracy when allowing negative companion masses, and a widely unappreciated, potential 10-minute ambiguity in the reported transit times. EXOFASTv2 is available at https://github.com/jdeast/EXOFASTv2 . The code is written in IDL, and includes an executable that can be run freely and legally without an IDL license or any knowledge of the language. Extensive documentation and tutorials are included in the distribution for a variety of example fits. Advanced amateurs and undergrads have successfully performed sophisticated global fits of complex planetary systems with EXOFASTv2. It is therefore a powerful tool for education and outreach as well as the broader professional community. △ Less

Submitted 22 July, 2019; originally announced July 2019.

Comments: 51 pages, 14 figures, 3 Tables. Submitted to PASP. Comments welcome

Detecting volcanically produced tori along orbits of exoplanets using UV spectroscopy

Authors: Kristina G. Kislyakova, Luca Fossati, Denis Shulyak, Eike Günther, Manuel Güdel, Colin P. Johnstone, Vladimir Airapetian, Sudeshna Boro Saikia, Allan Sacha Brun, Vera Dobos, Kevin France, Eric Gaidos, Maxim L. Khodachenko, Antonino F. Lanza, Helmut Lammer, Lena Noack, Rodrigo Luger, Antoine Strugarek, Aline Vidotto, Allison Youngblood

Abstract: We suggest to use the Hubble Space Telescople (HST) follow-up observations of the TESS targets for detecting possible plasma tori along the orbits of exoplanets orbiting M dwarfs. The source of the torus could be planetary volcanic activity due to tidal or electromagnetic induction heating. Fast losses to space for planets orbiting these active stars can lead to the lost material forming a torus a… ▽ More We suggest to use the Hubble Space Telescople (HST) follow-up observations of the TESS targets for detecting possible plasma tori along the orbits of exoplanets orbiting M dwarfs. The source of the torus could be planetary volcanic activity due to tidal or electromagnetic induction heating. Fast losses to space for planets orbiting these active stars can lead to the lost material forming a torus along the planetary orbit, similar to the Io plasma torus. We show that such torus would be potentially detectable by the HST in the UV. △ Less

Submitted 11 July, 2019; originally announced July 2019.

Comments: White paper submitted on the use of HST for TESS follow-up observations, 5 pages, 2 figures

On The XUV Luminosity Evolution of TRAPPIST-1

Authors: David P. Fleming, Rory Barnes, Rodrigo Luger, Jacob T. VanderPlas

Abstract: We model the long-term XUV luminosity of TRAPPIST-1 to constrain the evolving high-energy radiation environment experienced by its planetary system. Using Markov Chain Monte Carlo (MCMC), we derive probabilistic constraints for TRAPPIST-1's stellar and XUV evolution that account for observational uncertainties, degeneracies between model parameters, and empirical data of low-mass stars. We constra… ▽ More We model the long-term XUV luminosity of TRAPPIST-1 to constrain the evolving high-energy radiation environment experienced by its planetary system. Using Markov Chain Monte Carlo (MCMC), we derive probabilistic constraints for TRAPPIST-1's stellar and XUV evolution that account for observational uncertainties, degeneracies between model parameters, and empirical data of low-mass stars. We constrain TRAPPIST-1's mass to $m_{\star} = 0.089 \pm{0.001}$ M$_{\odot}$ and find that its early XUV luminosity likely saturated at $\log_{10}(L_{XUV}/L_{bol}) = -3.03^{+0.23}_{-0.12}$. From the posterior distribution, we infer that there is a ${\sim}40\%$ chance that TRAPPIST-1 is still in the saturated phase today, suggesting that TRAPPIST-1 has maintained high activity and $L_{XUV}/L_{bol} \approx 10^{-3}$ for several Gyrs. TRAPPIST-1's planetary system therefore likely experienced a persistent and extreme XUV flux environment, potentially driving significant atmospheric erosion and volatile loss. The inner planets likely received XUV fluxes ${\sim}10^3 - 10^4\times$ that of the modern Earth during TRAPPIST-1's ${\sim}1$ Gyr-long pre-main sequence phase. Deriving these constraints via MCMC is computationally non-trivial, so scaling our methods to constrain the XUV evolution of a larger number of M dwarfs that harbor terrestrial exoplanets would incur significant computational expenses. We demonstrate that approxposterior, an open source Python machine learning package for approximate Bayesian inference using Gaussian processes, accurately and efficiently replicates our analysis using $980\times$ less computational time and $1330\times$ fewer simulations than MCMC sampling using emcee. We find that approxposterior derives constraints with mean errors on the best fit values and $1σ$ uncertainties of $0.61\%$ and $5.5\%$, respectively, relative to emcee. △ Less

Submitted 17 February, 2020; v1 submitted 12 June, 2019; originally announced June 2019.

Comments: Accepted in ApJ, 16 pages, 7 figures

VPLanet: The Virtual Planet Simulator

Authors: Rory Barnes, Rodrigo Luger, Russell Deitrick, Peter Driscoll, Thomas R. Quinn, David P. Fleming, Hayden Smotherman, Diego V. McDonald, Caitlyn Wilhelm, Rodolfo Garcia, Patrick Barth, Benjamin Guyer, Victoria S. Meadows, Cecilia M. Bitz, Pramod Gupta, Shawn D. Domagal-Goldman, John Armstrong

Abstract: We describe a software package called VPLanet that simulates fundamental aspects of planetary system evolution over Gyr timescales, with a focus on investigating habitable worlds. In this initial release, eleven physics modules are included that model internal, atmospheric, rotational, orbital, stellar, and galactic processes. Many of these modules can be coupled simultaneously to simulate the evo… ▽ More We describe a software package called VPLanet that simulates fundamental aspects of planetary system evolution over Gyr timescales, with a focus on investigating habitable worlds. In this initial release, eleven physics modules are included that model internal, atmospheric, rotational, orbital, stellar, and galactic processes. Many of these modules can be coupled simultaneously to simulate the evolution of terrestrial planets, gaseous planets, and stars. The code is validated by reproducing a selection of observations and past results. VPLanet is written in C and designed so that the user can choose the physics modules to apply to an individual object at runtime without recompiling, i.e., a single executable can simulate the diverse phenomena that are relevant to a wide range of planetary and stellar systems. This feature is enabled by matrices and vectors of function pointers that are dynamically allocated and populated based on user input. The speed and modularity of VPLanet enables large parameter sweeps and the versatility to add/remove physical phenomena to assess their importance. VPLanet is publicly available from a repository that contains extensive documentation, numerous examples, Python scripts for plotting and data management, and infrastructure for community input and future development. △ Less

Submitted 27 August, 2019; v1 submitted 15 May, 2019; originally announced May 2019.

Comments: 75 pages, 34 figures, 10 tables, accepted to the Proceedings of the Astronomical Society of the Pacific. Source code, documentation, and examples available at https://github.com/VirtualPlanetaryLaboratory/vplanet

No Evidence for Lunar Transit in New Analysis of Hubble Space Telescope Observations of the Kepler-1625 System

Authors: Laura Kreidberg, Rodrigo Luger, Megan Bedell

Abstract: Observations of the Kepler-1625 system with the Kepler and Hubble Space Telescopes have suggested the presence of a candidate exomoon, Kepler-1625b I, a Neptune-radius satellite orbiting a long-period Jovian planet. Here we present a new analysis of the Hubble observations, using an independent data reduction pipeline. We find that the transit light curve is well fit with a planet-only model, with… ▽ More Observations of the Kepler-1625 system with the Kepler and Hubble Space Telescopes have suggested the presence of a candidate exomoon, Kepler-1625b I, a Neptune-radius satellite orbiting a long-period Jovian planet. Here we present a new analysis of the Hubble observations, using an independent data reduction pipeline. We find that the transit light curve is well fit with a planet-only model, with a best-fit $χ^2_ν$ equal to 1.01. The addition of a moon does not significantly improve the fit quality. We compare our results directly with the original light curve from Teachey & Kipping (2018), and find that we obtain a better fit to the data using a model with fewer free parameters (no moon). We discuss possible sources for the discrepancy in our results, and conclude that the lunar transit signal found by Teachey & Kipping (2018) was likely an artifact of the data reduction. This finding highlights the need to develop independent pipelines to confirm results that push the limits of measurement precision. △ Less

Submitted 21 May, 2019; v1 submitted 23 April, 2019; originally announced April 2019.

Comments: 7 pages, 5 figures, accepted to ApJL

TESS Photometric Mapping of a Terrestrial Planet in the Habitable Zone: Detection of Clouds, Oceans, and Continents

Authors: Rodrigo Luger, Megan Bedell, Roland Vanderspek, Christopher J. Burke

Abstract: To date, a handful of exoplanets have been photometrically mapped using phase-modulated reflection or emission from their surfaces, but the small amplitudes of such signals have limited previous maps almost exclusively to coarse dipolar features on hot giant planets. In this work, we uncover a signal using recently released data from the Transiting Exoplanet Survey Satellite (TESS), which we show… ▽ More To date, a handful of exoplanets have been photometrically mapped using phase-modulated reflection or emission from their surfaces, but the small amplitudes of such signals have limited previous maps almost exclusively to coarse dipolar features on hot giant planets. In this work, we uncover a signal using recently released data from the Transiting Exoplanet Survey Satellite (TESS), which we show corresponds to time-variable reflection from a terrestrial planet with a rotation period of 0.9972696 days. Using a spherical harmonic-based reflection model developed as an extension of the STARRY package, we are able to reconstruct the surface features of this rocky world. We recover a time-variable albedo map of the planet including persistent regions which we interpret as oceans and cloud banks indicative of continental features. We argue that this planet represents the most promising detection of a habitable world to date, although the potential intelligence of any life on it is yet to be determined. △ Less

Submitted 28 March, 2019; originally announced March 2019.

Comments: 22 pages, 9 figures. To be submitted to AAS journals after slight modifications

The Pointing Limits of Transiting Exoplanet Light Curve Characterization with Pixel Level De-correlation

Authors: Nicholas Saunders, Rodrigo Luger, Rory Barnes

Abstract: We present scope (Simulated CCD Observations for Photometric Experimentation), a Python package to create a forward model of telescope detectors and simulate stellar targets with motion relative to the CCD. The primary application of this package is the simulation of the Kepler Space Telescope detector to predict and characterize increased instrumental noise in the spacecraft's final campaigns of… ▽ More We present scope (Simulated CCD Observations for Photometric Experimentation), a Python package to create a forward model of telescope detectors and simulate stellar targets with motion relative to the CCD. The primary application of this package is the simulation of the Kepler Space Telescope detector to predict and characterize increased instrumental noise in the spacecraft's final campaigns of observation. As the fuel powering the spacecraft's stabilizing thrusters ran out and thruster fires began to sputter and fail, stellar Point Spread Functions (PSFs) experienced more extreme and less predictable motion relative to regions of varied sensitivity on the spacecraft detector, generating more noise in transiting exoplanet light curves. Using our simulations, we demonstrate that current de-trending techniques effectively capture and remove systematics caused by sensitivity variation for spacecraft motion as high as about ten times that typically experienced by K2. The scope package is open-source and has been generalized to allow custom detector and stellar target parameters. Future applications include simulating observations made by the Transiting Exoplanet Survey Satellite (TESS) and ground based observations with synthetic atmospheric interference as testbeds for noise-removal techniques. △ Less

Submitted 22 March, 2019; originally announced March 2019.

Comments: 11 pages, 4 figures

eleanor: An open-source tool for extracting light curves from the TESS Full-Frame Images

Authors: Adina D. Feinstein, Benjamin T. Montet, Daniel Foreman-Mackey, Megan E. Bedell, Nicholas Saunders, Jacob L. Bean, Jessie L. Christiansen, Christina Hedges, Rodrigo Luger, Daniel Scolnic, Jose Vinicius de Miranda Cardoso

Abstract: During its two year prime mission the Transiting Exoplanet Survey Satellite (TESS) will perform a time-series photometric survey covering over 80% of the sky. This survey comprises observations of 26 24 x 96 degree sectors that are each monitored continuously for approximately 27 days. The main goal of TESS is to find transiting planets around 200,000 pre-selected stars for which fixed aperture ph… ▽ More During its two year prime mission the Transiting Exoplanet Survey Satellite (TESS) will perform a time-series photometric survey covering over 80% of the sky. This survey comprises observations of 26 24 x 96 degree sectors that are each monitored continuously for approximately 27 days. The main goal of TESS is to find transiting planets around 200,000 pre-selected stars for which fixed aperture photometry is recorded every two minutes. However, TESS is also recording and delivering Full-Frame Images (FFIs) of each detector at a 30 minute cadence. We have created an open-source tool, eleanor, to produce light curves for objects in the TESS FFIs. Here, we describe the methods used in eleanor to produce light curves that are optimized for planet searches. The tool performs background subtraction, aperture and PSF photometry, decorrelation of instrument systematics, and cotrending using principal component analysis. We recover known transiting exoplanets in the FFIs to validate the pipeline and perform a limited search for new planet candidates in Sector 1. Our tests indicate that eleanor produces light curves with significantly less scatter than other tools that have been used in the literature. Cadence-stacked images, and raw and detrended eleanor light curves for each analyzed star will be hosted on MAST, with planet candidates on ExoFOP-TESS as Community TESS Objects of Interest (CTOIs). This work confirms the promise that the TESS FFIs will enable the detection of thousands of new exoplanets and a broad range of time domain astrophysics. △ Less

Submitted 13 June, 2019; v1 submitted 21 March, 2019; originally announced March 2019.

Comments: 21 pages, 13 figures, 2 tables, Accepted to PASP

Rotation Period Evolution in Low-Mass Binary Stars: The Impact of Tidal Torques and Magnetic Braking

Authors: David P. Fleming, Rory Barnes, James R. A. Davenport, Rodrigo Luger

Abstract: We examine how tides, stellar evolution, and magnetic braking shape the rotation period (P$_{rot}$) evolution of low-mass stellar binaries up to orbital periods (P$_{orb}$) of 100 d across a wide range tidal dissipation parameters using two common equilibrium tidal models. We find that many binaries with P$_{orb} < 20$ d tidally lock, and most with $P_{orb} < 4$ d tidally lock into synchronous rot… ▽ More We examine how tides, stellar evolution, and magnetic braking shape the rotation period (P$_{rot}$) evolution of low-mass stellar binaries up to orbital periods (P$_{orb}$) of 100 d across a wide range tidal dissipation parameters using two common equilibrium tidal models. We find that many binaries with P$_{orb} < 20$ d tidally lock, and most with $P_{orb} < 4$ d tidally lock into synchronous rotation on circularized orbits. At short P$_{orb}$, tidal torques produce a population of fast rotators that single-star only models of magnetic braking fail to produce. In many cases, we show that the competition between magnetic braking and tides produces a population of subsynchronous rotators that persists for Gyrs, even in short P$_{orb}$ binaries, qualitatively reproducing the subsynchronous eclipsing binaries (EBs) discovered in the Kepler field by Lurie et al. (2017). Both equilibrium tidal models predict that binaries can tidally-interact out to P$_{orb} \approx 80$ d, while the Constant Phase Lag tidal model predicts that binaries can tidally lock out to P$_{orb} \approx 100$ d. Tidal torques often force the P$_{rot}$ evolution of stellar binaries to depart from the long-term magnetic braking-driven spin down experienced by single stars, revealing that P$_{rot}$ is not be a valid proxy for age in all cases, i.e. gyrochronology can underpredict ages by up to $300\%$ unless one accounts for binarity. We suggest that accurate determinations of orbital eccentricties and P$_{rot}$ can be used to discriminate between which equilibrium tidal models best describes tidal interactions in low-mass binary stars. △ Less

Submitted 1 July, 2019; v1 submitted 13 March, 2019; originally announced March 2019.

Comments: Accepted, ApJ

Detecting Ocean Glint on Exoplanets Using Multiphase Mapping

Authors: Jacob Lustig-Yaeger, Victoria S. Meadows, Guadalupe Tovar Mendoza, Edward W. Schwieterman, Yuka Fujii, Rodrigo Luger, Tyler D. Robinson

Abstract: Rotational mapping and glint are two proposed methods to directly detect liquid water on the surface of habitable exoplanets. However, false positives for both methods may prevent the unambiguous detection of exoplanet oceans. We use simulations of Earth as an exoplanet to introduce a combination of multiwavelength, multiphase, time-series direct-imaging observations and accompanying analyses that… ▽ More Rotational mapping and glint are two proposed methods to directly detect liquid water on the surface of habitable exoplanets. However, false positives for both methods may prevent the unambiguous detection of exoplanet oceans. We use simulations of Earth as an exoplanet to introduce a combination of multiwavelength, multiphase, time-series direct-imaging observations and accompanying analyses that may improve the robustness of exoplanet ocean detection by spatially mapping the ocean glint signal. As the planet rotates, the glint spot appears to "blink" as Lambertian scattering continents interrupt the specular reflection from the ocean. This manifests itself as a strong source of periodic variability in crescent-phase reflected light curves. We invert these light curves to constrain the longitudinal slice maps and apparent albedo of two surfaces at both quadrature and crescent phase. At crescent phase, the retrieved apparent albedo of ocean-bearing longitudinal slices is increased by a factor of 5, compared to the albedo at quadrature phase, due to the contribution from glint. The land-bearing slices exhibit no significant change in apparent albedo with phase. The presence of forward-scattering clouds in our simulated observation increases the overall reflectivity toward crescent, but clouds do not correlate with any specific surfaces, thereby allowing for the phase-dependent glint effect to be interpreted as distinct from cloud scattering. Retrieving the same longitudinal map at quadrature and crescent phases may be used to tie changes in the apparent albedo with phase back to specific geographic surfaces, although this requires ideal geometries. We estimate that crescent-phase time-dependent glint detections are feasible for between 1-10 habitable zone exoplanets orbiting the nearest G, K, and M dwarfs using a space-based, high-contrast, direct-imaging telescope with a diameter >6 m. △ Less

Submitted 15 January, 2019; originally announced January 2019.

Comments: 21 pages, 9 figures, published in AJ

Journal ref: AJ, 156, 301 (2018)

Wobble: A Data-driven Analysis Technique for Time-series Stellar Spectra

Authors: Megan Bedell, David W. Hogg, Daniel Foreman-Mackey, Benjamin T. Montet, Rodrigo Luger

Abstract: In recent years, dedicated extreme-precision radial velocity (EPRV) spectrographs have produced vast quantities of high-resolution, high-signal-to-noise time-series spectra for bright stars. These data contain valuable information for the dual purposes of planet detection via the measured RVs and stellar characterization via the co-added spectra. However, considerable data analysis challenges exis… ▽ More In recent years, dedicated extreme-precision radial velocity (EPRV) spectrographs have produced vast quantities of high-resolution, high-signal-to-noise time-series spectra for bright stars. These data contain valuable information for the dual purposes of planet detection via the measured RVs and stellar characterization via the co-added spectra. However, considerable data analysis challenges exist in extracting these data products from the observed spectra at the highest possible precision, including the issue of poorly-characterized telluric absorption features and the common use of an assumed stellar spectral template. In both of these examples, precision-limiting reliance on external information can be sidestepped using the data directly. Here we propose a data-driven method to simultaneously extract precise RVs and infer the underlying stellar and telluric spectra using a linear model (in the log of flux). The model employs a convex objective and convex regularization to keep the optimization of the spectral components fast. We implement this method in wobble, an open-source python package which uses TensorFlow in one of its first non-neural-network applications to astronomical data. In this work, we demonstrate the performance of wobble on archival HARPS spectra. We recover the canonical exoplanet 51 Pegasi b, detect the secular RV evolution of the M dwarf Barnard's Star, and retrieve the Rossiter-McLaughlin effect for the Hot Jupiter HD 189733b. The method additionally produces extremely high-S/N composite stellar spectra and detailed time-variable telluric spectra, which we also present here. △ Less

Submitted 30 September, 2019; v1 submitted 2 January, 2019; originally announced January 2019.

Comments: closer match to published version

Journal ref: AJ, 158, 164 (2019)

STARRY: Analytic Occultation Light Curves

Authors: Rodrigo Luger, Eric Agol, Daniel Foreman-Mackey, David P. Fleming, Jacob Lustig-Yaeger, Russell Deitrick

Abstract: We derive analytic, closed form, numerically stable solutions for the total flux received from a spherical planet, moon or star during an occultation if the specific intensity map of the body is expressed as a sum of spherical harmonics. Our expressions are valid to arbitrary degree and may be computed recursively for speed. The formalism we develop here applies to the computation of stellar trans… ▽ More We derive analytic, closed form, numerically stable solutions for the total flux received from a spherical planet, moon or star during an occultation if the specific intensity map of the body is expressed as a sum of spherical harmonics. Our expressions are valid to arbitrary degree and may be computed recursively for speed. The formalism we develop here applies to the computation of stellar transit light curves, planetary secondary eclipse light curves, and planet-planet/planet-moon occultation light curves, as well as thermal (rotational) phase curves. In this paper we also introduce STARRY, an open-source package written in C++ and wrapped in Python that computes these light curves. The algorithm in STARRY is six orders of magnitude faster than direct numerical integration and several orders of magnitude more precise. STARRY also computes analytic derivatives of the light curves with respect to all input parameters for use in gradient-based optimization and inference, such as Hamiltonian Monte Carlo (HMC), allowing users to quickly and efficiently fit observed light curves to infer properties of a celestial body's surface map. △ Less

Submitted 15 October, 2018; originally announced October 2018.

Comments: 55 pages, 20 figures. Accepted to the Astronomical Journal. Check out the code at https://github.com/rodluger/starry

Evolved Climates and Observational Discriminants for the TRAPPIST-1 Planetary System

Authors: Andrew P. Lincowski, Victoria S. Meadows, David Crisp, Tyler D. Robinson, Rodrigo Luger, Jacob Lustig-Yaeger, Giada N. Arney

Abstract: The TRAPPIST-1 planetary system provides an unprecedented opportunity to study terrestrial exoplanet evolution with the James Webb Space Telescope (JWST) and ground-based observatories. Since M dwarf planets likely experience extreme volatile loss, the TRAPPIST-1 planets may have highly-evolved, possibly uninhabitable atmospheres. We used a versatile, 1D terrestrial-planet climate model with line-… ▽ More The TRAPPIST-1 planetary system provides an unprecedented opportunity to study terrestrial exoplanet evolution with the James Webb Space Telescope (JWST) and ground-based observatories. Since M dwarf planets likely experience extreme volatile loss, the TRAPPIST-1 planets may have highly-evolved, possibly uninhabitable atmospheres. We used a versatile, 1D terrestrial-planet climate model with line-by-line radiative transfer and mixing length convection (VPL Climate) coupled to a terrestrial photochemistry model to simulate environmental states for the TRAPPIST-1 planets. We present equilibrium climates with self-consistent atmospheric compositions, and observational discriminants of post-runaway, desiccated, 10-100 bar O2- and CO2-dominated atmospheres, including interior outgassing, as well as for water-rich compositions. Our simulations show a range of surface temperatures, most of which are not habitable, although an aqua-planet TRAPPIST-1 e could maintain a temperate surface given Earth-like geological outgassing and CO2. We find that a desiccated TRAPPIST-1 h may produce habitable surface temperatures beyond the maximum greenhouse distance. Potential observational discriminants for these atmospheres in transmission and emission spectra are influenced by photochemical processes and aerosol formation, and include collision-induced oxygen absorption (O2-O2), and O3, CO, SO2, H2O, and CH4 absorption features, with transit signals of up to 200 ppm. Our simulated transmission spectra are consistent with K2, HST, and Spitzer observations of the TRAPPIST-1 planets. For several terrestrial atmospheric compositions, we find that TRAPPIST-1 b is unlikely to produce aerosols. These results can inform JWST observation planning and data interpretation for the TRAPPIST-1 system and other M dwarf terrestrial planets. △ Less

Submitted 20 September, 2018; originally announced September 2018.

Comments: 44 pages, 19 figures, accepted for publication in the Astrophysical Journal

On The Lack of Circumbinary Planets Orbiting Isolated Binary Stars

Authors: David P. Fleming, Rory Barnes, David E. Graham, Rodrigo Luger, Thomas R. Quinn

Abstract: We outline a mechanism that explains the observed lack of circumbinary planets (CBPs) via coupled stellar-tidal evolution of isolated binary stars. Tidal forces between low-mass, short-period binary stars on the pre-main sequence slow the stellar rotations, transferring rotational angular momentum to the orbit as the stars approach the tidally locked state. This transfer increases the binary orbit… ▽ More We outline a mechanism that explains the observed lack of circumbinary planets (CBPs) via coupled stellar-tidal evolution of isolated binary stars. Tidal forces between low-mass, short-period binary stars on the pre-main sequence slow the stellar rotations, transferring rotational angular momentum to the orbit as the stars approach the tidally locked state. This transfer increases the binary orbital period, expanding the region of dynamical instability around the binary, and destabilizing CBPs that tend to preferentially orbit just beyond the initial dynamical stability limit. After the stars tidally lock, we find that angular momentum loss due to magnetic braking can significantly shrink the binary orbit, and hence the region of dynamical stability, over time impacting where surviving CBPs are observed relative to the boundary. We perform simulations over a wide range of parameter space and find that the expansion of the instability region occurs for most plausible initial conditions and that in some cases, the stability semi-major axis doubles from its initial value. We examine the dynamical and observable consequences of a CBP falling within the dynamical instability limit by running N-body simulations of circumbinary planetary systems and find that typically, at least one planet is ejected from the system. We apply our theory to the shortest period {\it Kepler} binary that possesses a CBP, Kepler-47, and find that its existence is consistent with our model. Under conservative assumptions, we find that coupled stellar-tidal evolution of pre-main sequence binary stars removes at least one close-in CBP in $87\%$ of multi-planet circumbinary systems. △ Less

Submitted 10 April, 2018; originally announced April 2018.

Comments: 26 pages, 16 figures, accepted, ApJ

Planet-Planet Occultations in TRAPPIST-1 and Other Exoplanet Systems

Authors: Rodrigo Luger, Jacob Lustig-Yaeger, Eric Agol

Abstract: We explore the occurrence and detectability of planet-planet occultations (PPOs) in exoplanet systems. These are events during which a planet occults the disk of another planet in the same system, imparting a small photometric signal as its thermal or reflected light is blocked. We focus on the planets in TRAPPIST-1, whose orbital planes we show are aligned to within 0.3 degrees at 90% confidence.… ▽ More We explore the occurrence and detectability of planet-planet occultations (PPOs) in exoplanet systems. These are events during which a planet occults the disk of another planet in the same system, imparting a small photometric signal as its thermal or reflected light is blocked. We focus on the planets in TRAPPIST-1, whose orbital planes we show are aligned to within 0.3 degrees at 90% confidence. We present a photodynamical model for predicting and computing PPOs in TRAPPIST-1 and other systems for various assumptions of the planets' atmospheric states. When marginalizing over the uncertainties on all orbital parameters, we find that the rate of PPOs in TRAPPIST-1 is about 1.4 per day. We investigate the prospects for detection of these events with the James Webb Space Telescope, finding that ~10-20 occultations per year of b and c should be above the noise level at 12-15 microns. Joint modeling of several of these PPOs could lead to a robust detection. Alternatively, observations with the proposed Origins Space Telescope should be able to detect individual PPOs at high signal-to-noise. We show how PPOs can be used to break transit timing variation degeneracies, imposing strong constraints on the eccentricities and masses of the planets, as well as to constrain the longitudes of nodes and thus the complete three-dimensional structure of the system. We further show how modeling of these events can be used to reveal a planet's day/night temperature contrast and construct crude surface maps. We make our photodynamical code available on github. △ Less

Submitted 15 November, 2017; originally announced November 2017.

Comments: 36 pages, 25 figures. Accepted to ApJ. Multi-purpose photodynamical code available at github.com/rodluger/planetplanet

Linear models for systematics and nuisances

Authors: Rodrigo Luger, Daniel Foreman-Mackey, David W. Hogg

Abstract: The target of many astronomical studies is the recovery of tiny astrophysical signals living in a sea of uninteresting (but usually dominant) noise. In many contexts (i.e., stellar time-series, or high-contrast imaging, or stellar spectroscopy), there are structured components in this noise caused by systematic effects in the astronomical source, the atmosphere, the telescope, or the detector. Mor… ▽ More The target of many astronomical studies is the recovery of tiny astrophysical signals living in a sea of uninteresting (but usually dominant) noise. In many contexts (i.e., stellar time-series, or high-contrast imaging, or stellar spectroscopy), there are structured components in this noise caused by systematic effects in the astronomical source, the atmosphere, the telescope, or the detector. More often than not, evaluation of the true physical model for these nuisances is computationally intractable and dependent on too many (unknown) parameters to allow rigorous probabilistic inference. Sometimes, housekeeping data---and often the science data themselves---can be used as predictors of the systematic noise. Linear combinations of simple functions of these predictors are often used as computationally tractable models that can capture the nuisances. These models can be used to fit and subtract systematics prior to investigation of the signals of interest, or they can be used in a simultaneous fit of the systematics and the signals. In this Note, we show that if a Gaussian prior is placed on the weights of the linear components, the weights can be marginalized out with an operation in pure linear algebra, which can (often) be made fast. We illustrate this model by demonstrating the applicability of a linear model for the non-linear systematics in K2 time-series data, where the dominant noise source for many stars is spacecraft motion and variability. △ Less

Submitted 30 October, 2017; originally announced October 2017.

Comments: 5 pages, 1 figure. Accepted to Research Notes of the AAS

A seven-planet resonant chain in TRAPPIST-1

Authors: Rodrigo Luger, Marko Sestovic, Ethan Kruse, Simon L. Grimm, Brice-Olivier Demory, Eric Agol, Emeline Bolmont, Daniel Fabrycky, Catarina S. Fernandes, Valérie Van Grootel, Adam Burgasser, Michaël Gillon, James G. Ingalls, Emmanuël Jehin, Sean N. Raymond, Franck Selsis, Amaury H. M. J. Triaud, Thomas Barclay, Geert Barentsen, Steve B. Howell, Laetitia Delrez, Julien de Wit, Daniel Foreman-Mackey, Daniel L. Holdsworth, Jérémy Leconte , et al. (8 additional authors not shown)

Abstract: The TRAPPIST-1 system is the first transiting planet system found orbiting an ultra-cool dwarf star. At least seven planets similar to Earth in radius and in mass were previously found to transit this host star. Subsequently, TRAPPIST-1 was observed as part of the K2 mission and, with these new data, we report the measurement of an 18.77 d orbital period for the outermost planet, TRAPPIST-1h, whic… ▽ More The TRAPPIST-1 system is the first transiting planet system found orbiting an ultra-cool dwarf star. At least seven planets similar to Earth in radius and in mass were previously found to transit this host star. Subsequently, TRAPPIST-1 was observed as part of the K2 mission and, with these new data, we report the measurement of an 18.77 d orbital period for the outermost planet, TRAPPIST-1h, which was unconstrained until now. This value matches our theoretical expectations based on Laplace relations and places TRAPPIST-1h as the seventh member of a complex chain, with three-body resonances linking every member. We find that TRAPPIST-1h has a radius of 0.727 Earth radii and an equilibrium temperature of 173 K. We have also measured the rotational period of the star at 3.3 d and detected a number of flares consistent with a low-activity, middle-aged, late M dwarf. △ Less

Submitted 22 May, 2017; v1 submitted 12 March, 2017; originally announced March 2017.

Comments: 42 pages, 8 figures, 2 tables

Journal ref: Nature Astronomy 1, 0129 (2017)

An update to the EVEREST K2 pipeline: Short cadence, saturated stars, and Kepler-like photometry down to Kp = 15

Authors: Rodrigo Luger, Ethan Kruse, Daniel Foreman-Mackey, Eric Agol, Nicholas Saunders

Abstract: We present an update to the EVEREST K2 pipeline that addresses various limitations in the previous version and improves the photometric precision of the de-trended light curves. We develop a fast regularization scheme for third order pixel level decorrelation (PLD) and adapt the algorithm to include the PLD vectors of neighboring stars to enhance the predictive power of the model and minimize over… ▽ More We present an update to the EVEREST K2 pipeline that addresses various limitations in the previous version and improves the photometric precision of the de-trended light curves. We develop a fast regularization scheme for third order pixel level decorrelation (PLD) and adapt the algorithm to include the PLD vectors of neighboring stars to enhance the predictive power of the model and minimize overfitting, particularly for faint stars. We also modify PLD to work for saturated stars and improve its performance on extremely variable stars. On average, EVEREST 2.0 light curves have 10-20% higher photometric precision than those in the previous version, yielding the highest precision light curves at all Kp magnitudes of any publicly available K2 catalog. For most K2 campaigns, we recover the original Kepler precision to at least Kp = 14, and to at least Kp = 15 for campaigns 1, 5, and 6. We also de-trend all short cadence targets observed by K2, obtaining even higher photometric precision for these stars. All light curves for campaigns 0-8 are available online in the EVEREST catalog, which will be continuously updated with future campaigns. EVEREST 2.0 is open source and is coded in a general framework that can be applied to other photometric surveys, including Kepler and the upcoming TESS mission. △ Less

Submitted 17 February, 2017; originally announced February 2017.

Comments: 16 pages, 21 figures. Submitted to AJ. Source code and documentation available at https://github.com/rodluger/everest

The Pale Green Dot: A Method to Characterize Proxima Centauri b using Exo-Aurorae

Authors: Rodrigo Luger, Jacob Lustig-Yaeger, David P. Fleming, Matt A. Tilley, Eric Agol, Victoria S. Meadows, Russell Deitrick, Rory Barnes

Abstract: We examine the feasibility of detecting auroral emission from the potentially habitable exoplanet Proxima Centauri b. Detection of aurorae would yield an independent confirmation of the planet's existence, constrain the presence and composition of its atmosphere, and determine the planet's eccentricity and inclination, thereby breaking the mass-inclination degeneracy. If Proxima Centauri b is a te… ▽ More We examine the feasibility of detecting auroral emission from the potentially habitable exoplanet Proxima Centauri b. Detection of aurorae would yield an independent confirmation of the planet's existence, constrain the presence and composition of its atmosphere, and determine the planet's eccentricity and inclination, thereby breaking the mass-inclination degeneracy. If Proxima Centauri b is a terrestrial world with an Earth-like atmosphere and magnetic field, we estimate the power at the 5577Å OI auroral line is on the order of 0.1 TW under steady-state stellar wind, or ${\sim} 100 {\times}$ stronger than that on Earth. This corresponds to a planet-star contrast ratio of $10^{-6}-10^{-7}$ in a narrow band about the 5577Å line, although higher contrast ($10^{-4}-10^{-5}$) may be possible during periods of strong magnetospheric disturbance (auroral power $1-10$ TW). We searched the Proxima Centauri b HARPS data for the 5577Å line and for other prominent oxygen and nitrogen lines, but find no signal, indicating that the OI auroral line contrast must be lower than $2\times 10^{-2}$ (with power $\lesssim$ 3,000 TW), consistent with our predictions. We find that observations of 0.1 TW auroral emission lines are likely infeasible with current and planned telescopes. However, future observations with a space-based coronagraphic telescope or a ground-based extremely large telescope (ELT) with a coronagraph could push sensitivity down to terawatt oxygen aurorae (contrast $7\times 10^{-6}$) with exposure times of ${\sim} 1$ day. If a coronagraph design contrast of $10^{-7}$ can be achieved with negligible instrumental noise, a future concept ELT could observe steady-state auroral emission in a few nights. △ Less

Submitted 14 February, 2017; v1 submitted 28 September, 2016; originally announced September 2016.

Comments: 19 pages, 4 tables, 10 figures. Accepted to ApJ

The Habitability of Proxima Centauri b: II: Environmental States and Observational Discriminants

Authors: Victoria S. Meadows, Giada N. Arney, Edward W. Schwieterman, Jacob Lustig-Yaeger, Andrew P. Lincowski, Tyler Robinson, Shawn D. Domagal-Goldman, Rory K. Barnes, David P. Fleming, Russell Deitrick, Rodrigo Luger, Peter E. Driscoll, Thomas R. Quinn, David Crisp

Abstract: Proxima Centauri b provides an unprecedented opportunity to understand the evolution and nature of terrestrial planets orbiting M dwarfs. Although Proxima Cen b orbits within its star's habitable zone, multiple plausible evolutionary paths could have generated different environments that may or may not be habitable. Here we use 1D coupled climate-photochemical models to generate self-consistent at… ▽ More Proxima Centauri b provides an unprecedented opportunity to understand the evolution and nature of terrestrial planets orbiting M dwarfs. Although Proxima Cen b orbits within its star's habitable zone, multiple plausible evolutionary paths could have generated different environments that may or may not be habitable. Here we use 1D coupled climate-photochemical models to generate self-consistent atmospheres for evolutionary scenarios predicted in our companion paper (Barnes et al., 2016). These include high-O2, high-CO2, and more Earth-like atmospheres, with either oxidizing or reducing compositions. We show that these modeled environments can be habitable or uninhabitable at Proxima Cen b's position in the habitable zone. We use radiative transfer models to generate synthetic spectra and thermal phase curves for these simulated environments, and instrument models to explore our ability to discriminate between possible planetary states. These results are applicable not only to Proxima Cen b, but to other terrestrial planets orbiting M dwarfs. Thermal phase curves may provide the first constraint on the existence of an atmosphere, and JWST observations longward of 7 microns could characterize atmospheric heat transport and molecular composition. Detection of ocean glint is unlikely with JWST, but may be within the reach of larger aperture telescopes. Direct imaging spectra may detect O4, which is diagnostic of massive water loss and O2 retention, rather than a photosynthesis. Similarly, strong CO2 and CO bands at wavelengths shortward of 2.5 μm would indicate a CO2-dominated atmosphere. If the planet is habitable and volatile-rich, direct imaging will be the best means of detecting habitability. Earth-like planets with microbial biospheres may be identified by the presence of CH4 and either photosynthetically produced O2 or a hydrocarbon haze layer. △ Less

Submitted 30 August, 2016; originally announced August 2016.

Comments: 93 pages, 36 figures, 3 tables. Submitted to Astrobiology

The Habitability of Proxima Centauri b I: Evolutionary Scenarios

Authors: Rory Barnes, Russell Deitrick, Rodrigo Luger, Peter E. Driscoll, Thomas R. Quinn, David P. Fleming, Benjamin Guyer, Diego V. McDonald, Victoria S. Meadows, Giada Arney, David Crisp, Shawn D. Domagal-Goldman, Daniel Foreman-Mackey, Nathan A. Kaib, Andrew Lincowski, Jacob Lustig-Yaeger, Eddie Schwieterman

Abstract: We analyze the evolution of the potentially habitable planet Proxima Centauri b to identify environmental factors that affect its long-term habitability. We consider physical processes acting on size scales ranging from the galactic to the stellar system to the planet's core. We find that there is a significant probability that Proxima Centauri has had encounters with its companion stars, Alpha Ce… ▽ More We analyze the evolution of the potentially habitable planet Proxima Centauri b to identify environmental factors that affect its long-term habitability. We consider physical processes acting on size scales ranging from the galactic to the stellar system to the planet's core. We find that there is a significant probability that Proxima Centauri has had encounters with its companion stars, Alpha Centauri A and B, that are close enough to destabilize an extended planetary system. If the system has an additional planet, as suggested by the discovery data, then it may perturb planet b's eccentricity and inclination, possibly driving those parameters to non-zero values, even in the presence of strong tidal damping. We also model the internal evolution of the planet, evaluating the roles of different radiogenic abundances and tidal heating and find that magnetic field generation is likely for billions of years. We find that if planet b formed in situ, then it experienced 169 +/- 13 million years in a runaway greenhouse as the star contracted during its formation. This early phase could remove up to 5 times as much water as in the modern Earth's oceans, possibly producing a large abiotic oxygen atmosphere. On the other hand, if Proxima Centauri b formed with a substantial hydrogen atmosphere (0.01 - 1% of the planet's mass), then this envelope could have shielded the water long enough for it to be retained before being blown off itself. After modeling this wide range of processes we conclude that water retention during the host star's pre-main sequence phase is the biggest obstacle for Proxima b's habitability. These results are all obtained with a new software package called VPLANET. △ Less

Submitted 5 March, 2018; v1 submitted 24 August, 2016; originally announced August 2016.

Comments: 46 pages, 25 figures. Final version, includes corrected geophysical evolution and expanded analysis orbital/rotational evolution and atmosphere loss

EVEREST: Pixel Level Decorrelation of K2 Light curves

Authors: Rodrigo Luger, Eric Agol, Ethan Kruse, Rory Barnes, Andrew Becker, Daniel Foreman-Mackey, Drake Deming

Abstract: We present EVEREST, an open-source pipeline for removing instrumental noise from K2 light curves. EVEREST employs a variant of pixel level decorrelation (PLD) to remove systematics introduced by the spacecraft's pointing error and a Gaussian process (GP) to capture astrophysical variability. We apply EVEREST to all K2 targets in campaigns 0-7, yielding light curves with precision comparable to tha… ▽ More We present EVEREST, an open-source pipeline for removing instrumental noise from K2 light curves. EVEREST employs a variant of pixel level decorrelation (PLD) to remove systematics introduced by the spacecraft's pointing error and a Gaussian process (GP) to capture astrophysical variability. We apply EVEREST to all K2 targets in campaigns 0-7, yielding light curves with precision comparable to that of the original Kepler mission for stars brighter than $K_p \approx 13$, and within a factor of two of the Kepler precision for fainter targets. We perform cross-validation and transit injection and recovery tests to validate the pipeline, and compare our light curves to the other de-trended light curves available for download at the MAST High Level Science Products archive. We find that EVEREST achieves the highest average precision of any of these pipelines for unsaturated K2 stars. The improved precision of these light curves will aid in exoplanet detection and characterization, investigations of stellar variability, asteroseismology, and other photometric studies. The EVEREST pipeline can also easily be applied to future surveys, such as the TESS mission, to correct for instrumental systematics and enable the detection of low signal-to-noise transiting exoplanets. The EVEREST light curves and the source code used to generate them are freely available online. △ Less

Submitted 2 July, 2016; originally announced July 2016.

Comments: 15 pages, 17 figures. Revised version submitted to AJ