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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Cavity-modified molecular dipole switching dynamics

Polaritonic states, which are formed by resonances between a molecular excitation and the photonic mode of a cavity, have a number of useful properties that offer new routes to control molecular photochemistry using electric fields. To provide a theoretical description of how polaritonic states affect the real-time electron dynamics in molecules, a new method is described where the effects of strong light–molecule coupling are implemented using real-time electronic structure theory. The coupling between the molecular electronic states and the cavity is described by the Pauli–Fierz Hamiltonian, and transitions between polaritonic states are induced via an external time-dependent electric field using time-dependent configuration interaction (TDCI) theory, producing quantum electrodynamics TDCI (QED-TDCI). This method is used to study laser-induced ultrafast charge transfer and dipole-switching dynamics of the LiCN molecule inside a cavity. The increase in cavity coupling strength is found to have a significant impact on the energies and transition dipole moments of the molecule–cavity system. The convergence of the polaritonic state energies as a function of the number of included electronic and photonic basis states is discussed.

Chemistry↗
DOE OSTI2024

Alternative CNDOL Fockians for fast and accurate description of molecular exciton properties

CNDOL is an a priori, approximate Fockian for molecular wave functions. In this study, we employ several modes of singly excited configuration interaction (CIS) to model molecular excitation properties by using four combinations of the one electron operator terms. Those options are compared to the experimental and theoretical data for a carefully selected set of molecules. The resulting excitons are represented by CIS wave functions that encompass all valence electrons in the system for each excited state energy. The Coulomb–exchange term associated to the calculated excitation energies is rationalized to evaluate theoretical exciton binding energies. This property is shown to be useful for discriminating the charge donation ability of molecular and supermolecular systems. Multielectronic 3D maps of exciton formal charges are showcased, demonstrating the applicability of these approximate wave functions for modeling properties of large molecules and clusters at nanoscales. This modeling proves useful in designing molecular photovoltaic devices. Our methodology holds potential applications in systematic evaluations of such systems and the development of fundamental artificial intelligence databases for predicting related properties.

Chemistry↗
DOE OSTI2024

Tracking ultrafast non-adiabatic dissociation dynamics of the deuterated water dication molecule

Here, we applied reaction microscopy to elucidate fast non-adiabatic dissociation dynamics of deuterated water molecules after direct photo-double ionization at 61 eV with synchrotron radiation. For the very rare D + + O + + D breakup channel, the particle momenta, angular, and energy distributions of electrons and ions, measured in coincidence, reveal distinct electronic dication states and their dissociation pathways via spin–orbit coupling and charge transfer at crossings and seams on the potential energy surfaces. Notably, we could distinguish between direct and fast sequential dissociation scenarios. For the latter case, our measurements reveal the geometry and orientation of the deuterated water molecule with respect to the polarization vector that leads to this rare 3-body molecular breakup channel. Aided by multi-reference configuration-interaction calculations, the dissociation dynamics could be traced on the relevant potential energy surfaces and particularly their crossings and seams. This approach also unraveled the ultrafast time scales governing these processes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Equation-of-motion internally contracted multireference unitary coupled-cluster theory

The accurate computation of excited states remains a challenge in electronic structure theory, especially for systems with a ground state that requires a multireference treatment. In this work, we introduce a novel equation-of-motion (EOM) extension of the internally contracted multireference unitary coupled-cluster framework (ic-MRUCC), termed EOM-ic-MRUCC. EOM-ic-MRUCC follows the transform-then-diagonalize approach, in analogy to its non-unitary counterpart. By employing a projective approach to optimize the ground state, the method retains additive separability and proper scaling with system size. We show that excitation energies are size-intensive if the EOM operator satisfies the “killer” and the projective conditions. Furthermore, we propose to represent changes in the reference state upon electron excitation via projected many-body operators that span the active orbitals and show that the EOM equations formulated in this way are invariant with respect to active orbital rotations. We test the EOM-ic-MRUCC method truncated to single and double excitations by computing the potential energy curves for several excited states of a BeH2 model system, the HF molecule, and water undergoing symmetric dissociation. Across these systems, our method delivers accurate excitation energies and potential energy curves within 5 mE h (∼0.14 eV) from full configuration interaction. Here, we find that truncating the Baker–Campbell–Hausdorff series to fourfold commutators contributes negligible errors (on the order of 10 −5 E h or less), offering a practical route to highly accurate excited-state calculations with reduced computational overhead.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Ground and excited state gradients with end-to-end differentiable semiempirical quantum chemistry

Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.

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

Active space selection with self-healing diffusion Monte Carlo algorithms for periodic solids

Multideterminant Diffusion Monte Carlo (DMC) displays improved accuracy over single determinant DMC. Self-Healing Diffusion Monte Carlo (SHDMC) is a DMC based method that iteratively improves a multideterminant trial wavefunction. Although configuration interaction or complete active space (CAS) methods are very accurate and computationally feasible for many systems, they are not optimal for application to solids. SHDMC is accurate and designed for application to solids, so developing SHDMC based active space selection algorithms is a worthy endeavor. Here, we present and compare active space selection algorithms that are designed for use in conjunction with SHDMC, without relying on external approaches. For benchmarking, we calculated the ground state energy of a small unit cell of graphene and compared the results with a complete basis set extrapolated selected CI and a reference SHDMC trajectory. We found that systematically expanding the active space using an “auto-branching” algorithm optimally balances accuracy with computational practicality. To the best of our knowledge, this is the first work that demonstrates completely self-contained DMC-based active space selection algorithms that do not depend on external methods for determinant selection.

Spanedda, Nicole [ORNL]↗
DOE OSTI2019

Resonant inelastic X-ray scattering determination of the electronic structure of oxyhemoglobin and its model complex

Hemoglobin and myoglobin are oxygen-binding proteins with S = 0 heme {FeO 2 } 8 active sites. The electronic structure of these sites has been the subject of much debate. This study uses Fe K-edge X-ray absorption spectroscopy (XAS) and 1s2p resonant inelastic X-ray scattering (RIXS) to study oxyhemoglobin and a related heme {FeO 2 } 8 model compound, [(pfp)Fe(1-MeIm)(O 2 )] (pfp = meso-tetra(α,α,α,α- o -pivalamido-phenyl)porphyrin, or TpivPP, 1-MeIm = 1-methylimidazole) (pfpO 2 ), which was previously analyzed using L-edge XAS. The K-edge XAS and RIXS data of pfpO 2 and oxyhemoglobin are compared with the data for low-spin Fe II and Fe III [Fe(tpp)(Im) 2 ] 0/+ (tpp = tetra-phenyl porphyrin) compounds, which serve as heme references. The X-ray data demonstrate that pfpO 2 is similar to Fe II , while oxyhemoglobin is qualitatively similar to Fe III , but with significant quantitative differences. Density-functional theory (DFT) calculations show that the cjontrast between pfpO 2 and oxyhemoglobin is due to a distal histidine H bond to O 2 and the less hydrophobic environment in the protein, which lead to more backbonding into the O 2 . A valence bond configuration interaction multiplet model is used to analyze the RIXS data and show that pfpO 2 is dominantly Fe II with 6–8% Fe III character, while oxyhemoglobin has a very mixed wave function that has 50–77% Fe III character and a partially polarized Fe–O 2 π-bond.

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

Wave function analysis with a maximum flow algorithm

An efficient algorithm for computing the maximum-flow path in a network is applied to the identification of the dominant configuration state functions (CSFs) in a graphically contracted function (GCF), configuration interaction, wave function. The flow network is a space of spin-adapted CSFs represented by a Shavitt graph, wherein the nodes correspond to orbital occupations and spin quantum numbers. The graph nodes are connected by arcs, and an arc density is defined as sums of the associated squared CSF coefficients. A max-min approach determines an upper bound to the maximum possible incoming flow for each graph node. A backtracking step generates a candidate walk and is followed by a limited search of alternative branching paths for the dominant CSF. The arc density contributions are removed from the graph, and the algorithm is reapplied to the updated graph. This list of generated walks can be partitioned in order to guarantee that the dominant CSFs have been identified. All of the steps in this algorithm are computationally efficient and do not depend on the potentially large dimension of the underlying linear CSF expansion space. An analysis of low-lying valence states of C-2 illustrates the method.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Formation of silicon monoxide by radiative association: the impact of resonances

In this work, detailed quantum chemistry calculations within the multireference configuration interaction approximation with the Davidson correction are presented using an aug-cc-pV6Z basis set, for the potential energy curves and transition dipole moments between low lying molecular states of singlet spin symmetry for the SiO molecule. The high quality molecular data are used to obtain radiative association cross sections and rate coefficients for collisions between ground state Si and O atoms. Quantal calculations are compared with semiclassical results. Using a quantum kinetic theory of radiative association in which quasibound levels are assumed to be in local thermodynamic equilibrium, we find that resonances play an important role in enhancing the rate coefficients at low temperatures by several orders of magnitude from that predicted by standard quantum scattering formulations. Finally, these new molecular formation rates may have important implications for applications in astrophysics.

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

Parametric CI+MBPT calculations of Th I energies and g-factors for even states

In this study, we study energies and g-factors for a large number of levels of Th I. We found that the accuracy of configuration-interaction many-body perturbation theory can be substantially improved by using multiple adjustment parameters, which are introduced to regulate the second-order MBPT corrections for single-valence energies and Coulomb interaction. The results for energies and g-factors are in excellent agreement with experiment, with accuracy sufficient to positively assign theoretical states of given angular momentum and parity to experimental levels. Finally, this theory will be further developed in the future to treat other, even more complicated, actinides such as U I.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

CI-MBPT line strengths and atomic probabilities for some transitions of neutral iodine

Iodine is a promising alternative to expensive xenon as a propellant for electric propulsion devices. Iodine optical spectroscopy is a valuable diagnostic tool for optimization of the iodine propellant. Recently, some measurements have been conducted with cells containing iodine plasma, and iodine atomic absorption at several wavelengths (911, 906, 206 and 1315 nm) was measured with a Ti:sapphire laser. Motivated by these experiments and future applications, here we have calculated line strengths and atomic probabilities for the transitions between different types of 5s 2 5p 5 , 5s 2 5p 4 6s, 5s 2 5p 4 6p and 5s 2 5p 4 5d states of atomic iodine, using configuration-interaction many-body perturbation theory. Since experimental transition data for atomic iodine are only available for a relatively small number of transitions, we also performed calculations for the bromine atom, to further validate our approach. We found good agreement between our theory and available experimental data for both elements in most cases. Therefore, we believe that our predicted line strengths and transition probabilities are quite reliable and will be useful for iodine spectroscopy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Resonant x-ray absorption of strong-field-ionized CF 3 Br

We report on an experimental and theoretical study of strong-field laser ionization of CF 3 Br followed by resonant x-ray absorption at the Br K-edge. Distinct 1s -> 4p, 5p Rydberg transitions of Br q+ (q = 1-4) atomic ions are observed and identified with Hartree-Fock-Slater and relativistic configuration interaction calculations. Time-dependent density functional theory and ab initio molecular dynamics calculations were performed to simulate the dissociative ionization process and the molecular orbitals for the q = 1-4 charge states. Measurements were made with both parallel and perpendicular linear polarizations of the laser and x-rays, but dichroism was not observed, indicating negligible alignment by the laser ionization process. This result is explained by calculations on atomic Br and the molecular simulation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

A comprehensive study of the radiative properties of NO—a first step toward a complete air opacity

Here, in this paper we propose a methodology to calculate the radiative properties of the diatomic molecular constituents of air, and utilize the present approach to model the radiative properties of nitrogen monoxide, NO. We also investigate the important physics involved in calculating accurate radiative quantities for air, such as the Planck and Rosseland mean opacities, and emission and absorption coefficients, as well as the couplings accounted for in rovibrational calculations. Complete active space self-consistent field multi-reference configuration interaction (CAS-MRCI) calculations were performed in order to model the NO X 2 Π, a 4 Π, b 4 Σ - , 1 2 Σ + , 2 2 Σ + , 3 2 Σ + , G 2 Σ - , B' 2 Δ, (C, B) 2 2 Π, (H', L) 3 2 Π, and 1 2 Φ adiabatic states, and calculate the respective molecular data. The γ X 2 Π - A 2 Σ + , ε X 2 Π - D 2 Σ + , β' X 2 Π - B' 2 Δ, '11 000 Å' A 2 Σ + - D 2 Σ + , 'infrared' X 2 Π - X 2 Π and X 2 Π - (C, B) 2 2 Π (δ and β) band systems are investigated in monochromatic spectra calculations, as well as the Ogawa a 4 Π - b 4 Σ - band and several other band systems. Several conclusions are drawn, such as the importance of including the Ogawa band, which has not been included in previous air radiative models or comprehensive line-list calculations, as well as the importance of performing coupled rovibrational line-list calculations in order to accurately calculate the Rosseland means. We also found that the additional band systems modeled here contribute significantly to the total Planck and Rosseland means.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Exact sum rules with approximate ground states

Electromagnetic and weak transitions tell us a great deal about the structure of atomic nuclei. Yet modeling transitions can be difficult: it is often easier to compute the ground state, if only as an approximation, than excited states. One alternative is through transition sum rules, in particular the non-energy-weighted and energy-weighted sum rules, which can be computed as expectation values of operators. We investigate by computing sum rules for a variety of nuclei, comparing the numerically exact full configuration-interaction shell model, as a reference, to Hartree-Fock, projected Hartree-Fock, and the nucleon pair approximation. Furthermore, these approximations yield reasonable agreement, which we explain by prior work on the systematics of transition moments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

How effective is the Brink–Axel hypothesis for astrophysical weak rates?

We explore the effectiveness of the Brink–Axel hypothesis (BAH) for the computation of stellar electron capture (EC) and β-decay (BD) rates, namely that the transition strength function depends only upon the transition energy and not upon the details of the initial state. For this purpose, we calculated Gamow–Teller (GT) strength distributions for a selection of sd-shell nuclides, using two different microscopic models, namely the proton–neutron quasiparticle random phase approximation and the full configuration-interaction shell model, taking into account the first 100 states of both the initial and final nuclides. The GT transition strengths among these levels evolve with initial state energy. These transition strength functions we folded into weak-interaction mediated rates in stellar matter, specifically EC and BD rates, for a range of densities 10 g cm –3 ≤ ρ ≤ 10 11 g cm –3 and range of temperatures 1 GK ≤ T ≤ 30 GK. When transitions from excited states were approximated using the BAH, augmented by so-called 'back-resonance' transitions, the rates were affected by up to three orders of magnitude or more at high temperatures and densities. Furthermore the BAH is not a reliable approximation for the calculation of stellar rates, especially in high temperature–density environments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Nucleon–nucleon correlations in the extreme oxygen isotopes

There has been an upsurge of interest in two-nucleon decays thanks to the studies of nucleon–nucleon correlations. In our previous work (2021 Phys. Rev. Lett. 126 142501), based on a novel time-dependent three-body approach, we demonstrated that the energy and angular correlations of the emitted nucleons can shed light on the structure of nucleonic pairs formed inside the nucleus. In this work, we apply the new framework to study the decay dynamics and properties of some extreme proton-rich and neutron-rich oxygen isotopes, including two-proton (2p) decays of 11,12 O and two-neutron (2n) decay of 26 O. Here we show that the low- ℓ components of 11,12 O wave functions, which are affected by continuum and configuration-interaction effects, strongly impact decay dynamics and asymptotic correlations. In the calculated wave functions of 11,12 O, diproton and cigarlike structures merge together during the tunneling process and the resulting energy and angular correlations are very consistent with the experimental data. The asymptotic correlations of the 2n decay of 26 O dramatically change as the 2n decay energy approaches the zero-energy threshold. The small reported value of Q 2n suggests that the 2n decay of this nucleus can be understood in terms of the universal phase-space limit.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Benchmark of correlation matrix renormalization method in molecule calculations

In this paper, we report benchmark calculations of the correlation matrix renormalization (CMR) approach for 23 molecules in the well-established G2 molecule set. This subset represents molecules with spin-singlet ground state in a variety of chemical bonding and coordination environments. The QUAsi-atomic minimal basis-set orbitals (QUAMBOs) are used as local orbitals in both CMR and full configuration interaction (FCI) calculations for comparison. The results obtained from the calculations are also compared with available experimental data. It is shown that the CMR method produces binding and dissociation energy curves in good agreement with the QUAMBO-FCI calculations as well as experimental results. Finally, the CMR benchmark calculations yield a standard deviation of 0.09 Å for the equilibrium bond length and 0.018 Hartree/atom for the formation energy, with a gain of great computational efficiency which scales like Hartree–Fock method.

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

Accurate CI-MBPT calculation of radiative lifetimes and transition probabilities of neutral lanthanum (La I) odd states with J = 3/2

In this work, an accurate calculation of radiative lifetimes and transition probabilities of the first 21 odd states with a total angular momentum J = 3/2 was performed for La I using configuration-interaction many-body perturbation theory (CI-MBPT). A technique to obtain correct transition amplitudes was devised for pairs of states with substantial mixing that improved the theoretical accuracy. The results were compared to known experimental and theoretical data for the considered transitions. It has been shown that the transition data calculated with CI-MBPT and with corrections for pairs of mixed states are in good agreement with the latest experimental data of Den Hartog et al, 2015 J. Phys. B: At. Mol. Phys. 48, 155 001, better than other theoretical results. The transition data can find many applications in astrophysics, and the CI-MBPT formalism with mixing correction is applicable to other atoms similar to or more complex than La I.

36 MATERIALS SCIENCE↗