Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear combination of atomic orbitals”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Photoionization of the O 2 molecule

Ground state total and partial photoionization cross sections of the O 2 molecule have been produced with the configuration-average distorted-wave (CADW) approximation, and an energy-corrected pseudo independent-atom (PIA) method. The PIA approach is based on molecular orbitals expressed as a linear combination of atomic orbitals, and utilizes the corresponding atomic-orbital photoionization cross sections. Here we use an energy-corrected PIA approach with atomic photoionization cross sections computed using a CADW method. Our results are found to be in good agreement with photoionization measurements across a wide photon energy range, and in reasonably good agreement with the CADW calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

A new computational framework for spinor-based relativistic exact two-component calculations using contracted basis functions

Here, a new computational framework for spinor-based relativistic exact two-component (X2C) calculations is developed using contracted basis sets with a spin–orbit contraction scheme. Generally contracted, j-adapted basis sets of p-block elements using primitive functions in the correlation-consistent basis sets are constructed for the X2C Hamiltonian with atomic mean-field spin–orbit integrals (the X2CAMF scheme). The contraction coefficients are taken from atomic X2CAMF Hartree–Fock spinors, thereby following the simple concept of a linear combination of atomic orbitals. Benchmark calculations of spin–orbit splittings, equilibrium bond lengths, and harmonic vibrational frequencies demonstrate the accuracy and efficacy of the j-adapted spin–orbit contraction scheme.

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

Orbital design of flat bands in non-line-graph lattices via line-graph wave functions

Line-graph (LG) lattices are known for having flat bands (FBs) from the destructive interference of Bloch wave functions encoded in only lattice symmetry. In this work, we develop a generic atomic/molecular orbital design principle for FBs in non-LG lattices. Based on linear combination of atomic orbital theory, we demonstrate that the underlying wave-function symmetry of FBs in a LG lattice can be transformed into the atomic/molecular orbital symmetry in a non-LG lattice. We illustrate such orbital-designed topological FBs in three 2D non-LG, square, trigonal, and hexagonal lattices, where the designed orbitals faithfully reproduce the corresponding lattice symmetries of checkerboard, kagome, and diatomic-kagome lattices, respectively. Interestingly, systematic design of FBs with a high Chern number is also achieved based on the same principle. Fundamentally our theory enriches the FB physics; practically, it significantly expands the scope of FB materials, since most materials have multiple atomic/molecular orbitals at each lattice site, rather than a single s orbital mandated in graph theory and generic lattice models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

First Principle Investigation of Electronic, Transport, and Bulk Properties of Zinc-Blende Magnesium Sulfide

We have studied electronic, structural, and transport properties of zinc-blende magnesium sulfide (zb-MgS). We employed a local density approximation (LDA) potential and the linear combination of atomic orbitals (LCAO) method. Our computational method is able to reach the ground state of a material, as dictated by the second theorem of density functional theory (DFT). Consequently, our findings have the physical content of DFT and agree with available, corresponding experimental ones. The calculated band gap of zb-MgS, a direct gap equal to 4.43 eV, obtained at the experimental lattice constant of 5.620 Å, completely agrees with the experimental band gap of 4.45 ± 0.2 eV. We also report total (DOS) and partial (pDOS) densities of states, electron and hole effective masses, the equilibrium lattice constant, and the bulk modulus. The calculated pDOS also agree with the experiment for the description of the states at the top and the bottom of the valence and conduction bands, respectively.

42 ENGINEERING↗
DOE OSTI2023

Understanding the evolution of double perovskite band structure upon dimensional reduction

Recent investigations into the effects of dimensional reduction on halide double perovskites have revealed an intriguing change in band structure when the three-dimensional (3D) perovskite is reduced to a two-dimensional (2D) perovskite with inorganic sheets of monolayer thickness (n = 1). The indirect bandgap of 3D Cs 2 AgBiBr 6 becomes direct in the n = 1 perovskite whereas the direct bandgap of 3D Cs 2 AgTlBr 6 becomes indirect at the n = 1 limit. Here, we apply a linear combination of atomic orbitals approach to uncover the orbital basis for this bandgap symmetry transition with dimensional reduction. We adapt our previously established method for predicting band structures of 3D double perovskites for application to their 2D congeners, emphasizing new considerations required for the 2D lattice. In particular, we consider the inequivalence of the terminal and bridging halides and the consequences of applying translational symmetry only along two dimensions. The valence and conduction bands of the layered perovskites can be derived from symmetry adapted linear combinations of halide p orbitals propagated across the 2D lattice. The dispersion of each band is then determined by the bonding and antibonding interactions of the metal and halide orbitals, thus affording predictions of the essential features of the band structure. We demonstrate this analysis for 2D Ag–Bi and Ag–Tl perovskites with sheets of mono- and bilayer thickness, establishing a detailed understanding of their band structures, which enables us to identify the key factors that drive the bandgap symmetry transitions observed at the n = 1 limit. Importantly, these insights also allow us to make the general prediction that direct → indirect or indirect → direct bandgap transitions in the monolayer limit are most likely in double perovskite compositions that involve participation of metal d orbitals at the band edges or that have no metal-orbital contributions to the valence band, laying the groundwork for the targeted realization of this phenomenon.

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

Discontinuous Galerkin discretization for quantum simulation of chemistry

All-electron electronic structure methods based on the linear combination of atomic orbitals method with Gaussian basis set discretization offer a well established, compact representation that forms much of the foundation of modern correlated quantum chemistry calculations—on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growth of the number of integrals. Recent results have shown that, for some quantum and classical algorithms, moving to representations with diagonal two-body operators can result in dramatically lower asymptotic costs, even if the number of functions required increases significantly. We introduce a way to interpolate between the two regimes in a systematic and controllable manner, such that the number of functions is minimized while maintaining a block-diagonal structure of the two-body operator and desirable properties of an original, primitive basis. Techniques are analyzed for leveraging the structure of this new representation on quantum computers. Empirical results for hydrogen chains suggest a scaling improvement from O(N 4.5 ) in molecular orbital representations to O(N 2.6 ) in our representation for quantum evolution in a fault-tolerant setting, and exhibit a constant factor crossover at 15 to 20 atoms. Moreover, we test these methods using modern density matrix renormalization group methods classically, and achieve excellent accuracy with respect to the complete basis set limit with a speedup of 1–2 orders of magnitude with respect to using the primitive or Gaussian basis sets alone. These results suggest our representation provides significant cost reductions while maintaining accuracy relative to molecular orbital or strictly diagonal approaches for modest-sized systems in both classical and quantum computation for correlated systems.

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

Atomistic Studies of Individual Impurities and Impurity Complexes in III-V Semiconductors. Final report

The central hypothesis of our investigation is that reduced symmetries significantly influence the properties of individual impurities and impurity complexes in semiconductors. To disentangle the properties of the local environments of these impurities and impurity complexes from the behavior of an ensemble of impurities, at U. Iowa we performed theoretical calculations of specific local environments and developed the theoretical technology to calculate spin and charge transport through those local environments, in which the contact is a localized probe (scanning tunneling microscopy tip). These calculations were designed to be brought into comparison with experimental measurements made by Prof. Gupta’s group at Ohio State University. The principal environment that we focused on is the (110) surface of GaAs and other III-V semiconductors, because pristine surfaces can be revealed through a natural cleavage plane. We studied how the p-d hybridization of transition-metal dopants is influenced by the energetic location and the local environment’s symmetry-breaking of the d levels, and how the 4f-shell states of rare earth impurities couple to the surroundings. We calculated long-distance properties using tight-binding (linear combination of atomic orbitals, LCAO) Green’s function methods. The spin character of the dopants is often coherent on the time-scale of electronic transport through the region, and so correlations of the spin character of a dopant with that of a spinpolarized STM tip, or another dopant in the vicinity, can influence the measured conductivity from STM. We constructed a robust theory for these effects that can be used to describe quantitatively the measurements made by spin-polarized STM.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Ground State Properties of the Wide Band Gap Semiconductor Beryllium Sulfide (BeS)

We report the results from self-consistent calculations of electronic, transport, and bulk properties of beryllium sulfide (BeS) in the zinc-blende phase, and employed an ab-initio local density approximation (LDA) potential and the linear combination of atomic orbitals (LCAO). We obtained the ground state properties of zb-BeS with the Bagayoko, Zhao, and Williams (BZW) computational method, as enhanced by Ekuma and Franklin (BZW-EF). Our findings include the electronic energy bands, the total (DOS) and partial (pDOS) densities of states, electron and hole effective masses, the equilibrium lattice constant, and the bulk modulus. The calculated band structure clearly shows that zb-BeS has an indirect energy band gap of 5.436 eV, from Γ to a point between Γ and X, for an experimental lattice constant of 4.863 Å. This is in excellent agreement with the experiment, unlike the findings of more than 15 previous density functional theory (DFT) calculations that did not perform the generalized minimization of the energy functional, required by the second DFT theorem, which is inherent to the implementation of our BZW-EF method.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Sparsity of the electron repulsion integral tensor using different localized virtual orbital representations in local second-order Møller–Plesset theory

Utilizing localized orbitals, local correlation theory can reduce the unphysically high system-size scaling of post-Hartree–Fock (post-HF) methods to linear scaling in insulating molecules. The sparsity of the four-index electron repulsion integral (ERI) tensor is central to achieving this reduction. For second-order Møller–Plesset theory (MP2), one of the simplest post-HF methods, only the (ia|jb) ERIs are needed, coupling occupied orbitals i, j and virtuals a, b. In this paper, we compare the numerical sparsity (called the “ragged list”) and two other approaches revealing the low-rank sparsity of the ERI. The ragged list requires only one set of (localized) virtual orbitals, and we find that the orthogonal valence virtual-hard virtual set of virtuals originally proposed by Subotnik et al. gives the sparsest ERI tensor. To further compress the ERI tensor, the pair natural orbital (PNO) type representation uses different sets of virtual orbitals for different occupied orbital pairs, while the occupied-specific virtual (OSV) approach uses different virtuals for each occupied orbital. Here, our results indicate that while the low-rank PNO representation achieves significant rank reduction, it also requires more memory than the ragged list. The OSV approach requires similar memory to that of the ragged list, but it involves greater algorithmic complexity. An approximation (called the “fixed sparsity pattern”) for solving the local MP2 equations using the numerically sparse ERI tensor is proposed and tested to be sufficiently accurate and to have highly controllable error. A low-scaling local MP2 algorithm based on the ragged list and the fixed sparsity pattern is therefore promising.

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

Machine Learned Hückel Theory: Interfacing Physics and Deep Neural Networks

The Hückel Hamiltonian is an incredibly simple tight-binding model known for its ability to capture qualitative physics phenomena arising from electron interactions in molecules and materials. Part of its simplicity arises from using only two types of empirically fit physics-motivated parameters: the first describes the orbital energies on each atom and the second describes electronic interactions and bonding between atoms. By replacing these empirical parameters with machine-learned dynamic values, we vastly increase the accuracy of the extended Hückel model. The dynamic values are generated with a deep neural network, which is trained to reproduce orbital energies and densities derived from density functional theory. The resulting model retains interpretability, while the deep neural network parameterization is smooth and accurate and reproduces insightful features of the original empirical parameterization. Altogether, this work shows the promise of utilizing machine learning to formulate simple, accurate, and dynamically parameterized physics models.

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

Landscape of coexisting excitonic states in the insulating single-layer cuprates and nickelates

We present an ab initio study of the excitonic states of a prototypical high-temperature superconductor La 2 CuO 4 and compare them to the isostructural single-layer nickelate La 2 NiO 4 . A key difference in the low-energy electronic structure leads to very different excitonic behavior. Excitons in La 2 CuO 4 are delocalized and can freely move in the CuO 2 plane without disturbing the antiferromagnetic order. In contrast, in La 2 NiO 4 , we find the low-lying excitonic states to be extremely localized, producing a nearly flat dispersion. The theoretically obtained excitonic dispersion and behavior are in excellent agreement with RIXS observations. To classify the excitons, we project the electron-hole coupling onto each atomic site including the full manifold of atomic orbitals. We find the excitons to be composed of a linear combination of exciton classes, including Mott-Hubbard, d-d, and charge-transfer. The implication of these excitations to the high-T c pairing mechanism is also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Coupling of Nondegenerate Topological Modes in Nitrogen Core-Doped Graphene Nanoribbons

Nitrogen core-doping of graphene nanoribbons (GNRs) allows trigonal planar carbon atoms along the backbone of GNRs to be substituted by higher-valency nitrogen atoms. The excess valence electrons are injected into the π-orbital system of the GNR, thereby changing not only its electronic occupation but also its topological properties. We have observed this topological change by synthesizing dilute nitrogen core-doped armchair GNRs with a width of five atoms (N 2 -5-AGNRs). The incorporation of pairs of trigonal planar nitrogen atoms results in the emergence of topological boundary states at the interface between doped and undoped segments of the GNR. These topological boundary states are offset in energy by approximately ΔE = 300 meV relative to the topological end states at the termini of finite 5-AGNRs. Scanning tunneling microscopy (STM) and spectroscopy (STS) reveal that for finite GNRs the two types of topological states can interact through a linear combination of orbitals, resulting in a pair of asymmetric hybridized states. This behavior is captured by an effective Hamiltonian of nondegenerate diatomic molecules, where the analogous interatomic hybridization interaction strength is tuned by the distance between GNR topological modes.

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

Spin-orbit entanglement in the Color Glass Condensate

We compute the spin-orbit correlations of quarks and gluons at small-x and show that the helicity and the orbital angular momentum of individual partons are strongly anti-aligned even in unpolarized or spinless hadrons and nuclei. Combined with the fact that gluons in the Color Glass Condensate are linearly polarized, our finding indicates that the helicity and the orbital angular momentum of single gluons are maximally entangled in a quantum mechanical sense.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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

Optically Induced Static Magnetization in Metal Halide Perovskite for Spin‐Related Optoelectronics

Abstract Understanding the feasibility to couple semiconducting and magnetic properties in metal halide perovskites through interface design opens new opportunities for creating the next generation spin‐related optoelectronics. In this work, a fundamentally new phenomenon of optically induced magnetization achieved by coupling photoexcited orbital magnetic dipoles with magnetic spins at perovskite/ferromagnetic interface is discovered. The depth‐sensitive polarized neutron reflectometry combined with in situ photoexcitation setup, constitutes key evidence of this novel effect. It is demonstrated that a circularly polarized photoexcitation induces a stable magnetization signal within the depth up to 7.5 nm into the surface of high‐quality perovskite (MAPbBr 3 ) film underneath a ferromagnetic cobalt layer at room temperature. In contrast, a linearly polarized light does not induce any detectable magnetization in the MAPbBr 3 . The observation reveals that photoexcited orbital magnetic dipoles at the surface of perovskite are coupled with the spins of the ferromagnetic atoms at the interface, leading to an optically induced magnetization within the perovskite’s surface. The finding demonstrates that perovskite semiconductor can be bridged with magnetism through optically controllable method at room temperature in this heterojunction design. This provides the new concept of utilizing spin and orbital degrees of freedom in new‐generation spin‐related optoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The FMO-DFTB Method

Although the fragment molecular orbital (FMO) method enables electronic structure calculations with near-linear scaling behavior with respect to system size, the computational cost of ab initio methods typically employed in conjunction with FMO is still prohibitive for routine calculations of very large systems or long timescale molecular dynamics simulations. We, therefore, combined the FMO and density-functional tight-binding (DFTB) method, which is one of the emerging semi-empirical quantum chemical methods, and have demonstrated that FMO-DFTB is capable of performing geometry optimizations for systems containing up to one million atoms using limited computational resources. In this chapter, we will review the basics of the DFTB method first before introducing FMO-DFTB, focusing on the relationship with density functional theory and other FMO methodologies. We also demonstrate the latest scalings of FMO-DFTB2 and DFTB3 using three-dimensional water clusters, showing that the most favorable scaling is O(N1.16). Applications of FMO-DFTB to various systems are briefly summarized, and an outlook to future applications is provided.

Nishimoto, Yoshio↗
DOE OSTI2022

Complex Fermi–Löwdin orbital self-interaction correction

Here, this paper introduces the use of complex Fermi orbital descriptors (FODs) in the Fermi–Löwdin self-interaction-corrected density functional theory (FLOSIC). With complex FODs, the Fermi–Löwdin orbitals (FLOs) that are used to evaluate the SIC correction to the total energy become complex. Complex FLO-SIC (cFLOSIC) calculations based on the local spin density approximation produce total energies that are generally lower than the corresponding energies found with FLOSIC restricted to real orbitals (rFLOSIC). The cFLOSIC results are qualitatively similar to earlier Perdew–Zunger SIC (PZ-SIC) calculations using complex orbitals. The energy lowering stems from the exchange–correlation part of the self-interaction correction. The Hartree part of the correction is more negative in rFLOSIC. The energy difference between real and complex solutions is greater for more strongly hybridized FLOs in atoms and for FLOs corresponding to double and triple bonds in molecules. The case of N2 is examined in detail to show the differences between the real and complex FLOs. We show that the complex triple-bond orbitals are simple, and physically appealing combinations of π and σ g orbitals that have not been discussed before. Consideration of complex FODs, and resulting unitary transformations, underscores the fact that FLO centroids are not necessarily good guesses for FOD positions in a FLOSIC calculation.

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

Persistence of deuterium and tritium nuclear spin-polarization in presence of high-frequency plasma waves

We present first-principles numerical calculations of the depolarization rate of spin-polarized deuterium and tritium nuclei in realistic tokamak plasmas, driven by resonant interactions with plasma waves. Backed up by first-of-a-kind linear and nonlinear simulations, we find that alpha particle-driven Alfvénic modes cause only negligible depolarization, which is contrary to expectations in prior literature. Other Alfvénic instabilities can in principle degrade polarization, but only under conditions unlikely to be realized on transport timescales. By combining full-orbit particle tracing with a dedicated depolarization solver, we demonstrate that wave-driven depolarization is surprisingly weak in SPARC and ITER-scale devices. These results provide strong evidence that spin-polarized fuel can maintain its polarization long enough to boost fusion reactivity, opening a viable path toward substantially enhanced performance in magnetic confinement fusion power plants.

Cook, James William Setchfield [UK Atomic Energy A↗