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

Removing Basis Set Incompleteness Error in Finite-Temperature Electronic Structure Calculations: Two-Electron Systems

We investigate the basis-set-size dependence for quantities related to interacting electrons in the canonical ensemble. Calculations are performed using exact diagonalization (finite temperature full configuration interaction method) on two-electron model systems–the uniform electron gas (UEG) and the helium atom. Our data reproduce previous observations of a competition for how the internal energy converges between the ground-state correlation energy and the high-temperature kinetic energy. We explore how this can be related to component parts of the internal energy including kinetic, exchange, and correlation energies and show there is surprising nuance in how this can be broken down into mostly monotonically converging quantities. We also show that separation of the free energy into a free energy with/without correlation allows for monotonic convergence with basis set size due to the variational principle. We find that the free energy convergence matches the previously observed convergence properties of the internal energy. We discuss the free energy divergence that happens when converging a finite basis analytical hydrogen atom to the complete basis set limit and compare this to the energies of a helium atom in a large periodic box. Reducing the box size, we saw convergence trends for the helium atom that were similar to the UEG.

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

The many-body expansion for metals. I. The alkaline earth metals Be, Mg, and Ca

We examine the many-body expansion (MBE) for alkaline earth metal clusters, Be n , Mg n , Ca n ( n = 4, 5, 6), at the Møller–Plesset second order perturbation theory, coupled-cluster singles and doubles with perturbative triples, multi-reference perturbation theory, and multi-reference configuration interaction levels of theory. The magnitude of each term in the MBE is evaluated for several geometrical configurations. We find that the behavior of the MBE for these clusters depends strongly on the geometrical arrangement and, to a lesser extent, on the level of theory used. Another factor that affects the MBE is the in situ (ground or excited) electronic state of the individual atoms in the cluster. For most geometries, the three-body term is the largest, followed by a steady decrease in absolute energy for subsequent terms. Though these systems exhibit non-negligible multi-reference effects, there was little qualitative difference in the MBE when employing single vs multi-reference methods. Useful insights into the connectivity and stability of these clusters have been drawn from the respective potential energy surfaces and quasi-atomic orbitals for the various dimers, trimers, and tetramers. Through these analyses, we investigate the similarities and differences in the binding energies of different-sized clusters for these metals.

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

Measures of complexity and entanglement in many-fermion systems

There is no unique and widely accepted definition of the complexity measure (CM) of a many-fermion wave function in the presence of interactions. The simplest many-fermion wave function is a Slater determinant. In shell-model or configuration interaction (CI) and other related methods, the state is represented as a superposition of a large number of Slater determinants, which in the case of CI calculations reaches about 20 billion terms [Johnson, arXiv:1809.07869]. Although in practice this number has been used as a CM for decades, it is ill defined: it is not unique, and it depends on the particular type and the number of single-particle wave functions used to construct the Slater determinants. Further, the canonical wave functions and/or natural orbitals [Löwdin, Adv. Phys. 5, 1 (1956); Löwdin and Shull, Phys. Rev. 101, 1730 (1956); Bardeen et al., Phys. Rev. 108, 1175 (1957); N. N. Bogoljubov, Il Nuovo Cimento 7, 794 (1958); Valatin, Il Nuovo Cimento 7, 843 (1958); de Gennes, Superconductivity of Metals and Alloys (CRC Press, Boca Raton, FL, 1999); Ring and Schuck, The Nuclear Many-Body Problem, 1st ed. (Springer-Verlag, Berlin, 2004)] and their corresponding occupation probabilities are intrinsic properties of any many-body wave function, irrespective of the representation, and they provide a unique solution to characterize the CM. The non-negative orbital entanglement entropy, which vanishes for a Slater determinant, provides the simplest CM, while a more complete measure of complexity is the entanglement spectrum. We illustrate these aspects in the case of a complex nonequilibrium time-dependent process, induced nuclear fission described within a real-time density functional theory framework extended to superfluid systems, which can describe simultaneously the long-range and the short-range correlations between fermions. The orbital entanglement entropy of the fissioning nucleus illustrates the localization mechanism of the many-body wave function in Fock and/or Hilbert space. The (minimal) number of Slater determinants required to represent such a complex many-body wave function with a well-defined number of particles in the case presented here is about 10 500 . The realistic case of the highly nonequilibrium nuclear fission process illustrated here is equivalent to a system of 23.328×10 9 interacting quantum spin-1/2 particles, a very large system for the study of quantum entanglement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Electromagnetic moments of 215,217 Bi: Probing shell evolution beyond N = 126

The nuclear properties of bismuth isotopes (Z = 83) , with just one valence proton above the closed spherical shell at (Z = 82) , are expected to be governed by a single unpaired proton. However, already in semimagic 209 Bi (Z = 83, N = 126) , , the magnetic moment (μ) strongly deviates from the single-particle Schmidt value. A near linear decrease in μ with the increase of N after the N = 126 magic number was observed up to N = 130 . In order to test whether this trend is kept at N > 130 and to reveal the underlying mechanisms, an investigation of 215,217 Bi (N = 132, 134) has been undertaken. The magnetic dipole and electric quadrupole moments of the I π = 9/2 - nuclear ground states in these isotopes have been measured for the first time using the in-source resonance-ionization spectroscopy technique at ISOLDE (CERN). It has been shown that the linearly decreasing trend of μ( 209,211,213 Bi g ) is broken in 215,217 Bi with a nearly constant value of μ observed. Experimental data have been compared to calculations in the framework of the configuration-interaction shell model with the monopole-based universal V MU +LS interaction. The peculiarities in the behavior of μ(Bi, 9/2 - ) with increasing neutron number are explained as being due to the shell evolution, change of the neutron orbitals occupancies and strong configuration mixing beyond N = 130 . Also, the difference in the μ trends for bismuth (Z = 83) and astatine (Z = 85) isotopes with N > 126 are reproduced by the shell-model calculations. It is shown that monopole interaction plays noticeable role in the description of the peculiarities of the μ behaviour. Additionally, the extension of the application of the V MU interaction to the μ isotopic trends for heavy nuclei is important for further study of the capabilities of this promising version of the shell-model calculations.

Dipole magnetic moments↗
DOE OSTI2021

Quantitative comparison of opacities calculated using the R -matrix and distorted-wave methods: Fe xvii

ABSTRACT We present here a detailed calculation of opacities for Fe xvii at the physical conditions corresponding to the base of the Solar convection zone. Many ingredients are involved in the calculation of opacities. We review the impact of each ingredient on the final monochromatic and mean opacities (Rosseland and Planck). The necessary atomic data were calculated with the R-matrix and the distorted-wave (DW) methods. We study the effect of broadening, of resolution, of the extent of configuration sets and of configuration interaction to understand the differences between several theoretical predictions as well as the existing large disagreement with measurements. New Dirac R-matrix calculations including all configurations up to the n = 4, 5, and 6 complexes have been performed as well as corresponding Breit–Pauli DW calculations. The DW calculations have been extended to include autoionizing initial levels. A quantitative contrast is made between comparable DW and R-matrix models. We have reached self-convergence with n = 6 R-matrix and DW calculations. Populations in autoionizing initial levels contribute significantly to the opacities and should not be neglected. The R-matrix and DW results are consistent under the similar treatment of resonance broadening. The comparison with the experiment shows a persistent difference in the continuum while the filling of the windows shows some improvement. This study defines our path to the next generation of opacities and opacity tables for stellar modelling.

Delahaye, F.↗
DOE OSTI2023

Solving the nuclear pairing model with neural network quantum states

In this work, we present a variational Monte Carlo method that solves the nuclear many-body problem in the occupation number formalism exploiting an artificial neural network representation of the groundstate wave function. A memory-efficient version of the stochastic reconfiguration algorithm is developed to train the network by minimizing the expectation value of the Hamiltonian. We benchmark this approach against widely used nuclear many-body methods by solving a model used to describe pairing in nuclei for different types of interaction and different values of the interaction strength. Despite its polynomial computational cost, our method outperforms coupled-cluster and provides energies that are in excellent agreement with the numerically-exact full configuration interaction values.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Definitive Assessment of the Accuracy, Variationality, and Convergence of Relativistic Coupled Cluster and Density Matrix Renormalization Group in 100-Orbital Space

Accuracy, variationality, and convergence underpin the reliability of modern electronic structure methods, yet definitive benchmarks in the relativistic regime remain elusive due to the absence of numerically exact full configuration interaction (CI) references. Recent algorithmic advances in the CI framework, enabled by the small-tensor-product (STP) decomposition approach, have dramatically extended the tractable size of the configuration space, making numerically exact CI calculations feasible in large active spaces previously beyond reach. In this paper, we employ the recently developed STP-CI framework to perform large-scale numerically exact CI calculations and directly benchmark relativistic coupled cluster and density matrix renormalization group methods. Definitive benchmarking of approximate relativistic electronic structure methods is ensured through the application of the gap theorem, which provides rigorous error bounds on the CI reference and establishes a controlled standard for assessing accuracy, variationality, and convergence.

Chemical calculations↗
DOE OSTI2024

WFOT: A Wave Function Overlap Tool between Single- and Multi-Reference Electronic Structure Methods for Spectroscopy Simulation

We report the development of a novel diagnostic tool, named wave function overlap tool (WFOT), designed to evaluate the overlap between wave functions computed at single-reference [i.e., time-dependent density functional theory or configuration interaction singles (CIS)] and multireference (i.e., CASSCF/CASPT2) electronic structure levels of theory. It relies on truncating the single- and multireference WFs to CIS-like expansions spanning the same configurational space and maximizing the molecular orbital overlap by means of a unitary transformation. To demonstrate the functionality of the tool, we calculate the transient spectrum of acetylacetone by evaluating excited state absorption signals with multireference quality on top of single-reference on-the-fly dynamics simulations. Semiautomatic spectra generation is facilitated by interfacing the tool with the COBRAMM package, which also allows one to use WFOT with several quantum chemistry codes such as Gaussian, NWChem, and OpenMolcas. Finally, other exciting possibilities for the utilization of the code beyond the simulation of transient absorption spectroscopy are eventually discussed.

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

Unusual behavior of Cooper minima of n s subshells in high- Z atoms

Here, a study of Cooper minima (CM) arising from the photoionization of 6$\textit{s}$, 5$\textit{s}$, and 4$\textit{s}$ subshells of high-$\textit{Z}$ atoms has been performed using Dirac-Fock (DF), two-channel relativistic-random-phase approximation (RRPA), and fully coupled RRPA. The results show huge splittings between $ns → εp_{3/2}$ and $ns → εp_{1/2}$ CM which increase with $\textit{Z}$ owing primarily to the relativistic interactions (spin orbit) that are attractive for the $εp_{1/2}$ final state but repulsive for the corresponding $εp_{3/2}$. In addition, it was found that correlation in the form of interchannel coupling (essentially configuration interaction in the final continuum states) plays a huge role in determining the location of the CM. For 6$\textit{s}$ photoionization, the 6$s → εp_{3/2}$ and 6$s → εp_{1/2}$ CM behave completely different as a function of $\textit{Z}$ . It was also found that for 5$\textit{s}$ and 4$\textit{s}$ photoionization, the CM move below the threshold, with increasing $\textit{Z}$ , and, at high enough $\textit{Z}$ , the 5s → εp3/2 and 4s → εp3/2 CM re-emerge into the continuum. The calculations have been carried out for the $\textit{ns}$ subshells of Hg ($\textit{Z}$ = 80), Rn ($\textit{Z}$ = 86), Ra ($\textit{Z}$ = 88), No ($\textit{Z}$ = 102), Cn ($\textit{Z}$ = 112), and Og ($\textit{Z}$ = 118).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Diabatic Valence-Hole States in the C 2 Molecule: “Putting Humpty Dumpty Together Again”

Despite the long history of spectroscopic studies of the C$_2$ molecule, fundamental questions about its chemical bonding are still being hotly debated. The complex electronic structure of C$_2$ is a consequence of its dense manifold of near-degenerate, low-lying electronic states. A global multi-state diabatic model is proposed here to disentangle the numerous configuration interactions within four symmetry manifolds of C$_2$ ($^{1}\Pi_g$, $^{3}\Pi_g$, $^{1}\Sigma_u^+$, and $^{3}\Sigma_u^+$). The key concept of our model is the existence of two "valence-hole" configurations, $2\sigma_g^22\sigma_u^11\pi_{u}^33\sigma_g^2$ for $^{1,3}\Pi_g$ states and $2\sigma_g^22\sigma_u^11\pi_{u}^43\sigma_g^1$ for $^{1,3}\Sigma_u^+$ states that derive from $3\sigma_g\leftarrow2\sigma_u$ electron promotion. The lowest-energy state from each of the four C$_2$ symmetry species is dominated by this type of valence-hole configuration at its equilibrium internuclear separation. Finally, as a result of their large binding energy (nominal bond order of 3) and correlation with the 2s$^2$2p$^2$+2s2p$^3$ separated-atom configurations, the presence of these valence-hole configurations has a profound impact on the $global$ electronic structure and unimolecular dynamics of C$_2$.

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

Modeling Strong Light-Matter Coupling in Correlated Systems: State-Averaged Cavity Quantum Electrodynamics Complete Active Space Self-Consistent Field Theory

The description of strongly correlated systems interacting with quantized cavity modes poses significant theoretical challenges due to the combinatorial scaling of electronic and photonic degrees of freedom. Recent advances addressing this complexity include cavity quantum electrodynamics (QED) generalizations of complete active space configuration interaction and density matrix renormalization group methods. In this work, we introduce a QED extension of state-averaged complete active space self-consistent field theory, which incorporates cavity-induced correlations through a second-order orbital optimization framework with robust convergence properties. The method is implemented using both photon number state and coherent state representations, with the latter showing robust origin invariance in the energies regardless of the completeness of the photonic Fock space. The implementation enables symmetry-free orbital relaxations to account for photon-mediated symmetry breaking in polaritonic systems. Numerical validation on lithium hydride, hydroxide anion, and magnesium hydride cation demonstrates that this method achieves significantly improved accuracy in modeling ground-state and polariton potential energy surfaces compared to QED-CASCI in a fixed orbital basis. In these studies, we reach sub-kcal/mol accuracy in potential energy surface in much smaller active spaces than are required for QED-CASCI. This advancement provides a more robust approach for studying cavity-altered chemical landscapes for ground and exited strongly coupled systems.

CASSCF↗
DOE OSTI2024

Barrier penetration in a discrete-basis formalism

The dynamics of a many-particle system are often modeled by mapping the Hamiltonian into a Schrödinger equation. An alternative approach is to solve the Hamiltonian equations directly in a model space of many-body configurations. In a previous paper the numerical convergence of the two approaches was compared with a simplified treatment of the Hamiltonian representation. Here, we extend the comparison to the nonorthogonal model spaces that would be obtained by the generator-coordinate method. With a suitable choice of the collective-variable grid, a configuration-interaction Hamiltonian can reproduce the Schrödinger dynamics very well. Furthermore, the method as implemented here requires that the barrier height is not much larger than the zero-point energy in the collective coordinates of the configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Accelerating eigenvalue computation for nuclear structure calculations via perturbative corrections

Subspace projection methods utilizing perturbative corrections have been proposed for computing the lowest few eigenvalues and corresponding eigenvectors of large Hamiltonian matrices. In this paper, we build upon these methods and introduce the term Subspace Projection with Perturbative Corrections (SPPC) method to refer to this approach. We tailor the SPPC for nuclear many-body Hamiltonians represented in a truncated configuration interaction subspace, i.e., the no-core shell model (NCSM). We use the hierarchical structure of the NCSM Hamiltonian to partition the Hamiltonian as the sum of two matrices. The first matrix corresponds to the Hamiltonian represented in a small configuration space, whereas the second is viewed as the perturbation to the first matrix. Eigenvalues and eigenvectors of the first matrix can be computed efficiently. Because of the split, perturbative corrections to the eigenvectors of the first matrix can be obtained efficiently from the solutions of a sequence of linear systems of equations defined in the small configuration space. These correction vectors can be combined with the approximate eigenvectors of the first matrix to construct a subspace from which more accurate approximations of the desired eigenpairs can be obtained. We show by numerical examples that the SPPC method can be more efficient than conventional iterative methods for solving large-scale eigenvalue problems such as the Lanczos, block Lanczos and the locally optimal block preconditioned conjugate gradient (LOBPCG) method. The method can also be combined with other methods to avoid convergence stagnation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Evaluation of the excitation spectra with diffusion Monte Carlo on an auxiliary bosonic ground state

We aim to improve upon the variational Monte Carlo (VMC) approach for excitations replacing the Jastrow factor by an auxiliary bosonic (AB) ground state and multiplying it by a fermionic component factor. The instantaneous change in imaginary time of an arbitrary excitation in the original interacting fermionic system is obtained by measuring observables via the ground-state distribution of walkers of an AB system that is subject to an auxiliary effective potential. The effective potential is used to (i) drive the AB system’s ground-state configuration space toward the configuration space of the excitations of the original fermionic system and (ii) subtract from a diffusion Monte Carlo (DMC) calculation contributions that can be included in conventional approximations, such as mean-field and configuration interaction (CI) methods. In this novel approach, the AB ground state is treated statistically in DMC, whereas the fermionic component of the original system is expanded in a basis. The excitation energies of the fermionic eigenstates are obtained by sampling a fermion–boson coupling term on the AB ground state. We show that this approach can take advantage of and correct for approximate eigenstates obtained via mean-field calculations or truncated interactions. We demonstrate that the AB ground-state factor incorporates the correlations missed by standard Jastrow factors, further reducing basis truncation errors. Relevant parts of the theory have been tested in soluble model systems and exhibit excellent agreement with exact analytical data and CI and VMC approaches. In particular, for limited basis set expansions and sufficient statistics, AB approaches outperform CI and VMC in terms of basis size for the same systems. The implementation of this method in current codes, despite being demanding, will be facilitated by reusing procedures already developed for calculating ground-state properties with DMC and excitations with VMC.

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

Origin of Ferromagnetic Exchange Coupling in Donor–Acceptor Biradical Analogues of Charge-Separated Excited States

A new donor–acceptor biradical complex, Tp Cum,Me Zn(SQ-VD) (Tp Cum,Me Zn + = zinc(II) hydro-tris(3-cumenyl-5-methylpyrazolyl)borate complex cation; SQ = orthosemiquinone; VD = oxoverdazyl), which is a ground-state analogue of a charge-separated excited state, has been synthesized and structurally characterized. Here, the magnetic exchange interaction between the S = 1/2 SQ and the S = 1/2 VD within the SQ-VD biradical ligand is observed to be ferromagnetic, with $J_{SQ-VD}$ = +77 cm –1 (H = -2$J_{SQ-VD}\bar{S}_{SQ}·\bar{S}_{VD}$) determined from an analysis of the variable-temperature magnetic susceptibility data. The pairwise biradical exchange interaction in Tp Cum,Me Zn(SQ-VD) can be compared with that of the related donor–acceptor biradical complex Tp Cum,Me Zn(SQ-NN) (NN = nitronyl nitroxide, S = 1/2), where $J_{SQ-NN}$ ≅ +550 cm –1 . This represents a dramatic reduction in the biradical exchange by a factor of ~7, despite the isolobal nature of the VD and NN acceptor radical SOMOs. Computations assessing the magnitude of the exchange were performed using a broken-symmetry density functional theory (DFT) approach. These computations are in good agreement with those computed at the CASSCF NEVPT2 level, which also reveals an S = 1 triplet ground state as observed in the magnetic susceptibility measurements. A combination of electronic absorption spectroscopy and CASSCF computations has been used to elucidate the electronic origin of the large difference in the magnitude of the biradical exchange coupling between Tp Cum,Me Zn(SQ-VD) and Tp Cum,Me Zn(SQ-NN). A Valence Bond Configuration Interaction (VBCI) model was previously employed to highlight the importance of mixing an SQ SOMO → NNLUMO charge transfer configuration into the electronic ground state to facilitate the stabilization of the high-spin triplet (S = 1) ground state in Tp Cum,Me Zn(SQ-NN). Here, CASSCF computations confirm the importance of mixing the pendant radical (e.g., VD, NN) LUMO (VD LUMO and NN LUMO ) with the SOMO of the SQ radical (SQ SOMO ) for stabilizing the triplet, in addition to spin polarization and charge transfer contributions to the exchange. An important electronic structure difference between Tp Cum,Me Zn(SQ-VD) and Tp Cum,Me Zn(SQ-NN), which leads to their different exchange couplings, is the reduced admixture of excited states that promote ferromagnetic exchange into the Tp Cum,Me Zn(SQ-VD) ground state, and the intrinsically weaker mixing between the VD LUMO and the SQ SOMO compared to that observed for Tp Cum,Me Zn(SQ-NN), where this orbital mixing is significant. The results of this comparative study contribute to a greater understanding of biradical exchange interactions, which are important to our understanding of excited-state singlet–triplet energy gaps, electron delocalization, and the generation of electron spin polarization in both the ground and excited states of (bpy)Pt(CAT-radical) complexes.

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

Atomistic and phase field simulations of three dimensional interactions of {$10\bar12$} twins with grain boundaries in Mg: twin transmission and dislocation emission

In polycrystals, the interaction of dislocations and twins with grain boundaries (GBs) plays a role in hardening and formability during plastic deformation. While dislocation-GB interactions are relatively well-understood, twin-GB interactions remain mostly unknown. In this work, an approach using molecular dynamics and phase-field simulations is followed to study the forward and lateral interactions between {$10\bar12$} twins and tilt grain boundaries in Mg. Molecular dynamics results show that the resolved shear stress on slip/twinning modes of the neighboring grain, not the geometric alignment, is the dominant factor in determining the outcome of the twin-GB interactions. For some lateral interaction configurations, as the misorientation angle increases, the resolved shear stress on the same {$10\bar12$} twin variant of the neighboring grain reduces while it increases for slip or I 2 stacking fault emissions or other twin modes such as {$11\bar12$} and {$10\bar11$}, explaining why twin transmission is not seen at high misorientation angles. Furthermore, lateral and forward interactions of the twin with tilt grain boundaries whose misorientation axes are normal to the coherent twin boundary show significantly different outcomes. For the forward interaction, the twin is absorbed and stacking faults are emitted when interacting for low misorientation angles (up to 30°) while the lateral interaction results in twin transmission, nucleation of a {$11\bar12$} twin, and emission of I 2 stacking faults. Finally, comparisons between twin interactions with symmetric and asymmetric tilt GBs with different GB structures show similar outcomes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

I ntera C hem : Exploring Excited States in Virtual Reality with Ab Initio Interactive Molecular Dynamics

InteraChem is an ab initio interactive molecular dynamics (AI-IMD) visualizer that leverages recent advances in virtual reality hardware and software, as well as the graphical processing unit (GPU)-accelerated TeraChem electronic structure package, in order to render quantum chemistry in real time. We introduce the exploration of electronically excited states via AI-IMD using the floating occupation molecular orbital-complete active space configuration interaction method. The optimization tools in InteraChem enable identification of excited state minima as well as minimum energy conical intersections for further characterization of excited state chemistry in small- to medium-sized systems. We demonstrate that finite-temperature Hartree–Fock theory is an efficient method to perform ground state AI-IMD. InteraChem allows users to track electronic properties such as molecular orbitals and bond order in real time, resulting in an interactive visualization tool that aids in the interpretation of excited state chemistry data and makes quantum chemistry more accessible for both research and educational purposes.

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

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear–electronic coupling in addition to accurately predicting Franck–Condon vertical excitations. Here, in this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗