Abstract We present a solution for the light curve of two bodies mutually transiting a star with polynomial limb darkening. The term “mutual transit” in this work refers to a transit of the star during which overlap occurs between the two transiting bodies. These could be an exoplanet with an exomoon companion, two exoplanets, an eclipsing binary and a planet, or two stars eclipsing a third in a triple-star system. We include analytic derivatives of the light curve with respect to the positions and radii of both bodies. We provide code that implements a photodynamical model for a mutual transit. We include two dynamical models, one for hierarchical systems in which a secondary body orbits a larger primary (e.g., an exomoon system) and a second for confocal systems in which two bodies independently orbit a central mass (e.g., two planets in widely separated orbits). Our code is fast enough to enable inference with Markov Chain Monte Carlo algorithms, and the inclusion of derivatives allows for the use of gradient-based inference methods such as Hamiltonian Monte Carlo. While applicable to a variety of systems, this work was undertaken primarily with exomoons in mind. It is our hope that making this code publicly available will reduce barriers for the community to assess the detectability of exomoons, conduct searches for exomoons, and attempt to validate existing exomoon candidates. We also anticipate that our code will be useful for studies of planet–planet transits in exoplanetary systems, transits of circumbinary planets, and eclipses in triple-star systems.
more »
« less
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
- PAR ID:
- 10353632
- 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
-
-
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
-
Abstract Some hydrogen-rich core-collapse supernovae (SNeIIP) exhibit evidence of a sustained energy source powering their light curves, resulting in a brighter and/or longer-lasting hydrogen recombination plateau phase. We present a semi-analytic SNIIP light-curve model that accounts for the effects of an arbitrary internal heating source, considering as special cases56Ni/56Co decay, a central engine (magnetar or accreting compact object), and shock interaction with a dense circumstellar disk. While a sustained internal power source can boost the plateau luminosity commensurate with the magnitude of the power, the duration of the recombination plateau can typically be increased by at most a factor of ∼2–3 compared to the zero-heating case. For a given ejecta mass and initial kinetic energy, the longest plateau duration is achieved for a constant heating rate at the highest magnitude that does not appreciably accelerate the ejecta. This finding has implications for the minimum ejecta mass required to explain particularly long-lasting SNe, such as iPTF14hls, and for confidently identifying rare explosions of the most massive hydrogen-rich (e.g., Population III) stars. We present a number of analytic estimates that elucidate the key features of the detailed model.more » « less
-
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
-
Abstract One of the unique aspects of Earth is that it has a fractionally large Moon, which is thought to have formed from a Moon-forming disk generated by a giant impact. The Moon stabilizes the Earth’s spin axis at least by several degrees and contributes to Earth’s stable climate. Given that impacts are common during planet formation, exomoons, which are moons around planets in extrasolar systems, should be common as well, but no exomoon has been confirmed. Here we propose that an initially vapor-rich moon-forming disk is not capable of forming a moon that is large with respect to the size of the planet because growing moonlets, which are building blocks of a moon, experience strong gas drag and quickly fall toward the planet. Our impact simulations show that terrestrial and icy planets that are larger than ~1.3−1.6 R ⊕ produce entirely vapor disks, which fail to form a fractionally large moon. This indicates that (1) our model supports the Moon-formation models that produce vapor-poor disks and (2) rocky and icy exoplanets whose radii are smaller than ~1.6 R ⊕ are ideal candidates for hosting fractionally large exomoons.more » « less
An official website of the United States government

