Exomoon
An exomoon is a natural satellite orbiting an exoplanet, a planet located beyond our Solar System, analogous to the moons that accompany planets within our own system but adapted to diverse extrasolar environments.[1] These bodies are expected to form through mechanisms such as giant impacts, accretion in circumplanetary disks, or capture from surrounding space, potentially influencing the host planet's atmosphere, rotation, and habitability.[1] Despite over a decade of targeted searches using space-based telescopes, no exomoons have been definitively confirmed as of November 2025, though candidates continue to emerge from transit photometry data.[2] The primary detection methods for exomoons rely on indirect signatures during planetary transits across their host stars, including transit timing variations (TTVs) caused by the gravitational tug of a moon on its planet, transit duration variations (TDVs) from the combined motion of the planet-moon system, and transit radius variations (TRVs) due to the moon's separate shadow.[1] Instruments like NASA's Kepler Space Telescope, Hubble Space Telescope, and James Webb Space Telescope (JWST) have been instrumental in these efforts, with JWST's infrared capabilities enabling spectroscopy to probe atmospheric compositions potentially influenced by moons.[1] Notable candidates include Kepler-1625b-i, a Neptune-sized moon proposed around a Jupiter-like exoplanet in 2018 based on TTVs and corroborated by Hubble observations, though subsequent analyses have questioned its signal amid noise.[1] Another is Kepler-1708b-i, identified in 2022 as a potential Earth-sized moon around a cold gas giant via reanalysis of Kepler data, representing one of the few candidates around non-hot-Jupiter worlds.[1] Recent advancements have spotlighted even more intriguing possibilities, such as a hypothesized volcanic exomoon orbiting the hot Saturn WASP-39b, approximately 700 light-years away, inferred from variable sodium, potassium, and sulfur dioxide signatures in JWST transmission spectra analyzed in 2025.[3] This candidate suggests tidally heated outgassing akin to Jupiter's moon Io, with simulated mass loss rates indicating a stable orbit within the planet's Roche limit.[3] Such findings underscore exomoons' potential role in enriching planetary atmospheres with volatiles, which could enhance prospects for subsurface oceans or biosignatures on otherwise barren worlds.[1] Exomoons hold profound implications for understanding planetary system architecture and the prevalence of life in the universe, as they may stabilize axial tilts for seasonal climates, drive tidal heating for liquid water, or serve as independent habitable niches.[1] Theoretical models predict that massive exomoons are more likely around giant planets in the outer regions of systems, where disk instabilities favor their formation, contrasting with the inward migration of hot Jupiters that might disrupt satellites. Ongoing surveys with JWST and future missions like the Habitable Worlds Observatory aim to resolve current candidates and detect dozens more, potentially revolutionizing our view of extrasolar diversity.Definition and Fundamentals
Definition of an Exomoon
An exomoon is defined as a natural satellite that orbits an exoplanet, a planet located outside the Solar System. Unlike exoplanets, which are primary bodies orbiting stars, exomoons are secondary objects gravitationally bound to their host exoplanets, distinguishing them from smaller asteroids or comets that may orbit stars independently or temporarily interact with planets.[4] They also differ from rogue planets or hypothetical free-floating moons, as exomoons maintain a stable orbital relationship with a planetary-mass host rather than drifting unbound through interstellar space.[5] The concept of exomoons emerged in the late 1990s, shortly after the discovery of the first exoplanets around main-sequence stars in 1995, when astronomers began considering the potential for satellite systems in extrasolar environments. The initial theoretical proposal for detecting such bodies via transit methods was outlined in 1999, building on the rapid expansion of exoplanet research.[6] Exomoons presuppose the existence of exoplanets as parent bodies, which are diverse in size, composition, and orbital configurations but share the fundamental characteristic of orbiting a star other than the Sun.[7] This hierarchical structure—stars hosting planets, which in turn host moons—mirrors the architecture of our Solar System, providing a foundational framework for understanding extrasolar satellite populations.[4]Nomenclature and Designation
The nomenclature for exomoons extends the established conventions for exoplanets, appending a lowercase "i" to the host planet's designation to denote the first potential satellite, as exemplified by the candidates Kepler-1625 b-i and Kepler-1708 b-i.[8][7] This practice mirrors the International Astronomical Union (IAU) guidelines for exoplanets, which assign provisional designations based on the host star's catalog name followed by a sequential lowercase letter (e.g., b, c) for the planet in order of discovery.[9] Due to the absence of confirmed exomoons, the IAU has not issued specific guidelines for their permanent naming, leaving all current designations provisional and tied to the discovery context rather than standardized approval.[9] Provisional names are typically introduced in peer-reviewed papers reporting candidate signals, such as those from transit photometry surveys, and may vary slightly in notation (e.g., Roman numeral I versus lowercase i) across sources until consensus emerges.[8][7] In astronomical databases, exomoons lack dedicated catalog entries owing to their unconfirmed status. NASA's Exoplanet Archive documents confirmed exoplanets and their hosts but does not include separate parameters or designations for exomoon candidates, referencing them instead through linked publications.[10] Similarly, the SIMBAD database catalogs exoplanet systems under their stellar or planetary identifiers but omits distinct entries for provisional exomoons, treating potential satellites as annotations within the host system's bibliography. The development of exomoon nomenclature traces back to the early 2010s, when theoretical models for exomoon formation and detection began incorporating hypothetical designations in simulations, though without specific candidates.[11] Practical adoption accelerated with the first observational candidate in 2018, establishing the host-planet-plus-i convention in subsequent literature.[8]Physical and Orbital Characteristics
General Properties
Exomoons are expected to exhibit a wide range of physical properties based on theoretical models that draw analogies from Solar System satellites, with masses ranging from approximately 0.01 Earth masses (M⊕) for small icy bodies to up to 1 M⊕ for larger, potentially habitable examples. Sizes may span from dwarf moon scales, around 0.25 Earth radii (R⊕), to Earth-sized objects at 1 R⊕, depending on formation environments around gas giants or super-Earths. Compositions are predicted to vary from predominantly icy mixtures of water, silicates, and minor volatiles for outer exomoons, akin to those orbiting Jupiter and Saturn, to more rocky interiors with metallic cores for inner or larger satellites. These properties provide baselines through comparisons to Solar System moons such as Ganymede and Titan, which serve as archetypes for exomoon expectations. Ganymede, with a mass of 0.025 M⊕, radius of 0.41 R⊕, and density of 1.94 g cm⁻³, features a differentiated internal structure including an iron-rich core, silicate mantle, and thick water-ice layer, offering a model for mid-sized, icy exomoons capable of generating internal magnetic fields. Similarly, Titan, at 0.023 M⊕, 0.40 R⊕, and 1.88 g cm⁻³ density, exemplifies a composition of water ice and rock with a dense nitrogen-methane atmosphere, suggesting that comparable exomoons could retain substantial volatile envelopes under certain orbital conditions. Theoretical estimates for exomoon densities range from 1.5–2.0 g cm⁻³ for icy types to 5.5 g cm⁻³ for rocky ones, with internal structures likely layered—potentially including subsurface oceans between rocky cores and icy mantles—driven by differentiation processes observed in Solar System analogues. Surface features of exomoons are inferred to include heavily cratered terrains from impacts, grooved or ridged regions from tectonic stresses, and possible cryovolcanic deposits, mirroring Ganymede's ancient, modified icy crust or Titan's organic-rich dunes and lakes. Atmospheres, if present, could consist of outgassed volatiles like H₂O vapor (30–460 bar partial pressure) and CO₂ (7–100 bar), particularly for masses above 0.5 M⊕, enabling greenhouse effects or surface oceans on larger bodies while smaller ones might lose atmospheres to space. Geological activity is modeled to arise from tidal heating and radiogenic decay, potentially sustaining cryovolcanism or plate-like tectonics on active exomoons, similar to Europa's subsurface dynamics but scaled to exoplanetary host influences.Orbital Dynamics and Inclination
The orbital dynamics of exomoons are governed by the gravitational interaction between the moon and its host exoplanet, analogous to the Keplerian motion observed in solar system moons. Key parameters include the semi-major axis a_m, which defines the average distance from the moon to the planet's center of mass, typically ranging from just beyond the Roche limit (approximately 2.2 times the moon's radius scaled by the mass ratio (M_p / M_m)^{1/3}) to about 0.49 times the planet's Hill radius for prograde orbits in stable configurations.[12] Eccentricity e_m is generally low to ensure dynamical stability, often constrained such that the periastron distance exceeds the Roche limit and apoastron remains within the Hill sphere, with models assuming values near zero or up to 0.01 for long-term survival.[12] The orbital period P_m follows accordingly, with lower bounds around 3-5 hours for dense rocky moons and upper limits scaled to the planet's orbital period, typically P_m / P_p < 0.2 to avoid perturbations from the host star.[12] These parameters are interrelated through Kepler's third law adapted for the moon-planet system, where the period relates to the semi-major axis via the combined masses of the planet M_p and moon M_m: T^2 \propto \frac{a^3}{M_p + M_m} This formulation, derived from the two-body problem, approximates the moon's orbit as M_p P_m^2 / M_* P_p^2 \approx a_m^3 / a_p^3, linking it to the planet's orbit around the star with semi-major axis a_p and period P_p.[12] In practice, M_m \ll M_p, so the moon's mass has minimal impact unless the satellite is unusually massive relative to the planet. Orbital inclination relative to the planet's equator plays a critical role in the dynamics, distinguishing prograde (co-aligned) from retrograde (opposed) orbits. Prograde orbits, common in the solar system, limit stable semi-major axes to roughly 0.3-0.35 times the Hill radius due to corotation resonances that destabilize wider orbits.[13] Retrograde orbits, by contrast, enhance dynamical stability, extending viable semi-major axes to 0.58-0.61 times the Hill radius—nearly double that of prograde cases—owing to reduced torque from the planet's oblateness and weaker resonance overlaps.[13] This difference implies that retrograde exomoons may persist at greater distances, potentially influencing their tidal interactions and migration paths, though high inclinations can trigger Lidov-Kozai oscillations that excite eccentricity.[14] Tidal forces between the exomoon and its host planet drive secular evolution of these orbits, causing gradual changes in semi-major axis and eccentricity over millions of years. For synchronously rotating planets, tides typically migrate moons outward if the planet spins faster than the orbital motion, increasing a_m while damping e_m toward zero; conversely, inward migration occurs for slower-spinning planets, potentially leading to Roche lobe overflow.[15] Orbital resonances, such as 2:1 mean-motion resonances among multiple moons, further modulate this evolution by inducing eccentricity oscillations (e.g., amplitudes up to 0.15) that librate around equilibrium values, even from initially circular orbits.[16] These resonant interactions, modeled via N-body simulations like REBOUND, can sustain temporary equilibria but often lead to enhanced tidal dissipation and orbital reconfiguration on timescales of 1-100 million years, depending on planetary mass and composition.[16]Stability and Theoretical Constraints
Stability Within Hill Spheres
The Hill sphere represents the gravitational domain surrounding an exoplanet where its influence dominates over that of the host star, allowing potential exomoons to remain bound. This region is approximated by the radius r_H \approx a_p \left( \frac{M_p}{3 M_*} \right)^{1/3}, where a_p is the semi-major axis of the planet's orbit around the star, M_p is the mass of the planet, and M_* is the mass of the star. Long-term stability of exomoons requires their orbits to lie well within this Hill sphere to avoid ejection due to stellar perturbations. Numerical investigations of hierarchical triple systems—consisting of the star, planet, and moon—reveal that circular, coplanar, prograde orbits remain stable up to roughly 0.4895 r_H, beyond which chaotic ejections occur on timescales of $10^4 to $10^5 orbital periods of the planet.[17] This critical fraction arises from the dynamics of the restricted three-body problem and has been confirmed through extensive N-body simulations integrating over millions of initial conditions. Several factors modulate the extent and duration of stable exomoon orbits within the Hill sphere. A higher planetary mass M_p expands r_H, permitting larger stable orbital radii for moons, while closer stellar proximity (smaller a_p) contracts the sphere, compressing the viable region. Additionally, multi-body perturbations from co-orbiting planets can narrow the stable zone by inducing secular resonances or close encounters that accelerate instabilities.[18] Simulations tailored to diverse exoplanet architectures demonstrate that exomoons positioned inside approximately 0.5 r_H can endure for gigayears, aligning with typical system ages. For instance, in temperate giant planet systems like those surveyed around Kepler targets, tidal evolution models over 10 Gyr indicate stable orbits for moons without significant eccentricity growth or loss, provided they avoid the outer Hill boundary.[7] In the specific case of the candidate exomoon around Kepler-1625b, orbital stability analyses confirm lifetimes exceeding 4.5 Gyr for semi-major axes up to 0.4 r_H, highlighting the robustness of inner Hill sphere regions against perturbations.[19]Absence of Close-In Exomoons
Close-in exomoons, those orbiting planets in tight orbits around their host stars such as hot Jupiters, face significant challenges to their stability due to intense stellar tidal forces. These tides exert a disruptive influence that can strip moons from their parent planets over time, as the strong gravitational gradient from the nearby star overcomes the planet's hold on its satellites. Theoretical models indicate that for planets with orbital periods less than about 10 days, the stellar tidal torque accelerates the outward migration of moons until they exceed the planet's Hill sphere and are lost to the star. This tidal disruption is exacerbated by the smaller Hill spheres of close-in planets, which limit the region where moons can remain bound—a concept rooted in the balance between planetary and stellar gravity as discussed in general orbital stability analyses. Observational biases further compound this issue, as the reduced Hill radius around hot Jupiters confines potential moons to orbits too compact for retention against tidal forces, leading models to predict low retention rates of substantial moons over gigayear timescales. These predictions arise from simulations incorporating tidal evolution, where prograde moons are particularly vulnerable to rapid ejection during the planet's post-formation phase. A 2025 study suggests that massive retrograde moons greater than 10 Earth masses may survive in approximately 6% of simulated hot Jupiter migration scenarios, potentially allowing for rare close-in exomoons.[20] Empirical evidence supports this theoretical scarcity, with no confirmed exomoon candidates identified among the more than 100 transiting hot Jupiters scrutinized in exoplanet surveys using transit timing variations and other indirect signals. For instance, analyses of Kepler data for hot Jupiters exhibiting potential perturbations have attributed variations to other causes, such as additional planets, rather than moons, reinforcing the absence of detectable exomoons in these systems as of November 2025. The rarity of close-in exomoons has broader implications for understanding planetary system architecture, particularly in the context of high-eccentricity migration models for hot Jupiters. During inward migration driven by disk interactions or other mechanisms, moons are often detached and scattered, either contributing to debris disks or becoming free-floating objects, which may explain the observed isolation of many close-in giants without satellite companions. This process highlights how dynamical evolution shapes the final configuration of mature planetary systems.Formation and Evolution
Formation Mechanisms
The primary mechanisms proposed for exomoon formation draw analogies from Solar System moon origins but adapt to the diverse environments of exoplanetary systems, where giant planets may form farther from their stars or undergo significant dynamical evolution. These include accretion within circumplanetary disks, debris generation from giant impacts, and gravitational capture of external bodies. Models from the 2010s onward emphasize how these processes can yield moons up to Mars mass or larger, potentially detectable around super-Jovian exoplanets.[21][11] Circumplanetary disk accretion, the dominant pathway for regular moons in the Solar System, involves the in situ formation of satellites from gas and dust disks surrounding nascent gas giant planets. As the planet accretes its envelope, material is captured into a subnebula where solids coagulate and gas condenses, leading to moon formation via core accretion or disk instabilities. Numerical models predict satellite-to-planet mass ratios around 10^{-4}, sufficient for Mars-sized exomoons around planets exceeding Jupiter's mass, with examples including the Galilean satellites of Jupiter formed in a gas-starved disk. This mechanism is most efficient beyond 5 AU, where protoplanetary disks supply ample material.[21][11] Giant impacts offer an alternative for forming massive, rocky exomoons, particularly around terrestrial or super-Earth exoplanets, through collisions that eject debris into orbit. Simulations of oblique impacts between protoplanets with masses 0.25 to 10 Earth masses, at velocities near escape speed, generate circumplanetary disks containing 0.1 to several Earth masses of material, which can coalesce into one or more moons. This process mirrors the formation of Earth's Moon from a Mars-sized impactor and could produce exomoons up to 0.5 Earth masses, with higher efficiencies for impacts involving 10-50% mass ratios at 10-15 km/s. Such events are plausible during the late stages of planet formation in dynamically active systems.[22][21] Capture mechanisms enable the acquisition of exomoons from external sources, such as interstellar objects or planetesimals, without requiring local disk material. In the "pull-down" scenario, a pre-formed protoplanet or rogue body is drawn into a stable orbit during the rapid gas accretion phase of a giant planet's envelope, aided by gas drag for energy dissipation. Binary exchange, where a captured pair is disrupted and one component retained, can also operate, as inferred for Neptune's Triton. These processes favor irregular, distant exomoons and may incorporate material from the interstellar medium, contrasting with the more ordered disk accretion.[23][21] Exomoon formation differs from Solar System patterns due to the prevalence of planetary migration in exoplanet systems, which can populate outer orbits with moons via dynamical scattering or enhanced capture opportunities. Unlike the relatively static formation of Jupiter's moons, migrating gas giants may sweep up additional material during inward or outward excursions, leading to hybrid systems with both accreted and captured satellites; population synthesis models from the 2010s indicate this can increase exomoon frequencies around super-Jovians by factors of 2-5 compared to non-migrating cases. Overall, these mechanisms predict a broader range of exomoon sizes and compositions than observed locally.[11][21]Evolutionary Processes
Exomoons undergo significant tidal evolution over billions of years, primarily driven by interactions between the moon, its host planet, and the central star. These tides cause orbital migration, where the moon's semi-major axis can expand or contract depending on the relative strengths of planetary and stellar tides. For close-in exomoons around giant planets, inward migration toward the planet is common due to dominant planetary tides, potentially leading to orbital decay and eventual disruption if unchecked. However, the evolving physical properties of the host planet—such as contraction of its radius and changes in internal structure—can reverse this trend, inducing outward migration and preventing collision or tidal breakup. Additionally, tidal torques synchronize the moon's rotation to its orbital period around the planet, resulting in synchronous rotation where one hemisphere permanently faces the host, a process that typically completes within a few million years for Earth-mass moons orbiting gas giants.[24][25] Atmospheric evolution in exomoons is shaped by stellar radiation and the host planet's magnetosphere, leading to potential loss or, in some cases, secondary gain through outgassing. Hydrodynamic escape, powered by extreme ultraviolet (XUV) radiation from young stars, is the primary loss mechanism, stripping hydrogen-rich envelopes at rates of approximately 10^{31}–10^{33} atoms per second for moons of 0.1–1 Earth masses at 1 AU from a Sun-like star. Smaller exomoons (≤0.1 M_⊕) lose their entire atmospheres within 1 million years under high XUV fluxes, while more massive ones (≥0.5 M_⊕) may retain substantial envelopes, evolving toward Earth-like compositions. In cooler stellar environments like those around M-dwarfs, prolonged high activity extends loss periods, but planetary magnetic fields can shield exomoons, reducing sputtering and ion pickup losses. Atmospheric gain occurs via volcanic or cryovolcanic outgassing, replenishing lost volatiles and potentially forming secondary atmospheres dominated by CO_2 or N_2 in mature systems.[26][25] Geological evolution of exomoons involves differentiation, volcanism, and cryovolcanism, largely powered by tidal and radiogenic heating. For icy exomoons with masses between 0.1 and 0.5 Earth masses, internal heating drives differentiation into iron-rich cores, rocky mantles, and icy shells, similar to Ganymede's structure with a moment of inertia factor of about 0.31. Tidal heating in close orbits (<10 planetary radii) induces volcanism on rocky or mixed-composition moons, producing heat fluxes exceeding 2.4 W/m² for Mars-mass bodies, sustaining activity for up to 100 million years and leading to silicate resurfacing akin to Io. Cryovolcanism predominates on icy exomoons, where tidal friction melts subsurface oceans, enabling eruptions of water-ammonia mixtures through geysers or flows, as inferred from analogs like Enceladus; this process can persist for billions of years if eccentricities are maintained by resonances, fostering global resurfacing and potential subsurface habitability. Radiogenic decay in larger moons supplements tidal energy, enhancing differentiation and long-term geological activity. Host planet migration profoundly impacts exomoon orbits, often accelerating decay through altered tidal regimes and dynamical instabilities. During inward planet migration, such as in high-eccentricity or disk-driven scenarios, exomoons within the shrinking Hill sphere experience enhanced stellar tides, hastening inward orbital migration and reducing stability lifetimes to 10^8–10^9 years for initial semi-major axes of 5–20 planetary radii. Simulations of planet-planet scattering reveal that up to 50% of exomoons are ejected or collide with the planet during chaotic phases, while survivors may expand outward if the planet's spin synchronizes with the moon's orbit, countering decay. In synchronized systems, exomoons can extend the planet's survival against stellar tidal inspiral by dissipating energy more efficiently, with decay timescales increased by orders of magnitude compared to moonless planets. These effects underscore the fragility of exomoon systems during planetary dynamical evolution.[27]Detection Methods
Imaging and Direct Techniques
Direct imaging techniques seek to spatially resolve exomoons from their host exoplanets using high-contrast observations, primarily in the mid-infrared where thermal emission from the moon can outshine reflected visible light.[28] This approach relies on coronagraphs to suppress starlight and adaptive optics to mitigate atmospheric distortion, aiming to detect the moon-planet separation directly.[28] The method is particularly suited to tidally heated exomoons (THEMs), whose internal heating from orbital eccentricity produces elevated temperatures and luminosities.[28] A major challenge is the extreme contrast ratio required—up to 10^6 or more—between the moon's faint signal and the host planet's glare, compounded by small angular separations typically below 0.1 arcseconds for nearby systems.[28] Current ground-based telescopes like the Very Large Telescope (VLT) with the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument have demonstrated contrasts sufficient for Jupiter-sized planets but fall short for resolving sub-Earth-sized moons without extended integrations.[28] Space-based platforms avoid atmospheric limitations, yet even they struggle with variable moon emission due to hotspots or partial eclipses by the planet.[28] Expected observational signatures include broadband thermal emission peaks in the 5–20 μm range for THEMs at 300–1000 K, potentially reaching luminosities of 0.1% of the host star for low-mass stars.[28] For an Earth-radius moon at 600 K orbiting a gas giant 5 pc away, detection at 5σ significance requires about 1 hour of exposure with instruments like Spitzer, though signal-to-noise drops rapidly for cooler or smaller bodies.[28] Reflected light signatures are weaker and harder to isolate without multi-wavelength color-magnitude analysis.[29] No exomoons have been directly imaged to date, with historical attempts focused on theoretical simulations rather than targeted observations due to instrumental constraints. Early proposals emphasized THEMs around young, wide-orbit gas giants, but no candidates have yielded resolved detections despite surveys of directly imaged exoplanets.[29] The James Webb Space Telescope (JWST), operational since 2022, enhances feasibility through its Mid-Infrared Instrument (MIRI); pre-launch predictions suggested it could image Earth-sized THEMs warmer than 300 K around the 24 nearest stars (within 4 pc) at 5σ in roughly 10,000 seconds of integration, with coronagraphic mode achieving contrasts of 10^{-5} to 10^{-6} at 0.5–1 arcseconds, sufficient for separations of 10–30 AU in nearby systems.[28] However, as of November 2025, JWST/MIRI searches, such as those around the rogue planet WISE 0855, have not yielded statistically significant exomoon detections.[30] Future ground-based extremely large telescopes, such as the Thirty Meter Telescope, could extend this to slightly more distant targets with adaptive secondary mirrors. Detection remains limited to wide-separation orbits (>10 AU) where contrast is more favorable and tidal heating sustains moon luminosity, primarily in young systems (<100 Myr) to avoid cooling over gigayear timescales.[28] Close-in exomoons are infeasible due to unresolved blending with the planet and reduced heating efficiency.[29]Spectroscopic and Timing Methods
Spectroscopic methods for detecting exomoons primarily rely on measuring subtle perturbations in the radial velocity of the host exoplanet caused by the gravitational pull of an orbiting moon. In Doppler spectroscopy, the moon induces a periodic wobble in the planet's motion around their common barycenter, manifesting as variations in the planet's radial velocity as observed from Earth. The amplitude of this velocity shift, \Delta v, is approximated by \Delta v \approx \left( \frac{M_m}{M_p} \right) v_p \sin i, where M_m is the moon's mass, M_p is the planet's mass, v_p is the planet's orbital velocity around the star, and i is the inclination of the planet-moon system relative to the line of sight.[31] This signal is particularly promising for directly imaged young, self-luminous gas giants, where high-resolution spectrographs can resolve the moon-induced perturbations on top of the planet's own stellar reflex motion.[32] For instance, observations of the planet HR 7672 B using the Keck Planet Imager and Characterizer have set limits on potential exomoons by constraining radial velocity semi-amplitudes below 100 m/s, highlighting the method's sensitivity to moon-to-planet mass ratios as low as 0.01.[31] Another spectroscopic approach involves detecting radio emissions generated by interactions between the exomoon and the host planet's magnetosphere, analogous to Jupiter's moon Io producing auroral radio bursts through plasma torus interactions. In gas giant systems, an exomoon can supply material to the planet's magnetosphere via volcanism or atmospheric stripping, leading to modulated decametric radio emissions that vary with the moon's orbital phase.[33] These emissions, detectable at frequencies around 10-40 MHz, could reveal exomoons around mature Jupiter-like exoplanets within 10-30 parsecs using ground-based radio telescopes like the Low Frequency Array (LOFAR).[33] Recent searches with the Giant Metrewave Radio Telescope (GMRT) on candidate "exo-Io" systems such as WASP-49 and HAT-P-12 have not detected such emissions as of 2022.[34] For multiple-exomoon systems, overlapping plasma tori may produce complex emission patterns, but single large moons like exo-Ios are the most promising targets due to their stronger Alfvén wave generation. Timing methods complement spectroscopy by analyzing variations in the host planet's transit events across a star's light curve. Transit timing variations (TTVs) arise from the gravitational tug of the moon, causing the planet to arrive early or late at the transit midpoint by amounts proportional to the moon's mass and semi-major axis, with the TTV amplitude scaling as \propto M_m a_m, where a_m is the moon's orbital distance.[35] These periodic deviations, typically on the order of minutes to hours for Earth-sized moons around Neptune-mass planets, can be extracted from high-precision photometry like that from the Kepler or TESS missions, though aliasing from sparse sampling often confines detectable signals to a narrow "exomoon corridor" in frequency space.[36] Complementary transit depth variations occur when the moon's shadow aligns with the planet's during transit, altering the combined occultation depth by a factor related to the moon's radius relative to the planet's, potentially detectable at parts-per-million levels in ultra-precise light curves.[37] For example, fits to Kepler data have used these depth perturbations as indicators of exomoons by revealing non-constant ingress/egress shapes inconsistent with a lone planet model.[37]Gravitational and Indirect Methods
Gravitational microlensing offers a method to detect exomoons by observing temporary increases in brightness when a planet-moon system aligns with a background star from Earth's perspective, where the moon's gravity can cause additional lensing anomalies distinct from the planet's primary event.[38] These anomalies manifest as secondary peaks or caustics in the light curve, allowing inference of the moon's mass and separation if the sampling is sufficiently dense and precise.[38] Simulations indicate that exomoons with masses around 0.01 Earth masses or larger could be detectable in events monitored by ground- or space-based surveys, particularly for wide-orbit systems where the moon's independent lensing is resolvable.[38] Upcoming joint surveys like Euclid and Roman are projected to enhance sensitivity, potentially identifying multiple exomoons per event through high-cadence photometry that resolves sub-percent deviations. Astrometric methods detect exomoons by measuring the positional wobble of the host exoplanet around the planet-moon barycenter projected on the sky, using high-precision astrometry or optical interferometry. The astrometric signature scales with the moon-to-planet mass ratio and orbital separation, enabling detection of Earth-mass exomoons around Saturn-mass planets at distances of several parsecs with ground-based telescopes observing over several years, or even closer systems with space-based astrometry. This approach is particularly suited for directly imaged exoplanets and avoids degeneracies in radial velocity methods, with recent proposals highlighting its promise for future facilities like the Habitable Worlds Observatory.[39] Pulsar timing provides another indirect avenue for exomoon detection, exploiting the precise measurement of pulse arrival times from millisecond pulsars to identify perturbations caused by orbiting planet-moon systems. In such systems, the barycenter motion induced by the moon alters the planet's orbital path around the pulsar, producing periodic residuals in the time-of-arrival (TOA) data with amplitudes scaling as the moon-to-planet mass ratio times the orbital separation. For known pulsar planets like those around PSR B1257+12, a moon with a separation of about 1% of the planet's orbit and a mass ratio of 10^{-3} could induce detectable TOA variations of order 10 microseconds, feasible with current radio telescope arrays. This method is particularly sensitive to close-in, massive exomoons but is limited by the rarity of pulsar-planet systems and requires long-term monitoring to distinguish moon signals from other astrophysical noise. Orbital sampling effects (OSE) enable the statistical detection of exomoons through subtle deviations in phase-folded transit light curves of exoplanets, where the moon's photocentric motion causes non-uniform sampling of the planet's transit depth across its orbital phases.[40] This results in an asymmetric or modulated transit profile in stacked observations, with the effect's amplitude proportional to the moon's radius relative to the planet and inversely to the orbital distance, allowing estimation of moon properties without resolving individual moon transits.[40] Analyses of Kepler data have searched for OSE signatures in thousands of exoplanet light curves, revealing potential signals in systems with transit depths varying by up to 0.1% over the planet's orbit, though confirmation requires modeling to rule out stellar variability or instrumental artifacts.[41] The method's efficacy depends on high-precision photometry and large sample sizes, making it suitable for surveys like TESS to infer exomoon prevalence indirectly via population-level statistics.[41] Indirect evidence for past exomoons around white dwarfs can be gleaned from perturbations in circumstellar debris disks, where destabilized moons contribute to irregular structures or variable accretion rates observed in infrared excess or atmospheric pollution.[42] N-body simulations show that exomoons, upon planetary disruption during post-main-sequence evolution, can scatter into the disk, inducing density waves or gaps that manifest as asymmetric thermal emission profiles detectable by telescopes like Spitzer or JWST.[42] For white dwarfs with cooling ages of 100-500 million years, such perturbations from moon masses exceeding 0.1 Earth masses could explain observed disk variabilities on timescales of years, distinct from pure planetesimal collisions.[43] This approach traces historical exomoon presence rather than current ones, with detectability enhanced in systems showing both polluted atmospheres and warped disks, as modeled in dynamical studies of tidal evolution.[43] Stability constraints on exomoon orbits around evolving hosts further influence the timing and amplitude of these disk perturbations.[42]Known Candidates
List of Candidates
As of November 2025, no exomoons have been definitively confirmed, but approximately 10 candidates have been proposed based on anomalous signals in observational data. These candidates meet a minimum threshold of statistical significance greater than 3σ, often derived from deviations in expected planetary transit profiles or additional photometric or spectroscopic features attributable to a companion moon. The vast majority originate from NASA's Kepler mission, which surveyed thousands of transiting exoplanets for moon-like signatures, supplemented by emerging data from the Transiting Exoplanet Survey Satellite (TESS) and the James Webb Space Telescope (JWST). Ground-based microlensing surveys have contributed marginal cases as well.[7][8][44] The following table summarizes the most prominent candidates, focusing on those with the strongest evidence:| Candidate Name | Host Planet | Discovery Year | Detection Method | Statistical Significance | Reference |
|---|---|---|---|---|---|
| Kepler-1625b-i | Kepler-1625 b (Jupiter-sized gas giant, ~0.7 AU orbit) | 2018 | Transit timing variations (TTVs) and transit depth anomalies in Kepler photometry, confirmed with Hubble Space Telescope follow-up | ~3σ | Teachey & Kipping (2018) |
| Kepler-1708b-i | Kepler-1708 b (Jupiter-sized gas giant, ~1.6 AU orbit) | 2022 | Moon-induced "photoeccentric" signal and TTVs in Kepler light curves | 4.8σ | Kipping et al. (2022) |
| Unnamed (WASP-49 b candidate) | WASP-49 b (Hot Jupiter, ~0.4 AU orbit) | 2024 | Spectroscopic detection of a persistent sodium cloud with orbital velocity suggesting an 8-hour period moon, observed with ESO's Very Large Telescope | >3σ | Mankovich et al. (2024) |
| WASP-39b-i | WASP-39 b (hot Saturn-mass gas giant, ~0.29 AU orbit) | 2025 | Variable Na, K, SO₂ signatures in JWST transmission spectra indicating tidally heated outgassing | Preliminary evidence (inferred >3σ) | Oza et al. (2025, preprint) |
| MOA-2011-BLG-262 Lb | MOA-2011-BLG-262 L (Free-floating gas giant, mass ~4 Jupiter masses) | 2013 | Perturbations in microlensing light curve from MOA survey, indicating a sub-Earth-mass moon | ~2σ (marginal, requires follow-up) | Bennett et al. (2014) |