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 469 records · Page 26

Stochastic tensor contraction for quantum chemistry

Many computational methods in ab initio quantum chemistry are formulated in terms of high-order tensor contractions, whose cost determines the size of system that can be studied. We introduce stochastic tensor contraction to perform such operations with greatly reduced cost, and present its application to the gold-standard quantum chemistry method, coupled cluster theory with up to perturbative triples. For total energy errors more stringent than chemical accuracy, we reduce the computational scaling to that of mean-field theory, while starting to approach the mean-field absolute cost, thereby challenging the existing cost-to-accuracy landscape. Benchmarks against state-of-the-art local correlation approximations further show that we achieve an order-of-magnitude improvement in both total computation time and error, with significantly reduced sensitivity to system dimensionality and electron delocalization. We conclude that stochastic tensor contraction is a powerful computational primitive to accelerate a wide range of quantum chemistry.

Chemical Physics (physics.chem-ph)↗

Benchmarking quantum trial wavefunctions for phaseless auxiliary-field quantum Monte Carlo

The phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) method is a stochastic imaginary-time projection technique for computing ground-state properties of strongly correlated quantum systems, with accuracy that depends critically on the choice of trial wavefunction. Here, we investigate ph-AFQMC with trial states prepared using parameterized quantum circuits. In this work, we present a comprehensive benchmarking study of quantum trial wavefunctions spanning unitary coupled-cluster, Hamiltonian-informed, Jastrow-inspired, and adaptively constructed ansatze. The benchmarking evaluates accuracy, expressibility, and scalability of these ansatze within the QC-AFQMC framework. We test these ansatze on linear hydrogen chains under bond stretching and find that several ansatz families produce chemically accurate ph-AFQMC energies across the dissociation curve. We have performed simulations using the CUDA-Q quantum development platform on the GPU partition of the Perlmutter supercomputer. When comparing ansatze at similar numbers of variational parameters, we find that different ansatz families yield comparable ph-AFQMC results despite exhibiting substantially different variational energies, optimization costs, and circuit depths. Our results indicate that the variational energy of an ansatz is not always a reliable indicator of its quality for ph-AFQMC and reveal instances of over-parameterization. In the strongly correlated regime, trial wavefunctions obtained from adaptive ansatze, exemplified here by ADAPT-VQE with the UCCSD operator pool, can outperform their fixed-ansatz counterparts (UCCSD) in terms of projected energies while using substantially more compact circuits, providing a flexible route to optimize quantum resources within the ph-AFQMC framework.

Rofougaran, Rod [LBNL, Berkeley; Columbia U.; PNL,↗

Determinants for proton location and electron coupled proton transfer in hydrogen bonded pentafluorophenol-anion clusters

This work aims to reveal the determinant factors for proton locations and electron coupled proton transfer (ECPT) in biologically relevant hydrogen bonded systems. Pentafluorophenol-anion clusters [C6F5O-•H+•A-]- are chosen to model active sites of biological functions, with the anion A- being systematically varied to ensure the proton affinities (PAs) of the anions well cover the referenced PA of C6F5O- from being appreciably smaller, similar, and to significantly larger. Negative ion photoelectron spectroscopy of these clusters provides spectroscopic evidence showing that proton location in the anionic state is largely following the PA prediction, while ECPT is observed only for the clusters with the anion possessing the electron binding energy (EBE) significantly larger than that of the referenced C6F5O-. Theoretical calculations suggest these clusters are stabilized by forming single strong hydrogen bond between donor and acceptor, and the associated charge and MO analyses fully support the experimental observations. The current holistic cluster model study indicates that PA is the right determinant that can be used to predict the proton location and describe hydrogen bonding structures, while both PA and EBE of anionic groups play important roles in facilitating ECPT process. This work was supported by U.S. Department of Energy (DOE), Office of Basic Energy Sciences, Division of Chemical Science, Geosciences, and Biosciences (X.-B.W.), and performed using EMSL, a national scientific user facility sponsored by DOE’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory, which is operated by Battelle Memorial Institute for the DOE. Part of this work was supported by the National Natural Science Foundation of China (No-11604046 for JZ; No. 11727810 for ZRS), as well as “Chenguang Program” by Shanghai Education Development Foundation, Shanghai Municipal Education Commission (No. 16CG38) and Natural Science Foundation of Shanghai (No. 16ZR1448100) (JZ).

Zhang, Jian↗

Isotope Effects in the Zundel–Eigen Isomerization of H + (H 2 O) 6

The isomerization pathway between the energetically low-lying Zundel and Eigen isomers of the protonated water hexamer was investigated using high-level ab initio calculations including a treatment of zero-point corrections. On the basis of these calculations, the Zundel–Eigen isomerization was found to proceed through a stable intermediate isomer, which consists of a four-membered ring with two single acceptor water molecules. Further, the inclusion of vibrational zero-point energy is shown to be important for accurately establishing the relative energies of the three relevant isomers involved in the Zundel–Eigen isomerization. Diffusion Monte Carlo calculations including anharmonic vibrational effects show that all three isomers of H + (H 2 O) 6 and D + (D 2 O) 6 have well-defined structures. The energetic ordering of the three isomers changes upon deuteration. The implications of these results for the vibrational spectra of H + (H 2 O) 6 and D + (D 2 O) 6 are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the prospects of optical cycling in diatomic cations: effects of transition metals, spin–orbit couplings, and multiple bonds

Molecules with optical cycling centers (OCCs) are highly desirable in the context of fundamental studies as well as applications (e.g., quantum computing) because they can be effectively cooled to very low temperatures by repeated absorption and emission (hence, cycling). Charged species offer additional advantages for experimental control and manipulation. Here, we present a systematic computational study of a series of diatomic radical-cations made of a d-block metal and a p-block ligand, that are isoelectronic (in their valence shell) to the successfully laser-cooled neutral molecules. Using high-level electronic structure methods, we characterize state and transition properties of low-lying electronic states and compute Franck-Condon factors. The computed branching ratios and radiative lifetimes reveal that the electronic transitions analogous to those successfully used in the laser cooling of neutral molecules are less than optimal in the cations. We propose alternative transitions suitable for optical cycling and highlight trends that could assist future designs of OCCs in charged or neutral molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗

Identifying Sample Provenance From SEM/EDS Automated Particle Analysis via Few-Shot Learning Coupled With Similarity Graph Clustering

Automated particle analysis (APA) provides a vast amount of compositional data via energy-dispersive X-ray spectroscopy along with size and shape data via scanning electron microscopy for individual particles in a sample. In many instances, APA data are leveraged to support identification of the source of a sample based on the detection of particles of a specific composition. Often, the particles that provide context make up a minuscule portion of the sample. Additionally, the interpretation of complex samples can be difficult due to the diversity of compositions both in the mixture and within a particle. In this work, we demonstrate a method to compute and cluster similarity graphs that describe inter-particle relationships within a sample using a multi-modal few-shot learning neural network. Here, as a proof-of-concept, we show that samples known to have been exposed to gunshot residue can be distinguished from samples occasionally mistaken for gunshot residue. Our workflow builds upon standard APA techniques and data processing methods to unveil additional information in a readily interpretable and quantitatively comparable format.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Reproducibility of calculations on Li species with correlation-consistent basis sets

Here we have found that two different basis set series, both labeled cc-pVxZ, are widely available and widely used for calculation on Li compounds such that calculations that should be the same are different, causing a reproducibility problem. Here we study the magnitude of the disagreement of these nominally-the-same basis sets for wave function and density functional calculations of Li clusters. We find deviations among the basis sets of as large as 22 millihartrees. In most cases one obtains lower energies with the noncanonical choices. Therefore, it is important to specify which version of the basis set is used for calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanistic implications of excited high-spin states, spin–spin coupling, and differential [2Fe–2S] + cluster temperature relaxations in the electron-bifurcating NfnABC from Thermococcus sibiricus

Electron bifurcation (EB) is a mechanism of biological energy transduction in which multiple oxidation–reduction (redox) reactions are thermodynamically coupled within a single enzyme, enabling the enzyme to harness the excess free energy from an exergonic process to drive an endergonic process. Because of this unprecedented chemistry, there is interest to translate EB principles to artificial and bioengineered systems, but a hurdle is that knowledge pertaining to the fundamental design principles of EB enzymes remains scarce. Here, we investigated the fundamental physical and electronic properties of electron transfer sites in a spectroscopically uncharacterized member of the BfuABC family of EB enzymes, the NADH-dependent reduced-ferredoxin:NADP + oxidoreductase from Thermococcus sibiricus (Tsi NfnABC). Cryo-EM structures of Tsi NfnABC previously demonstrated that it contains twelve redox cofactors: two flavins (one FAD and one FMN), eight [4Fe–4S] clusters, and two [2Fe–2S] clusters. The FMN, one [4Fe–4S] cluster, and one [2Fe–2S] cluster comprise the bifurcating active site termed the electron-bifurcating flavobicluster (BF-FBC), which is found in all BfuABC family members. By using electron paramagnetic resonance spectroscopy, we identified spectral signatures originating from interactions between the FMN radical and [4Fe–4S] + cluster in the BF-FBC and observed temperature dependent behavior of the BF-FBC's [2Fe–2S] + cluster indicative of moderately slow spin–lattice relaxation. Additionally, we uncovered numerous spectral features corresponding to half-integer, S > ½ spin states of [4Fe–4S] + clusters, including one attributable to the consequences of lysine-ligation of a [4Fe–4S] cluster unique to NfnABC. We contextualize these findings to electron transfer theory and NfnABC's structure. Our insights further the understanding of how enzymes are designed to exert control over electron transfer to conduct thermodynamically challenging reactions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Constructing a square-like copper cluster to boost C–C coupling for CO 2 electroreduction to ethylene

The CO 2 electroreduction reaction (CO 2 ER) to ethylene (C 2 H 4 ) offers the dual promise of lowering CO 2 emission while storing energy from renewable electricity, for which the development of highly efficient electrocatalysts is of great significance. Herein, by means of density functional theory (DFT) computations, we designed an electrocatalyst for CO 2 -to-C 2 H 4 conversion by anchoring a Cu 5 cluster supported on a MoS 2 monolayer with an S monovacancy (Cu 5 @MoS 2 ). Our results revealed that one Cu atom of the Cu 5 cluster was embedded into the framework of the defective MoS 2 monolayer, while the other four Cu atoms form a square-like island over the substrate surface. Interestingly, the C–C coupling between two *CO species can easily occur on the unique square-like active site with a low kinetic barrier of 0.56 eV to form the key *C 2 O 2 intermediate, which can then be hydrogenated to the C 2 H 4 product with a very low limiting potential (–0.32 eV). Significantly, alkaline conditions (pH = 13) are beneficial to further promote C 2 H 4 synthesis. Finally, our work may offer a new avenue to precisely modulate the structures of Cu clusters for converting CO 2 into high-value target products.

30 DIRECT ENERGY CONVERSION↗

Symmetry-projected cluster mean-field theory applied to spin systems

We introduce S z spin-projection based on cluster mean-field theory and apply it to the ground state of strongly correlated spin systems. In cluster mean-fields, the ground state wavefunction is written as a factorized tensor product of optimized cluster states. In previous work, we have focused on unrestricted cluster mean-field, where each cluster is S z symmetry adapted. We here remove this restriction by introducing a generalized cluster mean-field (GcMF) theory, where each cluster is allowed to access all S z sectors, breaking S z symmetry. In addition, a projection scheme is used to restore global S z , which gives rise to the S z spin-projected generalized cluster mean-field (S z GcMF). Both of these extensions contribute to accounting for inter-cluster correlations. We benchmark these methods on the 1D, quasi-2D, and 2D J 1 – J 2 and XXZ Heisenberg models. Furthermore, our results indicate that the new methods (GcMF and S z GcMF) provide a qualitative and semi-quantitative description of the Heisenberg lattices in the regimes considered, suggesting them as useful references for further inter-cluster correlations, which are discussed in this work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Binding of noble gas atoms by superhalogens

Because of their closed shells, noble gas (Ng) atoms (Ng = Ne, Ar, Kr, and Xe) seldom take part in chemical reactions, yet finding such mechanisms not only is of scientific interest but also has practical significance. Following a recent work by Mayer et al. [Proc. Natl. Acad. Sci. U. S. A. 116, 8167–8172 (2019)] on the room temperature binding of Ar to a superelectrophilic boron site embedded in a negative ion complex, B 12 (CN) 11 - , we have systematically studied the effect of cluster size and terminal ligands on the interaction of Ng by focusing on B 12 X 11 (Ng) (X = H, CN, and BO) and B 12 X 10 (Ng) 2 (X = CN and BO) whose stabilities are governed by the Wade–Mingos rule and on C 5 BX 5 (Ng) (X = H, F, and CN) and C 4 B 2 (CN) 4 (Ng) 2 whose stabilities are governed by the Huckel’s aromaticity rule. Our conclusion, based on density functional theory, is that both the cluster size and the terminal ligands matter—the interaction between the cluster and the Ng atoms becomes stronger with increasing cluster size and the electron affinity of the terminal ligands. Our studies also led to a counter-intuitive finding—removing multiple terminal ligands can enable electrophilic centers to bind multiple Ng atoms simultaneously without compromising their binding strength.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)↗

Breaking covalent bonds in the context of the many-body expansion (MBE). I. The purported “first row anomaly” in XH n (X = C, Si, Ge, Sn; n = 1–4)

We present a new, novel implementation of the Many-Body Expansion (MBE) to account for the breaking of covalent bonds, thus extending the range of applications from its previous popular usage in the breaking of hydrogen bonds in clusters to molecules. A central concept of the new implementation is the in situ atomic electronic state of an atom in a molecule that casts the one-body term as the energy required to promote it to that state from its ground state. The rest of the terms correspond to the individual diatomic, triatomic, etc., fragments. Its application to the atomization energies of the XH n series, X = C, Si, Ge, Sn and n = 1–4, suggests that the (negative, stabilizing) 2-B is by far the largest term in the MBE with the higher order terms oscillating between positive and negative values and decreasing dramatically in size with increasing rank of the expansion. The analysis offers an alternative explanation for the purported “first row anomaly” in the incremental H n–1 X–H bond energies seen when these energies are evaluated with respect to the lowest energy among the states of the XH n molecules. Due to the “flipping” of the ground/first excited state between CH 2 ( 3 B 1 ground state, 1 A 1 first excited state) and XH 2 , X = Si, Ge, Sn ( 1 A 1 ground state, 3 B 1 first excited state), the overall picture does not exhibit a “first row anomaly” when the incremental bond energies are evaluated with respect to the molecular states having the same in situ atomic states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synthesis of Co‐Doped MoS 2 Monolayers with Enhanced Valley Splitting

Abstract Internal magnetic moments induced by magnetic dopants in MoS 2 monolayers are shown to serve as a new means to engineer valley Zeeman splitting (VZS). Specifically, successful synthesis of monolayer MoS 2 doped with the magnetic element Co is reported, and the magnitude of the valley splitting is engineered by manipulating the dopant concentration. Valley splittings of 3.9, 5.2, and 6.15 meV at 7 T in Co‐doped MoS 2 with Co concentrations of 0.8%, 1.7%, and 2.5%, respectively, are achieved as revealed by polarization‐resolved photoluminescence (PL) spectroscopy. Atomic‐resolution electron microscopy studies clearly identify the magnetic sites of Co substitution in the MoS 2 lattice, forming two distinct types of configurations, namely isolated single dopants and tridopant clusters. Density functional theory (DFT) and model calculations reveal that the observed enhanced VZS arises from an internal magnetic field induced by the tridopant clusters, which couples to the spin, atomic orbital, and valley magnetic moment of carriers from the conduction and valence bands. The present study demonstrates a new method to control the valley pseudospin via magnetic dopants in layered semiconducting materials, paving the way toward magneto‐optical and spintronic devices.

Zhou, Jiadong↗

Hydrogen Bond Benchmark: Focal‐Point Analysis and Assessment of DFT Functionals

We performed a hierarchical, convergent ab initio benchmark study and systematically analyzed the performance of density functional approximations for describing hydrogen bonds in small neutral, cationic, and anionic complexes, as well as in larger systems involving amide, urea, deltamide, and squaramide moieties. Focal point analyses (FPA), extrapolating to the ab initio limit, were carried out using correlated wave function methods up to CCSDT(Q) for the small complexes and CCSD(T) for the larger systems, together with correlation-consistent Gaussian basis sets up to the complete basis set limit. Optimized geometries and vibrational frequencies were obtained at the CCSD(T) level. The resulting FPA hydrogen-bond energies converge within a few tenths of a kcal mol −1 . These reference data were used to evaluate 60 density functionals (including 12 dispersion-corrected), spanning the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. Overall, the meta-hybrid M06-2X provides the best performance for both hydrogen bond energies and geometries, while the dispersion-corrected GGAs BLYP-D3(BJ) and BLYP-D4 also yield accurate hydrogen-bond data and can serve as cost-effective options for studying large and complex systems.

coupled cluster theory↗

Biased Spin-State Energetics of Fe(II) Molecular Complexes within Density-Functional Theory and the Linear-Response Hubbard U Correction

The spin-state energetics of six Fe(II) molecular complexes are computed using the linear-response Hubbard U approach within DFT. The adiabatic energy differences, ΔE H-L , between the high-spin (S = 2) and the low-spin (S = 0) states are computed and compared with accurate-coupled cluster-corrected CASPT2 results. We show that DFT+U fails in correctly capturing the ground state for strong-field ligands yielding ΔE H-L that are almost constant throughout the molecular series. This bias toward high spin together with the metal/ligand charge transfer upon U correction are here quantified and explained using molecular orbital diagrams involving both σ- and π-bonding interactions. Furthermore, with increasing ligand-field strengths this bias also increases owing to the stronger molecular character of the metal/ligand Kohn–Sham orbitals thus resulting in large deviations from the reference larger than 4 eV. Smaller values of U can be employed to mitigate this effect and recover the right energetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Classical Preoptimization Approach for ADAPT-VQE: Maximizing the Potential of High-Performance Computing Resources to Improve Quantum Simulation of Chemical Applications

The ADAPT-VQE algorithm is a promising method for generating a compact ansatz based on derivatives of the underlying cost function, and it yields accurate predictions of electronic energies for molecules. In this work, we report the implementation and performance of ADAPT-VQE with our recently developed sparse wave function circuit solver (SWCS) in terms of accuracy and efficiency for molecular systems with up to 52 spin orbitals. The SWCS can be tuned to balance computational cost and accuracy, which extends the application of ADAPT-VQE for molecular electronic structure calculations to larger basis sets and a larger number of qubits. Using this tunable feature of the SWCS, we propose an alternative optimization procedure for ADAPT-VQE to reduce the computational cost of the optimization. Furthermore, by preoptimizing a quantum simulation with a parametrized ansatz generated with ADAPT-VQE/SWCS, we aim to utilize the power of classical high-performance computing in order to minimize the work required on noisy intermediate-scale quantum hardware, which offers a promising path toward demonstrating quantum advantage for chemical applications.

ADAPT-VQE↗

Mapping the Electronic Structure of the Uranium(VI) Dinitride Molecule, UN 2

A combined anion photoelectron spectroscopic and relativistic electronic-structure computational study of the electronic structure of the UN 2 molecule is reported. Since the photoelectron spectrum of uranium dinitride’s negative ion, UN 2 - , (UN 2 plus an additional electron) directly reflects the electronic structure of neutral UN 2 , we have measured and relied upon the photoelectron spectrum of the UN 2 – anion as a means of mapping the electronic structure of neutral UN 2 . In addition to the electron affinity of the UN 2 ground state (the energy released when attaching an electron to UN 2 molecule), energy levels of the UN 2 excited states (electronic states with higher energies than the ground state) were well characterized by the close interplay between experiment and high-level theory. We found that both electron attachment and electronic excitation significantly bend the UN 2 molecule and elongate its U≡N bond. Implications for the activation of UN 2 are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗