Search NASA⌕ Search

SEARCH · Search NASA

Results for “configuration-interaction method”

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 37 records · Page 2

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↗

Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes

We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Triple and quadruple excitation contributions to the binding in Be clusters - Calibration calculations on Be3

The contribution of connected triple and quadruple excitations to the binding in Be3 is investigated by comparing various coupled-cluster (CC) and truncated configuration-interaction (CI) treatments with multireference CI (MRCI) and full CI (FCI) calculations. The CC method with single and double excitations (CCSD) produces results that differ substantially from more elaborate treatments, but most extensions to CCSD that account approximately for connected triple excitations perform very well. In contrast, good agreement with CFI for Be2 can be achieved only with the highest level CC and MRCI methods.

Watts, John D.↗

Full configuration-interaction benchmark calculations for AlH

Full CI and CASSCF/SOCI theoretical computations are presented for the X 1Sigma(+), A 1Pi, and C 1Sigma(+) states of AlH. The aim is to establish benchmark data to calibrate methods which describe valence-Rydberg transitions in other molecules. The results are presented in extensive tables and graphs and compared with published theoretical and experimental values, and the conditions under which CASSCF/SOCI calculations accurately reproduce the FCI results are defined.

Bauschlicher, Charles W., Jr.↗

The study of molecular spectroscopy by ab initio methods

This review illustrates the potential of theory for solving spectroscopic problems. The accuracy of approximate techniques for including electron correlation have been calibrated by comparison with full configuration-interaction calculations. Examples of the application of ab initio calculations to vibrational, rotational, and electronic spectroscopy are given. It is shown that the state-averaged, complete active space self-consistent field, multireference configuration-interaction procedure provides a good approach for treating several electronic states accurately in a common molecular orbital basis.

Bauschlicher, Charles W., Jr.↗

Electron density-sensitive line ratios for Fe XII in the solar corona

The collision strengths for electron-impact excitation of fine-structure transitions from the ground 3s23p3 configuration to the excited 3s3p4 configuration in Fe XII are presented. These are calculated by transforming the K-matrices obtained in a nonrelativistic R-matrix method to K-matrices in intermediate coupling using term-coupling coefficients which take into account the effect of breakdown of target LS coupling. Extensive configuration-interaction wave functions are used to represent the target states. The effective collision strengths are determined assuming a Maxwellian distribution of electron energies. Rydberg series of resonances converging to the excited state thresholds are explicitly included in the calculation. The relative level populations and emission-line strengths for Fe XII are calculated using the new atomic data presented in this paper and that published by Tayal, Henry, and Pradhan (1987) and Tayal and Henry (1986) over a wide range of electron densities and temperatures (1-2 x 10 to the 6th K). Present results for collision strengths and emission-line strengths are compared with the previous calculation of Flower (1977). Some significant differences are noted.

Tayal, S. S.↗

Cross Sections for Electron Impact Excitation of Ions Relevant to Planetary Atmospheres Observation

The goal of this research grant was to calculate accurate oscillator strengths and electron collisional excitation strengths for inelastic transitions in atomic species of relevance to Planetary Atmospheres. Large scale configuration-interaction atomic structure calculations have been performed to obtain oscillator strengths and transition probabilities for transitions among the fine-structure levels and R-matrix method has been used in the calculations of electron-ion collision cross sections of C II, S I, S II, S III, and Ar II. A number of strong features due to ions of sulfur have been detected in the spectra of Jupiter satellite Io. The electron excitation cross sections for the C II and S II transitions are studied in collaboration with the experimental atomic physics group at the Jet Propulsion Laboratory. There is excellent agreement between experiment and theory which provide an accurate and broad-base test of the ability of theoretical methods used in the calculation of atomic processes. Specifically, research problems have been investigated for: electron impact excitation cross sections of C II: electron impact excitation cross sections of S III; energy levels and oscillator strengths for transitions in S III; collision strengths for electron collisional excitation of S II; electron impact excitation of inelastic transitions in Ar II; oscillator strengths of fine-structure transitions in neutral sulfur; cross sections for inelastic scattering of electrons from atomic nitrogen; and excitation of atomic ions by electron impact.

Tayal, Swaraj S.↗

Theoretical study of the dipole moment function of OH(X 2Pi)

A complete theoretical electric dipole moment function mu(r) is computed for OH(X 2Pi) at the full configuration-interaction (FCI) level for DZP + diffuse Gaussian basis set. The energy and dipole moment are also determined with 2d functions on oxygen at seven r values. The following methods are compared to the FCI for both the oxygen 1d and 2d basis sets: single-reference singles-plus-doubles configuration interaction (SDCI), SDCI with Davidson's correction, the coupled pair functional approach (CPF) and a modified form of CPF, CASSCF, and CASSCF-MRSDCI, and the external contracted CI(CCI) method. The dipole moments are evaluated both as an expectation value and by finite-field methods. The results support the superiority of evaluating the dipole moment as an energy derivative rather than as an expectation value.

Langhoff, Stephen R.↗

The computation of C-C and N-N bond dissociation energies for singly, doubly, and triply bonded systems

The bond dissociation energies (De) of C2H2, C2H4, C2H6, N2, N2H2, and N2H4 are studied at various levels of correlation treatment. The convergence of De with respect to the one-particle basis is studied at the single-reference modified coupled-pair-functional (MCPF)level. At all levels of correlation treatment, the errors in the bond dissociation energies increase with the degree of multiple bond character. The multireference configuration-interaction (MRCI) De values, corrected for an estimate of higher excitations, are in excellent agreement with those determined using the size-extensive averaged-coupled-pair-functional (ACPF) method. The full-valence complete-active-space self-consistent-field (CASSCF)/MRCI calculations are reproduced very well by MRCI calculations based on a CASSCF calculation that includes in the active space only those electrons involved in the C-C or N-N bonds.

Langhoff, Stephen R.↗

Low-energy enhancement in the magnetic dipole γ-ray strength functions of heavy nuclei

A low-energy enhancement (LEE), which was observed experimentally in the gamma-ray strength function (γSF) describing the decay of compound nuclei, would have profound effects on r-process nucleosynthesis if it persists in heavy neutron-rich nuclei. The LEE was shown to be a feature of the magnetic dipole (M1) strength function in configuration-interaction shell-model calculations in medium-mass nuclei. However, its existence in heavy open-shell nuclei remains an open question. Here, using a combination of many-body methods, we identify a LEE in the M1 γSFs of heavy samarium nuclei. In particular, we use the static-path plus random-phase approximation (SPA+RPA), which includes static and small-amplitude quantal fluctuations beyond the mean field. Using the SPA+RPA strength as a prior, we apply the maximum-entropy method (MEM) to obtain finite-temperature M1 γSFs from exact imaginary-time response functions calculated with the shell model Monte Carlo (SMMC) method. We find that the slope of the LEE in samarium isotopes is roughly independent of the average initial energy over a wide range below the neutron separation energy. As the neutron number increases, strength transfers to a low-energy excitation, which we interpret as the scissors mode built on top of excited states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Core-core and core-valence correlation

The effect of 1s core correlation on properties and energy separations are analyzed using full configuration-interaction (FCI) calculations. The Be1S - 1P, the C 3P - 5S,m and CH(+) 1Sigma(+) - 1Pi separations, and CH(+) spectroscopic constants, dipole moment, and 1Sigma(+) - 1Pi transition dipole moment have been studied. The results of the FCI calculations are compared to those obtained using approximate methods.

Bauschlicher, Charles W., Jr.↗

Theoretical study of the hyperfine parameters of OH

In the present study of the hyperfine parameters of O-17H as a function of the one- and n-particle spaces, all of the parameters except oxygen's spin density, b sub F(O), are sufficiently easily tractable to allow concentration on the computational requirements for accurate determination of b sub F(O). Full configuration-interaction (FCI) calculations in six Gaussian basis sets yield unambiguous results for (1) the effect of uncontracting the O s and p basis sets; (2) that of adding diffuse s and p functions; and (3) that of adding polarization functions to O. The size-extensive modified coupled-pair functional method yields b sub F values which are in fair agreement with FCI results.

Chong, Delano P.↗

Core-core and core-valence correlation

The effect of (1s) core correlation on properties and energy separations was analyzed using full configuration-interaction (FCI) calculations. The Be 1 S - 1 P, the C 3 P - 5 S and CH+ 1 Sigma + or - 1 Pi separations, and CH+ spectroscopic constants, dipole moment and 1 Sigma + - 1 Pi transition dipole moment were studied. The results of the FCI calculations are compared to those obtained using approximate methods. In addition, the generation of atomic natural orbital (ANO) basis sets, as a method for contracting a primitive basis set for both valence and core correlation, is discussed. When both core-core and core-valence correlation are included in the calculation, no suitable truncated CI approach consistently reproduces the FCI, and contraction of the basis set is very difficult. If the (nearly constant) core-core correlation is eliminated, and only the core-valence correlation is included, CASSCF/MRCI approached reproduce the FCI results and basis set contraction is significantly easier.

Bauschlicher, Charles W., Jr.↗

An improved guess for the variational calculation of charge-transfer excitations in large systems

Ab initio quantum-chemical methods that perform well for computing the electronic ground state are not straightforwardly transferable to electronically excited states, particularly in large molecular systems. Wave function theory offers high accuracy, but is often prohibitively expensive. Methods based on time-dependent density functional theory (TD-DFT) are crucially sensitive to the chosen exchange-correlation functional (XCF) parameterization, and system-specific tuning protocols were therefore proposed to address the method's robustness. Methods based on the variational relaxation of the excited-state electron density showcased promising results for the calculation of charge-transfer excitations, but the complex shape of the electronic hypersurface makes convergence to a specific excited state much more difficult than for the ground state when standard variational techniques are applied. We address the latter aspect by providing suitable initial guesses, which we obtain by two separate constrained algorithms. Combined with the squared-gradient minimization algorithm for all-electrons relaxation in a freeze-and-release scheme (FRZ-SGM), we demonstrate that orbital-optimized density functional theory (OO-DFT) calculations can reliably converge to the charge-transfer states of interest even for large molecular systems. We test the FRZ-SGM method on a phenothiazine-anthraquinone CT excitation in a supramolecular Pd(II) coordination cage complex as a function of the cage conformation. This compound has been studied experimentally prior to our work. We compare this freeze-and-release scheme to two XCF reparameterizations, which were recently proposed as low-cost TD-DFT-based alternatives to variational methods. Two dye-semiconductor complexes, which were previously investigated in the context of photovoltaic applications, serve as a second example to investigate the convergence and stability of the FRZ-SGM approach. Our results demonstrate that FRZ-SGM provides reliable convergence for charge-transfer excited states and avoids variational collapse to lower-lying electronic states, whereas time-dependent DFT calculations with an adequate tuning procedure for the range-separation parameter provide a computationally efficient initial estimate of the corresponding energies, with a computational cost comparable to that of configuration-interaction singles (CIS) calculations.

Bogo, Nicola↗

Collision Strengths for Electron Collisional Excitation of S II

Electron collisional excitation strengths for inelastic transitions in S II are calculated using the R-matrix method in a 19-state (3s(sup 2)3p(sup 3)(sup 4)S(sup o), (sup 2)D(sup o), (sup 2)p(sup o), 3s3p(sup 4)(sup 4)P, (sup 2)D, (sup 2)S, 3S(sup 2)3p(sup 2)3d(sup 2)P, (sup 4)F, (sup 4)D, (sup 2)F, (sup 4)P, 3s(sup 2)3p(sup 2)4s(sup 4)P, (sup 2)P, 3s(sup 2)3p(sup 2)4p(sup 2)S(s o), (sup 4)D(sup o), (sup 4)P(sup o), (sup 2)D(sup o), (sup 4)S(sup o), (sup 2)P(sup o)) close-coupling approximation. These target states are represented by extensive configuration-interaction wave functions that give excitation energies and oscillator strengths that are usually in good agreement with the experimental values and the available accurate calculations. The present results for collision strengths are in very good agreement with the recent merged beams energy loss measurement of Liao et al. and agree reasonably well with the 18-state R-matrix calculation of Ramsbottom, Bell, & Stafford, but show significant differences from the 12-state R-matrix calculation of Cai & Pradhan.

Tayal, S. S.↗

FEM R-Matrix Calculations of Electron-H2 Collisions

Electron-H2 Collisions have been studied using the finite-element R-matrix method. Our approach uses a mixed finite-element and Gaussian basis to describe the continuum electron. In the inner region, the multi-centered nature of the Gaussian basis provides an efficient representation in the regions of space near the nuclei, whereas the piecemeal and energy-independent nature of the FEM basis is particular well suited for R-matrix calculations. Fixed nuclei calculations have been carried out at 1.4 a(sub 0), the equilibrium internuclear distance of the ground state. Up to six target states are included: the X(sup 1)E+(sub g), b(sup 3)E+(sub u), a(sup 3)E+(sub g), B(sup 1)E+(sub u), c(sup 3)II(sub u), and C(sup 1)II(sub u) states, and configuration-interaction functions are used to describe the target states. The results will be compared with available theoretical and experimental data.

Huo, Winifred M.↗

Collision Strengths for Electron Impact Excitation of Inelastic Transitions in Ar II

We have calculated collision strengths for electron impact excitation of inelastic transitions in Ar II using the R-matrix method in two independent nine- and 19-state close-coupling approximations. In the nine-state calculation the 3s(sup 2)3p(sup 5)2p(sup 0), 3S(sup 3)p(sup 6)S-2, 3p(sup 4)(P-3)4s(sup 2)P, 3p(sup 4)(P-3)3d(sup 2)P, 3p(sup 4)(D-1)4s(sup 2)D, 3p(sup 4)(P-3)3d(sup 2)D, 3p(sup 4)(S-1)4s(sup 2)S, 3p(sup 4)(D-1)3d(sup 2)S and 3p(sup 4)(D-1)4d(sup 2)S states are included, while in the 19-state calculation these states plus an additional ten states 3p(sup 4)(P-3)3d(sup 2)F, 3p(sup 4)(P-3)4p(sup 2)D(sup 0), 2p(sup 0), 2S(sup 0), 3p(sup 4)(D-1)4p(sup 2)P(sup 0), 2D(sup O), 2F(sup 0), 3p(sup 4)(D-1)3d(sup 2)D, p-2 and 3p(sup 4)(S-1)4p(sup 2)P(sup 0) are considered. These target states are represented by fairly extensive configuration-interaction wavefunctions which yield excitation energies and oscillator strengths that are generally in good agreement with the available most accurate calculations and the experimental values. Rydberg series of resonances converging to the excited state thresholds are included in the calculation. The effective collision strengths are obtained assuming a Maxwellian distribution of electron energies which are tabulated over the temperature range (0.5-20) x 10(exp 4) K.

Tayal, S. S.↗

High-Performance Semiempirical Excited-State Molecular Dynamics Powered by Graphics Processing Units

Here, this Letter introduces excited-state molecular dynamics in PYSEQM, a GPU-accelerated semiempirical quantum chemistry engine implemented in PyTorch. The new module enables Born–Oppenheimer molecular dynamics (BOMD) using configuration-interaction singles and random phase approximation for excited states, allowing long trajectories and large statistical ensembles to be simulated efficiently on a single GPU. We also implement an extended Lagrangian excited-state BOMD (XL-ESMD) scheme that propagates auxiliary electronic variables, enabling relaxed ground and excited-state convergence thresholds without compromising energy conservation. The excited-state BOMD implementation scales smoothly from small chromophores to a nearly 900-atom dendrimer (taking 6.5 s per MD step). PYSEQM also supports batched execution, allowing many geometries or trajectories to be evaluated in a single GPU launch, substantially increasing throughput and making ensemble-based protocols routine. As a demonstration, we compute absorption, emission, and infrared spectra from trajectories propagated on the ground and first excited states. The XL-ESMD scheme yields identical spectra at significantly lower computational cost, establishing the role of extended Lagrangian based dynamics for efficient excited-state BOMD simulations. Beyond raw performance, PYSEQM’s PyTorch foundation provides automatic differentiation for forces, efficient GPU batching, and seamless interfacing with machine learning models. These capabilities position PYSEQM as a practical platform for machine learning-augmented excited-state dynamics and lay the foundation for future data-driven nonadiabatic excited-state dynamics modeling of ultrafast spectroscopic probes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗