Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Quantum Benchmarking”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Two-Component Multireference Restricted Active Space Configuration Interaction for the Computation of L-Edge X-ray Absorption Spectra

X-ray absorption spectroscopy is a powerful probe of local electronic and nuclear structure, providing insights into chemical processes. The theoretical prediction and interpretation of metal L-edge X-ray absorption spectra is complicated by both relativistic effects, including spin-orbit coupling, and by the multiconfigurational nature of the states involved. This work details an exact two-component multireference restricted active space configuration interaction (X2C-MRRASCI) scheme that uses an exact two-component state average complete active space self-consistent-field method (X2C-SA-CASSCF), which includes the spin-orbit coupling in a variational manner, for the accurate description of the electronic structure before using a restricted active space configuration interaction method to describe the core excited states of the X-ray spectrum. Benchmark calculations are presented for a series of iron containing complexes, with results showing key features of the spectrum being reproduced, including ligand to metal charge transfer and shake-up excitations.

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

Optical second harmonic generation in anisotropic multilayers with complete multireflection of linear and nonlinear waves using ♯SHAARP.ml package

Abstract Optical second harmonic generation (SHG) is a nonlinear optical effect widely used for nonlinear optical microscopy and laser frequency conversion. Closed-form analytical solution of the nonlinear optical responses is essential for evaluating materials whose optical properties are unknown a priori. A recent open-source code, ♯SHAARP. si , can provide such closed form solutions for crystals with arbitrary symmetries, orientations, and anisotropic properties at a single interface. However, optical components are often in the form of slabs, thin films on substrates, and multilayer heterostructures with multiple reflections of both the fundamental and up to ten different SHG waves at each interface, adding significant complexity. Many approximations have therefore been employed in the existing analytical approaches, such as slowly varying approximation, weak reflection of the nonlinear polarization, transparent medium, high crystallographic symmetry, Kleinman symmetry, easy crystal orientation along a high-symmetry direction, phase matching conditions and negligible interference among nonlinear waves, which may lead to large errors in the reported material properties. To avoid these approximations, we have developed an open-source package named Second Harmonic Analysis of Anisotropic Rotational Polarimetry in Multilayers (♯SHAARP. ml ). The reliability and accuracy are established by experimentally benchmarking with both the SHG polarimetry and Maker fringes using standard and commonly used nonlinear optical materials as well as twisted 2-dimensional heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

On symmetry-resolved generalized entropies

Symmetry-resolved entanglement, capturing the refined structure of quantum entanglement in systems with global symmetries, has attracted a lot of attention recently. In this manuscript, introducing the notion of symmetry-resolved generalized entropies, we aim to develop a computational framework suitable for the study of excited state symmetry-resolved entanglement as well as the dynamical evolution of symmetry-resolved entanglement in symmetry-preserving out-of-equilibrium settings. We illustrate our framework using the example of (1+1)-d free massless compact boson theory, and benchmark our results using lattice computation in the XX chain. As a byproduct, our computational framework also provides access to the probability distribution of the symmetry charge contained within a subsystem and the corresponding full counting statistics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Towards exact finite temperature electronic structure in solids and molecules (Final Technical Report)

This report describes the University of Iowa portion of a project that is now continuing at Michigan State University. We are developing novel methods, algorithms, and software to enable simulations of molecules and materials at high temperature. This is a key challenge in chemistry and materials science. By refining an approach called Density Matrix Quantum Monte Carlo (DMQMC), we developed faster and more accurate ways to conduct these simulations. These advances will help us understand how temperature affects the behavior of electrons, chemical bonds, and phase transitions in solids and molecules. These breakthroughs are especially important for applications where light and heat drive chemical reactions, superconductivity, and materials used in energy and sensing. In addition, this project involved the development of the open-source HANDE-QMC software package, supporting the broader community in benchmarking and developing finite-temperature electronic structure methods.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Meta-Generalized Gradient Approximation for the Cavity-Dependent Exchange-Correlation Interaction in Strongly Coupled Light–Matter Systems

Strong light–matter coupling in optical cavities enables the manipulation of chemical and physical properties without altering molecular composition. Theoretical modeling of such phenomena requires exchange-correlation (XC) functionals that account for both electron–electron and electron–photon (ep) interactions within quantum electrodynamical density functional theory (QEDFT). In this work, we develop a meta-generalized gradient approximation (meta-GGA) specifically targeting the cavity-dependent XC interaction in strongly coupled light–matter systems. This novel approximation is built upon a new semilocal polarizability approximation, which draws from the jellium-with-a-gap model, and can be extended to a “global hybrid” variant that goes beyond the isotropic model from previous approximations. The polarizability model yields significantly improved dispersion coefficients and benchmark calculations with the cavity-dependent XC functional demonstrate improved agreement with QED Hartree–Fock (QED-HF) reference energies. Application to the regioselectivity of brominated nitrobenzene intermediates reveals the functional’s capacity to capture cavity-induced energetic shifts. In conclusion, our results advance the Jacob’s ladder of functionals for QEDFT and provide a practical tool for modeling polaritonic chemistry.

Approximation↗
DOE OSTI2023

Solvent Dynamics Are Critical to Understanding Carbon Dioxide Dissolution and Hydration in Water

Simulations of carbon dioxide (CO 2 ) in water may aid in understanding the impact of its accumulation in aquatic environments and help advance technologies for carbon capture and utilization (via, e.g., mineralization). Quantum mechanical (QM) simulations based on static molecular models with polarizable continuum solvation poorly reproduce the energetics of CO 2 hydration to form carbonic acid in water, independent of the level of QM theory employed. Only with density-functional-theory-based molecular dynamics and rare-event sampling, followed by energy corrections based on embedded correlated wavefunction theory (in conjunction with density functional embedding theory), can a close agreement between theory and experiment be achieved. Such multilevel simulations can serve as benchmarks for simpler, less costly models, giving insight into potential errors of the latter. The strong influence of sampling/averaging over dynamical solvent configurations on the energetics stems from the difference in polarity of both the transition state and product (both polar) versus the reactant (nonpolar). When a solute undergoes a change in polarity during reaction, affecting its interaction with the solvent, careful assessment of the energetic contribution of the solvent response to this change is critical. Here, we show that static models (without structural sampling) that incorporate three explicit water molecules can yield far superior results than models with more explicit water molecules because fewer water molecules yield less configurational artifacts. Static models intelligently incorporating both explicit (molecules directly participating in the reaction) and implicit solvation, along with a proper QM theory, e.g., CCSD(T) for closed-shell systems, can close the accuracy gap between static and dynamic models.

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

Li1−xNiO2 Many-body DMC Benchmark Dataset

The dataset contains all numerical data generated in support of the manuscript “Many‑body Benchmark of Electronic Charge and Spin Densities for Li1–xNiO2​” (Journal of Chemical Theory and Computation, DOI: 10.1021/acs.jctc.5c02097, URL: https://pubs.acs.org/doi/10.1021/acs.jctc.5c02097). The materials included in this repository are: 1. Data files used to produce all figures and tables in the main manuscript and supporting information. 2. Benchmark density‑functional theory (DFT) datasets used for the charge‑ and spin‑density analyses. 3. Reference many‑body diffusion Monte Carlo (DMC) calculations and associated input/output files.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Multireference Quantum Krylov Algorithm for Strongly Correlated Electrons

We introduce a multireference selected quantum Krylov (MRSQK) algorithm suitable for quantum simulation of many-body problems. MRSQK is a low-cost alternative to the quantum phase estimation algorithm that generates a target state as a linear combination of non-orthogonal Krylov basis states. This basis is constructed from a set of reference states via real-time evolution; thus, avoiding the numerical optimization of parameters. An efficient algorithm for the evaluation of the off-diagonal matrix elements of the overlap and Hamiltonian matrices is discussed and a selection procedure is introduced to identify a basis of orthogonal references that ameliorates the linear dependency problem. Preliminary benchmarks on linear H 6 , H 8 , and BeH 2 indicate that MRSQK can predict the energy of these systems accurately using very compact Krylov bases.

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

𝐹 𝐾 ⁡/𝐹 𝜋 from Möbius domain-wall fermions solved on gradient-flowed HISQ ensembles

We report the results of a lattice quantum chromodynamics calculation of 𝐹 𝐾 ⁡/𝐹 𝜋 using Möbius domain-wall fermions computed on gradient-flowed 𝑁 𝑓 =2 +1 +1 highly improved staggered quark (HISQ) ensembles. The calculation is performed with five values of the pion mass ranging from 130 ≲ 𝑚 𝜋 ≲ 400 MeV , four lattice spacings of 𝑎 ∼ 0.15, 0.12, 0.09 and 0.06 fm and multiple values of the lattice volume. The interpolation/extrapolation to the physical pion and kaon mass point, the continuum, and infinite volume limits are performed with a variety of different extrapolation functions utilizing both the relevant mixed-action effective field theory expressions as well as discretization-enhanced continuum chiral perturbation theory formulas. We find that the 𝑎 ∼ 0.06 fm ensemble is helpful, but not necessary to achieve a subpercent determination of 𝐹 𝐾 ⁡/𝐹 𝜋 . We also include an estimate of the strong isospin breaking corrections and arrive at a final result of 𝐹 $\hat{K}$ + ⁡ /𝐹 $\hat{𝜋}$ + = 1.1942⁢(45) with all sources of statistical and systematic uncertainty included. This is consistent with the Flavour Lattice Averaging Group average value, providing an important benchmark for our lattice action. Combining our result with experimental measurements of the pion and kaon leptonic decays leads to a determination of |𝑉 𝑢⁢𝑠 |/|𝑉 𝑢⁢𝑑 | = 0.2311⁢(10).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

ON-OFF neuromorphic ISING machines using Fowler-Nordheim annealers

We introduce NeuroSA, a neuromorphic architecture specifically designed to ensure asymptotic convergence to the ground state of an Ising problem using a Fowler-Nordheim quantum mechanical tunneling based threshold-annealing process. The core component of NeuroSA consists of a pair of asynchronous ON-OFF neurons, which effectively map classical simulated annealing dynamics onto a network of integrate-and-fire neurons. The threshold of each ON-OFF neuron pair is adaptively adjusted by an FN annealer and the resulting spiking dynamics replicates the optimal escape mechanism and convergence of SA, particularly at low-temperatures. To validate the effectiveness of our neuromorphic Ising machine, we systematically solved benchmark combinatorial optimization problems such as MAX-CUT and Max Independent Set. Across multiple runs, NeuroSA consistently generates distribution of solutions that are concentrated around the state-of-the-art results (within 99%) or surpass the current state-of-the-art solutions for Max Independent Set benchmarks. Furthermore, NeuroSA is able to achieve these superior distributions without any graph-specific hyperparameter tuning. For practical illustration, we present results from an implementation of NeuroSA on the SpiNNaker2 platform, highlighting the feasibility of mapping our proposed architecture onto a standard neuromorphic accelerator platform.

42 ENGINEERING↗
DOE OSTI2026

Fermionic mean-field dynamics for spin systems beyond free fermions

We introduce the fermionized time-dependent Hartree–Fock (fTDHF), a real-time quantum dynamics method for spin-1/2 Hamiltonians following their mapping to fermions via the Jordan-Wigner transformation. fTDHF is formally equivalent to exact dynamics in the case of free fermions, and can efficiently handle non-local string operators arising from long-range interactions via transition matrix elements between non-orthogonal Slater determinants. We show that the fTDHF method can be implemented on a classical computer with a cost that scales polynomially with system size, and linearly with the time steps. We benchmark fTDHF against exact dynamics on three separate spin-1/2 models, representing adiabatic preparation of states with long-range correlations, disorder-driven observation of many-body localization, and particle production in the Schwinger model. For each of these systems, fTDHF is shown to reproduce the qualitative dynamics generated by the exact evolutions, while maintaining a simple physical picture due to its mean-field nature.

Dutta, Rishab↗
DOE OSTI2022

Route to Chemical Accuracy for Computational Uranium Thermochemistry

Benchmark spinor-based relativistic coupled-cluster calculations for the ionization energies of the uranium atom, the uranium monoxide molecule (UO), and the uranium dioxide molecule (UO 2 ) and for the bond dissociation energies of UO and UO 2 are reported. The accuracy of these calculations in the treatments of relativistic, electron correlation, and basis-set effects is analyzed. Furthermore, the intrinsic convergence of the computed results and the favorable comparison with the experimental values demonstrate the unique applicability of the spinor representation of quantum-chemical methods to open-shell uranium-containing atomic and molecular species with uranium oxidation states ranging from U(0) to U(V).

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

Towards efficient light emitters via computational design of molecules with inverted singlet-triplet gaps

To move toward rational design of efficient organic light emitting diodes based on the radical idea of inverted singlet-triplet gap (INVEST) systems, we propose a set of novel quantum chemical approaches, predictive but low-cost, to unveil a set of structural-property relationships. We perform a computational study of a series of substituted molecules based on a small set of known INVEST molecules. Our study demonstrates a high degree of correlation between the intramolecular charge transfer and the singlet-triplet energy gap and hints towards the use of a quantitative estimate of charge transfer to predict and modulate these energy gaps. We aim to create a database of INVEST molecules that includes accurate benchmarks of singlet-triplet energy gaps. Furthermore, we aim to link structural features and molecular properties, enabling a control knob for rational design.

42 ENGINEERING↗
DOE OSTI2022

Role of Nematicity in Controlling Spin Fluctuations and Superconducting Tc in Bulk FeSe

FeSe undergoes a transition from a tetragonal to a slightly orthorhombic phase at 90 K and becomes a superconductor below 8 K. The orthorhombic phase is sometimes called a nematic phase because quantum oscillation, neutron, and other measurements detect a significant asymmetry in x and y. How nematicity affects superconductivity has recently become a matter of intense speculation. Here, we employ an advanced ab initio Green's function description of superconductivity and show that bulk tetragonal FeSe would, in principle, superconduct with almost the same critical temperature Tc as the nematic phase. The mechanism driving the observed nematicity is not yet understood. Since the present theory underestimates it, we simulate the full nematic asymmetry by artificially enhancing the orthorhombic distortion. For benchmarking, we compare theoretical spin susceptibilities against experimentally observed data over all energies and relevant momenta. When the orthorhombic distortion is adjusted to correlate with observed nematicity in spin susceptibility, the enhanced nematicity causes spectral weight redistribution in the Fe-3dxz and Fe-dyz orbitals, but it leads to at most a 10-15% increment in Tc. This is because the dxy orbital always remains the most strongly correlated and provides most of the source of the superconducting glue. Nematicity suppresses the density of states at the Fermi level; nevertheless, Tc increases, in contradiction to both BCS theory and the theory of Bose-Einstein condensation. We show how the increase is connected to the structure of the particle-particle vertex. Our results suggest that while nematicity may be an intrinsic property of bulk FeSe, it is not the primary force driving the superconducting pairing.

bulk FeSe↗
DOE OSTI2022

$t\overline{t}t\overline{t}$ signatures through the lens of color-octet scalars

We reinterpret two recent LHC searches for events containing four top quarks ($t\overline{t}t\overline{t}$) in the context of supersymmetric models with Dirac gauginos and color-octet scalars (sgluons). We explore whether sgluon contributions to the four-top production cross section σ ($pp$ → $t\overline{t}t\overline{t}$) can accommodate an excess of four-top events recently reported by the ATLAS collaboration. We also study constraints on these models from an ATLAS search for new phenomena with high jet multiplicity and significant missing transverse energy ($E^{miss}_{T}$) sensitive to signals with four top quarks. We find that these two analyses provide complementary constraints, with the jets + $E^{miss}_{T}$ search exceeding the four-top cross section measurement in sensitivity for sgluons heavier than about 800 GeV. We ultimately find that either a scalar or a pseudoscalar sgluon can currently fit the ATLAS excess in a range of reasonable benchmark scenarios, though a pseudoscalar in minimal Dirac gaugino models is ruled out. We finally offer sensitivity projections for these analyses at the HL-LHC, mapping the 5σ discovery potential in sgluon parameter space and computing exclusion limits at 95% CL in scenarios where no excess is found.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Bootstrap embedding for interacting electrons in phonon coherent-state mean field

Here, we develop a Fermi–Bose bootstrap embedding framework for the ground state of interacting electrons coupled to a phonon mean field. The method combines bootstrap embedding for correlated electrons with a self-consistent coherent-state mean-field treatment for phonons. This method models the interacting electron–phonon problem as a system of correlated electrons traveling in a self-consistently specified potential landscape, allowing for efficient treatment of large lattice systems. Convergence of the methods for fragment size and total system size is demonstrated for the one-dimensional Hubbard–Holstein model for up to 350 sites. Finite-size scaling is performed to extrapolate to the infinite system size. Benchmarking against the density matrix renormalization group for a small 8-site system at half- and quarter-filling shows an orders-of-magnitude runtime advantage. The comparison further reveals that the method performs best in regimes dominated by localization, such as the Mott insulating phase and the strong-coupling tiny polaron regime, where the local embedding ansatz is still valid. However, due to the mean-field treatment for phonons, we find limitations of our methods in the weakly coupled delocalized region and at the Peierls transition, where quantum phonon fluctuations and long-range kinetic correlations become substantial.

Islam, Shariful [North Carolina State University, ↗
DOE OSTI2026

Auxiliary-state-facilitated phase synchronization phenomena in isolated spin systems

Extending classical synchronization to the quantum domain is of great interest both from the fundamental physics point of view and with a view toward quantum technology applications. This work characterizes phase synchronization of an effective spin-1 system, which is realized by coupling three quantum states with infinite lifetime to auxiliary excited states that have a finite lifetime. Integrating out the excited states, the effective spin-1 model features coherent and incoherent effective couplings. The following are our key findings. (i) Phase synchronization can be controlled by adjusting the phases of the couplings to the excited states. (ii) Unlike in the paradigmatic spin-1 system studied in the literature, where the dissipative couplings describe decay into the limit-cycle state, the effective spin-1 model investigated in this work is governed by a competition between dissipative decay into and out of the limit-cycle state, with the dissipative decay out of the limit-cycle state playing a critical role. (iii) We identify a parameter regime where phase synchronization of the effective spin-1 system is, in the absence of coherent effective couplings, governed entirely by the effective dissipators. The effective spin-1 model is benchmarked through comparisons with master-equation calculations for the full Hilbert space. Physical insights are gained through analytical perturbation theory calculations. In conclusion, our findings, which are expected to hold for a broad class of energy-level and coupling schemes, are examined using hyperfine states of 87 Rb as an example system.

Nonlinear optics↗
DOE OSTI2026

Hydride and Seek: Comparing Crystallographic Hydride Placement Techniques with an Open-Shell Cobalt Complex

Locating hydrides is crucial in organometallic chemistry but difficult to do accurately using X-ray diffraction. Electron diffraction has been proposed as a way to overcome this problem but has not been systematically compared to neutron diffraction and to quantum crystallography (Hirshfeld atom refinement, HAR) to test this hypothesis. Here, we present a comparative analysis of methods for a terminal cobalt hydride complex by comparing a single-crystal neutron diffraction reference structure to results from single-crystal X-ray diffraction with and without Hirshfeld atom refinement (HAR, NoSpherA2), density functional theory (DFT), and electron diffraction (3D-ED/MicroED) refined under kinematical and dynamical formalisms. Conventional X-ray diffraction gives lower precision than neutron diffraction as expected. Despite expected improvements, HAR gives systematic deviation from the neutron benchmark. Interestingly, optimized DFT equilibrium geometries are closer to the neutron value than the value from HAR. On the other hand, electron diffraction with a high-quality data set coupled with dynamical refinement localizes the hydride in difference maps and gives excellent agreement with the neutron data. Dynamical refinement is crucial, as kinematical refinement does not allow assignment of a hydride peak. This cross-modal comparison defines the conditions under which 3D-ED/MicroED delivers high-precision metal–hydride distances for this open-shell cobalt hydride.

anions↗