Disintegrating Exoplanets: Creating Size Constraints by Statistically Peering Through the Debris

We study two intriguing exoplanets, Kepler-1520b and K2-22b, that are being disintegrated by their host stars, producing dust and debris pulled from their surface into tails that trail and precede the exoplanets in their orbits, making it difficult to discern the true nature of the objects. Our goal is to constrain the radius of the underlying objects, and while previous studies have done this in the past by selecting shallow transit events, we attempt a new statistical approach to model the intrinsic astrophysical and photon noise distributions simultaneously. We assume that the lightcurve flux distribution is distributed as a convolution of a Gaussian photon noise component and a Raleigh astrophysical component. The Raleigh curve has a finite flux maximum, which we fit with a Hamiltonian Markov Chain. With these methods, a more accurate flux maximum may be estimated, producing a better final value for the radius of these exoplanets. To determine statistical significance, we used the python package PyMC3 to find the posterior distribution for our data with Gaussian, Rayleigh, and joint function curves and plotting it against our collected flux. After completing this analysis, we were able to place a constraint on the maximum radius of K2-22b at approximately 15,000km. However, we were unable to constrain the radius of Kepler-1520b due to either transits obscured by forward scattering of dust or a grazing transit.