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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Heat of formation determination of the ground and excited state of cyanomethylene (HCCN) radical

Ab initio electronic structure theory has been used to characterize the structure of the ground triplet and lowest singlet excited states of cyanomethylene. The geometries, vibrational frequencies, and heats of formation have been determined using second-order Moller-Plesset perturbation, single and double excitation configuration interaction, and quadratic configuration interaction theory. The heat of formation is predicted with isodesmic reaction and Gaussian-2 theory (G2) for the ground triplet and first excited singlet states of cyanomethylene. For the ground state Delta-H(sub 0)(sup f,0) is 114.8+/-2 kcal/mol while for the excited single state it is 126.5+/-2 kcal/mol.

Francisco, Joseph S.↗
DOE OSTI2017

Atomic structure considerations for the low-temperature opacity of Sn

Here, we have begun a preliminary investigation into the opacity of Sn at low temperatures (< 50 eV). The emissivity and opacity of Sn is a crucial factor in determining the utility of Sn in EUV lithography, with numerous industrial implications. To this end, we have been exploring the accuracy of some approximations used in opacity models for the relevant ion stages of Sn (neutral through ~ 18 times ionized). We also find that the use of intermediate-coupling, as compared to full configuration-interaction, is not adequate to obtain accurate line positions of the important bound-bound transitions in Sn. One requires full configuration-interaction to properly describe the strong mixing between the various n=4 sub-shells that give rise to the Δn= 0 transitions that dominate the opacity spectrum at low temperatures. Furthermore, since calculations that include full configuration-interaction for large numbers of configurations quickly become computationally prohibitive, we have explored hybrid calculations, in which full configuration-interaction is retained for the most important transitions, while intermediate-coupling is employed for all other transitions. After extensive exploration of the atomic structure properties, local-thermodynamic-equilibrium (LTE) opacities are generated using the ATOMIC code at selected temperatures and densities and compared to experiment.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Approximate singly excited states from a two-component Hartree-Fock reference

For many molecules, relaxing the spin symmetry constraint on the wave function results in the lowest energy mean-field solution. The two-component Hartree-Fock (2cHF) method relaxes all spin symmetry constraints, and the wave function is no longer an eigenfunction of the total spin, spin projection, or time-reversal symmetry operators. For ground state energies, 2cHF is a superior mean-field method for describing spin-frustrated molecules. For excited states, the utility of 2cHF is uncertain. Here, we implement the 2cHF extensions of two single-reference excited state methods, the two-component configuration interaction singles and time-dependent Hartree-Fock. We compare the results to the analogous methods based off of the unrestricted Hartree-Fock approximation, as well as the full configuration interaction for three small molecules with distinct 2cHF solutions, and discuss the nature of the 2cHF excited state solutions.

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

Radiative rates for E1, E2, M1, and M2 transitions in Br-like ions with 43 ≤ Z ≤ 50

Energies and lifetimes are reported for the eight Br-like ions with 43≤Z≤50, namely Tc IX, Ru X, Rh XI, Pd XII, Ag XIII, Cd XIV, In XV, and Sn XVI. Results are listed for the lowest 375 levels, which mostly belong to the 4s{sup 2}4p{sup 5}, 4s{sup 2}4p{sup 4}4ℓ, 4s4p{sup 6},4s{sup 2}4p{sup 4}5ℓ, 4s{sup 2}4p{sup 3}4d{sup 2}, 4s4p{sup 5}4ℓ, and 4s4p{sup 5}5ℓ configurations. Extensive configuration interaction among 39 configurations (generating 3990 levels) has been considered and the general-purpose relativistic atomic structure package (grasp) has been adopted for the calculations. Radiative rates are listed for all E1, E2, M1, and M2 transitions involving the lowest 375 levels. Previous experimental and theoretical energies are available for only a few levels of three, namely Ru X, Rh XI and Pd XII. Differences with the measured energies are up to 4% but the present results are an improvement (by up to 0.3 Ryd) in comparison to other recently reported theoretical data. Similarly for radiative rates and lifetimes, prior results are limited to those involving only 31 levels of the 4s{sup 2}4p{sup 5}, 4s{sup 2}4p{sup 4}4d, and 4s4p{sup 6} configurations for the last four ions. Moreover, there are generally no discrepancies with our results, although the larger calculations reported here differ by up to two orders of magnitude for a few transitions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Weak entanglement approximation for nuclear structure

The interacting shell model, a configuration-interaction method, is a venerable approach for low-lying nuclear structure calculations, but it is hampered by the exponential growth of its basis dimension as one increases the single-particle space and/or the number of active particles. Recent, quantum-information-inspired work has demonstrated that the proton and neutron sectors of a nuclear wave function are weakly entangled. Furthermore, the entanglement is smaller for nuclides away from N = Z, such as heavy, neutron-rich nuclides. Here, in this study, we implement a weak entanglement approximation to bipartite configuration-interaction wave functions, approximating low-lying levels by coupling a relatively small number of many-proton and many-neutron states. This truncation scheme, which we present in the context of past approaches, reduces the basis dimension by many orders of magnitude while preserving essential features of nuclear spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Quantum Package 2.0: An Open-Source Determinant-Driven Suite of Programs

Quantum chemistry is a discipline which relies heavily on very expensive numerical computations. The scaling of correlated wave function methods lies, in their standard implementation, between O(N-5) and O(e(N)), where N is proportional to the system size. Therefore, performing accurate calculations on chemically meaningful systems requires (i) approximations that can lower the computational scaling and (ii) efficient implementations that take advantage of modern massively parallel architectures. QUANTUM PACKAGE is an open-source programming environment for quantum chemistry specially designed for wave function methods. Its main goal is the development of determinant-driven selected configuration interaction (sCI) methods and multireference second-order perturbation theory (PT2). The determinant-driven framework allows the programmer to include any arbitrary set of determinants in the reference space, hence providing greater methodological freedom. The sCI method implemented in QUANTUM PACKAGE is based on the CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) algorithm which complements the variational sCI energy with a PT2 correction. Additional external plugins have been recently added to perform calculations with multireference coupled cluster theory and range-separated density-functional theory. All the programs are developed with the IRPF90 code generator, which simplifies collaborative work and the development of new features. QUANTUM PACKAGE strives to allow easy implementation and experimentation of new methods, while making parallel computation as simple and efficient as possible on modern supercomputer architectures. Currently, the code enables, routinely, to realize runs on roughly 2 000 CPU cores, with tens of millions of determinants in the reference space. Moreover, we have been able to push up to 12 288 cores in order to test its parallel efficiency. In the present manuscript, we also introduce some key new developments: (i) a renormalized second-order perturbative correction for efficient extrapolation to the full CI limit and (ii) a stochastic version of the CIPSI selection performed simultaneously to the PT2 calculation at no extra cost.

Garniron, Yann↗
DOE OSTI2024

Time propagation of electronic wavefunctions using nonorthogonal determinant expansions

The use of truncated configuration interaction in real-time time-dependent simulations of electron dynamics provides a balance of computational cost and accuracy, while avoiding some of the failures associated with real-time time-dependent density functional theory. However, low-order truncated configuration interaction also has limitations, such as overestimation of polarizability in configuration interaction singles, even when perturbative doubles are included. Increasing the size of the determinant expansion may not be computationally feasible, and so, in this work, we investigate the use of nonorthogonality in the determinant expansion to establish the extent to which higher-order substitutions can be recovered, providing an improved description of electron dynamics. Model systems are investigated to quantify the extent to which different methods accurately reproduce the (hyper)polarizability, including the high-harmonic generation spectrum of H 2 , water, and butadiene.

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

Atomic structure considerations for the low-temperature opacity of Xenon

Here, we have begun a preliminary investigation into the opacity of Xe at low temperatures (<50 eV). The emissivity and opacity of Xe is a crucial factor in determining the utility of Xe in EUV lithography, with numerous industrial applications. To this end, we have been exploring the accuracy of some approximations used in opacity models for the relevant ion stages of Xe (6 times ionized through 20 times ionized). Due to the complexity of Xe atomic structure, one needs to use full configuration-interaction to properly describe the strong mixing in the various n=4 sub-shells that give rise to the Δn=0 and Δn=1 transitions that dominate the opacity spectrum at low temperatures. Since calculations that include full configuration-interaction for large numbers of configurations quickly become computationally prohibitive, we have explored hybrid calculations, in which full configuration-interaction is retained for the most important transitions, while intermediate coupling is employed for all other transitions. After verifying atomic structure properties, local-thermodynamic-equilibrium (LTE) opacities are generated using the ATOMIC code at selected temperatures and densities and compared to various experimental results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Perspective on Kramers symmetry breaking and restoration in relativistic electronic structure methods for open-shell systems

Without rigorous symmetry constraints, solutions to approximate electronic structure methods may artificially break symmetry. In the case of the relativistic electronic structure, if time-reversal symmetry is not enforced in calculations of molecules not subject to a magnetic field, it is possible to artificially break Kramers degeneracy in open shell systems. This leads to a description of excited states that may be qualitatively incorrect. Despite this, different electronic structure methods to incorporate correlation and excited states can partially restore Kramers degeneracy from a broken symmetry solution. For single-reference techniques, the inclusion of double and possibly triple excitations in the ground state provides much of the needed correction. Formally, however, this imbalanced treatment of the Kramers-paired spaces is a multi-reference problem, and so methods such as complete-active-space methods perform much better at recovering much of the correct symmetry by state averaging. Using multi-reference configuration interaction, any additional corrections can be obtained as the solution approaches the full configuration interaction limit. A recently proposed “Kramers contamination” value is also used to assess the magnitude of symmetry breaking.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A well-scaling natural orbital theory

Here, we introduce an energy functional for ground-state electronic structure calculations. Its variables are the natural spin-orbitals of singlet many-body wave functions and their joint occupation probabilities deriving from controlled approximations to the two-particle density matrix that yield algebraic scaling in general, and Hartree–Fock scaling in its seniority-zero version. Results from the latter version for small molecular systems are compared with those of highly accurate quantum-chemical computations. The energies lie above full configuration interaction calculations, close to doubly occupied configuration interaction calculations. Their accuracy is considerably greater than that obtained from current density-functional theory approximations and from current functionals of the oneparticle density matrix.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS1983

The dissociation energy of Cu2 - Do we want to perform multi-reference singles and doubles CIs on many-electron systems?

Reference is made to a study by Hay and Dunning (1976), who showed that, for most problems, it is not necessary to perform full FOCI (first order configuration interaction); instead, only that subset of configuration state functions (CSFs) that differ by only one or two spin orbitals from any of a list of reference configurations need be used. This is referred to as POLCI. Calculations of the D(e) of Cu2 using both SDCI (singles and doubles configuration interaction) and POLCI models are presented. It is found that the size-consistency problems associated with SDCI cause an error in D(e), whereas the POLCI results are in excellent agreement with experiment. It is proposed that the orbitals be classified as inactive and active and that the inactive-inactive correlation be eliminated to reduce the size-consistency error. In addition, it is recommended that SDCI treatments be applied with caution to systems having a large number of electrons.

Bauschlicher, C. W., Jr.↗
DOE OSTI2019

Theoretical level energies and transition data for ion W{sup 28+}

The ab initio quasirelativistic approach was used to derive transition data for Pd-like tungsten W{sup 28+} ion. This approach was developed specifically for the calculation of spectral parameters for highly charged ions. The relativistic effects were taken into account in the Breit–Pauli approximation. The configuration interaction method was applied to include electron correlation effects. The ground configuration 4p{sup 6}4d{sup 10} and the lowest excited 4p{sup 6}4d{sup 9}4f were considered. The calculations for W{sup 28+} were performed with different sets of interacting configurations in order to evaluate reliability and accuracy of produced results. Level energies, radiative lifetimes τ, Landé g-factors were determined for the excited configuration 4p{sup 6}4d{sup 9}4f of W{sup 28+} levels. The radiative transition wavelengths λ and emission transition probabilities A for the electric and magnetic transitions of various multipole order (up to E5 and M4) among the levels of these two configurations were produced. The uncertainties of computed spectroscopic parameters were evaluated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

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

Theoretical study of the dipole moment of oxygen monofluoride (OF)

The ground-state potential curve and dipole-moment function of OF are calculated theoretically using the complete active-space self-consistent-field levels, externally contracted configuration-interaction levels, or multireference (singles plus doubles) configuration-interaction levels. Both an extended Gaussian basis set and a double-zeta-plus-polarization basis set were applied. The results are presented in extensive tables and graphs. Best results are achieved using a large Gaussian basis set and taking the valence-correlation energy into account. It is suggested that OF may best be detected by its laser-magnetic-resonance spectrum in the IR.

Langhoff, S. R.↗
NASA NTRS1993

Structure and thermochemistry of ClO2 radicals

The structure of ClO2 has been calculated for the X 2A-double prime ground state using unrestricted Hartree-Fock (UHF), unrestricted second-order Moller-Plesset perturbation (UMP2), configuration interaction employing single and double excitation (CISD), and quadratic configuration interaction (QCI) ab initio molecular orbital methods. Calculations using UMP2 and CISD wave functions predict a ClO bond length of 1.728 +/- 0.01 A. The single-configuration-based QCI in the singles and doubles space with perturbation inclusion of triple substitutions, denoted QCISD(T), yield a ClO bond length of 2.205 A. The QCI results are consistent with results of Jensen (1990) who showed that the ClO bond length is 2.181 A using annihilated self-consistent methods (AUMP2). The thermochemistry of ClO2 radical has been calculated using MP2 and QCI methods employing an isodesmic scheme. Our scheme predicts the heat of formation for ClO2 at 0 K to be 24.6 +/- 2 kcal/mol.

Francisco, J. S.↗
NASA NTRS2005

Research on Spectroscopy, Opacity, and Atmospheres

I propose to continue providing observers with basic data for interpreting spectra from stars, novas, supernovas, clusters, and galaxies. These data will include allowed and forbidden line lists, both laboratory and computed, for the first five to ten ions of all atoms and for all relevant diatomic molecules. I will eventually expand to all ions of the first thirty elements to treat far UV and X-ray spectra, and for envelope opacities. I also include triatomic molecules provided by other researchers. I have also made CDs with Partridge and Schwenke's water data for work on UV stars. The line data also serve as input to my model atmosphere and synthesis programs that generate energy distributions, photometry, limb darkening, and spectra that can be used for planning observations and for fitting observed spectra. The spectrum synthesis programs produce detailed plots with the lines identified. Grids of stellar spectra can be used for radial velocity-, rotation-, or abundance templates and for population synthesis. I am fitting spectra of bright stars to test the data and to produce atlases to guide observers. For each star the whole spectrum is computed from the UV to the far IR. The line data, opacities, models, spectra, and programs are freely distributed on CDs and on my Web site and represent a unique resource for many NASA programs. I am now in full production of new line lists for atoms. I am computing all ions of all elements from H to Zn and the first 5 ions of all the heavier elements, about 800 ions. For each ion I treat as many as 61 even and 61 odd configurations, computing all energy levels and eigenvectors. The Hamiltonian is determined from a scaled-Hartree-Fock starting guess by least squares fitting the observed energy levels. The average energy of each configuration is used in computing scaled-Thomas-Fermi-Dirac wavefunctions for each configuration which in turn are used to compute allowed and forbidden transition integrals. These are multiplied into the LS allowed and forbidden transition arrays. The transition arrays are transformed to the observed coupling to yield the allowed and forbidden line lists. Results are put on the web as they are finished. Provided I get funding,there will be more than 500 million lines. I will then compare ion by ion, to all the laboratory and computed data in the literature and make up a working line list for spectrum synthesis and opacity calculations with the best available data. As the laboratory spectrum analyses are improved, I will redo the calculations with the new energy levels. My original plan when I started the new calculations was to run through all the atoms using my old Cray programs from the 1980's that were limited to 1100 x 1100 arrays in the Hamiltonian for each J. Then I would go back and rerun the more complicated cases with 3000 x 3000 arrays so that I could include many more configurations and more configuration interactions. At present I am limited to 61 even and 61 odd configurations and I try to include everything up through n = 9. The current program runs on Alpha workstations. I decided to test the big program on Fe I and Fe II to see whether there was any great difference in the low configurations compared to those from the Cray program. Besides increasing the number of E1 lines by a factor of 6 to 7.7 million, there was an unexpected result: the electric quadrupole transitions were 10 times stronger than before because the transition integrals are weighted by r(exp 2) ---they become very large for high n, and because there are numerous configuration interactions that mix the low and high configurations. As a check I was able to reproduce Carstang's (1962) lower results by running his three configurations with my program. Since my model atom is still only a subset of a real Fe II ion, the true quadrupole A values are probably larger than mine. The magnetic dipole lines are affected by the mixing but the overall scale does not change. Because of this scovery I decided that there was no point in computing the small array cases. I have been running with as many configurations as I can and with thousands of parameters in the Hamiltonian. The computer runs take much longer to set up and produce than I had expected. I have concentrated on redoing the low iron group spectra, especially to get data for supernova modelers. I have done only Ca I -- Zn I, Ca II -- Zn II, CU I -- Cu XXIX, Zn I - Zn XXX, for practice at high stages of ionization, C I, C II, S I, and CL I and Ag I for people who were working on the laboratory spectra. Check my web site kurucz.harvard.edu for current additions. My latest calculations have been for carbon I and sulphur I, and silicon I is under way using the same elaborate approach as for C I, which took many months to do. These line lists greatly increase the number of lines in the ultraviolet, in the visible, and especially in the infrared. They will increase the opacity in A, F, and G stars. They will account for many unidentified lines in the sun.

Kurucz, Robert L.↗
DOE OSTI2019

Relativistic large scale CI calculations of energies, transition rates and lifetimes in Ca-like ions between Co VIII and Zn XI

Energies, E1, E2 and M1 transition rates and lifetimes are reported for the lowest 134 fine-structure levels of the 3p{sup 6}3d{sup 2}, 3p{sup 6}3d4s, 3p{sup 5}3d{sup 3} and 3p{sup 6}3d4p configurations in all Ca-like ions between Co VIII and Zn XI. The calculations have been performed using the multiconfiguration Dirac–Hartree–Fock (MCDHF) and relativistic configuration interaction (RCI) methods. Important core–valance and core–core electron correlation effects are systematically accounted with large scale configuration interaction expansion. The Breit interaction and leading quantum-electrodynamical (QED) contributions such as the vacuum polarization and self-energy corrections are also included in the RCI calculations. The results are compared with the compiled data from the NIST ASD and other calculations to evaluate the accuracy of the present calculations. Our results have good agreement with previous relativistic many-body perturbation theory (RMBPT) calculations and the data from the NIST ASD, except for Zn XI. Present calculations provide some missing data, especially for Cu X and Zn XI, in the NIST ASD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Communication: Towards ab initio self-energy embedding theory in quantum chemistry

The self-energy embedding theory (SEET), in which the active space self-energy is embedded in the self-energy obtained from a perturbative method treating the non-local correlation effects, was recently developed in our group. In SEET, the double counting problem does not appear and the accuracy can be improved either by increasing the perturbation order or by enlarging the active space. This method was first calibrated for the 2D Hubbard lattice showing promising results. In this paper, we report an extension of SEET to quantum chemical ab initio Hamiltonians for applications to molecular systems. The self-consistent second-order Green’s function method is used to describe the non-local correlations, while the full configuration interaction method is carried out to capture strong correlation within the active space. Using few proof-of-concept examples, we show that SEET yields results of comparable quality to n-electron valence state second-order perturbation theory with the same active space, and furthermore, the full active space can be split into smaller active spaces without further implementation. Moreover, SEET avoids intruder states and does not require any high-order reduced density matrices. These advantages show that SEET is a promising method to describe physical and chemical properties of challenging molecules requiring large active spaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗