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↗
NASA NTRS1990

General contraction of Gaussian basis sets. II - Atomic natural orbitals and the calculation of atomic and molecular properties

A recently proposed scheme for using natural orbitals from atomic configuration interaction wave functions as a basis set for linear combination of atomic orbitals (LCAO) calculations is extended for the calculation of molecular properties. For one-electron properties like multipole moments, which are determined largely by the outermost regions of the molecular wave function, it is necessary to increase the flexibility of the basis in these regions. This is most easily done by uncontracting the outermost Gaussian primitives, and/or by adding diffuse primitives. A similar approach can be employed for the calculation of polarizabilities. Properties which are not dominated by the long-range part of the wave function, such as spectroscopic constants or electric field gradients at the nucleus, can generally be treated satisfactorily with the original atomic natural orbital sets.

Almlof, Jan↗
NASA NTRS1989

General contraction of Gaussian basis sets. Part 2: Atomic natural orbitals and the calculation of atomic and molecular properties

A recently proposed scheme for using natural orbitals from atomic configuration interaction (CI) wave functions as a basis set for linear combination of atomic orbitals (LCAO) calculations is extended for the calculation of molecular properties. For one-electron properties like multipole moments, which are determined largely by the outermost regions of the molecular wave function, it is necessary to increase the flexibility of the basis in these regions. This is most easily done by uncontracting the outmost Gaussian primitives, and/or by adding diffuse primitives. A similar approach can be employed for the calculation of polarizabilities. Properties which are not dominated by the long-range part of the wave function, such as spectroscopic constants or electric field gradients at the nucleus, can generally be treated satisfactorily with the original atomic natural orbital (ANO) sets.

Almloef, Jan↗
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↗
NASA NTRS1971

Study of different approximation in the calculation of g tensors - H2/+/

The theory of the g tensor in one-electron systems is briefly reviewed and calculations are performed in several ways for the hydrogen molecular ion in order to test approximations which must be made for larger systems. Approximate ground state wavefunctions are determined variationally. The first-order wavefunction with respect to the orbit-field perturbation is calculated, and the second-order g tensor is determined. The results of the various approximate calculations are compared and discussed. It is found that the linear combination of atomic orbitals method is rather poor and that two center integrals cannot be neglected.

De Montgolfier, P.↗
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↗
NASA NTRS1989

Charge exchange transition probability for collisions between unlike ions and atoms within the adiabatic approximation

A simple formula for the transition probability for electron exchange between unlike ions and atoms is established within the adiabatic approximation by employing the Linear Combination of Atomic Orbitals (LCAO) method. The formula also involves an adiabatic parameter, introduced by Massey, and thus the difficulties arising from the internal energy defect and the adiabatic approximation are avoided. Specific reactions Li(+++) + H to Li(++) + H(+) and Be(4+) + H to Be(3+) + H(+) are considered as examples. The calculated capture cross section results of the present work are compared with the experimental data and with the calculation of other authors over the velocity range of 10(7) cm/sec to 10(8) cm/sec.

Xu, Y. J.↗
NASA NTRS2001

Current Induced Forces on Adsorbates at Conducting Carbon Nanotubes

The problem addressed is force on an impurity due to electric current. In order to determine the size and direction of the force, as well as its dependence on the voltage, LCAO (Linear Combination of Atomic Orbitals) and the Green's Function approach were implemented. The resultant induced mean force was directed outwards in addition to a smaller component in the electron flow direction. At low bias, force is larger for larger radius, and crosses over the opposite trend when increasing the bias. The crossover current is approximately inversely proportional to the radius.

Mingo, N.↗
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↗