We present the application of Magnus based methods to the solution of first order coupled ordinary differential equations in High Energy Density (HED) physics applications. Our focus is on the application to quantum mechanical methods, specifically on the solution of the radial Dirac equation for real and complex energies. HED applications require accurate solutions across a wide range of spatial and energy domains, including regimes where the solutions exhibit pronounced oscillatory behavior. Such cases pose significant computational challenges. We demonstrate that Magnus-based integrators can efficiently and accurately address these challenges. We discuss the implementation of the Magnus method for the solution of the radial Dirac equation, including practical considerations such as the evaluation of matrix exponentials, numerical integration, error estimation, and adaptive step size control. We also discuss the application of these methods to complex energy Green’s function techniques and the efficient approximation of integrals of the solutions relevant to HED electronic structure calculations. We demonstrate the accuracy and robustness of the resulting method in applications to the free-particle case, for which analytic solutions are available for comparison, as well as the challenging case of gold at HED conditions.

We examine the performance of pure boron, boron carbide, high density carbon, and boron nitride ablators in the polar direct drive exploding pusher (PDXP) platform. The platform uses the polar direct drive configuration at the National Ignition Facility to drive high ion temperatures in a room temperature capsule and has potential applications for plasma physics studies and as a neutron source. The higher tensile strength of these materials compared to plastic enables a thinner ablator to support higher gas pressures, which could help optimize its performance for plasma physics experiments, while ablators containing boron enable the possibility of collecting additional data to constrain models of the platform. Applying recently developed and experimentally validated equation of state models for the boron materials, we examine the performance of these materials as ablators in 2D simulations, with particular focus on changes to the ablator and gas areal density, as well as the predicted symmetry of the inherently 2D implosion.

Boron carbide (B_{4}C) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. We report the results of a comprehensive computational study of the equation of state (EOS) of B_{4}C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B_{4}C over a broad range of temperatures (∼6×10^{3}-5×10^{8} K) and densities (0.025-50 g/cm^{3}). We assess that the largest discrepancies between theoretical predictions are ≲5% near the compression maximum at 1-2×10^{6} K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ∼18% and occur at temperatures between 6×10^{3}-2×10^{5} K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

The equation of state (EOS) of materials at warm dense conditions poses significant challenges to both theory and experiment. We report a combined computational, modeling, and experimental investigation leveraging new theoretical and experimental capabilities to investigate warm-dense boron nitride (BN). The simulation methodologies include path integral Monte Carlo (PIMC), several density functional theory (DFT) molecular dynamics methods [plane-wave pseudopotential, Fermi operator expansion (FOE), and spectral quadrature (SQ)], activity expansion (actex), and all-electron Green's function Korringa-Kohn-Rostoker (mecca), and compute the pressure and internal energy of BN over a broad range of densities and temperatures. Our experiments were conducted at the Omega laser facility and the Hugoniot response of BN to unprecedented pressures (1200–2650 GPa). The EOSs computed using different methods cross validate one another in the warm-dense matter regime, and the experimental Hugoniot data are in good agreement with our theoretical predictions. By comparing the EOS results from different methods, we assess that the largest discrepancies between theoretical predictions are ≲4% in pressure and ≲3% in energy and occur at 106K, slightly below the peak compression that corresponds to the K-shell ionization regime. At these conditions, we find remarkable consistency between the EOS from DFT calculations performed on different platforms and using different exchange-correlation functionals and those from PIMC using free-particle nodes. This provides strong evidence for the accuracy of both PIMC and DFT in the high-pressure, hightemperature regime. Moreover, the recently developed SQ and FOE methods produce EOS data that have significantly smaller statistical error bars than PIMC, and so represent significant advances for efficient computation at high temperatures. The shock Hugoniot predicted by PIMC, actex, and mecca shows a maximum compression ratio of 4.55±0.05 for an initial density of 2.26g/cm3, higher than the Thomas-Fermi predictions by about 5%. In addition, we construct tabular EOS models that are consistent with the firstprinciples simulations and the experimental data. Our findings clarify the ionic and electronic structure of BN over a broad range of temperatures and densities and quantify their roles in the EOS and properties of this material. The tabular models may be utilized for future simulations of laser-driven experiments that include BN as a candidate ablator material. Disciplines Condensed Matter Physics | Materials Science and Engineering | Plasma and Beam Physics Authors Shuai Zhang, Amy Lazicki, Burkhard Militzer, Lin H. Yang, Kyle Caspersen, Jim A. Gaffney, Markus W. Däne, John E. Pask, Walter R. Johnson, Abhiraj Sharma, Phanish Suryanarayana, Duane D. Johnson, Andrey V. Smirnov, Philip A. Sterne, David Erskine, Richard A. London, Federica Coppari, Damian Swift, Joseph Nilsen, Art J. Nelson, and Heather D. Whitley PHYSICAL REVIEW B 99, 165103 (2019) Equation of state of boron nitride combining computation, modeling, and experiment Shuai Zhang,1,* Amy Lazicki,1,† Burkhard Militzer,2,3,‡ Lin H. Yang,1 Kyle Caspersen,1 Jim A. Gaffney,1 Markus W. Däne,1 John E. Pask,1 Walter R. Johnson,4 Abhiraj Sharma,5 Phanish Suryanarayana,5 Duane D. Johnson,6,7 Andrey V. Smirnov,6 Philip A. Sterne,1 David Erskine,1 Richard A. London,1 Federica Coppari,1 Damian Swift,1 Joseph Nilsen,1 Art J. Nelson,1 and Heather D. Whitley1,§ 1Lawrence Livermore National Laboratory, Livermore, California 94550, USA 2Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA 3Department of Astronomy, University of California, Berkeley, California 94720, USA 4Department of Physics, 225 Nieuwland Science Hall, University of Notre Dame, Notre Dame, Indiana 46556, USA 5College of Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA 6Division of Materials Science & Engineering, Ames Laboratory, Ames, Iowa 50011, USA 7Department of Materials Science & Engineering, Iowa State University, Ames, Iowa 50011, USA (Received 2 February 2019; published 3 April 2019) The equation of state (EOS) of materials at warm dense conditions poses significant challenges to both theory and experiment. We report a combined computational, modeling, and experimental investigation leveraging new theoretical and experimental capabilities to investigate warm-dense boron nitride (BN). The simulation methodologies include path integral Monte Carlo (PIMC), several density functional theory (DFT) molecular dynamics methods [plane-wave pseudopotential, Fermi operator expansion (FOE), and spectral quadrature (SQ)], activity expansion (ACTEX), and all-electron Green’s function Korringa-Kohn-Rostoker (MECCA), and compute the pressure and internal energy of BN over a broad range of densities and temperatures. Our experiments were conducted at the Omega laser facility and the Hugoniot response of BN to unprecedented pressures (1200–2650 GPa). The EOSs computed using different methods cross validate one another in the warm-dense matter regime, and the experimental Hugoniot data are in good agreement with our theoretical predictions. By comparing the EOS results from different methods, we assess that the largest discrepancies between theoretical predictions are 4% in pressure and 3% in energy and occur at 106 K, slightly below the peak compression that corresponds to the K-shell ionization regime. At these conditions, we find remarkable consistency between the EOS from DFT calculations performed on different platforms and using different exchange-correlation functionals and those from PIMC using free-particle nodes. This provides strong evidence for the accuracy of both PIMC and DFT in the high-pressure, high-temperature regime. Moreover, the recently developed SQ and FOE methods produce EOS data that have significantly smaller statistical error bars than PIMC, and so represent significant advances for efficient computation at high temperatures. The shock Hugoniot predicted by PIMC, ACTEX, and MECCA shows a maximum compression ratio of 4.55 ± 0.05 for an initial density of 2.26 g/cm3, higher than the Thomas-Fermi predictions by about 5%. In addition, we construct tabular EOS models that are consistent with the first-principles simulations and the experimental data. Our findings clarify the ionic and electronic structure of BN over a broad range of temperatures and densities and quantify their roles in the EOS and properties of this material. The tabular models may be utilized for future simulations of laser-driven experiments that include BN as a candidate ablator material. DOI: 10.1103/PhysRevB.99.165103

The equation of state (EOS) of materials at warm dense conditions poses significant challenges to both theory and experiment. We report a combined computational, modeling, and experimental investigation leveraging new theoretical and experimental capabilities to investigate warm-dense boron nitride (BN). The simulation methodologies include path integral Monte Carlo (PIMC), several density functional theory (DFT) molecular dynamics methods [plane-wave pseudopotential, Fermi operator expansion (FOE), and spectral quadrature (SQ)], activity expansion (ACTEX), and all-electron Green's function Korringa-Kohn-Rostoker (MECCA), and compute the pressure and internal energy of BN over a broad range of densities and temperatures. Our experiments were conducted at the Omega laser facility and the Hugoniot response of BN to unprecedented pressures (1200-2650 GPa). The EOSs computed using different methods cross validate one another in the warm-dense matter regime, and the experimental Hugoniot data are in good agreement with our theoretical predictions. By comparing the EOS results from different methods, we assess that the largest discrepancies between theoretical predictions are less than or similar to 4% in pressure and less than or similar to 3% in energy and occur at 10(6) K, slightly below the peak compression that corresponds to the K-shell ionization regime. At these conditions, we find remarkable consistency between the EOS from DFT calculations performed on different platforms and using different exchange-correlation functionals and those from PIMC using free-particle nodes. This provides strong evidence for the accuracy of both PIMC and DFT in the high-pressure, high-temperature regime. Moreover, the recently developed SQ and FOE methods produce EOS data that have significantly smaller statistical error bars than PIMC, and so represent significant advances for efficient computation at high temperatures. The shock Hugoniot predicted by PIMC, ACTEX, and MECCA shows a maximum compression ratio of 4.55 +/- 0.05 for an initial density of 2.26 g/cm(3), higher than the Thomas-Fermi predictions by about 5%. In addition, we construct tabular EOS models that are consistent with the first-principles simulations and the experimental data. Our findings clarify the ionic and electronic structure of BN over a broad range of temperatures and densities and quantify their roles in the EOS and properties of this material. The tabular models may be utilized for future simulations of laser-driven experiments that include BN as a candidate ablator material.

Based on a computational procedure for determining the functional derivative with respect to the density of any antisymmetric N-particle wave function for a non-interacting system that leads to the density, we devise a test as to whether or not a wave function known to lead to a given density corresponds to a solution of a Schrödinger equation for some potential. We examine explicitly the case of non-interacting systems described by Slater determinants. Numerical examples for the cases of a one-dimensional square-well potential with infinite walls and the harmonic oscillator potential illustrate the formalism.

Magneto-electric coupling offers a new pathway to information storage in magnetic memory devices. This phenomenon has been observed in various materials ranging from insulators to semiconductors. In bulk metallic systems, magneto-electric coupling has been disregarded as the electric field cannot enter bulk metals. In this work, we show that a substantial magneto-electric coupling exists in metallic Fe nano-islands grown on Cu(111). Using the electric field in the tunnel junction of a scanning tunneling microscope, the magnetic order parameter and the crystal structure of Fe was changed on the nanometer scale. This allows high density nonvolatile information storage by means of magneto-electric coupling in a simple metallic system.

An ab initio study of magnetic exchange interactions in antiferromagnetic and strongly correlated 3d transition metal monoxides is presented. Their electronic structure is calculated using the local self-interaction correction approach, implemented within the Korringa-Kohn-Rostoker band structure method, which is based on multiple scattering theory. The Heisenberg exchange constants are evaluated with the magnetic force theorem. Based on these the corresponding Neel temperatures T_N and spin wave dispersions are calculated. The Neel temperatures are obtained using mean field approximation, random phase approximation and Monte Carlo simulations. The pressure dependence of T_N is investigated using exchange constants calculated for different lattice constants. All the calculated results are compared to experimental data.