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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
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 OSTI2019

Dirac fermions and flat bands in the ideal kagome metal FeSn

A kagome lattice of 3d transition metal ions is a versatile platform for correlated topological phases hosting symmetry-protected electronic excitations and magnetic ground states. However, the paradigmatic states of the idealized two-dimensional kagome lattice-Dirac fermions and flat bands-have not been simultaneously observed. Here, we use angle-resolved photoemission spectroscopy and de Haas-van Alphen quantum oscillations to reveal coexisting surface and bulk Dirac fermions as well as flat bands in the antiferromagnetic kagome metal FeSn, which has spatially decoupled kagome planes. Our band structure calculations and matrix element simulations demonstrate that the bulk Dirac bands arise from in-plane localized Fe-3d orbitals, and evidence that the coexisting Dirac surface state realizes a rare example of fully spin-polarized two-dimensional Dirac fermions due to spin-layer locking in FeSn. Finally, the prospect to harness these prototypical excitations in a kagome lattice is a frontier of great promise at the confluence of topology, magnetism and strongly correlated physics.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Singlet fission: the role of molecular packing [Book Chapter]

The crucial role of molecular packing in singlet fission in dimers, aggregates, and solids is summarized and the key roles of Davydov splitting and the size of the electronic matrix element are pointed out. Here, simple rules and computer algorithms for these effects are described and their shortcomings are listed.

Buchanan, Eric A.↗
DOE OSTI2019

Optimal arrangements of 1,3-diphenylisobenzofuran molecule pairs for fast singlet fission

A simplified version of the frontier orbital model has been applied to pairs of C 2 , C 2v , C s , and C 1 symmetry 1,3-diphenylisobenzofuran rotamers to determine their best packing for fast singlet fission (SF). For each rotamer the square of the electronic matrix element for SF was calculated at 2.2 × 10 9 pair geometries and a few thousand most significant physically accessible local maxima were identified in the six-dimensional space of mutual arrangements. At these pair geometries, SF energy balance was evaluated, relative SF rate constants were approximated using Marcus theory, and the SF rate constant k SF was maximized by further optimization of the geometry of the molecular pair. The process resulted in 142, 67, 214, and 291 unique geometries for the C 2 , C 2v , C s , and C 1 symmetry molecular pairs, respectively, predicted to be superior to the C 2 symmetrized known crystal pair structure. Furthermore, these optimized pair geometries and their triplet biexciton binding energies are reported as targets for crystal engineering and/or covalent dimer synthesis, and as possible starting points for high-level pair geometry optimizations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2019

Efficient construction of generalized master equation memory kernels for multi-state systems from nonadiabatic quantum-classical dynamics

Techniques derived from the generalized quantum master equation (GQME) framework have provided the basis for elucidating energy and charge transfer in systems ranging from molecular solids to photosynthetic complexes. Recently, the nonperturbative combination of the GQME with quantum-classical methods has resulted in approaches whose accuracy and efficiency exceed those of the original quantumclassical schemes while offering significant accuracy improvements over perturbative expansions of the GQME. In this work, we show that, while the non-Markovian memory kernel required to propagate the GQME scales quartically with the number of subsystem states, the number of trajectories required scales at most quadratically when using quantum-classical methods to construct the kernel. We then present an algorithm that allows further acceleration of the quantum-classical GQME by offering a way to selectively sample the kernel matrix elements that are most important to the process of interest. We demonstrate the utility of these advances by applying the combination of Ehrenfest mean field theory with the GQME (MF-GQME) to models of the Fenna-Matthews-Olson (FMO) complex and the light harvesting complex II (LHCII), with 7 and 14 states, respectively. This allows us to show that the MF-GQME is able to accurately capture all the relevant dynamical time scales in LHCII: the initial nonequilibrium population transfer on the femtosecond time scale, the steady state-type trapping on the picosecond time scale, and the long time population relaxation. Remarkably, all of these physical effects spanning tens of picoseconds can be encoded in a memory kernel that decays only after ~65 fs.

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

Influence of the polarization anisotropy on the linewidth enhancement factor and reflection sensitivity of 1.55- μ m InP-based InAs quantum dash lasers

This work investigates the effect of the polarization anisotropy on the linewidth enhancement factor and the reflection sensitivity of InAs/InP quantum dash semiconductor lasers. The results show that the small linewidth enhancement factor and high stability against external optical feedback are obtained for nanostructures oriented perpendicular to the cavity axis as opposed to those oriented parallel to the cavity axis. Effective simulations on the critical feedback level of these two lasers are also in agreement with experimental results. Such anisotropy is attributed to the polarization dependence of the transition matrix element in these quantum nanostructures.

Dong, B. (ORCID:0000000158266723)↗
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↗
DOE OSTI2023

Electronic heat generation in semiconductors: Non-equilibrium excitation and evolution of zone-edge phonons via electron–phonon scattering in photo-excited germanium

We investigate experimentally and using first-principles theory the generation of phonons and the relaxation of carriers on picosecond timescales across the Brillouin zone of photo-excited Ge by inter-valley electron–phonon scattering. The phonons generated are typical of those generated in semiconductor devices, contributing to the accumulation of heat within the material. We simulate the time-evolution of phonon populations, based on first-principles band structure and electron–phonon and phonon–phonon matrix elements, and compare them to data from time-resolved x-ray diffuse scattering experiments, performed at the Linac Coherent Light Source x-ray free-electron laser facility, following photo-excitation by a 50 fs near-infrared optical pulse. We show that the intensity of the non-thermal x-ray diffuse scattering signal, which is observed to grow substantially near the L-point of the Brillouin zone over 3–5 ps, is due to phonons generated by scattering of carriers between the Δ and L valleys. These phonons have low group velocities, resulting in a heat bottleneck. With the inclusion of phonon decay through 3-phonon processes, the simulations also account for other non-thermal features observed in the x-ray diffuse scattering intensity, which are due to anharmonic phonon–phonon scattering of the phonons initially generated by electron–phonon scattering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Spin–orbit effects on the electronic and optical properties of lead iodide

Lead iodide (PbI2) has gained much interest due to its direct electronic gap in the visible range and layered crystal structure. It has thereby been considered as a promising material for applications in atomically thin optoelectronic devices. In this work, we present a detailed investigation of the effect of spin–orbit coupling (SOC) that arises from the presence of heavy atoms on the electronic and optical properties of PbI2 using first-principles calculations based on density-functional theory and many-body perturbation theory. We find that SOC not only alters the bandgap but also induces the mixing of orbital characters, resulting in a significant change in the overall band structure and charge carrier effective masses. Moreover, the band orbital mixing caused by SOC results in the dramatic change in optical transition matrix elements and, correspondingly, the absorption spectrum. Our experimentally measured absorption spectra validate the calculation results and demonstrate the importance of SOC in the optical processes of PbI2. Our findings provide insights that are important for the potential use of PbI2 as a material platform for visible optoelectronic devices.

Physics↗
DOE OSTI2023

Vibrational heat-bath configuration interaction with semistochastic perturbation theory using harmonic oscillator or VSCF modals

Vibrational heat-bath configuration interaction (VHCI)—a selected configuration interaction technique for vibrational structure theory—has recently been developed in two independent works [J. H. Fetherolf and T. C. Berkelbach, J. Chem. Phys. 154, 074104 (2021); A. U. Bhatty and K. R. Brorsen, Mol. Phys. 119, e1936250 (2021)], where it was shown to provide accuracy on par with the most accurate vibrational structure methods with a low computational cost. Here, we eliminate the memory bottleneck of the second-order perturbation theory correction using the same (semi)stochastic approach developed previously for electronic structure theory. This allows us to treat, in an unbiased manner, much larger perturbative spaces, which are necessary for high accuracy in large systems. Stochastic errors are easily controlled to be less than 1 cm−1. We also report two other developments: (i) we propose a new heat-bath criterion and an associated exact implicit sorting algorithm for potential energy surfaces expressible as a sum of products of one-dimensional potentials; (ii) we formulate VHCI to use a vibrational self-consistent field (VSCF) reference, as opposed to the harmonic oscillator reference configuration used in previous reports. Our tests are done with quartic and sextic force fields, for which we find that with VSCF, the minor improvements to accuracy are outweighed by the higher computational cost associated the matrix element evaluations. We expect VSCF-based VHCI to be important for more general potential representations, for which the harmonic oscillator basis function integrals are no longer analytic.

Chemistry↗
DOE OSTI2023

Tailoring light-induced charge transfer and intersystem crossing in FeCO using time-dependent spin–orbit configuration interaction

Real-time (RT) electronic structure methods provide a natural framework for describing light–matter interactions in arbitrary time-dependent electromagnetic fields (EMF). Optically induced excited state transitions are of particular interest, which require tuned EMF to drive population transfer to and from the specific state(s) of interest. Intersystem crossing, or spin-flip, may be driven through shaped EMF or laser pulses. These transitions can result in long-lived “spin-trapped” excited states, which are especially useful for materials requiring charge separation or protracted excited state lifetimes. Time-dependent configuration interaction (TDCI) is unique among RT methods in that it may be implemented in a basis of eigenstates, allowing for rapid propagation of the time-dependent Schrödinger equation. The recent spin–orbit TDCI (TD-SOCI) enables a real-time description of spin-flip dynamics in an arbitrary EMF and, therefore, provides an ideal framework for rational pulse design. The present study explores the mechanism of multiple spin-flip pathways for a model transition metal complex, FeCO, using shaped pulses designed to drive controlled intersystem crossing and charge transfer. These results show that extremely tunable excited state dynamics can be achieved by considering the dipole transition matrix elements between the states of interest.

Chemistry↗
DOE OSTI2024

Electron transfer in a crystalline cytochrome with four hemes

Diffusion of electrons over distances on the order of 100 μm has been observed in crystals of a small tetraheme cytochrome (STC) from Shewanella oneidensis [J. Huang et al . J. Am. Chem. Soc. 142, 10459–10467 (2020)]. Electron transfer between hemes in adjacent subunits of the crystal is slower and more strongly dependent on temperature than had been expected based on semiclassical electron-transfer theory. Here, in this work, we explore explanations for these findings by molecular-dynamics simulations of crystalline and monomeric STC. New procedures are developed for including time-dependent quantum mechanical energy differences in the gap between the energies of the reactant and product states and for evaluating fluctuations of the electronic-interaction matrix element that couples the two hemes. Rate constants for electron transfer are calculated from the time- and temperature-dependent energy gaps, coupling factors, and Franck–Condon-weighted densities of states using an expression with no freely adjustable parameters. Back reactions are considered, as are the effects of various protonation states of the carboxyl groups on the heme side chains. Interactions with water are found to dominate the fluctuations of the energy gap between the reactant and product states. The calculated rate constant for electron transfer from heme IV to heme Ib in a neighboring subunit at 300 K agrees well with the measured value. However, the calculated activation energy of the reaction in the crystal is considerably smaller than observed. We suggest two possible explanations for this discrepancy. The calculated rate constant for transfer from heme I to II within the same subunit of the crystal is about one-third that for monomeric STC in solution.

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

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗