Search for: All records

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.

 

Stellar variability is a limiting factor for planet detection and characterization, particularly around active M-type stars. Here we revisit one of the most active stars from the Kepler mission, the M4 star GJ 1243, and use a sample of 414 flare events from 11 months of 1-minute cadence light curves to study the empirical morphology of white-light stellar flares. We use a Gaussian process detrending technique to account for the underlying starspots. We present an improved analytic, continuous flare template that is generated by stacking the flares onto a scaled time and amplitude and uses a Markov Chain Monte Carlo analysis to fit the model. Our model is defined using classical flare events but can also be used to model complex, multipeaked flare events. We demonstrate the utility of our model using TESS data at the 10-minute, 2-minute, and 20 s cadence modes. Our new flare model code is made publicly available on GitHub.⁵⁵

https://github.com/lupitatovar/Llamaradas-Estelares

 

Abstract We describe a new transit-detection algorithm designed to detect single-transit events in discontinuous Perkins INfrared Exosatellite Survey (PINES) observations of L and T dwarfs. We use this algorithm to search for transits in 131 PINES light curves and identify two transit candidates: 2MASS J18212815+1414010 (2MASS J1821+1414) and 2MASS J08350622+1953050 (2MASS J0835+1953). We disfavor 2MASS J1821+1414 as a genuine transit candidate due to the known variability properties of the source. We cannot rule out the planetary nature of 2MASS J0835+1953's candidate event and perform follow-up observations in an attempt to recover a second transit. A repeat event has yet to be observed, but these observations suggest that target variability is an unlikely cause of the candidate transit. We perform a Markov Chain Monte Carlo simulation of the light curve and estimate a planet radius ranging from 4.2 − 1.6 + 3.5 R ⊕ to 5.8 − 2.1 + 4.8 R ⊕ , depending on the host’s age. Finally, we perform an injection and recovery simulation on our light-curve sample. We inject planets into our data using measured M-dwarf planet occurrence rates and attempt to recover them using our transit-search algorithm. Our detection rates suggest that, assuming M-dwarf planet occurrence rates, we should have roughly a 1% chance of detecting a candidate that could cause the transit depth we observe for 2MASS J0835+1953. If 2MASS J0835+1953 b is confirmed, it would suggest an enhancement in the occurrence of short-period planets around L and T dwarfs in comparison to M dwarfs, which would challenge predictions from planet formation models. 

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 

Abstract We describe the Perkins INfrared Exosatellite Survey (PINES), a near-infrared photometric search for short-period transiting planets and moons around a sample of 393 spectroscopically confirmed L- and T-type dwarfs. PINES is performed with Boston University’s 1.8 m Perkins Telescope Observatory, located on Anderson Mesa, Arizona. We discuss the observational strategy of the survey, which was designed to optimize the number of expected transit detections, and describe custom automated observing procedures for performing PINES observations. We detail the steps of the PINES Analysis Toolkit ( PAT ), software that is used to create light curves from PINES images. We assess the impact of second-order extinction due to changing precipitable water vapor on our observations and find that the magnitude of this effect is minimized in Mauna Kea Observatories J band. We demonstrate the validity of PAT through the recovery of a transit of WASP-2 b and known variable brown dwarfs, and use it to identify a new variable L/T transition object: the T2 dwarf WISE J045746.08-020719.2. We report on the measured photometric precision of the survey and use it to estimate our transit-detection sensitivity. We find that for our median brightness targets, assuming contributions from white noise only, we are sensitive to the detection of 2.5 R ⊕ planets and larger. PINES will test whether the increase in sub-Neptune-sized planet occurrence with decreasing host mass continues into the L- and T-dwarf regime.