Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Pseudopotential”

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.

107 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Radial Two-Dimensional Ion Crystals in a Linear Paul Trap

We experimentally study two-dimensional (2D) Coulomb crystals in the “radial-2D” phase of a linear Paul trap. This phase is identified by a 2D ion lattice aligned entirely with the radial plane and is created by imposing a large ratio of axial to radial trapping potentials. Using arrays of up to 19 171 Yb + ions, we demonstrate that the structural phase boundaries of such crystals are well described by the pseudopotential approximation, despite the time-dependent ion positions driven by intrinsic micromotion. We further observe that micromotion-induced heating of the radial-2D crystal is confined to the radial plane. Finally, we verify that the transverse motional modes, which are used in most ion-trap quantum simulation schemes, are well-predictable numerically and remain decoupled and cold in this geometry. Our results establish radial-2D ion crystals as a robust experimental platform for realizing a variety of theoretical proposals in quantum simulation and computation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Quantifying uncertainty in high-throughput density functional theory: A comparison of AFLOW, Materials Project, and OQMD

A central challenge in high-throughput density functional theory (HT-DFT) calculations is selecting a combination of input parameters and postprocessing techniques that can be used across all materials classes, while also managing accuracy-cost tradeoffs. To investigate the effects of these parameter choices, we consolidate three large HT-DFT databases: Automatic-FLOW (AFLOW), the Materials Project (MP), and the Open Quantum Materials Database (OQMD), and compare reported properties across each pair of databases for materials calculated using the same initial crystal structure. We find that HT-DFT formation energies and volumes are generally more reproducible than band gaps and total magnetizations; for instance, a notable fraction of records disagree on whether a material is metallic (up to 7%) or magnetic (up to 15%). The variance between calculated properties is as high as 0.105 eV/atom (median relative absolute difference, or MRAD, of 6%) for formation energy, 0.65 Å 3 / atom (MRAD of 4%) for volume, 0.21 eV (MRAD of 9%) for band gap, and 0.15 μB/ formula unit (MRAD of 8%) for total magnetization, comparable to the differences between DFT and experiment. Here, we trace some of the larger discrepancies to choices involving pseudopotentials, the DFT + U formalism, and elemental reference states, and argue that further standardization of HT-DFT would be beneficial to reproducibility.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Band structures and Z 2 invariants of two-dimensional transition metal dichalcogenide monolayers from fully relativistic Dirac-Kohn-Sham theory using Gaussian-type orbitals

Two-dimensional (2D) materials exhibit a wide range of remarkable phenomena, many of which owe their existence to the relativistic spin-orbit coupling (SOC) effects. To understand and predict properties of materials containing heavy elements, such as the transition-metal dichalcogenides (TMDs), relativistic effects must be taken into account in first-principles calculations. We present an all-electron method based on the four-component Dirac Hamiltonian and Gaussian-type orbitals (GTOs) that overcomes complications associated with linear dependencies and ill-conditioned matrices that arise when diffuse functions are included in the basis. Until now, there has been no systematic study of the convergence of GTO basis sets for periodic solids either at the nonrelativistic or the relativistic level. Here, in this study, we provide such a study of relativistic band structures of the 2D TMDs in the hexagonal (2H), tetragonal (1T), and distorted tetragonal (1T') structures, along with a discussion of their SOC-driven properties (Rashba splitting and Z 2 topological invariants). We demonstrate the viability of our approach even when large basis sets with multiple basis functions involving various valence orbitals (denoted triple- and quadruple-ζ) are used in the relativistic regime. Our method does not require the use of pseudopotentials and provides access to all electronic states within the same framework. Our study paves the way for direct studies of material properties, such as the parameters in spin Hamiltonians, that depend heavily on the electron density near atomic nuclei where relativistic and SOC effects are the strongest.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Solving the electronic structure problem for over 100000 atoms in real space

Using a real-space high-order finite-difference approach, we investigate the electronic structure of large spherical silicon nanoclusters. Within Kohn-Sham density functional theory and using pseudopotentials, we report the self-consistent field convergence of a system with over 100000 atoms: a Si 107,641 ⁢H 9,084 nanocluster with a diameter of 16 nm. Our approach uses Chebyshev-filtered subspace iteration to speed up the convergence of the eigenspace, and blockwise Hilbert space-filling curves to speed up sparse matrix-vector multiplications, all of which are implemented in the parsec code. For the largest system, we utilized 2048 nodes (114 688 cores) on the Frontera machine in the Texas Advanced Computing Center. Our quantitative analysis of the electronic structure shows how it gradually approaches its bulk counterpart as a function of nanocluster size. The band gap is enlarged due to quantum confinement in nanoclusters, but decreases as the system size increases, as expected. In conclusion, our work serves as a proof of concept for the capacity of the real-space approach in efficiently parallelizing very large calculations using high-performance computer platforms, which can straightforwardly be replicated in other systems with more than 10 5 atoms.

0-dimensional systems↗
DOE OSTI2021

Kohn-Sham Solver (KSSOLV) v2.0

KSSOLV is a MATLAB toolbox for solving Kohn-Sham density functional theory based electronic structure eigenvalue problems. It uses an object oriented features of MATLAB to represent atom, molecules, wavefunctions and Hamiltonians and their operations. It is designed to make it easier for users to prototype and test new algorithms for solving the Kohn-Sham problem. KSSOLV2.0 contains significant improvement over the original KSSOLV described in a paper published in ACM Transaction on Mathematical Software (attached). In addition to performing ground state calculation for small molecules, it can also perform geometry optimization for both molecules and solids. It uses standard pseudopotentials and implements local density approximation, generalized gradient approximation and hybrid functionals. Future releases will also include time-dependent DFT and post DFT calculations such as the GW quasi-particle energy calculation and Bethe-Salpeter equation solver for optical absorption.

Yang, Chao↗
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 OSTI2026

Custom-trained Machine-learning Interatomic Potentials: ZnCl2 Aqueous Solution

This dataset was generated using an iterative active-learning strategy implemented in the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials for aqueous ZnCl2 solutions. Each active-learning cycle consisted of three stages: training, exploration, and labeling. The initial training set combined configurations generated in this work from enhanced-sampling ab initio molecular dynamics simulations with configurations from a previously reported neural-network-potential study of aqueous ZnCl2. The enhanced-sampling ab initio molecular dynamics simulations involved Zn–Cl separation and the chloride coordination number around Zn²? as collective variables. These configurations served as the seed dataset. Subsequent active-learning cycles expanded the training set by identifying and labeling configurations that were poorly represented by the current models, thereby improving coverage of ion-association states and changes in local coordination and charge-state environments relevant to the solution free-energy landscape. For all selected configurations, single-point calculations of the total energies and atomic forces were performed within density functional theory using the CP2K Quickstep module. Reference calculations employed the revPBE-D3 and r2SCAN exchange-correlation functionals. Motivated by recent work on aqueous Zn²?, the main revPBE calculations omitted D3 dispersion contributions involving Zn²?, while retaining the D3 correction for water and chloride. For comparison, fully dispersion-corrected revPBE-D3 reference calculations were also performed, with D3 applied to all species, including Zn²?. Valence electrons were treated explicitly, while core electrons were represented using norm-conserving Goedecker–Teter–Hutter pseudopotentials. The wave functions were expanded using the mixed Gaussian-and-plane-wave scheme with TZV2P-MOLOPT basis sets for all elements and a 600 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent-field convergence was accelerated using the orbital-transformation and Direct Inversion in the Iterative Subspace algorithms, with a convergence threshold of 10?6. All single-point calculations were performed in periodic orthorhombic cells. The CELL_REF keyword in CP2K was used to define a fixed reference cell with a box length of 25 Å. This treatment ensured a consistent reference for configurations extracted from NpT trajectories with fluctuating cell dimensions. The resulting DFT energies and atomic forces constitute the ground-truth labels used to train the MLIPs. The resulting MLIP was trained for aqueous ZnCl2 solutions spanning concentrations from 0 to 30 molal and a broad pH range, from strongly acidic to strongly basic conditions. Representative examples of configurations included in the MLIP training dataset are provided below. These include 1) Representative configurations from the dataset labeled at the revPBE-D3 level, with D3 dispersion interactions involving Zn2+ excluded (revPBE-wo-D3). 2) Representative configurations from the dataset labeled at the fully dispersion-corrected revPBE-D3 level, with D3 interactions applied to all species, including Zn2+ (revPBE-D3). 3) Representative configurations from the dataset labeled at the r2SCAN level of theory (r2SCAN).

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2023

High-Temperature Thermodynamics of Uranium from Ab Initio Modeling

We present high-temperature thermodynamic properties for uranium in its γ phase (γ-U) from first-principles, relativistic, and anharmonic theory. The results are compared to CALPHAD modeling. The ab initio electronic structure is obtained from density-functional theory (DFT) that includes spin–orbit coupling and an added self-consistent orbital-polarization (OP) mechanism for more accurate treatment of magnetism. The first-principles method is coupled to a lattice dynamics scheme that is used to model anharmonic lattice vibrations, namely, Self-Consistent Ab Initio Lattice Dynamics (SCAILD). The methodology can be summarized in the acronym DFT + OP + SCAILD. Upon thermal expansion, γ-U develops non-negligible magnetic moments that are included for the first time in thermodynamic theory. The all-electron DFT approach is shown to model γ-U better than the commonly used pseudopotential method. In addition to CALPHAD, DFT + OP + SCAILD thermodynamic properties are compared with other ab initio and semiempirical modeling and experiments. Our first-principles approach produces Gibbs free energy that is essentially identical to CALPHAD. The DFT + OP + SCAILD heat capacity is close to CALPHAD and most experimental data and is predicted to have a significant thermal dependence due to the electronic contribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization

In this work, we present DFT-FE 1.0, building on DFT-FE 0.6 [Comput. Phys. Commun. 246, 106853 (2020)], to conduct fast and accurate large-scale density functional theory (DFT) calculations (reaching ~ 100,000 electrons) on both many-core CPU and hybrid CPU-GPU computing architectures. This work involves improvements in the real-space formulation—via an improved treatment of the electrostatic interactions that substantially enhances the computational efficiency—as well high-performance computing aspects, including the GPU acceleration of all the key compute kernels in DFT-FE. We demonstrate the accuracy by comparing the ground-state energies, ionic forces and cell stresses on a wide-range of benchmark systems against those obtained from widely used DFT codes. Further, we demonstrate the numerical efficiency of our implementation, which yields ~ 20× CPU-GPU speed-up by using GPU acceleration on hybrid CPU-GPU nodes. Notably, owing to the parallel-scaling of the GPU implementation, we obtain wall-times of 80–140 seconds for full ground-state calculations, with stringent accuracy, on benchmark systems containing ~ 6, 000 – 15,000 electrons.

pseudopotential↗
DOE OSTI2022

Investigation of the Failure of Marcus Theory for Hydrated Electron Reactions

Reactions of the hydrated electron with a wide variety of substrates have been found to exhibit unusually similar activation energies in a manner incompatible with Marcus electron transfer theory. Given the fundamental linear response assumption of Marcus theory, one possible explanation for this apparent failure is that the underlying free energy surfaces governing the reactions are not harmonic, i.e., hydrated electron structural fluctuations exhibit non-Gaussian behavior. In this work, we test this hypothesis by using simulations to calculate the hydrated electron vertical detachment energy distribution. Furthermore, we consider both cavity and non-cavity models for the hydrated electron, between which the actual hydrated electron behavior is expected to lie. Our results identify a possible origin for non-Gaussian behavior of the hydrated electron, but show that it is not of sufficient magnitude to explain the failure of Marcus theory to describe its reactions. Thus, other explanations must be sought.

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

Role of Superlattice Phonons in Charge Localization Across Quantum Dot Arrays

Understanding charge transport in semiconductor quantum dot (QD) assemblies is important for developing the next generation of solar cells and light-harvesting devices based on QD technology. One of the key factors that governs the transport in such systems is related to the hybridization between the QDs. Recent experiments have successfully synthesized QD molecules, arrays, and assemblies by directly fusing the QDs, with enhanced hybridization leading to high carrier mobilities and coherent band-like electronic transport. Here, in this work, we theoretically investigate the electron transfer dynamics across a finite CdSe-CdS core-shell QD array, considering up to seven interconnected QDs in one dimension. We find that, even in the absence of structural and size disorder, electron transfer can become localized by the emergent low-frequency superlattice vibrational modes when the connecting neck between QDs is narrow. On the other hand, we also identify a regime where the same vibrational modes facilitate coherent electron transport when the connecting necks are wide. Overall, we elucidate the crucial effects of electronic and superlattice symmetries and their couplings when designing high-mobility devices based on QD superlattices.

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

Unraveling Exciton Trap Dynamics and Nonradiative Loss Pathways in Quantum Dots via Atomistic Simulations

Surface defects in colloidal quantum dots are a major source of nonradiative losses, yet the microscopic mechanisms underlying exciton trapping and recombination remain elusive. Here, we develop a model Hamiltonian based on atomistic electronic calculations to investigate exciton dynamics in CdSe/CdS core/shell QDs containing a single hole trap introduced by an unpassivated sulfur atom. By systematically varying the defect depth and reorganization energy, we uncover how defect-induced excitonic states mediate energy relaxation pathways. Our simulations reveal that a single localized defect can induce a rich spectrum of excitonic states, leading to multiple dynamical regimes, from slow, energetically off-resonant trapping to fast, cascaded relaxation through in-gap defect states. Crucially, we quantify how defect-induced polaron shifts and exciton-phonon couplings govern the balance between efficient radiative emission and rapid nonradiative decay. These insights clarify the microscopic origin of defect-assisted loss channels and suggest pathways for tailoring QD optoelectronic properties via surface and defect engineering.

Defects↗
DOE OSTI2021

Induced anomalous Hall effect of massive Dirac fermionsin ZrTe 5 and HfTe 5 thin flakes

Research on the anomalous Hall effect (AHE) has been lasting for a century to make clear the underlying physical mechanism. Generally, the AHE appears in magnetic materials, in which the extrinsic process related to scattering effects and intrinsic contribution connected with Berry curvature are crucial. Recently, AHE has been counterintuitively observed in nonmagnetic topological materials and attributed to the existence of Weyl points. However, the Weyl point scenario would lead to unsaturated AHE even in large magnetic fields and contradicts the saturation of AHE in several tesla (T) in experiments. In this work, we investigate the Hall effect of ZrTe 5 and HfTe 5 thin flakes in static ultrahigh magnetic fields up to 33 T. We find the AHE saturates to 55(70)Ω 1 cm 1 for ZrTe 5 (HfTe 5 ) thin flakes above ~10T. Combining detailed magnetotransport experiments and Berry curvature calculations, we clarify that the splitting of massive Dirac bands without Weyl points can be responsible for AHE in nonmagnetic topological materials ZrTe 5 and HfTe 5 thin flakes. This model can identify our thin flake samples to be weak topological insulators and serve as a tool to probe the band structure topology in topological materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Structural Principles in Liquids and Glasses: Bottom-Up or Top-Down

The conventional approach to elucidate the atomic structure of liquid and glass is to start with local structural units made of several atoms, and to use them as building blocks to form a global structure, the bottom-up approach. We propose to add an alternative top-down approach in which we start with a global high-temperature gas state and then apply interatomic potentials to all atoms at once. This causes collective density wave instability in all directions with the same wavelength. These two driving forces, local and global, are in competition and are mutually frustrated. The final structure is determined through the compromise of frustration between these two, which creates the medium-range-order. This even-handed approach on global and local potential energy landscapes explains the distinct natures of short-range order and medium-range order, and strong temperature dependence of various properties of liquid.

36 MATERIALS SCIENCE↗