Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled cluster”

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

Studying excited-state-specific perturbation theory on the Thiel set

Here, we explore the performance of a recently introduced N 5 -scaling excited-state-specific second order perturbation theory (ESMP2) on the singlet excitations of the Thiel benchmarking set. We find that, without regularization, ESMP2 is quite sensitive to π system size, performing well in molecules with small π systems but poorly in those with larger π systems. With regularization, ESMP2 is far less sensitive to π system size and shows a higher overall accuracy on the Thiel set than CC2, equation of motion-coupled cluster with singles and doubles, CC3, and a wide variety of time-dependent density functional approaches. Unsurprisingly, even regularized ESMP2 is less accurate than multi-reference perturbation theory on this test set, which can, in part, be explained by the set’s inclusion of some doubly excited states but none of the strong charge transfer states that often pose challenges for state-averaging. Beyond energetics, we find that the ESMP2 doubles norm offers a relatively low-cost way to test for doubly excited character without the need to define an active space.

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

Collisional excitation of C2H− by H2: New interaction potential and scattering calculations

Interstellar anions play an important role in astrochemistry as being tracers of the physical and chemical conditions in cold molecular clouds and circumstellar gas. The local thermodynamic equilibrium is generally not fulfilled in media where anions are detected and radiative and collisional data are required to model the observed lines. The C2H− anion has not yet been detected in the interstellar medium; however, collisional data could be used for non-LTE models that would help in identifying the most intense lines. For this purpose, we have computed the first 4D potential energy surface (PES) of the C2H−–H2 complex using an explicitly correlated coupled-cluster approach. The PES is characterized by a single deep minimum with a well-depth of 924.96 cm−1. From this interaction potential, we derived excitation cross sections and rate coefficients of C2H− induced by collisions with para- and ortho-H2. The results obtained for collisions with para-H2 are compared to previous calculations performed using a 2D-PES obtained from an average over H2 rotations.

Chemistry↗
DOE OSTI2023

A parallel, distributed memory implementation of the adaptive sampling configuration interaction method

The many-body simulation of quantum systems is an active field of research that involves several different methods targeting various computing platforms. Many methods commonly employed, particularly coupled cluster methods, have been adapted to leverage the latest advances in modern high-performance computing. Selected configuration interaction (sCI) methods have seen extensive usage and development in recent years. However, the development of sCI methods targeting massively parallel resources has been explored only in a few research works. Here, we present a parallel, distributed memory implementation of the adaptive sampling configuration interaction approach (ASCI) for sCI. In particular, we will address the key concerns pertaining to the parallelization of the determinant search and selection, Hamiltonian formation, and the variational eigenvalue calculation for the ASCI method. Load balancing in the search step is achieved through the application of memory-efficient determinant constraints originally developed for the ASCI-PT2 method. The presented benchmarks demonstrate near optimal speedup for ASCI calculations of Cr 2 (24e, 30o) with 10 6 , 10 7 , and 3 × 10 8 variational determinants on up to 16 384 CPUs. Importantly, to the best of the authors’ knowledge, this is the largest variational ASCI calculation to date.

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

Repartitioned Brillouin-Wigner perturbation theory with a size-consistent second-order correlation energy

Second-order Møller-Plesset perturbation theory (MP2) often breaks down catastrophically in small-gap systems, leaving much to be desired in its performance for myriad chemical applications such as noncovalent interactions, thermochemistry, and dative bonding in transition metal complexes. This divergence problem has reignited interest in Brillouin-Wigner perturbation theory (BWPT), which is regular at all orders but lacks size consistency and extensivity, severely limiting its application to chemistry. In this work, we propose an alternative partitioning of the Hamiltonian that leads to a regular BWPT perturbation series that, through the second order, is size-extensive, size-consistent (provided its Hartree–Fock reference is also), and orbital invariant. Here, our second-order size-consistent Brillouin-Wigner (BW-s2) approach can describe the exact dissociation limit of H 2 in a minimal basis set, regardless of the spin polarization of the reference orbitals. More broadly, we find that BW-s2 offers improvements relative to MP2 for covalent bond breaking, noncovalent interaction energies, and metal/organic reaction energies, although rivaling coupled-cluster with single and double substitutions for thermochemical properties.

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

A new generation of effective core potentials: Selected lanthanides and heavy elements

Here, we construct correlation-consistent effective core potentials (ccECPs) for a selected set of heavy atoms and ƒ-elements that are of significant current interest in materials and chemical applications, including Y, Zr, Nb, Rh, Ta, Re, Pt, Gd, and Tb. As is customary, ccECPs consist of spin orbit averaged relativistic effective potential (AREP) and effective spin-orbit (SO) terms. For the AREP part, our constructions are carried out within a relativistic coupled-cluster framework while also taking into account objective function one-particle characteristics for improved convergence in optimizations. The transferability is adjusted using binding curves of hydride and oxide molecules. We address the difficulties encountered with ƒ-elements, such as the presence of large cores and multiple near-degeneracies of excited levels. For these elements, we construct ccECPs with core valence partitioning that includes 4ƒ-subshell in the valence space. The developed ccECPs achieve an excellent balance between accuracy, size of the valence space, and transferability and are also suitable to be used in plane wave codes with reasonable energy cutoffs.

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

Multitask methods for predicting molecular properties from heterogeneous data

Data generation remains a bottleneck in training surrogate models to predict molecular properties. We demonstrate that multitask Gaussian process regression overcomes this limitation by leveraging both expensive and cheap data sources. In particular, we consider training sets constructed from coupled-cluster (CC) and density functional theory (DFT) data. We report that multitask surrogates can predict at CC-level accuracy with a reduction in data generation cost by over an order of magnitude. Of note, our approach allows the training set to include DFT data generated by a heterogeneous mix of exchange–correlation functionals without imposing any artificial hierarchy on functional accuracy. More generally, the multitask framework can accommodate a wider range of training set structures—including the full disparity between the different levels of fidelity—than existing kernel approaches based on Δ-learning although we show that the accuracy of the two approaches can be similar. Consequently, multitask regression can be a tool for reducing data generation costs even further by opportunistically exploiting existing data sources.

Chemistry↗
DOE OSTI2024

Improving gas adsorption modeling for MOFs by local calibration of Hubbard U parameters

While computational screening with density functional theory (DFT) is frequently employed for the screening of metal–organic frameworks (MOFs) for gas separation and storage, commonly applied generalized gradient approximations (GGAs) exhibit self-interaction errors, which hinder the predictions of adsorption energies. We investigate the Hubbard U parameter to augment DFT calculations for full periodic MOFs, targeting a more precise modeling of gas molecule–MOF interactions, specifically for N2, CO2, and O2. We introduce a calibration scheme for the U parameter, which is tailored for each MOF, by leveraging higher-level calculations on the secondary building unit (SBU) of the MOF. When applied to the full periodic MOF, the U parameter calibrated against hybrid HSE06 calculations of SBUs successfully reproduces hybrid-quality calculations of the adsorption energy of the periodic MOF. The mean absolute deviation of adsorption energies reduces from 0.13 eV for a standard GGA treatment to 0.06 eV with the calibrated U, demonstrating the utility of the calibration procedure when applied to the full MOF structure. Furthermore, attempting to use coupled cluster singles and doubles with perturbative triples calculations of isolated SBUs for this calibration procedure shows varying degrees of success in predicting the experimental heat of adsorption. It improves accuracy for N2 adsorption for cases of overbinding, whereas its impact on CO2 is minimal, and ambiguities in spin state assignment hinder consistent improvements of O2 adsorption. Our findings emphasize the limitations of cluster models and advocate the use of full periodic MOF systems with a calibrated U parameter, providing a more comprehensive understanding of gas adsorption in MOFs.

Chemistry↗
DOE OSTI2024

Benchmarking the exponential ansatz for the Holstein model

Polarons are quasiparticles formed as a result of lattice distortions induced by charge carriers. The single-electron Holstein model captures the fundamentals of single polaron physics. We examine the power of the exponential ansatz for the polaron ground-state wavefunction in its coupled cluster, canonical transformation, and (canonically transformed) perturbative variants across the parameter space of the Holstein model. Our benchmark serves to guide future developments of polaron wavefunctions beyond the single-electron Holstein model.

Chemistry↗
DOE OSTI2024

Extending GPU-accelerated Gaussian integrals in the TeraChem software package to f type orbitals: Implementation and applications

Here, the increasing availability of graphics processing units (GPUs) for scientific computing has prompted interest in accelerating quantum chemical calculations through their use. However, the complexity of integral kernels for high angular momentum basis functions often limits the utility of GPU implementations with large basis sets or for metal containing systems. In this work, we report the implementation of f function support in the GPU-accelerated TeraChem software package through the development of efficient kernels for the evaluation of Hamiltonian integrals. The high efficiency of the resulting code is demonstrated through density functional theory (DFT) calculations on increasingly large organic molecules and transition metal complexes, as well as coupled cluster singles and doubles calculations on water clusters. Preliminary investigations into Ni(I) catalysis with DFT and the photochemistry of MnH(CH 3 ) with complete active space self-consistent field are also carried out. Overall, our GPU-accelerated software appears to be well-suited for fast simulation of large transition metal containing systems, as well as organic molecules.

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

Electronic structure and photophysics of a supermolecular iron complex having a long MLCT-state lifetime and panchromatic absorption

Exploiting earth-abundant iron-based metal complexes as high-performance photosensitizers demands long-lived electronically excited metal-to-ligand charge-transfer (MLCT) states, but these species suffer typically from femtosecond timescale charge-transfer (CT)-state quenching by low-lying nonreactive metal-centered (MC) states. Here, we engineer supermolecular Fe(II) chromophores based on the bis(tridentate-ligand)metal(II)-ethyne-(porphinato)zinc(II) conjugated framework, previously shown to give rise to highly delocalized low-lying 3 MLCT states for other Group VIII metal (Ru, Os) complexes. Electronic spectral, potentiometric, and ultrafast pump–probe transient dynamical data demonstrate that a combination of a strong σ-donating tridentate ligand and a (porphinato)zinc(II) moiety with low-lying π*-energy levels, sufficiently destabilize MC states and stabilize supermolecular MLCT states to realize Fe(II) complexes that express 3 MLCT state photophysics reminiscent of their heavy-metal analogs. The resulting Fe(II) chromophore archetype, FeNHCPZn, features a highly polarized CT state having a profoundly extended 3MLCT lifetime (160 ps), 3 MLCT phosphorescence, and ambient environment stability. Density functional and domain-based local pair natural orbital coupled cluster [DLPNO-CCSD(T)] theory reveal triplet-state wavefunction spatial distributions consistent with electronic spectroscopic and excited-state dynamical data, further underscoring the dramatic Fe metal-to-extended ligand CT character of electronically excited FeNHCPZn. This design further prompts intense panchromatic absorptivity via redistributing high-energy absorptive oscillator strength throughout the visible spectral domain, while maintaining a substantial excited-state oxidation potential for wide-ranging photochemistry––highlighted by the ability of FeNHCPZn to photoinject charges into a SnO 2 /FTO electrode in a dye-sensitized solar cell (DSSC) architecture. Concepts enumerated herein afford opportunities for replacing traditional rare-metal–based emitters for solar-energy conversion and photoluminescence applications.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Methanediol from cloud-processed formaldehyde is only a minor source of atmospheric formic acid

Atmospheric formic acid is severely underpredicted by models. A recent study proposed that this discrepancy can be resolved by abundant formic acid production from the reaction ( 1 ) between hydroxyl radical and methanediol derived from in-cloud formaldehyde processing and provided a chamber-experiment-derived rate constant, k 1 = 7.5 × 10 −12 cm 3 s −1 . High-level accuracy coupled cluster calculations in combination with E,J -resolved two-dimensional master equation analyses yield k 1 = (2.4 ± 0.5) × 10 −12 cm 3 s −1 for relevant atmospheric conditions ( T = 260–310 K and P = 0–1 atm). We attribute this significant discrepancy to HCOOH formation from other molecules in the chamber experiments. More importantly, we show that reversible aqueous processes result indirectly in the equilibration on a 10 min. time scale of the gas-phase reaction HCHO + H 2 O ⇌ HOCH 2 OH (2) with a HOCH 2 OH to HCHO ratio of only ca . 2%. Although HOCH 2 OH outgassing upon cloud evaporation typically increases this ratio by a factor of 1.5–5, as determined by numerical simulations, its in-cloud reprocessing is shown using a global model to strongly limit the gas-phase sink and the resulting production of formic acid. Based on the combined findings in this work, we derive a range of 1.2–8.5 Tg/y for the global HCOOH production from cloud-derived HOCH 2 OH reacting with OH. The best estimate, 3.3 Tg/y, is about 30 times less than recently reported. The theoretical equilibrium constant K eq (2) determined in this work also allows us to estimate the Henry’s law constant of methanediol (8.1 × 10 5 M atm −1 at 280 K).

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

New reactions of diazene and related species for modelling combustion of amine fuels

Potential energy surfaces for reactions involving N 2 H 2 isomers of diazene (diimide) have been explored using density functional theory, with energies based on coupled-cluster theory. A focus is on processes that create or consume these species, and isomerisation between the E (trans) and Z (cis) forms of HNNH. Furthermore, these include isomerisation and dissociation pathways for HNNH, addition of H atoms to form N 2 H 3 , abstraction by H atoms yielding short-lived NNH, and abstraction reactions of H with N 2 H 3 . Transition state and capture theories are applied for high-pressure-limiting behaviour, while low-pressure and falloff regions are characterised via the methods of Troe and coworkers. Rate constants and thermochemistry are provided to improve models of diamine chemistry, relevant to the combustion of NH 3 especially at high concentrations, high pressures or under reducing conditions. Results indicate that amine radical recombination mainly yields the E HNNH isomer, while H-abstraction from N 2 H 3 results in E HNNH and H 2 NN. However, at elevated temperature E → Z isomerisation becomes competitive, and Z HNNH, being more reactive, acts to enhance the diazene consumption rate.

10 SYNTHETIC FUELS↗
DOE OSTI2020

From synchrotrons for XFELs: the soft x-ray near-edge spectrum of the ESCA molecule

A predictive understanding of soft x-ray near-edge absorption spectra of small molecules is an enduring theoretical challenge and of current interest for x-ray probes of molecular dynamics. We report the experimental absorption spectrum for the electron spectroscopy for chemical analysis (ESCA) molecule (ethyl trifluoroacetate) near the carbon 1s absorption edge between 285-300 eV. Furthermore, the ESCA molecule with four chemically distinct carbon sites has previously served as a theoretical benchmark for photoelectron spectra and now for photoabsorption spectra. We report a simple edge-specific approach for systematically expanding standard basis sets to properly describe diffuse Rydberg orbitals and the importance of triple excitations in equation-of-motion coupled-cluster calculations of the energy interval between valence and Rydberg excitations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Resonant x-ray emission across the L 3 edge of uranium compounds

A narrow bandwidth x-ray beamline and a multi-crystal von Hamos spectrometer were used to record x-ray absorption and x-ray emission (XES) across the uranium L 3 edge of UO 2 , UO 3 , Cs 2 UO 2 Cl 4 , and Cs 2 UCl 6 . Measurements were made over 17150–17250 eV with an instrumental resolution of ∼2 eV. This resolution allowed the use of Lorentzian peak fits to the L 3 N 4 and L 3 N 5 characteristic x-ray fluorescence lines that have ∼12–13 eV lifetime widths. The fluorescence yields of the four compounds display strong white lines near threshold and multiple scattering features at higher energies. The XES spectra were fit with three Lorentzians, two for the L 3 N 4 and L 3 N 5 lines and one for inelastic scattering by 2p–6d excitations and 4d–6d final states. The intensities of the white lines were largely due to the inelastic scattering component. To support the measurements, relativistic equation-of-motion coupled-cluster and restricted active space configuration interaction calculations were performed for uranium core-excitation energies in Cs 2 UO 2 Cl 4 and Cs 2 UCl 6 . These calculations with rigorous treatments of relativistic effects are shown to provide accurate uranium L 3 -edge binding energies, L 3 N 5 emission energies, as well as the energy losses in the inelastic scattering process. The measured and calculated results show variations between compounds with U(IV) oxidation states that contain 6d 0 5f 2 electrons in their ground configurations (UO 2 and Cs 2 UCl 6 ) compared with U(VI) compounds (UO 3 and Cs 2 UO 2 Cl 4 ) with empty 5f and 6d configurations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Merging ab initio theory and few-body approach for (d, p) reactions

A framework for A(d, p)B reactions is introduced by merging the microscopic approach to computing the properties of the nucleon-target systems and the three-body n + p + A reaction formalism, thus providing a consistent link between the reaction cross sections and the underlying microscopic structure. In this first step toward a full microscopic description, we focus on the inclusion of the neutron-target microscopic properties. The properties of the neutron-target subsystem are encapsulated in the Green's function which is computed with the coupled Cluster theory using a chiral nucleon-nucleon and three-nucleon interactions. Subsequently, this many-body information is introduced in the few-body Green's Function Transfer approach to (d, p) reactions. Our benchmarks on stable targets 40,48 Ca show an excellent agreement with the data. We then proceed to make specific predictions for (d, p) on neutron rich 52,54 Ca isotopes. These predictions are directly relevant to testing the new magic numbers N = 32, 34 and are expected to be feasible in the first campaign of the projected FRIB facility.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Adaptive pruning-based optimization of parameterized quantum circuits

Abstract Variational hybrid quantum–classical algorithms are powerful tools to maximize the use of noisy intermediate-scale quantum devices. While past studies have developed powerful and expressive ansatze, their near-term applications have been limited by the difficulty of optimizing in the vast parameter space. In this work, we propose a heuristic optimization strategy for such ansatze used in variational quantum algorithms, which we call ‘parameter-efficient circuit training (PECT)’. Instead of optimizing all of the ansatz parameters at once, PECT launches a sequence of variational algorithms, in which each iteration of the algorithm activates and optimizes a subset of the total parameter set. To update the parameter subset between iterations, we adapt the Dynamic Sparse Reparameterization scheme which was originally proposed for training deep convolutional neural networks. We demonstrate PECT for the Variational Quantum Eigensolver, in which we benchmark unitary coupled-cluster ansatze including UCCSD and k -UpCCGSD, as well as the Low-Depth Circuit Ansatz (LDCA), to estimate ground state energies of molecular systems. We additionally use a layerwise variant of PECT to optimize a hardware-efficient circuit for the Sycamore processor to estimate the ground state energy densities of the one-dimensional Fermi-Hubbard model. From our numerical data, we find that PECT can enable optimizations of certain ansatze that were previously difficult to converge and more generally can improve the performance of variational algorithms by reducing the optimization runtime and/or the depth of circuits that encode the solution candidate(s).

Physics↗
DOE OSTI2026

Ionization potential of radium monofluoride

The ionization potential (IP) of radium monofluoride (RaF) was measured to be 4.969(2)[10] eV, revealing a relativistic enhancement in the series of alkaline earth monofluorides. The results are in agreement with a relativistic coupled-cluster prediction of 4.981(7) eV, incorporating up to quantum electrodynamics corrections. Using the same computational methodology, an improved calculation for the dissociation energy (𝐷 0 ) of 5.54(5) eV is presented. This confirms that RaF joins the group of diatomic molecules for which 𝐷 0 > IP, paving the way for precision control and interrogation of its Rydberg states.

Wilkins, S. G. [Massachusetts Inst. of Technology ↗
DOE OSTI2026

Toward scalable quantum computations of atomic nuclei

We solve the nuclear two-body and three-body bound states via quantum simulations of pionless effective field theory on a lattice in position space. While the employed lattice remains small, the usage of local Hamiltonians including two- and three-body forces ensures that the number of Pauli terms scales linearly with increasing numbers of lattice sites. We use an adaptive ansatz grown from unitary coupled cluster theory to parametrize the ground states of the deuteron and 3 He, compute their corresponding energies, and analyze the scaling of the required computational resources. Our quantum simulations reproduce exact benchmarks for 2 H and 3 He within 100 keV, requiring at most 30 layers in the ansatz and thus resulting in modest circuit depths. Additionally, we find the number of shots required to reach a given precision scales linearly in the lattice size and more mildly in the system size. Furthermore, based on the agreement with exact benchmarks and mild scaling, we conclude that this can be an efficient, scalable approach for quantum computations of nuclear ground states, particularly to prepare initial states for quantum phase estimation or other filtering algorithms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗