Unified Model of Heated Plasma Expansion

Motivated by the need to predict plasma density and temperature distributions created in the early stages of high-intensity laser-plasma interactions, we develop a fluid model of plasma expansion into vacuum that incorporates external heating. We propose a new three-parameter family of self-similar solutions for plasma expansion that models a wide range of spatiotemporal variations of the electron temperature. Depending on the relative scales of the heated plasma domain $L$, the Debye length $λ_D$ and an emergent ion-acoustic correlation length $λ_s$, characterized by the parameters $\frac{λ_s}{λ_D}$ and $\frac{L}{λ_s}$, a spectrum of dynamical behaviors for the expanding plasma are identified. The behavior is classified into five dynamical regimes, ranging from nearly quasineutral expansion to the formation of bare ion slabs susceptible to Coulomb explosion. Self-similar solutions with spatially uniform electron temperature with exponential time dependence are analyzed, and the dynamics in the five asymptotic limits in the parameter space are detailed. Scaling relations for the length scales and energies of the expanding plasma are proposed. The self-similar framework is applied to laser-plasma interactions, specifically addressing the plasma dynamics at a target surface as the target interacts with the rising laser intensity envelope. The results offer insights into the expansion behavior based on the laser-plasma parameters, and scaling relations for optimizing laser-plasma schemes and guiding experimental designs in high-intensity laser experiments.

💡 Research Summary

The paper presents a comprehensive fluid‑based theory for the expansion of a laser‑heated plasma into vacuum, explicitly incorporating external electron heating. Starting from the collisionless two‑fluid (electron and ion) equations coupled through Poisson’s equation, the authors identify three fundamental length scales that govern the dynamics: the overall plasma size L, the electron Debye length λ_D, and an ion‑acoustic correlation length λ_s = C_s/γ, where C_s = √(Z T_e/m_i) is the ion‑sound speed and γ is the characteristic expansion rate. By forming the dimensionless ratios λ_s/λ_D and L/λ_s, the authors map the entire parameter space onto a two‑dimensional diagram that contains five distinct dynamical regimes:

1. Quasi‑neutral expansion (λ_s ≪ λ_D ≪ L) – electrons and ions remain tightly coupled; the electric field is confined to a Debye‑scale sheath at the vacuum interface.
2. Coulomb‑explosion precursor (λ_D ≫ L ≫ λ_s) – electrons outrun the ions, creating a highly charged ion slab that converts electrostatic potential energy into ion kinetic energy.
3. Ablation‑like regime (λ_s ≈ L ≫ λ_D) – the ion‑acoustic scale matches the plasma size, leading to rapid ion acceleration and a thin, high‑velocity ion front.
4. Expanding hot‑electron cloud (λ_s ≫ L ≫ λ_D) – a hot electron layer expands while ions remain essentially stationary; the electron cloud dominates the field structure.
5. Mixed regime (λ_s ≈ λ_D ≈ L) – all three scales are comparable, producing a hybrid of charge‑separation and quasi‑neutral behavior.

To capture these behaviors analytically, the authors employ a self‑similar Ansatz. They introduce a time‑dependent global length X(t) that grows much larger than its initial value, define the similarity variable ξ = x/X(t), and assume power‑law scalings for density, velocity, temperature, and electric field:

• n_α = n_α0 (X₀/X)^m N_α(ξ)
• v_α = Ẋ V_α(ξ)
• T_e = (m_i X₀² γ²/Z) (X₀/X)^{m‑2} Θ(ξ)
• eE = (m_i X₀ γ²/Z) (X₀/X)^{m‑1} E(ξ).

The requirement that λ_D, λ_s, and L remain proportional during the evolution forces the scaling exponent m to be linked to the heating law (e.g., exponential temperature rise). This yields “incomplete” self‑similar solutions because the intrinsic Debye length introduces an explicit dependence on X(t). Nevertheless, the solutions are universal attractors: they become independent of the detailed initial density or temperature profiles once X(t)≫X₀.

From the self‑similar equations the authors derive explicit scaling relations for each regime. For instance, the expansion rate obeys γ ∝ T_e^{1/2} (λ_s/λ_D)^{‑1/2}, showing that stronger heating (larger λ_s) accelerates the plasma, while a larger Debye length (greater charge separation) slows it. The ion front velocity can be sub‑sonic or supersonic depending on the position in the (λ_s/λ_D, L/λ_s) plane, and shock‑like structures appear near the leading edge in the Coulomb‑explosion limit.

Energy partition is analyzed by integrating the fluid equations: the externally supplied electron heating power Q is divided among electron thermal energy, electrostatic field energy, and ion kinetic energy. In the quasi‑neutral regime most of the energy stays as electron heat, whereas in the Coulomb‑explosion regime a substantial fraction is converted into ion kinetic energy. The authors provide analytic expressions for these fractions as functions of the two dimensionless parameters, offering a tool for optimizing ion‑acceleration efficiency.

Finally, the model is applied to a realistic laser‑target interaction where the laser intensity rises with time. By prescribing a heating function Q(t) that follows the laser envelope, the self‑similar solution predicts the formation of an electron sheath, its growth, and the subsequent ion acceleration. The analysis predicts the timing of sheath formation, the peak electric field, and the maximum ion energy as explicit functions of laser intensity, pulse duration, and target material properties (initial density, ion charge state). This framework therefore bridges the gap between early‑time plasma formation and the later high‑energy ion acceleration stage, providing a unified description that can guide experimental design for laser‑driven ion sources, high‑harmonic generation, and relativistic transparency studies.

In summary, the paper introduces a unified, self‑similar theory that captures five distinct plasma‑expansion regimes, derives scaling laws for length, velocity, and energy, and demonstrates how these results can be directly employed to predict and optimize the early dynamics of laser‑heated plasmas in high‑intensity experiments.