skip to main content


Title: Analytic Light Curves in Reflected Light: Phase Curves, Occultations, and Non-Lambertian Scattering for Spherical Planets and Moons
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, 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 ∼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. 6 6 https://github.com/rodluger/starrynight  more » « less
Award ID(s):
1907342
NSF-PAR ID:
10353632
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
The Astronomical Journal
Volume:
164
Issue:
1
ISSN:
0004-6256
Page Range / eLocation ID:
4
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    We present simulated optical light curves of super-Eddington tidal disruption events (TDEs) using the ZEro-BeRnoulli Accretion (ZEBRA) flow model, which proposes that during the super-Eddington phase, the disc is quasi-spherical, radiation-pressure dominated, and accompanied by the production of strong jets. We construct light curves for both on- and off-axis (with respect to the jet) observers to account for the anisotropic nature of the jetted emission. We find that at optical wavelengths, emission from the accretion flow is orders of magnitude brighter than that produced by the jet, even with boosting from synchrotron self-Compton. Comparing to the observed jetted TDE Swift J2058.4+0516, we find that the ZEBRA model accurately captures the time-scale for which accretion remains super-Eddington and reproduces the luminosity of the transient. However, we find the shape of the light curves deviate at early times and the radius and temperature of our modelled ZEBRA are ∼2.7–4.1 times smaller and ∼1.4–2.3 times larger, respectively, than observed. We suggest that this indicates the ZEBRA inflates more, and more rapidly, than currently predicted by the model, and we discuss possible extensions to the model to account for this. Such refinements, coupled with valuable new data from upcoming large-scale surveys, could help to resolve the nature of super-Eddington TDEs and how they are powered.

     
    more » « less
  2. Abstract During the Zwicky Transient Facility (ZTF) Phase I operations, 78 hydrogen-poor superluminous supernovae (SLSNe-I) were discovered in less than 3 yr, constituting the largest sample from a single survey. This paper (Paper I) presents the data, including the optical/UV light curves and classification spectra, while Paper II in this series will focus on the detailed analysis of the light curves and modeling. Our photometry is primarily taken by ZTF in the g , r , and i bands, and with additional data from other ground-based facilities and Swift. The events of our sample cover a redshift range of z = 0.06 − 0.67, with a median and 1 σ error (16% and 84% percentiles) of z med = 0.265 − 0.135 + 0.143 . The peak luminosity covers −22.8 mag ≤ M g ,peak ≤ −19.8 mag, with a median value of − 21.48 − 0.61 + 1.13 mag. The light curves evolve slowly with a mean rest-frame rise time of t rise = 41.9 ± 17.8 days. The luminosity and timescale distributions suggest that low-luminosity SLSNe-I with a peak luminosity ∼−20 mag or extremely fast-rising events (<10 days) exist, but are rare. We confirm previous findings that slowly rising SLSNe-I also tend to fade slowly. The rest-frame color and temperature evolution show large scatters, suggesting that the SLSN-I population may have diverse spectral energy distributions. The peak rest-frame color shows a moderate correlation with the peak absolute magnitude, i.e., brighter SLSNe-I tend to have bluer colors. With optical and UV photometry, we construct the bolometric luminosity and derive a bolometric correction relation that is generally applicable for converting g , r -band photometry to the bolometric luminosity for SLSNe-I. 
    more » « less
  3. Abstract

    Nearly one million light curves from the TESS Year 1 southern hemisphere extracted from Full Field Images with the DIAmante pipeline are processed through the AutoRegressive Planet Search statistical procedure. ARIMA models remove lingering autocorrelated noise, the Transit Comb Filter identifies the strongest periodic signal in the light curve, and a Random Forest machine-learning classifier is trained and applied to identify the best potential candidates. Classifier training sets are based on injections of planetary transit signals, eclipsing binaries, and other variable stars. The optimized classifier has a True Positive Rate of 92.5% and a False Positive Rate of 0.43% from the labeled training set. The result of this DIAmante TESS autoregressive planet search of the southern ecliptic hemisphere analysis is a list of 7377 potential exoplanet candidates. The classifier had a 64% recall rate for previously confirmed exoplanets and a 78% negative recall rate for known False Positives. The completeness map of the injected planetary signals shows high recall rates for planets with 8–30Rradii and periods 0.6–13 days and poor completeness for planets with radii <2Ror periods <1 day. The list has many False Alarms and False Positives that need to be culled with multifaceted vetting operations (Paper II).

     
    more » « less
  4. ABSTRACT

    Eclipsing binaries are important benchmark objects to test and calibrate stellar structure and evolution models. This is especially true for binaries with a fully convective M-dwarf component for which direct measurements of these stars’ masses and radii are difficult using other techniques. Within the potential of M-dwarfs to be exoplanet host stars, the accuracy of theoretical predictions of their radius and effective temperature as a function of their mass is an active topic of discussion. Not only the parameters of transiting exoplanets but also the success of future atmospheric characterization relies on accurate theoretical predictions. We present the analysis of five eclipsing binaries with low-mass stellar companions out of a subsample of 23, for which we obtained ultra-high-precision light curves using the CHEOPS satellite. The observation of their primary and secondary eclipses are combined with spectroscopic measurements to precisely model the primary parameters and derive the M-dwarfs mass, radius, surface gravity, and effective temperature estimates using the PYCHEOPS data analysis software. Combining these results to the same set of parameters derived from TESS light curves, we find very good agreement (better than 1 per cent for radius and better than 0.2 per cent for surface gravity). We also analyse the importance of precise orbits from radial velocity measurements and find them to be crucial to derive M-dwarf radii in a regime below 5 per cent accuracy. These results add five valuable data points to the mass–radius diagram of fully convective M-dwarfs.

     
    more » « less
  5. ABSTRACT

    A new window to observing individual stars and other small sources at cosmological distances was opened recently, with the detection of several caustic-crossing events in galaxy cluster fields. Many more such events are expected soon from dedicated campaigns with the Hubble Space Telescope and the James Webb Space Telescope. These events can not only teach us about the lensed sources themselves, such as individual high-redshift stars, star clusters, or accretion discs, but through their light curves they also hold information about the point-mass function of the lens, and thus, potentially, the composition of dark matter. We present here a simple method for simulating light curves of such events, i.e. the change in apparent magnitude of the source as it sweeps over the net of caustics generated by microlenses embedded around the critical region of the lens. The method is recursive and so any reasonably sized small source can be accommodated, down to sub-solar scales, in principle. We compare the method, which we dub Adaptive Boundary Method, with other common methods such as simple inverse ray shooting, and demonstrate that it is significantly more efficient and accurate in the small-source and high-magnification regime of interest. A python version of the code is made publicly available in an open-source fashion for simulating future events.

     
    more » « less