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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

No Core CI calculations for light nuclei with chiral 2- and 3-body forces

The atomic nucleus is a self-bound system of strongly interacting nucleons. In No-Core Configuration Interaction calculations, the nuclear wavefunction is expanded in Slater determinants of single-nucleon wavefunctions (Configurations), and the many-body Schrödinger equation becomes a large sparse matrix problem. The challenge is to reach numerical convergence to within quantified numerical uncertainties for physical observables using finite truncations of the infinite-dimensional basis space. We discuss strategies for constructing and solving the resulting large sparse matrices for a set of low-lying eigenvalues and eigenvectors on current multicore computer architectures. Several of these strategies have been implemented in the code MFDn, a hybrid MPI/OpenMP Fortran code for ab initio nuclear structure calculations that scales well to over 200,000 cores. We discuss how the similarity renormalization group can be used to improve the numerical convergence. We present results for excitation energies and other selected observables for 8 Be and 12 C using realistic 2- and 3-body forces obtained from chiral perturbation theory. Finally, we demonstrate that collective phenomena such as rotational band structures can emerge from these microscopic calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Molecular electric moments calculated by using natural orbital functional theory

The molecular electric dipole, quadrupole, and octupole moments of a selected set of 21 spin-compensated molecules are determined employing the extended version of the Piris natural orbital functional 6 (PNOF6), using the triple-ζ Gaussian basis set with polarization functions developed by Sadlej, at the experimental geometries. The performance of the PNOF6 is established by carrying out a statistical analysis of the mean absolute errors with respect to the experiment. The calculated PNOF6 electric moments agree satisfactorily with the corresponding experimental data and are in good agreement with the values obtained by accurate ab initio methods, namely, the coupled-cluster single and doubles and multi-reference single and double excitation configuration interaction methods.

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

Role of momentum in the generator-coordinate method applied to barrier penetration

Nuclear fission at barrier-top energies is conventionally modeled by a one-dimensional Schrödinger equation applied to internal fission channels, but that treatment is hard to justify in the configuration-interaction approach to nuclear Hamiltonians. Here we show that inclusion of states of finite momentum by the generator coordinate method (GCM) considerably extends the range of energies at which GCM-based Hamiltonians could reproduce the Schrödinger treatment. Furthermore, the transmission probabilities for crossing the barrier are calculated by a discrete version of Kohn's variational method, which may also be useful for other systems of interacting fermions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Understanding the chemical bonding of ground and excited states of HfO and HfB with correlated wavefunction theory and density functional approximations

Knowledge of the chemical bonding of HfO and HfB ground and low-lying electronic states provides essential insights into a range of catalysts and materials that contain Hf–O or Hf–B moieties. Here, we carry out high-level multi-reference configuration interaction theory and coupled cluster quantum chemical calculations on these systems. We compute full potential energy curves, excitation energies, ionization energies, electronic configurations, and spectroscopic parameters with large quadruple-ζ and quintuple-ζ quality correlation consistent basis sets. We also investigate equilibrium chemical bonding patterns and effects of correlating core electrons on property predictions. Differences in the ground state electron configuration of HfB(X 4 Σ - ) and HfO(X 1 Σ + ) lead to a significantly stronger bond in HfO than HfB, as judged by both dissociation energies and equilibrium bond distances. We extend our analysis to the chemical bonding patterns of the isovalent HfX (X = O, S, Se, Te, and Po) series and observe similar trends. We also note a linear trend between the decreasing value of the dissociation energy (D e ) from HfO to HfPo and the singlet–triplet energy gap (ΔE S–T ) of the molecule. Finally, we compare these benchmark results to those obtained using density functional theory (DFT) with 23 exchange–correlation functionals spanning multiple rungs of “Jacob’s ladder.” When comparing DFT errors to coupled cluster reference values on dissociation energies, excitation energies, and ionization energies of HfB and HfO, we observe semi-local generalized gradient approximations to significantly outperform more complex and high-cost functionals.

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

Energies and E1, M1, E2, M2 transition rates for states of the 2s{sup 2}2p, 2s2p{sup 2}, and 2p{sup 3} configurations in boron-like ions between N III and Zn XXVI

Energies, E1, M1, E2, M2 transition rates, line strengths, oscillator strengths, and lifetimes from relativistic configuration interaction calculations are reported for the states of the (1s{sup 2})2s{sup 2}2p, 2s2p{sup 2}, and 2p{sup 3} configurations in all boron-like ions between N III and Zn XXVI. Valence, core-valence, and core-core correlation effects were accounted for through single-double multireference (SD-MR) expansions to increasing sets of active orbitals.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Energy levels and radiative transition rates for Ba XLVIII

Highlights: • Lowest 431 levels of Ba XLVIII are calculated. • Radiative data for E1, E2, M1 and M2 transitions from the lowest 3 levels to 431 levels are presented. • Lifetimes for lowest 431 fine structure levels are provided. Energy levels and radiative rates are reported for transitions in F-like Ba XLVIII. Configuration interaction has been included among 27 configurations (generating 431 levels) over a wide energy range up to 618 Rydbergs, and the fully relativistic multi-configurational Dirac–Fock method adopted for the calculations. To assess the accuracy, calculations have also been performed with the flexible atomic code, FAC. Radiative rates, oscillator strengths and line strengths are reported for all electric dipole, magnetic dipole, electric quadrupole, and magnetic quadrupole transitions from the lowest 3 levels, although calculations have been performed for a much larger number of levels. We have made comparisons of our results with existing available results and a good agreement has been achieved. Additionally, lifetimes for all 431 levels are listed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Monte Carlo Shell Model for ab initio nuclear structure

We report on our recent application of the Monte Carlo Shell Model to no-core calculations. At the initial stage of the application, we have performed benchmark calculations in the p-shell region. Results are compared with those in the Full Configuration Interaction and No-Core Full Configuration methods. These are found to be consistent with each other within quoted uncertainties when they could be quantified. The preliminary results in N shell = 5 reveal the onset of systematic convergence pattern.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

Recent development of Monte Carlo shell model and its application to no-core calculations

One of the major challenges in nuclear theory is to reproduce and to predict nuclear structure from ab initio calculations with realistic nuclear forces. As the current limitation of direct diagonalization of Hamiltonian matrices by Lanczos iteration method is around the order of matrix dimensionality 10 10 in shell-model calculations, it is difficult to access heavier nuclei beyond the p shell with sufficiently large basis spaces. It is possible to overcome this difficulty by utilizing efficient approximate methods to reproduce full ab initio solutions with good precision and quantified uncertainties. Following the major success of the Monte Carlo shell model (MCSM) with an assumed inert core in the sd- and pf-shell regions and also by recent developments in the MCSM algorithm, the no-core MCSM is expected to be one of the most powerful tools to meet these conditions. We have performed benchmark calculations in the p-shell region. Results of energies are compared with those in the full configuration interaction and no-core full configuration methods. These are found to be consistent with each other within quoted uncertainties when they could be quantified. We also compare and discuss the radial density of the helium-4 ground state extracted from the MCSM and FCI many-body wave functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Optimized coordinates in vibrational coupled cluster calculations

The use of variationally optimized coordinates, which minimize the vibrational self-consistent field (VSCF) ground state energy with respect to orthogonal transformations of the coordinates, has recently been shown to improve the convergence of vibrational configuration interaction (VCI) towards the exact full VCI [K. Yagi, M. Keçeli, and S. Hirata, J. Chem. Phys. 137, 204118 (2012)]. The present paper proposes an incorporation of optimized coordinates into the vibrational coupled cluster (VCC), which has in the past been shown to outperform VCI in approximate calculations where similar restricted state spaces are employed in VCI and VCC. An embarrassingly parallel algorithm for variational optimization of coordinates for VSCF is implemented and the resulting coordinates and potentials are introduced into a VCC program. The performance of VCC in optimized coordinates (denoted oc-VCC) is examined through pilot applications to water, formaldehyde, and a series of water clusters (dimer, trimer, and hexamer) by comparing the calculated vibrational energy levels with those of the conventional VCC in normal coordinates and VCI in optimized coordinates. For water clusters, in particular, oc-VCC is found to gain orders of magnitude improvement in the accuracy, exemplifying that the combination of optimized coordinates localized to each monomer with the size-extensive VCC wave function provides a supreme description of systems consisting of weakly interacting sub-systems.

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

Electronic couplings for singlet fission: Orbital choice and extrapolation to the complete basis set limit

Abstract For the search for promising singlet fission candidates, the calculation of the effective electronic coupling, which is required to estimate the singlet fission rate between the initially excited state (S 0 S 1 ) and the multiexcitonic state ( 1 TT, two triplets on neighboring molecules, coupled into a singlet), should be sufficiently reliable and fast enough to explore the configuration space. We propose here to modify the calculation of the effective electronic coupling using a nonorthogonal configuration interaction approach by: (a) using only one set of orbitals, optimized for the triplet state of the molecules, to describe all molecular electronic states, and (b) only taking the leading configurations into consideration. Furthermore, we also studied the basis set convergence of the electronic coupling, and we found, by comparison to the complete basis set limit obtained using the cc‐pV n Z series of basis sets, that both the aug‐cc‐pVDZ and 6–311++G** basis sets are a good compromise between accuracy and computational feasibility. The proposed approach enables future work on larger clusters of molecules than dimers.

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

Radiative rates for E1, E2, M1, and M2 transitions in F-like ions with 55≤Z≤73

Energy levels, radiative rates and lifetimes are reported for 19 F-like ions with 55≤Z≤73, among 113 levels of the 2s{sup 2}2p{sup 5}, 2s2p{sup 6}, 2s{sup 2}2p{sup 4}3ℓ, 2s2p{sup 5}3ℓ, and 2p{sup 6}3ℓ configurations. The general-purpose relativistic atomic structure package (grasp) has been adopted for the calculations, and radiative rates (and other associated parameters, such as oscillator strengths and line strengths) are listed for all E1, E2, M1, and M2 transitions of the ions. Comparisons are made with earlier available theoretical and experimental energies, especially for Ba XLVIII. Nevertheless, calculations have also been performed with the flexible atomic code (fac), and with a much larger configuration interaction with up to 38 089 levels, for further accuracy assessments, particularly for energy levels.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Pair extended coupled cluster doubles

The accurate and efficient description of strongly correlated systems remains an important challenge for computational methods. Doubly occupied configuration interaction (DOCI), in which all electrons are paired and no correlations which break these pairs are permitted, can in many cases provide an accurate account of strong correlations, albeit at combinatorial computational cost. Recently, there has been significant interest in a method we refer to as pair coupled cluster doubles (pCCD), a variant of coupled cluster doubles in which the electrons are paired. This is simply because pCCD provides energies nearly identical to those of DOCI, but at mean-field computational cost (disregarding the cost of the two-electron integral transformation). Here, we introduce the more complete pair extended coupled cluster doubles (pECCD) approach which, like pCCD, has mean-field cost and reproduces DOCI energetically. We show that unlike pCCD, pECCD also reproduces the DOCI wave function with high accuracy. Moreover, pECCD yields sensible albeit inexact results even for attractive interactions where pCCD breaks down.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Energy decomposition analysis for exciplexes using absolutely localized molecular orbitals

An energy decomposition analysis (EDA) scheme is developed for understanding the intermolecular interaction involving molecules in their excited states. The EDA utilizes absolutely localized molecular orbitals to define intermediate states and is compatible with excited state methods based on linear response theory such as configuration interaction singles and time-dependent density functional theory. The shift in excitation energy when an excited molecule interacts with the environment is decomposed into frozen, polarization, and charge transfer contributions, and the frozen term can be further separated into Pauli repulsion and electrostatics. These terms can then be added to their counterparts obtained from the ground state EDA to form a decomposition of the total interaction energy. The EDA scheme is applied to study a variety of systems, including some model systems to demonstrate the correct behavior of all the proposed energy components as well as more realistic systems such as hydrogen-bonding complexes (e.g., formamide-water, pyridine/pyrimidine-water) and halide (F - ,Cl - )-water clusters that involve charge-transfer-to-solvent excitations.

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

Photodissociation spectroscopy of the dysprosium monochloride molecular ion

We have performed a combined experimental and theoretical study of the photodissociation cross section of the molecular ion DyCl{sup +}. The photodissociation cross section for the photon energy range 35 500 cm{sup −1} to 47 500 cm{sup −1} is measured using an integrated ion trap and time-of-flight mass spectrometer; we observe a broad, asymmetric profile that is peaked near 43 000 cm{sup −1}. The theoretical cross section is determined from electronic potentials and transition dipole moments calculated using the relativistic configuration-interaction valence-bond and coupled-cluster methods. The electronic structure of DyCl{sup +} is extremely complex due to the presence of multiple open electronic shells, including the 4f{sup 10} configuration. The molecule has nine attractive potentials with ionically bonded electrons and 99 repulsive potentials dissociating to a ground state Dy{sup +} ion and Cl atom. We explain the lack of symmetry in the cross section as due to multiple contributions from one-electron-dominated transitions between the vibrational ground state and several resolved repulsive excited states.

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

N-representability-driven reconstruction of the two-electron reduced-density matrix for a real-time time-dependent electronic structure method

Propagating the equations of motion (EOM) for the one-electron reduced-density matrix (1-RDM) requires knowledge of the corresponding two-electron RDM (2-RDM). We show that the indeterminacy of this expression can be removed through a constrained optimization that resembles the variational optimization of the ground-state 2-RDM subject to a set of known N-representability conditions. Electronic excitation energies can then be obtained by propagating the EOM for the 1-RDM and following the dipole moment after the system interacts with an oscillating external electric field. For simple systems with well-separated excited states whose symmetry differs from that of the ground state, excitation energies obtained from this method are comparable to those obtained from full configuration interaction computations. Although the optimized 2-RDM satisfies necessary N-representability conditions, the procedure cannot guarantee a unique mapping from the 1-RDM to the 2-RDM. This deficiency is evident in the mean-field-quality description of transitions to states of the same symmetry as the ground state, as well as in the inability of the method to describe Rabi oscillations.

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

Theoretical level energies and transition data for 4p{sup 6}4d{sup 5}, 4p{sup 5}4d{sup 6} and 4p{sup 6}4d{sup 4}4f configurations of W{sup 33+} ion

The ab initio quasirelativistic approach developed specifically for the calculation of spectral parameters of highly charged ions was used to derive transition data for the Nb-like tungsten ion W{sup 33+}. The configuration interaction method was applied to include electron correlation effects. The relativistic effects were taken into account in the Breit–Pauli approximation. The level energies, radiative lifetimes τ, Landé g-factors were determined for the ground configuration 4p{sup 6}4d{sup 5} and two excited configurations 4p{sup 5}4d{sup 6} and 4p{sup 6}4d{sup 4}4f. The radiative transition wavelengths λ and emission transition probabilities A for the electric dipole, electric quadrupole, electric octupole, magnetic dipole, and magnetic quadrupole transitions among the levels of these configurations were produced. The uncertainties of computed spectroscopic parameters were evaluated.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Quantifying charge resonance and multiexciton character in coupled chromophores by charge and spin cumulant analysis

We extend excited-state structural analysis to quantify the charge-resonance and multi-exciton character in wave functions of weakly interacting chromophores such as molecular dimers. The approach employs charge and spin cumulants which describe inter-fragment electronic correlations in molecular complexes. We introduce indexes corresponding to the weights of local, charge resonance, and biexciton (with different spin structure) configurations that can be computed for general wave functions thus allowing one to quantify the character of doubly excited states. The utility of the approach is illustrated by applications to several small dimers, e.g., He-H{sub 2}, (H{sub 2}){sub 2}, and (C{sub 2}H{sub 4}){sub 2}, using full and restricted configuration interaction schemes. In addition, we present calculations for several systems relevant to singlet fission, such as tetracene, 1,6-diphenyl-1,3,5-hexatriene, and 1,3-diphenylisobenzofuran dimers.

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

Energies and E1, M1, E2, and M2 transition rates for states of the 2s{sup 2}2p{sup 3}, 2s2p{sup 4}, and 2p{sup 5} configurations in nitrogen-like ions between F III and Kr XXX

Based on relativistic wavefunctions from multiconfiguration Dirac–Hartree–Fock and configuration interaction calculations, E1, M1, E2, and M2 transition rates, weighted oscillator strengths, and lifetimes are evaluated for the states of the (1s{sup 2})2s{sup 2}2p{sup 3},2s2p{sup 4}, and 2p{sup 5} configurations in all nitrogen-like ions between F III and Kr XXX. The wavefunction expansions include valence, core–valence, and core–core correlation effects through single–double multireference expansions to increasing sets of active orbitals. The computed energies agree very well with experimental values, with differences of only 300–600 cm{sup −1} for the majority of the levels and ions in the sequence. Computed transitions rates are in close agreement with available data from MCHF-BP calculations by Tachiev and Froese Fischer [G.I. Tachiev, C. Froese Fischer, A and A 385 (2002) 716].

74 ATOMIC AND MOLECULAR PHYSICS↗