Search NASA⌕ Search

SEARCH · Search NASA

Results for “THEORIES - IMPACT”

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 343 records · Page 19

Steam Reforming on Transition-metal Carbides from Density-functional Theory

A screening study of the steam reforming reaction on clean and oxygen covered early transition-metal carbides surfaces is performed by means of density-functional theory calculations. It is found that carbides provide a wide spectrum of reactivities, from too reactive via suitable to too inert. Several molybdenum-based systems are identified as possible steam reforming catalysts. The findings suggest that carbides provide a playground for reactivity tuning, comparable to the one for pure metals.

36 MATERIALS SCIENCE↗

Gauge invariance and Ward identities in nonlinear response theory

We present a formal analysis of nonlinear response functions in terms of correlation functions in real- and imaginary-time domains. In particular, we show that causal nonlinear response functions, expressed in terms of nested commutators in real time, can be obtained from the analytic continuation of time-ordered response functions, which are more easily amenable to diagrammatic calculation. This generalizes the well-known result of linear response theory. We then use gauge invariance arguments to derive exact relations between second-order response functions in density and current channels. These identities, which are non-perturbative in the strength of inter-particle interactions, allow us to establish exact connections between nonlinear optics calculations done in different electromagnetic gauges.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Probing the Structures and Electronic Properties of Dual-Phosphorus-Doped Gold Cluster Anions (AunP-2, n = 1–8): A Density functional Theory Investigation

The geometries of gold clusters doped with two phosphorus atoms, (AunP-2, n = 1–8) were investigated using density functional theory (DFT) methods. Various two-dimensional (2D) and three-dimensional (3D) structures of the doped clusters were studied. The results indicate that the structures of dual-phosphorus-doped gold clusters exhibit large differences from those of pure gold clusters with small cluster sizes. In our study, as for Au6P-2, two cis–trans isomers were found. The global minimum of Au8P-2 presents a similar configuration to that of Au-20, a pyramid-shaped unit, and the potential novel optical and catalytic properties of this structure warrant further attention. The higher stability of AunP-2 clusters relative to Au-n+2 (n = 1–8) clusters was verified based on various energy parameters, and the results indicate that the phosphorus atom can improve the stabilities of the gold clusters. We then explored the evolutionary path of (n = 1–8) clusters. We found that AunP-2 clusters exhibit the 2D–3D structural transition at n = 6, which is much clearer and faster than that of pure gold clusters and single-phosphorus-doped clusters. The electronic properties of AunP-2 (n = 1–8) were then investigated. The photoelectron spectra provide additional fundamental information on the structures and molecular orbitals shed light on the evolution of AunP-2 (n = 1–8). Natural bond orbital (NBO) described the charge distribution in stabilizing structures and revealed the strong relativistic effects of the gold atoms.

Environmental Molecular Sciences Laboratory↗

Fast plane wave density functional theory molecular dynamics calculations on multi-GPU machines

Plane wave pseudopotential (PWP) density functional theory (DFT) calculation is the most widely used method for material simulations, but its absolute speed stagnated due to the inability to use large scale CPU based computers. By a drastic redesign of the algorithm, and moving all the major computation parts into GPU, we have reached a speed of 12 s per molecular dynamics (MD) step for a 512 atom system using 256 GPU cards. This is about 20 times faster than the CPU version of the code regardless of the number of CPU cores used. Our tests and analysis on different GPU platforms and configurations shed lights on the optimal GPU deployments for PWP-DFT calculations. An 1800 step MD simulation is used to study the liquid phase properties of GaInP.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigation of Structural Evolution of Li 1.1 V 3 O 8 by In Situ X-ray Diffraction and Density Functional Theory Calculations

Li 1.1 V 3 O 8 is a promising cathode material for Li-ion batteries, because of its high theoretical capacity (362 mAh g-1) and good rate performance. In this study, the structural evolution of Li 1.1 V 3 O 8 material during electrochemical dis(charge) processes was investigated using a combination of theoretical calculations and experimental data. Density functional theory (DFT) was used to predict the intermediate structures at various lithiation states, as well as the stability of major phases. In order to validate these predictions, in situ X-ray diffraction (XRD) data was collected operando, allowing for the phase transformations to be monitored under current load and eliminating the possibilities of structural relaxation processes and environmental oxidation. Rietveld refinement was performed to fit the diffraction data with the DFT-derived structures and to analyze the fractions of major phases as a function of dis(charge). The DFT calculations identified three stable states that were validated by the in situ XRD result: a Li-poor α-phase (Li1), a Li-rich α-phase (Li2.5), and a β-phase (Li4). The DFT-predicted particle shape based on the surface energy of the (100), (001), and (010) planes rationalized the preferential orientation of Li 1.1 V 3 O 8 particles along the [010] direction in the electrode. Furthermore, the onset and offset of the α → β transition, as well as the phase fractions of α and β determined via in situ XRD, related well with the DFT-derived relative stability of each phase. Thus, by integrating DFT calculations with experimental work, this work provides a thorough understanding of the structural transformations in Li 1.1 V 3 O 8 during electrochemical dis(charge).

25 ENERGY STORAGE↗

Dipole Moment Calculations Using Multiconfiguration Pair-Density Functional Theory and Hybrid Multiconfiguration Pair-Density Functional Theory

The dipole moment is the molecular property that most directly indicates molecular polarity. The accuracy of computed dipole moments depends strongly on the quality of the calculated electron density, and the breakdown of single-reference methods for strongly correlated systems can lead to poor predictions of the dipole moments in those cases. Here, we derive the analytical expression for obtaining the electric dipole moment by multiconfiguration pair density functional theory (MCPDFT), and we assess the accuracy of MC-PDFT for predicting dipole moments at equilibrium and nonequilibrium geometries. We show that MC-PDFT dipole moment curves have reasonable behavior even for stretched geometries, and they significantly improve upon the CASSCF results by capturing more electron correlation. Herein, the analysis of a dataset consisting of 18 first-row transition metal diatomics and 6 main-group polyatomic molecules with multireference character suggests that MC-PDFT and its hybrid extension (HMC-PDFT) perform comparably to CASPT2 and MRCISD+Q methods and have a mean unsigned deviation of 0.2–0.3 D with respect to the best available dipole moment reference values. We explored the dependence of the predicted dipole moments upon the choice of the on-top density functional and active space, and we recommend the tPBE and hybrid tPBE0 on-top choices for the functionals combined with the moderate correlated participating orbital scheme for selecting the active space. With these choices, the mean unsigned deviations (in debyes) of the calculated equilibrium dipole moments from the best estimates are 0.77 for CASSCF, 0.29 for MC-PDFT, 0.24 for HMCPDFT, 0.28 for CASPT2, and 0.25 for MRCISD+Q. These results are encouraging because the computational cost of MC-PDFT or HMC-PDFT is largely reduced compared to the CASPT2 and MRCISD+Q methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pseudodiagonalization Method for Accelerating Nonlinear Subspace Diagonalization in Density Functional Theory

In density functional theory, each self-consistent field (SCF) nonlinear step updates the discretized Kohn-Sham orbitals by solving a linear eigenvalue problem. The concept of pseudodiagonalization is to solve this linear eigenvalue problem approximately, and specifically utilizing a method involving a small number of Jacobi rotations that takes advantage of the good initial guess to the solution given by the approximation to the orbitals from the previous SCF iteration. The approximate solution to the linear eigenvalue problem can be very rapid, particularly for those steps near SCF convergence. Here, we adapt pseudodiagonalization to finite-temperature and metallic systems, where partially-occupied orbitals must be individually resolved with some accuracy. We apply pseudodiagonalization to the subspace eigenvalue problem that arises in Chebyshev-filtered subspace iteration. In tests on metallic and other systems for a range of temperatures, we show that pseudodiagonalization achieves similar rates of SCF convergence to exact diagonalization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multicomponent Cholesky Decomposition: Application to Nuclear–Electronic Orbital Theory

The Cholesky decomposition technique is commonly used to reduce the memory requirement for storing two-particle repulsion integrals in quantum chemistry calculations that use atomic orbital bases. However, when quantum methods use multicomponent bases, such as nuclear–electronic orbitals, additional challenges are introduced due to asymmetric two-particle integrals. This work proposes several multicomponent Cholesky decomposition methods for calculations using nuclear–electronic orbital density functional theory. To analyze the errors in different Cholesky decomposition components, benchmark calculations using water clusters are carried out. The largest benchmark calculation is a water cluster (H 2 O) 27 where all 54 protons are treated quantum mechanically. Furthermore, this study provides energetic and complexity analyses to demonstrate the accuracy and performance of the proposed multicomponent Cholesky decomposition method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Extending Orbital-Optimized Density Functional Theory to L-Edge XPS and Beyond: Spin–Orbit Coupling via Nonorthogonal Quasi-Degenerate Perturbation Theory

Quantum mechanical calculations of core electron binding energies (CEBEs) are relevant to interpreting X-ray photoelectron spectroscopy (XPS). Orbital-optimized density functional theory (OO-DFT) accurately predicts K-edge CEBEs but is challenged by the presence of significant spin–orbit coupling (SOC) at L- and higher edges involving inner-shell orbitals with nonzero angular momentum. Here, to extend OO-DFT to L-edges and higher, our method utilizes scalar-relativistic, spin-restricted open-shell OO-DFT to construct a minimal, quasi-degenerate basis of core-hole states corresponding to a chosen inner-shell (e.g., ionizing all six possible 2p spin orbitals). Nonorthogonal configuration interaction (NOCI) is then used to obtain the matrix elements of the full Hamiltonian including SOC in this quasi-degenerate model space of determinants. Using a screened 1-electron SOC operator parametrized with the Dirac-Coulomb-Breit (DCB) Hamiltonian results in doublet splitting (DS) values for third row elements that are nearly in quantitative agreement with experiment. The resulting NOCI eigenvalues are shifted by the average of the (scalar) OO-DFT CEBEs to yield CEBEs (split by SOC) corrected for dynamic correlation. Comparing calculations on gas phase molecules with experimental results establishes that NO-QDPT with the SCAN functional (NO-QDPT/SCAN), using the DCB screened 1-electron SOC operator is accurate to about 0.2 eV for L-edge CEBEs of molecules containing third row atoms. However, this NO-QDPT approach becomes less accurate for fourth-row elements starting in the middle of the 3d transition metal series, with errors increasing as atomic number increases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solid Oganesson via a Many-Body Interaction Expansion Based on Relativistic Coupled-Cluster Theory and from Plane-Wave Relativistic Density Functional Theory

In this work, many-body potentials up to fourth order in the expansion are constructed using non-relativistic, scalar-relativistic and relativistic coupled-cluster theory to accurately describe the interaction between superheavy oganesson atoms. The obtained distance dependent energy values were fitted to extended two-body Lennard-Jones and three-body Axilrod-Teller-Muto potentials, with the fourth-order term treated through a classical long-range Drude dipole interaction model. From these interac tion potentials, spectroscopic constants for the oganesson dimer and solid state properties were obtained. Furthermore, these high-level results are compared to scalar relativistic and two-component plane-wave DFT calculations based on a tailor-made projector augmented-wave pseudopotential (PAW-PP) and newly derived parameters for Grimme’s dispersion correction. It is shown that the functionals PBE-D3(BJ), PBEsol and in particular SCAN provide excellent agreement with the many-body ref erence for solid Og. Finally, the results for Og are compared and related to the lighter rare gas elements and periodic trends are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗