Stability and instability of a one-dimensional MHD model

1. Problem Statement and Context

The paper investigates the stability and instability properties of a specific one-dimensional magnetohydrodynamics (MHD) model introduced by Dai, Vasseur, and Zhang (2023). This model serves as a simplified proxy for the 3D ideal incompressible MHD equations, designed to capture the interplay between vortex stretching and transport mechanisms.

• The Model: The system is defined on the torus $\mathbb{T} = [-\pi, \pi]$ with variables $\omega^\pm$ (related to Elsässer variables) and parameters $a, p, q$.
  $$\partial_t \omega^\pm + a u^\pm \partial_\theta \omega^\pm = p \omega^\pm H \omega^\mp + q \omega^\mp H \omega^\pm$$
  where $H$ is the Hilbert transform and $u^\pm$ are determined by $\partial_\theta u^\pm = H \omega^\pm$.
• Connection to Known Models: In the absence of a magnetic field (or when $\omega^+ = \omega^-$), the system reduces to the De Gregorio model ($a=1, p+q=1$), which is a critical model for studying finite-time blow-up vs. global regularity in 3D Euler equations.
• Specific Focus: The authors analyze the system near the first excited state (steady state) on the torus:
  $$\Omega^\pm_2(\theta) = -\sin(2\theta), \quad U^\pm_2(\theta) = \frac{1}{2}\sin(2\theta).$$
  This extends recent work by Guo and Jiu (2025) on the De Gregorio model to the coupled MHD system.

2. Methodology

The authors employ a rigorous functional analytic approach centered on perturbation analysis and spectral theory.

A. Functional Framework

• Weighted Sobolev Space ($H^2$): The analysis is conducted in a specific Hilbert space $H^2$ inspired by previous works on the De Gregorio model. The space consists of odd functions with a weighted norm involving the singular weight $\rho_2 = (4\pi \sin^2 \theta)^{-1}$:
  $$\|\eta\|_{H^2} = \left( \int_{-\pi}^\pi \rho_2 |\partial_\theta \eta|^2 d\theta \right)^{1/2}.$$
• Orthonormal Basis: A complete orthonormal basis $\{e_{2,k}\}_{k \ge 1}$ is constructed for $H^2$, defined as:
  $$e_{2,k} = \frac{\sin((k+2)\theta)}{k+2} - \frac{\sin(k\theta)}{k}.$$
  This basis is crucial because it diagonalizes the linear operators in a way that eliminates cross-terms in the energy estimates.

B. Linearization and Operator Decomposition

The system is linearized around the excited state $\Omega_2$. The perturbation equations for $\eta^\pm$ are written as:
$$\partial_t \eta^+ = \mathcal{L}_+ \eta^+ + N_1, \quad \partial_t \eta^- = \mathcal{L}_- \eta^- + N_2.$$

• Key Insight: The linear operator $\mathcal{L}_+$ governing the $\eta^+$ component is structurally identical to the linearized operator of the De Gregorio model.
• Spectral Analysis: The authors expand the solution in the basis $\{e_{2,k}\}$. They analyze the coefficients $d_{2,k}$ arising from the action of the linear operators on the basis functions.
  • For $\mathcal{L}_-$, the coefficients are strictly negative, implying exponential decay (stability).
  • For $\mathcal{L}_+$, the coefficients change sign. Specifically, the term corresponding to $k=1$ is positive, while higher modes are negative. This sign change is the source of instability.

C. Proof Techniques

1. Galerkin Method: Used to establish global well-posedness for the linearized system and local well-posedness for the nonlinear system.
2. Second-Order ODE Analysis: To prove instability, the authors derive a second-order ordinary differential inequality for the $H^2$ norm of the perturbation. By analyzing a positive definite quadratic form associated with the Fourier coefficients, they show that the norm grows exponentially for specific initial data.
3. Contradiction Argument: For nonlinear instability, they construct a family of solutions to the nonlinear system and compare them to the linearized solution using a scaling argument, proving that the nonlinear solution cannot remain bounded if the linear part is unstable.
4. Bootstrap Argument: Used for the nonlinear stability result, showing that if initial data is sufficiently small and lies in a specific subspace (only even modes), the nonlinear terms are dominated by the linear decay.

3. Key Contributions and Results

A. Global Well-Posedness (Linear)

• Result: The linearized system is globally well-posed in the weighted Sobolev space $H^2$ for initial data with sufficient regularity ($m > 3$).
• Significance: This establishes a solid foundation for analyzing the nonlinear dynamics near the excited state.

B. Linear Instability

• Result: The linearized system is unstable. For a broad class of initial data (specifically those with a non-zero projection onto the first mode $e_{2,1}$ and satisfying a spectral condition), the $H^2$ norm of the solution grows exponentially in time.
• Mechanism: The instability is driven entirely by the $\eta^+$ component (the "plus" mode), which mimics the instability mechanism of the De Gregorio model. The $\eta^-$ component remains stable.

C. Nonlinear Instability

• Result: The full nonlinear system is nonlinearly unstable under the Lipschitz structure.
• Definition: There exist initial data arbitrarily close to the excited state such that the solution grows to a macroscopic size in finite time, violating the Lipschitz continuity of the flow map with respect to initial data.
• Method: Proven via a contradiction argument showing that if the nonlinear solution were stable, it would contradict the known exponential growth of the linearized solution.

D. Nonlinear Stability (Conditional)

• Result: The system is nonlinearly stable (globally well-posed and exponentially decaying) if the initial perturbation $\eta^+_0$ contains only even Fourier modes (i.e., $\eta^+_0 = \sum \eta_{2,2k} e_{2,2k}$).
• Significance: This identifies a specific subspace of the phase space where the destabilizing mechanism (the odd mode $k=1$) is absent, allowing the stabilizing transport effects to dominate.

4. Significance and Implications

1. Extension of De Gregorio Results: The paper successfully generalizes the recent stability/instability dichotomy found in the De Gregorio model (a single scalar equation) to a coupled MHD system. This confirms that the magnetic field coupling does not suppress the fundamental instability mechanism inherent in the vortex stretching term.
2. Mechanism of Instability: It rigorously identifies that the instability near excited states in 1D MHD models is driven by the same spectral properties as the Euler/De Gregorio models. The "plus" mode acts as the primary driver of blow-up, while the "minus" mode acts as a stable component.
3. Subspace Stability: The discovery that restricting initial data to even modes restores global stability provides a deeper understanding of the geometric constraints required to prevent singularity formation in these models.
4. Methodological Advancement: The use of the specific weighted basis $\{e_{2,k}\}$ to diagonalize the linear operator and the derivation of the second-order ODE for the energy norm offer a robust template for analyzing similar nonlocal PDEs with mixed stability properties.

5. Conclusion

The paper establishes a sharp dichotomy for the 1D MHD model near the first excited state: the system is linearly and nonlinearly unstable for generic perturbations due to the dominance of the vortex stretching term in the $\eta^+$ component, but it is globally stable for perturbations restricted to a specific subspace (even modes). These results deepen the mathematical understanding of how magnetic fields interact with fluid instabilities in simplified 1D settings, offering insights relevant to the broader question of singularity formation in 3D MHD.