Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Wave potential theory”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Identifying and Tuning the In Situ Oxygen-Rich Surface of Molybdenum Nitride Electrocatalysts for Oxygen Reduction

Rigorous in situ studies of electrocatalysts are required to enable the design of higher performing materials. Nonplatinum group metals for oxygen reduction reaction (ORR) catalysis containing light elements such as O, N, and C are known to be susceptible to both ex situ and in situ oxidation, leading to challenges associated with ex situ characterization methods. We have previously shown that the bulk O content plays an important role in the activity and selectivity of Mo–N catalysts, but further understanding of the role of composition and morphological changes at the surface is needed. Here, we report the measurement of in situ surface changes to a molybdenum nitride (MoN) thin film under ORR conditions using grazing incidence X-ray absorption and reflectivity. We show that the half-wave potential of MoN can be improved by ~90 mV by potential conditioning up to 0.8 V versus RHE. Utilizing electrochemical analysis, dissolution monitoring, and surface-sensitive X-ray techniques, we show that under moderate polarization (0.3–0.7 V vs RHE) there is local ligand distortion, O incorporation, and amorphization of the MoN surface, without changes in roughness. Furthermore, with a controlled potential hold procedure, we show that the surface changes concurrent with potential conditioning are stable under ORR relevant potentials. Conversely, at higher potentials (≥0.8 V vs RHE), the film incorporates O, dissolves, and roughens, suggesting that in this higher potential regime, the performance enhancements are due to increased access to active sites. Density functional theory calculations and Pourbaix analysis provide insights into film stability and O incorporation as a function of potential. These findings coupled with in situ electrochemical surface-sensitive X-ray techniques demonstrate an approach to studying nontraditional surfaces in which we can leverage our understanding of surface dynamics to improve performance with the rational, in situ tuning of active sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Superconductivity at carrier density 10 17 cm – 3 in quasi-one-dimensional Li 0.9 Mo 6 O 17

Superconductivity is a fascinating phenomenon that involves an attractive interaction by which electrons are paired and can move without electrical resistance below a critical temperature T c . A paradigm for understanding superconductive pairing, the successful “BCS” theory, developed more than 65 years ago, is challenged by certain classes of materials for which parameters (e.g. the value of T c , the density of electrons, etc.) fall outside the range for which BCS theory is applicable. The study of such materials is important because new physics may be revealed and because superconductivity has many potential applications in sensing and electronics. Superconductivity at very low carrier density is a rare and fascinating phenomenon that is of considerable importance and interest in condensed matter physics. The conditions of high density of states and effective screening of Coulomb repulsion, both key ingredients of stable Cooper pairing in the BCS theory, may not be met at low carrier density. Thus, such systems may also harbor unconventional electron pairing (e.g. non-phononic) or pairing without phase coherence as a precursor to Bose-Einstein condensation. An intriguing route toward extreme low-density (and possibly unconventional) superconductivity may be found in quasi-one-dimensional metals for which density-wave instabilities compete for the ground state and can lead to gapping of substantial portions of the Fermi surface (FS), leaving a residual (reconstructed) FS comprising a very small carrier density. Our work highlights a novel material in this unusual category: Li 0.9 Mo 6 O 17 (T c =2 K) (also known as "lithium purple bronze" or "LiPB"). We present both electrical and thermoelectric magnetotransport measurements that reveal superconductivity to occur at carrier density 2 x10 17 cm -3 , among the lowest known for any superconductor, and four orders of magnitude lower than is expected from its chemical valence and band structure. With its exceptionally large ratio of superconducting to Fermi temperature, T c /TF ≳ 0.1, LiPB is properly placed at the top of a very small group of superconductors (including cuprates, heavy-fermion and Fe-based superconductors) for which this ratio falls in the BCS-BEC crossover regime and for which pairing is believed to be spin mediated.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Gravitational wave imprints of left-right symmetric model with minimal Higgs sector

We study the gravitational wave imprints of left-right symmetric model equipped with universal seesaw mechanism allowing for the natural generation of hierarchical masses of the Standard Model fermions. The scalar sector of this model is the minimal one, consisting of only two Higgs doublets. Following the construction of the full thermal potential for this model, we perform a scan of the entire parameter space and identify the region in which the cosmic phase transition associated with the left-right symmetry breaking gives gravitational wave signals detectable by a variety of planned space-based interferometers. Then we also discuss the relevant collider implications of this beyond the Standard Model scenario.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Beyond $\rho ^{2/3}$ scaling: microscopic origins and multimessengers of high-density nuclear symmetry energy

The nature and Equation of State (EOS) of dense neutron-rich matter are still very poorly known, while they have broad impacts on many interesting issues in both astrophysics and nuclear physics. In particular, the nuclear symmetry energy $E_{\textrm{sym}}(\rho )$ encoding the cost to make nuclear matter more neutron-rich has been the most uncertain component of the EOS of dense neutron-rich nucleonic matter. It significantly affects the radii, tidal deformations, cooling rates, and frequencies of various oscillation modes of isolated neutron stars as well as the strain amplitude and frequencies of gravitational waves from their mergers, besides its many effects on the structures of nuclei as well as the dynamics and observables of their collisions. Siemens (1970s) observed that $E_{\textrm{sym}}(\rho )$ scales as $(\rho /\rho _0)^{2/3}$ near the saturation density $\rho _0$ of nuclear matter, since both the kinetic part and the potential contribution (quadratic in momentum) exhibit this dependence. The scaling holds if: (1) the nucleon isoscalar potential is quadratic in momentum, and (2) the isovector interaction is weakly density-dependent. After examining many empirical evidences and understanding theoretical findings in the literature, we conclude that: (1) Siemens’ $\rho ^{2/3}$ scaling is robust and serves as a valuable benchmark for both nuclear theories and experiments up to $2\rho _0$ but breaks down at higher densities, (2) Experimental and theoretical findings about $E_{\textrm{sym}}(\rho )$ up to $2\rho _0$ are broadly consistent, but uncertainties remain large for its curvature $K_{\textrm{sym}}(\rho )$ and higher-order parameters, (3) Above $2\rho _0$ , uncertainties grow due to poorly constrained spin-isospin-dependent tensor and three-body forces as well as the resulting nucleon short-range correlations. Looking forward, combining signals from both observations of neutron stars and terrestrial heavy-ion reaction experiments is the most promising path to finally constraining the high-density $E_{\textrm{sym}}(\rho )$ and the EOS of supradense neutron-rich matter. Multiple examples of community efforts to further constrain the high-density $E_{\textrm{sym}}(\rho )$ using both real and mocked data of present and future high-precision observations of neutron stars, as well as heavy-ion collisions involving high-energy rare isotopes, are briefly reviewed.

Li, Bao-An [East Texas A & M Univ., Commerce, TX (↗
DOE OSTI2020

Analysis of multi-electron, multi-step homogeneous catalysis by rotating disc electrode voltammetry: theory, application, and obstacles

Rotating disc electrode (RDE) voltammetry has been widely adopted for the study of heterogenized molecular electrocatalysts for multi-step fuel-forming reactions but this tool has never been comprehensively applied to their homogeneous analogues. Here, the utility and limitations of RDE techniques for mechanistic and kinetic analysis of homogeneous molecular catalysts that mediate multi-electron, multi-substrate redox transformations are explored. Using the ECEC' reaction mechanism as a case study, two theoretical models are derived based on the Nernst diffusion layer model and the Hale transformation. Current–potential curves generated by these computational strategies are compared under a variety of limiting conditions to identify conditions under which the more minimalist Nernst Diffusion Layer approach can be applied. Based on this theoretical treatment, strategies for extracting kinetic information from the plateau current and the foot of the catalytic wave are derived. RDEV is applied to a cobaloxime hydrogen evolution reaction (HER) catalyst under non-aqueous conditions in order to experimentally validate this theoretical framework and explore the feasibility of RDE as a tool for studying homogeneous catalysts. Importantly, analysis of the foot-of-the-wave via this theoretical framework provides rate constants for elementary reaction steps that agree with those extracted from stationary voltammetric methods, supporting the application of RDE to study homogeneous fuel-forming catalysts. Lastly, obstacles encountered during the kinetic analysis of cobaloxime, along with the voltammetric signatures used to diagnose this reactivity, are discussed with the goal of guiding groups working to improve RDE set-ups and help researchers avoid misinterpretation of RDE data.

08 HYDROGEN↗
DOE OSTI2021

Considering nonlocality in the optical potentials within eikonal models

Background: For its simplicity, the eikonal method is the tool of choice to analyze nuclear reactions at high energies (E > 100 MeV/nucleon), including knockout reactions. However, so far, the effective interactions used in this method are assumed to be fully local. Purpose: Given the recent studies on nonlocal optical potentials, in this work we assess whether nonlocality in the optical potentials is expected to impact reactions at high energies and then explore different avenues for extending the eikonal method to include nonlocal interactions. Method: We compare angular distributions obtained for nonlocal interactions (using the exact R-matrix approach for elastic scattering and the adiabatic distorted wave approximation for transfer) with those obtained using their local-equivalent interactions. Results: Our results show that transfer observables are significantly impacted by nonlocality in the high-energy regime. Because knockout reactions are dominated by stripping (transfer to inelastic channels), nonlocality is expected to have a large effect on knockout observables too. Three approaches are explored for extending the eikonal method to nonlocal interactions, including an iterative method and a perturbation theory. Conclusions: None of the derived extensions of the eikonal model provide a good description of elastic scattering. Here, this paper suggests that nonlocality removes the formal simplicity associated with the eikonal model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Stability of the modulator in a plasma-modulated plasma accelerator

We explore the regime of operation of the modulator stage of a recently proposed laser-plasma accelerator scheme [Phys. Rev. Lett. 127, 184801 (2021)], dubbed the plasma-modulated plasma accelerator (P-MoPA). The P-MoPA scheme offers a potential route to high-repetition-rate, GeV-scale plasma accelerators driven by picosecond-duration laser pulses from, for example, kilohertz thin-disk lasers. The first stage of the P-MoPA scheme is a plasma modulator in which a long, high-energy “drive” pulse is spectrally modulated by copropagating in a plasma channel with the low-amplitude plasma wave driven by a short, low-energy “seed” pulse. The spectrally modulated drive pulse is converted to a train of short pulses, by introducing dispersion, which can resonantly drive a large wakefield in a subsequent accelerator stage with the same on-axis plasma density as the modulator. In this paper we derive the 3D analytic theory for the evolution of the drive pulse in the plasma modulator and show that the spectral modulation is independent of transverse coordinate, which is ideal for compression into a pulse train. We then identify a transverse mode instability (TMI), similar to the TMI observed in optical fiber lasers, which sets limits on the energy of the drive pulse for a given set of laser-plasma parameters. We compare this analytic theory with particle-in-cell (PIC) simulations and find that even higher energy drive pulses can be modulated than those demonstrated in the original proposal.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Ab initio leading order effective potential for elastic proton scattering based on the symmetry-adapted no-core shell model

Calculating microscopic optical potentials for elastic scattering at intermediate energies from light nuclei in an ab initio fashion within the Watson expansion has been established within the last few years. Based on the Watson expansion of the multiple scattering series, we employ a nonlocal translationally invariant nuclear density derived within the symmetry-adapted no-core shell model (SA-NCSM) framework from a chiral next-to-next-to-leading order (NNLO) nucleon-nucleon interaction and the very same interaction for a consistent full-folding calculation of the effective (optical) potential for nucleon-nucleus scattering for medium-heavy nuclei. The leading order effective (optical) folding potential is computed by integrating over a translationally invariant SA-NCSM one-body scalar density, spin-projected momentum distribution, and the Wolfenstein amplitudes 𝐴, 𝐶, and 𝑀. The resulting nonlocal potentials serve as input for a momentum space Lippmann-Schwinger equation. In the SA-NCSM, the model space is systematically up-selected using Sp⁡(3,ℝ) symmetry considerations. For the light nucleus of 6 He, we establish a systematic selection scheme in the SA-NCSM for scattering observables. Then, we apply this scheme to calculations of scattering observables, such as differential cross sections, analyzing powers, and spin rotation functions for elastic proton scattering from 20 Ne and 40 Ca in the energy regime between 65 and 200 MeV, and compare to available data. Furthermore, our calculations show that the leading order effective nucleon-nucleus potential in the Watson expansion of multiple scattering theory obtained from an up-selected SA-NCSM model space describes 40 Ca elastic scattering observables reasonably well to about 60 degrees in the center-of-mass frame, which coincides roughly with the validity of the NNLO chiral interaction used to calculate both the nucleon-nucleon amplitudes and the one-body scalar and spin nuclear densities.

Ab initio calculations↗
DOE OSTI2022

Twistronics of Janus transition metal dichalcogenide bilayers

Twisted multilayers of two-dimensional (2D) materials are an increasingly important platform for investigating quantum phases of matter, and in particular, strongly correlated electrons. The moiré pattern introduced by the relative twist between layers creates effective potentials of long wavelength, leading to electron localization. However, in contrast to the abundance of 2D materials, few twisted heterostructures have been studied until now. Here, in this study, we develop a first-principles continuum theory to study the electronic bands introduced by moire patterns of twisted Janus transition metal dichalcogenide (TMD) homo- and heterobilayers. The model includes lattice relaxation, stacking-dependent effective mass, and Rashba spin-orbit coupling. We then perform a high-throughput generation and characterization of DFT-extracted continuum models for more than a hundred possible combinations of materials and stackings. Our model predicts that the moiré physics and emergent symmetries depend on chemical composition, vertical layer orientation, and twist angle, so that the miniband wave functions can form triangular, honeycomb, and kagome networks. Rashba spin-orbit effects, peculiar of these systems, can dominate the moiré bandwidth at small angles. Our work enables detailed investigation of Janus twisted heterostructures, allowing the discovery and control of novel electronic phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The complex heavy-quark potential in an anisotropic quark-gluon plasma — Statics and dynamics

We generalize a complex heavy-quark potential model from an isotropic QCD plasma to an anisotropic one by replacing the Debye mass m> D with an anisotropic screening mass depending on the quark pair alignment with respect to the direction of anisotropy. Such an angle-dependent mass is determined by matching the perturbative contributions in the potential model to the exact result obtained in the Hard-Thermal-Loop resummed perturbation theory. An advantage of the resulting potential model is that its angular dependence can be effectively described by using a set of angle-averaged screening masses as proposed in our previous work. Consequently, one could solve a one-dimensional Schrödinger equation with a potential model built by changing the anisotropic screening masses into the corresponding angle-averaged ones, and reproduce the full three-dimensional results for the binding energies and decay widths of low-lying quarkonium bound states to very high accuracy. Finally, turning to dynamics, we demonstrate that the one-dimensional effective potential can accurately describe the time evolution of the vacuum overlaps obtained using the full three-dimensional anisotropic potential. This includes the splitting of different p-wave polarizations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Datasets for Custom-trained Machine-learning Interatomic Potentials: Nitric Acid Aqueous Solution

This dataset was generated using an iterative active learning strategy with the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials (MLIPs) for aqueous nitric acid. Each active-learning cycle consisted of three stages: (1) training, (2) exploration, and (3) labeling. The initial training set comprised approximately 800 randomly selected configurations from a previous study by Lewis et al. (https://doi.org/10.1021/jp205510q), which investigated nitric acid solutions at 2, 3, 4, and 5 mol/L. For all configurations, single-point calculations of atomic forces and total energies were performed at the quantum density functional theory BLYP-D2 and PBE-D3 levels of theory using the CP2K Quickstep module. Valence electrons were treated explicitly, while core electrons on all atoms were represented by norm-conserving Goedecker–Teter–Hutter (GTH) pseudopotentials. Long-range dispersion interactions were accounted for using Grimme dispersion corrections. Wave functions were expanded in a mixed Gaussian-and-plane-wave scheme using TZV2P-MOLOPT basis sets for all elements and an 800 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent field convergence was accelerated using orbital transformation and Direct Inversion in the Iterative Subspace, with a convergence threshold of 10^{-6}. All single-point calculations were carried out in periodic orthorhombic cells whose dimensions match those of the molecular configurations sampled from earlier trajectories. The CELL_REF keyword in CP2K was used to define a fixed reference cell, ensuring consistency in the reference data used for MLIP training, particularly when cell fluctuations are present in NpT simulations. The resulting high-fidelity energies and forces constitute the ground-truth labels used to train the MLIPs contained in this dataset.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2025

Prototypes of Nonrelativistic Spin Splitting and Polarization in Symmetry Broken Antiferromagnets

Antiferromagnets that break both space-time reversal and translation-spin-rotation symmetries were recently predicted [L.-D. Yuan, Z. Wang, J.-W. Luo, E. I. Rashba, and A. Zunger, Phys. Rev. B 102, 014422 (2020)] to possess splitting between the otherwise spin-degenerate energy bands even without the relativistic spin-orbit coupling (SOC). Here, we point out that such nonrelativistic spin splitting (NRSS)—in particular, “spin splitting type 4” (SST-4) symmetry-broken antiferromagnets—can be divided into subgroups having distinct patterns of spin splitting and spin textures, depending on additional auxiliary symmetries of spin interconversion and polarity. These SST-4 subgroups include the 𝛼-type (no spin-interconverting symmetry) having spin splitting at the Brillouin zone center, as well as the 𝛽 subgroup in which a rotation symmetry is applied and determines the alternating spin texture and the 𝛾 subgroup having exclusively reflection spin-interconverting symmetry. Unlike ferrimagnets, the 𝛼-type compounds are shown to have tiny net magnetization at finite temperature and thus avoid the adverse effect of the stray field. The 𝛼 and 𝛽 subgroups can be either polar or nonpolar, whereas the 𝛾 subgroup is polar only, providing a basis for possible switching by external fields. The combination of NRSS-enabling and auxiliary symmetries is used here as a filter for identifying previously synthesized compounds as specific prototypes. Their characteristic splitting and spin polarization are calculated by density functional theory to the benefit of potential future experimental testing. Interesting results are as follows: (i) SOC-independent NRSS can exceed the magnitude of the SOC-induced Rashba and Dresselhaus spin splitting in semiconductors. (ii) Examples of predicted 𝛼-type insulating compounds include BiCrO 3 (nonpolar) and Mn 2 ⁢ScSbO 6 (polar), the latter having spin splitting of 158 meV and 160 meV in the valence and conduction bands, respectively. (iii) The 𝛽-type (Cu 2 ⁢Y 2 ⁢O 5 and FeF 2 ) and 𝛾-type compounds (Mn 4 ⁢Nb 2 ⁢O 9 and FeScO 3 ) are distinguished both by their auxiliary symmetries and polarity. The spin textures of 𝛾-type compounds are mirror reflected with spin degeneracy of the wave vectors on that mirror. These observations will likely broaden the experimental playing field of NRSS physics significantly.

antiferromagnets↗
DOE OSTI2024

Analytical model and dynamical phase-field simulation of terahertz transmission across ferroelectrics

Here, we theoretically investigate the steady-state transmission of continuous terahertz (THz) wave across a freestanding ferroelectric slab. Based on the Landau-Ginzburg-Devonshire theory of ferroelectrics and the coupled equations of motion for polarization and electromagnetic (EM) waves, we derive the analytical expressions of the frequency- and thickness-dependent dielectric susceptibility and transmission coefficient at the thin slab limit in the harmonic excitation regime. When the slab thickness is much smaller than the THz wavelength in the ferroelectric, the analytical predictions agree well with the numerical simulations from a dynamical phase-field model that incorporates the coupled dynamics of strain, polarization, and EM wave in multiphase systems. At larger thicknesses, the transmission is mainly determined by the frequency-dependent attenuation of THz waves in the ferroelectric and the formation of a standing polarization/THz wave. Our results advance the understanding of the interaction between THz wave and ferroelectrics and suggest the potential of exploiting ferroelectrics to achieve low-heat-dissipation, nonvolatile voltage modulation of THz transmission for high-data-rate wireless communication.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Three-dimensional higher-order saddle-point-induced flatbands in Co-based kagome metals

The saddle point (Van Hove singularity) exhibits a divergent density of states in two-dimensional systems, leading to fascinating phenomena such as strong correlations and unconventional superconductivity, yet it is seldom observed in three-dimensional (3D) systems. In this work we find two types of 3D higher-order saddle points (HOSPs) in emerging 3D kagome metals YbCo 6⁢ Ge 6 and MgCo 6 ⁢Ge 6 . Both HOSPs exhibit a singularity in their density of states, which is significantly enhanced compared to the ordinary saddle point. The HOSP near the Fermi energy generates a flatband extending a large area in the Brillouin zone, potentially amplifying the correlation effect and fostering electronic instabilities. Two types of HOSPs exhibit distinct robustness upon element substitution and lattice distortions in these kagome compounds. Our work paves the way for engineering exotic band structures, such as saddle points and flatbands, and exploring interesting phenomena in Co-based kagome materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The inflated Chern-Simons number in spectator chromo-natural inflation

The chromo-natural inflation (CNI) scenario predicts a potentially detectable chiral gravitational wave signal, generated by a Chern-Simons coupling between a rolling scalar axion field and an SU(2) gauge field with an isotropy-preserving classical background during inflation. However, the generation of this signal requires a very large integer Chern-Simons level, which can be challenging to explain or embed in a UV-complete model. We show that this challenge persists in the phenomenologically viable spectator field CNI (S-CNI) model. Furthermore, we show that a clockwork scenario giving rise to a large integer as a product of small integers can never produce a Chern-Simons level large enough to have successful S-CNI phenomenology. We briefly discuss other constraints on the model, both in effective field theory based on partial-wave unitarity bounds and in quantum gravity based on the Weak Gravity Conjecture, which may be relevant for further explorations of alternative UV completions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Band structure and Fermi surface nesting in LaSb 2

Here we use high-resolution angle-resolved photoemission spectroscopy (ARPES) and density functional theory (DFT) to investigate the electronic structure of the charge density wave (CDW) system LaSb 2 . This compound is among an interesting group of materials that manifests both a CDW transition and lower temperature superconductivity. We find the DFT calculations to be in good agreement with our ARPES data. The Fermi surface of LaSb 2 consists of two small hole pockets close to Γ and four larger pockets near the Brillouin zone boundary. The key features of the Fermi surface do not vary significantly with temperature. A saddle point is present at -0.19 eV below the Fermi level at Γ. Saddle points in the band structure have more pronounced effects on materials properties when they are located closer to the Fermi level, making doped LaSb 2 compounds a potential interesting subject of future research. Multiple peaks are present in the generalized, electronic susceptibility calculations, indicating the presence of three possible nesting vectors. We were not able to detect any signatures of the CDW transition at 355 K down to the lowest temperature of 7 K, pointing to the subtle nature of this transition. This is unusual, given that such a high transition temperature is expected to be associated with the presence of a large CDW gap. This is confirmed through investigation of the Fermi surface and through analysis of momentum distribution curves. It is possible that changes are subtle and occur below current sensitivity of our measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Experimental Measurements of Ion Diffusion Coefficients and Heating in a Multi-Ion-Species Plasma Shock

Collisional plasma shocks generated from supersonic flows are an important feature in many astrophysical and laboratory high-energy-density plasmas. Compared to single-ion-species plasma shocks, plasma shock fronts with multiple ion species contain additional structure, including interspecies ion separation driven by gradients in species concentration, temperature, pressure, and electric potential. We present time-resolved density and temperature measurements of two ion species in collisional plasma shocks produced by head-on merging of supersonic plasma jets, allowing determination of the ion diffusion coefficients. Here our results provide the first experimental validation of the fundamental inter-ion-species transport theory. The temperature separation, a higher-order effect reported here, is valuable for advancements in modeling HED and ICF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

NbV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-vanadium (Nb-V). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb and V. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name case-N where N ranges between 11 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal-amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. and Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15-50 (1996) (4) Kresse, G. and Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗