Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupled dynamics”

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 181 records · Page 10

High-throughput spin-bath characterization of spin defects in semiconductors

Detailed knowledge of the local environments of spin defects in semiconductors, such as nitrogenvacancy (NV) centers in diamond or divacancies in silicon carbide, is crucial for optimizing control and entanglement protocols in quantum sensing and information applications. However, at present a direct experimental characterization of individual defect environments is not scalable, as conventional spin-bath measurements are time consuming and difficult to automate. Achieving high-throughput characterization requires short experiments to probe the spin bath. However, with fewer and noisier measurements, the inverse problem of recovering spin-bath properties from measured data becomes ill posed, with multiple spin baths having a high likelihood of yielding the same data. In this work, we present a set of computational tools to resolve the ill-posed inverse problem of recovering the atomic positions and hyperfine couplings of random nuclei surrounding spin defects from sparse, noisy experimental coherence data, which can be obtained in hours. Here, we use a trans-dimensional Bayesian approach that incorporates ab initio data to yield full posterior distributions over nuclear spin environments, enabling robust recovery from limited data. We also provide practical tools and guidelines to determine the limits of detectability for hyperfine couplings under specific dynamical decoupling sequences and sampling conditions. In addition, we demonstrate how the tools developed here, in combination with ab initio simulations of spin baths, can guide the design of efficient experimental protocols for application-specific high-throughput screening. To showcase the utility of our approach, we apply it to design fast dynamical decoupling experiments to characterize the spin baths often individual NV centers in diamond. While the primary focus is on accelerating spin-bath characterization of spin defects, this Bayesian approach also lays the foundation for digital-twin studies of spin defects, where a virtual model of the spin-defect system evolves in real time with ongoing experimental measurements. Together, the set of tools we designed and applied paves the way for scalable deployment of spin defects in semiconductors for quantum sensing and information applications.

Bayesian methods↗

Interplay of 𝑑- and 𝑝-states in RbTi3⁢Bi5 and CsTi3⁢Bi5 flat-band kagome metals

Shifting the Fermi level of the celebrated 𝐴⁢𝑀3⁢𝑋5 (135) compounds into proximity of flat bands strongly enhances electronic correlations and severely affects the formation of density waves and superconductivity. Our broadband infrared spectroscopy measurements of RbTi3⁢Bi5 and CsTi3⁢Bi5 combined with density-functional band-structure calculations reveal that the correlated Ti 𝑑-states are intricately coupled with the Bi 𝑝-states that form a tilted Dirac crossing. Electron-phonon coupling manifests itself in the strong damping of itinerant carriers and in the anomalous shape of the phonon line in RbTi3⁢Bi5. An anomaly in these spectral features around 150 K can be paralleled to the onset of nematicity detected by low-temperature probes. Our findings show that the materials with low band filling open unexplored directions in the physics of kagome metals and involve electronic states of different nature strongly coupled with lattice dynamics.

Wenzel, Matthew [ORNL]↗

Robust resonator-assisted ZZ cancellation in superconducting quantum processors

Strong qubit interactions are essential for faster two-qubit gates, but they often come with undesirable ZZ interactions that limit gate fidelity. Existing methods to mitigate these interactions, such as flux-tunable couplers, can introduce additional noise and complexity. In contrast, our work presents a simpler, more robust approach using a driven resonator to cancel the static ZZ interaction between qubits [1]. The experiment was performed on a revised 9-qubit quantum processing unit from Rigetti, developed in collaboration with SQMS scientists. We validate the resonator-induced-phase (RIP) interaction, where an off-resonant drive on the resonator dynamically cancels ZZ coupling. This marks an important step toward high gate fidelities. We also explore the entangling gates enabled by this coupling scheme, focusing on minimizing gate duration and qubit decoherence while maintaining effective ZZ cancellation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Development of an Improved RELAP5-3D Model for the High Temperature Test Facility

High-temperature gas-cooled reactors (HTGRs) are rapidly approaching deployment. Confidence in transient analysis of these systems for design, optimization, and licensing calculations requires modeling and simulation tools that have been validated against data relevant to HTGR conditions. The High Temperature Test Facility (HTTF) is an integral effects thermal hydraulics test facility for prismatic HTGRs. In spring and summer of 2019, HTTF was used for a series of experiments that now serve as the basis for the OECD/NEA Thermal Hydraulic Code Validation Benchmark for High Temperature Gas-Cooled Reactors using HTTF Data (HTGR T/H Benchmark). This benchmark contains problems for systems code, computational fluid dynamics (CFD), and coupled systems code/CFD modeling representing lower plenum mixing and both the depressurized and pressurized conduction cooldown (DCC and PCC respectively) transients. Benchmark problems include exercises for code-to-code and code-to-data comparisons as well as an exercise for error scaling between HTTF and the Modular High Temperature Gas-Cooled Reactor, which serves as the basis for the HTTF design. Previous analysis as part of the HTGR T/H benchmark used a RELAP5-3D model developed at Idaho National Laboratory (INL) and demonstrated an ability to reproduce trends in the measured data but difficulties reproducing experimental values within their uncertainty. These difficulties were largely attributed to assumptions made during the development of the initial RELAP5-3D model, which predated the HTTF experiments. A significant cause of difficulty reproducing the measured temperatures may be the radial nodalization of the previous RELAP5-3D model. The new model provides a finer nodalization to assess the impact of radial nodalization and allows for asymmetric heating within the core, which was a feature of multiple HTTF experiments. In this paper, we present the new RELAP5-3D model of HTTF. In addition to describing the new model, this paper compares the new and old models and provides results for a full-power steady state, a DCC, and a PCC in HTTF. These analyses are based on the code-to-code comparison exercises for the DCC and PCC problems of the HTGR T/H benchmark. We present the results of these exercises from the new model and compare them to the results of the old model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Modeling Radiolysis and Chemical Reactions during Dry Storage of Aluminum-clad Spent Nuclear Fuel

After aluminum-clad spent nuclear fuel (ASNF) is removed from the reactor, it is initially stored in spent fuel pools, which are specially designed water-filled basins that provide temporary cooling to reduce the temperature of the fuel assemblies and provide radiation shielding. ASNF continues to generate heat due to the radioactive decay of elements within the fuel, which persists for many years post-shutdown as the residual radioactive products decay into more stable elements. During the wet storage period, an oxyhydroxide layer composed of boehmite/bayerite forms on the surfaces of the aluminum cladding from exposure to water in the pools. Road-ready packaging for long-term disposition of the ASNF involves dry storage in helium backfilled DOE standard canisters (DSCs). When the ASNF is removed from water storage and dried, most of the water is removed, but some physisorbed and chemisorbed water remains in the oxyhydroxide layers. This residual water can produce hydrogen when exposed to radiation from the ASNF during dry storage. Predicting hydrogen accumulation over time in the DSCs is critical for long-term storage considerations. Previous modeling efforts have developed coupled computational fluid dynamics (CFD)-chemical models to simulate temperature, pressure, and gas phase concentrations within the DSCs. These models use the thermal field predicted by CFD as input to a radiolysis model for the gas phase and the surface oxyhydroxide layer chemistry. Given the long storage period of the DSCs and the impracticality of long-term experiments, a simulation-based approach is necessary to assess chemical evolution within the canisters. This study advances the development of a modeling framework designed to simulate the chemical evolution of spent fuel canisters. Both thermal and radiation-driven reactions are considered, with radiation kinetics quantified using G-values. Sensitivity analysis identifies key parameters influencing species composition. Reaction pathway diagrams offer insight into dominant species formation routes, enabling more effective comparisons between model predictions and experimental observations, particularly regarding the production of hydrogen. Results show that the model predicts significant hydrogen gas production with minimal oxygen generation, primarily due to hydrogen formation via boehmite pathways. These findings underscore the importance of accurately characterizing surface-bound species and radiolysis kinetics. A deeper understanding of these mechanisms is critical for evaluating the long-term safety of nuclear waste storage.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗

Multiscale Nuclear-Electronic Orbital Quantum Dynamics in Complex Environments

Many renewable energy conversion processes rely on the movement of protons as well as electrons through either electrocatalysis or photoexcitation. The simulation of such processes requires a quantum mechanical description of coupled nuclear-electronic dynamics in a solvent or heterogeneous chemical environment. The overall objective of this project is the development of theoretical and computational capabilities for simulating nuclear-electronic quantum dynamics in complex environments and the creation of high-performance, open-source software. This multiscale framework will enable simulations of the real-time dynamics of nonequilibrium excited state proton-coupled electron transfer, quantum decoherence, vibronic energy transfer, and ultrafast radiolysis, as well as their associated time-resolved multidimensional spectroscopies. An important outcome of this project will be a sustainable, reusable, and interoperable open-source software ecosystem. This software will be designed for emerging exascale and future national leadership computers. Another key outcome will be a multiscale quantum dynamics method and software enabling simulations of nonequilibrium nuclear-electronic quantum dynamics in complex environments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

pH Regulates Ion Dynamics in Carboxylated Mixed Conductors

Coupled ionic and electronic transport underpins processes as diverse as electrochemical energy conversion, biological signaling, and soft adaptive electronics. Yet, how chemical environments such as pH modulate this coupling at the molecular scale remains poorly understood. Here, we show that the protonation state of carboxylated polythiophenes provides precise chemical control over ion dynamics, doping efficiency, solvent uptake, and mechanical response. Using a suite of multimodal operando techniques, supported by simulations, we reveal that pH dictates the balance of cation/anion uptake during electrochemical doping. Mapping across pH uncovers a quasi-nonswelling regime (≈pH 3–3.5) where charge compensation proceeds with minimal volumetric change yet pronounced stiffening. These findings establish molecular acidity as a general strategy to program ionic preference and mechanical stability, offering design principles for pH-responsive mixed conductors and soft electronic materials that couple ionic, electronic, and mechanical functionality.

PH↗

Exploring the effect of ELM and code-coupling frequencies on plasma and material modeling of dynamic recycling in divertors

Abstract Integrated modeling of plasma-surface interactions provides a comprehensive and self-consistent description of the system, moving the field closer to developing predictive and design capabilities for plasma facing components. One such workflow, including descriptions for the scrape-off-layer plasma, ion-surface interactions and the sub-surface evolution, was previously used to address steady-state scenarios and has recently been extended to incorporate time-dependence and two-way information flow. The new model can address dynamic recycling in transient scenarios, such as the application presented in this paper: the evolution of W samples pre-damaged by helium and exposed to ELMy H-mode plasmas in the DIII-D DiMES. A first set of simulations explored the effect of ELM frequency. This study was discussed in detail in this conference’s proceedings and is summarized here. The 2nd set of simulations, which is the focus of this paper, explores the effect of code-coupling frequency. These simulations include initial SOLPS solutions converged to the inter-ELM state, ion impact energy ( E in ) and angles ( A in ) calculated by hPIC2, and an improved heat transfer description in Xolotl. The model predicts increases in particle fluxes and decreases in heat fluxes by 10%–20% with the coupling time-step. Compared with the first set of simulations, the less shallow impact angle leads to smaller reflection rates and significant D implantation. The higher fraction of implanted flux (and deeper), in particular during ELMs, increases the accumulated D content in the W near-surface region. Future expansion of the workflow includes coupling to hPIC2 and GITR to ensure accurate descriptions of E in and A in , and W impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamical Generation of Higher-Order Spin-Orbit Coupling, Topology, and Persistent Spin Texture in Light-Irradiated Altermagnets

Altermagnets have been identified as the third category of magnetic materials, exhibiting momentum-dependent spin splitting characterized by even powers of momentum. Here, in this Letter, we show that when subjected to elliptically polarized light, these materials serve as an exemplary framework for the dynamic generation of topological bands featuring higher-order spin-orbit coupling (SOC). Notably, while the generated Zeeman field remains invariant to the particular altermagnetic ordering, the induced higher-order SOCs are related to the magnitude and symmetry of the altermagnetic order. Specifically, we show that an altermagnet exhibiting 𝐤 2⁢𝑛 -spin splitting can generate odd-in-𝐤 spin-orbit couplings up to 𝐤 2⁢𝑛−1 . In the limit of circularly polarized light, the only correction is 𝑘 2⁢𝑛−1 , with all lower-order contributions being nullified. Interestingly, light-induced SOCs significantly impact the low-energy band topology, where their Chern numbers change by Δ⁢𝐶 = ±1, 2, 3 for 𝑑, 𝑔, 𝑖-wave altermagnets. Finally, we find a critical field in which a persistent spin texture is realized, a highly desirable state with predicted infinite spin lifetime. Our Letter showcases light as a powerful, controllable tool for engineering complex and exciting phenomena in altermagnets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Reconceptualizing the Ir III Role in Metallaphotoredox Catalysis: From Strong Photooxidant to Potent Energy Donor

Dual Ir III /L n Ni II metallaphotoredox catalyzed C(sp 3 )–C(sp 2 ) cross-coupling reactions are widely assumed to proceed by photoinduced single electron transfer steps due to the highly oxidizing Ir III * excited state (Ir III = [Ir(dF(CF 3 )ppy) 2 (dtbbpy)] + [PF 6 ] – ; dF(CF 3 )ppy = 2-(2,4-difluorophenyl)-5-(trifluoromethyl)pyridine; L n = dtbbpy = 4,4'-di-tert-butyl-2,2'-bipyridine). Using time-resolved absorption and emission spectroscopy, we reveal that energy transfer between Ir III * and various LnNi II precatalysts and intermediates with k q ≥ 10 8 M –1 s –1 also drives catalysis. Specifically, the excited states of L n Ni II dihalide precatalysts/organometallic intermediates accessible by energy transfer appear to drive bond homolysis, halogen radical elimination, and reductive elimination reactions that facilitate formation of cross-coupled products. Energy transfer dynamics consequently circumvent the need for photoinduced electron transfer, thereby extending substrate scopes to coupling partners that cannot be oxidized by Ir III *. Within a cross-electrophile coupling model reaction between 4-bromobenzotrifluoride and bromocyclohexane, energy transfer activates the L n Ni II precatalyst at early reaction times before nucleophilic reductants are present. In the absence of Ir III , direct excitation of L n Ni II (Br) 2 also activates the precatalyst to form a L n Ni II (Br)(Aryl) intermediate. To compare energy transfer and electron transfer kinetics, we determined rate constants for reductive quenching by Br – (k SET = 4.1 × 10 8 M –1 s –1 ) and for the subsequent electron transfer from reduced Ir III•– to L n Ni II (Br) 2 (k SET = 4.1 × 10 7 M –1 s –1 ) using Stern-Volmer analysis and pulse radiolysis, respectively. Energy transfer rate constants are competitive with the electron transfer rate constants and energy transfer is a parallel pathway within metallaphotoredox catalysis. Exploiting the energy transfer mechanism, we demonstrate highly selective cross-electrophile coupling between 4-chlorobenzotrifluoride and bromocyclohexane to form exclusively cross-coupled product. Here, with alkyl-trifluoroborate nucleophiles that do not reductively quench IrIII* emission, transmetalation with L n Ni II (Br/Cl)(Aryl) followed by energy transfer also drives excited state reductive elimination to form C(sp 3 )–C(sp 2 ) cross-coupled product. Similarly, energy transfer rather than Ni II oxidation drives C(sp 2 )–OR reductive elimination, despite the strongly oxidizing ability of Ir III *. In total, these reactions demonstrate energy transfer processes from Ir III * to L n Ni II in metallaphotoredox catalysis that can unlock alternative reactive pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamical Inference and 3D Imaging of Magnetized Dusty Plasmas

This is the final technical report for this DOE award. The motion of dust particles in a laboratory plasma has been studied for nearly 30 years and has led to the discovery of strongly-coupled crystalline structures and nonequilibrium dynamics governed by gravitational, hydrodynamic, and electrostatic forces. Many of the inter-particle forces are complex; they can be non-reciprocal, and non-additive. Deciphering the mélange of particle interactions has been a challenging task, nevertheless, magnetized dusty plasmas remains an open challenge with limited understanding. With applications ranging from magnetically-confined plasmas for fusion to near-surface planetary environments, this research field is ripe for well-controlled, laboratory experiments and new theoretical tools. In collaboration with the MDPX experiment at Auburn University, this proposal aims to tease apart both the known and unknown forces that drive dusty plasmas in magnetized environments using high-resolution, three-dimensional imaging and particle tracking coupled with modern dynamical inference and machine learning techniques.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling graphene sheet growth and dynamical matrix calculations using molecular dynamics

Molecular dynamics (MD) has been an incredibly useful tool to model physical processes that were synthesized experimentally but not fully understood. MD, through the use of semi-empirical inter-atomic potentials, has allowed understanding of different physical processes in materials science. Yet as well as providing useful insights into materials science, molecular dynamics has a wider range of usability. In this report, I will be detailing how MD can be used to study graphene formation from a carbon liquid which requires high temperatures and pressures. Beyond this, I will describe the usefulness of MD for understanding the physics for phonon transport quantum sensors. To do this, MD was employed to determine the dynamical matrix by treating atoms as coupled oscillators. An accurate understanding of the dynamical matrix of a system is required to calculate the non-equilibrium Green’s function used to describe the phonon transport within phonon wave-guides. I found that, across multiple pressures and temperatures, randomly placed carbon atoms will show evidence of pent-first formation with semi-empirical models. Density functional theory (DFT), on the other hand, was too computationally expensive to use for full scale MD simulations, but we have the possibility of training a machine learned interatomic potential to approximate DFT for carbon in the environments being studied for pent-first graphene sheet formation.

36 MATERIALS SCIENCE↗

Control Co-Design Framework for Joint Optimization of PID Governor and Hydraulic Dynamics for a Hydropower Plant

This paper presents a control co-design (CCD) framework tailored to the hydraulic and mechanical subsystems of hydropower systems, aimed at enhancing their dynamic performance and responsiveness to grid demands. Conventional sequential design approaches often fall short in capturing the coupled interactions between plant dynamics and control objectives. The proposed CCD methodology enables simultaneous optimization of key hydro plant parameters - such as governor settings, penstock characteristics, and valve control - with the control system architecture. Results indicate that both CCD and sequential design approaches achieve comparable transient and steady-state performance, with minor discrepancies arising from differences in their optimized parameters. However, CCD demonstrates a distinct advantage by tripling the water time constant (Tw) while maintaining similar dynamic performance to the sequential approach. This increase in Tw enables optimization of penstock dimensions, as it is directly influenced by penstock design, thereby contributing to overall cost optimization.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Multiscale electro-chemo-mechanical model of high-capacity amorphous silicon anode material in sodium-ion batteries

The growing interest in sodium-ion batteries (SIBs) is fueled by their abundant resources and environmentally friendly nature, with amorphous silicon (a-Si) emerging as a promising anode material for enhancing capacity. However, the key challenge lies in sustaining reversible capacity during cycling. Here, in this work, we developed a multiscale electrochemical model incorporating an a-Si anode to elucidate the performance parameters of SIBs. Additionally, we integrated an electro-chemo-mechanical model at the particle level to capture stress generation, an essential factor in the degradation of high-capacity electrodes. Unlike existing models, our approach accounts for large-deformation chemo-mechanics at the particle scale and includes simulations under varying charge rates to explore multiscale behavior. The results reveal that coupled sodiation significantly prolongs complete cycling times and influences discharge dynamics, indicating that neglecting this coupling leads to an underestimation of actual capacity. Furthermore, we observed pronounced polarization effects at higher charge rates, resulting in heterogeneous stress distributions across the anode. With the identification of critical failure parameters for both active particles and binder materials, offering novel insights for mitigating degradation in high-capacity electrode systems.

Amorphous silicon anode↗

Experimental determination of energy-level structures of diamond silicon-vacancy centers in off-axis magnetic fields

Here, we report the development of an experimental approach to efficiently determine the energy level structure of an individual silicon vacancy (SiV) center in a magnetic field along an arbitrary direction. This approach uses two coupling rates (one each for the ground and the excited states) to characterize the combined effects of static strain and dynamic Jahn-Teller coupling and exploits the fact that orbital Zeeman effects vanish when the magnetic field is normal to the SiV axis. With an analytical expression for the energy level structure of the SiV under a transverse magnetic field, the two coupling rates can be directly derived from two measurements: one on the frequency separation between two spin-conserved transitions and the other on the coherent population trapping resonance of the SiV ground spin states. A detailed comparison between the numerical calculation and the experimental result on the dependence of the spin-conserved splitting on both the amplitude and direction of the magnetic field further reveals unequal orbital magnetic coupling for the ground and excited states, indicating that unequal orbital quenching factors are needed for an accurate description of the SiV energy level structure in a magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Data-driven equation-free dynamics applied to many-protein complexes: The microtubule tip relaxation

Microtubules (MTs) constitute the largest components of the eukaryotic cytoskeleton and play crucial roles in various cellular processes, including mitosis and intracellular transport. The property allowing MTs to cater to such diverse roles is attributed to dynamic instability, which is coupled to the hydrolysis of GTP (guanosine-5'-triphosphate) to GDP (guanosine-5'-diphosphate) within the β-tubulin monomers. Understanding the equilibrium dynamics and the structural features of both GDP- and GTP-complexed MT tips, especially at an all-atom level, remains challenging for both experimental and computational methods because of their dynamic nature and the prohibitive computational demands of simulating large, many-protein systems. This study employs the “equation-free” multiscale computational method to accelerate the relaxation of all-atom simulations of MT tips toward their putative equilibrium conformation. Using large MT lattice systems (14 protofilaments × 8 heterodimers) comprising ~21-38 million atoms, we applied this multiscale approach to leapfrog through time and nearly double the computational efficiency in realizing relaxed all-atom conformations of GDP- and GTP-complexed MT tips. Commencing from an initial 4 μs unbiased all-atom simulation, we interleave coarse projective “equation-free” jumps with short bursts of all-atom molecular dynamics simulation to realize an additional effective simulation time of 1.875 μs. Our 5.875 μs of effective simulation trajectories for each system expose the subtle yet essential differences in the structures of MT tips as a function of whether β-tubulin monomer is complexed with GDP or GTP, as well as the lateral interactions within the MT tip, offering a refined understanding of features underlying MT dynamic instability. Furthermore, the approach presents a robust and generalizable framework for future explorations of large biomolecular systems at atomic resolution.

Wu, Jiangbo [University of Chicago, IL (United Sta↗

Coupled Induction Machine and HVAC Models for Simulating HVAC Performance Considering Grid Dynamics in Buildings

This paper presents the development of novel models that integrate induction machines with HVAC equipment, such as pumps, heat pumps, and chillers, to analyze the impact of electrical parameters on the operational performance of thermo-fluid systems. The proposed model employs a coupling technique that captures the dynamic interactions between induction machines and HVAC systems. By integrating electrical, thermal, and mechanical dynamics, the models provide a comprehensive framework for simulating real-world scenarios, including interactions with the electrical grid. This achievement was made possible through the development of a Computationally Efficient and Accurate Induction Machine (CEAIM) model. Implemented using the equation-based Modelica language, the CEAIM model has been validated against experimental results, manufacturer data sheets, and various operating conditions. Its performance has been compared with existing induction machine models in the Modelica Standard Library (MSL), demonstrating superior accuracy and computational efficiency. The CEAIM model predicts torque, speed, and power consumption with a coefficient of determination (R 2 ) ranging from 0.98 to 1 and a coefficient of variation of root mean square error (CVRMSE) between 0.27% and 6.67%. Additionally, CEAIM scales more efficiently than conventional MSL models, with a slower computational growth rate in large-scale simulations. After thorough validation of the CEAIM model, it was coupled with HVAC equipment as this approach provides a detailed multi-dimensional view of capturing electrical transients and mechanical performance. To support this, a case study was conducted to showcase its capabilities.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Efficient Quantum Gibbs Samplers with Kubo–Martin–Schwinger Detailed Balance Condition

Lindblad dynamics and other open-system dynamics provide a promising path towards efficient Gibbs sampling on quantum computers. In these proposals, the Lindbladian is obtained via an algorithmic construction akin to designing an artificial thermostat in classical Monte Carlo or molecular dynamics methods, rather than being treated as an approximation to weakly coupled system-bath unitary dynamics. Recently, Chen, Kastoryano, and Gilyén (arXiv:2311.09207) introduced the first efficiently implementable Lindbladian satisfying the Kubo–Martin–Schwinger (KMS) detailed balance condition, which ensures that the Gibbs state is a fixed point of the dynamics and is applicable to non-commuting Hamiltonians. This Gibbs sampler uses a continuously parameterized set of jump operators, and the energy resolution required for implementing each jump operator depends only logarithmically on the precision and the mixing time. In this work, we build upon the structural characterization of KMS detailed balanced Lindbladians by Fagnola and Umanità, and develop a family of efficient quantum Gibbs samplers using a finite set of jump operators (the number can be as few as one), akin to the classical Markov chain-based sampling algorithm. Compared to the existing works, our quantum Gibbs samplers have a comparable quantum simulation cost but with greater design flexibility and a much simpler implementation and error analysis. Moreover, it encompasses the construction of Chen, Kastoryano, and Gilyén as a special instance.

97 MATHEMATICS AND COMPUTING↗