LUNA: An algorithm for generating dynamic planet-moon transits

Abstract: It has been previously shown that moons of extrasolar planets may be detectable with the Kepler Mission, for moon masses above ~0.2 Earth masses Kipping et al. 2009c. Transit timing effects have been formerly identified as a potent tool to this end, exploiting the dynamics of the system. In this work, we explore the simulation of transit light curves of a planet plus a single moon including not only the transit timing effects but also the light curve signal of the moon itself. We introduce our new algorithm, LUNA, which produces transit light curves for both bodies, analytically accounting for shadow overlaps, stellar limb darkening and planet-moon dynamical motion. By building the dynamics into the core of LUNA, the routine automatically accounts for transit timing/duration variations and ingress/egress asymmetries for not only the planet, but also the moon. We then generate some artificial data for two feasibly detectable hypothetical systems of interest: a i) prograde and ii) retrograde Earth-like moon around a habitable-zone Neptune for a M-dwarf system. We fit the hypothetical systems using LUNA and demonstrate the feasibility of detecting these cases with Kepler photometry.

• The paper presents LUNA, an analytic algorithm designed for rapidly simulating dynamic transit light curves of planet-moon systems, incorporating key effects like TTV/TDV and limb darkening.
• LUNA's analytic nature allows for rapid, large-scale simulation and parameter space exploration, proving effective even when tested with noisy data comparable to Kepler observations.
• This algorithm provides researchers a powerful tool to search archival data for exomoons, enhancing understanding of planetary system architecture and dynamics.

Overview of the Exomoon Transit Algorithm

David M. Kipping's paper presents a comprehensive analytic algorithm designed to simulate transits of exoplanetary systems with moons, a topic of significant interest in exoplanetary research. The algorithm aims to produce dynamic transit light curves that account for both the planetary and lunar transits, incorporating various physical and observational effects such as transit timing, duration variations, and stellar limb darkening.

Analytical Approach and Methodological Contributions

The proposed algorithm stands out for its analytic nature; it eschews numerical components in favor of closed-form solutions, thereby offering rapid computations suitable for large-scale data analysis in transit searches. This analytic framework is critical, as potential exomoon detections require extensive parameter space exploration, often through computationally demanding techniques like Monte Carlo methods.

In terms of methodology, the paper details the transit generation process for a star-planet-moon system, emphasizing:

• Dynamic Integration: By integrating the dynamics of the planet-moon system, the algorithm inherently accounts for interactions such as positional and velocity wobbles that induce transit timing variations (TTV) and transit duration variations (TDV). This dynamic approach marks an improvement over static approximations, traditionally used in transit modeling.
• Limb Darkening: The algorithm meticulously incorporates stellar limb darkening, supporting non-linear laws which are essential for analyzing high-precision data such as that from the Kepler mission.
• Orbital Elements: It rigorously includes all relevant orbital elements, allowing simulations to reflect real-world complexities such as orbital eccentricity and inclination effects on transit light curves.

Numerical Experiments and Feasibility Studies

The paper provides examples of simulated data for hypothetical exomoons in distinct scenarios, such as a prograde or retrograde Earth-like moon orbiting a Neptune-sized planet within an M-dwarf system's habitable zone. These examples are integral as they validate the algorithm's capacity to retrieve parameters accurately in noisy environments comparable to those experienced with Kepler data.

These simulations highlight that the algorithm can distinguish between different moon orbital configurations based on the photometric signatures in the light curves. Moreover, the paper demonstrates feasibility in detecting such systems with current instrumentation, reinforcing the algorithm's potential utility in exomoon searches.

Implications and Future Directions

The creation of this algorithm has significant theoretical and practical implications for the field of exoplanetary science. Practically, it equips researchers with a robust tool for mining archival data for exomoons, potentially broadening our understanding of such systems' frequency and characteristics. Theoretically, the inclusion of an exomoon detection capability will enhance our understanding of planetary system dynamics and the potential habitability of moons.

Future research could focus on extending the algorithm to model binary-planet systems or incorporate more complex dynamical interactions not covered in the current framework. Additionally, as observational techniques advance, the algorithm could be adjusted to account for direct moon transits detectable with next-generation telescopes.

In conclusion, Kipping's analytic algorithm represents a significant step forward in exoplanetary transit modeling, offering an efficient and comprehensive approach to simulate and analyze the complex interplay of planet-moon-star systems. This not only elevates the prospect of exomoon discoveries but also deepens the broader understanding of the architecture and dynamics of planetary systems beyond our own solar system.