Search NASA⌕ Search

Engineering topics

Correa, Alfredo A.

Publications and source records attributed to Correa, Alfredo A..

Software stewardship and advancement of a high-performance computing scientific application: QMCPACK

Here, we provide an overview of the software engineering efforts and their impact in QMCPACK, a production-level ab-initio Quantum Monte Carlo open-source code targeting high-performance computing (HPC) systems. Aspects included are: (i) strategic expansion of continuous integration (CI) targeting CPUs, using GitHub Actions own runners, and NVIDIA and AMD GPUs used in pre-exascale systems, (ii) incremental reduction of memory leaks using sanitizers, (iii) incorporation of Docker containers for CI and reproducibility, and (iv) refactoring efforts to improve maintainability, testing coverage, and memory lifetime management. We quantify the value of these improvements by providing metrics to illustrate the shift towards a predictive, rather than reactive, maintenance approach. Our goal, in documenting the impact of these efforts on QMCPACK, is to contribute to the body of knowledge on the importance of research software engineering (RSE) for the stewardship and advancement of community HPC codes to enable scientific discovery at scale.

97 MATHEMATICS AND COMPUTING↗

Direct observation of strong momentum-dependent electron-phonon coupling in a metal

Phonon scattering in metals is one of the most fundamental processes in materials science. However, understanding such processes has remained challenging and requires detailed information on interactions between phonons and electrons. We use an ultrafast electron diffuse scattering technique to resolve the nonequilibrium phonon dynamics in femtosecond–laser-excited tungsten in both time and momentum. We determine transient populations of phonon modes which show strong momentum dependence initiated by electron-phonon coupling. For phonons near Brillouin zone border, we observe a transient rise in their population on a timescale of approximately 1 picosecond driven by the strong electron-phonon coupling, followed by a slow decay on a timescale of approximately 8 picosecond governed by the weaker phonon-phonon relaxation process. We find that the exceptional harmonicity of tungsten is needed for isolating the two processes, resulting in long-lived nonequilibrium phonons in a pure metal. Our finding highlights that electron-phonon scattering can be the determinant factor in the phonon thermal transport of metals.

36 MATERIALS SCIENCE↗

Ultrafast lattice disordering can be accelerated by electronic collisional forces

In the prevalent picture of ultrafast structural phase transitions, the atomic motion occurs in a slowly varying potential energy surface determined adiabatically by the fast electrons. However, this ignores non-conservative forces caused by electron-lattice collisions, which can significantly influence atomic motion. Most ultrafast techniques only probe the average structure and are less sensitive to random displacements, and therefore do not detect the role played by non-conservative forces in phase transitions. Here we show that the lattice dynamics of the prototypical insulator-to-metal transition of VO 2 cannot be described by a potential energy alone. We use the sample temperature to control the preexisting lattice disorder before ultrafast photoexcitation across the phase transition and our ultrafast diffuse scattering experiments show that the fluctuations characteristic of the rutile metal develop equally fast (120 fs) at initial temperatures of 100 K and 300 K. This indicates that additional non-conservative forces are responsible for the increased lattice disorder. In conclusion, these results highlight the need for more sophisticated descriptions of ultrafast phenomena beyond the Born-Oppenheimer approximation as well as ultrafast probes of spatial fluctuations beyond the average unit cell measured by diffraction.

36 MATERIALS SCIENCE↗

Electron dynamics in extended systems within real-time time-dependent density-functional theory

Due to a beneficial balance of computational cost and accuracy, real-time time-dependent density-functional theory has emerged as a promising first-principles framework to describe electron real-time dynamics. Here we discuss recent implementations around this approach, in particular in the context of complex, extended systems. Results include an analysis of the computational cost associated with numerical propagation and when using absorbing boundary conditions. Here we extensively explore the shortcomings for describing electron–electron scattering in real time and compare to many-body perturbation theory. Modern improvements of the description of exchange and correlation are reviewed. In this work, we specifically focus on the Qb@ll code, which we have mainly used for these types of simulations over the last years, and we conclude by pointing to further progress needed going forward.

2D materials↗

Radioactive decay of $\mathrm {{}^{90}Sr}$ in cement: a non-equilibrium first-principles investigation

Cement is an inexpensive and relatively easily manageable material that is used as a last barrier for nuclear waste disposal. Under these conditions, the cement is in contact with low radiation doses, but there is a distinct possibility of being contaminated with radioactive products. Of particular concern is the medium lived half-life product 90 Sr (28.8 years) due to its ability to replace Ca. 90 Sr undergoes β-decay to 90 Y which, in turn, β-decays to stable 90 Zr. In this work, we discuss systematically the chain of non-equilibrium processes that result as a consequence of β-decay events in cement. We first use density functional-based methods to study the consequences of the sudden increase of the nuclear charge from Z to Z+1, a possible induced ionization and the perturbation of the surrounding electronic charge. Secondly, we use molecular dynamics simulations to study the recoil of the daughter nucleus. Finally, we discuss the damage caused by the ionization cascade produced during the propagation of the β-electron and the resulting chemical and structural perturbation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantum Non-Equilibrium Dynamics Prediction of Electronic Transport Coefficients in Non-Linear Regimes

The goal of this work was to develop unique LLNL capabilities for the computational prediction of electrical, optical and thermal transport properties from quantum simulations. Using a recently developed theory of microscopic quantum transport, we were able to model compressed metallic hydrogen and it’s non-linear and out-of-equilibrium contributions in electrical and thermal transport. This theory was implemented in the open-source, real-space time-dependent density functional theory code Octopus. This model is general and widely applicable to high energy density systems that are important to LLNL and NNSA mission areas.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A measurement of the equation of state of carbon envelopes of white dwarfs

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. In this paper 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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗