Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “charge density waves”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Dual-cone variational calculation of the two-electron reduced density matrix

The computation of strongly correlated quantum systems is challenging because of its potentially exponential scaling in the number of electron configurations. Variational calculation of the two-electron reduced density matrix (2-RDM) without the many-electron wave function exploits the pairwise nature of the electronic Coulomb interaction to compute a lower bound on the ground-state energy with polynomial computational scaling. Recently, a dual-cone formulation of the variational 2-RDM calculation was shown to generate the ground-state energy, albeit not the 2-RDM, at a substantially reduced computational cost, especially for higher $\textit{N}$-representability conditions such as the T2 constraint. Here we generalize the dual-cone variational 2-RDM method to compute not only the ground-state energy but also the 2-RDM. The central result is that we can compute the 2-RDM from a generalization of the Hellmann-Feynman theorem. Specifically, we prove that in the Lagrangian formulation of the dual-cone optimization the 2-RDM is the Lagrange multiplier. Further, we apply the method to computing the energies and properties of strongly correlated electrons—including atomic charges, electron densities, dipole moments, and orbital occupations—in an illustrative hydrogen chain and the nitrogen-fixation catalyst FeMoco. The dual variational computation of the 2-RDM with T2 or higher N -representability conditions provides a polynomially scaling approach to strongly correlated molecules and materials with significant applications in atomic and molecular and condensed-matter chemistry and physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The role of high-order electron correlation effects in a model system for non-valence correlation-bound anions

The diffusion Monte Carlo (DMC), auxiliary field quantum Monte Carlo (AFQMC), and equation-of-motion coupled cluster (EOM-CC) methods are used to calculate the electron binding energy (EBE) of the non-valence anion state of a model (H2O)4 cluster. Two geometries are considered, one at which the anion is unbound and the other at which it is bound in the Hartree–Fock (HF) approximation. It is demonstrated that DMC calculations can recover from the use of a HF trial wave function that has collapsed onto a discretized continuum solution, although larger EBEs are obtained when using a trial wave function for the anion that provides a more realistic description of the charge distribution and, hence, of the nodal surface. For the geometry at which the cluster has a non-valence correlation-bound anion, both the inclusion of triples in the EOM-CC method and the inclusion of supplemental diffuse d functions in the basis set are important. DMC calculations with suitable trial wave functions give EBE values in good agreement with our best estimate EOM-CC result. AFQMC using a trial wave function for the anion with a realistic electron density gives a value of the EBE nearly identical to the EOM-CC result when using the same basis set. For the geometry at which the anion is bound in the HF approximation, the inclusion of triple excitations in the EOM-CC calculations is much less important. The best estimate EOM-CC EBE value is in good agreement with the results of DMC calculations with appropriate trial wave functions.

Upadhyay, Shiv (ORCID:0000000285010501)↗
DOE OSTI2022

Mixed excitonic nature in water-oxidized BiVO 4 surfaces with defects

BiVO 4 is a promising photocatalyst for efficient water oxidation, with surface reactivity determined by the structure of active catalytic sites. Surface oxidation in the presence of oxygen vacancies induces electron localization, suggesting an atomistic route to improve the charge transfer efficiency within the catalytic cycle. In this paper, we study the effect of oxygen vacancies on the electronic and optical properties at BiVO 4 surfaces upon water oxidation. We use density functional theory and many-body perturbation theory to explore the change in the electronic and quasiparticle energy levels and to evaluate the electron-hole coupling as a function of the underlying structure. Further we show that while the presence of defects alters the atomic structure and largely modifies the wave-function nature, leading to defect-localized states at the quasiparticle gap region, the optical excitations remain largely unchanged due to the substantial hybridization of defect and nondefect electron-hole transitions. Our findings suggest that defect-induced surface oxidation supports improved electron transport, both through bound and tunable electronic states and via a mixed nature of the optical transitions, expected to reduce electron-hole defect trapping.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nuclear-level effective theory of 𝜇→𝑒 conversion: Formalism and applications

Over the next decade new 𝜇→𝑒 conversion searches at Fermilab (Mu2e) and J-PARC (COMET, DeeMe) are expected to advance limits on charged lepton flavor violation (CLFV) by more than four orders of magnitude. Here, by considering the consequence of 𝑃 and 𝐶⁢𝑃 on elastic 𝜇→𝑒 conversion and the structure of possible charge and current densities, we show that rates are governed by six nuclear responses and a single scale, 𝑞/𝑚 𝑁 , where 𝑞 ≈ 𝑚 𝜇 is the momentum transferred from the leptons to the nucleus. To relate this result to microscopic formulations of CLFV, we construct in nonrelativistic effective theory (NRET) the CLFV nucleon-level interaction, pointing out the relevance of the dimensionless scales 𝑦=($^{𝑞⁢𝑏}_2$) 2 ⁢ > |$\overrightarrow{𝑣}_N$⁢|>|$\overrightarrow{𝑣}_μ$⁢|>|$\overrightarrow{𝑣}_T$| , where 𝑏 is the nuclear size, $\overrightarrow{𝑣}_N$ and $\overrightarrow{𝑣}_μ$ are the nucleon and muon intrinsic velocities, and $\overrightarrow{𝑣}_T$ is the target recoil velocity. We discuss previous work, noting the lack of a systematic treatment of the various small parameters. Because the parameter 𝑦 is not small, a proper calculation of 𝜇→𝑒 conversion requires a full multipole expansion of the nuclear response functions, an apparently daunting task with Coulomb-distorted electron partial waves. We demonstrate that the multipole expansion can be carried out to high precision by introducing a simplifying local momentum 𝑞 eff for the electron. Previous work has been limited to simple charge or spin interactions, thereby treating the nucleus effectively as a point particle. We show that such formulations are not compatible with the general form of the 𝜇→𝑒 conversion rate, failing to generate three of the six allowed nuclear response functions. The inclusion of the nucleon velocity $\overrightarrow{𝑣}_N$ yields an NRET with 16 operators and a rate of the general form. Consequently, in the current discovery era for CLFV, it provides the most sensible starting point for experimental analysis, defining what can and cannot be determined about CLFV from the highly exclusive process of 𝜇→𝑒 conversion. Finally, we expand the NRET operator basis to account for the effects of $\overrightarrow{𝑣}_μ$, associated with the muon's lower component, generating corrections to the CLFV coefficients of the point-nucleus response functions. Using advanced shell-model methods, we compute 𝜇→𝑒 conversion rates for a series of experimental targets, deriving bounds on the coefficients of the CLFV operators. These calculations are the first to include a general basis of CLFV operators, full evaluation of the associated nuclear response functions, and an accurate treatment of electron and muon Coulomb effects. We discuss target selection as an experimental “knob” that can be turned to probe the microscopic origins of CLFV. We describe two types of coherence that enhance certain CLFV operators and selection rules that blind elastic 𝜇→𝑒 conversion to others. We discuss the matching of the NRET onto higher level effective field theories, such as those constructed at the light quark level, noting opportunities to build on existing work in direct detection of dark matter. We discuss the relation of 𝜇→𝑒 conversion to 𝜇→𝑒+𝛾 and 𝜇→3⁢𝑒, showing how MEG II and Mu3e results will complement those of Mu2e and COMET. Finally we describe a accompanying script—in Mathematica and Python versions—that can be used to compute 𝜇→𝑒 conversion rates in various nuclear targets for the full set of NRET operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

On anomalous transport of multi-species plasma associated with the resistive ballooning and resistive drift waves driven turbulence

Anomalous transport of multi-species plasma related to the resistive ballooning and resistive drift wave turbulence is considered in a “cold” ion approximation. It is found that similar to the resistive drift wave turbulence [see A. R. Knyazev and S. I. Krasheninnikov, Phys. Plasmas 31, 012502 (2024); and S. I. Krasheninnikov and R. D. Smirnov, Phys. Plasmas (to be published)] the addition of the ballooning drive does not change the main features of anomalous transport of the multi-species plasma: (i) The transport of all ion species is described as a transport of the passive scalars in the turbulent field of the electrostatic potential and electron density perturbation; (ii) the density of ion species with a larger ratio of the mass to charge has the tendency to the accumulation/depletion in the vortices of plasma flow; and (iii) the cross-field transport of all plasma species (including electrons and ions) is described by the same anomalous transport coefficient.

Physics↗
DOE OSTI2024

GPAW: An open Python package for electronic structure calculations

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe–Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn–Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

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

Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell

Development of particle-in-cell (PIC) methods using finite element based methods (FEMs) have been a topic of renewed interest; this has largely been driven by (a) the ability of finite element methods to better model geometry, (b) better understanding of function spaces that are necessary to represent all Maxwell quantities, and (c) more recently, the fundamental rubrics that should be obeyed in space and time so as to satisfy Gauss' laws and the equation of continuity. In that vein, methods have been developed recently that satisfy these equations and are agnostic to time stepping methods. While this development is indeed a significant advance, it should be noted that implicit FEM transient solvers support an underlying null space that corresponds to a gradient of a scalar potential ∇Φ(r) (or t∇Φ(r) in the case of wave equation solvers). While explicit schemes do not suffer from this drawback, they are only conditionally stable with time step sizes that are mesh dependent and very small. Furthermore, the null space produces spurious charge that can corrupt the desired physics of a PIC simulation. The way to overcome this bottleneck, and indeed, satisfy all four Maxwell's equation is to use a quasi-Helmholtz formulation on a tessellation. In the re-formulation presented, we strictly satisfy the equation of continuity and Gauss' laws for both the electric and magnetic flux densities. Results illustrating the efficacy of this scheme will be demonstrated by analyzing non-neutral beams with space-charge and the adiabatic expansion of a neutral plasma with realistic parameters (Debye length and real mass ratios).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Overdriven-detonation states produced by spherically diverging waves

A series of experiments are currently underway at the Detonation Science and Technology group, within Los Alamos National Laboratory, to study the overdriven-detonation states achievable in energetic materials as a result of detonation-wave interactions. A multi-component, energetic-material array was designed to study the amplification of velocity and pressure states produced by spherically diverging detonation waves in pentaerythritol tetranitrate (PETN) output charges. The unique geometry provides a low-jitter, highly controlled series of interactions between three independent-detonation inputs. Streak-camera imaging was performed on the output face of PETN pellets ranging in thickness from 2.5 – 10 mm to characterize the resulting breakout profile. Additionally, photonic Doppler velocimetry (PDV) measurements were collected at the output-pellet surface to determine simultaneity within the system. Detonation-wave velocities upwards of 16 mm/µs were measured, as compared to a steady-state detonation velocity of 7.9 mm/µs for the PETN pressing density investigated. Finally, additional experiments are being conducted to measure the pressure amplification generated at key areas of interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Ground-state phase diagram of the t-t ′ -J model

Significance Theoretical efforts to understand the superconducting cuprates have focused on a few simplified but challenging models: notably the Hubbard and the t - t ′ - J models. Using density matrix renormalization group (DMRG) simulations on width-8 cylinders, we study the ground-state phase diagram of the t - t ′- J model. The next nearest-neighbor hopping t ′ allows one to distinguish between hole ( t ′ < 0 ) and electron doping ( t ′ > 0 ). We find for the magnetic, charge, and single-particle properties that this model roughly captures the differences between the electron- and hole-doped cuprates. However, the superconducting properties differ. For hole doping ( t ′ < 0 ), there is an absence of superconductivity, while for electron doping ( t ′ > 0), we find d -wave superconductivity coexisting with antiferromagnetism and a secondary triplet p -wave pairing component.

Jiang, Shengtao↗
DOE OSTI2023

Microscopic calculation of fission product yields for odd-mass nuclei

Fission data are essential inputs to reaction networks involved in nucleosynthesis simulations and nuclear forensics. In such applications as well as in the description of multichance fission, the characteristics of fission for odd-mass nuclei are just as important as those for even-even nuclei. The fission theories that aim at explicitly describing fission dynamics are typically based on some variant of the nuclear mean-field theory. In such cases, the treatment of systems with an odd number of particles is markedly more involved, both formally and computationally. In this article, we use the blocking prescription of the Hartree-Fock-Bogoliubov theory with Skyrme energy functionals to compute the deformation properties of odd-mass uranium isotopes. We show that the resulting fission fragment distributions depend quite significantly on the spin of the odd neutron. By direct calculation of the spin distribution of the fissioning nucleus, we propose a methodology to rigorously predict the charge and mass distributions in odd-mass nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Charged particle tracking with quantum annealing optimization

Abstract At the High Luminosity Large Hadron Collider (HL-LHC), traditional track reconstruction techniques that are critical for physics analysis will need to be upgraded to scale with track density. Quantum annealing has shown promise in its ability to solve combinatorial optimization problems amidst an ongoing effort to establish evidence of a quantum speedup. As a step towards exploiting such potential speedup, we investigate a track reconstruction approach by adapting the existing geometric Denby-Peterson (Hopfield) network method to the quantum annealing framework for HL-LHC conditions. We develop additional techniques to embed the problem onto existing and near-term quantum annealing hardware. Results using simulated annealing and quantum annealing with the D-Wave 2X system on the TrackML open dataset are presented, demonstrating the successful application of a quantum annealing algorithm to the track reconstruction challenge. We find that combinatorial optimization problems can effectively reconstruct tracks, suggesting possible applications for fast hardware-specific implementations at the HL-LHC while leaving open the possibility of a quantum speedup for tracking.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Identifying Majorana vortex modes via nonlocal transport

The combination of two-dimensional Dirac surface states with s-wave superconductivity is expected to generate localized topological Majorana zero modes in vortex cores. Putative experimental signatures of these modes have been reported for heterostructures of proximitized topological insulators, iron-based superconductors or certain transition metal dichalcogenides. Despite these efforts, the Majorana nature of the observed excitation is still under debate. In this work, we propose to identify the presence of Majorana vortex modes using a nonlocal transport measurement protocol originally employed for one-dimensional settings. In the case of an isolated subgap state, the protocol provides a spatial map of the ratio of local charge- and probability-density which offers a clear distinction between Majorana and ordinary fermionic modes. We show that these distinctive features survive in the experimentally relevant case of hybridizing vortex core modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

NbZr_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-zirconium (Nb-Zr). 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 Zr. 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. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10, 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. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of 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. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system 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. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & 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. & 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. & 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. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2024

TaZr_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 tantalum-zirconium (Ta-Zr). 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 Ta and Zr. 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 80, 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. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10 and between 81 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. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of 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. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system 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. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & 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. & 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. & 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. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2020

Excluded-volume model for quarkyonic matter: Three-flavor baryon-quark mixture

The single-flavor excluded-volume model based on the effective size of baryons reproduces the hard-soft density evolution of the equation of state (EoS) required by the recent studies of GW170817. This phenomenological model basically realizes the concept of quarkyonic matter which is introduced from large-$N_c$ gauge theory for dense matter. Further, enhanced nucleon interactions and dynamically generated quark degrees of freedom can reproduce the hard-soft evolution of the EoS. In this paper, we extend the excluded-volume model to a three-flavor system by considering electromagnetic charge and possible weak equilibrium in order to obtain a proper description for the hard-soft behavior of the EoS inferred from the gravitational waves observations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Positive Energy Functional for Massless Scalars in Rotating Black Hole Backgrounds of Maximal Ungauged Supergravity

Here, we outline a proof of the stability of a massless neutral scalar field ψ in the background of a wide class of four dimensional asymptotically flat rotating and “electrically charged” solutions of supergravity, and the low energy limit of string theory, known as STU metrics. Despite their complexity, we find it possible to circumvent the difficulties presented by the existence of ergo regions and the related phenomenon of superradiance in the original metrics by following a strategy due to Whiting, and passing to an auxiliary metric admitting an everywhere lightlike Killing field and constructing a scalar field ψ (related to a possible unstable mode ψ by a nonlocal transformation) which satisfies the massless wave equation with respect to the auxiliary metric. By contrast with the case for ψ, the associated energy density of ψ is not only conserved but is also non-negative.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Ultrafast Molecular Frame Quantum Tomography

Here, we develop and experimentally demonstrate a methodology for a full molecular frame quantum tomography (MFQT) of dynamical polyatomic systems. We exemplify this approach through the complete characterization of an electronically nonadiabatic wave packet in ammonia (NH 3 ). The method exploits both energy and time-domain spectroscopic data, and yields the lab frame density matrix (LFDM) for the system, the elements of which are populations and coherences. The LFDM fully characterizes electronic and nuclear dynamics in the molecular frame, yielding the time- and orientation-angle dependent expectation values of any relevant operator. For example, the time-dependent molecular frame electronic probability density may be constructed, yielding information on electronic dynamics in the molecular frame. In NH 3 , we observe that electronic coherences are induced by nuclear dynamics which nonadiabatically drive electronic motions (charge migration) in the molecular frame. Here, the nuclear dynamics are rotational and it is nonadiabatic Coriolis coupling which drives the coherences. Interestingly, the nuclear-driven electronic coherence is preserved over longer timescales. In general, MFQT can help quantify entanglement between electronic and nuclear degrees of freedom, and provide new routes to the study of ultrafast molecular dynamics, charge migration, quantum information processing, and optimal control schemes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

The elphbolt ab initio solver for the coupled electron-phonon Boltzmann transport equations

elphbolt is a modern Fortran (2018 standard) code for efficiently solving the coupled electron–phonon Boltzmann transport equations from first principles. Using results from density functional and density functional perturbation theory as inputs, it can calculate the effect of the non-equilibrium phonons on the electronic transport (phonon drag) and non-equilibrium electrons on the phononic transport (electron drag) in a fully self-consistent manner and obeying the constraints mandated by thermodynamics. It can calculate the lattice, charge, and thermoelectric transport coefficients for the temperature gradient and electric fields, and the effect of the mutual electron–phonon drag on these transport properties. The code fully exploits the symmetries of the crystal and the transport-active window to allow the sampling of extremely fine electron and phonon wave vector meshes required for accurately capturing the drag phenomena. The corray feature of modern Fortran, which offers native and convenient support for parallelization, is utilized. The code is compact, readable, well-documented, and extensible by design.

electron drag↗