Deriving the Continuous β-FPUT Model for Decayless Solar Coronal Oscillations Directly from Resistive MHD

Technical Summary: Deriving the Continuous $\beta$-FPUT Model for Decayless Solar Coronal Oscillations Directly from Resistive MHD

Problem Statement
The paper addresses a long-standing paradox in solar coronal seismology: the observation of persistent, low-amplitude transverse oscillations in coronal loops that do not exhibit the rapid damping predicted by classical linear Magnetohydrodynamics (MHD). Traditional models, relying on phase mixing and resonant absorption in transversely inhomogeneous media, predict an asymptotic decay law proportional to $\exp(-t^3)$. However, high-resolution observations from the Solar Dynamics Observatory (SDO) reveal "decayless" oscillations that persist for dozens of periods without measurable attenuation. Prevailing explanations often invoke external driving mechanisms (e.g., background flows or red noise) to sustain these waves, effectively treating the corona as an open system. This work proposes an alternative, self-contained analytical framework that eliminates the need for external forcing by demonstrating that intrinsic non-linear waveguide dynamics can naturally arrest phase mixing.

Methodology
The author derives a continuous fluid model directly from the governing vector equations of resistive MHD, including mass continuity, momentum, and induction equations under zero gravity. The derivation employs a multi-scale perturbation expansion up to third order ($O(\epsilon^3)$) for a line-tied coronal loop modeled as a 2D plasma slab with a uniform longitudinal magnetic field ($B_0 \hat{z}$) and a smooth cross-field density gradient ($\rho_0(x)$).

Key methodological steps include:

1. Symmetry Analysis: The study assumes out-of-plane polarization ($v_y, b_y$) and spatial translation invariance along the displacement coordinate ($\partial/\partial y \equiv 0$). This symmetry causes all quadratic ($O(\epsilon^2)$) advective and Maxwell tension terms to vanish identically.
2. Non-linear Cascade: The analysis tracks a two-step process where the quadratic longitudinal ponderomotive force drives field-aligned density modulations. These modulations back-react on the primary transverse wave, generating a strictly cubic ($O(\epsilon^3)$) self-interaction term.
3. Reduction to mKdV: By projecting the wave field onto the fundamental geometric eigenmode of the line-tied loop and applying a long-wavelength asymptotic expansion, the governing equations are reduced to a Modified Korteweg-de Vries (mKdV) system. This system is identified as the exact fluid-continuum analog of the cubic $\beta$-Fermi-Pasta-Ulam-Tsingou ($\beta$-FPUT) lattice.
4. Numerical Validation: The analytical predictions are tested using numerical simulations of a discrete lattice of $N=32$ coupled magnetic strands. The simulations utilize a second-order symplectic velocity Verlet scheme to integrate the coupled ordinary differential equations. Crucially, to prioritize the uncompromised conservation properties of the mechanical FPUT analogy, the discrete lattice representation explicitly ignores plasma resistivity and physical non-uniformities that would lead to phase mixing along the boundary layers.

Key Contributions and Results

• Derivation of the $\beta$-FPUT Continuum: The paper establishes that the continuous resistive MHD equations, under specific polarization and symmetry constraints, map directly to the mKdV equation. This provides a rigorous physical basis for the $\beta$-FPUT model in a fluid context, where the macro-geometric waveguide gradient ($L_x$) acts as the dispersive mechanism rather than discrete lattice spacing.
• Mechanism for Decayless Oscillations: The study proves that the cubic non-linearity induces an amplitude-dependent frequency self-focusing. This mechanism phase-locks adjacent fluid elements, suppressing the development of cross-field spatial gradients that lead to classical phase mixing. The non-linear coupling timescale ($\tau_{NL} \sim 10^3$ s) is shown to be two orders of magnitude faster than the near-ideal phase-mixing timescale ($\tau_{PM} \sim 10^5$ s), effectively arresting wave attenuation without external driving.
• Closed-Form Coronal Seismology: A significant contribution is the derivation of closed-form expressions linking macroscopic observables to unresolved micro-structural parameters. By inverting the recurrence time of the non-linear energy cascade, the paper provides formulas to calculate:
  • The cross-field density gradient scale length ($L_x$).
  • The internal volumetric filling factor ($f_V$).
    These calculations rely solely on the observed duration of the decayless phase ($\Delta \tau$), the loop length ($L$), the macro-radius ($R$), and the initial displacement amplitude ($a_0$). Crucially, the derivation shows that the internal Alfvén speed and local plasma density cancel out, removing the need for traditional emission measure or spectroscopic line-ratio diagnostics.
• Numerical Confirmation: Simulations of the discrete lattice confirm that the system resists multi-scale thermalization. Instead, the strands execute a highly synchronized, periodic energy exchange, validating the recurrence dynamics predicted by the analytical model within the idealized, non-dissipative limit of the FPUT analogy.

Significance
The paper claims to establish a new paradigm for coronal seismology by shifting from phenomenological open-system models to a self-contained physical framework. By demonstrating that intrinsic non-linear dynamics can sustain oscillations, it resolves the damping paradox without invoking unverified external forcing mechanisms. Furthermore, the framework offers a practical diagnostic tool: it allows observers to infer fine-scale, otherwise unresolvable waveguide properties (specifically $L_x$ and $f_V$) by monitoring the macroscopic duration of oscillation decay. This approach is likened to Enrico Fermi's use of macro-scale displacements to estimate atomic yields, positioning the monitoring of oscillation durations as a method to extract detailed internal structural data of the solar atmosphere. The work provides a theoretical basis for interpreting decayless oscillations as stable, self-sustaining fluid phenomena governed by the continuous $\beta$-FPUT system.