Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CASSCF”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

On the electron affinity of the oxygen atom

The electron affinity (EA) of oxygen is computed to be 1.287 eV, using 2p electron full configuration-interaction (CI) wave functions expanded in a 6s5p3d2f Slater-type orbital basis. The best complete active space self-consistent field - multireference CI (CASSCF-MRCI) result including only 2p correlation is 1.263 eV. However, inclusion of 2s intrashell and 2s2p intershell correlation increases the computed EA to 1.290 at the CASSCF-MRCI level. At the full CI basis set limit, the 2s contribution to the electron affinity is estimated to be as large as 0.1 eV. This study clearly establishes the synergistic effect between the higher excitations and basis set completeness on the electron affinity when the 2s electrons are correlated.

Bauschlicher, C. W., Jr.↗
NASA NTRS1986

On the electron affinity of the oxygen atom

The electron affinity of oxygen is computed to be 1.287 eV, at the full CI level using a 6s5p3d 2f Slater-type orbital basis and correlating only the 2p electrons. The best CASSCF-MRCI result including only 2p correlation is 1.263 eV. However, inclusion of 2s intrashell and 2s2p intershell coorelation increases the computed EA to 1.290 eV at the CASSCF-MRCI level. At the full CI basis set limit, the 2s contribution to the electron affinity is estimated to be as large as 0.1 eV. The higher excitation contribution to the electron affinity is found to increase substantially with basis set completeness, especially when the 2s electrons are correlated. Relativistic effects are shown to make a small (less than 0.01 eV) change in the EA.

Bauschlicher, C. W., Jr.↗
NASA NTRS1987

Full CI benchmark calculations for several states of the same symmetry

Full CI (FCI) wave functions are used to compute energies for several electronic states of the same symmetry for SiH2, CH2, and CH2(+). It is found that CASSCF/multireference CI wave functions yield results very similar to FCI, irrespective of whether the CASSCF MOs are optimized independently for each state or using an average of the CASSCF energies for all desired states. The ionization potentials and excitation energies obtained from the FCI calculations should help calibrate methods (such as Green's function approaches, equations of motion and propagator methods, and cluster expansions) in which energy differences are computed directly.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1987

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.↗
NASA NTRS1985

Theoretical studies of transition metal dimers

The CASSCF approach was used to perform the MCSCF calculations for a number of transition metal dimers, including the Sc2, Ti2, Cr2, Cu2, TiV, Y2, Nb2, and Mo2 molecules; in addition, CASSCF/CI calculations were carried out for Sc2, Ti2, Cu2, and Y2. The CASSCF procedure is shown to provide a consistent set of calculations for these molecules, from which trends and a simple qualitative picture of the electronic structure may be derived. In particular, the calculations confirmed the ground states of the Sc2 and the TiV, and led to predictions for other molecules in this series. In addition to specific predictions, the study provides a simple qualitative picture of the bonding in these dimers.

Walch, Stephen P.↗
NASA NTRS1991

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.↗
NASA NTRS1989

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

The bond dissociation energies (D sub e) of C2H2, C2H4, C2H6, N2, N2H2, and N2H4 are studied at various levels of correlation treatment. The convergence of D sub e 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) D sub e 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. It was found that 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. To achieve chemical accuracy (1 kcal/mole) for the D sub e values of the doubly bonded species C2H4 and N2H2 requires one particle basis sets including up through h angular momentum functions (l = 5) and a multireference treatment of electron correlation: still higher levels of calculation are required to achieve chemical accuracy for the triply bonded species C2H2 and N2.

Langhoff, Stephen R.↗
NASA NTRS1994

Theoretical Study of the Nitric Oxide epsilon and 1100 A Bands

Potential-energy curves, transition-moment functions, Einstein Alpha coefficients and radiative lifetimes associated with the Rydberg state of nitric oxide are calculated in large Gaussian basis sets employing internally contracted MRCI from state-averaged CASSCF wave functions. Computationally, a strong sensitivity of the results to the choice of the CASSCF active space is seen, and careful examination of this behavior gives insight into matters related to the treatment of weakly occupied orbitals, which is a consideration in any study Rydberg excited states. The epsilon system transition-moment function and radiative lifetimes are in disagreement with another calculation in the literature, but reasons for the discrepancies are examined, and the present results are determined to be more accurate.

Sheehy, Jeffrey A.↗
NASA NTRS1995

The MCH2+ Systems: Do ScCH2+ and TiCH2+ Have C(sub s) or C(sub 2v) Symmetry and A Comparison of the B3LYP Method to Other Approaches

ScCH2(+) and TiCH2(+) are found to have C(sub s) symmetry, while the remaining systems have the expected C(sub 2v) symmetry. The C(sub s) symmetry structure is favored for ScCH2(+) and TiCH2(+) because it allows donation from one of the CH bonds into an empty 3d orbital. This distortion stabilizes the systems by less than 2 kcal/mol at the B3LYP, CASSCF, and CCSD(T) levels of theory. The SCF and B3LYP approaches find the C2(sub v) structure to be a saddle point, while the CASSCF approach finds the C(sub 2v) structure to be a stationary point. For V-Cu there are no empty 3d orbitals of the correct symmetry, and therefore the distortion to C(sub s), is unfavorable. The B3LYP binding energies are in good agreement with experiment and our previous best estimates.

Ricca, Alessandra↗
DOE OSTI2025

Perspective on Many-Body Methods for Molecular Polaritonic Systems

Recent advances in strong light–matter interactions have revealed a wealth of new physical phenomena in molecules embedded in optical cavities, including modified chemical reactivity, altered excitation spectra, and novel quantum correlations. To describe these effects from first-principles, the field of ab initio quantum electrodynamics (QED) has emerged as a compelling extension of quantum chemistry that treats electronic and photonic degrees of freedom on equal footing. In this Perspective, we review the growing landscape of many-body QED methods, including Hartree–Fock, density functional theory (QEDFT), time-dependent DFT (QED-TDDFT), configuration interaction (QED-CI), complete active space (QED-CASSCF), coupled cluster (QED-CC), quantum Monte Carlo (QED-QMC), and density matrix renormalization group (QED-DMRG), highlighting recent developments and implementations. We further explore real-time methods, gradient and Hessian formalisms, and the integration of nonadiabatic nuclear dynamics. Applications range from benchmark simulations of polaritonic chemistry to quantum simulations on emerging quantum hardware. We conclude by outlining future directions for theory development and interdisciplinary efforts at the interface of quantum chemistry, condensed matter, and quantum optics.

36 MATERIALS SCIENCE↗
DOE OSTI2026

NOCI-F Electronic Couplings in Assemblies of Indolonaphthyridine Molecules: From Dimers to the Full Stack

Key electronic processes related to molecular excitonic states of finite stacks of indolonaphthyridine molecules are analyzed via the non-orthogonal configuration interaction with fragments (NOCI-F) method. Indolonaphthyridine is an organic chromophore that can undergo several electronic photoexcitation-related intermolecular processes, such as exciton and electron transfer. The structures studied here are noncrystalline arrangements built as either ordered stacks of indolonaphthyridine or stacks extracted from molecular dynamics simulations including thermal disorder. Taking dimers or trimers from either model, we performed CASSCF and NOCI-F calculations to quantify the intermolecular electronic couplings governing singlet fission, excited singlet and triplet diffusion, and hole and electron diffusion processes. Also, comparing the results for the different models, we studied the effect of structural disorder and distortion on these couplings. Finally, we present a newly developed, advanced postanalysis tool. It takes the NOCI-F data as input to carry out a multifragment full Hamiltonian procedure that involves the complete stack, providing physical information not available from the dimer/trimer models, hence giving access to additional insight into the material’s properties.

coupling reactions↗
DOE OSTI2025

Nonunitary Variational Quantum Eigensolver with the Localized Active Space Method and Cost Mitigation

Accurately describing strongly correlated systems with affordable quantum resources remains a central challenge for quantum chemistry applications on near and intermediate term quantum computers. The localized active space self-consistent field (LASSCF) approximates the complete active space self-consistent field (CASSCF) by generating active space-based wave functions within specific fragments while treating interfragment correlation with mean-field approach, hence is computationally less expensive. Hardware-efficient ansatzes (HEA) offer affordable and shallower circuits, yet they often fail to capture the necessary correlation. Previously, Jastrow-factor-inspired nonunitary qubit operators were proposed to use with HEA for variational quantum eigensolver (VQE) calculations (so-called nuVQE), as they do not increase circuit depths and recover correlation beyond the mean-field level for Hartree–Fock initial states. Here, in this study, we explore running nuVQE with LASSCF as the initial state. The method, named LAS-nuVQE, is shown to recover interfragment correlations, reach chemical accuracy with a small number of gates (<70) in both H 4 and square cyclobutadiene (C 4 H 4 ), and produces more accurate energetics than its HEA counterparts at all circuit depths. To further address the inherent symmetry-breaking in HEA, we implemented spin-constrained LAS-nuVQE to extend the capabilities of HEA further and show spin-pure results for square cyclobutadiene. We also mitigate the increased measurement overhead of nuVQE via Pauli grouping and shot-frugal sampling, reducing measurement costs by up to 2 orders of magnitude compared to ungrouped operator, and show that one can achieve better accuracy with a small number of shots (10 3–4 ) per one expectation value calculation compared to noiseless simulations with one or two orders of magnitude more shots. Finally, wall clock time estimates show that, with our measurement mitigation protocols, nuVQE becomes a cheaper and more accurate alternative than vanilla VQE with HEA. Taken together, these developments illustrate a practical pathway toward performing multireference chemical simulations with accuracy and affordable resources on today’s quantum hardware, achieving both accuracy and affordability in challenging correlated systems.

Wang, Qiaohong [Univ. of Chicago, IL (United State↗
DOE OSTI2024

Exploring the Computational Aspects of Propylene Oligomerization Catalysis Using M′ 2 M Type Trimetallic MOF Nodes

Metal-organic frameworks have emerged as promising materials in the field of catalysis. They offer an optimal ground for screening catalysts and tailoring their catalytic properties. Here, in this work, via density functional theory (DFT) calculations, we investigated the catalytic activity of the trimetallic MOF nodes, M' 2 M for propylene oligomerization, by varying the active metal M from Sc to Cu with M' being Fe, aiming to grasp the impact of altering the active atoms on the catalyst's activity. Additionally, we examined how substituting the spectator atom, M', with other transition metals, i.e., from Sc to Cu, affects these energy barriers, keeping Ni as the active metal. We proposed several cases with lower or comparable energy barriers to the experimentally reported Fe 2 Ni trimetallic MOF node. In addition, we found a correlative relationship between spin-density from natural population analysis and energy barriers in the realm of C-C bond formation, whereby an elevation in spin-density is found to be inversely proportional to the magnitude of the energy barriers. Moreover, we calculated the energy barriers for C-C coupling and beta-hydride elimination using multireference NEVPT2 calculations on top of the CASSCF wave function to validate the rate-determining step that is predicted by DFT.

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

Conical Intersection Accessibility Dictates Brightness in Red Fluorescent Proteins

Red fluorescent protein (RFP) variants are highly sought after for in-vivo imaging since longer wavelengths improve depth and contrast in fluorescence imaging. However, the lower energy emission wavelength usually correlates with a lower fluorescent quantum yield compared to their green emitting counterparts. To guide the rational design of bright variants, we have theoretically assessed two variants (mScarlet and mRouge) which are reported to have very different brightness. Using an α-CASSCF QM/MM framework (chromophore and all protein residues within 6 Å of it in the QM region, for a total of more than 450 QM atoms), we identify key points on the ground and first excited state potential energy surfaces. The brighter variant mScarlet has a rigid scaffold, and the chromophore stays largely planar on the ground state. The dimmer variant mRouge shows more flexibility and can accommodate a pre-twisted chromophore conformation which provides easier access to conical intersections. Notably, the main difference between the variants lies in the intersection seam regions, which appear largely inaccessible in mScarlet but partially accessible in mRouge. This observation is mainly related with changes in the cavity charge distribution, the hydrogen-bonding network involving the chromophore and a key ARG/THR mutation (which changes both charge and steric hindrance).

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

Using a multistate mapping approach to surface hopping to predict the ultrafast electron diffraction signal of gas-phase cyclobutanone

Using the recently developed multistate mapping approach to surface hopping (multistate MASH) method combined with SA(3)-CASSCF(12,12)/aug-cc-pVDZ electronic structure calculations, the gas-phase isotropic ultrafast electron diffraction (UED) of cyclobutanone is predicted and analyzed. After excitation into the n-3s Rydberg state (S2), cyclobutanone can relax through two S2/S1 conical intersections, one characterized by compression of the CO bond and the other by dissociation of the α–CC bond. Subsequent transfer into the ground state (S0) is then achieved via two additional S1/S0 conical intersections that lead to three reaction pathways: α ring-opening, ethene/ketene production, and CO liberation. The isotropic gas-phase UED signal is predicted from the multistate MASH simulations, allowing for a direct comparison to the experimental data. This work, which is a contribution to the cyclobutanone prediction challenge, facilitates the identification of the main photoproducts in the UED signal and thereby emphasizes the importance of dynamics simulations for the interpretation of ultrafast experiments.

Chemistry↗
DOE OSTI2025

Excited electronic states of Na 2 and K 2 : The potential for long-lived “reservoir” states leading to collision induced population inversions

Potential energy curves (PECs) for the spin-free (ΛS) and spin–orbit (Ω) states associated with the four lowest-lying dissociation channels of Na 2 and K 2 were calculated at the SA-CASSCF/SO-CASPT2/aug-cc-pwCVQZ-DK level. The PECs of Na 2 were consistent with the experimental data and with the FS-CCSD (2,0) calculations, reproducing the double-well and the “shelf” character for some of the potentials of the excited states. For K 2 , the PECs behaved in a similar way and the spectroscopic parameters for the ground and the excited states are in good agreement with the available experimental values. The dissociation energy of K 2 was predicted to be D e = 4454 cm −1 , within an agreement of 5 cm −1 with the experiments. For Na 2 , D e = 5789 cm −1 compared to the experimental value of 6022 cm −1 . The inclusion of spin–orbit coupling effects resulted in avoided crossings, which affect the PECs. Spin–orbit changes the predicted curves for some excited Ω states arising from ΛS states that overlap each other, affecting their associated vibrational frequencies and bond distances. Here, the current studies of the low-lying states in K 2 reveal a similar structure to those of Na 2 , which suggests the accessibility of long-lived energy storing reservoir states and possible population inversions in K 2 following prior experimental work on the reaction of halogen atoms with Na 3 to produce excited states of Na 2 .

Ab-initio methods↗
DOE OSTI2025

Fast methods for multisite charge transfer. Processes II. Analytic nuclear gradients and nonadiabatic dynamics for cCASSCF(1,n) and cCASSCF(2n-1,n) wavefunctions

In this work we derive and implement analytic nuclear gradients and derivative couplings for a constrained complete active space self-consistent field with a small active space designed to model electron or hole transfer. Using a Lagrangian formalism, we are able to differentiate both the CASSCF energy and the constraint (which is required for smooth surfaces over a wide range of parameter space), and the resulting efficient algorithm can be immediately applied to nonadiabatic dynamics simulations of charge transfer processes. Here, we run initial surface-hopping simulations of a proton coupled electron transfer event for a phenoxyl–phenol system.

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

Many roads to the seam: How conformational flexibility drives nonadiabatic relaxation in a prototypical tetrapyrrolic chromophore

Large and structurally flexible chromophores pose challenges for in silico modeling of photodeactivation due to the many vibrational modes that can funnel the system toward energy degeneracy. In this work, we examine how the multiple degrees of freedom in biliverdin, a prototypical tetrapyrrolic chromophore, cooperate to drive access to the S 1 /S 0 intersection seam in vacuo. We begin by mapping the ground-state potential energy surface to identify representative biliverdin conformers relevant to photoexcitation. We then use a CASSCF-based framework to map the excited-state landscape and characterize the intersection seam, identifying distinct conical-intersection types. Finally, we employ ab initio multiple spawning to resolve the dynamical pathways by which the system accesses these regions. DFT potential-energy and free-energy mappings indicate that, although several conformers are relevant, the “locked-helix” ZsZsZs conformer predominates in the ground state. The intersection seam comprises numerous geometrically distinct regions characterized by varying degrees and combinations of dihedral torsion, pyramidalization, and bond-length alternation. Yet only select regions lie within energetic reach, and moderate barriers separate them from the S 1 minimum. Nonadiabatic dynamics combined with multivariate analyses show that, despite extensive mode coupling during deactivation that guides the system toward multiple regions of the seam, a single dihedral torsion, together with bond-length alternation, predominantly drives energy degeneracy. This work offers new insight into biliverdin’s intrinsic photochemical response and underscores a general feature of flexible chromophores: many modes may participate during photorelaxation, but only a limited subset ultimately dictates seam accessibility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗