SEARCH · Search NASA

Results for “Quantum computing models”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

DOE OSTI2024

HamLib: A library of Hamiltonians for benchmarking quantum algorithms and hardware

In order to characterize and benchmark computational hardware, software, and algorithms, it is essential to have many problem instances on-hand. This is no less true for quantum computation, where a large collection of real-world problem instances would allow for benchmarking studies that in turn help to improve both algorithms and hardware designs. To this end, here we present a large dataset of qubit-based quantum Hamiltonians. The dataset, called HamLib (for Hamiltonian Library), is freely available online and contains problem sizes ranging from 2 to 1000 qubits. HamLib includes problem instances of the Heisenberg model, Fermi-Hubbard model, Bose-Hubbard model, molecular electronic structure, molecular vibrational structure, MaxCut, Max- k -SAT, Max- k -Cut, QMaxCut, and the traveling salesperson problem. The goals of this effort are (a) to save researchers time by eliminating the need to prepare problem instances and map them to qubit representations, (b) to allow for more thorough tests of new algorithms and hardware, and (c) to allow for reproducibility and standardization across research studies.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Capturing Surface Coverage Effects in Heterogeneous Catalysis

Adsorbate–adsorbate lateral interactions at relevant surface coverages have a significant effect on chemical kinetics, thereby influencing the activity of a heterogeneous catalyst. Coverage-dependent kinetic and thermodynamic parameters therefore must be included in studies of such complex systems to properly predict the turnover frequencies and kinetic trends. Thus, it becomes extremely important to accurately capture the strength of lateral interactions between neighboring species under realistic reaction conditions. In this Perspective, we discuss the various existing computational and experimental methods for determining adspecies coverage and configurational effects. The choice of the tools and methods employed in such studies depends on factors such as time, length scales, computational cost, the presence of solvents, and reaction conditions. The applications of each method and the respective challenges are also discussed here. As a result, we discuss the recent developments and future of the state-of-the-art for inclusion of surface coverage and configuration into a holistic picture for accurate predictions of catalytic behavior.

09 BIOMASS FUELS↗

DOE OSTI2024

U(1) fields from qubits: An approach via D-theory algebra

A new quantum link microstructure was proposed for the lattice quantum chromodynamics (QCD) Hamiltonian, replacing the Wilson gauge links with a bilinear of fermionic qubits, later generalized to D-theory. This formalism provides a general framework for building lattice field theory algorithms for quantum computing. We focus mostly on the simplest case of a quantum rotor for a single compact U(1) field. We also make some progress for non-Abelian setups, making it clear that the ideas developed in the U(1) case extend to other groups. These in turn are building blocks for 1 + 0 -dimensional ( 1 + 0 -D) matrix models, 1 + 1 -D sigma models and non-Abelian gauge theories in 2 + 1 and 3 + 1 dimensions. By introducing multiple flavors for the U(1) field, where the flavor symmetry is gauged, we can efficiently approach the infinite-dimensional Hilbert space of the quantum O(2) rotor with increasing flavors. The emphasis of the method is on preserving the symplectic algebra exchanging fermionic qubits by sigma matrices (or hard bosons) and developing a formal strategy capable of generalization to a SU ( 3 ) field for lattice QCD and other non-Abelian 1 + 1 -D sigma models or 3 + 1 -D gauge theories. For U(1), we discuss briefly the qubit algorithms for the study of the discrete 1 + 1 -D sine-Gordon equation. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

DOE OSTI2024

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS2003

Interfacial Reaction Studies Using ONIOM

In this report, we focus on the calculations of the energetics and chemical kinetics of heterogeneous reactions for Organometallic vapor phase epitaxy (OMVPE). The work described in this report builds upon our own previous thermochemical and chemical kinetics studies. The first of these articles refers to the prediction of thermochemical properties, and the latter one deals with the prediction of rate constants for gaseous homolytic dissociation reactions. The calculations of this investigation are at the microscopic level. The systems chosen consisted of a gallium nitride (GaN) substrate, and molecular nitrogen (N2) and ammonia (NH3) as adsorbants. The energetics for the adsorption and the adsorbant dissociation processes were estimated, and reaction rate constants for the dissociation reactions of free and adsorbed molecules were predicted. The energetics for substrate decomposition was also computed. The ONIOM method, implemented in the Gaussian98 program, was used to perform the calculations. This approach has been selected since it allows dividing the system into two layers that can be treated at different levels of accuracy. The atoms of the substrate were modeled using molecular mechanics6 with universal force fields, whereas the adsorbed molecules were approximated using quantum mechanics, based on density functional theory methods with B3LYP functionals and 6-311G(d,p) basis sets. Calculations for the substrate were performed in slabs of several unit cells in each direction. The N2 and NH3 adsorbates were attached to a central location at the Ga-lined surface.

Cardelino, Beatriz H.↗

DOE OSTI2025

Multi-scale computational screening and mechanistic insights of cyclic amines as solvents for improved lignocellulosic biomass processing

A computational screening workflow for the efficient deconstruction of cellulose, lignin and hemicellulose fractions of lignocellulosic biomass using cyclic amines as solvents. Lignocellulosic biomass is a promising feedstock for production of affordable fuels and chemicals from renewable resources. Effective solubilization and subsequent deconstruction of its cellulose, hemicellulose, and lignin fractions is essential for the viability of future biorefineries. This study used quantum chemistry-based equilibrium thermodynamics methods to evaluate the potential of 650 cyclic amines to solubilize cellulose, hemicellulose, and lignin. The activity coefficients of solvent - biopolymer interactions were predicted using the COSMO-RS (COnductor-like Screening MOdel for Real Solvents) method and used to identify cyclic amines that can efficiently dissolve and extract selective fractions of biopolymers during biomass pretreatment. Among the 650 cyclic amines, 1-piperazineethanmaine was predicted to be an effective solvent for extracting all three polymers and was experimentally shown to achieve the highest lignin removal (97.1%). Non-covalent interaction, reduced density gradient and quantum chemical calculations were performed to elucidate the dissolution mechanism of lignin, cellulose and hemicellulose and gain further molecular level insights into the interactions between the cyclic amines and biomass polymers that promote efficient solubilization and extraction. These analyses indicated that 1-piperazineethanmaine and 1-methylimidazole make noncovalent van der Waals, electrostatic interactions and hydrogen bonding with lignin, leading to enhanced lignin removal, while the strong intramolecular hydrogen bonding interactions in cellulose and hemicellulose result in weaker solvent-biopolymer interactions. Overall, the computational approach provided an efficient method for identifying cyclic amines tailored for optimal biomass pretreatment and resulted in the identification of a potential new class of solvents for effective biomass pretreatment.

Kumar, Nikhil↗

DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Simulating spin dynamics of supersolid states in a quantum Ising magnet

Motivated by a recent experimental study on the quantum Ising magnet K 2 Co(SeO 3 ) 2 that presented spectroscopic evidence of zero-field supersolidity (Chen et al., arXiv:2402.15869), we simulate the excitation spectrum of the corresponding microscopic XXZ model for the compound, using the recently developed excitation ansatz for infinite projected entangled-pair states. Here, we map out the ground state phase diagram and compute the dynamical spin structure factors across a range of magnetic field strengths, focusing especially on the two supersolid phases found near zero and saturation fields. Our simulated excitation spectra for the zero-field supersolid “Y” phase are in excellent agreement with the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2009

Development of a GCR Event-based Risk Model

A goal at NASA is to develop event-based systems biology models of space radiation risks that will replace the current dose-based empirical models. Complex and varied biochemical signaling processes transmit the initial DNA and oxidative damage from space radiation into cellular and tissue responses. Mis-repaired damage or aberrant signals can lead to genomic instability, persistent oxidative stress or inflammation, which are causative of cancer and CNS risks. Protective signaling through adaptive responses or cell repopulation is also possible. We are developing a computational simulation approach to galactic cosmic ray (GCR) effects that is based on biological events rather than average quantities such as dose, fluence, or dose equivalent. The goal of the GCR Event-based Risk Model (GERMcode) is to provide a simulation tool to describe and integrate physical and biological events into stochastic models of space radiation risks. We used the quantum multiple scattering model of heavy ion fragmentation (QMSFRG) and well known energy loss processes to develop a stochastic Monte-Carlo based model of GCR transport in spacecraft shielding and tissue. We validated the accuracy of the model by comparing to physical data from the NASA Space Radiation Laboratory (NSRL). Our simulation approach allows us to time-tag each GCR proton or heavy ion interaction in tissue including correlated secondary ions often of high multiplicity. Conventional space radiation risk assessment employs average quantities, and assumes linearity and additivity of responses over the complete range of GCR charge and energies. To investigate possible deviations from these assumptions, we studied several biological response pathway models of varying induction and relaxation times including the ATM, TGF -Smad, and WNT signaling pathways. We then considered small volumes of interacting cells and the time-dependent biophysical events that the GCR would produce within these tissue volumes to estimate how GCR event rates mapped to biological signaling induction and relaxation times. We considered several hypotheses related to signaling and cancer risk, and then performed simulations for conditions where aberrant or adaptive signaling would occur on long-duration space mission. Our results do not support the conventional assumptions of dose, linearity and additivity. A discussion on how event-based systems biology models, which focus on biological signaling as the mechanism to propagate damage or adaptation, can be further developed for cancer and CNS space radiation risk projections is given.

Cucinotta, Francis A.↗

DOE OSTI2026

LDRD FY25 Program Overview

As Lawrence Livermore National Laboratory’s (LLNL’s) Laboratory Directed Research and Development (LDRD) program enters its fifth decade of leading-edge research and development, its impact and importance have never been stronger. The program continues to advance strategic investments in pioneering science, technology, and engineering, ensuring LLNL will be ready to deliver on our mission as it evolves over the coming decades. Investing in LDRD research, and the people who perform this critical work, gives LLNL the ability to sustain our role as a leader in the Department of Energy and National Nuclear Security Administration enterprise. The LDRD program enables high-risk, high-payoff research that anticipates emerging threats and future mission needs. By nurturing the ingenuity of the Lab’s greatest asset, its people, LDRD funding advances not only our research but also grows and nurtures our workforce: engaging future innovators with student mentoring, challenging postdoctoral researchers to apply their skills to support national security, and strengthening the leadership skills of early career staff. This annual report documents how LDRD investments advance LLNL’s science, technology, and engineering across our mission space. To assess LDRD’s impact we track both short and long-term metrics such as peer-reviewed publications, number of students, or professional fellows. In addition to reviewing these metrics, I encourage you to delve deeper into the breadth of science and technology that illustrate the strategic value of this research portfolio. For instance, a recent exploratory research project used advanced manufacturing to construct miniaturized three-dimensional ion traps for a quantum computer with reduced quantum error rates to enable applications that address national security missions and support basic science. Another project has delved into studying detonation by examining deflagration to enhance the safety and security of the nuclear weapons stockpile. LDRD researchers are also deploying AI agents on two of the world’s most powerful supercomputers to automate and accelerate inertial confinement fusion experiments. Other teams are delivering more accurate optical constants to enable improved validation for aluminum to advance atomic and molecular physics models. LDRD-driven discoveries of how metals deform under extreme conditions strengthen our ability to model and design materials for demanding national security environments. National security challenges are increasingly complex and continuously evolving. LDRD focuses our most innovative science and technology on these challenges, ensuring the Laboratory is developing creative, forward-leaning solutions for our nation and the world. The following pages feature highlights of published scientific advances, patents, and honors that stem from LDRD investments. As you read this report, I hope you will understand how these investments position the Laboratory, and our partners, to meet the demands of the decades ahead.

36 MATERIALS SCIENCE↗

NASA NTRSDate not supplied

Scalar Atomic Defect-Based Solid-State Self-calibrating Magnetometer (3SM) for Space Plasma Analysis

The Earth’s magnetosphere is a system of multiple, co-located particle populations interacting via plasma waves. Things to understand are driving processes, radiation belt and ring current issues, auroral physics, internal plasma processes, and magnetosphere-ionosphere mapping issues. These plasma-physics processes enable the Earth’s magnetosphere to evolve in response to temporal changes in the solar wind and they underlie the phenomena of space weather, which impacts spacecraft systems, astronauts, radio communications, and ground based electric-power grids. The objective of 3SM is to measure magnetic field strength and to calibrate the Vector Magnetometer (VM) device to maintain absolute accuracy during a mission. The 3SM features make the instrument preferably suited not only for the traditional role of scalar magnetometers as absolute references for the calibration of the on-board vector instruments, but also for extended operational capacities, such as higher frequency scalar measurements (of potential interest for magnetosphere studies for the low frequency part of the spectrum) or autonomous scalar / vector operations. Diamond has been the solid-state platform of choice for quantum device technologies for some time, however, it suffers from difficulties such as scalability, integration, and cost. While the diamond platform is very useful for quantum technologies, further development is needed to make it viable. The material platform of choice for NASA Glenn’s Quantum Sensing And Spin Physics (Q-SASP) is silicon carbide (SiC). This is due to the much higher industry development of the SiC material platform for high-power and high-temperature electronics. It leverages both the decades-long SiC development expertise and infrastructure at NASA Glenn and its growing capabilities in quantum metrology. To make SiC devices usable for quantum technologies such as quantum sources, a much deeper understanding of defects is needed. Q-SASP is developing quantum metrology capabilities to evaluate the energy structure, defect formation energy, band structure augmentation, generation/recombination rates, and limits of dipole-dipole coupling in non-metal implanted SiC devices. This can be achieved by analysis of zero-field splitting, low-field resonance, and singlet-triplet mixing through various forms of Electrically Detectable Magnetic Resonance (EDMR) and Near-Zero Field Magnetic Resonance (NZFMR) spectroscopy. This work will discuss recent system developments, device developments, computational modeling, and spectroscopy results and analysis of defects created by non-metal implantations in SiC devices. The defect formation energies of V Si ,V C , V C V Si , N C V Si , N Si , N C in 4H-SiC are previously reported values in other research [2]-[3]. The defect formation energies of PSi and PC were calculated in GPAW [fig 1A]. The basic underlying mechanism of the zero-field phenomenon is the mixing of singlet and triplet states [4]-[5]. In most spin-dependent transport, two electron spins are involved, and thus one must consider each of their interactions with the field. We investigated the electronic and magnetic properties of 4H-SiC and 6H-SiC. The defect formation energy helps us determine what types of defects we are observing in the SiC EDMR experiment. They have very low formation energy (it is negative). The phosphorus substitution in 4H-SiC is a very stable defect. The band diagrams provide us with vital information about how the electronic properties of SiC (such as band gap) change as we add non-metal defects. The zero-field splitting parameters allow us to study the inflection point in the NZFMR [fig 1B]. We clearly observed zero-field splitting. We also noted that the zero-field splitting remained constant with changing bias. We aspect it zero-field splitting to remain constant while the hyperfine and exchange interaction perturbations shift under the influence of an external magnetic field. This is the essence of quantum magnetometry and self-calibration.

space plasma↗

DOE OSTI2024

Colossal Strain Tuning of Ferroelectric Transitions in KNbO 3 Thin Films

Strong coupling between polarization ( P ) and strain (ɛ) in ferroelectric complex oxides offers unique opportunities to dramatically tune their properties. Here colossal strain tuning of ferroelectricity in epitaxial KNbO 3 thin films grown by sub‐oxide molecular beam epitaxy is demonstrated. While bulk KNbO 3 exhibits three ferroelectric transitions and a Curie temperature ( T c ) of ≈676 K, phase‐field modeling predicts that a biaxial strain of as little as −0.6% pushes its T c > 975 K, its decomposition temperature in air, and for −1.4% strain, to T c > 1325 K, its melting point. Furthermore, a strain of −1.5% can stabilize a single phase throughout the entire temperature range of its stability. A combination of temperature‐dependent second harmonic generation measurements, synchrotron‐based X‐ray reciprocal space mapping, ferroelectric measurements, and transmission electron microscopy reveal a single tetragonal phase from 10 K to 975 K, an enhancement of ≈46% in the tetragonal phase remanent polarization ( P r ), and a ≈200% enhancement in its optical second harmonic generation coefficients over bulk values. These properties in a lead‐free system, but with properties comparable or superior to lead‐based systems, make it an attractive candidate for applications ranging from high‐temperature ferroelectric memory to cryogenic temperature quantum computing.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Superconductor to exciton condensate transition in a model copper-oxide material

Abstract Superconductivity and exciton condensation are fundamental phenomena in condensed matter physics, associated with the condensation of electron–electron and electron–hole pairs, respectively, into coherent quantum states. In this study, we present evidence of a superconductor to exciton condensate transition within the context of the three-band Hubbard model of copper-oxide-like materials. As the electron–electron repulsion increases, the superconducting phase is superseded by exciton condensation. In support of theoretical predictions—not yet realized experimentally—we observe the coexistence of the two condensates in the vicinity of the transition where the quantum states become a superposition of electron–electron and electron–hole condensates. Coexistence is rigorously computed from large eigenvalues and their eigenvectors in both the two-electron reduced density matrix (2-RDM) and the particle-hole RDM, which we obtain from a direct variational ground-state energy minimization with respect to the 2-RDM by semidefinite programming. We further discern that adjacent d orbitals and intervening p orbitals facilitate electron–electron pairing between copper orbitals, thereby supporting the superexchange mechanism for superconductivity. These observations suggest the feasibility of witnessing a superconductor to exciton condensate transition in copper-oxide analogs, bearing significant implications for identifying materials conducive to efficient transport processes.

Schouten, Anna O. (ORCID:0000000264255219)↗

NASA NTRS1997

Quantum Chemical Calculations

The current methods of quantum chemical calculations will be reviewed. The accent will be on the accuracy that can be achieved with these methods. The basis set requirements and computer resources for the various methods will be discussed. The utility of the methods will be illustrated with some examples, which include the calculation of accurate bond energies for SiF$_n$ and SiF$_n^+$ and the modeling of chemical data storage.

Bauschlicher, Charles W.↗

DOE OSTI2025

Impact of dynamics, entanglement and Markovian noise on the fidelity of few-qubit digital quantum simulation

Quantum algorithms have been proposed to accelerate the simulation of the chaotic dynamical systems that are ubiquitous in the physics of plasmas. Quantum computers without error correction might even use noise to their advantage to calculate the Lyapunov exponent by measuring the Loschmidt echo fidelity decay rate. For the first time, digital Hamiltonian simulations of the quantum sawtooth map, performed on the IBM-Q quantum hardware platform, show that the fidelity decay rate of a digital quantum simulation increases during the transition from dynamical localization to chaotic diffusion in the map. The observed error per CNOT gate increases by $1.5{\times }$ as the dynamics varies from localized to diffusive, while only changing the phases of virtual RZ gates and keeping the overall gate count constant. A gate-based Lindblad noise model that captures the effective change in relaxation and dephasing errors during gate operation qualitatively explains the effect of dynamics on fidelity as being due to the localization and entanglement of the states created. Specifically, highly delocalized states that are entangled with random phases show an increased sensitivity to dephasing and, on average, a similar sensitivity to relaxation as localized states. In contrast, delocalized unentangled states show an increased sensitivity to dephasing but a lower sensitivity to relaxation. This gate-based Lindblad model is shown to be a useful benchmarking tool by estimating the effective Lindblad coherence times during CNOT gates and finding a consistent $2\unicode{x2013}3{\times }$ shorter $T_2$ time than reported for idle qubits. Thus, the interplay of the dynamics of a simulation with the noise processes that are active can strongly influence the overall fidelity decay rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Iodine recombination in xenon solvent: Clusters in the gas to liquid-like state transition

Supercritical fluids (SCFs) have attracted significant attention as solvents for chemical reactions due to their unique properties, such as high diffusivity, low viscosity, and tunable solvation properties. These properties profoundly influence reaction kinetics and are often attributed to the formation of molecular clusters within SCFs. To study the effect of supercritical solvent on chemical reactivity and dynamics of reactions, one needs to understand the dynamics of clusters in supercritical fluid. Extensive experiments on the photodissociation and recombination of iodine in supercritical fluids served as a model system for understanding these effects. Experimental studies have been complemented by theoretical and computational investigations, which mostly employ Monte Carlo or empirical molecular dynamics simulations. However, computational studies using non-reactive force fields and ab initio approaches present challenges in capturing reactive processes at larger scales within supercritical fluids. Here, in this work, we developed the ReaxFF parameters by training against quantum mechanics data. ReaxFF reactive force field based molecular dynamics simulations were performed, studying the dynamics of a xenon solvent and cage effect at different thermodynamic conditions for the iodine recombination reaction. We show that the conditions near the critical point are the optimal conditions to study the cage effect. We show that the average lifetime of xenon clusters ranging between 5 and 11 ps is comparable to iodine geminate recombination. Our simulation results of iodine recombination in xenon solvent demonstrate the higher probability of iodine molecule formation in the presence of xenon clusters. Finally, we show that the supercritical condition exhibits the highest recombination rate for iodine atoms.

Cage effect↗

DOE OSTI2024

Modeling athermal phonons in novel materials using the G4CMP simulation toolkit

Understanding phonon and charge propagation in superconducting devices plays an important role in both performing low-threshold dark matter searches and limiting correlated errors in superconducting qubits. The Geant4 Condensed Matter Physics (G4CMP) package, originally developed for the Cryogenic Dark Matter Search (CDMS) experiment, models charge and phonon transport within silicon and germanium detectors and has been validated by experimental measurements of phonon caustics, mean charge-carrier drift velocities, and heat pulse propagation times. Here, in this work, we present a concise framework for expanding the capabilities for phonon transport to a number of other novel substrate materials of interest to the dark matter and quantum computing communities, including sapphire (Al 2 O 3 ), gallium arsenide (GaAs), lithium fluoride (LiF), calcium tungstate (CaWO 4 ), and calcium fluoride (CaF 2 ). We demonstrate the use of this framework in generating phonon transport properties of these materials and compare these properties with experimentally-determined values where available.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2026

The critical role of intrinsic defects and many-body interactions on the stability of MnBi2Te4

Intrinsic antisite defects pose a major challenge to understanding and predicting the exotic properties of the layered topological magnetic insulator MnBi2Te4 (MBT). In this work, we study the origin of the abundance of intrinsic defects in MBT, including many-body defect–defect interactions and many-body electronic correlations. Until now, ab initio methods have struggled to explain thermodynamic stability and properties influenced by defect behavior in MBT. We model native Mn–Bi antisite defects in MBT at finite temperatures using a cluster expansion that includes defect–defect interactions. To overcome the limitations of conventional density functional theory (DFT), we introduce a hybrid approach that incorporates high-accuracy quantum Monte Carlo (QMC) calculations, introducing missing correlations. This strategy allows for accurate estimation of defect energetics and finite-temperature properties. We compute the configurational free energy, defect concentration, and configurational heat capacity, revealing a second-order order–disorder phase transition near the experimental synthesis temperature. Our study provides the first theoretical insight into the thermodynamics of intrinsic defects in MBT. The negative free energy relative to pristine MBT at synthesis temperatures indicates that Mn–Bi antisite formation is thermodynamically spontaneous. We also present a broadly applicable general framework for correcting low-level theoretical theories using highly accurate many-body corrections from QMC.

Ghaffar, Abdul [ORNL] (ORCID:0000000241190168)↗