Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum mechanical calculations”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Phase diagram of the TIP4P/Ice water model by enhanced sampling simulations

In this work, we studied the phase diagram for the TIP4P/Ice water model using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice–liquid free energy differences from biased coexistence simulations that reversibly sample the melting and growth of layers of ice. We computed a total of 19 melting points for five different ice polymorphs, which are in excellent agreement with the melting lines obtained from the integration of the Clausius–Clapeyron equation. For proton-ordered and fully proton-disordered ice phases, the results are in very good agreement with previous calculations based on thermodynamic integration. For the partially proton-disordered ice III, we find a large increase in stability that is in line with previous observations using direct coexistence simulations for the TIP4P/2005 model. This issue highlights the robustness of the approach employed here for ice polymorphs with diverse degrees of proton disorder. Our approach is general and can be applied to the calculation of other complex phase diagrams.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Strong influence of nonmagnetic ligands on the momentum-dependent spin splitting in antiferromagnets

Here, recent studies show the nonrelativistic antiferromagnetic ordering could generate momentum-dependent spin splitting analogous to the Rashba effect but free from the requirement of relativistic spin-orbit coupling. Whereas the classification of such compounds can be illustrated by different spin-splitting prototypes (SSTs) from symmetry analysis and density-functional-theory calculations, the huge variation in chemical bonding and structures of these diverse compounds possibly clouds the issue of how much of the variation in spin splitting can be traced back to the symmetry-defined characteristics, rather to the underlining chemical and structural diversity. The alternative model Hamiltonian approaches do not confront the issues of chemical and structural complexity but often consider only the magnetic sublattice, dealing with the all-important effects of the nonmagnetic ligands via renormalizing the interactions between the magnetic sites. To this end, we constructed a DFT model Hamiltonian that allows us to study SSTs at constant chemistry while retaining the realistic atomic-scale structure including ligands. This is accomplished by using a single, universal magnetic skeletal lattice (Ni 2+ ions in rocksalt NiO) and designing small displacements of the nonmagnetic (oxygen) sublattice which produce, by design, the different SST magnetic symmetries. We show that (i) even similar crystal structures having very similar band structures can lead to contrasting behavior of spin splitting vs momentum, and (ii) even subtle deformations of the nonmagnetic ligand sublattice could cause a giant spin splitting in AFM-induced SST. This is a paradigm shift relative to the convention of modeling magnets without considering the nonmagnetic ligand that mediates indirect magnetic interaction (e.g., superexchange).

36 MATERIALS SCIENCE↗
DOE OSTI2023

An over-approximation of entropy for elemental multiconfigurational ground state electronic structures

The ability to predict the chemical and physical properties of a material is directly related to the structure and interactions of its electrons. For materials comprised of f-block elements (the lanthanides and actinides found in the last two rows of the periodic table), the complexity of electronic structure has presented great difficulty in understanding, modeling, and predicting material properties. The complexity of multiconfigurational ground state electronic structures is illustrated herein by the combinatorics of electron permutations within individual and cumulative occupancy configurations. A non-integer orbital occupancy representation of multiconfigurational ground states is described for superposition mixing between multiple near-energy degenerate occupancy configurations and generalized in such a way that established ground states are returned by approximation for elements with less-complex electronic structures. By considering the occupancy configurations as statistical mechanics macrostates, and the permutations of electrons as statistical mechanics microstates within those macrostates, an over-approximation of entropy for multiconfigurational elemental ground state electronic structures has been calculated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Quantum Information for Particle Theorists

Lectures given at the Theoretical Advanced Study Institute (TASI 2020), 1-26 June 2020. The topics covered include quantum circuits, entanglement, quantum teleportation, Bell inequalities, quantum entropy and decoherence, classical versus quantum measurement, the area law for entanglement entropy in quantum field theory, and simulating quantum field theory on a quantum computer. Along the way we confront the fundamental sloppiness of how we all learned (and some of us taught) quantum mechanics in college. Links to a Python notebook and Mathematica notebooks will allow the reader to reproduce and extend the calculations, as well as perform five experiments on a quantum simulator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Salt-flux synthesis, crystal structure and theoretical characterization of Rb 0.74 Ga 6.62 Ti 0·38 O 11

Here single crystals of Rb 0.74 Ga 6.62 Ti 0·38 O 11 (RGTO) were grown from a mixed RbCl–RbF flux at 850 °C. The compound crystallizes in the RbGa 7 O 11 structure type, which is reminiscent of the hollandite and β-Ga 2 O 3 structure types. RGTO crystallizes in the monoclinic space group P2/m with lattice parameters a = 8.3355 (8) Å, b = 3.0286 (3) Å, c = 9.5028 (9) Å, and β = 114.620 (3)°. The crystal structure of RGTO is comprised of GaO 6 and mixed (Ga/Ti)O 6 octahedra and GaO 4 tetrahedra connected in a complex three-dimensional, anionic framework exhibiting eight-sided channels that are occupied by disordered Rb cations required for charge balance. First-principles calculations in the form of density functional theory were performed, which indicated the complex to be a charge transfer semiconductor.

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

Wormhole calculus, replicas, and entropies

We investigate contributions of spacetime wormholes, describing baby universe emission and absorption, to calculations of entropies and correlation functions, for example those based on the replica method. We find that the rules of the “wormhole calculus”, developed in the 1980s, together with standard quantum mechanical prescriptions for computing entropies and correlators, imply definite rules for limited patterns of connection between replica factors in simple calculations. These results stand in contrast with assumptions that all topologies connecting replicas should be summed over, and call into question the explanation for the latter. In a “free” approximation baby universes introduce probability distributions for coupling constants, and we review and extend arguments that successive experiments in a “parent” universe increasingly precisely fix such couplings, resulting in ultimately pure evolution. Once this has happened, the nontrivial question remains of how topology-changing effects can modify the standard description of black hole information loss.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Computational Spectroscopy of the Cr–Cr Bond in Coordination Complexes

In this work, we report the accurate computational vibrational analysis of the Cr–Cr bond in dichromium complexes using second-order multireference complete active space methods (CASPT2), allowing direct comparison with experimental spectroscopic data both to facilitate interpreting the low-energy region of the spectra and to provide insights into the nature of the bonds themselves. Recent technological development by the authors has realized such computation for the first time. Accurate simulation of the vibrational structure of these compounds has been hampered by their notorious multiconfigurational electronic structure that yields bond distances that do not correlate with bond order. Some measured Cr–Cr vibrational stretching modes, ν(Cr 2 ), have suggested weaker bonding, even for so-called ultrashort Cr–Cr bonds, while others are in line with the bond distance. Here, we optimize geometries and compute ν(Cr 2 ) with CASPT2 for three well-characterized complexes, Cr 2 (O 2 CCH 3 ) 4 (H 2 O) 2 , Cr 2 (mhp) 4 , and Cr 2 (dmp) 4 . We obtain CASPT2 harmonic ν(Cr 2 ) modes in good agreement with experiment at 282 cm –1 for Cr 2 (mhp) 4 and 353 cm –1 for Cr 2 (dmp) 4 , compute 50 Cr and 54 Cr isotope shifts, and demonstrate that the use of the so-called IPEA shift leads to improved Cr–Cr distances. Additionally, normal mode sampling was used to estimate anharmonicity along ν(Cr 2 ), leading to an anharmonic mode of 272 cm –1 for Cr 2 (mhp) 4 and 333 cm –1 for Cr 2 (dmp) 4 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Effect of Geometry, Spin, and Orbital Optimization in Achieving Accurate, Correlated Results for Iron–Sulfur Cubanes

Iron–sulfur clusters comprise an important functional motif in the catalytic centers of biological systems, capable of enabling important chemical transformations at ambient conditions. This remarkable capability derives from a notoriously complex electronic structure that is characterized by a high density of states that is sensitive to geometric changes. The spectral sensitivity to subtle geometric changes has received little attention from correlated, large active space calculations, owing partly to the exceptional computational complexity for treating these large and correlated systems accurately. Here, to provide insight into this aspect, we report the first Complete Active Space Self Consistent Field (CASSCF) calculations for different geometries of the [Fe(II/III) 4 S 4 (SMe) 4 ] -2 clusters using two complementary, correlated solvers: spin-pure Adaptive Sampling Configuration Interaction (ASCI) and Density Matrix Renormalization Group (DMRG). We find that the previously established picture of a double-exchange driven magnetic structure, with minute energy gaps (<1 mHa) between consecutive spin states, has a weak dependence on the underlying geometry. However, the spin gap between the singlet and the spin state 2S + 1 = 19, corresponding to a maximal number of Fe-d electrons being unpaired and of parallel spin, is strongly geometry dependent, changing by a factor of 3 upon slight deformations that are still within biologically relevant parameters. The CASSCF orbital optimization procedure, using active spaces as large as 86 electrons in 52 orbitals, was found to reduce this gap compared to typical mean-field orbital approaches. Our results show the need for performing large active space calculations to unveil the challenging electronic structure of these complex catalytic centers and should serve as accurate starting points for fully correlated treatments upon inclusion of dynamical correlation outside the active space.

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

Study of Self-Interaction Errors in Density Functional Calculations of Magnetic Exchange Coupling Constants Using Three Self-Interaction Correction Methods

We examine the role of self-interaction errors (SIE) removal on the evaluation of magnetic exchange coupling constants. In particular we analyze the effect of scaling down the self-interaction-correction (SIC) for three non-empirical density functional approximations (DFAs) namely, the local spin density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation, and recent SCAN family of meta-GGA functionals. To this end, we employ three one-electron SIC methods: Perdew-Zunger [Perdew, J. P.; Zunger, A. Phys. Rev. B, 1981, 23, 5048] SIC, the orbitalwise scaled SIC method [Vydrov, O. A. et al., J. Chem. Phys. 2006, 124, 094108], and the recent local scaling method [Zope, R. R. et al., J. Chem. Phys. 2019, 151, 214108]. We compute the magnetic exchange coupling constants using the spin projection and non projection approaches for sets of molecules composed of dinuclear and polynuclear H-He models, organic radical molecules, and chlorocuprate, and compare these results against accurate theories and experiment. Our results show that for the systems that mainly consist of single electron regions, PZSIC performs well but for more complex organic systems and the chlorcuprates, an overcorrecting tendency of PZSIC combined with the DFAs utilized in this work is more pronounced, and in such cases LSIC with kinetic energy density ratio performs better than PZSIC. Furthermore, analysis of the results in terms of SIC corrections to the density and to the total energy shows that both density and energy correction are required to obtain an improved prediction of magnetic exchange couplings.

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

Realistic Ion Dynamics through Charge Renormalization in Nonaqueous Electrolytes

While many practically important electrolytes contain lithium ions, interactions of these ions are particularly difficult to probe experimentally because of their small X-ray and neutron scattering cross sections and large neutron absorption cross sections. Molecular dynamics (MD) is a powerful tool for understanding the properties of nonaqueous electrolyte solutions from the atomic level, but the accuracy of this computational method crucially depends on the physics built into the classical force field. In this study, we demonstrate that several force fields for lithium bistriflimide (LiTFSI) in acetonitrile yield a solution structure that is consistent with the neutron scattering experiments, yet these models produce dramatically different ion dynamics in solution. Such glaring discrepancies indicate that inadequate representation of long-range interactions leads to excessive ionic association and ion-pair clustering. We show that reasonable agreement with the experimental observations can be achieved by renormalization of the ion charges using a “titration” method suggested herewith. This simple modification produces realistic concentration dependencies for ionic diffusion and conductivity in <2 M solutions, without loss in quality for simulation of the structure.

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

Exact-Factorization-Based Surface Hopping for Multistate Dynamics

A surface-hopping algorithm recently derived from the exact factorization approach, SHXF, introduces an additional term in the electronic equation of surface hopping that couples electronic states through the quantum momentum. Furthermore, this term not only provides a first-principles description of decoherence, but here we show it is crucial to accurately capture nonadiabatic dynamics when more than two states are occupied at any given time. Using a vibronic coupling model of the uracil cation, we show that the lack of this term in traditional surface-hopping methods, including those with decoherence corrections, leads to failure to predict the dynamics through a three-state intersection, while SHXF performs similarly to the multiconfiguration time-dependent Hartree quantum dynamics benchmark.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A first-principles study of hydrogen surface coverage on δ-Pu (100), (111), and (110) surfaces

Hydriding corrosion of plutonium leads to surface cracking, pitting, and ultimately structural failure. Laboratory experiments demonstrate that hydriding begins on the surface or near the subsurface of plutonium. However, there has not yet been a systematic evaluation of hydrogen surface coverage on plutonium. In this work, we compute the surface energies of the low facet surfaces of face-centered cubic δ-Pu. The adsorption free energies of expected hydrogen structures at low and high coverage are presented along with the likely progression for filling sites as the H 2 partial pressure increases. Implications for near-equilibrium pressure hydride nucleation and non-equilibrium millibar pressure hydriding are discussed.

36 MATERIALS SCIENCE↗
NASA NTRS1966

Pictures of molecular orbitals.

Molecular orbital electron charge density pictures, noting representation close to Hartree-Fock calculation method

MOLECULAR ORBIT↗
DOE OSTI2020

Shell-model study of calcium isotopes toward their drip line

Here, we report in this paper a study in terms of the nuclear shell model about the location of the calcium isotopes drip line. The starting point is considering the realistic two-body potential derived by Entem and Machleidt within chiral perturbation theory at next-to-next-to-next-to-leading order (N 3 LO), as well as a chiral three-body force at next-to-next-to-leading order (N 2 LO) whose structure and low-energy constants are consistent with the two-body potential. Then we construct the effective single-particle energies and residual interaction needed to diagonalize the shell-model Hamiltonian. The calculated two-neutron separation energies agree nicely with experiment until 56 Ca, which is the heaviest isotope whose mass has been measured, and do not show any sign of two-neutron emission until 70 Ca . We discuss the role of the choice of the model space in determining the neutron drip line, and also the dependence of the results on the parameters of the shell-model Hamiltonian.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Laser-free trapped-ion entangling gates with simultaneous insensitivity to qubit and motional decoherence

The dominant error sources for state-of-the-art laser-free trapped-ion entangling gates are decoherence of the qubit state and the ion motion. The effect of these decoherence mechanisms can be suppressed with additional control fields or with techniques that have the disadvantage of reducing gate speed. In this study, we propose using a near-motional-frequency magnetic field gradient to implement a laser-free gate that is simultaneously resilient to both types of decoherence, does not require additional control fields, and has a relatively small cost in gate speed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Direct and indirect spin current generation and spin-orbit torques in ferromagnet/nonmagnet/ferromagnet trilayers

Spin-orbit torques in ferromagnet/nonmagnet/ferromagnet trilayers are studied using a combination of symmetry analysis, circuit theory, semiclassical simulations, and first-principles calculations using the nonequilibrium Green's function method with supercell disorder averaging. Here, we focus on unconventional processes involving the interplay between the two ferromagnetic layers, which are classified into direct and indirect mechanisms. The direct mechanism involves spin current generation by one ferromagnetic layer and its subsequent absorption by the other. In the indirect mechanism, the in-plane spin-polarized current from one ferromagnetic layer “leaks” into the other layer, where it is converted into an out-of-plane spin current and reabsorbed by the original layer. The direct mechanism results in a predominantly dampinglike torque, which damps the magnetization towards a certain direction 𝐬 𝑑 . The indirect mechanism results in a predominantly fieldlike torque with respect to a generally different direction 𝐬 𝑓 . Similarly to the current-in-plane giant magnetoresistance, the indirect mechanism is only active if the thickness of the nonmagnetic spacer is smaller than or comparable to the mean free path. Numerical calculations for a semiclassical model based on the Boltzmann equation confirm the presence of both direct and indirect mechanisms of spin current generation. First-principles calculations reveal sizable unconventional spin-orbit torques in Co/Cu/Co, Py/Cu/Py, and Co/Pt/Co trilayers and provide strong evidence of indirect spin current generation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗