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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Efficient Pathways of Excitation Energy Transfer from Delocalized S 2 Excitons in the Peridinin–Chlorophyll a –Protein Complex

Excitation energy transfer (EET) in peridinin–chlorophyll–protein (PCP) complexes is dominated by the S 1 → Q y pathway, but the high efficiencies cannot be solely explained by this one pathway. We postulate that EET from peridinin S 2 excitons may also be important. We use complete active space configuration interaction calculations and the transition density cube method to calculate Coulombic couplings between peridinin and chlorophyll a in the PCP complex of the dinoflagellate Amphidinium carterae and compare monomeric and dimeric delocalized peridinin S 2 excited states. Our calculations show that the S 2 → Q y EET pathway from peridinin to chlorophyll a is the dominant energy transfer pathway from the S 2 excited state in PCP, with several values in the sub-picosecond range. Here, this result suggests that the S 2 → Q y EET pathway may be responsible for the reported chlorophyll a bleaching signature seen in experiment at around 200 fs, and not the S 2 → Q x pathway as previously suggested.

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

Contextual subspace variational quantum eigensolver calculation of the dissociation curve of molecular nitrogen on a superconducting quantum computer

Abstract We present an experimental demonstration of the Contextual Subspace Variational Quantum Eigensolver on superconducting hardware. Calculating the potential energy curve of molecular nitrogen proves challenging for many conventional quantum chemistry techniques, since static correlation dominates in the dissociation limit. Our quantum simulations retain good agreement with the Full Configuration Interaction energy, outperforming all benchmarked single-reference wavefunction techniques in capturing the bond-breaking appropriately. Moreover, our methodology is competitive with multiconfigurational approaches but at a saving of quantum resource, meaning larger active spaces can be treated for a fixed qubit allowance. To achieve this result, we deploy an error mitigation/suppression strategy comprised of Dynamical Decoupling, Measurement-Error Mitigation and Zero-Noise Extrapolation. Circuit parallelization also provides passive noise-averaging and improves the effective shot yield to reduce the measurement overhead. Furthermore, we introduce a modified adaptive ansatz construction algorithm that incorporates hardware awareness into our variational circuits, minimizing the transpilation cost for the target qubit topology.

Physics↗
DOE OSTI2025

Dissipative ground state preparation in ab initio electronic structure theory

Dissipative engineering is a powerful tool for quantum state preparation, and has drawn significant attention in quantum algorithms and quantum many-body physics in recent years. In this work, we introduce a novel approach using the Lindblad dynamics to efficiently prepare the ground state for general ab initio electronic structure problems on quantum computers, without variational parameters. These problems often involve Hamiltonians that lack geometric locality or sparsity structures, which we address by proposing two generic types of jump operators for the Lindblad dynamics. Type-I jump operators break the particle number symmetry and should be simulated in the Fock space. Type-II jump operators preserves the particle number symmetry and can be simulated more efficiently in the full configuration interaction space. For both types of jump operators, we prove that in a simplified Hartree-Fock framework, the spectral gap of our Lindbladian is lower bounded by a universal constant. For physical observables such as energy and reduced density matrices, the convergence rate of our Lindblad dynamics with Type-I jump operators remains universal, while the convergence rate with Type-II jump operators only depends on coarse grained information such as the number of orbitals and the number of electrons. To validate our approach, we employ a Monte Carlo trajectory-based algorithm for simulating the Lindblad dynamics for full ab initio Hamiltonians, demonstrating its effectiveness on molecular systems amenable to exact wavefunction treatment.

Quantum chemistry↗
DOE OSTI2023

The role of electron correlations in the electronic structure of putative Chern magnet TbMn6Sn6

Abstract A member of the RMn 6 Sn 6 rare-earth family materials, TbMn 6 Sn 6 , recently showed experimental signatures of the realization of a quantum-limit Chern magnet. In this work, we use quantum Monte Carlo (QMC) and density functional theory with Hubbard U (DFT + U ) calculations to examine the electronic structure of TbMn 6 Sn 6 . To do so, we optimize accurate, correlation-consistent pseudopotentials for Tb and Sn using coupled-cluster and configuration–interaction (CI) methods. We find that DFT + U and single-reference QMC calculations suffer from the same overestimation of the magnetic moments as meta-GGA and hybrid density functional approximations. Our findings point to the need for improved orbitals/wavefunctions for this class of materials, such as natural orbitals from CI, or for the inclusion of multi-reference effects that capture the static correlations for an accurate prediction of magnetic properties. DFT + U with Mn magnetic moments adjusted to the experiment predict the Dirac crossing in bulk to be close to the Fermi level, within ~120 meV, in agreement with the experiments. Our non-stoichiometric slab calculations show that the Dirac crossing approaches even closer to the Fermi level, suggesting the possible realization of Chern magnetism in this limit.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Ab initio pressure-dependent reaction kinetics of methyl propanoate radicals

The unimolecular dissociation and isomerization kinetics of the three methyl propanoate (MP) radicals, CH 3 CH 2 C(=O)O$\dot{C}$H 2 (MP-m), CH 3 $\dot{C}$HC(=O)OCH 3 (MP-α), and $\dot{C}$H 2 CH 2 C(=O)OCH 3 (MP-β), are theoretically investigated using high-level ab initio methods and the Rice–Ramsperger–Kassel–Marcus (RRKM)/master equation (ME) theory. Stationary-point energies are obtained using the coupled cluster singles and doubles with perturbative triples correction (CCSD(T)), multi-reference singles and doubles configuration interaction (MRSDCI) with the Davidson-Silver (DS) correction, and multi-reference averaged coupled pair functional (MRACPF2) theories. The isomerization barriers between the three radicals are predicted to be generally lower than the corresponding bond dissociation channels, leading to a strongly coupled reaction system in subsequent kinetics studies. The phenomenological temperature- and pressure-dependent rate coefficients are computed using the RRKM/ME theory over a temperature range of 500 to 2000 K and at a pressure range of 0.01 atm to the high-pressure limit, which are then fitted to modified Arrhenius expressions. The β-scission rate coefficients of MP-α to CH 3 CHC(=O) and CH 3 $\dot{O}$ are predicted to be the smallest because of its highest activation energy among all studied unimolecular reactions channels. Analysis of branching fractions shows that both MP-m and MP-α radicals mainly decompose to the bimolecular products CH 3 CH 2 $\dot{C}$(=O) and H 2 CO, whereas the MP-β radical primarily decomposes via cleavage of a C–C bond to form C 2 H 4 and CH 3 O$\dot{C}$(=O). The isomerization channels dominate at low temperatures, the branching fractions of which decrease with increasing temperature and become very minor at about 2000 K. Our accurate rate coefficients and branching fractions help to illuminate the unique combustion properties of MP.

09 BIOMASS FUELS↗
DOE OSTI2017

Controlling the excited-state dynamics of low band gap, near-infrared absorbers via proquinoidal unit electronic structural modulation

While the influence of proquinoidal character upon the linear absorption spectrum of low optical bandgap π-conjugated polymers and molecules is well understood, its impact upon excited-state relaxation pathways and dynamics remains obscure. We report the syntheses, electronic structural properties, and excited-state dynamics of a series of model highly conjugated near-infrared (NIR)-absorbing chromophores based on a (porphinato)metal(II)-proquinoidal spacer-(porphinato)metal(II) (PM-Sp-PM) structural motif. A combination of excited-state dynamical studies and time-dependent density functional theory calculations: (i) points to the cardinal role that excited-state configuration interaction (CI) plays in determining the magnitudes of S 1 → S 0 radiative (k r ), S 1 → T 1 intersystem crossing (k ISC ), and S 1 → S 0 internal conversion (k IC ) rate constants in these PM-Sp-PM chromophores, and (ii) suggests that a primary determinant of CI magnitude derives from the energetic alignment of the PM and Sp fragment LUMOs (ΔE L ). These insights not only enable steering of excited-state relaxation dynamics of high oscillator strength NIR absorbers to realize either substantial fluorescence or long-lived triplets (τ$_ {T_1}$ > μs) generated at unit quantum yield (Φ ISC = 100%), but also crafting of those having counter-intuitive properties: for example, while (porphinato)platinum compounds are well known to generate non-emissive triplet states (Φ ISC = 100%) upon optical excitation at ambient temperature, diminishing the extent of excited-state CI in these systems realizes long-wavelength absorbing heavy-metal fluorophores. In conclusion, this work highlights approaches to: (i) modulate low-lying singlet excited-state lifetime over the picosecond-to-nanosecond time domain, (ii) achieve NIR fluorescence with quantum yields up to 25%, (iii) tune the magnitude of S 1 –T 1 ISC rate constant from 10 9 to 10 12 s -1 and (iv) realize T 1 -state lifetimes that range from ~0.1 to several μs, for these model PM-Sp-PM chromophores, and renders new insights to evolve bespoke photophysical properties for low optical bandgap π-conjugated polymers and molecules based on proquinoidal conjugation motifs.

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

An ab initio based full-dimensional potential energy surface for OH + O 2 ⮂ HO 3 and low-lying vibrational levels of HO 3

To offer an in-depth understanding of the HO 3 radical and its dissociation to OH + O 2 , a six-dimensional potential energy surface (PES) has been constructed by fitting 2087 energy points for the electronic ground state of HO 3 (X 2 A") using the permutation invariant polynomial-neural network (PIP-NN) approach. The energy points were calculated using an explicitly-correlated and Davidson-corrected multi-reference configuration interaction approach with the correlation-consistent polarized valence double zeta basis (MRCI( Q )-F12/VDZ-F12). On the PES, the trans -HO 3 isomer is found to be the global minimum, 33.0 cm -1 below the cis -HO 3 conformer, which is consistent with previous high-level theoretical investigations. The dissociation to the OH + O 2 asymptote from both conformers is shown to be barrierless. As a benchmark from a recently developed high-accuracy thermochemistry protocol, D 0 for trans -HO 3 is calculated to be 2.29 ± 0.36 kcal mol -1 , only slightly deeper than the value of 2.08 kcal mol -1 obtained using the PES, and in reasonable agreement with the experimentally estimated value of 2.93 ± 0.07 kcal mol -1 . Using this PES, low-lying vibrational energy levels of HO 3 are determined using an exact quantum Hamiltonian and compared with available experimental results.

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

Beyond the Coulson–Fischer point: characterizing single excitation CI and TDDFT for excited states in single bond dissociations

Linear response time dependent density functional theory (TDDFT), which builds upon configuration interaction singles (CIS) and TD-Hartree–Fock (TDHF), is the most widely used class of excited state quantum chemistry methods and is often employed to study photochemical processes. This paper studies the behavior of the resulting excited state potential energy surfaces beyond the Coulson–Fischer (CF) point in single bond dissociations, when the optimal reference determinant is spin-polarized. Many excited states exhibit sharp kinks at the CF point, and connect to different dissociation limits via a zone of unphysical concave curvature. In particular, the unrestricted M S = 0 lowest triplet T 1 state changes character, and does not dissociate into ground state fragments. The unrestricted M S = ±1 T 1 CIS states better approximate the physical dissociation limit, but their degeneracy is broken beyond the CF point for most single bond dissociations. On the other hand, the M S = ±1 T 1 TDHF states reach the asymptote too soon, by merging with the ground state from the CF point onwards. Use of local exchange–correlation functionals causes M S = ±1 T 1 TDDFT states to resemble their unphysical M S = 0 counterpart. The 2 orbital, 2-electron model system of minimal basis H 2 is analytically treated to understand the origin of these issues, revealing that the lack of double excitations is at the root of these remarkable observations. In conclusion, the behavior of excited state surfaces is also numerically examined for species like H 2 , NH 3 , C 2 H 6 and LiH in extended basis sets.

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

On the nonadiabatic collisional quenching of OH(A) by H 2 : a four coupled quasi-diabatic state description

A four-state diabatic potential energy matrix (DPEM), H d , for the description of the nonadiabatic quenching of OH(A 2 Σ + ) by collisions with H 2 is reported. The DPEM is constructed as a fit to adiabatic energies, energy gradients, and derivative couplings obtained exclusively from multireference configuration interaction wave functions. A four-adiabatic-electronic-state representation is used in order to describe all energetically accessible regions of the nuclear coordinate space. Partial permutation-inversion symmetry is incorporated into the representation. Here, the fit is based on electronic structure data at 42 882 points, described by over 1.6 million least squares equations with a root mean square (mean unsigned) error of 178(83) cm –1 . Comparison of ab initio and H d determined minima, saddle points, and energy minimized points on C 2v , C s , C ∞v , and C 1 (noncoplanar) portions of two conical intersection seams are used to establish the accuracy of the H d .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Dissociative electron attachment to 5-bromo-uracil: non-adiabatic dynamics on complex-valued potential energy surfaces

Electron induced dissociation reactions are relevant to many fields, ranging from prebiotic chemistry to cancer treatments. However, the simulation of dissociation electron attachment (DEA) dynamics is very challenging because the auto-ionization widths of the transient negative ions must be accounted for. We propose an adaptation of the ab initio multiple spawning (AIMS) method for complex-valued potential energy surfaces, along the lines of recent developments based on surface hopping dynamics. Our approach combines models for the energy dependence of the auto-ionization widths, obtained from scattering calculations, with survival probabilities computed for the trajectory basis functions employed in the AIMS dynamics. The method is applied to simulate the DEA dynamics of 5-bromouracil in full dimensionality, i.e., taking all the vibrational modes into consideration. The propagation starts on the $π_{2}^{*}$ resonance state and describes the formation of Br - anions mediated by non-adiabatic couplings. The potential energies, gradients and non-adiabatic couplings were computed with the fractional-occupancy molecular orbital complete-active-space configuration-interaction method, and the calculated DEA cross section are consistent with the observed DEA intensities.

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

Application of saturation absorption spectroscopy to study the hyperfine structure of 235U and accurate 235U/238U isotope ratio determinations at 861.031 nm

The uranium transition 5f 3 6d7s 2 ( 5 L 6 ) → 5f 3 7s 2 7p ( 5 K 5 ) at 861.031 nm was studied using saturation absorption spectroscopy. The hyperfine structure of the 235 U isotope was determined by measuring the Doppler-suppressed absorption spectrum, and the hyperfine structure constants were obtained by employing transition-specific line strength theory and spectral fitting. The magnetic-dipole and electric-quadrupole constants were determined to be A = -3.2887 ± 0.1957 mK and B = 31.8712 ± 8.0789 mK for the upper state 5f 3 7s 2 7p ( 5 K 5 ). A theoretical calculation for the hyperfine structure constants of the 235 U isotope was also carried out by the relativistic configuration-interaction (RCI) method, and the obtained values were found to be in agreement with the measured values. In addition, the line profile of 238 U was experimentally investigated for different excitation intensities and modeled using the RCI method for estimating the transition probability and oscillator strength, with various residual broadening mechanisms considered beyond the radiative decay. Finally, we show that these results allow us to make improvements on the precision of spectral measurements using our previously developed technique for isotope ratio determinations in atomic beams.

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

On the choice of reference orbitals for linear-response calculations of solution-phase K-edge X-ray absorption spectra

The simplest response theory methods for computing vertical excitation spectra in condensed-phase are configuration interaction with single excitations (CIS) and linear-response time-dependent density functional theory (TDDFT) within the Tamm–Dancoff approximation. In applications to X-ray absorption spectroscopy (XAS), methods like CIS and TDDFT that codify only single excitations into the wave function are prone to catastrophic errors in main-edge and post-edge features whose shapes act as a crucial fingerprint in structural analyses of liquids. Here, we show that these errors manifest primarily due to a lack of orbital relaxation in conventional linear-response theories and that core-ionized (n - 1-electron) references, like those of electron-affinity TDDFT, can eliminate the errors in the spectral profile, even in the highest-energy parts of the post-edge. Crucially, we find that single excitations atop core-ionized references are sufficient to elucidate liquid-phase XAS spectra with semi-quantitative accuracy, opening the door for methods like electron-affinity CIS/TDDFT to be used as efficient alternatives to higher-order wave function approaches.

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

Electronic structure of strongly correlated systems: recent developments in multiconfiguration pair-density functional theory and multiconfiguration nonclassical-energy functional theory

Strong electron correlation plays an important role in transition-metal and heavy-metal chemistry, magnetic molecules, bond breaking, biradicals, excited states, and many functional materials, but it provides a significant challenge for modern electronic structure theory. The treatment of strongly correlated systems usually requires a multireference method to adequately describe spin densities and near-degeneracy correlation. However, quantitative computation of dynamic correlation with multireference wave functions is often difficult or impractical. Multiconfiguration pair-density functional theory (MC-PDFT) provides a way to blend multiconfiguration wave function theory and density functional theory to quantitatively treat both near-degeneracy correlation and dynamic correlation in strongly correlated systems; it is more affordable than multireference perturbation theory, multireference configuration interaction, or multireference coupled cluster theory and more accurate for many properties than Kohn–Sham density functional theory. This perspective article provides a brief introduction to strongly correlated systems and previously reviewed progress on MC-PDFT followed by a discussion of several recent developments and applications of MC-PDFT and related methods, including localized-active-space MC-PDFT, generalized active-space MC-PDFT, density-matrix-renormalization-group MC-PDFT, hybrid MC-PDFT, multistate MC-PDFT, spin–orbit coupling, analytic gradients, and dipole moments. We also review the more recently introduced multiconfiguration nonclassical-energy functional theory (MC-NEFT), which is like MC-PDFT but allows for other ingredients in the nonclassical-energy functional. We discuss two new kinds of MC-NEFT methods, namely multiconfiguration density coherence functional theory and machine-learned functionals.

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

The Li + CaF → Ca + LiF chemical reaction under cold conditions

The calcium monofluoride (CaF) molecule has emerged as a promising candidate for precision measurements, quantum simulation, and ultracold chemistry experiments. Inelastic and reactive collisions of laser cooled CaF molecules in optical tweezers have recently been reported and collisions of cold Li atoms with CaF are of current experimental interest. In this paper, we report ab initio electronic structure and full-dimensional quantum dynamical calculations of the Li + CaF → LiF + Ca chemical reaction. The electronic structure calculations are performed using the internally contracted multi-reference configuration-interaction method with Davidson correction (MRCI + Q). An analytic fit of the interaction energies is obtained using a many-body expansion method. Furthermore, a coupled-channel quantum reactive scattering approach implemented in hyperspherical coordinates is adopted for the scattering calculations under cold conditions. Results show that the Li + CaF reaction populates several low-lying vibrational levels and many rotational levels of the product LiF molecule and that the reaction is inefficient in the 1–100 mK regime allowing sympathetic cooling of CaF by collisions with cold Li atoms.

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

Atomic autoionization in the photo-dissociation of super-excited deuterated water molecules fragmenting into D + + O + + D

We present the relaxation dynamics of deuterated water molecules via autoionization, initiated by the absorption of a 61 eV photon, producing the very rare D + + O + + D breakup channel. We employ the COLd target recoil ion momentum spectroscopy method to measure the 3D momenta of the ionic fragments and emitted electrons from the dissociating molecule in coincidence. We interpret the results using the potential energy surfaces extracted from multi-reference configuration interaction calculations. The measured particle energy distributions can be related to a super-excited monocationic state located above the double ionization threshold of D 2 O. The autoionized electron energy shows a sharp distribution centered around 0.5 eV, which is a signature of the atomic oxygen autoionization occurring in the direct and sequential dissociation processes of D 2 O + * at a large internuclear distance. In this way, an O + radical fragment and a low-energy electron are created, both of which can trigger secondary reactions in their environment.

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

Rapidly convergent quantum Monte Carlo using a Chebyshev projector

The multireference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multireference problem. These include a logarithmically scaling metric-tree-based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants with essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.

Zhao, Zijun↗
DOE OSTI2025

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↗