Researcher profile

Damian C. Swift

Damian C. Swift contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 published item(s)

preprint2022arXiv

Development of slurry targets for high repetition-rate XFEL experiments

Combining an x-ray free electron laser (XFEL) with high power laser drivers enables the study of phase transitions, equation-of-state, grain growth, strength, and transformation pathways as a function of pressure to 100s GPa along different thermodynamic compression paths. Future high-repetition rate laser operation will enable data to be accumulated at >1 Hz which poses a number of experimental challenges including the need to rapidly replenish the target. Here, we present a combined shock-compression and X-ray diffraction study on vol% epoxy(50)-crystalline grains(50) (slurry) targets, which can be fashioned into extruded ribbons for high repetition-rate operation. For shock-loaded NaCl-slurry samples, we observe pressure, density and temperature states within the embedded NaCl grains consistent with observations for shock-compressed single-crystal NaCl.

preprint2022arXiv

Dislocation-based strength model for high energy density conditions

We derive a continuum-level plasticity model for polycrystalline materials in the high energy density regime, based on a single dislocation density and single mobility mechanism, with an evolution model for the dislocation density. The model is formulated explicitly in terms of quantities connected closely with equation of state (EOS) theory, in particular the shear modulus and Einstein temperature, which reduces the number of unconstrained parameters while increasing the range of applicability. The least constrained component is the Peierls barrier $E_P$, which is however accessible by atomistic simulations. We demonstrate an efficient method to estimate the variation of $E_P$ with compression, constrained to fit a single flow stress datum. The formulation for dislocation mobility accounts for some or possibly all of the stiffening at high strain rates usually attributed to phonon drag. The configurational energy of the dislocations is accounted for explicitly, giving a self-consistent calculation of the conversion of plastic work to heat. The configurational energy is predicted to contribute to the mean pressure, and may reach several percent in the terapascal range, which may be significant when inferring scalar EOS data from dynamic loading experiments. The bulk elastic strain energy also contributes to the pressure, but appears to be much smaller. Although inherently describing the plastic relaxation of elastic strain, the model can be manipulated to estimate the flow stress as a function of mass density, temperature, and strain rate, which is convenient to compare with other models and inferences from experiment. The deduced flow stress reproduces systematic trends observed in elastic waves and instability growth experiments, and makes testable predictions of trends versus material and crystal type over a wide range of pressure and strain rate.

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.

preprint2021arXiv

Atom-in-jellium predictions of the shear modulus at high pressure

Atom-in-jellium calculations of the Einstein frequency in condensed matter and of the equation of state were used to predict the variation of shear modulus from zero pressure to ~$10^7$ g/cm$^3$, for several elements relevant to white dwarf (WD) stars and other self-gravitating systems. This is by far the widest range reported electronic structure calculation of shear modulus, spanning from ambient through the one-component plasma to extreme relativistic conditions. The predictions were based on a relationship between Debye temperature and shear modulus which we assess to be accurate at the o(10%) level, and is the first known use of atom-in-jellium theory to calculate a shear modulus. We assessed the overall accuracy of the method by comparing with experimental measurements and more detailed electronic structure calculations at lower pressures.

preprint2021arXiv

Equation of state and strength of diamond in high pressure ramp loading

Diamond is used extensively as a component in high energy density experiments, but existing equation of state (EOS) models do not capture its observed response to dynamic loading. In particular, in contrast with first principles theoretical EOS models, no solid-solid phase changes have been detected, and no general-purpose EOS models match the measured ambient isotherm. We have performed density functional theory (DFT) calculations of the diamond phase to ~10TPa, well beyond its predicted range of thermodynamic stability, and used these results as the basis of a Mie-Greuneisen EOS. We also performed DFT calculations of the elastic moduli, and calibrated an algebraic elasticity model for use in simulations. We then estimated the flow stress of diamond by comparison with the stress-density relation measured experimentally in ramp-loading experiments. The resulting constitutive model allows us to place a constraint on the Taylor-Quinney factor (the fraction of plastic work converted to heat) from the observation that diamond does not melt on ramp compression.

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.

preprint2019arXiv

High pressure melt locus of iron from atom-in-jellium calculations

Although usually considered as a technique for predicting electron states in dense plasmas, atom-in-jellium calculations can be used to predict the mean displacement of the ion from its equilibrium position in colder matter, as a function of compression and temperature. The Lindemann criterion of a critical displacement for melting can then be employed to predict the melt locus, normalizing for instance to the observed melt temperature or to more direct simulations such as molecular dynamics (MD). This approach reproduces the high pressure melting behavior of Al as calculated using the Lindemann model and thermal vibrations in the solid. Applied to Fe, we find that it reproduces the limited-range melt locus of a multiphase equation of state (EOS) and the results of ab initio MD simulations, and agrees less well with a Lindemann construction using an older EOS. The resulting melt locus lies significantly above the older melt locus for pressures above 1.5\,TPa, but is closer to recent ab initio MD results and extrapolations of an analytic fit to them. This study confirms the importance of core freezing in massive exoplanets, predicting that a slightly smaller range of exoplanets than previously assessed would be likely to exhibit dynamo generation of magnetic fields by convection in the liquid portion of the core.

preprint2019arXiv

High temperature ion-thermal behavior from average-atom calculations

Atom-in-jellium calculations of the Einstein frequency were used to calculate the mean displacement of an ion over a wide range of compression and temperature. Expressed as a fraction of the Wigner-Seitz radius, the displacement is a measure of the asymptotic freedom of the ion at high temperature, and thus of the change in heat capacity from 6 to 3 quadratic degrees of freedom per atom. A functional form for free energy was proposed based on the Maxwell-Boltzmann distribution as a correction to the Debye free energy, with a single free parameter representing the effective density of potential modes to be saturated. This parameter was investigated using molecular dynamics simulations, and found to be ~0.2 per atom. In this way, the ion-thermal contribution can be calculated for a wide-range equation of state (EOS) without requiring a large number of molecular dynamics simulations. Example calculations were performed for carbon, including the sensitivity of key EOS loci to ionic freedom.

preprint2018arXiv

Non-Iterative Characteristics Analysis for High-Pressure Ramp Loading

In the canonical ramp compression experiment, a smoothly-increasing load is applied to the surface of the sample, and the particle velocity history is measured at two or more different distances into the sample, at interfaces where the surface of the sample can be probed. The velocity histories are used to deduce a stress-density relation, usually using iterative Lagrangian analysis to account for the perturbing effect of the impedance mismatch at the interface. In that technique, a stress- density relation is assumed in order to correct for the perturbation, and is adjusted until it becomes consistent with the deduced stress-density relation. This process is subject to the usual difficulties of nonlinear optimization, such as the existence of local minima (sensitivity to the initial guess), possible failure to converge, and relatively large computational effort. We show that, by considering the interaction of successive characteristics reaching the interfaces, the stress-density relation can be deduced directly by recursion rather than iteration. This calculation is orders of magnitude faster than iterative analysis, and does not require an initial guess. Direct recursion may be less suitable for very noisy data, but it was robust when applied to trial data. The stress-density relation deduced was identical to the result from iterative Lagrangian analysis.