NePTune: Neptune & Exoplanet Dynamics

• NePTune is a research concept that examines Neptune’s unique orbital dynamics and exoplanetary migration with precise geometric and simulation techniques.
• It studies the formation of the Neptune ridge and desert through f-mode oscillations, tidal interactions, and hydrodynamical shock processes.
• The findings inform planetary formation theories and enhance methods for measuring orbital parameters and atmospheric properties with high accuracy.

NePTune designates a set of distinct research concepts and methodologies related to Neptune, the eighth planet in the Solar System, and analogous exoplanetary phenomena involving Neptune-mass objects. The term is most prominently featured in contexts such as (i) the "Neptune ridge" and "Neptune desert" in exoplanet orbital distributions, as explained by the role of tides and high-eccentricity migration in dynamical astrophysics, and (ii) methods for precise determination of Neptune's orbital and physical parameters using geometric and observational techniques.

1. Neptune Ridge and Desert: Dynamical Origins

The period distribution of hot gaseous exoplanets reveals pronounced mass-dependent features: a "Neptune ridge" (clustering of sub-Saturns/Neptunes at periods $P \sim 3$–$6$ days) and a "Neptune desert," a sharp deficit of such planets at $P \lesssim 3$ days, in contrast to the continuous presence of hot Jupiters at shorter periods. This structure is explained in the framework of high-eccentricity migration (HEM), where planets with initially large semi-major axes experience eccentricity excitation via secular or scattering interactions, driving repeated pericenter passages near their Roche radii (Zanazzi et al., 18 Jun 2026).

During each periastron passage, the star's tidal gravitational potential excites the planet's $\ell=2$, $m=2$ fundamental oscillation mode (f-mode). The f-mode, evolving on the dynamical timescale $t_p = (G M_p / R_p^3)^{-1/2}$, acts as a reservoir for orbital energy. This oscillation can attain amplitudes sufficient to drive supersonic radial velocities, generating shocks in the gaseous envelope.

Whether the f-mode energy is dissipated by radiative diffusion or by driving a hydrodynamical wind is critical to the outcome. If the deposited energy cools radiatively, orbital circularization deposits the planet at $P \sim 3$–$6$ days ("ridge"/"pile-up"). For deep shocks with deposited energy exceeding the local Eddington luminosity, a super-Eddington wind unbinds the envelope, exposing the solid core and leading to a deficit ("desert") at shorter periods.

2. Hydrodynamical Shocks, Cooling Regimes, and Mass Loss

Three-dimensional Athena++ simulations demonstrate that when pericenter passages occur at $r_p = 1.8$–$2.2\, r_t$ (with $6$0, the tidal radius), f-mode excitation rapidly leads to shock formation wherever the surface radial velocity $6$1. The shock-heated envelope dissipates energy either via radiative diffusion (for shocks at moderate optical depth) or enters a wind-driven regime (for deeper, more energetic shocks where power exceeds $6$2). The boundary between these outcomes in parameter space directly maps onto the edges between the Neptune ridge and Neptune desert (Zanazzi et al., 18 Jun 2026).

The process can be formalized in an iterative map over successive pericenter passages:

• Amplitude growth and damping of the f-mode,
• Stepwise updates to semi-major axis, eccentricity, and planetary mass,
• Thresholds for mass loss set by the "wave-breaking" radius and mass above,
• Post-mass-loss structure recalculated using MESA+GYRE tabulated models.

Repeated super-Eddington outflows strip the atmosphere in $6$3 orbits, leaving bare cores on short-period orbits characteristic of the Neptune desert.

3. Quantitative Predictions and Observational Signatures

For a $6$4 star and a $6$5 planetary core with a $6$6–$6$7 H/He envelope, the inner boundary for envelope loss, $6$8, corresponds to orbital periods $6$9 days—matching the observed inner edge of the Neptune desert. The outer "ridge" at $P \lesssim 3$0 aligns with $P \lesssim 3$1 days, and the mechanism naturally produces a deficit of sub-Saturns at $P \lesssim 3$2 days with clustering at $P \lesssim 3$3–$P \lesssim 3$4 days. Hot Jupiters, with higher densities and different f-mode thresholds, display comparable but shifted features.

Desert Neptunes are predicted to exhibit large spin–orbit misalignments, consistent with HEM, and may produce observable optical/infrared "tidal flares" reaching $P \lesssim 3$5–$P \lesssim 3$6 during envelope loss events (Zanazzi et al., 18 Jun 2026).

4. Geometric Determination of Neptune's Orbital Parameters and Mass

A direct geometric method exploiting Newtonian perturbation theory allows recovery of Neptune’s six orbital elements and mass solely from high-precision Cartesian position and velocity data for the Sun and the remaining planets (Bhatnagar et al., 2021). The approach:

• Analyzes the net acceleration vector imposed by Neptune on Uranus,
• Extracts conjunctions and oppositions using meridian-plane geometry,
• Derives the synodic and sidereal periods, and from these, the semi-major axis ($P \lesssim 3$7 AU) and eccentricity ($P \lesssim 3$8),
• Determines inclination, longitude of ascending node, argument of perihelion, and mean anomaly from the direction vector sequence,
• Computes Neptune’s mass as $P \lesssim 3$9 kg by inverting the perturbation equation over multiple epochs.

All computations are fundamentally geometric and require only Newtonian vector mechanics and root-finding, with 1–2% accuracy in orbital distance and $\ell=2$0 in mass. This methodology reinforces the physical link between planetary perturbations and orbital reconstruction (Bhatnagar et al., 2021).

5. Constraints on Neptune's Ionospheric H$\ell=2$1 Emission

Observations targeting the molecular ion H$\ell=2$2 in Neptune's upper atmosphere using the IRTF/iSHELL instrument yielded no detected emission in the $\ell=2$3–$\ell=2$4 µm region, after 15.4 h of integration (Melin et al., 2017). The derived $\ell=2$5 upper limit on the vertical column density, $\ell=2$6 at $\ell=2$7 K, is at least a factor of 5 lower than predicted by prevalent photochemical models.

This discrepancy is attributable to efficient quenching of H$\ell=2$8 by polyatomic neutrals (e.g., CH$\ell=2$9, CO$m=2$0, H$m=2$1O) lofted to thermospheric altitudes by strong vertical mixing. The resulting destruction of H$m=2$2 is expressed by the loss reaction $m=2$3, imposing a lower limit on the eddy diffusion coefficient $m=2$4–$m=2$5 cm$m=2$6 s$m=2$7. Subsequent JWST/NIRSpec observations are poised to lower detection thresholds an order of magnitude further, probing fundamental questions regarding Neptune’s ionospheric structure and vertical transport.

6. Implications for Planet Formation and Demographic Structure

The physical processes underlying the Neptune ridge and desert illuminate the fate of sub-Saturn planets during migration. The interplay of f-mode energy injection, dissipation mechanisms (radiative cooling vs. wind loss), and planetary structural response accounts for the sharply defined period-mass substructure observed in exoplanet surveys. Efficient formation of bare cores at short periods via envelope stripping reconciles the paucity of Neptune-mass planets at $m=2$8 days and the observed accumulation at longer periods ("ridge"), offering a direct, quantitative connection between dynamical tides and planetary demographics (Zanazzi et al., 18 Jun 2026).