White dwarfs represent the final state of evolution for most stars
. Certain classes of white dwarfs pulsate
, leading to observable brightness variations, and analysis of these variations with ...theoretical stellar models probes their internal structure. Modelling of these pulsating stars provides stringent tests of white dwarf models and a detailed picture of the outcome of the late stages of stellar evolution
. However, the high-energy-density states that exist in white dwarfs are extremely difficult to reach and to measure in the laboratory, so theoretical predictions are largely untested at these conditions. Here we report measurements of the relationship between pressure and density along the principal shock Hugoniot (equations describing the state of the sample material before and after the passage of the shock derived from conservation laws) of hydrocarbon to within five per cent. The observed maximum compressibility is consistent with theoretical models that include detailed electronic structure. This is relevant for the equation of state of matter at pressures ranging from 100 million to 450 million atmospheres, where the understanding of white dwarf physics is sensitive to the equation of state and where models differ considerably. The measurements test these equation-of-state relations that are used in the modelling of white dwarfs and inertial confinement fusion experiments
, and we predict an increase in compressibility due to ionization of the inner-core orbitals of carbon. We also find that a detailed treatment of the electronic structure and the electron degeneracy pressure is required to capture the measured shape of the pressure-density evolution for hydrocarbon before peak compression. Our results illuminate the equation of state of the white dwarf envelope (the region surrounding the stellar core that contains partially ionized and partially degenerate non-ideal plasmas), which is a weak link in the constitutive physics informing the structure and evolution of white dwarf stars
.
Full text
Available for:
FZAB, GEOZS, IJS, IMTLJ, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
The existence of liquid carbon as an intermediate phase preceding the formation of novel carbon materials has been a point of contention for several decades. Experimental observation of such a liquid ...state requires nonthermal melting of solid carbon materials at various laser fluences and pulse properties. Reflectivity experiments performed in the mid-1980s reached opposing conclusions regarding the metallic or insulating properties of the purported liquid state. Time-resolved X-ray absorption studies showed shortening of C–C bonds and increasing diffraction densities, thought to evidence a liquid or glassy carbon state, respectively. Nevertheless, none of these experiments provided information on the electronic structure of the proposed liquid state. Herein, we report the results of time-resolved resonant inelastic X-ray scattering (RIXS) and time-resolved X-ray emission spectroscopy (XES) studies on amorphous carbon (a-C) and ultrananocrystalline diamond (UNCD) as a function of delay time between the irradiating pulse and X-ray probe. For both a-C and UNCD, we attribute decreases in RIXS or XES signals to transition blocking, relaxation, and finally, ablation. Increased signal at 20 ps following the irradiation of the UNCD is attributed to the probable formation of nanoscale structures in the ablation plume. Differences in the amount of signal observed between a-C and UNCD are explained by the difference in sample thickness and, specifically, incomplete melting of the UNCD film. Comparisons to spectral simulations based on MD trajectories at extreme conditions indicate that the carbon state in our experiments is crystalline. Normal mode analysis confirmed that symmetrical bending or stretching of the C–C bonds in the diamond lattice results in XES spectra with small intensity differences. Overall, we observed no evidence of melting to a liquid state, as determined by the lack of changes in the spectral properties for up to 100 ps delays following the melting pulses.
Full text
Available for:
IJS, KILJ, NUK, PNG, UL, UM
We compute electrical and thermal conductivities of hydrogen plasmas in the nondegenerate regime using Kohn-Sham density functional theory (DFT) and an application of the Kubo-Greenwood response ...formula, and demonstrate that for thermal conductivity, the mean-field treatment of the electron-electron (e-e) interaction therein is insufficient to reproduce the weak-coupling limit obtained by plasma kinetic theories. An explicit e-e scattering correction to the DFT is posited by appealing to Matthiessen's Rule and the results of our computations of conductivities with the quantum Lenard-Balescu (QLB) equation. Further motivation of our correction is provided by an argument arising from the Zubarev quantum kinetic theory approach. Significant emphasis is placed on our efforts to produce properly converged results for plasma transport using Kohn-Sham DFT, so that an accurate assessment of the importance and efficacy of our e-e scattering corrections to the thermal conductivity can be made.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
We present density functional theory calculations of phosphorus dopants in bulk silicon and of several properties relating to their use as spin qubits for quantum computation. Rather than a mixed ...pseudopotential or a Heitler-London approach, we have used an explicit treatment for the phosphorus donor and examined the detailed electronic structure of the system as a function of the isotropic doping fraction, including lattice relaxation due to the presence of the impurity. Doping electron densities (rho sub(doped) - rho sub(bulk)) and spin densities (rho sub(arrowup) - rho sub(arrowdown)) are examined in order to study the properties of the dopant electron as a function of the isotropic doping fraction. Doping potentials (V sub(doped) - V sub(bulk)) are also calculated for use in calculations of the scattering cross sections of the phosphorus dopants, which are important in the understanding of electrically detected magnetic resonance experiments. We find that the electron density around the dopant leads to nonspherical features in the doping potentials, such as trigonal lobes in the (001) plane at energy scales of +12 eV near the nucleus and of -700 meV extending away from the dopants. These features are generally neglected in effective mass theory and will affect the coupling between the donor electron and the phosphorus nucleus. Our density functional calculations reveal detail in the densities and potentials of the dopants which are not evident in calculations that do not include explicit treatment of the phosphorus donor atom and relaxation of the crystal lattice. These details can also be used to parametrize tight-binding models for simulation of large-scale devices.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
A new method for the determination of clay swelling thermodynamics from computer simulation is discussed and evaluated. This method allows for the determination of temperature, pressure, and water ...chemical potential dependence of clay swelling from simulations at a single thermodynamic state point. The temperature dependence and pressure dependence of clay swelling are shown to be directly related to the composite system entropy and volume change, respectively, that accompany swelling. Expressions for the chemical potential dependence of clay swelling are used to determine constant pressure layer spacing and adsorption isotherms, quantities that are well suited for comparison with experimental measurements. This method is evaluated through grand isoshear ensemble simulations of Na-montmorillonite, a prototypical swelling clay. Approximations associated with all expressions are discussed with explicit calculations used to demonstrate their regimes of validity.
Full text
Available for:
IJS, KILJ, NUK, PNG, UL, UM
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