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

Strategies to Obtain Reliable Energy Landscapes from Embedded Multireference Correlated Wavefunction Methods for Surface Reactions

Embedded correlated wavefunction (ECW) theory is a powerful tool for studying ground- and excited-state reaction mechanisms and associated energetics in heterogeneous catalysis. Several factors are important to obtaining reliable ECW energies, critically the construction of consistent active spaces (ASs) along reaction pathways when using a multireference correlated wavefunction (CW) method that relies on a subset of orbital spaces in the configuration interaction expansion to account for static electron correlation, e.g., complete AS self-consistent field theory, in addition to the adequate partitioning of the system into a cluster and environment, as well as the choice of a suitable basis set and number of states included in excited-state simulations. Here, in this work, we conducted a series of systematic studies to develop best-practice guidelines for ground- and excited-state ECW theory simulations, utilizing the decomposition of NH 3 on Pd(111) as an example. We determine that ECW theory results are relatively insensitive to cluster size, the aug-cc-pVDZ basis set provides an adequate compromise between computational complexity and accuracy, and that a fixed-clean-surface approximation holds well for the derivation of the embedding potential. Additionally, we demonstrate that a merging approach, which involves generating ASs from the molecular fragments at each configuration, is preferable to a creeping approach, which utilizes ASs from adjacent structures as an initial guess, for the generation of consistent potential energy curves involving open-d-shell metal surfaces, and, finally, we show that it is essential to include bands of excited states in their entirety when simulating excited-state reaction pathways.

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

Cross-shell excited configurations in the structure of 34 Si

The cross-shell excited states of 34 Si have been investigated via β decays of the 4 − ground state and the 1 + isomeric state of 34 Al. Since the valence protons and valence neutrons occupy different major shells in the ground state as well as the intruder 1 + isomeric state of 34 Al, intruder levels of 34 Si are populated via allowed β decays. Spin assignments to such intruder levels of 34 Si were established through γ-γ angular correlation analysis for the negative-parity states. The configurations of such intruder states play crucial roles in our understanding of the N = 20 shell gap evolution. A configuration interaction model derived from the FSU Hamiltonian was utilized in order to interpret the intruder states in 34 Si. Furthermore, shell model interaction derived from a more fundamental theory with the valence space in medium similarity renormalization group method was also employed to interpret the structure of 34 Si.

20 ≤ A ≤ 38↗
DOE OSTI2016

Adsorbate-Induced Changes in Magnetic Interactions in Fe 2 (dobdc) with Adsorbed Hydrocarbon Molecules

Metal–organic frameworks, which are a special case of coordination polymers, form a class of materials with numerous applications due to their high porosities and large internal surface areas and in some cases also due to paramagnetic metal ions. A family of such materials of general formula M 2 (dobdc) (where M is a divalent metal ion, and dobdc 4 – is 2,5-dioxido-1,4-benzenedicarboxylate) has attracted considerable attention for gas separation, catalysis, and magnetism. Here, in this work, we explore the magnetic properties of a member of this class, Fe 2 (dobdc), both in the activated form (bare MOF after solvent removal) and when hydrocarbons are bound to the open coordination sites of the metal. We report quantum mechanical electronic structure calculations using both cluster models and a periodic model, and we compare our results to previously reported theoretical studies. We find that hydrocarbon adsorption only mildly affects the isotropic couplings, but that the isotropic magnetic couplings obtained with hybrid exchange–correlation functionals for the cluster models were found to be in good agreement with the experimental values, and the ones obtained using local exchange–correlation functionals with empirical Coulomb and exchange integrals for the periodic models were also found to agree well with experiments. Furthermore, hybrid density functionals for cluster models and local density functionals with empirical Coulomb and exchange integrals for periodic models are both found to give good agreement with the experimental result that the adsorption of ethylene changes the magnetic ordering. We also used second-order n-electron valence state perturbation theory and contracted spin-orbit configuration interaction to study the role of adsorption on the crystal-field splitting of the quintet manifold and on the single-ion anisotropy of the iron center.

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

Simulating the absorption spectra of helium clusters (N = 70, 150, 231, 300) using a charge transfer correction to superposition of fragment single excitations

Simulations of the n = 2 absorption spectra of HeN(N = 70, 150, 231, 300) clusters are reported, with nuclear configurations sampled by path integral molecular dynamics. The electronic structure is treated by a new approach, ALMO-CIS+CT, which is a formulation of configuration interaction singles (CIS) based on absolutely localized molecular orbitals (ALMOs). The method generalizes the previously reported ALMO-CIS model [K. D. Closser et al. J. Chem. Theory Comput. 11, 5791 (2015)] to include spatially localized charge transfer (CT) effects. It is designed to recover large numbers of excited states in atomic and molecular clusters, such as the entire n = 2 Rydberg band in helium clusters. ALMO-CIS+CT is shown to recover most of the error caused by neglecting charge transfer in ALMO-CIS and has comparable accuracy to standard CIS for helium clusters. For the n = 2 band, CT stabilizes states towards the blue edge by up to 0.5 eV. ALMO-CIS+CT retains the formal cubic scaling of ALMO-CIS with respect to system size. With improvements to the implementation over that originally reported for ALMO-CIS, ALMO-CIS+CT is able to treat helium clusters with hundreds of atoms using modest computing resources. Here, a detailed simulation of the absorption spectra associated with the 2s and 2p bands of helium clusters up to 300 atoms is reported, using path integral molecular dynamics with a spherical boundary condition to generate atomic configurations at 3 K. The main features of experimentally reported fluorescence excitation spectra for helium clusters are reproduced.

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

Understanding the nature of mean-field semiclassical light-matter dynamics: An investigation of energy transfer, electron-electron correlations, external driving, and long-time detailed balance

Semiclassical electrodynamics (with quantum matter plus classical electrodynamics fields) is an appealing approach for studying light-matter interactions, especially for realistic molecular systems. However, there is no unique semiclassical scheme. On the one hand, intermolecular interactions can be described instantaneously by static two-body interactions connecting two different molecules, while a classical transverse E field acts as a spectator at short distance; we will call this Hamiltonian no. I. On the other hand, intermolecular interactions can also be described as effects that are mediated exclusively through a classical one-body Ε field without any quantum effects at all (assuming we ignore electronic exchange); we will call this Hamiltonian no. II. Moreover, one can also mix these two different Hamiltonians into a third, hybrid Hamiltonian, which preserves quantum electron-electron correlations for lower excitations but describes higher excitations in a mean-field way. To investigate which semiclassical scheme is most reliable for practical use, here we study the real-time dynamics of a minimalistic many-site model—a pair of identical two-level systems (TLSs)—undergoing either resonance energy transfer (RET) or collectively driven dynamics. While both approaches (no. 1 and no. 2) perform reasonably well when there is no strong external excitation, we find that no single approach is perfect for all conditions (and all methods fail when a strong external field is applied). Each method has its own distinct problems: Hamiltonian no. I performs best for RET but behaves in a complicated manner for driven dynamics; Hamiltonian no. II is always stable, but obviously fails for RET at short distances. One key finding is that, for externally driven dynamics, a full configuration-interaction description of Hamiltonian no. I strongly overestimates the long-time electronic energy, highlighting the not obvious fact that, if one plans to merge quantum molecules with classical light, a full, exact treatment of electron-electron correlations can actually lead to worse results than a simple mean-field electronic structure treatment. Future work will need to investigate (i) how these algorithms behave in the context of more than a pair of TLSs and (ii) whether or not these algorithms can be improved in general by including crucial aspects of spontaneous emission.

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

Large Charge‐Transfer Energy in LiFePO 4 Revealed by Full‐Multiplet Calculation for the Fe L 3 ‐edge Soft X‐ray Emission Spectra

Abstract We analyzed the Fe 3d electronic structure in LiFePO 4 /FePO 4 (LFP/FP) nanowire with a high cyclability by using soft X‐ray emission spectroscopy (XES) combined with configuration‐interaction full‐multiplet (CIFM) calculation. The ex situ Fe L 2,3 ‐edge resonant XES (RXES) spectra for LFP and FP are ascribed to oxidation states of Fe 2+ and Fe 3+ , respectively. CIFM calculations for Fe 2+ and Fe 3+ states reproduced the Fe L 3 RXES spectra for LFP and FP, respectively. In the calculations for both states, the charge‐transfer energy was considerably larger than those for typical iron oxides, indicating very little electron transfer from the O 2p to Fe 3d orbitals and a weak hybridization on the Fe−O bond during the charge–discharge reactions.

Asakura, Daisuke↗
DOE OSTI2018

Double-buffered, heterogeneous CPU + GPU integral digestion algorithm for single-excitation calculations involving a large number of excited states

The most widely used quantum-chemical models for excited states are single-excitation theories, a category that includes configuration interaction with single substitutions, time-dependent density functional theory, and also a recently developed ab initio exciton model. When a large number of excited states are desired, these calculations incur a significant bottleneck in the “digestion” step in which two-electron integrals are contracted with density or density-like matrices. We present an implementation that moves this step onto graphical processing units (GPUs), and introduce a double-buffer scheme that minimizes latency by computing integrals on the central processing units (CPUs) concurrently with their digestion on the GPUs. An automatic code generation scheme simplifies the implementation of high-performance GPU kernels. For the exciton model, which requires separate excited-state calculations on each electronically coupled chromophore, the heterogeneous implementation described here results in speedups of 2–6× versus a CPU-only implementation. For traditional time-dependent density functional theory calculations, we obtain speedups of up to 5× when a large number of excited states is computed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Electronic structure of Li 1,2,3 +,0,– and nature of the bonding in Li 2,3 +,0,–

Abstract The current study of the small lithium molecules Li 2 +,0,− and Li 3 +,0,− focuses on the nature of the bonding in these molecules as well as their structures and energetics (bond energies, ionization energies, and electron affinities). Valence CASSCF (2s,2p) calculations incorporate nondynamical electron correlation in the calculations, while the corresponding multireference configuration interaction and coupled cluster calculations incorporate dynamical electron correlation. Treatment of nondynamical correlation is critical for properly describing the Li 2,3 +,0,− molecules as well as the Li − anion with dynamical correlation, in general, only fine‐tuning the predictions. All lithium molecules and ions are bound, with the Li 3 + and Li 2 + ions being the most strongly bound, followed by Li 3 − , Li 2 , Li 2 − and Li 3 . The minimum energy structures of Li 3 +,0,− are, respectively, an equilateral triangle, an isosceles triangle, and a linear structure. The results of SCGVB calculations are analyzed to obtain insights into the nature of the bonding in these molecules. An important finding of this work is that interstitial orbitals, a concept first put forward by McAdon and Goddard in 1985, play an essential role in the bonding of all lithium molecules considered here except for Li 2 . The interstitial orbitals found in the Li 3 +,0 molecules likely give rise to the non‐nuclear attractors/maxima observed in these molecules.

Chemistry↗
DOE OSTI2023

A promising intersection of excited‐state‐specific methods from quantum chemistry and quantum Monte Carlo

Abstract We present a discussion of recent progress in excited‐state‐specific quantum chemistry and quantum Monte Carlo alongside a demonstration of how a combination of methods from these two fields can offer reliably accurate excited state predictions across singly excited, doubly excited, and charge transfer states. Both of these fields have seen important advances supporting excited state simulation in recent years, including the introduction of more effective excited‐state‐specific optimization methods, improved handling of complicated wave function forms, and ways of explicitly balancing the quality of wave functions for ground and excited states. To emphasize the promise that exists at this intersection, we provide demonstrations using a combination of excited‐state‐specific complete active space self‐consistent field theory, selected configuration interaction, and state‐specific variance minimization. These demonstrations show that combining excited‐state‐specific quantum chemistry and variational Monte Carlo can be more reliably accurate than either equation of motion coupled cluster theory or multi‐reference perturbation theory, and that it can offer new clarity in cases where existing high‐level methods do not agree. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Software > Quantum Chemistry

Otis, Leon↗
DOE OSTI2022

Formation of phosphorus monoxide through the $$\mathbf {P}(^{4}S)+\mathbf {O}_{\mathbf {2}}(^{3}\Sigma ^{-})\rightarrow \mathbf {O}(^{3}P)+\mathbf {PO}(^{2}\Pi )$$ reaction

Phosphorus is a key and vital element for a diverse set of important biological molecules, being indispensable for life as we know. A deeper comprehension of its role in astrochemistry and atmospheric chemistry may aid in finding answers to how this element became available on Earth. The PO molecule is one of the main reservoirs of phosphorus in the interstellar medium (ISM), and a better understanding of the mechanisms and rate coefficients for its formation in the ISM is important for modelling its abundances. Here in this work, we perform multireference configuration interaction calculations on the formation of PO via the P( 4 S) + O 2 (3 Σ - ) reaction, analyzing its potential energy surface and rate coefficients for the global reaction on both doublet and quartet states. We also perform DFT (M06-2X) and CCSD(T) calculations, in order to compare the results. We found that the OPO system possesses a high multiconfigurational character, making DFT and CCSD methodologies not suitable for its potential energy landscape calculation. The rate coefficients have been calculated using the master equation system solver (MESS) package, and the results compared to recent experimental data. It is shown that the quartet state contributes for temperatures higher than 700K. The computed rate coefficient can be described by a modified Arrhenius equation [α(T/300) β exp (-γ/T)] with α = 1.44 x 10 -12 cm 3 s -1 , β= -1.66 and γ= 704 K.

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

Natural orbital description of the halo nucleus 6 He

Ab initio calculations of nuclei face the challenge of simultaneously describing the strong short-range internucleon correlations and the long-range properties of weakly bound halo nucleons. Natural orbitals, which diagonalize the one-body density matrix, provide a basis which is better matched to the physical structure of the many-body wave function. Here, we demonstrate that the use of natural orbitals significantly improves convergence for ab initio no-core configuration interaction calculations of the neutron halo nucleus 6 He, relative to the traditional oscillator basis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Electron impact excitation of N-like ions from the ICFT R-matrix calculation

We carried out intermediate coupling frame transformation (ICFT) R-matrix calculations on N-like ions (11≤Z≤20), which included a total of 272 fine structure levels in both the configuration interaction (CI) and close-coupling (CC) collision expansions. Comparisons with previous theoretical values were made in terms of energy levels, radiative transition rates, and oscillator strengths. The effective collision strengths (ϒ) were estimated in the temperature range of 10{sup 4}–10{sup 7} K. An excellent agreement with available computations was reached for dipole-allowed transitions, while differences in ϒ were found for non-dipole allowed and forbidden transitions. This discrepancy is mainly attributed to the difference in CC expansions. Our calculations provide a good supplement to the electron impact excitation data on N-like ions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Energy levels and radiative rates for transitions in Cr-like Kr XIII, Tc XX and Xe XXXI

Energy levels, radiative transition probabilities, Landé g-factors, and radiative lifetimes are reported for 144 levels of the 3p{sup 6}3d{sup 6} and 3p{sup 5}3d{sup 7} configurations of the Cr-like ions Kr XIII, Tc XX and Xe XXXI. Since to our knowledge no prior data for energy levels or A-values, experimental or theoretical, are available in the literature for these ions, calculations have been performed with three independent atomic structure codes, namely GRASP, FAC and a code implementing the QR approach, with increasing configuration interaction (CI) basis with each code. This is to establish the accuracy of the calculated results so that data can be confidently applied in the modeling of plasmas. Our recommended data are determined within the QR approach, because not only has the largest CI basis been included in this calculation, the reliability has already been established in our earlier work for other Cr-like ions, for which some experimental and theoretical data were available for comparison. Results for radiative rates (A-values) and transition wavelengths (λ) are listed for all strong E1, E2 and M1 transitions among the 144 levels of each ion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Level density of the sd-nuclei—Statistical shell-model predictions

Accurate knowledge of the nuclear level density is important both from a theoretical viewpoint as a powerful instrument for studying nuclear structure and for numerous applications. For example, astrophysical reactions responsible for the nucleosynthesis in the universe can be understood only if we know the nuclear level density. We use the configuration-interaction nuclear shell model to predict nuclear level density for all nuclei in the sd-shell, both total and for individual spins (only with positive parity). To avoid the diagonalization in large model spaces we use the moments method based on statistical properties of nuclear many-body systems. In the cases where the diagonalization is possible, the results of the moments method practically coincide with those from the shell-model calculations. Using the computed level densities, we fit the parameters of the Constant Temperature phenomenological model, which can be used by practitioners in their studies of nuclear reactions at excitation energies appropriate for the sd-shell nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

K- and L-edge X-ray absorption spectroscopy (XAS) and resonant inelastic X-ray scattering (RIXS) determination of differential orbital covalency (DOC) of transition metal sites

Continual advancements in the development of synchrotron radiation sources have resulted in X-ray based spectroscopic techniques capable of probing the electronic and structural properties of numerous systems. This review gives an overview of the application of metal K-edge and L-edge X-ray absorption spectroscopy (XAS), as well as Kα resonant inelastic X-ray scattering (RIXS), to the study of electronic structure in transition metal sites with emphasis on experimentally quantifying 3d orbital covalency. The specific sensitivities of K-edge XAS, L-edge XAS, and RIXS are discussed emphasizing the complementary nature of the methods. L-edge XAS and RIXS are sensitive to mixing between 3d orbitals and ligand valence orbitals, and to the differential orbital covalency (DOC), that is, the difference in the covalencies for different symmetry sets of the d orbitals. Both L-edge XAS and RIXS are highly sensitive to and enable separation of σ and π donor bonding and π back bonding contributions to bonding. Applying ligand field multiplet simulations, including charge transfer via valence bond configuration interactions, DOC can be obtained for direct comparison with density functional theory calculations and to understand chemical trends. Here, the application of RIXS as a probe of frontier molecular orbitals in a heme enzyme demonstrates the potential of this method for the study of metal sites in highly covalent coordination sites in bioinorganic chemistry.

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

Impact of electron density on the fixed-node errors in Quantum Monte Carlo of atomic systems

We analyze the effect of increasing charge density on the fixed-node errors in diffusion Monte Carlo (FNDMC). We compare FNDMC ground state energies of 4-electron atomic systems with nuclear charge using Hartree–Fock and restricted Configuration Interaction trial wave functions. We perform these calculations for several different values of the nuclear charge, and we analyze the systematic trends of the recovered correlation energies as a function of . For this testing system we demonstrate that the fixed-node error of the Hartree–Fock trial wave functions increases linearly with increasing the nuclear charge.

Chemistry↗
DOE OSTI2015

Relativistic opacities for astrophysical applications

Here, we report on the use of the Los Alamos suite of relativistic atomic physics codes to generate radiative opacities for the modeling of astrophysically relevant plasmas under local thermodynamic equilibrium (LTE) conditions. The atomic structure calculations are carried out in fine-structure detail, including full configuration interaction. Three example applications are considered: iron opacities at conditions relevant to the base of the solar convection zone, nickel opacities for the modeling of stellar envelopes, and samarium opacities for the modeling of light curves produced by neutron star mergers. In the first two examples, comparisons are made between opacities that are generated with the fully and semi-relativistic capabilities in the Los Alamos suite of codes. As expected for these highly charged, iron-peak ions, the two methods produce reasonably similar results, providing confidence that the numerical methods have been correctly implemented. However, discrepancies greater than 10% are observed for nickel and investigated in detail. In the final application, the relativistic capability is used in a preliminary investigation of the complicated absorption spectrum associated with cold lanthanide elements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Low-lying states and total internal partition sums of CH

The electronic structure and spin-orbit states of the CH radical have been systematically investigated using multi-reference configuration interaction (MRCI) and single-reference coupled-cluster (CC) methods. These calculations were performed in conjunction with large correlation-consistent basis sets of quadruple-, quintuple-, and sextuple-ζ quality. To achieve high accuracy, electronic energies for all states were extrapolated to the complete basis set (CBS) limit, enabling the detailed construction of potential energy curves and determination of reliable spectroscopic constants. Spin-orbit coupling effects were explicitly incorporated, and vibrational energy levels were computed via Numerov analysis. Furthermore, the resulting values exhibit good to excellent agreement with available experimental data. Dipole moment and transition dipole moment curves were evaluated to assess the opacity characteristics of CH, revealing that transitions such as Χ 2 Π (u′′ = 0) → Α 2 Δ (u′ = 0), Χ 2 Π (u′′ = 0) → Β 2 Σ − (u′ = 0), Χ 2 Π (u′′ = 0) → C 2 Σ + (u′ = 0), and Χ 2 Π (u′′ = 0) → D 2 Σ + (u′ = 3) are particularly probable. Finally, the total internal partition function sum (TIPS) of CH was computed over a broad temperature range (10–30,000 K) based on our high-accuracy ab initio results.

74 ATOMIC AND MOLECULAR PHYSICS↗