Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CONFIGURATION INTERACTION”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Optimization of the Reax force field for the lithium–oxygen system using a high fidelity charge model

Studies using molecular dynamics (MD) have long struggled to simulate the failure modes of materials, predicting unrealistically high ductility and failing to capture brittle fracture. The primary cause of this shortcoming is an inadequate description of bond breaking. While reactive force fields such as ReaxFF show improvements compared to traditional force fields, the charge models used yield unphysical partial charges, especially during dissociation of ionic bonds. This flaw may be remedied by using the atom-condensed Kohn–Sham density functional theory (DFT) approximated to a second order (ACKS2) charge model for determining partial charges. In this work, we present a new ACKS2-enabled Reax force field for fracture simulations of lithium oxide systems, which was obtained by training against an extensive set of DFT, multireference configuration interaction (MRCI), and MRCI+Q reference data using genetic optimization techniques. This new force field significantly improves the bond breaking behavior, but still cannot fully capture the brittle fracture in MD simulations, suggesting more research is needed to improve simulation of brittle fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Electronic structure and correlated wave functions of a few electron quantum dots

The energy spectra and wave functions of a few electrons confined by a quasi-one-dimensional harmonic and anharmonic potentials have been studied by using a full configuration interaction method employing a Cartesian anisotropic Gaussian basis set. The energy spectra are classified into three regimes of the strength of confinement, namely, large, medium and small. The polyad quantum number defined by a total number of nodes in the wave functions is shown to be a key ingredient to interpret the energy spectra for the whole range of the confinement strength. The nodal pattern of the wave functions exhibits normal modes for the harmonic confining potential, indicating collective motions of electrons. These normal modes are shown to undergo a transition to local modes for an anharmonic potential with large anharmonicity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

A global ab initio potential energy surface for the X{sup 2}A{sup ′} ground state of the Si + OH → SiO + H reaction

We report the first global potential energy surface (PES) for the X{sup 2}A{sup ′} ground electronic state of the Si({sup 3}P) + OH(X{sup 2}Π) → SiO(X{sup 1}Σ{sub g}{sup +}) + H({sup 2}S) reaction. The PES is based on a large number of ab initio energies obtained from multireference configuration interaction calculations plus Davidson correction (MRCI+Q) using basis sets of quadruple zeta quality. Corrections were applied to the ab initio energies in the reactant channel allowing a proper description of long-range interactions between Si({sup 3}P) and OH(X{sup 2}Π). An analytical representation of the global PES has been developed by means of the reproducing kernel Hilbert space method. The reaction is found barrierless. Two minima, corresponding to the SiOH and HSiO isomers, and six saddle points, among which the isomerization transition state, have been characterized on the PES. The vibrational spectra of the SiOH/HSiO radicals have been computed from second-order perturbation theory and quantum dynamics methods. The structural, energetic, and spectroscopic properties of the two isomers are in good agreement with experimental data and previous high quality calculations.

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

Computing x-ray absorption spectra from linear-response particles atop optimized holes

State specific orbital optimized density functional theory (OO-DFT) methods, such as restricted open-shell Kohn–Sham (ROKS), can attain semiquantitative accuracy for predicting x-ray absorption spectra of closed-shell molecules. OO-DFT methods, however, require that each state be individually optimized. In this Communication, we present an approach to generate an approximate core-excited state density for use with the ROKS energy ansatz, which is capable of giving reasonable accuracy without requiring state-specific optimization. Herein, this is achieved by fully optimizing the core-hole through the core-ionized state, followed by the use of electron-addition configuration interaction singles to obtain the particle level. This hybrid approach can be viewed as a DFT generalization of the static-exchange (STEX) method and can attain ~0.6 eV rms error for the K-edges of C–F through the use of local functionals, such as PBE and OLYP. This ROKS(STEX) approach can also be used to identify important transitions for full OO ROKS treatment and can thus help reduce the computational cost of obtaining OO-DFT quality spectra. ROKS(STEX), therefore, appears to be a useful technique for the efficient prediction of x-ray absorption spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2015

The Structure, Anharmonic Vibrational Frequencies, and Intensities of NNHNN+

A semi-global potential energy surface (PES) and quartic force field (QFF) based on fitting high-level electronic structure energies are presented to describe the structures and spectroscopic properties of NNHNN+. The equilibrium structure of NNHNN+ is linear with the proton equidistant between the two nitrogen groups and thus of D(sub ∞h) symmetry. Vibrational second-order perturbation theory (VPT2) calculations based on the QFF fails to describe the proton “rattle” motion, i.e., the antisymmetric proton stretch, due to the very flat nature of PES around the global minimum, but performs properly for other modes with sharper potential wells. Vibrational self-consistent field/virtual state configuration interaction (VSCF/VCI) calculations using a version of MULTIMODE without angular momentum terms successfully describe this motion and predict the fundamental to be at 759 cm(exp -1). This is in good agreement with the value of 746 cm(exp -1) from a fixed-node diffusion Monte Carlo calculation and the experimental Ar-tagged result of 743 cm(exp -1). Other VSCF/VCI energies are in good agreement with other experimentally reported ones. Both double-harmonic intensity and rigorous MULTIMODE intensity calculations show the proton transfer fundamental has a very strong intensity.

Astrochemistry↗
DOE OSTI2013

Electron correlation dynamics in atoms and molecules

In this paper, we present quantum dynamical calculations on electron correlation dynamics in atoms and molecules using explicitly time-dependent ab initio configuration interaction theory. The goals are (i) to show that in which cases it is possible to switch off the electronic correlation by ultrashort laser pulses, and (ii) to understand the temporal evolution and the time scale on which it reappears. We characterize the appearance of correlation through electron-electron scattering when starting from an uncorrelated state, and we identify pathways for the preparation of a Hartree-Fock state from the correlated, true ground state. Exemplary results for noble gases, alkaline earth elements, and selected molecules are provided. For Mg we show that the uncorrelated state can be prepared using a shaped ultrashort laser pulse.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Evaluating the gradients of localized diabatic state energies and couplings at minimum cost

We calculate the gradients of Boys diabatic state energies and couplings when the electronic vector space is generated by configuration interaction singles. Our approach follows the Lagrangian approach of Paz and Glover (rather than direct differentiation of the adiabatic-to-diabatic approaches that have been published previously). The result is that we achieve a dramatic increase in savings over previous approaches, and the present approach should be immediately useful to scientists focused on electronic relaxation, especially chemists studying electron transfer who wish to go beyond the Condon approximation. Here, a future extension to time-dependent density functional theory in the Tamm–Dancoff approximation is clear.

Chemical physics↗
DOE OSTI2013

Full dimensional potential energy surface for the ground state of H{sub 4}{sup +} system based on triatomic-in-molecules formalism

In this work, we present a global potential energy surface for the ground electronic state of the H{sub 4}{sup +} based on ab initio calculations. The final fit is based on triatomics-in-molecules (TRIM) approximation and it includes extra four-body terms for the better description of some discrepancies found on the TRIM model. The TRIM method itself allows a very accurate description of the asymptotic regions. The global fit uses more than 19 000 multireference configuration interaction ab initio points. The global potential energy surface has an overall root mean square error of 0.013 eV for energies up to 2 eV above the global minimum. This work presents an analysis of the stationary points, reactant and product channels, and crossing between the two lowest TRIM adiabatic states. It is as well included a brief description of the two first excited states of the TRIM matrix, concluding that TRIM method is a very good approximation not only for the ground state but also for at least two of the excited states of H{sub 4}{sup +} system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Ground state property calculations of LiH n complexes using IBM Qiskit’s quantum simulator

In this study, the variational quantum eigensolver (VQE) on a quantum simulator is used in calculating ground state electronic structure properties of the LiH n , n = 1–3, complexes including their singly charged ions. Results calculated using classical electronic structure algorithms are also included. We investigate the use of the unitary coupled cluster with singles and doubles (UCCSD) Ansatz using VQE within Qiskit and compare results to full configuration interaction (FCI) calculations. Computed ground state energies, electron affinities, ionization potentials, and dipole moments are considered. We report the first-of-its-kind simulated quantum computing results of selected LiH n species and use the parity orbital to qubit mapping scheme. We find that VQE/UCCSD results are comparable to classical coupled clusters with singles and doubles for all considered systems with respect to FCI. A VQE calculation cost evaluation is included in which we evaluate performance using both Jordan–Wigner and parity orbital to qubit mapping schemes. We also discuss some of the current limitations of utilizing VQE for the study of chemical systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

A diabatic representation of the two lowest electronic states of Li{sub 3}

Using the Multi-Reference Configuration Interaction method, the adiabatic potential energy surfaces of Li{sub 3} are computed. The two lowest electronic states are bound and exhibit a conical intersection. By fitting the calculated potential energy surfaces to the cubic E ⊗ ε Jahn-Teller model we extract the effective Jahn-Teller parameters corresponding to Li{sub 3}. These are used to set up the transformation matrix which transforms from the adiabatic to a diabatic representation. This diabatization method gives a Hamiltonian for Li{sub 3} which is free from singular non-adiabatic couplings and should be accurate for large internuclear distances, and it thereby allows for bound dynamics in the vicinity of the conical intersection to be explored.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Conjugate Fermi holes and its manifestation in He-like systems

The structure of genuine and conjugate Fermi holes in two-electron atomic systems, namely He and He-like atomic ions, has been studied relying on accurate full configuration interaction wave functions. The standard Fermi hole exists in the vicinity of region in the two-electron coordinate space satisfying the well-known condition, r{sub 1} = r{sub 2}, while the conjugate Fermi hole exists in the vicinity of region close to this genuine Fermi hole but satisfying r{sub 1} ≠ r{sub 2} instead of r{sub 1} = r{sub 2}. Existence of these holes has shown to give an insightful interpretation of the origin of the first Hund rule and of the anomalously strong angular correlation manifested in the series of the singlet-triplet pair of singly-excited states of the aforementioned systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

ATOMIC DATA OF Cu I FOR THE INVESTIGATION OF ELEMENT ABUNDANCE

A complete set of atomic data of Cu I, including the energy levels, oscillator strengths, and photoionization cross sections, is theoretically studied to investigate element abundance including nonlocal thermodynamic equilibrium (NLTE) effects. The calculations are carried out by using the R-matrix method in the LS-coupling scheme. Twenty terms of Cu II are utilized as target states, and extensive configuration interactions are included to properly delineate the quantum states of Cu II and Cu I. One hundred thirteen bound states and 1699 oscillator strengths for E1 transitions between these states are obtained. Photoionization cross sections for all bound states are calculated in a photon energy range covering 1.28 Ry from the threshold of the respective state. Resonances shown in the photoionization cross sections are identified, and some strong resonances are expected to play an important role in NLTE modeling. The atomic data in this work represent the first complete data set for copper abundance studies. Our results are compared with the experimental and other theoretical data wherever available.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Dramatic changes in electronic structure revealed by fractionally charged nuclei

Discontinuous changes in the electronic structure upon infinitesimal changes to the Hamiltonian are demonstrated. These are revealed in one and two electron molecular systems by full configuration interaction (FCI) calculations when the realm of the nuclear charge is extended to be fractional. FCI electron densities in these systems show dramatic changes in real space and illustrate the transfer, hopping, and removal of electrons. This is due to the particle nature of electrons seen in stretched systems and is a manifestation of an energy derivative discontinuity at constant number of electrons. Dramatic errors of density functional theory densities are seen in real space as this physics is missing from currently used approximations. The movements of electrons in these simple systems encapsulate those in real physical processes, from chemical reactions to electron transport and pose a great challenge for the development of new electronic structure methods.

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

Excitonic optical properties of wurtzite ZnS quantum dots under pressure

By means of atomistic empirical pseudopotentials combined with a configuration interaction approach, we have studied the optical properties of wurtzite ZnS quantum dots in the presence of strong quantum confinement effects as a function of pressure. We find the pressure coefficients of quantum dots to be highly size-dependent and reduced by as much as 23% in comparison to the bulk value of 63 meV/GPa obtained from density functional theory calculations. The many-body excitonic effects on the quantum dot pressure coefficients are found to be marginal. The absolute gap deformation potential of quantum dots originates mainly from the energy change of the lowest unoccupied molecular orbital state. Finally, we find that the exciton spin-splitting increases nearly linearly as a function of applied pressure.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Probing ab initio emergence of nuclear rotation

Structural phenomena in nuclei, from shell structure and clustering to superfluidity and collective rotations and vibrations, reflect emergent degrees of freedom. Ab initio theory describes nuclei directly from a fully microscopic formulation. We can therefore look to ab initio theory as a means of exploring the emergence of effective degrees of freedom in nuclei. For the illustrative case of emergent rotational bands in the $Be$ isotopes, we establish an understanding of the underlying oscillator space and angular momentum (orbital and spin) structure. We consider no-core configuration interaction (NCCI) calculations for 7, 9, 11 Be with the Daejeon16 internucleon interaction. Although shell model or rotational degrees of freedom are not assumed in the ab initio theory, the NCCI results are suggestive of the emergence of effective shell model degrees of freedom (0$\hbar$ω and 2$\hbar$ω excitations) and LS-scheme rotational degrees of freedom, consistent with an Elliott–Wilsdon SU(3) description. Finally, these results provide some basic insight into the connection between emergent effective collective rotational and shell model degrees of freedom in these light nuclei and the underlying ab initio microscopic description.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Effective bond orders from two-step spin–orbit coupling approaches: The I{sub 2}, At{sub 2}, IO{sup +}, and AtO{sup +} case studies

The nature of chemical bonds in heavy main-group diatomics is discussed from the viewpoint of effective bond orders, which are computed from spin–orbit wave functions resulting from spin–orbit configuration interaction calculations. The reliability of the relativistic correlated wave functions obtained in such two-step spin–orbit coupling frameworks is assessed by benchmark studies of the spectroscopic constants with respect to either experimental data, or state-of-the-art fully relativistic correlated calculations. The I{sub 2}, At{sub 2}, IO{sup +}, and AtO{sup +} species are considered, and differences and similarities between the astatine and iodine elements are highlighted. In particular, we demonstrate that spin–orbit coupling weakens the covalent character of the bond in At{sub 2} even more than electron correlation, making the consideration of spin–orbit coupling compulsory for discussing chemical bonding in heavy (6p) main group element systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Study of the electronic and rovibronic structure of the X {sup 2}Σ{sup +}, A {sup 2}Π, and B {sup 2}Σ{sup +} states of AlO

The electronic structure of the X {sup 2}Σ{sup +}, A {sup 2}Π, and B {sup 2}Σ{sup +} states of aluminum monoxide (AlO) are studied via ab initio multi-reference configuration interaction calculations. Core correlation corrections, several basis sets, and active space choices are considered. Angular momentum and spin-orbit coupling terms are obtained at different levels of theory. The resulting ab initio curves are used to solve the associated rovibronic problem for the total angular momentum J up to 112.5 and then also refined by fitting to the experimental wavenumbers available in the literature, reproducing them with the root-mean-square error of 0.07 cm{sup −1}. Theoretical rovibronic energy levels of AlO in its X {sup 2}Σ{sup +}, A {sup 2}Π, and B {sup 2}Σ{sup +} electronic states are presented including those from the X − B blue-green system.

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

Calculation of spin–orbit couplings using RASCI spinless one-particle density matrices: Theory and applications

Here we present the formalism and implementation for calculations of spin–orbit couplings (SOCs) using the Breit–Pauli Hamiltonian and non-relativistic wave functions described by the restricted active space configuration interaction (RASCI) method with general excitation operators of spin-conserving spin-flipping, ionizing, and electron-attaching types. The implementation is based on the application of the Wigner–Eckart theorem within the spin space, which enables the calculation of the entire SOC matrix based on the explicit calculation of just one transition between the two spin multiplets. Numeric results for a diverse set of atoms and molecules highlight the importance of a balanced treatment of correlation and adequate basis sets and illustrate the overall robust performance of RASCI SOCs. The new implementation is a useful addition to the methodological toolkit for studying spin-forbidden processes and molecular magnetism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗