Search NASA⌕ Search

SEARCH · Search NASA

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 289 records · Page 16

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Toward scalable many-body calculations for nuclear open quantum systems using the Gamow Shell Model

Drip-line nuclei have very different properties from those of the valley of stability, as they are weakly bound and resonant. Therefore, the models devised for stable nuclei can no longer be applied therein. Hence, a new theoretical tool, the Gamow Shell Model (GSM), has been developed to study the many-body states occurring at the limits of the nuclear chart. GSM is a configuration interaction model based on the use of the so-called Berggren basis, which contains bound, resonant and scattering states, so that inter-nucleon correlations are fully taken into account and the asymptotes of extended many-body wave functions are precisely handled. However, large complex symmetric matrices must be diagonalized in this framework, therefore the use of very powerful parallel machines is needed therein. So, in order to fully take advantage of their power, a 2D partitioning scheme using hybrid MPI/OpenMP parallelization has been developed in our GSM code. The specificities of the 2D partitioning scheme in the GSM framework will be described and illustrated with numerical examples. It will then be shown that the introduction of this scheme in the GSM code greatly enhances its capabilities.

2D partitioning↗

The electron affinity of the uranium atom

The results of a combined experimental and computational study of the uranium atom are presented with the aim of determining its electron affinity. Experimentally, the electron affinity of uranium was measured via negative ion photoelectron spectroscopy of the uranium atomic anion, U – . Computationally, the electron affinities of both thorium and uranium were calculated by conducting relativistic coupled-cluster and multi-reference configuration interaction calculations. The experimentally determined value of the electron affinity of the uranium atom was determined to be 0.309 ± 0.025 eV. The computationally predicted electron affinity of uranium based on composite coupled cluster calculations and full four-component spin–orbit coupling was found to be 0.232 eV. Predominately due to a better convergence of the coupled cluster sequence for Th and Th – , the final calculated electron affinity of Th, 0.565 eV, was in much better agreement with the accurate experimental value of 0.608 eV. In both scenarios, the ground state of the anion corresponds to electron attachment to the 6d orbital.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Transition matrices and orbitals from reduced density matrix theory

In this contribution, we report two different methodologies for characterizing the electronic structure reorganization occurring when a chromophore undergoes an electronic transition. For the first method, we start by setting the theoretical background necessary to the reinterpretation through simple tensor analysis of (i) the transition density matrix and (ii) the natural transition orbitals in the scope of reduced density matrix theory. This novel interpretation is made more clear thanks to a short compendium of the one-particle reduced density matrix theory in a Fock space. The formalism is further applied to two different classes of excited states calculation methods, both requiring a single-determinant reference, that express an excited state as a hole-particle mono-excited configurations expansion, to which particle-hole correlation is coupled (time-dependent Hartree-Fock/time-dependent density functional theory) or not (configuration interaction single/Tamm-Dancoff approximation). For the second methodology presented in this paper, we introduce a novel and complementary concept related to electronic transitions with the canonical transition density matrix and the canonical transition orbitals. Their expression actually reflects the electronic cloud polarisation in the orbital space with a decomposition based on the actual contribution of one-particle excitations from occupied canonical orbitals to virtual ones. This approach validates our novel interpretation of the transition density matrix elements in terms of the Euclidean norm of elementary transition vectors in a linear tensor space. A proper use of these new concepts leads to the conclusion that despite the different principles underlying their construction, they provide two equivalent excited states topological analyses. This connexion is evidenced through simple illustrations of (in)organic dyes electronic transitions analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Towards a spectroscopically accurate set of potentials for heavy hydride laser cooling candidates: Effective core potential calculations of BaH

BaH (and its isotopomers) is an attractive molecular candidate for laser cooling to ultracold temperatures and a potential precursor for the production of ultracold gases of hydrogen and deuterium. The theoretical challenge is to simulate the laser cooling cycle as reliably as possible and this paper addresses the generation of a highly accurate ab initio {sup 2}Σ{sup +} potential for such studies. The performance of various basis sets within the multi-reference configuration-interaction (MRCI) approximation with the Davidson correction is tested and taken to the Complete Basis Set (CBS) limit. It is shown that the calculated molecular constants using a 46 electron effective core-potential and even-tempered augmented polarized core-valence basis sets (aug-pCVnZ-PP, n = 4 and 5) but only including three active electrons in the MRCI calculation are in excellent agreement with the available experimental values. The predicted dissociation energy D{sub e} for the X{sup 2}Σ{sup +} state (extrapolated to the CBS limit) is 16 895.12 cm{sup −1} (2.094 eV), which agrees within 0.1% of a revised experimental value of <16 910.6 cm{sup −1}, while the calculated r{sub e} is within 0.03 pm of the experimental result.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic structure theory with molecular point group symmetries on quantum annealers

Quantum computation has the potential to revolutionize quantum chemistry through major speedups in computation times and an exponential reduction in computational resources. Here, we combine the symmetry-adapted Jordan–Wigner encoding based on the full Boolean symmetry group $\mathbb{Z}$$^{k}_{2}$ with our new implementation of the Xia–Bian–Kais (XBK) method for improving the efficiency of electronic structure theory calculations on quantum annealers, particularly by reducing the number of qubits needed to achieve the same accuracy. By providing a more extensive symmetry-adapted encoding (SAE) than previous work, we are able to simulate molecules larger than those previously reported that have been studied using methods developed for quantum annealers and without using an active space. We calculated the potential energy surfaces of H 2 , LiH, He 2 , H 2 O, O 2 , N 2 , Li 2 , F 2 , CO, BH 3 , NH 3 , and CH 4 , with the largest molecule in the STO-6G basis set requiring 16 qubits with our SAE, and compared them with full configuration interaction results. The application of SAE to the XBK method provides an exponential reduction in the size of the Hilbert space and scales well with the size of the problem. It does not introduce significant additional errors for even or large values of a key variational parameter that determines the number of ancilla qubits used in the XBK method’s Hamiltonian embedding, or for certain molecules such as He 2 and H 2 O. Here, we provide an explanation for this behavior and a recommendation on the usage of our method. In addition, we briefly discuss the potential of extracting electronic excited states from our method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ab initio calculations of the electronic structure of the low-lying states for the ultracold LiYb molecule

Ab initio techniques have been applied to investigate the electronic structure of the LiYb molecule. The potential energy curves have been computed in the Born–Oppenheimer approximation for the ground and 29 low-lying doublet and quartet excited electronic states. Complete active space self-consistent field, multi-reference configuration interaction, and Rayleigh Schrödinger perturbation theory to second order calculations have been utilized to investigate these states. The spectroscopic constants, ω{sub e}, R{sub e}, B{sub e}, …, and the static dipole moment, μ, have been investigated by using the two different techniques of calculation with five different types of basis. The eigenvalues, E{sub v}, the rotational constant, B{sub v}, the centrifugal distortion constant, D{sub v}, and the abscissas of the turning points, R{sub min} and R{sub max}, have been calculated by using the canonical functions approach. The comparison between the values of the present work, calculated by different techniques, and those available in the literature for several electronic states shows a very good agreement. Twenty-one new electronic states have been studied here for the first time.

74 ATOMIC AND MOLECULAR PHYSICS↗

Seniority number description of potential energy surfaces: Symmetric dissociation of water, N{sub 2}, C{sub 2}, and Be{sub 2}

The present study further explores the concept of the seniority number (Ω) by examining different configuration interaction (CI) truncation strategies in generating compact wave functions in a systematic way. While the role of Ω in addressing static (strong) correlation problem has been addressed in numerous previous studies, the usefulness of seniority number in describing weak (dynamic) correlation has not been investigated in a systematic way. Thus, the overall objective in the present work is to investigate the role of Ω in addressing also dynamic electron correlation in addition to the static correlation. Two systematic CI truncation strategies are compared beyond minimal basis sets and full valence active spaces. One approach is based on the seniority number (defined as the total number of singly occupied orbitals in a determinant) and another is based on an excitation-level limitation. In addition, molecular orbitals are energy-optimized using multiconfigurational-self-consistent-field procedure for all these wave functions. The test cases include the symmetric dissociation of water (6-31G), N{sub 2} (6-31G), C{sub 2} (6-31G), and Be{sub 2} (cc-pVTZ). We find that the potential energy profile for H{sub 2}O dissociation can be reasonably well described using only the Ω = 0 sector of the CI wave function. For the Be{sub 2} case, we show that the full CI potential energy curve (cc-pVTZ) is almost exactly reproduced using either Ω-based (including configurations having up to Ω = 2 in the virtual-orbital-space) or excitation-based (up to single-plus-double-substitutions) selection methods, both out of a full-valence-reference function. Finally, in dissociation cases of N{sub 2} and C{sub 2}, we shall also consider novel hybrid wave functions obtained by a union of a set of CI configurations representing the full valence space and a set of CI configurations where seniority-number restriction is imposed for a complete set (full-valence-space and virtual) of correlated molecular orbitals, simultaneously. We discuss the usefulness of the seniority number concept in addressing both static and dynamic electron correlation problems along dissociation paths.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Communication: CDFT-CI couplings can be unreliable when there is fractional charge transfer

Constrained density functional theory with configuration interaction (CDFT-CI) is a useful, low-cost tool for the computational prediction of electronic couplings between pseudo-diabatic constrained electronic states. Such couplings are of paramount importance in electron transfer theory and transition state theory, among other areas of chemistry. Unfortunately, CDFT-CI occasionally fails significantly, predicting a coupling that does not decay exponentially with distance and/or overestimating the expected coupling by an order of magnitude or more. In this communication, we show that the eigenvalues of the difference density matrix between the two constrained states can be used as an a priori metric to determine when CDFT-CI are likely to be reliable: when the eigenvalues are near 0 or ±1, transfer of a whole electron is occurring, and CDFT-CI can be trusted. We demonstrate the utility of this metric with several illustrative examples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-photon absorption cross sections within equation-of-motion coupled-cluster formalism using resolution-of-the-identity and Cholesky decomposition representations: Theory, implementation, and benchmarks

The equation-of-motion coupled-cluster (EOM-CC) methods provide a robust description of electronically excited states and their properties. Here, we present a formalism for two-photon absorption (2PA) cross sections for the equation-of-motion for excitation energies CC with single and double substitutions (EOM-CC for electronically excited states with single and double substitutions) wave functions. Rather than the response theory formulation, we employ the expectation-value approach which is commonly used within EOM-CC, configuration interaction, and algebraic diagrammatic construction frameworks. In addition to canonical implementation, we also exploit resolution-of-the-identity (RI) and Cholesky decomposition (CD) for the electron-repulsion integrals to reduce memory requirements and to increase parallel efficiency. The new methods are benchmarked against the CCSD and CC3 response theories for several small molecules. We found that the expectation-value 2PA cross sections are within 5% from the quadratic response CCSD values. The RI and CD approximations lead to small errors relative to the canonical implementation (less than 4%) while affording computational savings. RI/CD successfully address the well-known issue of large basis set requirements for 2PA cross sections calculations. The capabilities of the new code are illustrated by calculations of the 2PA cross sections for model chromophores of the photoactive yellow and green fluorescent proteins.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

CAS without SCF—Why to use CASCI and where to get the orbitals

The complete active space self-consistent field (CASSCF) method has seen broad adoption due to its ability to describe the electronic structure of both the ground and excited states of molecules over a broader swath of the potential energy surface than is possible with the simpler Hartree–Fock approximation. However, it also has a reputation for being unwieldy, computationally costly, and un-black-box. Here, we discuss a class of alternatives, complete active space configuration interaction (CASCI) methods, paying particular attention to their application to electronic excited states. The goal of this Perspective is fourfold. First, we argue that CASCI is not merely an approximation to CASSCF, in that it can be designed to have important qualitative advantages over CASSCF. Second, we present several insights drawn from our experience experimenting with different schemes for computing orbitals to be employed in CASCI. Third, we argue that CASCI is well suited for application to nanomaterials. Lastly, we reason that, with the rise in new low-scaling approaches for describing multireference systems, there is a greater need than ever to develop new methods for defining orbitals that provide an efficient and accurate description of both static correlation and electronic excitations in a limited active space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multireference linearized coupled cluster theory for strongly correlated systems using matrix product states

We propose a multireference linearized coupled cluster theory using matrix product states (MPSs-LCC) which provides remarkably accurate ground-state energies, at a computational cost that has the same scaling as multireference configuration interaction singles and doubles, for a wide variety of electronic Hamiltonians. These range from first-row dimers at equilibrium and stretched geometries to highly multireference systems such as the chromium dimer and lattice models such as periodic two-dimensional 1-band and 3-band Hubbard models. The MPS-LCC theory shows a speed up of several orders of magnitude over the usual Density Matrix Renormalization Group (DMRG) algorithm while delivering energies in excellent agreement with converged DMRG calculations. Also, in all the benchmark calculations presented here, MPS-LCC outperformed the commonly used multi-reference quantum chemistry methods in some cases giving energies in excess of an order of magnitude more accurate. As a size-extensive method that can treat large active spaces, MPS-LCC opens up the use of multireference quantum chemical techniques in strongly correlated ab initio Hamiltonians, including two- and three-dimensional solids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗