Finite-time stability of two-soliton solutions of the Hartree equation with elliptic trajectories

Abstract: We prove the finite-time stability of two-soliton solutions of the three-dimensional gravitational Hartree equation whose centers remain close to an elliptic solution of the Kepler two-body problem.

• The paper demonstrates finite-time H¹ stability of two-soliton configurations with centers following elliptic Keplerian orbits in the 3D Hartree equation.
• It employs high-order accelerated approximations and quantitative modulation analysis to achieve ultra-small error bounds over extended time scales.
• Localized energy methods combined with harmonic cancellation techniques ensure that parameter dynamics remain tightly controlled despite persistent nonlocal interactions.

Finite-Time Stability of Two-Soliton Solutions of the Hartree Equation with Elliptic Trajectories

Introduction and Context

This work rigorously addresses the finite-time stability of two-soliton solutions of the three-dimensional gravitational Hartree equation, focusing specifically on the challenging case where the soliton centers execute elliptic Keplerian orbits. The Hartree equation,

$$iu_t+\Delta u - \phi_{|u|^2}u=0,\quad \phi_{|u|^2} = -\frac{1}{4\pi|x|}*|u|^2,$$

is a canonical nonlocal dispersive model appearing as an effective description in quantum many-body systems. Its soliton solutions, generated from the ground state $Q$ by symmetry transformations, are exponentially localized, but the nonlocal nonlinearity induces only polynomial decay in the self-consistent potential, leading to persistent long-range interactions atypical of local NLS-type models.

Previous results in the field established the existence and stability of multisoliton states in the Hartree framework only for dynamically expansive (hyperbolic or parabolic) configurations, where soliton centers recede and interaction errors decay in time [Krieger–Martel–Raphaël 2009; Wu 2026]. Contrastingly, the elliptic case, characterized by bounded inter-soliton separation and persistent strong interaction, had not been resolved: the error analysis therein cannot exploit temporal decay. This distinguishes the present work both in technical approach and in the dynamical regime analyzed.

Main Results and Methodology

The principal achievement is the finite-time $H^1$-stability of two-soliton Hartree solutions whose centers follow approximate elliptic Keplerian trajectories. More precisely, given any prescribed small $\delta>0$, uniform eccentricity bound $C_e$, and scaling parameters, it is proved that initial data sufficiently close (up to $O(r_0^{-3})$ in $H^1$, with $r_0$ the characteristic scale of the ellipse) to a carefully constructed two-soliton profile generates a solution remaining within $\delta$ of a reference two-soliton configuration, with centers tracking an elliptic orbit, for times up to $O(r_0^2)$. This timescale notably exceeds the orbital period $Q$0, allowing control over many cycles.

The argument has several substantial components:

• High-Order Accelerated Approximation: The construction of initial data exploits an $Q$1-th order approximate two-soliton profile, where $Q$2, engineered so that the local error is exponentially small in $Q$3 and uniform in time over $Q$4. This is crucial due to the absence of decay with time in the elliptic regime; errors must be controlled for all times considered, relying on expansion order.
• Quantitative Modulation Analysis: The evolution of the soliton parameters (center location, velocity, phase, scale) is governed by a system that, up to high order, realizes the effective two-body Kepler problem. Technical advances include the precision of the error estimates and demonstrating that lower-order dipole effects are exactly canceled due to a nontrivial interaction phase adjustment (termed interaction shift).
• Localized Energy Methods: Stability is routed through a localized coercive functional, refined to accommodate the nonlocality and lack of spatial separation between the solitons. The error growth in time is sublinear and can be bounded using Gronwall-type arguments due to the error's exponential smallness.

The accumulated error, energy growth, and parameter drift remain of order $Q$5 over the full timespan due to the combination of high-order approximation, coercivity, and intricate ODE analysis of the parameter flow around the elliptic solution.

Strong Claims and Quantitative Results

Explicit uniform time control: The theorem establishes that for initial separation $Q$6, two-soliton configurations with centers at elliptic separation remain $Q$7-close in $Q$8 norm for all $Q$9, with the profile error on initial data matching the approximation order, i.e., $H^1$0.

Precision of Parameter Dynamics: The approximate two-soliton dynamics deviates from true elliptic motion only by $H^1$1 in the effective equations for center-of-mass and velocity. This is nontrivial: the analysis shows that all lower-order terms vanish, owing to nontrivial harmonic cancellations and the careful design of the interaction phase.

Necessity of High-Order Approximation: The analysis demonstrates that, due to factorial growth in the expansion coefficients, the approximation order must be $H^1$2; otherwise, the error is insufficiently small to control the solution beyond $H^1$3one period.

Technical Innovations

Interaction Shift and Harmonic Cancellation

The main technical innovation is the introduction and systematic utilization of the interaction shift in the phase evolution equation, ensuring all first-order (and many higher-order) terms in the parameter ODE expansions are cancelled, leaving only genuinely high-order residuals. This leverages subtle representations via spherical harmonics, allowing certain terms in the error expansion to be shown identically zero due to orthogonality properties, thereby reducing the effective error to an ultra-small remainder.

Cutoff and Quantitative Admissibility

A further refinement involves introducing spatial cutoffs in the Taylor expansion of the nonlocal potential to ensure the uniformity of the error with respect to the parameters and spatial localization—a step essential for the elliptic regime where the centers do not separate.

Finite-Time Versus Asymptotic Stability

While previous work could consider the infinite-time stability due to asymptotic spatial separation in the non-elliptic regime, here the analysis is necessarily finite-time; persistent interaction prohibits extending these results to global-in-time stability, and achieving long but finite timescales requires fundamentally new approximation and error-propagation approaches.

Implications and Future Directions

This work closes a significant gap in the theory of nonlinear, nonlocal dispersive PDEs by showing that bounded soliton complexes (in the three-dimensional Hartree/mean-field gravitational model) can realize stable, long-lived 'molecular' states even under significant nonlinear, nonlocal interaction—a sharp contrast with the local NLS case.

From a dynamical systems perspective, the result illuminates the interplay between nonlocality, symmetry, and multipole cancellations in PDE effective dynamics. It raises questions about whether full asymptotic (infinite time) stability in the elliptic regime is ever possible, and it provides both new tools and motivations for considering more general, possibly many-soliton (with $H^1$4) bound states in nonlocal dispersive models.

For applications in mathematical physics—such as models for boson stars, galaxy dynamics, or quantum mean-field theory—these findings confirm the plausibility of persistent, strongly-coupled 'bound' solitary structures beyond the asymptotic regime. On a methodological level, the approach outlined in this work, combining high-order approximation, parameter modulation, and norm localization, is expected to generalize to other nonlocal PDEs and could be useful for studying strongly interacting coherent structures in higher dimensions and more complex nonlinearities.

Conclusion

In summary, the paper delivers the first rigorous construction and finite-time stability result for elliptic (non-expanding) two-soliton states in the 3D Hartree equation, overcoming the principal challenge posed by persistent nonlocal interaction. Through delicate harmonic analysis, modulation theory, and localized energy estimates, the study sets a new standard for precision and rigor in the treatment of multi-soliton dynamics in nonlocal dispersive systems. Future advancements are likely to focus on extending these techniques to global stability, higher soliton numbers, and broader classes of nonlocal nonlinear PDEs.

References

• Krieger, Martel, and Raphaël, "Two-soliton solutions to the three-dimensional gravitational Hartree equation" [KMR2bodyHartree, math/0606094]
• Y. Wu, "Existence of multisoliton solutions of the gravitational Hartree equation in three dimensions" [Hartree3Dmultisoliton, (Daghero et al., 2022)]
• E. Lenzmann, "Uniqueness of ground states for pseudorelativistic Hartree equations" [Lenzmanngroundstate, math/0702148]