Whether ion–ion correlations can materially raise dense-plasma opacities — and thereby contribute to resolving the solar opacity problem — is disputed: calculations with the Rozsnyai-type ion-correlation model give enhancements of the Rosseland mean opacity of roughly 10% for the solar mixture near the convection-zone boundary (1.5–2% at the solar center) and roughly 15% for iron at Z-facility conditions — with the STAR opacity values quoted by Ovechkin et al. corresponding to a 20% increase for iron at Z-machine conditions — while the Starrett–Saumon average-atom framework, which treats electronic and ionic structure self-consistently, gives only a few percent but suffers from an unresolved defect: there is no physically consistent prescription for the asymptotic free-electron chemical potential μe\mu_e, because the natural variational condition derived from a first-order cluster expansion violates the Debye–Hückel limit. The open problem is to construct a self-consistent prescription for μe\mu_e within this framework (or an improved cluster/variational treatment accounting for the dependence of the pair-correlation function g(r)g(r) on the asymptotic free-electron density), verify its Debye–Hückel limit, and remove the resulting ambiguity in the mean ion charge Z0Z_0, so that the ion-correlation change of Rosseland mean opacities at solar-interior and laboratory conditions can be evaluated consistently — bearing on whether ion correlations can materially contribute to resolving the solar opacity discrepancy.

Background

Opacity and the solar opacity problem. The radiative opacity κ(ω)\kappa(\omega) of a plasma is its photon absorption coefficient per unit mass as a function of photon energy ω\omega, comprising bound–bound (spectral line), bound–free (photoionization), and free–free absorption. Stellar-interior codes use the Rosseland mean opacity ϰR\varkappa_R — a harmonic mean of 1/κ(ω)1/\kappa(\omega) weighted by the energy derivative of the Planck distribution — which controls radiative energy transport through the Sun. The downward revision of solar photospheric abundances broke the long-standing agreement between solar models and helioseismic observations; restoring it appears to require a significant opacity increase near the base of the solar convection zone. Independently, the laboratory iron-opacity measurement at Sandia's Z machine (Bailey et al., Nature 517, 56 (2015)) found Fe opacity 30–400% above leading opacity models at T=182T=182 eV, ρ=0.17\rho=0.17 g/cm³, while chromium and nickel measured with the same technique agree with models (Nagayama et al., Phys. Rev. Lett. 122, 235001 (2019)). Identifying missing opacity microphysics is therefore an actively contested problem, and ion–ion correlations are one of the proposed candidate mechanisms whose true size is uncertain.

Average-atom models and ion correlations. An average-atom (AA) model treats a single ion embedded in a self-consistent electron distribution representing the surrounding plasma, solving Kohn–Sham density-functional-theory equations in an atomic cell. In the ion-sphere (IS) picture the ion sits at the center of a neutral Wigner–Seitz sphere of radius r0=(3/4πni)1/3r_0=(3/4\pi n_i)^{1/3} (with nin_i the ion number density), and the surrounding ions are taken as uniformly distributed, so the ion–ion pair-correlation function is the step function g(r)=θ(r−r0)g(r)=\theta(r-r_0). Ion-correlation (IC) models replace this by a realistic g(r)g(r) obtained from the Ornstein–Zernike integral equation with the hypernetted-chain (HNC) closure, using an effective interionic potential screened by the free electrons; the resulting spatial ordering of ions shifts bound-electron energies and wavefunctions, alters shell populations, and hence changes opacities. The classic Rozsnyai model (Phys. Rev. A 43, 3035 (1991)) approximates the screened interionic potential as V~II(r)=rV2(r)\tilde V_{II}(r)=rV^2(r) in terms of the self-consistent electron potential V(r)V(r). The Starrett–Saumon model (Phys. Rev. E 87, 013104 (2013); High Energy Density Physics 10, 35 (2014)) instead derives the screening from the pseudoatom free-electron density and couples the ion structure to the average-atom electronic structure consistently; it is the basis of the RESEOS code's IC opacity calculations, which use the super-transition-array (STA) method — a grouping of the astronomically many electronic configurations of mid- and high-ZZ ions into superconfigurations that makes detailed opacity computation tractable.

The chemical-potential problem. In an AA model the free-electron chemical potential μe\mu_e fixes the asymptotic free-electron density ne0n_e^0 and thereby the ionization balance — the mean ion charge Z0Z_0, the bound-shell occupation numbers, and hence the entire opacity. Without ion correlations, μe\mu_e follows from charge conservation inside the atomic cell. With ion correlations, the number of electrons per ion acquires an additional term ∫(g(r)−1)d3r\int(g(r)-1)d^3r, and the variational condition for μe\mu_e derived by Chihara (High Energy Density Physics 19, 38 (2016)) from a first-order cluster expansion of the electron number and free energy, when applied within the Starrett–Saumon model, is inconsistent with the Debye–Hückel limit — the exact weak-coupling, weak-degeneracy benchmark in which the interionic potential and pair-correlation function are known analytically. Ovechkin, Loboda, and Falkov (High Energy Density Physics 30, 29 (2019)) demonstrated this inconsistency in their Appendix A: in the limit Γei=βZ0/r0≪1\Gamma_{ei}=\beta Z_0/r_0\ll1, Γii=βZ02/r0≪1\Gamma_{ii}=\beta Z_0^2/r_0\ll1, eβμe≪1e^{\beta\mu_e}\ll1 (with Z0≫1Z_0\gg1 assumed), the model's interionic potentials and pair-correlation functions already reproduce the Debye–Hückel result, but the two sides of Chihara's variational relation reduce to −1/(β(Z0+1)ni)-1/(\beta(Z_0+1)n_i) versus μe/((Z0+1)ni)\mu_e/((Z_0+1)n_i), which cannot be reconciled; the authors attribute the defect to the use of the first-order cluster expansion for the electron number and, lacking a consistent alternative, kept the ion-sphere value of μe\mu_e, which makes the ion-correlation effect on Z0Z_0 — and therefore on opacity predictions — ambiguous. Their calculations nevertheless showed that ion correlations raise the iron Rosseland mean by only 5.5% at Z-machine conditions, far below the ~20% obtained by Krief et al. (Astrophys. J. 856, 135 (2018)) with the Rozsnyai model in the STAR code (whose own abstract reports ~15% for iron at Z-facility conditions and ~10% for the solar mixture near the convection zone); they attributed the difference to the Rozsnyai model's underestimation of ion–ion repulsion. Resolving which magnitude is right — and whether ion correlations can matter at all for the solar opacity problem — requires the chemically consistent μe\mu_e prescription that the framework still lacks.

Problem Statement

The problem, as posed by the source literature, is to remove the obstacle that currently prevents a definitive quantitative assessment of ion-correlation effects on plasma opacities within the Starrett–Saumon average-atom framework. The source's Conclusion requires a single open item — a more accurate method to calculate the electron chemical potential — while its Section 3 acknowledges several limitations of the current implementation that bear on the same quantification. The required deliverable is the item-(1) prescription — a self-consistent μe\mu_e with its Debye–Hückel limiting behavior and the removal of the Z0Z_0 ambiguity — together with the item-(3) determination, carried out consistently with that prescription, of how ion correlations change the Rosseland mean opacity at solar-interior and Z-experiment conditions; the item-(2) limitations are supplied as physical context and are not independent requirements, though a solution must state clearly where its determination depends on them.

(1) Core — chemical potential: construct a physically consistent, self-consistent method to determine the asymptotic free-electron chemical potential μe\mu_e (equivalently the asymptotic free-electron density ne0n_e^0) in the Starrett–Saumon ion-correlation model — for instance an improved cluster-expansion or variational treatment that accounts for the dependence of the pair-correlation function g(r)g(r) on ne0n_e^0 — which (i) reduces correctly to the Debye–Hückel result in the limit of weak non-ideality and weak electron degeneracy (Γei=βZ0/r0≪1\Gamma_{ei}=\beta Z_0/r_0\ll1, Γii=βZ02/r0≪1\Gamma_{ii}=\beta Z_0^2/r_0\ll1, eβμe≪1e^{\beta\mu_e}\ll1, with Z0≫1Z_0\gg1 assumed in the source analysis), resolving the inconsistency demonstrated in Appendix A of Ovechkin et al. (2019), where the model's structural quantities already reproduce the Debye–Hückel limit but the two sides of Chihara's variational condition reduce to −1/(β(Z0+1)ni)-1/(\beta(Z_0+1)n_i) versus μe/((Z0+1)ni)\mu_e/((Z_0+1)n_i), and (ii) removes the resulting ambiguity in the mean ion charge Z0Z_0 and hence in opacity predictions.

(2) Acknowledged limitations of the current implementation (described in Section 3 of the source, not posed there as open requirements), relevant as context for the quantification: (a) the ion-correlation calculations are done in the first-order two-stage approximation (ion-sphere solve, single evaluation of g(r)g(r), ion-correlation solve), though the source reports that this approach alone yields results very similar to those obtained with fully self-consistent iteration of g(r)g(r) and the electronic structure; (b) electron-collisional broadening uses frequency-independent collisional widths, which overestimate far-wing photoabsorption — without collisional broadening, the ion-correlation increase of the iron Rosseland mean at the Z-machine benchmark drops from 5.5% to 2.7%; (c) the structural properties of plasmas calculated with the Starrett–Saumon model agree well with ab initio calculations and experimental data, which the source regards as an appropriate verification of the model.

(3) Unifying objective: with a consistent μe\mu_e prescription in hand, determine defensibly how the ion-correlation change of Rosseland mean opacities at solar-interior and Z-experiment conditions is affected — in particular whether the enhancement is at the few-percent level found by Ovechkin et al. (5.5% for iron at T=182T=182 eV, ρ=0.17\rho=0.17 g/cm³) or at the larger level reported by Krief et al. with the Rozsnyai model (~10% for the solar mixture near the convection-zone boundary and ~15% for iron at Z-facility conditions per their own abstract, with the STAR opacity values quoted by Ovechkin et al. giving a 20% increase at that same iron point) — and thereby assess, in the conditional terms of the source's Conclusion, whether the account of ion correlations is important in the context of the solar opacity problem.

A resolution of the core problem would supply the corrected μe\mu_e prescription with its Debye–Hückel limiting behavior and would clarify the consequent ion-correlation-induced changes in Z0Z_0, g(r)g(r), V(r)V(r), and ϰR\varkappa_R. The source demonstrates the model at benchmark conditions including iron at T=182T=182 eV, ρ=0.17\rho=0.17 g/cm³ and T=180T=180 eV, ρ=0.473\rho=0.473 g/cm³; aluminum at T=100T=100 eV, ρ=0.0045\rho=0.0045 g/cm³; iron at T=125T=125 eV, ρ=4.46\rho=4.46 g/cm³; and bromine at T=270T=270 eV, ρ=0.026\rho=0.026 g/cm³, which illustrate the conditions at which the correction bears on opacity predictions.

References

  1. A. A. Ovechkin, P. A. Loboda, A. L. Falkov, 'Plasma opacity calculations using the Starrett and Saumon average-atom model with ion correlations', High Energy Density Physics 30, 29–40 (2019). DOI: 10.1016/j.hedp.2019.01.003, https://doi.org/10.1016/j.hedp.2019.01.003
  2. C. E. Starrett, D. Saumon, 'Electronic and ionic structures of warm and hot dense matter', Physical Review E 87, 013104 (2013). DOI: 10.1103/PhysRevE.87.013104, https://doi.org/10.1103/PhysRevE.87.013104
  3. C. E. Starrett, D. Saumon, 'A simple method for determining the ionic structure of warm dense matter', High Energy Density Physics 10, 35–42 (2014). DOI: 10.1016/j.hedp.2013.12.001, https://doi.org/10.1016/j.hedp.2013.12.001
  4. B. F. Rozsnyai, 'Photoabsorption in hot plasmas based on the ion-sphere and ion-correlation models', Physical Review A 43, 3035–3042 (1991). DOI: 10.1103/PhysRevA.43.3035, https://doi.org/10.1103/PhysRevA.43.3035
  5. M. Krief, Y. Kurzweil, A. Feigel, D. Gazit, 'The Effect of Ionic Correlations on Radiative Properties in the Solar Interior and Terrestrial Experiments', The Astrophysical Journal 856, 135 (2018). DOI: 10.3847/1538-4357/aab353, https://doi.org/10.3847/1538-4357/aab353
  6. J. E. Bailey, T. Nagayama, G. P. Loisel, G. A. Rochau, C. Blancard, J. Colgan, Ph. Cosse, G. Faussurier, C. J. Fontes, F. Gilleron, I. Golovkin, S. B. Hansen, C. A. Iglesias, D. P. Kilcrease, J. J. MacFarlane, R. C. Mancini, S. N. Nahar, C. Orban, J.-C. Pain, A. K. Pradhan, M. Sherrill, B. G. Wilson, 'A higher-than-predicted measurement of iron opacity at solar interior temperatures', Nature 517, 56–59 (2015). DOI: 10.1038/nature14048, https://doi.org/10.1038/nature14048
  7. T. Nagayama, J. E. Bailey, G. P. Loisel, G. S. Dunham, G. A. Rochau, C. Blancard, J. Colgan, Ph. Cossé, G. Faussurier, C. J. Fontes, F. Gilleron, S. B. Hansen, C. A. Iglesias, I. E. Golovkin, D. P. Kilcrease, J. J. MacFarlane, R. C. Mancini, R. M. More, C. Orban, J.-C. Pain, M. E. Sherrill, B. G. Wilson, 'Systematic Study of L-Shell Opacity at Stellar Interior Temperatures', Physical Review Letters 122, 235001 (2019). DOI: 10.1103/PhysRevLett.122.235001, https://doi.org/10.1103/PhysRevLett.122.235001
  8. C. A. Iglesias, V. Sonnad, B. G. Wilson, J. I. Castor, 'Frequency dependent electron collisional widths for opacity calculations', High Energy Density Physics 5, 97–104 (2009). DOI: 10.1016/j.hedp.2009.03.003, https://doi.org/10.1016/j.hedp.2009.03.003
  9. J. Chihara, 'Average atom model based on Quantum Hyper-Netted Chain method', High Energy Density Physics 19, 38–47 (2016). DOI: 10.1016/j.hedp.2016.03.002, https://doi.org/10.1016/j.hedp.2016.03.002
  10. N. M. Gill, C. E. Starrett, 'Mean-force scattering potential for calculating optical properties of dense plasmas', High Energy Density Physics 31, 24–30 (2019). DOI: 10.1016/j.hedp.2019.02.001, https://doi.org/10.1016/j.hedp.2019.02.001
  11. N. M. Gill, C. J. Fontes, C. E. Starrett, 'Time-dependent density functional theory applied to average atom opacity', Physical Review E 103, 043206 (2021). DOI: 10.1103/PhysRevE.103.043206, https://doi.org/10.1103/PhysRevE.103.043206
  12. N. R. Shaffer, C. E. Starrett, 'Dense plasma opacity via the multiple-scattering method', Physical Review E 105, 015203 (2022). DOI: 10.1103/PhysRevE.105.015203, https://doi.org/10.1103/PhysRevE.105.015203 (arXiv:2201.02132, https://arxiv.org/abs/2201.02132)
  13. V. V. Karasiev, S. X. Hu, N. R. Shaffer, G. Miloshevsky, 'First-principles study of L-shell iron and chromium opacity at stellar interior temperatures', Physical Review E 106, 065202 (2022). DOI: 10.1103/PhysRevE.106.065202, https://doi.org/10.1103/PhysRevE.106.065202
  14. C. E. Starrett, N. R. Shaffer, D. Saumon, R. Perriot, T. Nelson, L. A. Collins, C. Ticknor, 'Model for the electrical conductivity in dense plasma mixtures', High Energy Density Physics 36, 100752 (2020). DOI: 10.1016/j.hedp.2020.100752, https://doi.org/10.1016/j.hedp.2020.100752
  15. A. A. Ovechkin, P. A. Loboda, V. V. Popova, E. Yu. Akulinina, M. E. Berezovskaya, A. S. Korolev, S. V. Kolchugin, 'Plasma ionization balance in chemical-picture and average-atom models', Physical Review E 108, 015207 (2023). DOI: 10.1103/PhysRevE.108.015207, https://doi.org/10.1103/PhysRevE.108.015207
  16. A. A. Ovechkin, N. A. Smirnov, P. A. Loboda, A. L. Falkov, 'Calculations of the electron–ion coupling factor using the Starrett and Saumon average-atom model', Physics of Plasmas 32, 122707 (2025). DOI: 10.1063/5.0301053, https://doi.org/10.1063/5.0301053
  17. G. Faussurier, 'Neutrality in the average-atom model', Physics of Plasmas 31, 042703 (2024). DOI: 10.1063/5.0202120, https://doi.org/10.1063/5.0202120
  18. J.-C. Pain, 'Solar Opacity Calculations: Recent Theoretical Advances Prompted by Laser and Z-Pinch Experiments', Solar Physics 299, 140 (2024). DOI: 10.1007/s11207-024-02374-z, https://doi.org/10.1007/s11207-024-02374-z
  19. J. L. Zeng, C. Gao, P. F. Liu, Y. J. Li, C. S. Meng, Y. Hou, D. D. Kang, J. M. Yuan, 'Electron localization enhanced photon absorption for the missing opacity in solar interior', Science China Physics, Mechanics & Astronomy 65, 233011 (2022). DOI: 10.1007/s11433-021-1812-1, https://doi.org/10.1007/s11433-021-1812-1
  20. C. A. Iglesias, 'On the transient spatial localization model', High Energy Density Physics 47, 101043 (2023). DOI: 10.1016/j.hedp.2023.101043, https://doi.org/10.1016/j.hedp.2023.101043
  21. F. Zhou, S. Niu, S. Lu, C. Li, X. Gao, Y. Wu, Y. Qu, J. Wang, 'Plasma screening and configuration interaction effects induced large enhancement on L-shell photoionization cross sections and opacity', arXiv:2607.21238 (2026), https://arxiv.org/abs/2607.21238