Relativistic calculations of electron impact excitation cross-sections of neutral tungsten

This paper presents the first comprehensive set of relativistic distorted-wave calculations for fine-structure-resolved electron impact excitation cross-sections and radiative transition probabilities of neutral tungsten (WI), emphasizing the critical role of metastable states in collisional-radiative modeling for plasma diagnostics.

The Gist

Imagine a giant, glowing pot of soup called a fusion reactor. Scientists are trying to cook this soup to create clean, limitless energy. To keep the soup from boiling over or getting too hot, they line the pot's walls with a special metal called Tungsten. It's like the ultimate non-stick pan: it can handle extreme heat without melting and doesn't react badly with the soup.

However, sometimes tiny bits of this Tungsten get knocked off the wall and float into the soup. This is a problem because Tungsten is heavy and acts like a "heat sink," sucking the energy out of the soup and cooling it down too fast. If the soup cools too much, the fusion reaction stops.

To fix this, scientists need to know exactly how much Tungsten is in the soup and where it is. They do this by looking at the light the Tungsten emits. Think of it like a barcode. Every time a Tungsten atom gets bumped by a fast-moving electron (a tiny particle zooming through the soup), it jumps to a higher energy level and then falls back down, flashing a specific color of light. By reading these colors, scientists can diagnose the health of the reactor.

The Problem:
To read this "barcode" correctly, scientists need a perfect instruction manual (a database) that tells them:

1. Exactly how much energy it takes to bump a Tungsten atom.
2. Exactly what color of light it will flash when it falls back down.

The problem is that Tungsten is a very complicated atom (it has 74 electrons, which is a lot!). It's like a crowded dance floor where everyone is bumping into each other. Previous manuals were incomplete or had errors, especially for the "lazy" Tungsten atoms that aren't in their most energetic state (called metastable states). Ignoring these lazy atoms is like trying to predict traffic by only looking at cars moving fast, ignoring the ones stuck in traffic jams.

What This Paper Does:
The authors of this paper are like super-accurate cartographers drawing a new, detailed map of the Tungsten atom's behavior.

1. The Method (The RDW Tool): They used a powerful computer program (called the Flexible Atomic Code) that uses a method called "Relativistic Distorted Wave."

   • Analogy: Imagine trying to predict how a ball bounces off a trampoline. If the trampoline is flat and simple, it's easy. But Tungsten's trampoline is wobbly, heavy, and distorted by gravity (relativity). This paper uses a super-sophisticated simulation to figure out exactly how the ball bounces off this messy, heavy trampoline.

2. The Discovery (The Metastable Secret): They didn't just look at the Tungsten atoms in their "ground state" (the most common, calm state). They also looked at the metastable states.

   • Analogy: Think of a crowd of people. Most are standing still (ground state). But some are pacing back and forth, waiting for a cue (metastable). The paper found that when you bump into these "pacing" people, they react much more violently and brightly than the calm ones.
   • Key Finding: They discovered that for many important light flashes, the "pacing" (metastable) atoms are actually the main source of the signal, not the calm ones. If you ignore them, your diagnosis of the reactor is wrong.

3. The Results (The New Map):

   • They calculated the "bump energy" and "light flash" for thousands of different scenarios, covering a wide range of speeds (energies) that electrons might have in a reactor.
   • They checked their work against the gold standard (NIST database) and found their map is very accurate, especially for the tricky, low-energy areas where previous maps were fuzzy.
   • They also calculated how likely these atoms are to flash light (radiative rates) and compared them to other experts' data, noting where the complexity of Tungsten makes things tricky.

Why This Matters:
This paper provides a highly accurate, updated instruction manual for Tungsten.

• For Fusion Reactors: It helps scientists build better computer models to predict how much Tungsten is eroding from the walls. This allows them to adjust the reactor to keep the "soup" hot and stable.
• For Diagnostics: It ensures that when they look at the light from the reactor, they aren't misinterpreting the signal. They can now tell, with much higher confidence, exactly how much Tungsten is floating around and where it came from.

In a Nutshell:
This paper is a massive upgrade to the "Tungsten Dictionary." It tells us that to understand the light from a fusion reactor, we can't just look at the calm atoms; we have to pay attention to the "pacing" ones too. With this new, detailed map, scientists can better manage the heat and safety of future clean-energy power plants.

Technical Summary

Here is a detailed technical summary of the paper "Relativistic calculations of electron impact excitation cross-sections of neutral tungsten".

1. Problem Statement

Neutral tungsten (W I) is a critical component in the edge and divertor regions of magnetic confinement fusion devices (such as ITER). Sputtered tungsten enters the plasma as neutral atoms, and their visible/near-UV emission lines are primary diagnostics for measuring tungsten influx and erosion. However, quantitative interpretation of these spectra via Collisional-Radiative Models (CRMs) requires accurate atomic data, specifically:

• Fine-structure resolved energy levels and radiative rates.
• Electron-impact excitation (EIE) cross-sections connecting ground and metastable states to excited configurations.

Challenges:

• Atomic Complexity: Tungsten (Z=74) has an open 5d shell, leading to strong relativistic effects, significant spin-orbit mixing, and dense configuration interaction (CI).
• Metastable Populations: Under edge plasma conditions, long-lived metastable states accumulate significant populations. Excitation from these metastables can dominate over ground-state excitation, yet many existing datasets focus only on the ground state.
• Data Gaps: Existing W I collision datasets are often sparse, limited to specific lines, or restricted to low energy ranges, lacking the comprehensive coverage needed for predictive modeling.

2. Methodology

The authors employed a two-step theoretical approach using the Flexible Atomic Code (FAC):

A. Atomic Structure Calculation (Target Description)

• Method: Multi-configurational Dirac-Fock (MCDF) approach.
• Configuration Interaction (CI): An extensive CI expansion was used, including 43,556 fine-structure resolved energy levels. The configurations included valence-valence and core-valence correlations (involving the 5p shell and deeper cores).
• Benchmarking: Calculated energy levels were compared against NIST recommended values and other theoretical works (HFR, GRASP, MDFGME).
• Correction: To improve accuracy, the calculated energy levels were corrected to match NIST values before calculating cross-sections.

B. Scattering Calculation (EIE Cross-sections)

• Method: Relativistic Distorted-Wave (RDW) method.
• Scope: Calculations were performed for excitations from the ground state ($5d^4 6s^2 \ ^5D_0$) and **six metastable states** ($^5D_{1-4}$,$^3P_0$, and$^7S_3$) to excited levels primarily belonging to the$5d^4 6s(6D)6p$and$5d^5(6S)6p$ configurations.
• Energy Range: Incident electron energies from the excitation threshold up to 500 eV.
• Low-Energy Correction: Since RDW is perturbative and less accurate near threshold, a correction factor was applied to the cross-sections for low incident energies.
• Radiative Data: Radiative transition probabilities (A-values) were calculated for selected prominent transitions.

3. Key Contributions

• Comprehensive Metastable Coverage: Unlike previous studies focusing mainly on the ground state, this work provides a complete dataset for excitations from seven initial states (1 ground + 6 metastables), addressing a critical gap in CRM inputs.
• Fine-Structure Resolution: The study reports cross-sections for specific fine-structure levels rather than just term-averaged values, which is essential for accurate spectral modeling.
• Extensive Configuration Set: The use of a massive CI expansion (43,556 levels) ensures that strong configuration mixing and relativistic effects in neutral tungsten are properly captured.
• New Data Generation: Most of the reported fine-structure-resolved cross-sections are presented for the first time, covering a broad energy range (0–500 eV) suitable for edge and laboratory plasmas.

4. Key Results

• Energy Levels: The calculated energy levels showed good agreement with NIST data for low-lying levels, though discrepancies increased for highly excited configurations. The MCDF+CI model successfully reproduced the key low-lying structure.
• Cross-Section Magnitudes:
  • Metastable Dominance: Excitation from the **$5d^5(6S)6s \ ^7S_3$** metastable state yielded the **largest cross-sections** (reaching$\sim 7 \times 10^{-16} \text{ cm}^2$ near threshold) for many channels, significantly exceeding those from the ground state.
  • Quintet Metastables: Excitation from the $^5D_{1-4}$ metastables also produced large cross-sections (order of $10^{-17} \text{ cm}^2$), comparable to or exceeding ground-state excitation for specific channels.
  • Selection Rules: The results confirmed that electric-dipole allowed transitions (opposite parity, $\Delta J = 0, \pm 1$) dominate. However, strong spin-orbit mixing in tungsten activates "forbidden" channels (e.g., spin-flips), redistributing excitation strength.
  • Statistical Dependence: For the $^7S_3$ initial state, cross-sections showed a clear $(2J_f + 1)$ statistical dependence, with higher-$J$ final states having larger cross-sections.
• Radiative Rates: Calculated A-values generally agreed with NIST and other theoretical data within the same order of magnitude, though deviations were noted for complex transitions due to strong CI effects.
• Comparison with R-Matrix: Comparisons with R-matrix calculations (limited to 30 eV) showed agreement for some transitions (e.g., 522.47 nm) but deviations for others (e.g., 400.88 nm) near the threshold, attributed to the perturbative nature of the RDW method missing resonance structures.

5. Significance

• Improved Diagnostics: This dataset allows for more accurate interpretation of tungsten spectroscopy in fusion devices. By including metastable excitation, models can better predict line intensities and inferred tungsten influx rates (using S/XB concepts).
• Modeling Accuracy: The data highlights that ignoring metastable populations leads to systematic errors in predicting plasma emissivities. The $^7S_3$ state, in particular, acts as a major "feeder" for excited levels.
• Resource for Fusion Research: The provided dataset (energies, A-values, and EIE cross-sections) is internally consistent and ready for implementation in CRMs, supporting the operation of current devices and future reactors like ITER where tungsten is the primary plasma-facing material.
• Future Work: The authors plan to integrate this atomic data into full collisional-radiative models to simulate spectral line intensities across a wide range of plasma conditions, further validating the impact of metastable populations.