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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Vibrational heat-bath configuration interaction with semistochastic perturbation theory using harmonic oscillator or VSCF modals

Vibrational heat-bath configuration interaction (VHCI)—a selected configuration interaction technique for vibrational structure theory—has recently been developed in two independent works [J. H. Fetherolf and T. C. Berkelbach, J. Chem. Phys. 154, 074104 (2021); A. U. Bhatty and K. R. Brorsen, Mol. Phys. 119, e1936250 (2021)], where it was shown to provide accuracy on par with the most accurate vibrational structure methods with a low computational cost. Here, we eliminate the memory bottleneck of the second-order perturbation theory correction using the same (semi)stochastic approach developed previously for electronic structure theory. This allows us to treat, in an unbiased manner, much larger perturbative spaces, which are necessary for high accuracy in large systems. Stochastic errors are easily controlled to be less than 1 cm−1. We also report two other developments: (i) we propose a new heat-bath criterion and an associated exact implicit sorting algorithm for potential energy surfaces expressible as a sum of products of one-dimensional potentials; (ii) we formulate VHCI to use a vibrational self-consistent field (VSCF) reference, as opposed to the harmonic oscillator reference configuration used in previous reports. Our tests are done with quartic and sextic force fields, for which we find that with VSCF, the minor improvements to accuracy are outweighed by the higher computational cost associated the matrix element evaluations. We expect VSCF-based VHCI to be important for more general potential representations, for which the harmonic oscillator basis function integrals are no longer analytic.

Chemistry↗
DOE OSTI2014

A hybrid configuration interaction treatment based on seniority number and excitation schemes

We present a configuration interaction method in which the Hamiltonian of an N-electron system is projected on Slater determinants selected according to the seniority-number criterion along with the traditional excitation-based procedure. This proposed method is especially useful to describe systems which exhibit dynamic (weak) correlation at determined geometric arrangements (where the excitation-based procedure is more suitable) but show static (strong) correlation at other arrangements (where the seniority-number technique is preferred). The hybrid method amends the shortcomings of both individual determinant selection procedures, yielding correct shapes of potential energy curves with results closer to those provided by the full configuration interaction method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

GronOR: Massively Parallel and GPU-Accelerated Non-Orthogonal Configuration Interaction for Large Molecular Systems

GronOR is a program package for non-orthogonal configuration interaction calculations for an electronic wave function built in terms of anti-symmetrized products of multi-configuration molecular fragment wave functions. The two-electron integrals that have to be processed may be expressed in terms of atomic orbitals or in terms of an orbital basis determined from the molecular orbitals of the fragments. The code has been specifically designed for execution on distributed memory massively parallel and Graphics Processing Unit (GPU)-accelerated computer architectures, using an MPI+OpenACC/OpenMP programming approach. The task-based execution model used in the implementation allows for linear scaling with the number of nodes on the largest pre-exascale architectures available, provides hardware fault resiliency, and enables effective execution on systems with distinct central processing unit-only and GPU-accelerated partitions. The code interfaces with existing multi-configuration electronic structure codes that provide optimized molecular fragment orbitals, configuration interaction coefficients, and the required integrals. Algorithm and implementation details, parallel and accelerated performance benchmarks, and an analysis of the sensitivity of the accuracy of results and computational performance to thresholds used in the calculations are presented.

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

Intersystem-crossing and phosphorescence rates in fac-Ir{sup III}(ppy){sub 3}: A theoretical study involving multi-reference configuration interaction wavefunctions

We have employed combined density functional theory and multi-reference configuration interaction methods including spin–orbit coupling (SOC) effects to investigate the photophysics of the green phosphorescent emitter fac-tris-(2-phenylpyridine)iridium (fac-Ir(ppy){sub 3}). A critical evaluation of our quantum chemical approaches shows that a perturbational treatment of SOC is the method of choice for computing the UV/Vis spectrum of this heavy transition metal complex while multi-reference spin–orbit configuration interaction is preferable for calculating the phosphorescence rates. The particular choice of the spin–orbit interaction operator is found to be of minor importance. Intersystem crossing (ISC) rates have been determined by Fourier transformation of the time correlation function of the transition including Dushinsky rotations. In the electronic ground state, fac-Ir(ppy){sub 3} is C{sub 3} symmetric. The calculated UV/Vis spectrum is in excellent agreement with experiment. The effect of SOC is particularly pronounced for the metal-to-ligand charge transfer (MLCT) band in the visible region of the absorption spectrum which does not only extend its spectral onset towards longer wavelengths but also experiences a blue shift of its maximum. Pseudo-Jahn-Teller interaction leads to asymmetric coordinate displacements in the lowest MLCT states. Substantial electronic SOC and a small energy gap make ISC an ultrafast process in fac-Ir(ppy){sub 3}. For the S{sub 1}↝T{sub 1} non-radiative transition, we compute a rate constant of k{sub ISC} = 6.9 × 10{sup 12} s{sup −1} which exceeds the rate constant of radiative decay to the electronic ground state by more than six orders of magnitude, in agreement with the experimental observation of a subpicosecond ISC process and a triplet quantum yield close to unity. As a consequence of the geometric distortion in the T{sub 1} state, the T{sub 1} → S{sub 0} transition densities are localized on one of the phenylpyridyl moieties. In our best quantum chemical model, we obtain phosphorescence decay times of 264 μs, 13 μs, and 0.9 μs, respectively, for the T{sub 1,I}, T{sub 1,II}, and T{sub 1,III} fine-structure levels in dichloromethane (DCM) solution. In addition to reproducing the correct orders of magnitude for the individual phosphorescence emission probabilities, our theoretical study gives insight into the underlying mechanisms. In terms of intensity borrowing from spin-allowed transitions, the low emission probability of the T{sub 1,I} substate is caused by the mutual cancellation of contributions from several singlet states to the total transition dipole moment. Their contributions do not cancel but add up in case of the much faster T{sub 1,III} → S{sub 0} emission while the T{sub 1,II} → S{sub 0} emission is dominated by intensity borrowing from a single spin-allowed process, i.e., the S{sub 2} → S{sub 0} transition.

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

Modern Approaches to Exact Diagonalization and Selected Configuration Interaction with the Adaptive Sampling CI Method

Recent advances in selected configuration interaction methods have made them competitive with the most accurate techniques available and, hence, creating an increasingly powerful tool for solving quantum Hamiltonians. In this work, we build on recent advances from the adaptive sampling configuration interaction (ASCI) algorithm. We show that a useful paradigm for generating efficient selected CI/exact diagonalization algorithms is driven by fast sorting algorithms, much in the same way iterative diagonalization is based on the paradigm of matrix vector multiplication. We present several new algorithms for all parts of performing a selected CI, which includes new ASCI search, dynamic bit masking, fast orbital rotations, fast diagonal matrix elements, and residue arrays. The ASCI search algorithm can be used in several different modes, which includes an integral driven search and a coefficient driven search. The algorithms presented here are fast and scalable, and we find that because they are built on fast sorting algorithms they are more efficient than all other approaches we considered. After introducing these techniques, we present ASCI results applied to a large range of systems and basis sets to demonstrate the types of simulations that can be practically treated at the full-CI level with modern methods and hardware, presenting double- and triple-ζ benchmark data for the G1 data set. The largest of these calculations is Si$_2$H$_6$ which is a simulation of 34 electrons in 152 orbitals. We also present some preliminary results for fast deterministic perturbation theory simulations that use hash functions to maintain high efficiency for treating large basis sets.

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

Efficient Treatment of Large Active Spaces through Multi-GPU Parallel Implementation of Direct Configuration Interaction

In this study, we have extended our graphical processing unit (GPU)-accelerated direct configuration interaction program to multiple devices, reducing iteration times for configuration spaces of 165 million determinants to only 3 s using NVIDIA P100 GPUs. Similar improvements in the one- and two-particle reduced density matrix formation allow for fast analytical energy gradients and electronic properties. Our parallel algorithm enables the calculation of arbitrarily large configuration spaces (limited only by available system memory), with iteration times of 13 min for an active space of 18 electrons in 18 orbitals (2.4 billion determinants) using six consumer grade NVIDIA 1080Ti GPUs. These advances enable routine molecular dynamics simulations, geometry optimizations, and absorption spectrum calculations for molecules with large configuration spaces, a task that has heretofore required massive computational effort. In this work, we demonstrate the utility of our program by generating the absorption spectrum for diphenyl acetylene at the floating occupation molecular orbital complete active space configuration interaction level of theory. Lastly, several active spaces were investigated to assess the dependence of spectral features on orbital space dimension.

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

Approaching exact hyperpolarizabilities via sum-over-states Monte Carlo configuration interaction

We propose using sum-over-states calculations with the compact wavefunctions of Monte Carlo configuration interaction to approach accurate values for higher-order dipole properties up to second hyperpolarizabilities in a controlled way. We apply the approach to small systems that can generally be compared with full configuration interaction (FCI) results. We consider hydrogen fluoride with a 6-31g basis and then look at results, including frequency dependent properties, in an aug-cc-pVDZ basis. We extend one calculation beyond FCI by using an aug-cc-pVTZ basis. The properties of an H{sub 4} molecule with multireference character are calculated in an aug-cc-pVDZ basis. We then investigate this method on a strongly multireference system with a larger FCI space by modelling the properties of carbon monoxide with a stretched geometry. The behavior of the approach with increasing basis size is considered by calculating results for the neon atom using aug-cc-pVDZ to aug-cc-pVQZ. We finally test if the unusual change in polarizability between the first two states of molecular oxygen can be reproduced by this method in a 6-31g basis.

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

Extended relativistic configuration interaction and many-body perturbation calculations of spectroscopic data for the N ≤ 6 configurations in Ne-like ions between Cr xv and Kr xxvii

Level energies, wavelengths, electric dipole, magnetic dipole, electric quadrupole, and magnetic quadrupole transition rates, oscillator strengths, and line strengths from combined relativistic configuration interaction and many-body perturbation calculations are reported for the 201 fine-structure states of the 2s{sup 2}2p{sup 6}, 2s{sup 2}2p{sup 5}3l, 2s2p{sup 6}3l, 2s{sup 2}2p{sup 5}4l, 2s2p{sup 6}4l, 2s{sup 2}2p{sup 5}5l, and 2s{sup 2}2p{sup 5}6l configurations in all Ne-like ions between Cr xv and Kr xxvii. Calculated level energies and transition data are compared with experiments from the National Institute of Standards and Technology (NIST) and CHIANTI databases, and other recent benchmark calculations. The mean energy difference with the NIST experiments is only 0.05%. The present calculations significantly increase the amount of accurate spectroscopic data for the n > 3 states in a number of Ne-like ions of astrophysical interest. A complete data set should be helpful for analyzing new observations from solar and other astrophysical sources, and is also likely to be useful for modeling and diagnosing a variety of plasmas, including astronomical and fusion plasma.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1988

Full configuration interaction benchmark calculations for TiH

Full configuration interaction (FCI) results for the 3F and 5F states of Ti and the 4Phi and 2Delta states of TiH are presented. While the coupled pair functional (CPF) or modified CPF approaches are found to work well for the 4Phi state of TiH, they do not perform as well for the 2Delta state. Although for mu, the CASSCF/MRCI methods do well for the 2Delta state, when the active space is expanded outside the normal valence definition, the dipole moment is only brought into agreement with the FCI results by natural orbital iterations.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2014

Fitting coupled potential energy surfaces for large systems: Method and construction of a 3-state representation for phenol photodissociation in the full 33 internal degrees of freedom using multireference configuration interaction determined data

A recently reported algorithm for representing adiabatic states coupled by conical intersections using a quasi-diabatic state Hamiltonian in four and five atom systems is extended to treat nonadiabatic processes in considerably larger molecules. The method treats all internal degrees of freedom and uses electronic structure data from ab initio multireference configuration interaction wave functions with nuclear configuration selection based on quasi-classical surface hopping trajectories. The method is shown here to be able to treat ∼30 internal degrees of freedom including dissociative and large amplitude internal motion. Two procedures are introduced which are essential to the algorithm, a null space projector which removes basis functions from the fitting process until they are needed and a partial diagonalization technique which allows for automated, but accurate, treatment of the vicinity of extended seams of conical intersections of two or more states. These procedures are described in detail. The method is illustrated using the photodissociaton of phenol, C{sub 6}H{sub 5}OH(X{sup ~1}A{sup ′}) + hv → C{sub 6}H{sub 5}OH(A{sup ~1}A{sup ′}, B{sup ~1}A{sup ′′}) → C{sub 6}H{sub 5}O(X{sup ~2}B{sub 1}, A{sup ~2}B{sub 2}) + H as a test case. Ab initio electronic structure data for the 1,2,3{sup 1}A states of phenol, which are coupled by conical intersections, are obtained from multireference first order configuration interaction wave functions. The design of bases to simultaneously treat large amplitude motion and dissociation is described, as is the ability of the fitting procedure to smooth the irregularities in the electronic energies attributable to the orbital changes that are inherent to nonadiabatic processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Analytic nuclear forces and molecular properties from full configuration interaction quantum Monte Carlo

Unbiased stochastic sampling of the one- and two-body reduced density matrices is achieved in full configuration interaction quantum Monte Carlo with the introduction of a second, “replica” ensemble of walkers, whose population evolves in imaginary time independently from the first and which entails only modest additional computational overheads. The matrices obtained from this approach are shown to be representative of full configuration-interaction quality and hence provide a realistic opportunity to achieve high-quality results for a range of properties whose operators do not necessarily commute with the Hamiltonian. A density-matrix formulated quasi-variational energy estimator having been already proposed and investigated, the present work extends the scope of the theory to take in studies of analytic nuclear forces, molecular dipole moments, and polarisabilities, with extensive comparison to exact results where possible. These new results confirm the suitability of the sampling technique and, where sufficiently large basis sets are available, achieve close agreement with experimental values, expanding the scope of the method to new areas of investigation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Configuration–Interaction Perturbation Theory Calculations of Pu II

Configuration–interaction perturbation theory (CI–PT) is applied to calculations of low-energy states of Pu II. This ion is quite challenging due to a large number of possible determinants arising from seven valence electrons and strong relativistic effects. The CI–PT calculations agree with experiments for the energies and g-factors for many low-energy states that allowed positive identification of the theoretical levels. Isotope shifts were also used to aid in identification, and, in case of the odd states, fitting with three independent parameters was used to match theoretical isotope shifts to the experimental values with good accuracy. The CI–PT approach tested here on the Pu II ion can be generally used to calculate properties of many complex atoms, including U I that can find application in fundamental and applied science.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Fixed-node diffusion Monte Carlo potential energy curve of the fluorine molecule F{sub 2} using selected configuration interaction trial wavefunctions

The potential energy curve of the F{sub 2} molecule is calculated with Fixed-Node Diffusion Monte Carlo (FN-DMC) using Configuration Interaction (CI)-type trial wavefunctions. To keep the number of determinants reasonable and thus make FN-DMC calculations feasible in practice, the CI expansion is restricted to those determinants that contribute the most to the total energy. The selection of the determinants is made using the CIPSI approach (Configuration Interaction using a Perturbative Selection made Iteratively). The trial wavefunction used in FN-DMC is directly issued from the deterministic CI program; no Jastrow factor is used and no preliminary multi-parameter stochastic optimization of the trial wavefunction is performed. The nodes of CIPSI wavefunctions are found to reduce significantly the fixed-node error and to be systematically improved upon increasing the number of selected determinants. To reduce the non-parallelism error of the potential energy curve, a scheme based on the use of a R-dependent number of determinants is introduced. Using Dunning’s cc-pVDZ basis set, the FN-DMC energy curve of F{sub 2} is found to be of a quality similar to that obtained with full configuration interaction/cc-pVQZ.

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

A combined DFT and restricted open-shell configuration interaction method including spin-orbit coupling: Application to transition metal L-edge X-ray absorption spectroscopy

A novel restricted-open-shell configuration interaction with singles (ROCIS) approach for the calculation of transition metal L-edge X-ray absorption spectra is introduced. In this method, one first calculates the ground state and a number of excited states of the non-relativistic Hamiltonian. By construction, the total spin is a good quantum number in each of these states. For a ground state with total spin S excited states with spin S Prime = S, S - 1, and S + 1 are constructed. Using Wigner-Eckart algebra, all magnetic sublevels with M{sub S}= S, Horizontal-Ellipsis , -S for each multiplet of spin S are obtained. The spin-orbit operator is represented by a mean-field approximation to the full Breit-Pauli spin-orbit operator and is diagonalized over this N-particle basis. This is equivalent to a quasi-degenerate treatment of the spin-orbit interaction to all orders. Importantly, the excitation space spans all of the molecular multiplets that arise from the atomic Russell-Saunders terms. Hence, the method represents a rigorous first-principles approach to the complicated low-symmetry molecular multiplet problem met in L-edge X-ray absorption spectroscopy. In order to gain computational efficiency, as well as additional accuracy, the excitation space is restricted to single excitations and the configuration interaction matrix is slightly parameterized in order to account for dynamic correlation effects in an average way. To this end, it is advantageous to employ Kohn-Sham rather than Hartree-Fock orbitals thus defining the density functional theory/ROCIS method. However, the method can also be used in an entirely non-empirical fashion. Only three global empirical parameters are introduced and have been determined here for future application of the method to any system containing any transition metal. The three parameters were carefully calibrated using the L-edge X-ray absorption spectroscopy spectra of a test set of coordination complexes containing first row transition metals. These parameters are universal and transferable. Hence, there are no adjustable parameters that are used to fit experimental X-ray absorption spectra. Thus, the new approach classifies as a predictive first-principles method rather than an analysis tool. A series of calculations on transition metal compounds containing Cu, Ti, Fe, and Ni in various oxidation and spin states is investigated and a detailed comparison to experimental data is reported. In most cases, the approach yields good to excellent agreement with experiment. In addition, the origin of the observed spectral features is discussed in terms of the electronic structure of the investigated compounds.

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

Ab initio molecular orbital-configuration interaction based quantum master equation (MOQME) approach to the dynamic first hyperpolarizabilities of asymmetric π-conjugated systems

In this study, we apply the ab initio molecular orbital - configuration interaction based quantum master equation (MOQME) approach to the calculation and analysis of the dynamic first hyperpolarizabilities (β) of asymmetric π-conjugated molecules. In this approach, we construct the excited state models by the ab initio configuration interaction singles method. Then, time evolutions of system reduced density matrix ρ(t) and system polarization p(t) are calculated by the QME approach. Dynamic β in the second harmonic generation is calculated based on the nonperturbative definition of nonlinear optical susceptibility, using the frequency domain system polarization p(ω). Spatial contributions of electrons to β are analyzed based on the dynamic hyperpolarizability density map, which visualizes the second-order response of charge density oscillating with a frequency of 2ω. We apply the present method to the calculation of the dynamic β of a series of donor/acceptor substituted polyene oligomers, and then discuss the applicability of the MOQME method to the calculation and analysis of dynamic NLO properties of molecular systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Improved Fast Randomized Iteration Approach to Full Configuration Interaction

In this work, we present three modifications to our recently introduced fast randomized iteration method for full configuration interaction (FCI-FRI) and investigate their effects on the method’s performance for Ne, H 2 O, and N 2 . The initiator approximation, originally developed for full configuration interaction quantum Monte Carlo, significantly reduces statistical error in FCI-FRI when few samples are used in compression operations, enabling its application to larger chemical systems. The semistochastic extension, which involves exactly preserving a fixed subset of elements in each compression, improves statistical efficiency in some cases but reduces it in others. We also developed a new approach to sampling excitations that yields consistent improvements in statistical efficiency and reductions in computational cost. We discuss possible strategies based on our findings for improving the performance of stochastic quantum chemistry methods more generally.

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

Decomposition of the configuration-interaction coefficients in the multiconfiguration time-dependent Hartree-Fock method

An approximate implementation of the multiconfiguration time-dependent Hartree-Fock method is proposed, in which the matrix of configuration-interaction coefficients is decomposed into a product of matrices of smaller dimension. The applicability of this method in which all the configurations are kept in the expansion of the wave function, while the configuration-interaction coefficients are approximately calculated, is discussed by showing the results on three model systems: a one-dimensional model of a beryllium atom, a one-dimensional model of a carbon atom, and a one-dimensional model of a chain of four hydrogen atoms. The time-dependent electronic dynamics induced by a few-cycle, long-wavelength laser pulse is found to be well described at a lower computational cost compared to the standard multiconfiguration time-dependent Hartree-Fock treatment. Drawbacks of the method are also discussed.

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

Highly correlated configuration interaction calculations on water with large orbital bases

A priori selected configuration interaction (SCI) with truncation energy error [C. F. Bunge, J. Chem. Phys. 125, 014107 (2006)] and CI by parts [C. F. Bunge and R. Carbó-Dorca, J. Chem. Phys. 125, 014108 (2006)] are used to approximate the total nonrelativistic electronic ground state energy of water at fixed experimental geometry with CI up to sextuple excitations. Correlation-consistent polarized core-valence basis sets (cc-pCVnZ) up to sextuple zeta and augmented correlation-consistent polarized core-valence basis sets (aug-cc-pCVnZ) up to quintuple zeta quality are employed. Truncation energy errors range between less than 1 μhartree, and 100 μhartree for the largest orbital set. Coupled cluster CCSD and CCSD(T) calculations are also obtained for comparison. Our best upper bound, −76.4343 hartree, obtained by SCI with up to sextuple excitations with a cc-pCV6Z basis recovers more than 98.8% of the correlation energy of the system, and it is only about 3 kcal/mol above the “experimental” value. Despite that the present energy upper bounds are far below all previous ones, comparatively large dispersion errors in the determination of the extrapolated energies to the complete basis set do not allow to determine a reliable estimation of the full CI energy with an accuracy better than 0.6 mhartree (0.4 kcal/mol)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗