Multiple scattering theory for dense plasmas

Predicting excitation energies in warm and hot dense matter

In a dense plasma environment, the energy levels of an ion shift relative to the isolated ion values. This shift is reflected in the optical spectrum of the plasma and can be measured in, for example, emission experiments. In this work we use a recently developed method of modeling electronic states in warm dense matter to predict these level energies. In this model excited state energies are calculated directly by enforcing constrained one-electron occupation factors, thus allowing the calculation of specific transition and ionization energies. This model includes plasma effects self-consistently, so the effect of continuum lowering is included in an ab initio sense. We use the model to calculate the K-edge and K-alpha energies of solid density magnesium, aluminum, and silicon over a range of temperatures, finding close agreement with experimental results. We also calculate the ionization potential depression to compare to widely used models and investigate the effects of temperature on the lowering of the continuum.

Excited states in warm and hot dense matter

Accurate modeling of warm and hot dense matter is challenging in part due to the multitude of excited states that must be considered. Here, we present a variational framework that models these excited states. In this framework an excited state is defined by a set of effective one-electron occupation factors, and the corresponding energy is defined by the effective one-body energy with an exchange and correlation term. The variational framework is applied to an atom-in-plasma model (a generalization of the so-called average atom model). Comparisons with a density functional theory based average atom model generally reveal good agreement in the calculated pressure, but our model also gives access to the excitation energies and charge state distributions.

Terahertz Spatiotemporal Wave Synthesis in Random Systems

Complex media have emerged as a powerful and robust framework to control light-matter interactions designed for task-specific optical functionalities. Studies on wavefront shaping through disordered systems have demonstrated optical wave manipulation capabilities beyond conventional optics, including aberration-free and subwavelength focusing. However, achieving arbitrary and simultaneous control over the spatial and temporal features of light remains challenging. In particular, no practical solution exists for field-level arbitrary spatiotemporal control of wave packets. A new paradigm shift has emerged in the terahertz frequency domain, offering methods for absolute time-domain measurements of the scattered electric field, enabling direct field-based wave synthesis. In this work, we report the experimental demonstration of field-level control of single-cycle terahertz pulses on arbitrary spatial points through complex disordered media.

Spectral-partitioned Kohn-Sham density functional theory

We introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

Dense plasma opacity via the multiple-scattering method

The calculation of the optical properties of hot dense plasmas with a model that has self-consistent plasma physics is a grand challenge for high energy density science. Here we exploit a recently developed electronic structure model that uses multiple scattering theory to solve the Kohn-Sham density functional theory equations for dense plasmas. We calculate opacities in this regime, validate the method, and apply it to recent experimental measurements of opacity for Cr, Ni, and Fe. Good agreement is found in the quasicontinuum region for Cr and Ni, while the self-consistent plasma physics of the approach cannot explain the observed difference between models and the experiment for Fe.

Effect of ionic disorder on the principal shock Hugoniot

The effect of ionic disorder on the principal Hugoniot is investigated using multiple scattering theory to very high pressure (Gbar). Calculations using molecular dynamics to simulate ionic disorder are compared to those with a fixed crystal lattice, for both carbon and aluminum. For the range of conditions considered here we find that ionic disorder has a relatively minor influence. It is most important at the onset of shell ionization and we find that, at higher pressures, the subtle effect of the ionic environment is overwhelmed by the larger number of ionized electrons with higher thermal energies.

Equation of state modeling with pseudoatom molecular dynamics

Using a modified version of the pseudoatom molecular-dynamics approach, the silicon and oxygen equations of state were generated and then employed to construct the equation of state of silicon dioxide. The results are supported by the close agreement with ab initio simulations of the silicon pressure and experimental shock Hugoniot of silicon dioxide. Ion thermal contributions to thermodynamic functions provided by the PAMD simulations are compared to their counterparts obtained with the one-component plasma and charged-hard-sphere approximations.