Source author record

Heather D. Whitley

Heather D. Whitley appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2022arXiv

Improved analysis of converging shock radiographs

We previously reported an experimental platform to induce a spherically-convergent shock in a sample using laser-driven ablation, probed with time-resolved x-ray radiography, and an analysis method to deduce states along the principal shock Hugoniot simultaneously with the x-ray opacity. We have now developed a modified method of analysis that is numerically better-conditioned and faster, and usually provides a better representation of the radiograph with correspondingly lower uncertainties. The previous approach was based on optimizing parameters in a model of the density distribution as a function of radius and time, warped to follow loci such as the shock and the outside of the sample. The converging shock configuration can be described more efficiently in terms of the shocked density and sound speed, expressed as functions of the shock speed Studies of the Hugoniot from various theoretical equations of state (EOS) indicate that, in the typical range of states explored by these experiments, these functions can be described by low-order polynomials. Similarly, few-parameter functions were found suitable for representing the variation of x-ray opacity with shock pressure. This approach was found to perform better in most cases than an alternative method based on parameterization of the EOS.

preprint2022arXiv

Shock Hugoniot of diamond from 3 to 80 TPa

The principal Hugoniot of carbon, initially diamond, was measured from 3 to 80 TPa (30 to 800 million atmospheres), the highest pressure ever achieved, using radiography of spherically-converging shocks. The shocks were generated by ablation of a plastic coating by soft x-rays in a laser-heated hohlraum at the National Ignition Facility (NIF). Experiments were performed with low and high drive powers, spanning different but overlapping pressure ranges. The radius-time history of the shock, and the profile of mass density behind, were determined by profile-matching from a time-resolved x-ray radiograph across the diameter of the sphere. Above ~50 TPa, the heating induced by the shock was great enough to ionize a significant fraction of K-shell electrons, reducing the opacity to the 10.2 keV probe x-rays. The opacity and mass density were deduced simultaneously using the constraint that the total mass of the sample was constant. The Hugoniot and opacity were consistent with density functional theory calculations of the electronic states and equation of state (EOS), and varied significantly from theoretical Hugoniots based on Thomas-Fermi theory. Theoretical models used to predict the compressibility of diamond ablator experiments at the NIF, producing the highest neutron yields so far from inertial confinement fusion experiments, are qualitatively consistent with our EOS measurements but appear to overpredict the compressibility slightly. These measurements help to evaluate theoretical techniques and constrain wide-range EOS models applicable to white dwarf stars, which are the ultimate evolutionary form of at least 97% of stars in the galaxy.

preprint2021arXiv

Atom-in-jellium equations of state and melt curves in the white dwarf regime

Atom-in-jellium calculations of the electron states, and perturbative calculations of the Einstein frequency, were used to construct equations of state (EOS) from around $10^{-5}$ to $10^7$g/cm$^3$ and $10^{-4}$ to $10^{6}$eV for elements relevant to white dwarf (WD) stars. This is the widest range reported for self-consistent electronic shell structure calculations. Elements of the same ratio of atomic weight to atomic number were predicted to asymptote to the same $T=0$ isotherm, suggesting that, contrary to recent studies of the crystallization of WDs, the amount of gravitational energy that could be released by separation of oxygen and carbon is small. A generalized Lindemann criterion based on the amplitude of the ion-thermal oscillations calculated using atom-in-jellium theory, previously used to extrapolate melt curves for metals, was found to reproduce previous thermodynamic studies of the melt curve of the one component plasma with a choice of vibration amplitude consistent with low pressure results. For elements for which low pressure melting satisfies the same amplitude criterion, such as Al, this melt model thus gives a likely estimate of the melt curve over the full range of normal electronic matter; for the other elements, it provides a useful constraint on the melt locus.

preprint2020arXiv

Phase transformation in boron under shock compression

Using first-principles molecular dynamics, we calculated the equation of state and shock Hugoniot of various boron phases. We find a large mismatch between Hugoniots based on existing knowledge of the equilibrium phase diagram and those measured by shock experiments, which could be reconciled if the $α$-B$_{12}$/$β\rightarrowγ$-B$_{28}$ transition is significantly over-pressurized in boron under shock compression. Our results also indicate that there exists an anomaly and negative Clapeyron slope along the melting curve of boron at 100 GPa and 1500--3000 Kelvin. These results enable in-depth understanding of matter under shock compression, in particular the significance of compression-rate dependence of phase transitions and kinetic effects in experimental measurements.

preprint2020arXiv

Simultaneous compression and opacity data from time-series radiography with a Lagrangian marker

Time-resolved radiography can be used to obtain absolute shock Hugoniot states by simultaneously measuring at least two mechanical parameters of the shock, and this technique is particularly suitable for one-dimensional converging shocks where a single experiment probes a range of pressures as the converging shock strengthens. However, at sufficiently high pressures, the shocked material becomes hot enough that the x-ray opacity falls significantly. If the system includes a Lagrangian marker, such that the mass within the marker is known, this additional information can be used to constrain the opacity as well as the Hugoniot state. In the limit that the opacity changes only on shock heating, and not significantly on subsequent isentropic compression, the opacity of shocked material can be determined uniquely. More generally, it is necessary to assume the form of the variation of opacity with isentropic compression, or to introduce multiple marker layers. Alternatively, assuming either the equation of state or the opacity, the presence of a marker layer in such experiments enables the non-assumed property to be deduced more accurately than from the radiographic density reconstruction alone. An example analysis is shown for measurements of a converging shock wave in polystyrene, at the National Ignition Facility.

preprint2013arXiv

Large-scale atomistic density functional theory calculations of phosphorus-doped silicon quantum bits

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 and spin densities are examined in order to study the properties of the dopant electron as a function of the isotropic doping fraction. Doping potentials 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 non-spherical 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 parameterize tight-binding models for simulation of large-scale devices.

preprint2003arXiv

Structure and energetics of helium adsorption on nanosurfaces

The ground and excited state properties of small helium clusters, 4He_N, containing nanoscale (~3-10 Angstroms) planar aromatic molecules have been studied with quantum Monte Carlo methods. Ground state structures and energies are obtained from importance-sampled, rigid-body diffusion Monte Carlo. Excited state energies due to helium vibrational motion are evaluated using the projection operator, imaginary time spectral evolution technique. We examine the adsorption of N helium atoms (N less than or equal to 24) on a series of planar aromatic molecules (benzene, naphthalene, anthracene, tetracene, phthalocyanine). The first layer of helium atoms is well-localized on the molecule surface, and we find well-defined localized excitations due to in-plane vibrational motion of helium on the molecule surface. We discuss the implications of these confined excitations for the molecule spectroscopy.