Characterisation of X- and O-points in Wendelstein 7-X with respect to coil currents

This paper presents a fast automated scheme to characterize X- and O-points in Wendelstein 7-X by analyzing how variations in superconducting coil currents influence magnetic field topology, island chain properties, and the stability metric Tr(M) across standard, high/low iota, and randomly sampled configurations.

The Gist

Imagine a Stellarator (like the Wendelstein 7-X, or W7-X) as a giant, futuristic donut-shaped oven designed to cook nuclear fusion fuel. The challenge? The fuel is a super-hot gas called plasma, and it wants to escape. If it touches the walls, the oven breaks. To keep it contained, scientists use powerful magnets to create an invisible "cage" of magnetic fields that twists and turns the plasma into a stable shape.

However, this magnetic cage isn't perfect. Sometimes, the lines of the magnetic field get tangled and form little loops or "islands" in the edge of the plasma. Some of these islands are bad (they let heat escape), while others are actually helpful (they act like a drain to remove waste heat).

This paper is like a master mechanic's manual for tweaking the knobs on this magnetic oven. The authors wanted to understand exactly how changing the electricity in the different magnets changes the shape of these magnetic islands.

Here is the breakdown using simple analogies:

1. The Three Types of "Knobs" (Coils)

The W7-X machine has three main types of magnets (coils) that act as the controls. Think of them as different tools in a toolbox:

• The "Shape-Shifters" (Non-Planar Coils): These are the complex, twisted magnets that give the plasma its unique, twisted donut shape.

  • Analogy: Imagine these are the architects who draw the blueprints. They decide the overall shape of the room.
  • What the paper found: Once the room is built, turning these knobs a little bit doesn't move the furniture much. They are great for setting the stage, but not for fine-tuning the details.

• The "Shifters" (Planar Coils): These are flat, simple magnets.

  • Analogy: These are like moving trucks. They can push the entire magnetic cage inward or outward, or tilt it up and down.
  • What the paper found: These are the most powerful tools for deciding where the magnetic islands appear. If you want the "drain" (the island) to hit a specific spot on the wall, you use these knobs to slide the whole pattern left, right, up, or down.

• The "Tuner" (Control Coil): This is a small, special magnet close to the plasma.

  • Analogy: This is like the fine-tuning knob on a radio or the volume dial. It doesn't move the station or the room; it just adjusts the clarity and size of the signal.
  • What the paper found: This knob is surprisingly powerful. It can change the size of the magnetic islands and even flip their nature. It can turn a "bad" island into a "good" one, or vice versa, without moving the whole pattern.

2. The "Fixed Points" (X and O Points)

To understand the islands, the scientists looked at specific spots called Fixed Points.

• O-Points: Think of these as the center of a whirlpool. The magnetic field lines swirl around this point. It's a calm, stable spot.
• X-Points: Think of these as the crossroads or the "saddle" between two hills. This is where the magnetic field lines cross and split.
  • Why it matters: The X-point is the "doorway" where heat and particles escape the main plasma and get directed to the divertor (the exhaust system).

3. The "Magic Number" (Tr(M))

The authors developed a fast computer program to calculate a specific number for every X-point and O-point. Let's call this the "Stability Score."

• If the score is 2, the island is perfectly balanced (or non-existent).
• If the score is above 2, it's an X-point (the doorway).
• If the score is below 2, it's an O-point (the whirlpool center).

The paper discovered a fascinating trick: The Control Coil can change the Stability Score of the X-point and the O-point differently.

• The Analogy: Imagine a seesaw. Usually, if you push one side down, the other goes up. But here, the Control Coil can push the X-point down (turning it into an O-point) while leaving the O-point almost untouched. This allows scientists to change the type of island chain without moving its location.

4. The "Greene's Residue" (The Size Predictor)

One of the coolest findings is a shortcut for measuring how big an island is.

• The Problem: Measuring the exact size of a magnetic island in a computer simulation is slow and hard.
• The Solution: The authors found that the "Stability Score" (specifically how far it is from 2) is a perfect proxy for the island's size.
• The Analogy: It's like judging the size of a storm by looking at the barometer. You don't need to measure the wind speed directly; if the pressure drops a certain amount, you know the storm is big.
• Why this helps: Now, scientists can quickly scan thousands of magnetic configurations and instantly know which ones will create islands of the perfect size to act as an exhaust system, without doing hours of complex calculations.

Summary: Why does this matter?

Fusion reactors need to get rid of "ash" (helium waste) and excess heat. The W7-X uses a clever system called an Island Divertor, which uses these magnetic islands to channel waste out of the machine.

This paper gives engineers a cheat sheet. It tells them:

1. Use the Planar Coils to move the island to the right spot.
2. Use the Control Coil to adjust the island's size and type.
3. Use the "Stability Score" to quickly guess if the island is the right size.

By understanding these rules, scientists can design better fusion reactors that don't overheat and can run for longer, bringing us closer to clean, limitless energy.

Technical Summary

1. Problem Statement

Commercially viable stellarator reactors require robust exhaust solutions to manage heat and particle flow (e.g., removing helium ash and preventing thermal overload of plasma-facing components). The Wendelstein 7-X (W7-X) utilizes an island divertor concept, which relies on magnetic island chains with low-order rational rotational transform ($\iota$) to divert plasma.

While the general roles of W7-X's electromagnetic coils are known, there is a lack of quantitative understanding regarding how specific variations in coil currents affect the topology of the edge magnetic field, specifically the location, stability, and size of fixed points (X-points and O-points) that define these island chains. The paper aims to quantify these relationships across the operational space of W7-X to better predict and control edge topology.

2. Methodology

The authors developed a fast, automated scheme to analyze vacuum magnetic field topology without relying on the calculation of unperturbed nested flux surfaces (a common but non-unique and computationally expensive approach in resonant perturbation models).

• Fixed Point Detection: The scheme locates fixed points of the field line map (points where field lines return to their starting position after $n$ field periods). It specifically targets periodicities $n=4, 5, 6$, corresponding to the 5/4, 5/5, and 5/6 island chains.
• Jacobian Analysis: For each fixed point, the code calculates the trace of the Jacobian matrix ($Tr(M)$) of the field line map.
  • $Tr(M) > 2$: X-point (hyperbolic, unstable).
  • $-2 < Tr(M) < 2$: O-point (elliptic, stable).
  • $Tr(M) = 2$: Rational surface (intact).
  • The quantity $|Tr(M) - 2|$ is used as a proxy for island size and relates to Greene's residue.
• Simulation Tools: The field line tracing is performed using EMC3-Lite (a Biot-Savart solver with an Adams-Bashforth-Moulton integrator), interfaced via a Python wrapper.
• Scan Strategy: Two sets of scans were performed:
  1. 1D Scans: Individual coil currents were varied around three paradigmatic configurations: "Standard" (5/5), "High $\iota$" (5/4), and "Low $\iota$" (5/6).
  2. Large Random Scan: Over $2 \times 10^5$ magnetic configurations were generated by randomly sampling coil currents within the normal operational range of W7-X.

3. Key Contributions

• Quantification of Coil Roles: The study provides a quantitative map of how Non-Planar Coils (NPCs), Planar Coils (PCs), and Control Coils (CC) specifically influence fixed point properties (location, $Tr(M)$, and island size).
• Asymmetry in Control Coil Effects: The paper identifies a previously under-quantified phenomenon where the control coil affects X-points and O-points within the same island chain unequally, potentially causing topological transitions (e.g., X-point to O-point) while leaving the other point largely unchanged.
• Proxy for Island Size: The authors establish a strong correlation between $|Tr(M) - 2|$ and the physical area of internal island chains, offering a computationally efficient metric for identifying experimental candidates with specific island widths.

4. Key Results

A. Roles of Specific Coil Types

• Non-Planar Coils (NPCs): Primarily responsible for plasma shaping and generating the rotational transform ($\iota$). They have a relatively minor effect on the radial location of the magnetic axis and edge fixed points compared to other coils.
• Planar Coils (PCs): These are the primary drivers for the radial location of the magnetic axis and fixed points.
  • The difference in current between Planar Coils A and B ($I_{PC,A} - I_{PC,B}$) controls the inward/outward shift of the plasma.
  • The sum of their currents ($I_{PC,A} + I_{PC,B}$) shifts the $\iota$ profile up or down, thereby controlling the periodicity of the island chains.
• Control Coil (CC):
  • Has negligible effect on the magnetic axis location or $\iota$.
  • Has the strongest influence on the stability ($Tr(M)$) of edge fixed points.
  • Can induce X-point to O-point transitions (changing the island chain topology, e.g., 5/5 to 10/10) by altering the phase and size of the islands.

B. Dependence on Configuration and Radius

• Sensitivity to Radius: The sensitivity of fixed points to coil currents (particularly the Control Coil) increases with the minor radius of the fixed point. Fixed points closer to the coils (as seen in the "High $\iota$" configuration) exhibit larger variations in $Tr(M)$ and location compared to those in the "Low $\iota$" configuration.
• Asymmetric Response: In the Standard configuration, the Control Coil can shift the X-point of a 5/5 chain to an O-point (creating a 10/10 chain) while the O-point remains largely unaffected. This asymmetry is less pronounced in other configurations.

C. Large Scan Statistics

• Phase Space Distribution: The distribution of fixed points in the phase space of Planar Coil currents forms a "diamond" shape.
  • Internal Islands: Mostly $n=5$ and $n=6$ periodicities.
  • Divertor-Relevant Islands: All periodicities ($n=4, 5, 6$) are possible.
  • Near-Coil Islands: Mostly $n=4$ and $n=5$.
• Tr(M) Trends:
  • $|Tr(M) - 2|$ generally increases with the minor radius of the fixed point.
  • For $n=5$ points, $Tr(M)$ increases with Control Coil current in the "divertor-relevant" range.
  • For $n=4$ points, the trend is opposite to $n=5$.
  • $n=6$ points show little dependence on Control Coil current.

D. Island Size Correlation

• A positive correlation was found between the island area and $|Tr(M) - 2|$ for internal island chains.
• This validates $|Tr(M) - 2|$ as a proxy for island size, allowing researchers to optimize for specific island widths without calculating the full geometric area for every configuration.

5. Significance

This work bridges the gap between theoretical coil control and practical magnetic topology in stellarators.

1. Experimental Optimization: The findings allow operators to predict how specific coil current adjustments will alter the edge topology, aiding in the selection of configurations that optimize exhaust (divertor performance) or suppress MHD activity.
2. Internal Island Control: By linking $|Tr(M) - 2|$ to island size, the study provides a tool to intentionally tailor internal island chains. Recent experiments suggest these chains can suppress turbulence and increase energy confinement; this method offers a way to target them.
3. Efficiency: The automated scheme is significantly faster than traditional methods requiring the reconstruction of unperturbed flux surfaces, enabling rapid scanning of the vast operational space of W7-X.

In conclusion, the paper confirms the distinct roles of W7-X's coil systems and provides a robust, quantitative framework for manipulating edge magnetic topology through precise coil current control.