Dual-moon forced dynamics and nonlinear aggregation in Saturn's F ring: From quasi-periodicity to modulated oscillations

Abstract

We develop a minimal nonlinear model to investigate the oscillatory dynamics of Saturn's F ring under dual-moon forcing from Prometheus and Pandora. The model extends classical predator--prey dynamics by incorporating both a nonlinear mass aggregation term $kM^n$ and explicit dual-frequency forcing, capturing how higher-order coagulation physics interacts with multi-moon perturbations. Through extensive numerical integration and dynamical systems analysis, including time-series, spectral, stroboscopic mapping, and rotation number diagnostics, we identify distinct dynamical regimes controlled by the parameters $n$ and $k$. For moderate nonlinearity $(n=1.28, k=0.54)$, the system exhibits regular quasi-periodic motion on a two-torus, characterized by smooth amplitude modulation and discrete spectral lines. As nonlinearity increases $(n=1.30, k=0.62)$, the dynamics transition to strongly modulated oscillations with intermittent phase slips, broadened Poincaré bands, and sideband-rich spectra. A rotation number heatmap reveals organized structures in parameter space, with smooth quasi-periodic regions bounded by near-locking bands analogous to Arnold tongues. Our results demonstrate that the F ring's complex morphology can emerge from deterministic multi-frequency dynamics rather than stochastic processes, with the system operating near critical boundaries where small parameter variations can trigger macroscopic reorganization. The model provides a framework for understanding pattern formation in other driven granular systems while offering testable predictions for ring observations.

Figures (5)

• Figure 1: Upper row: time series and phase portrait over a $10\,T_s$ window. Bottom row: stroboscopic section sampled every $T_s$ over a very long window $(10^4 \ T_S)$ and a periodogram of $M(t)$ over $10\,T_s$. Parameters from Table \ref{['tab:params']} with $(n=1.30, k=0.62)$.
• Figure 2: Same set of diagnostics as Fig. \ref{['fig:Times_Series1']}. Parameters: $(n=1.28, k=0.54)$.
• Figure 3: Unit–circle stroboscopic maps across $(n,k)$. Each panel shows $Z_j=z_j/|z_j|$ obtained by centering the stroboscopic samples $(M(t_j),X(t_j))$, mapping to the complex plane, and normalizing to the unit circle.
• Figure 4: Stroboscopic radial return maps across $(n,k)$. Each panel shows the relation between successive stroboscopic radii $r_{n}$ and $r_{n+1}$ computed from centered samples $(M(t_j),X(t_j))$ taken every $T_s$. The gray diagonal is $r_{n+1}=r_n$.
• Figure 5: Rotation number heatmap over $(k,n)$. Color encodes the stroboscopic rotation number $\rho(k,n)\bmod 1$ computed from long integrations of the two–moon model; grey cells (“No Solution”) indicate parameter pairs for which reliable stroboscopic sequences could not be obtained. Smooth light orange color variation signals quasi–periodic motion with an incommensurate rotation number, whereas the patches of blue-brown regions indicate near-locking frequencies