Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Pseudopotentials”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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 OSTI2024

VZr_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 vanadium-zirconium (V-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 V 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.

36 MATERIALS SCIENCE↗
DOE OSTI2025

TaTi_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-titanium (Ta-Ti). 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. 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 Hf and Nb. 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 Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HfNb_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 hafnium-niobium (Hf-Nb). 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. 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 Hf and Nb. 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,029 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
NASA NTRS1982

An investigation of impurity centers in semiconductors of variable composition. Part 1: General theory and some applications

A theory of deep point defects imbedded in otherwise perfect semiconductor crystals is developed with the aid of pseudopotentials. The dominant short-range forces engendered by the impurity are sufficiently weakened in all cases where the cancellation theorem of the pseudopotential formalism is operative. Thus, effective-mass-like equations exhibiting local effective potentials derived from nonlocal pseudopotentials are shown to be valid for a large class of defects. A two-band secular determinant for the energy eigenvalues of deep defects is also derived from the set of integral equations which corresponds to the set of differential equations of the effective-mass type. Subsequently, the theory in its simplest form, is applied to the system Al(x)Ga(1-x)As:Se. It is shown that the one-electron donor level of Se within the forbidden gap of Al(x)Ga(1-x)As as a function of the AlAs mole fraction x reaches its maximum of about 300 meV (as measured from the conduction band edge) at the cross-over from the direct to the indirect band-gap at x = 0.44 in agreement with experiments.

Vonroos, O. H.↗
DOE OSTI2024

Spectral scheme for atomic structure calculations in density functional theory

In this study, we present a spectral scheme for atomic structure calculations in pseudopotential Kohn-Sham density functional theory. In particular, after applying an exponential transformation of the radial coordinates, we employ global polynomial interpolation on a Chebyshev grid, with derivative operators approximated using the Chebyshev differentiation matrix, and integrations using Clenshaw-Curtis quadrature. We demonstrate the accuracy and efficiency of the scheme through spin-polarized and unpolarized calculations for representative atoms, while considering local, semilocal, and hybrid exchange-correlation functionals. In particular, we find that $\mathcal{O}$(200) grid points are sufficient to achieve an accuracy of 1 microhartree in the eigenvalues for optimized norm conserving Vanderbilt pseudopotentials spanning the periodic table from atomic number Ζ = 1 to 83.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1978

Adhesion of a bimetallic interface

The Hohenberg-Kohn and Kohn-Sham formalisms are used to examine binding (binding energy as a function of separation) for combinations of the simple metals Al(111), Zn(0001), Mg(0001), and Na(110) in contact. Similar metal contacts between Al, Zn, Mg, and Na are examined self-consistently in an ab initio calculation using the Kohn-Sham formalism. Crystallinity is included using the Aschroft pseudopotential via first order perturbation theory for the electron-ion interaction; and the ion-ion interaction is included exactly via a lattice sum. Binding energy was determined both in the local-density approximation and including gradient corrections to the exchange and correlation energy. Binding was found in all cases. In dissimilar metal contacts, interfacial bonding was greater than that in the weaker material predicting the possibility of metallic transfer. The nonzero position of the energy minimum in like metal contacts is explained in terms of consistency between the Ashcroft pseudopotential and the bulk charge density. Good agreement with experimental surface energies is obtained in the self-consistent calculation when nonlocal terms are included.

Ferrante, J.↗
NASA NTRS1985

Thermodynamic properties of liquid Si and Ge

Thermodynamic properties of liquid Si and Ge are calculated using standard variational techniques, with the hard-sphere system as a reference system. Third-order electron-ion pseudopotential contributions are included nonvariationally and are shown to have little effect on the free energy, although a larger one on the derivatives of the free energy. Two different model pseudopotentials are used and give similar results in good agreement with experiment.

Shih, W.-H.↗
DOE OSTI2025

Modelling Thomson Scattering in a Hydrogen Plasma at Stellar Interior Conditions Using the Hypernetted‐Chain Approach

Under the extreme conditions found in small stars, where electron degeneracy and Coulomb coupling are significant, accurate modeling of Thomson scattering is crucial for determining opacity, a primary quantity for stellar energy transport. We use hypernetted-chain calculations, incorporating quantum pseudopotentials and electron-exchange effects to obtain the electron–electron static structure factor to calculate the Thomson scattering transport cross-section for conditions prevailing in the interior of small stars. These results are compared to those from average-atom simulations and analytical calculations. Our findings support laboratory astrophysics experiments aimed at benchmarking opacity models for stellar interiors, particularly for red dwarf stars, and help to bridge theoretical models with observations.

ICF↗
DOE OSTI2025

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗
DOE OSTI2024

Prediction of Redox Potentials for U, Np, Pu, and Am in Aqueous Solution

The redox properties of the actinides in aqueous solution are important for fuel production/reprocessing and understanding the environmental impact of nuclear waste. The redox potentials for U, Np, Pu, and Am in oxidation states from 0 up to VII (as appropriate) in aqueous solutions have been predicted at the density functional theory level with the B3LYP functional, Stuttgart small core pseudopotential basis sets for the actinides, and explicit (30H 2 O molecules)/implicit treatment of the aqueous solvent using the self-consistent reaction field COSMO and SMD approaches for the implicit solvation. The predictions of the structural parameters of clusters incorporating first and second solvation shells are consistent with the available experimental data., Our results are typically within 0.2 V of the available experimental data using two explicit solvation shells with an implicit solvent model. The use of the PW91 functional substantially improved the prediction of the Pu(VI/V) redox couple. The redox couples for An(VI/IV) and An(V/IV) which involve the addition of protons and removal of the actinyl oxygens led to slightly larger differences from experiment. Here, the An(IV/0) and An(III/0) couples were reliably predicted with our approach. Predictions of the unknown An(II/I) redox potentials were negative, consistent with expectations, and predictions for unknown An(VII/VI), An(III/II), and An(II/0) redox couples improve prior estimates.

Actinides↗
DOE OSTI2024

Prediction of Redox Potentials for Ac, Th, and Pa in Aqueous Solution

Density functional theory in conjunction with small core pseudopotentials and the associated basis sets was used to calculate potentials for multiple redox couples, covering a range of oxidation states for Ac (0 to III), Th (0 to IV), and Pa (0 to V) in aqueous solution. Solvation effects were incorporated using a supermolecule-continuum approach, with 30 water molecules representing two solvation shells, and the COSMO and SMD implicit solvation models. The calculated geometries for Ac(III), Th(IV), and Pa(V) were in reasonable agreement with the available experimental data. Using the COSMO model with the B3LYP functional, the calculated redox potentials were within ± 0.2 V from experiment for most redox couples. Several pathways were explored for the Pa(V/IV) redox couple for different forms of Pa(V) and Pa(IV). Most Pa(V/IV) redox couples have very similar potentials, ranging from 0 to -0.4 V up to a pH of 1.4. At pH = 1.4, the potentials shift to values that are more negative than -0.7 V, reflecting the growing unfavorable nature of the redox process at higher pH levels. The calculated values for An(III/II) potentials were consistent with prior estimates and the available experimental data. The predicted redox potentials for An(II/I) were highly negative, as expected. For An(I/0) potentials, Th and Pa exhibited positive values, contrasting with the negative values calculated for Ac. Furthermore, the An +m /An(0) potentials agreed better with the experimental data when using the COSMO solvation model as compared to the SMD model.

Chemical calculations↗
DOE OSTI2025

Prediction of Redox Potentials for Different Oxidation States of U, Np, Pu, and Am in Alkaline Aqueous Solution

The redox potentials for U, Np, Pu, and Am for oxidation states +III up to +VIII in alkaline aqueous solutions were predicted using density functional theory (DFT) and small-core pseudopotentials and their basis sets, with a hybrid explicit/implicit solvent model using SHE = 4.28 V. For each oxidation state, various oxo/hydroxo complexes were evaluated resulting in a variety of one-electron redox pathways. For An(VIII/VII) couples, the predicted redox potentials for the [An(VIII)O 5 (OH)] -3 /[An(VII)O 4 (OH) 2 ] -3 , or [An(VIII)O 4 (OH) 2 ] -2 /[An(VII)O 4 (OH) 2 ] -3 couples are in good agreement with existing estimates. For An(VII/VI) redox couples, all couples, particularly [An(VII)O 4 (OH) 2 ] -3 /[An(VI)O 2 (OH) 4 ] -2 , were in agreement with experimental values for U, Np, and Pu, but the results for Am showed larger differences from the estimated potentials. The An(VI/V) couples were consistent with experiment for dioxo/tetrahydroxo couples, and the An(V/IV) couples showed acceptable agreement based on actinide-specific couples, with neutral hydroxides often favored in the +IV state. The An(IV/III) couples were consistent with literature values when modeled as soluble neutral hydroxides. The use of our approach yielded calculated redox potentials that were within ±0.2 V from experimental or estimated values consistent with our prior calculations on redox potentials of actinides from Ac to Am in acidic aqueous solutions. Finally, as a result, this supports the robustness of our DFT-based methodology for predicting actinide redox potentials, offering valuable insights for actinide chemistry in aqueous solutions.

Actinides↗
DOE OSTI2025

Prediction of Redox Potentials for the Late Actinides Cm to Lr Using Electronic Structure Methods

Our previously developed computational method for calculating the aqueous redox potentials of the early actinides has been extended to the later elements in the actinide series: Cm, Bk, Cf, Es, Fm, Md, No, and Lr in multiple oxidation states. These calculations were performed using density functional theory with small-core pseudopotentials and their associated basis sets. Solvation effects were considered via a supermolecule-continuum approach, with 30 water molecules representing two solvation shells. Both the COSMO and SMD implicit solvation models were utilized. The structural parameters and hydration numbers for Cm(III), Bk(III), Bk(IV), and Cf(III) are in reasonable agreement with the available experimental data. For redox processes involving atomic cations in solution, the B3LYP/COSMO approach predicted redox potentials to within ±0.2 V of experiment for most redox couples, consistent with our prior work. Inclusion of spin-orbit corrections in specific redox pairs, especially those with the later actinides in high oxidation states, yields improved results relative to calculations including only scalar-relativistic corrections. The An +m /An(0) redox potentials were calculated using a Born-Haber cycle incorporating sublimation, ionization, and hydration energies. Due to a lack of experimental data, three sets of ionization energies were used for the Born-Haber cycle. The calculated An(III/0) potentials showed better agreement with experimental data when using the COSMO solvation model and the test set comprising the NIST recommended ionization energies. Furthermore, the Md(II/0) potential was better described with the SMD model, whereas No(II/0) was not well described by all methods. Finally, the computational approach was able to predict redox potentials that for most cases agreed with the current available experimental or estimated data.

Actinides↗
DOE OSTI2024

Ligand-Induced Size-Dependent Circular Dichroism in Quantum Dots

Recent experiments have probed the chiral properties of semiconductor nanocrystal (NC) quantum dots (QDs), but understanding the circular dichroism line shape, excitonic features, and chirality induction mechanism remains a challenge. We propose an atomistic pseudopotential method to model chiral ligand passivated QDs, computing circular dichroism (CD) spectra for CdSe QDs (2.6-3.8 nm). We find strong agreement between calculated and measured line shapes, predicting consistent bisignate line shapes with decreasing CD magnitude as size increases. Further, our analysis reveals the origin of bisignate line shapes, arising from nondegenerate excitons with opposing angular momenta. We also explore the impact of chiral ligand orientation on QD surfaces, observing changes in the optical activity magnitude and sign. This orientation sensitivity offers the means to distinguish ordered from disordered ligand configurations, facilitating the study of order-disorder transitions at ligand-QD interfaces.

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

Low-energy pathways lead to self-healing defects in CsPbBr 3

Self-regulation of free charge carriers in perovskites via Schottky defect formation has been posited as the origin of the well-known defect-tolerance of metal halide perovskite materials. Understanding the mechanisms of self-regulation promises to lead to the fabrication of better performing solar cell materials with higher efficiencies. We investigate many such mechanisms here for CsPbBr 3 , a popular representative of a more commercially viable all-inorganic metal halide perovskite. We investigate different atomic-level mechanisms and pathways of the diffusion and recombination of neutral and charged interstitials and vacancies (Schottky pairs) in CsPbBr 3 . We use nudged elastic band calculations and ab initio-derived pseudopotentials within quantum ESPRESSO to determine energies of formation and migration and hence the activation energies for these defects. While halide vacancies are known to exhibit low formation energies, the migration of interstitials is less studied. Our calculations uncover interstitial defect pathways capable of producing an activation energy at, or below, the single experimental value of 0.53 eV observed for the slow, temperature-dependent recovery of light-induced conductivity in bulk CsPbBr 3 . Our work reveals the existence of a low-energy diffusion pathway involving a concerted “domino effect” of interstitials, with the net result that interstitials can diffuse more readily over longer distances than expected. This observation suggests that defect self-healing can be promoted if the “domino effect” strategy can be engaged.

14 SOLAR ENERGY↗
DOE OSTI2024

Tunable noninteracting free-energy density functionals for high-energy-density physics applications

In this work, we introduce the concept of a tunable noninteracting free-energy density functional and present two examples realized: (i) via a simple one-parameter convex combination of two existing functionals and (ii) via the construction of a generalized gradient approximation (GGA) enhancement factor that contains one free parameter and is designed to satisfy a set of incorporated constraints. Functional (i), constructed as a combination of the local Thomas–Fermi and a pseudopotential-adapted GGA for the noninteracting free-energy, has already demonstrated its practical usability for establishing the high temperature end of the equation of state of deuterium [Phys. Rev. B 104, 144104 (2021)] and CHON resin [Phys. Rev. E 106, 045207 (2022)] for inertial confinement fusion applications. Hugoniot calculations for liquid deuterium are given as another example of how the application of computationally efficient orbital-free density functional theory (OF-DFT) can be utilized with the employment of the developed functionals. Once the functionals have been tuned such that the OF-DFT Hugoniot calculation matches the Kohn–Sham solution at some low-temperature point, agreement with the reference Kohn–Sham results for the rest of the high temperature Hugoniot path is very good with relative errors for compression and pressure on the order of 2% or less.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗