Search NASA⌕ Search

SEARCH · Search NASA

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

The voltage threshold for arcing for solar cells in LEO: Flight and ground test results

Ground and flight results of solar cell arcing in low Earth orbit (LEO) conditions are compared and interpreted. It is shown that an apparent voltage threshold for arcing may be produced by a strong power law dependence of arc rate on voltage, combined with a limited observation time. The change in this apparent threshold with plasma density is a reflection of the density dependence of the arc rate. A nearly linear dependence of arc rate on density is inferred from the data. A real voltage threshold for arcing for 2 by 2 cm solar cells may exist however, independent of plasma density, near -230 V relative to the plasma. Here, arc rates may change by more than an order of magnitude for a change of only 30 V in array potential. For 5.9 by 5.9 solar cells, the voltage dependence of the arc rate is steeper, and the data are insufficient to indicate the existence of an arcing increased by an atomic oxygen plasma, as is found in LEO, and by arcing from the backs of welded-through substrates.

Ferguson, D. C.↗
NASA NTRS1986

The voltage threshold for arcing for solar cells in Leo - Flight and ground test results

Ground and flight results of solar cell arcing in low earth orbit (LEO) conditions are compared and interpreted. It is shown that an apparent voltage threshold for arcing may be produced by a storage power law dependence of arc rate on voltage, combined with a limited observation time. The change in this apparent threshold with plasma density is a reflection of the density dependence of the arc rate. A nearly linear dependence of arc rate on density is inferred from the data. A real voltage threshold for arcing for 2 by 2 cm solar cells may exist however, independent of plasma density, near -230 V relative to the plasma. Here, arc rates may change by more than an order of magnitude for a change of only 30 V in array potential. For 5.9 by 5.9 solar cells, the voltage dependence of the arc rate is steeper, and the data are insufficient to indicate the existence of an arcing increased by an atomic oxygen plasma, as is found in LEO, and by arcing from the backs of welded-through substrates.

Ferguson, Dale C.↗
NASA NTRS1982

Coupling between the nickel-carbon and carbon-oxygen stretch motion in NiCO

Linear-combination-of-Gaussian-type-orbital (LCGTO) X-alpha calculations are performed on NiCO for various linear Ni-CO and NiC-O bond distances. The basis sets used are moderately large, approximately double zeta, and thus are fairly accurate within the X-alpha approximation. The electronic structure and equilibrium bond distances are consistent with recent CI calculations, indicating a much stronger Ni-CO bond than occurs in Ni(CO)4. The calculations agree to within 100/cm with the experimentally observed NiC-O stretch frequency. The Ni-CO stretch is predicted to occur at 656/cm, although at an intensity that is 640 times less. Compared with experimental results for CO adsorbed on single crystal Ni(100), the computed vibrational frequencies and intensities are less appropriate than the relevant experimental values for Ni(CO)4. This is interpreted to mean that an accurate description of the electronic structure of the nickel atom participating directly in the surface bond requires proper accounting of the bonds to its other nearest neighbors.

Dunlap, B. I.↗
DOE OSTI2023

Synthesis, crystal and electronic structures, linear and nonlinear optical properties, and photocurrent response of oxyhalides CeHaVIO 4 (Ha = Cl, Br; VI = Mo, W)

Four heteroanionic oxyhalides, CeClMoO 4 , CeBrMoO 4 , CeClWO 4 , and CeBrWO 4 , have been studied as multifunctional materials, which show a combination of good second harmonic generation (SHG) response and photocurrent signals. Millimeter-sized CeHaVIO 4 (Ha = Cl, Br; VI = Mo, W) crystals were grown by halide salt flux. The crystal structure of CeHaVIO 4 crystals was accurately determined by single-crystal X-ray diffraction. CeClMoO 4 , CeBrMoO 4 , and CeBrWO 4 are isostructural to each other, and crystallize in the acentric LaBrMoO 4 structure type. CeClWO 4 crystallizes in a new structure type with unit cell parameters of a = 19.6059(2) Å, b = 5.89450(10) Å, c = 7.80090(10) Å, and β = 101.4746(8)°. The bandgaps of CeHaVIO 4 fall into the range of 2.8(1)–3.1(1) eV, which are much smaller than those of isotypic LaHaVIO 4 (Ha = Cl, Br; VI = Mo, W) in the range of 3.9(1)–4.3(1) eV. The narrowing of bandgaps in CeHaVIO 4 originates from the presence of partially filled 4f orbitals of cerium atoms, which was confirmed by density functional theory (DFT) calculations. The moderate bandgaps make CeHaVIO 4 suitable for infrared nonlinear optical (IR NLO) applications. CeBrMoO 4 and CeBrWO 4 exhibit moderate SHG responses of 0.58× AGS and 0.46× AGS, respectively, and are both type-I phase-matching materials. Moderate SHG response, easy growth of crystals, high ambient stability, and type-I phase-matching behavior make CeBrMoO 4 and CeBrWO 4 great materials for IR NLO applications. CeHaVIO 4 films also exhibited good photocurrent response upon light radiation. Furthermore, this work demonstrates the rich structural chemistry of the REHaVIO 4 (RE = Y, La–Lu; Ha = Cl, Br; VI = Mo, W) family and the potential presence of more multifunctional materials.

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

Graph-based quantum response theory and shadow Born–Oppenheimer molecular dynamics

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born–Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. Further, the proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗