Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix elements”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Local Pair Natural Orbital-Based Coupled-Cluster Theory through Full Quadruples (DLPNO–CCSDTQ)

In this work, we implement a local pair natural orbitalbased coupled-cluster method through the full treatment of quadruple excitations (CCSDTQ). The domain-based local pair natural orbital (DLPNO) approach, which has successfully been applied to lower levels of coupled-cluster theory, is utilized in our algorithm, and thus our algorithm is called DLPNO-CCSDTQ. For simplicity in the working equations and in the implementation, we t 1 -dress the twoelectron integrals as well as Fock matrix elements. Our method can recover CCSDTQ-CCSDT and CCSDTQ-CCSDT(Q) energy differences on the order of 0.01−0.05 kcal mol −1 , even at a loose quadruples natural orbital (QNO) occupation number cutoff of 3.33 × 10 −6 . To highlight the capabilities of our code and its potential future applications, we showcase computations that would be intractable with canonical CCSDTQ, such as the benzene dimer, (H 2 O) 17 , and adamantane. With sufficient computing resources, computations up to 15 heavy atoms (40 atoms overall) may be feasible for fully bonded 3D systems.

Cluster chemistry↗
DOE OSTI2020

A Multireference Quantum Krylov Algorithm for Strongly Correlated Electrons

We introduce a multireference selected quantum Krylov (MRSQK) algorithm suitable for quantum simulation of many-body problems. MRSQK is a low-cost alternative to the quantum phase estimation algorithm that generates a target state as a linear combination of non-orthogonal Krylov basis states. This basis is constructed from a set of reference states via real-time evolution; thus, avoiding the numerical optimization of parameters. An efficient algorithm for the evaluation of the off-diagonal matrix elements of the overlap and Hamiltonian matrices is discussed and a selection procedure is introduced to identify a basis of orthogonal references that ameliorates the linear dependency problem. Preliminary benchmarks on linear H 6 , H 8 , and BeH 2 indicate that MRSQK can predict the energy of these systems accurately using very compact Krylov bases.

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

Single-Photon-Induced Electron Spin Polarization of Two Exchange-Coupled Stable Radicals

Transient electron paramagnetic resonance spectroscopy has been used to probe photoinduced electron spin polarization of a stable exchange-coupled organic biradical in a Pt(II) complex comprising 4,4'-di-tert-butyl-2,2'-bipyridine (bpy) and 3,6-bis(ethynyl-para-phenyl-nitronyl nitroxide)-o-catecholate (CAT(o-C$\equiv$C-Ph-NN) 2 ). Photoexcitation results in four unpaired spins in excited states of this complex, with spins being localized on each of the two radicals, CAT •+ and bpy •– . The four spins are all exchange-coupled in these excited states, and an off-diagonal matrix element in the CAT •+ -NN exchange allows for exchange-enhanced intersystem crossing to the 3 T 1a state, which possesses (bpy •– )Pt(CAT •+ ) chromophoric triplet character. Further, fast mixing between this 3 T 1a state and other thermally accessible excited LL'CT state(s) followed by fast relaxation provides spin polarization of the exchange-coupled NN radicals in the 3 S 0 ground state of the complex. Our results demonstrate that well-defined quantum states of a ground-state biradical can be initialized with single-photon excitation and have the potential for further spin manipulation directed toward quantum information science applications.

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

Photoredox Catalysis with Spin Magnetic Field Effects: A Scheme for Chiral Resolution

Quantities that break both mirror symmetry and time-reversal symmetry, such as the orbital angular momentum, are known to connect molecular chirality with an applied magnetic field. This concept has led to observations such as magneto-chiral dichroism and chirality-induced spin selectivity (CISS). However, being small, these effects often require additional amplification procedures such as flow chemistry to achieve bulk enantioseparation. In this work, we demonstrate how the magnetic field effect on photogenerated radical pairs, which also breaks time-reversal symmetry, can be harnessed for enantiopurification. Fundamental to this process is the collective decay of the singlet and triplet radical-pair states made possible by an applied magnetic field. Because opposite enantiomers exhibit spin-orbit coupling matrix elements of opposite signs, the singlet and triplet decay channels interfere constructively in one enantiomer. Meanwhile, molecules of the other enantiomer are funnelled into the first enantiomer through excited-state chirality inversion, achieving (dynamic kinetic) chiral resolution. Using an axially chiral binaphthyl derivative and a borane photosensitizer as prototype, we predict an appreciable enantiomeric excess (e.e.) of 90% to be possible at steady state, attained within hundreds of milliseconds when irradiated by a laser. Importantly, our analytical results showcase regimes of perfect enantioselectivity (100% e.e.), accessible by further chemical optimization of the photosensitizer for which general strategies are discussed. Altogether, this work illustrates a so-far untapped but powerful control knob for photoredox catalysis based on spin chemistry principles.

Magnetic properties↗
DOE OSTI2021

van der Waals driven anharmonic melting of the 3D charge density wave in VSe 2

Understanding of charge-density wave (CDW) phases is a main challenge in condensed matter due to their presence in high-T c superconductors or transition metal dichalcogenides (TMDs). Among TMDs, the origin of the CDW in VSe 2 remains highly debated. Here, by means of inelastic x-ray scattering and first-principles calculations, we show that the CDW transition is driven by the collapse at 110 K of an acoustic mode at q CDW = (2.25 0 0.7) r.l.u. The softening starts below 225 K and expands over a wide region of the Brillouin zone, identifying the electron-phonon interaction as the driving force of the CDW. This is supported by our calculations that determine a large momentum-dependence of the electron-phonon matrix-elements that peak at the CDW wave vector. Our first-principles anharmonic calculations reproduce the temperature dependence of the soft mode and the T CDW onset only when considering the out-of-plane van der Waals interactions, which reveal crucial for the melting of the CDW phase. The nature of the charge density wave transition in VSe 2 is still debated. Here, the authors demonstrate that the transition is mainly driven by electron-phonon interactions, despite the presence of the Fermi-surface nesting, and that Wan-der-Waals forces are responsible for melting of the charge density wave order.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Bringing weak transitions to light

Abstract Weak transitions between quantum states are of fundamental importance for a broad range of phenomena from analytical biochemistry to precision physics, but generally challenge experimental detection. Due to their small cross sections scaling with the absolute square of their transition matrix elements, spectroscopic measurements often fail in particular in the presence of competing background processes. Here we introduce a general concept to break this scaling law and enhance the transition probability by exploiting a stronger laser-coupled pathway to the same excited state. We demonstrate the concept experimentally by attosecond transient absorption spectroscopy in helium atoms. The quasi-forbidden transitions from the ground state 1s 2 to the weakly coupled doubly excited 2p3dandsp 2,4− states are boosted by an order of magnitude. Enhancing single-photon-suppressed transitions can find widespread applicability, from spectral diagnostics of complex molecules in life and chemical sciences to precision spectroscopy of weak transitions in metastable atomic nuclei in the search for new physics.

Science & Technology - Other Topics↗
DOE OSTI2023

Substrate effects on spin relaxation in two-dimensional Dirac materials with strong spin-orbit coupling

Understanding substrate effects on spin dynamics and relaxation is of key importance for spin-based information technologies. However, the key factors that determine such effects, in particular for materials with strong spin-orbit coupling (SOC), have not been well understood. Here we performed first-principles real-time density-matrix dynamics simulations with SOC and the electron phonon and electron-impurity scattering for spin lifetimes (τ s ) of supported/free-standing germanene, a prototypical strong SOC 2D Dirac material. We show that the effects of different substrates on τ s can surprisingly differ by two orders of magnitude. We find that substrate effects on τ s are closely related to substrate-induced modifications of the SOC-field anisotropy, which changes the spin flip scattering matrix elements. We propose a new electronic quantity, named spin-flip angle θ ↑↓ , to characterize spin relaxation through intervalley spin-flip scattering. We find that $τ^{-1}_s$ is approximately proportional to the averaged value of sin 2 (θ ↑↓ /2), which serves as a guiding parameter of controlling spin relaxation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electron–phonon physics from first principles using the EPW code

EPW is an open-source software for ab initio calculations of electron–phonon interactions and related materials properties. The code combines density functional perturbation theory and maximally localized Wannier functions to efficiently compute electron–phonon coupling matrix elements, and to perform predictive calculations of temperature-dependent properties and phonon-assisted quantum processes in bulk solids and low-dimensional materials. Here, we report on significant developments in the code since 2016, namely: a transport module for the calculation of charge carrier mobility under electric and magnetic fields using the Boltzmann transport equation; a superconductivity module for calculations of phonon-mediated superconductors using the anisotropic multi-band Eliashberg theory; an optics module for calculations of phonon-assisted indirect transitions; a module for the calculation of small and large polarons without supercells; and a module for calculating band structure renormalization and temperature-dependent optical spectra using the special displacement method. For each capability, we outline the methodology and implementation and provide example calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Primitive to conventional geometry projection for efficient phonon transport calculations

Abstract The primitive Wigner-Seitz cell and corresponding first Brillouin zone (FBZ) are typically used in calculations of lattice vibrational and transport properties as they contain the smallest number of degrees of freedom and thus have the cheapest computational cost. However, in complex materials, the FBZ can take on irregular shapes where lattice symmetries are not apparent. Thus, conventional cells (with more atoms and regular shapes) are often used to describe materials, though dynamical and transport calculations are more expensive. Here we discuss an efficient anharmonic lattice dynamic method that maps conventional cell dynamics to primitive cell dynamics based on translational symmetries. Such symmetries have not been utilized in typical lattice dynamical calculations. This leads to phase interference conditions that act like conserved quantum numbers and a conservation rule for phonon scattering that is hidden in conventional dynamics which significantly reduces the computational cost. We demonstrate this method for phonon transport in a variety of materials with inputs from first-principles calculations and attribute its efficiency to reduced scattering phase space and fewer summations in scattering matrix element calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Downfolding from ab initio to interacting model Hamiltonians: comprehensive analysis and benchmarking of the DFT+cRPA approach

Abstract Model Hamiltonians are regularly derived from first principles to describe correlated matter. However, the standard methods for this contain a number of largely unexplored approximations. For a strongly correlated impurity model system, here we carefully compare a standard downfolding technique with the best possible ground-truth estimates for charge-neutral excited-state energies and wave functions using state-of-the-art first-principles many-body wave function approaches. To this end, we use the vanadocene molecule and analyze all downfolding aspects, including the Hamiltonian form, target basis, double-counting correction, and Coulomb interaction screening models. We find that the choice of target-space basis functions emerges as a key factor for the quality of the downfolded results, while orbital-dependent double-counting corrections diminish the quality. Background screening of the Coulomb interaction matrix elements primarily affects crystal-field excitations. Our benchmark uncovers the relative importance of each downfolding step and offers insights into the potential accuracy of minimal downfolded model Hamiltonians.

Chemistry↗
DOE OSTI2025

A quantum eigenvalue solver based on tensor networks

Electronic ground states are of central importance in chemical simulations, but have remained beyond the reach of efficient classical algorithms except in cases of weak electron correlation or one-dimensional spatial geometry. We introduce a hybrid quantum-classical eigenvalue solver that constructs a wavefunction ansatz from a linear combination of matrix product states in rotated orbital bases, enabling the characterization of strongly correlated ground states with arbitrary spatial geometry. The energy is converged via a gradient-free generalized sweep algorithm based on quantum subspace diagonalization, with a potentially exponential speedup in the off-diagonal matrix element contractions upon translation into compact quantum circuits of linear depth in the number of qubits. Chemical accuracy is attained in numerical experiments for both a stretched water molecule and an octahedral arrangement of hydrogen atoms, achieving substantially better correlation energies compared to a unitary coupled-cluster benchmark, with orders of magnitude reductions in quantum resource estimates and a surprisingly high tolerance to shot noise. This proof-of-concept study suggests a promising new avenue for scaling up simulations of strongly correlated chemical systems on near-term quantum hardware.

chemistry↗
DOE OSTI2021

Unveiling the orbital texture of 1T-TiTe 2 using intrinsic linear dichroism in multidimensional photoemission spectroscopy

The momentum-dependent orbital character in crystalline solids, referred to as orbital texture, is of capital importance in the emergence of symmetry-broken collective phases, such as charge density waves as well as superconducting and topological states of matter. By performing extreme ultraviolet multidimensional angle-resolved photoemission spectroscopy for two different crystal orientations linked to each other by mirror symmetry, we isolate and identify the role of orbital texture in photoemission from the transition metal dichalcogenide 1T-TiTe 2 . By comparing our experimental results with theoretical calculations based on both a quantitative one-step model of photoemission and an intuitive tight-binding model, we unambiguously demonstrate the link between the momentum-dependent orbital orientation and the emergence of strong intrinsic linear dichroism in the photoelectron angular distributions. Our results represent an important step towards going beyond band structure (eigenvalues) mapping and learning about electronic wavefunction and orbital texture of solids by exploiting matrix element effects in photoemission spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On-line installation of the Superallowed Transition Beta-Neutrino Decay Ion Coincidence Trap

The Cabibbo-Kobayashi-Maskawa quark mixing matrix currently does not satisfy unitarity at the 2σ-level. This could be the result of an inaccurate value of one or both of its largest matrix elements V us and V ud . In the case of V ud , the most precise measurement is obtained from the f t -value measurements of superallowed beta-transitions between 0 + states. The accuracy of this determination can, in turn, be tested by extracting V ud in other transitions including superallowed transitions between mirror nuclei. The Superallowed Transition Beta-Neutrino Decay Ion Coincidence Trap (St. Benedict) is currently under construction at the Nuclear Science Laboratory of the University of Notre Dame to perform such a determination, with the goal of shedding more light on this tension with unitarity. St. Benedict will take a radioactive ion beam produced by TwinSol , thermalize it in a large volume gas catcher, then transport it in two separate differentially-pumped volumes using a radio-frequency (RF) carpet and a radio-frequency quadrupole (RFQ) ion guide before injecting it in an RFQ trap to create cool ion bunches for injection in the measurement Paul trap. In this paper, we detail the installation of the beam preparation components of St. Benedict, and present the results of the first RIBs successfully stopped and extracted from its gas catcher.

Brodeur, M.↗
DOE OSTI2021

Bond order predictions using deep neural networks

Machine learning is an extremely powerful tool for the modern theoretical chemist since it provides a method for bypassing costly algorithms for solving the Schrödinger equation. Already, it has proven able to infer molecular and atomic properties such as charges, enthalpies, dipoles, excited state energies, and others. Most of these machine learning algorithms proceed by inferring properties of individual atoms, even breaking down total molecular energy into individual atomic contributions. In this paper, we introduce a modified version of the Hierarchically Interacting Particle Neural Network (HIP-NN) capable of making predictions on the bonds between atoms rather than on the atoms themselves. We train the modified HIP-NN to infer bond orders for a large number of small organic molecules as computed via the Natural Bond Orbital package. We demonstrate that the trained model is extensible to molecules much larger than those in the training set by studying its performance on the COMP6 dataset. This method has applications in cheminformatics and force field parameterization and opens a promising future for machine learning models to predict other quantities that are defined between atoms such as density matrix elements, Hamiltonian parameters, and molecular reactivities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A coupled cluster framework for electrons and phonons

Here, we describe a coupled cluster framework for coupled systems of electrons and harmonic phonons. Neutral and charged excitations are accessed via the equation-of-motion version of the theory. Benchmarks on the Hubbard–Holstein model allow us to assess the strengths and weaknesses of different coupled cluster approximations, which generally perform well for weak to moderate coupling. Finally, we report progress toward an implementation for ab initio calculations on solids and present some preliminary results on finite-size models of diamond with a linear electron–phonon coupling. We also report the implementation of electron–phonon coupling matrix elements from crystalline Gaussian type orbitals within the PySCF program package.

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

Analytical derivatives of the individual state energies in ensemble density functional theory. II. Implementation on graphical processing units (GPUs)

Conical intersections control excited state reactivity, and thus, elucidating and predicting their geometric and energetic characteristics are crucial for understanding photochemistry. Locating these intersections requires accurate and efficient electronic structure methods. Unfortunately, the most accurate methods (e.g., multireference perturbation theories such as XMS-CASPT2) are computationally challenging for large molecules. The state-interaction state-averaged restricted ensemble referenced Kohn–Sham (SI-SA-REKS) method is a computationally efficient alternative. The application of SI-SA-REKS to photochemistry was previously hampered by a lack of analytical nuclear gradients and nonadiabatic coupling matrix elements. We have recently derived analytical energy derivatives for the SI-SA-REKS method and implemented the method effectively on graphical processing units. We demonstrate that our implementation gives the correct conical intersection topography and energetics for several examples. Furthermore, our implementation of SI-SA-REKS is computationally efficient, with observed sub-quadratic scaling as a function of molecular size. This demonstrates the promise of SI-SA-REKS for excited state dynamics of large molecular systems.

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

Construction of linearly independent non-orthogonal AGP states

We show how to construct a linearly independent set of antisymmetrized geminal power (AGP) states, which allows us to rewrite our recently introduced geminal replacement models as linear combinations of non-orthogonal AGPs. Furthermore, this greatly simplifies the evaluation of matrix elements and permits us to introduce an AGP-based selective configuration interaction method, which can reach arbitrary excitation levels relative to a reference AGP, balancing accuracy and cost as we see fit.

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

A non-hierarchical correlation discrete variable representation

The correlation discrete variable representation (CDVR) facilitates (multi-layer) multi-configurational time-dependent Hartree (MCTDH) calculations with general potentials. It employs a layered grid representation to efficiently evaluate all potential matrix elements appearing in the MCTDH equations of motion. The original CDVR approach and its multi-layer extension show a hierarchical structure: the size of the grids employed at the different layers increases when moving from an upper layer to a lower one. In this work, a non-hierarchical CDVR approach, which uses identically structured quadratures at all layers of the MCTDH wavefunction representation, is introduced. Here, the non-hierarchical CDVR approach crucially reduces the number of grid points required, compared to the hierarchical CDVR, shows superior scaling properties, and yields identical results for all three representations showing the same topology. Numerical tests studying the photodissociation of NOCl and the vibrational states of CH 3 demonstrate the accuracy of the non-hierarchical CDVR approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗