A first-principles variational (Hamiltonian) theory by Ruiz [Phys. Plasmas 31, 122701 (2024)] describes an imploding spherical thin shell undergoing Rayleigh–Taylor instability (RTI), dynamically coupled to an adiabatically compressed interior hot spot, and yields quasilinear predictions of the degradation of inertial-confinement-fusion (ICF) performance metrics - stagnation pressure, residual kinetic energy, and areal density - as functions of initial axisymmetric perturbations and the convergence ratio. In its "Conclusions and Future Work" section, the paper proposes three extensions of this reduced theory: generalization to fully three-dimensional RTI eigenmodes (restoring invariance under rotations of the initial perturbation direction), inclusion of the acceleration phase with mass ablation (rocket-model physics and drive-asymmetry seeding of the later RTI), and replacement of the adiabatic hot-spot closure with one containing alpha heating, radiative losses, and ablated-mass exchange so that RTI can be connected to ignition thresholds and burn initiation. These three extensions have not been carried out, leaving the spherical ICF problem open.

Background

In inertial confinement fusion (ICF), a spherical fuel capsule is rapidly compressed by laser light or by X-rays generated in a surrounding hohlraum (a high-atomic-number cavity that converts laser light into X-rays). The imploding shell stagnates against the fuel it has compressed into a central hot spot, and ignition requires the hot spot to self-heat through alpha-particle deposition faster than it loses energy. The Rayleigh–Taylor instability (RTI) - the growth of perturbations at an interface between fluids of different density that is accelerated in the direction of the heavier fluid - develops at the shell's inner surface during deceleration and is recognized as a principal limiter of ICF performance. Its impact is quantified through metrics such as the stagnation pressure (peak hot-spot pressure), the areal density (column surface density, ∫ρ dr\int\rho\,dr), the residual kinetic energy (kinetic energy still present at stagnation in flows that do no compressive work), and the convergence ratio (initial-to-final hot-spot radius). Reduced analytical models matter here because fully three-dimensional (3D) radiation-hydrodynamic simulations of imploded capsules are computationally expensive.

Ruiz [Phys. Plasmas 31, 122701 (2024)] developed a first-principles variational theory of this problem: the dynamics of an imploding spherical thin shell (all shell mass concentrated on a surface) undergoing RTI is derived from an action principle, giving a Hamiltonian model in which the shell is dynamically coupled to an adiabatically compressed fluid in its interior. Perturbations are restricted to two-dimensional (2D) axisymmetric shapes expanded in Legendre modes Pℓ(cos⁡ϑ)P_\ell(\cos\vartheta) (axisymmetric angular shapes on the sphere labeled by the mode number ℓ\ell). A quasilinear analysis - a weakly nonlinear expansion in the perturbation amplitude keeping the leading correction to volume-averaged quantities - gives explicit degradation trends of the performance metrics as functions of the initial amplitude and mode number and of the convergence ratio, and these trends are confirmed qualitatively by symplectic numerical integration (a structure-preserving time-stepping method for Hamiltonian systems) of the fully nonlinear reduced equations. The model recovers, within one framework, effects previously studied separately: residual-kinetic-energy yield degradation [Woo et al., Phys. Plasmas 25, 052704 (2018)] and early pseudo-3D thin-shell stability calculations [Taguchi and Mima, Phys. Plasmas 2, 2790 (1995)]. Results on ignition-threshold sensitivity to areal-density asymmetries [Woo and Betti, Phys. Plasmas 28, 054503 (2021); Hurricane et al., Phys. Plasmas 28, 022704 (2021)] are not recovered by the adiabatic model; the source cites them as related prior work that the proposed non-adiabatic hot-spot extension would build upon.

The model's acknowledged limitations define the open target. It describes only the deceleration phase: the drive/acceleration stage and mass ablation at the shell surface (rocket-model physics, in which blowing-off surface material accelerates the shell like a rocket) are absent, so drive-asymmetry seeding of the later RTI cannot be followed. The hot spot is adiabatic: no alpha heating (energy deposition by fusion-born helium nuclei), no bremsstrahlung losses (radiation emitted by electrons decelerated in the field of ions), and no ablated mass entering the hot spot, so ignition cannot be addressed. Perturbations are 2D axisymmetric, and the model is not invariant under rotations of the initial perturbation direction - an invariance that a 3D description should restore. The 2D model predicts different degradation for prolate and oblate implosions (capsule shapes elongated versus flattened along the symmetry axis) with equal-amplitude perturbations; because a sign flip of the initial amplitude is not a rotation, part of that prolate/oblate difference is a physical geometric effect (a sign-dependent higher-order volume term) that persists in a correct 3D description, and only the component arising from the assumed axisymmetry is an artifact that 3D should remove. In addition, the thin-shell closure overestimates the compression work and the confinement time, and the quasilinear predictions break down at large amplitudes where spike structures (narrow fingers of dense shell material penetrating the lighter fluid) form and Lagrangian nodes (material points that move with the shell) coalesce on the symmetry axis.

Problem Statement

The problem, as posed in the "Conclusions and Future Work" section of Ruiz [Phys. Plasmas 31, 122701 (2024)], is to extend the Hamiltonian thin-shell theory of RTI-driven performance degradation in ICF implosions along the three directions proposed there:

  1. Three-dimensional RTI asymmetries. Generalize the model from axisymmetric Legendre modes to fully 3D perturbations of the imploding shell, and determine how the 3D RTI eigenmodes modify the modal growth, the partition of energy into residual kinetic energy, and the degradation of the stagnation conditions (pressure, areal density). The 3D model should restore invariance under rotations of the initial perturbation direction, removing the component of the 2D theory's prolate/oblate difference that arises from its assumed axisymmetry; invariance with respect to the sign of the initial perturbation amplitude should not be expected, since a sign flip is not a rotation. The 3D model is intended to provide theoretical insight at a cost far below 3D radiation-hydrodynamic simulation.

  2. Acceleration phase and mass ablation. Extend the model to the drive/acceleration stage of the implosion, incorporating mass ablation of the shell surface as in the rocket model, so that shell perturbations seeded by asymmetries in the radiation drive can be followed through their acceleration-phase (RTI) growth and their subsequent impact on the stagnation phase.

  3. Non-adiabatic hot-spot dynamics and ignition. Replace the adiabatic interior closure with a hot-spot model that includes internal heat sources (alpha heating), energy losses (bremsstrahlung x-ray emission), and mass ablation from the cold shell into the hot spot - in the spirit of the thermodynamic ignition condition of Hurricane et al. (2021) - and use it to establish how RTI-induced asymmetries shift the initiation of ignition and the subsequent burn, i.e., the ignition thresholds of the coupled shell–hot-spot system.

The extended theory should remain consistent with the published 2D adiabatic limit of the source theory and retain its reduced-order character, yielding analytic or controlled reduced-order predictions of the degradation of the ICF performance metrics and, for the non-adiabatic hot-spot extension, of ignition-threshold shifts as functions of the initial asymmetry parameters and implosion characteristics (e.g., the convergence ratio).

The three directions are independent extensions: a self-contained solution to any one of them constitutes a resolving contribution, while addressing all three constitutes the complete objective. Completing the first direction requires covering all three quantities listed there - the modal growth, the partition of energy into residual kinetic energy, and the degradation of the stagnation conditions.

References

  1. D. E. Ruiz, "Degradation of performance in ICF implosions due to Rayleigh–Taylor instabilities: A Hamiltonian perspective", Physics of Plasmas 31, 122701 (2024). DOI: 10.1063/5.0233454; arXiv:2411.12933. URL: https://doi.org/10.1063/5.0233454 (arXiv: https://arxiv.org/abs/2411.12933)
  2. D. E. Ruiz and M. Simmons, "Nonlinear analysis of the degradation effects on Z-pinch implosions due to azimuthal magneto-Rayleigh–Taylor instabilities", Physics of Plasmas 32, 082701 (2025). DOI: 10.1063/5.0266167. URL: https://doi.org/10.1063/5.0266167
  3. D. E. Ruiz, C. A. Williams, and R. A. Vesey, "Asymptotic scaling laws for the stagnation conditions of Z-pinch implosions", Physics of Plasmas 32, 012707 (2025). DOI: 10.1063/5.0242367; arXiv:2501.12284. URL: https://doi.org/10.1063/5.0242367 (arXiv: https://arxiv.org/abs/2501.12284)
  4. X. Li, Y. Wu, Z. Chen, M. Yang, and J. Zhang, "Bell Plesset Effects on Rayleigh Taylor Instability of Three Dimensional Spherical Geometry" (2025). arXiv:2511.19028. URL: https://arxiv.org/abs/2511.19028
  5. D. Liu, T. Tao, J. Li, Q. Jia, R. Yan, and J. Zheng, "A one-dimensional mixing model for the impact of ablative Rayleigh-Taylor instability on compression dynamics" (2024, revised 2025). arXiv:2412.00699. URL: https://arxiv.org/abs/2412.00699
  6. K. M. Woo and R. Betti, "Impact of areal-density asymmetries on the loss of confinement and ignition threshold in inertial confinement fusion capsules", Physics of Plasmas 28, 054503 (2021). DOI: 10.1063/5.0045317. URL: https://doi.org/10.1063/5.0045317
  7. T. Taguchi and K. Mima, "Pseudo-three-dimensional simulation on stability of thin shell target implosion", Physics of Plasmas 2, 2790 (1995). DOI: 10.1063/1.871177. URL: https://doi.org/10.1063/1.871177
  8. K. M. Woo, R. Betti, D. Shvarts, A. Bose, D. Patel, R. Yan, P.-Y. Chang, O. M. Mannion, R. Epstein, J. A. Delettrez, M. Charissis, K. S. Anderson, P. B. Radha, A. Shvydky, I. V. Igumenshchev, V. Gopalaswamy, A. R. Christopherson, J. Sanz, and H. Aluie, "Effects of residual kinetic energy on yield degradation and ion temperature asymmetries in inertial confinement fusion implosions", Physics of Plasmas 25, 052704 (2018). DOI: 10.1063/1.5026706. URL: https://doi.org/10.1063/1.5026706
  9. O. A. Hurricane, S. A. Maclaren, M. D. Rosen, J. H. Hammer, P. T. Springer, and R. Betti, "A thermodynamic condition for ignition and burn-propagation in cryogenic layer inertially confined fusion implosions", Physics of Plasmas 28, 022704 (2021). DOI: 10.1063/5.0035583. URL: https://doi.org/10.1063/5.0035583