Perspective: New directions in dynamical density functional theory

Classical dynamical density functional theory (DDFT) has become one of the central modeling approaches in nonequilibrium soft matter physics. Recent years have seen the emergence of novel and interesting fields of application for DDFT. In particular, there has been a remarkable growth in the amount of work related to chemistry. Moreover, DDFT has stimulated research on other theories such as phase field crystal models and power functional theory. In this perspective, we summarize the latest developments in the field of DDFT and discuss a variety of possible directions for future research.

Bridge functions of classical one-component plasmas

In a recent paper, Lucco Castello et al. [arXiv:2107.03537] performed systematic extractions of classical one-component plasma bridge functions from molecular dynamics simulations and provided an accurate parametrization that was incorporated in their isomorph-based empirically modified hypernetted chain approach for Yukawa one-component plasmas. Here the extraction technique and parametrization strategy are described in detail, while the deficiencies of earlier efforts are discussed. The structural and thermodynamic predictions of the updated version of the integral equation theory approach are compared with extensive available simulation results revealing a truly unprecedented level of accuracy in the entire dense liquid region of the Yukawa phase diagram.

Transport properties of a quasisymmetric binary nitrogen-oxygen mixture in the warm dense regime

Transport properties of mixtures in the warm dense matter (WDM) regime play an important role in natural astrophysics. However, a physical understanding of ionic transport properties in quasisymmetric liquid mixtures has remained elusive. Here, we present extensive ab initio molecular dynamics (AIMD) simulations on the ionic diffusion and viscosity of a quasisymmetric binary nitrogen-oxygen (N-O) mixture in a wide warm dense regime of 8-120 kK and 4.5-8.0 g/cm^{3}. Diffusion and viscosity of N-O mixtures with different compositions are obtained by using the Green-Kubo formula. Unlike asymmetric mixtures, the change of proportions in N-O mixtures slightly affects the viscosity and diffusion in the strong-coupling region. Furthermore, the AIMD results are used to build and verify a global pseudo-ion in jellium (PIJ) model for ionic transport calculations. The PIJ model succeeds in reproducing the transport properties of N-O mixtures where ionization has occurred, and provides a promising alternative approach to obtaining comparable results to AIMD simulations with relatively small computational costs. Our current results highlight the characteristic features of the quasisymmetric binary mixtures and demonstrate the applicability of the PIJ model in the WDM regime.

The five problems of irreversibility

Thermodynamics has a clear arrow of time, characterized by the irreversible approach to equilibrium. This stands in contrast to the laws of microscopic theories, which are invariant under time-reversal. Foundational discussions of this "problem of irreversibility" often focus on historical considerations, and do therefore not take results of modern physical research on this topic into account. In this article, I will close this gap by studying the implications of dynamical density functional theory (DDFT), a central method of modern nonequilibrium statistical mechanics not previously considered in philosophy of physics, for this debate. For this purpose, the philosophical discussion of irreversibility is structured into five problems, concerned with the source of irreversibility in thermodynamics, the definition of equilibrium and entropy, the justification of coarse-graining, the approach to equilibrium and the arrow of time. For each of these problems, it is shown that DDFT provides novel insights that are of importance for both physicists and philosophers of physics.

Sculpted ultracold neutral plasmas

Ultracold neutral plasma (UNP) experiments allow for careful control of plasma properties across Coulomb coupling regimes. Here, we examine how UNPs can be used to study heterogeneous, nonequilibrium phenomena, including nonlinear waves, transport, hydrodynamics, kinetics, stopping power, and instabilities. Through a series of molecular dynamics simulations, we have explored UNPs formed with spatially modulated ionizing radiation. We have developed a computational model for such sculpted UNPs that includes an ionic screened Coulomb interaction with a spatiotemporal screening length, and Langevin-based spatial ion-electron and ion-neutral collisions. We have also developed a hydrodynamics model and have extracted its field quantities (density, flow velocity, and temperature) from the molecular dynamics simulation data, allowing us to investigate kinetics by examining moment ratios and phase-space dynamics; we find that it is possible to create UNPs that vary from nearly perfect fluids (Euler limit) to highly kinetic plasmas. We have examined plasmas in three geometries: a solid rod, a hollow rod, and a gapped slab; we have studied basic properties of these plasmas, including the spatial Coulomb coupling parameter. By varying the initial conditions, we find that we can design experimental plasmas that would allow the exploration of a wide range of phenomena, including shock and blast waves, stopping power, two-stream instabilities, and much more. Using an evaporative cooling geometry, our results suggest that much larger Coulomb couplings can be achieved, possibly in excess of 10.

High-frequency response of classical strongly coupled plasmas

The dynamic structure factor (DSF) of the Yukawa system is here obtained with highly converged molecular dynamics (MD) over the entire liquid phase. The data provide a rigorous test of theoretical models of ion-acoustic wave-dispersion relations, the intermediate scattering function, and the high-frequency response. We compare our MD results with seven diverse models, finding good agreement among those that enforce the three basic sum rules for dispersion properties, although one of the models has previously unreported spurious peaks. The MD simulations reveal that at intermediate frequencies ω, the high-frequency response of the DSF follows a power law, going approximately as ω^{-p}, where p>0, and p shows nontrivial dependencies on the wave vector q and the plasma parameters κ and Γ. In contrast, among the seven comparison models, the predicted high-frequency response is found to be independent of {q,κ,Γ}. This high-frequency power suggests a useful fitting form. In addition, these results expose limitations of several models and, moreover, suggest that some approaches are difficult or impossible to extend because of the lack of finite moments. We also find the double-plasmon resonance peak in our MD simulations that none of the theoretical models predicts.

Excess pressure and electric fields in nonideal plasma hydrodynamics

Nonideal plasmas have nontrivial space and time correlations, which simultaneously impact both the excess thermodynamic quantities as well as the collision processes. However, hydrodynamics models for designing and interpreting nonideal plasma experiments, such as inertial-confinement fusion experiments, typically neglect electrodynamics, although some models include electric fields indirectly through a generalized Fick's law. However, because most transport models are not computed self-consistently with the equation of state, there is double counting of the forces in the excess thermodynamic quantities and the collision terms. Here we employ the statistical mechanical hydrodynamic theory of Irving and Kirkwood [J. Chem. Phys. 18, 817 (1950)JCPSA60021-960610.1063/1.1747782] to examine inhomogeneous, nonideal plasmas that contain electric fields. We show that it is not possible to simultaneously separate terms that correspond to electric fields and excess pressure; rather, these quantities arise from the same interparticle Coulomb forces. Moreover, new terms associated with nonlocality appear in the presence of strong inhomogeneities.