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, 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.26g/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.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
The electronic structure of Pb1–xSnxTe is studied by using the relativistic Korringa-Kohn-Rostoker Green function method in the framework of density functional theory. For all concentrations x, ...Pb1–xSnxTe is a direct semiconductor with a narrow band gap. In contrast to pure lead telluride, tin telluride shows an inverted band characteristic close to the Fermi energy. It will be shown that this particular property can be tuned, first, by alloying PbTe and SnTe and, second, by applying hydrostatic pressure or uniaxial strain. Furthermore, the magnitude of strain needed to switch between the regular and inverted band gap can be tuned by the alloy composition. In conclusion, there is a range of potential usage of Pb1–xSnxTe for spintronic applications.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
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.
We present and discuss density functional theory calculations of magnetic properties of the family of ferromagnetic compounds, ( )B, focusing specifically on the magnetocrystalline anisotropy energy ...(MAE). Using periodic supercells of various sizes (up to 96 atoms), it is shown that the general qualitative features of the composition dependence of the MAE is in agreement with experimental findings, while our predicted magnitudes are larger than those of experiment. We find that the use of small supercells (6 and 12-atom) favors larger MAE values relative to a statistical sample of configurations constructed with 96-atom supercells. The effect of lattice relaxations is shown to be small. Calculations of the Curie temperature for this alloy are also presented.
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.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
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.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
In this work, we report on a methodology for the treatment of the Coulomb energy and potential in Kohn–Sham density functional theory that is free from self-interaction effects. Specifically, we ...determine the Coulomb potential given as the functional derivative of the Coulomb energy with respect to the density, where the Coulomb energy is calculated explicitly in terms of the pair density of the Kohn–Sham orbitals. This is accomplished by taking advantage of an orthonormal and complete basis that is an explicit functional of the density that then allows for the functional differentiation of the pair density with respect to the density to be performed explicitly. This approach leads to a new formalism that provides an analytic, closed-form determination of the exchange potential. This method is applied to one-dimensional model systems and to the atoms Helium through Krypton based on an exchange only implementation. Comparison of our total energies (denoted SIF) to those obtained using the usual Hartree–Fock (HF) and optimized effective potential (OEP) methods reveals the hierarchy EHF≤EOEP≤ESIF that is indicative of the greater variation freedom implicit in the former two methods.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
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.