Mercier--Cotsaftis and Grad--Shafranov equations for anisotropic plasma

This brief review discusses the historical development of generalizing the Grad-Shafranov equation to describe anisotropic plasmas, with a specific focus on the contributions of Mercier and Cotsaftis.

The Gist

The Big Picture: Balancing a Hot, Wobbly Balloon

Imagine you are trying to hold a giant, super-hot balloon of gas (plasma) in place using invisible magnetic ropes. If the gas is calm and behaves the same way in every direction, it's relatively easy to predict how it will sit. This is the classic Grad–Shafranov equation—the "rulebook" for how these magnetic balloons behave in standard fusion reactors (like tokamaks).

However, in modern experiments, we blast this gas with powerful beams of energy. This makes the gas "wobbly" and anisotropic. Think of it like a crowd of people in a room:

• Isotropic (Normal): Everyone is jiggling randomly in all directions.
• Anisotropic (Wobbly): Everyone is suddenly running fast in one direction but standing still in another.

When the plasma acts like this "wobbly crowd," the old rulebook breaks down. Scientists need a new, more complex rulebook to keep the balloon from popping.

The Confusion: Too Many Names for One Thing

The author of this paper, Igor Kotelnikov, is acting like a librarian who is frustrated by a messy filing system. He noticed that scientists are using four or five different names for essentially the same new rulebook:

1. Generalized Grad–Shafranov Equation
2. Modified Grad–Shafranov Equation
3. Anisotropic Grad–Shafranov Equation
4. Mercier–Cotsaftis Equation

It's like if everyone called a "car" by different names depending on who invented it: "The Ford," "The Auto," "The Horseless Carriage," and "The Smith-Mobile." This makes it hard to find information and understand the history.

The History: Who Actually Wrote the Rulebook?

The paper digs into the history to set the record straight:

• The Famous Duo (Grad & Shafranov): In the 1950s, Harold Grad and V. Shafranov wrote the famous rulebook for the normal (isotropic) plasma. They are the "celebrity authors" everyone knows.
• The Forgotten Pioneers (Mercier & Cotsaftis): In 1961, two French scientists, Claude Mercier and Michel Cotsaftis, figured out how to write the rulebook for the wobbly (anisotropic) plasma. They did the hard math first.
• The "Re-invention": Six years later, Harold Grad (the famous one) published the same math again but didn't mention Mercier and Cotsaftis. Because Grad was more famous, the scientific community started calling it the "Generalized Grad–Shafranov Equation," effectively erasing the French pioneers from the title.

The Author's Analogy:
Kotelnikov compares this to the Hubble Constant (the speed at which the universe expands). For decades, it was just called "Hubble's Law." But it turns out a Belgian priest named Georges Lemaître figured it out mathematically two years earlier. It took until 2018 for the scientific world to officially rename it the "Hubble-Lemaître Law" to give credit where it's due.

Kotelnikov argues we should do the same for plasma physics: Call it the Mercier–Cotsaftis Equation.

The "Diamagnetic Trap" Problem: A New Twist

The paper starts by discussing a recent experiment involving a "diamagnetic trap."

• The Old Way: Imagine a magnetic cage holding the gas.
• The New Way: The gas is so energetic it pushes the magnetic cage out, creating a bubble where the magnetic field is almost gone.

In this new scenario, the math gets weird. The equation used to describe it (Equation 1 in the paper) isn't a standard math problem anymore; it's an integro-differential equation.

• Simple Analogy: A normal equation is like a recipe: "If you add 2 eggs, you get a cake."
• The New Equation: It's like a recipe that says, "If you add 2 eggs, the result depends on how many eggs you added yesterday, and how big the pan is tomorrow." It's a messy, non-local relationship.

The author points out that calling this messy new equation the "Grad–Shafranov equation" is technically incorrect because it's a different type of math problem entirely. He jokes that maybe it should be called the Beklemishev–Khristo Equation instead.

The "Middle Ground" Solution

Despite the messiness of the new "diamagnetic trap" equation, the author suggests that the Mercier–Cotsaftis Equation (the one from 1961) is actually the perfect tool for describing the middle stage of this process. It's the sweet spot where the magnetic field is being pushed out but hasn't vanished completely yet.

The Takeaway

1. The Math: We have a powerful equation to describe hot, wobbly plasma, but it's complicated.
2. The Names: We are using too many confusing names for it.
3. The History: We forgot the French scientists (Mercier and Cotsaftis) who discovered the anisotropic version first.
4. The Goal: The author wants us to stop calling it the "Generalized Grad–Shafranov" and start calling it the Mercier–Cotsaftis Equation to honor the original discoverers and clear up the confusion.

In short: It's a plea for historical justice and better organization in the world of fusion energy research.

Technical Summary

1. Problem Statement

The paper addresses a terminological and conceptual confusion in plasma physics regarding the equilibrium equations for anisotropic plasmas (where parallel pressure $p_\parallel$ differs from perpendicular pressure $p_\perp$).

• The Confusion: Recent literature (e.g., Khristo and Beklemishev) refers to equations describing diamagnetic traps as the "Grad–Shafranov equation" (GSE). However, the author argues this is a misnomer because the equation used in these specific cases is integro-differential (non-local) rather than the standard nonlinear partial differential equation (PDE) that defines the canonical GSE.
• The Gap: There is a lack of clarity in the scientific community regarding the historical lineage and distinct naming of the equations that generalize the GSE to anisotropic cases. Various terms are used interchangeably or inconsistently: Generalized GSE (GGSE), Modified GSE (MGSE), Anisotropic GSE (AGSE), and Mercier–Cotsaftis Equation (MCE).
• The Goal: To clarify the historical record, distinguish between the different mathematical forms of these equations, and advocate for the restoration of the "Mercier–Cotsaftis" name to honor the original discoverers of the anisotropic equilibrium formulation.

2. Methodology

The author employs a historical review and bibliometric analysis combined with theoretical derivation:

• Bibliometric Analysis: The author conducted an approximate count of scientific articles, monographs, and textbooks mentioning various terms (GGSE, MGSE, AGSE, MCE) across leading journals (Physics of Plasmas, Journal of Plasma Physics, Nuclear Fusion, Plasma Physics and Controlled Fusion).
• Theoretical Comparison: The paper compares the mathematical structure of:
  1. The Canonical Grad–Shafranov Equation (isotropic pressure).
  2. The Generalized Grad–Shafranov Equation (anisotropic pressure, derived by Grad in 1967).
  3. The Mercier–Cotsaftis Equation (original 1961 formulation by Mercier and Cotsaftis).
• Historical Reconstruction: The author traces the publication history of these equations, noting that while Mercier and Cotsaftis published the anisotropic equilibrium equation in 1961, Harold Grad reintroduced it in 1967 without citing them, leading to the current dominance of the "Grad–Shafranov" nomenclature in the West.

3. Key Contributions

A. Clarification of Terminology and History

• Historical Priority: The paper establishes that Claude Mercier and Michel Cotsaftis (1961) were the first to derive the equilibrium equation for anisotropic plasma, predating Harold Grad's 1967 work by six years.
• Nomenclature Statistics:
  • Terms containing "Grad–Shafranov" (GGSE, MGSE, AGSE) appear in 800–1,200 articles and 15–30 monographs.
  • The term "Mercier–Cotsaftis Equation" (MCE) appears in only 40–70 articles and 5–10 monographs.
  • The author notes that MCE is often used implicitly in stability analyses (e.g., deriving the Mercier criterion) without being explicitly named.

B. Mathematical Distinctions

The author rigorously distinguishes between the different forms of the equilibrium equation:

1. Canonical GSE (Isotropic):

   • A nonlinear PDE where the right-hand side depends only on the magnetic flux $\psi$.
   • Equation: $\Delta^* \psi = -\mu_0 R^2 \frac{dp}{d\psi} - \mu_0^2 i_p \frac{di_p}{d\psi}$.

2. Generalized GSE / MCE (Anisotropic):

   • For anisotropic pressure ($p_\parallel \neq p_\perp$), the equilibrium condition involves a pressure tensor.
   • Grad's Formulation (1967): Introduces an anisotropy parameter $\sigma = 1 - \mu_0(p_\parallel - p_\perp)/B^2$ and a modified current flux $F = \sigma i_p$.
   • Mercier–Cotsaftis Formulation: An alternative form where the Shafranov operator is modified to include the anisotropy term directly:
     $$\nabla \cdot \left( \frac{\sigma}{R^2} \nabla \psi \right) = -\mu_0 \frac{\partial p_\parallel}{\partial \psi} - \mu_0^2 \frac{F}{\sigma R^2} \frac{dF}{d\psi}$$
   • Significance: This form preserves the elliptic type only if $\sigma > 0$. Violation of this condition signals firehose instability. Thus, the MCE inherently links equilibrium structure to stability criteria.

C. Critique of Recent "Diamagnetic Trap" Equations

• The paper critiques recent equations used for diamagnetic traps (e.g., Eq. 1 in the paper).
• The Issue: In diamagnetic traps, the fast ion orbit diameter is comparable to the plasma column, making the relationship between current $J_\theta$ and flux $\psi$ non-local.
• Conclusion: These equations are integro-differential, not PDEs. Calling them "Grad–Shafranov" is mathematically incorrect, as the GSE class is strictly defined as nonlinear PDEs.

4. Results

• Bibliometric Findings: The "Grad–Shafranov" label is overwhelmingly dominant in modern literature (over 10,000 references when including implicit uses), while the "Mercier–Cotsaftis" label is rare and often restricted to French schools or historical contexts.
• Mathematical Equivalence: The paper confirms that the modern "Generalized Grad–Shafranov Equation" and the "Mercier–Cotsaftis Equation" are mathematically equivalent formulations of the same physical problem (anisotropic equilibrium), differing mainly in the arrangement of terms and the definition of the operator.
• Stability Connection: The MCE form is shown to be particularly useful for analyzing local stability (Mercier criterion) and identifying instability thresholds (firehose/mirror) directly from the equilibrium structure.

5. Significance

• Historical Justice: The paper argues for restoring the name "Mercier–Cotsaftis Equation" to the anisotropic equilibrium formulation, similar to the IAU's 2018 renaming of the "Hubble Law" to the "Hubble–Lemaître Law." This acknowledges the priority of Mercier and Cotsaftis (1961) over Grad (1967).
• Conceptual Clarity: By distinguishing between the PDE-based GSE (isotropic/generalized) and the integro-differential equations used in diamagnetic traps, the author prevents mathematical confusion in future modeling efforts.
• Practical Utility: The paper highlights that while the "Generalized GSE" is a "living tool" for modern tokamak simulation, the "Mercier–Cotsaftis" form offers specific advantages for stability analysis and understanding the physical limits of plasma equilibrium (e.g., firehose instability).
• Future Recommendations: The author suggests that future literature should adopt the "Mercier–Cotsaftis" name for the anisotropic equilibrium equation to reduce confusion and honor the original discoverers, while reserving "Grad–Shafranov" strictly for the isotropic or standard generalized PDE forms.