A Framework to Demonstrate a DNP3 Interface with a CIM-Based Data Integration Platform
This poster was presented at the 2024 IEEE Power & Energy Society General Meeting, July 21-25, 2024, Seattle, Washington.
Distributed Network Protocol (DNP3)↗
SEARCH · Engineering Papers
Results for “complex sample environment”
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.
High Semi-Volatile Organic Aerosol Contributions Associated with Ammonium Nitrate in the Coastal Urban Environment
Organic components often contribute 50% or more of the submicron aerosol mass in coastal urban environments, but their partitioning between the gas and particle phases is controlled by a complex mixture of unidentified organic compounds that are poorly constrained by observations. This study compares daily filter organic functional groups (OFG) with online organic mass fragments from La Jolla, California, as part of the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE), quantifying the contributions of four types of non-volatile (NV) organic emission sources to the submicron composition. Daily filters retained only 0.79–0.98 µg/m3 NV submicron organic mass concentration, even though 1.8–1.9 µg/m3 non-refractory (NR) submicron organic mass concentration was measured online. The 62–64% of measured NR submicron organic mass concentration that exceeded what remained on the filters after 23-hr sampling is interpreted as semi-volatile, consistent with the moderate correlation of the NR NV difference to NR ammonium, NR nitrate, and biomass burning-related NV and NR organic factors. The association between semi-volatile organic components and ammonium nitrate likely results from both co-emission and co-evaporation. Size-resolved filter analysis showed that NV organic mass concentration accounts for 68% of NR organic mass concentration for 0.5–1 µm dry diameter but account for 9.0% for 0.18–0.5 µm dry diameter showing the differences in volatility between particle modes. Importantly, volatility of organic components was size-dependent Information Classification: General and associated with ammonium nitrate and biomass burning, providing guidance for constraining atmospheric aerosol properties in global models.
Pelayo, Christian↗
Predicting oxidation damage of ultra high-temperature carbide ceramics in extreme environments using machine learning
Determining the oxidation resistance of UHTC carbides in extreme environments is challenging theoretically and experimentally due to the high dimensional complexity of influencing variables and intricate testing setups. Herein we demonstrate the use of machine learning (ML) models trained with experimental literature data to predict the oxide thickness of UHTC carbides exposed to air based on composition, mean grain size, relative densification, holding time, and temperature. A multi-dimensional database with 76 occurrences is created containing experimental results of Hf, Zr, and Ta carbides plus additives. In this study, the preprocessed database is then used to train ML models to predict their oxidation behavior. The trained model predicts the oxidation damage in the form of an average oxide thickness in UHTC carbides with a Mean Absolute Error (MAE) of ±65.45 μm for samples in the testing set that developed thicknesses up to 1000 μm. The model successfully predicted oxidation damage for a recession rate lower than 60 μm/min. It is noticed that the ensemble method MAE is increased to ±134.34 μm while forecasting the oxidation of samples with a recession rate higher than the threshold. The unprecedented approach is a novel way to predict the damage through the oxidation of carbide compounds before processing for a smarter design with room for improvement.
36 MATERIALS SCIENCE↗
Enabling machine learning-ready HPC ensembles with Merlin
With the growing complexity of computational and experimental facilities, many scientific researchers are turning to machine learning (ML) techniques to analyze large scale ensemble data. With complexities such as multi-component workflows, heterogeneous machine architectures, parallel file systems, and batch scheduling, care must be taken to facilitate this analysis in a high performance computing (HPC) environment. Here, we present Merlin, a workflow framework to enable large ML-friendly ensembles of scientific HPC simulations. By augmenting traditional HPC with distributed compute technologies, Merlin aims to lower the barrier for scientific subject matter experts to incorporate ML into their analysis. As a producer–consumer workflow model, Merlin enables multi-machine, cross-batch job, dynamically allocated yet persistent workflows capable of utilizing surge-compute resources. Key features of Merlin are a flexible HPC-centric interface, low per-task overhead, multi-tiered fault recovery, and a hierarchical sampling algorithm that allows for $\mathscr{O}$(N) task execution and $\mathscr{O}$(N ln N) task queuing to ensembles of millions of tasks. In addition to Merlin’s design, we test the algorithm’s performance in an HPC center and demonstrate the ability to enqueue 40 million simulations in 100 s, with a 30 millisecond per-task overhead that is independent of ensemble size. Finally, we describe some example applications that Merlin has enabled on leadership-class HPC resources, such as the ML-augmented optimization of nuclear fusion experiments and the calibration of infectious disease models to study the progression of and possible mitigation strategies for COVID-19.
97 MATHEMATICS AND COMPUTING↗
Studying Aerosol, Clouds, and Air Quality in the Coastal Urban Environment of Southeastern Texas
A multi-agency succession of field campaigns was conducted in southeastern Texas during July 2021 through October 2022 to study the complex interactions of aerosols, clouds and air pollution in the coastal urban environment. As part of the Tracking Aerosol Convection interactions Experiment (TRACER), the TRACER- Air Quality (TAQ) campaign the Experiment of Sea Breeze Convection, Aerosols, Precipitation and Environment (ESCAPE) and the Convective Cloud Urban Boundary Layer Experiment (CUBE), a combination of ground-based supersites and mobile laboratories, shipborne measurements and aircraft-based instrumentation were deployed. These diverse platforms collected high-resolution data to characterize the aerosol microphysics and chemistry, cloud and precipitation micro- and macro-physical properties, environmental thermodynamics and air quality-relevant constituents that are being used in follow-on analysis and modeling activities. We present the overall deployment setups, a summary of the campaign conditions and a sampling of early research results related to: (a) aerosol precursors in the urban environment, (b) influences of local meteorology on air pollution, (c) detailed observations of the sea breeze circulation, (d) retrieved supersaturation in convective updrafts, (e) characterizing the convective updraft lifecycle, (f) variability in lightning characteristics of convective storms and (g) urban influences on surface energy fluxes. The work concludes with discussion of future research activities highlighted by the TRACER model-intercomparison project to explore the representation of aerosol-convective interactions in high-resolution simulations.
54 ENVIRONMENTAL SCIENCES↗
A Benchmark to Test Generalization Capabilities of Deep Learning Methods to Classify Severe Convective Storms in a Changing Climate
Abstract This is a test case study assessing the ability of deep learning methods to generalize to a future climate (end of 21st century) when trained to classify thunderstorms in model output representative of the present‐day climate. A convolutional neural network (CNN) was trained to classify strongly rotating thunderstorms from a current climate created using the Weather Research and Forecasting model at high‐resolution, then evaluated against thunderstorms from a future climate and found to perform with skill and comparatively in both climates. Despite training with labels derived from a threshold value of a severe thunderstorm diagnostic (updraft helicity), which was not used as an input attribute, the CNN learned physical characteristics of organized convection and environments that are not captured by the diagnostic heuristic. Physical features were not prescribed but rather learned from the data, such as the importance of dry air at mid‐levels for intense thunderstorm development when low‐level moisture is present (i.e., convective available potential energy). Explanation techniques also revealed that thunderstorms classified as strongly rotating are associated with learned rotation signatures. Results show that the creation of synthetic data with ground truth is a viable alternative to human‐labeled data and that a CNN is able to generalize a target using learned features that would be difficult to encode due to spatial complexity. Most importantly, results from this study show that deep learning is capable of generalizing to future climate extremes and can exhibit out‐of‐sample robustness with hyperparameter tuning in certain applications.
54 ENVIRONMENTAL SCIENCES↗
Module-Fluidics: Building Blocks for Spatio-Temporal Microenvironment Control
Generating the desired solute concentration signal in micro-environments is vital to many applications ranging from micromixing to analyzing cellular response to a dynamic microenvironment. We propose a new modular design to generate targeted temporally varying concentration signals in microfluidic systems while minimizing perturbations to the flow field. The modularized design, here referred to as module-fluidics, similar in principle to interlocking toy bricks, is constructed from a combination of two building blocks and allows one to achieve versatility and flexibility in dynamically controlling input concentration. The building blocks are an oscillator and an integrator, and their combination enables the creation of controlled and complex concentration signals, with different user-defined time-scales. We show two basic connection patterns, in-series and in-parallel, to test the generation, integration, sampling and superposition of temporally-varying signals. All such signals can be fully characterized by analytic functions, in analogy with electric circuits, and allow one to perform design and optimization before fabrication. Such modularization offers a versatile and promising platform that allows one to create highly customizable time-dependent concentration inputs which can be targeted to the specific application of interest.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
High-temperature active oxidation of nanocrystalline silicon-carbide: A reactive force-field molecular dynamics study
Flexible woven SiC ceramics are prone to accelerated fiber embrittlement under high temperature oxidation in dynamic oxygen environments. The nanocrystalline structure of the constituent fibers impacts the reaction kinetics and phase transformations during active oxidation. However, fundamental understanding and quantification of grain boundary effects on oxidation behavior in nanocrystalline SiC remain elusive when temperatures exceed 1500 K. This study deploys large-scale molecular dynamics simulations with a reactive force-field to elucidate the complex roles of atomic oxygen reservoir conditions and grain size on oxidation kinetics and the nature of oxides produced in both monocrystalline and nanocrystalline 3C-SiC between 1100 K and 2000 K. The simulations with dynamically replenished oxygen provide good agreement with oxidation kinetics and activation energies for the monocrystalline Si(100) and C(100) orientations published in the available literature. This study reveals that, by contrast, nanocrystalline SiC samples exhibit two distinct oxidation kinetics with a transition point at 1500 K due to surface melting, which is supported by experimental evidence. The introduction of a grain-boundary network produces a two-fold decrease in oxidation activation energies compared to monocrystalline SiC below 1500 K. Above 1500 K, however, the activation energies rise substantially due to the formation of a liquid Si phase at the SiC/Si oxide interface. Further, it is shown that the stability of the interfacial liquid phase is promoted by incoherent grain boundaries in the crystalline SiC. These findings are important for the deployment of nanocrystalline SiC fibers in advanced thermal protection systems for high-temperature applications.
36 MATERIALS SCIENCE↗
Microwave-Based Synthesis of Functional Morphological Variants and Carbon Nanotube-Based Composites of VS 4 for Electrochemical Applications
A novel facile, fast, and efficient microwave-assisted method was developed to synthesize a number of diverse nanostructured motifs (ranging from nanorods to nanoflowers) of VS 4 along with its associated composite heterostructures, VS 4 /multi-walled carbon nanotube (MWNT; i.e., multi-walled carbon nanotubes). Specifically, we have probed and correlated the effects of a number of specific experimental variables, including primarily precursor, solvent, temperature, and time. We noted that nanorods formed more readily with VO(acac) 2 as the vanadium precursor and n-methyl-2-pyrrolidone (NMP) as a polar reaction solvent. By contrast, we determined that hierarchical three-dimensional (3D) nanoflower-like assemblies, ranging in size from 100 to 200 nm in average diameter, could be controllably synthesized by using Na 3 VO 4 as the vanadium precursor and an aqueous water: polar solvent mixture as the reaction medium. We also observed that VS 4 disintegrates, when in the presence of either air, solution, or a combination of these environments, and established that the extent of VS4 nanorod decomposition could be almost fully prevented by storage under nitrogen. From an application’s perspective, our VS 4 is electrochemically active and shows behavior, consistent with the literature. In particular, as compared with pristine VS 4 nanorods alone, we observed enhanced electrochemical activity with (i) 3D hierarchical flower-like motifs, (ii) unique necklace-like VS 4 nanorod–MWNT composites, and (iii) samples in which as-prepared VS4 nanorods had been annealed. Moreover, we found that the rational application of specific physical and chemical processing treatments, such as (i) thermal annealing to improve crystallinity, (ii) the addition of MWNTs to form conductive composites, and (iii) the evolution of morphology from one-dimensional (1D) nanorods to more complex 3D nanoflowers, was favorable to the resulting electrochemical performance with respect to increasing stability and reversibility.
25 ENERGY STORAGE↗
Actinide Molten Salts: A Machine-Learning Potential Molecular Dynamics Study
We know that actinide molten salts represent a class of important materials in nuclear energy. Understanding them at a molecular level is critical to proper and optimal design of relevant technological applications. Yet, owing to the complexity of electronic structure due to the 5f orbitals, computational studies of heavy elements in condensed phases using ab initio potentials to study the structure and dynamics of these elements embedded in molten salts are difficult. This lack of efficient computational protocols makes it difficult to obtain information on properties that require extensive statistical sampling like transport. To tackle this problem, we adopted a machine-learning approach to study ThCl 4 -NaCl and UCl 3 -NaCl binary systems. The machine-learning potential, with the density functional theory accuracy, allows us to obtain long molecular dynamics trajectories (ns) for large systems (10 3 atoms) at a considerably low computing cost, thereby efficiently gaining information about their bonding structures, thermodynamics, and dynamics at a range of temperature. We observed a considerable change in the coordination environments of actinide elements and their characteristic coordination-sphere lifetime. Our study also suggests that actinides in molten salts may not follow well known entropy-scaling laws.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Fiberboard Properties and Degradation Rates for Storage of 9975 Shipping Package in KAC (2020 Status Report)
Savannah River Site (SRS) stores plutonium materials within model 9975 shipping packages in the K-Area Complex (KAC). The 9975 shipping package consists of a 35 gallon stainless steel drum, Celotex fiberboard insulation, lead shield, and primary and secondary containment vessels. The 9975 shipping package design, performance, and analysis for safe transport of radioactive material are described in the Safety Analysis Report for Packaging (SARP). Celotex fiberboard provides three safety functions: thermal insulation to limit internal temperature during a fire, criticality control, and resistance to package crushing. The fiberboard material must retain its dimensions and density within certain ranges to provide the required impact resistance, criticality control, and fire resistance. The SRS Surveillance Program monitors material performance to establish a basis for service life and ensures the continued integrity of 9975 packages. Fiberboard samples, taken from multiple fiberboard assemblies fabricated from cane and softwood fiberboard, have been conditioning in elevated temperature environments since 2005. The samples are periodically examined to monitor thermal, mechanical, and physical properties, and assess degradation trends. Fiberboard properties of interest that are evaluated to demonstrate acceptable long-term performance include dimensional stability, density, compressive strength, thermal conductivity, and specific heat capacity. Duplicate samples from multiple package sources have been tested to identify the range of variability in fiberboard properties and degradation rates. Baseline and long-term testing of fiberboard material properties have been reported previously; reference 6 summarized experimental results of cane and softwood fiberboard through May 2019 and presented degradation models for the measured properties. This report presents the cumulative data collected through June 2020 and the corresponding updated aging models.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Macroscopic and microscopic behaviors of Mn(II) (ad)sorption to goethite with the effects of dissolved carbonates under anoxic conditions
Manganese and its compounds have been extensively researched because of their far-reaching roles in a wide range of biogeochemical processes in natural systems. The (ad)sorption behavior of Mn(II), however, is poorly understood despite its important role as the primary reaction step for surface-catalyzed Mn(II) oxidation that is the principal abiotic process forming various Mn (oxyhydr)oxides in nature. Here, we systematically examined Mn(II) (ad)sorption to one of the most common natural sorbents, goethite, in oxygen- and carbonate-free systems. Traditional sorption edge and isotherm experiments were conducted by varying sorbate-to-sorbent ratio (0.027–15 μmol∙m -2 ) and solution pH (pH 5.0–9.0). The effects of dissolved carbonates on Mn(II) sorption were also assessed in a range of naturally prevalent concentrations (0.5–10 mM NaHCO 3 ). The Mn(II) uptake on goethite followed a typical Langmuir isotherm in the absence of dissolved carbonates, with increasing maximum adsorption capacities (Γ max ) from 0.19 at pH 6.5 to 3.4 μmol∙m -2 at pH 9.0. The presence of dissolved carbonates raised the extent of Mn(II) adsorption, which appeared to be directly correspondent to that of the adsorption of dissolved carbonates. Extended X-ray absorption fine structure (EXAFS) analysis indicated that Mn(II) predominantly formed inner-sphere binuclear bidentate surface complexes. Mn(II) uptake became deviated from the Langmuir model and showed a clear indication of surface precipitation when the Mn(II) sorption density (Γ) exceeded a threshold value in a given solution composition. This secondary Mn(II) phase was identified as rhodochrosite using X-ray diffraction (XRD) and transmission electron microscope (TEM). Furthermore, it would be plausible that a minor fraction of adsorbed Mn(II) coexisted with the secondary rhodochrosite according to a linear combination fitting (LCF) of the X-ray absorption near edge structure (XANES) spectra of Mn reference and samples. These systematic investigations of the macroscopic and microscopic behaviors of Mn(II) (ad)sorption to goethite provide a critical avenue for disentangling surface-catalytic Mn(II) oxidation processes, which ultimately lead to the formation of diverse Mn (oxyhydr)oxides in the environment.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Operational Analytics Studies for ATLAS Distributed Computing: Data Popularity Forecast and Utilization of the WLCG Centers
Operational analytics is the direction of research related to the analysis of the current state of computing processes and the prediction of future states in order to anticipate imbalances and take timely measures to stabilize a complex system. There are two relevant areas in ATLAS Distributed Computing that are currently the focus of studies: user physics analysis including the forecast of popularity of data samples among users, and evaluating WLCG centers for their readiness to process user analysis payloads. Studying these areas is challenging due to the complexity involved, as it requires a comprehensive understanding of numerous boundary conditions typically found in large-scale distributed computing infrastructures. Forecasts of data popularity are problematic without the categorization of user tasks by their types (data transformation or physics analysis), which do not always appear on the surface but may induce noise, which introduces significant distortions for predictive analysis. Evaluating the WLCG resources by their analysis workloads is also a challenging task as it is necessary to find a balance between the workload of the resource, its performance, the waiting time for jobs on it, as well as the volume of jobs that it processes. This is especially difficult in a heterogeneous computing environment, where legacy resources are used along with modern high-performance machines. We will look at these areas of research in detail and discuss what tools and methods are used in our work, demonstrating results already obtained.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
High Energy Density Physics of Inertial Confinement Fusion Ablator Materials
The historic December 5, 2022 experiment at Lawrence Livermore National Lab’s (LLNL) National Ignition Facility (NIF) reached fusion energy ignition for the first time. This is the most important scientific breakthrough of the 21st century paves the way to future clean inertial fusion energy (IFE). The diamond (high density carbon (HDC)) ablator material used in this experiment displays detrimental effects due to the development of hydrodynamic instabilities at the diamond/fuel interface under shock compression. New alternatives to diamond ablators are required to step up the energy yield in ICF experiments. The unique combination of mechanical strength (approaching that of diamond), the ability to accommodate high-Z dopants (in contrast to diamond), and the tunability of the properties (through synthesis material with varying sp 3 content) make amorphous carbon (a-C) a promising material for next-generation IFE ablative capsules. However, despite its critical importance to the IFE program, the behavior of a-C carbon at extreme temperatures and pressures remains largely unexplored. The primary goals of this project were to perform groundbreaking dynamic compression experiments and predictive simulations to uncover the fundamental high-energy-density physics of amorphous carbon. Our goals were (1) to uncover the metastability range of amorphous carbon and probe phase transitions to diamond or metastable supercooled liquid carbon; (2) to acquire high-quality equation of state (EOS) data and develop an experimentally validated EOS from machine-learning MD simulations of the complex states of carbon; and (3) to uncover the complex behavior of carbon liquid in both thermodynamically stable and metastable supercooled states by accessing large areas of carbon phase diagram with amorphous samples with variable sp 3 content. Our proposed experimental program included measurements of equation of state and diffraction measurements using the Omega EP laser at the Laboratory of Laser Energetics at the University of Rochester. The theoretical/simulation program involved the development of machine-learning models of the complex response of amorphous carbon under dynamic compression by performing molecular dynamics simulations at experimental time and length scales using leadership class DOE supercomputers. Simulations guided experiments to observe predicted phenomena and acquire critical experimental data in specific pressure-temperature domains to validate theoretical models. This research delivered fundamental properties of novel amorphous carbon IFE ablator material, including phase diagram and EOS. These results will aid in IFE target design and implosion experiments. A unique combination of predictive simulations and dynamic and static experiments provided a highly inspirational intellectual environment for graduate students and postdocs involved in this project.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Characterizing Radiation-Induced Chromium Redox Chemistry
The redox behavior of Cr(III) and Cr(VI) in nuclear environments were studied by analyzing aqueous irradiated chromium ions using UV-Visible spectroscopy. The formation of different chromium oxidation states can occur through gamma radiolysis of chromium solutions. Characterization of radiation-induced chromium is not well understood, and yet important as chromium can make its way into the primary coolant of reactors because of the corrosion of stainless steel reactor components. Furthermore, Cr(VI) is toxic and not suitable for environmental release, therefore Cr speciation is important. Fricke dosimetry was performed to collect doses rates, specific extinction coefficients were obtained from various analytical methods for the two oxidation states of chromium. Chromium samples--conditioned at various pH's and concentration to understand the effects--were then irradiated using a Coblat-60 gamma irradiator. Three optical spectroscopy analytical methods were used , including direct oxidation state measurements and complexation by either EDTA or DPC. Complexation by EDTA was found to not be a reliable method. Radiolytic oxidation of Cr(III) to Cr(VI), and reduction of Cr(VI) to Cr(III) were clearly seen, except for Cr(III) samples were a pH less than or equal to 2, where no change occurred.
38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
Machine learning bridges local static structure with multiple properties in metallic glasses
A long-standing challenge in the metallic glass (MG) community has been how to quantitatively gauge the influence of the intricate local packing environment on the response (such as the propensity for atomic rearrangement) of the atomic configuration to external stimuli. Here we establish this structure–property relation by representing the complex amorphous structure using a single, flexibility-orientated structural quantity. This structural flexibility (SF) couples to a bona fide structural representation, the pair distribution function (PDF) of individual atoms, through a weighting function that reflects what matters in the static atomic configuration to dynamic responses. Machine learning is used, employing microscopic flexibility volume as the supervisory signal, to establish via direct regression an optimized weighting vector, which is proven robust for all quenching rates, deformation conditions, and different compositions in a given (e.g., Cu x Zr 100-x ) alloy system. Additionally, the SF is evaluated solely from the particle positions (PDF), for any structure variation, from the atomic scale up to sample average. Strong correlations are demonstrated between SF and a broad range of properties, including vibrational, diffusional, as well as elastic and plastic relaxation responses.
36 MATERIALS SCIENCE↗
Landau–Zener Transition Enhanced Quantum Sensing in Spin Defects of Hexagonal Boron Nitride
Negatively charged boron vacancies (V $^{–}_{B}$ ) in hexagonal boron nitride (hBN) comprise a promising quantum sensing platform, optically addressable at room temperature and transferable onto samples. However, broad hyperfine-split spin transitions of the ensemble pose challenges for quantum sensing with conventional resonant excitation due to limited spectral coverage. While V $^{–}_{B}$ in isotopically enriched hBN using 10 B and 15 N isotopes (h 10 B 15 N) exhibits sharper spectral features, significant inhomogeneous broadening persists. We show that, implemented via frequency modulation on an FPGA, a frequency-ramped microwave pulse achieves around 4-fold greater |0⟩→|−1⟩ spin-state population transfer and thus contrast than resonant microwave excitation and thus 16-fold shorter measurement time for spin relaxation-based quantum sensing. Quantum dynamics simulations reveal that an effective two-state Landau–Zener model captures the complex relationship between population inversion and pulse length with relaxations incorporated. Our approach is robust and valuable for quantum relaxometry with spin defects in hBN, especially in noisy environments.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Quantum Markov chain Monte Carlo with digital dissipative dynamics on quantum computers
Modeling the dynamics of a quantum system connected to the environment is critical for advancing our understanding of complex quantum processes, as most quantum processes in nature are affected by an environment. Modeling a macroscopic environment on a quantum simulator may be achieved by coupling independent ancilla qubits that facilitate energy exchange in an appropriate manner with the system and mimic an environment. This approach requires a large, and possibly exponential number of ancillary degrees of freedom which is impractical. In contrast, we develop a digital quantum algorithm that simulates interaction with an environment using a small number of ancilla qubits. By combining periodic modulation of the ancilla energies, or spectral combing, with periodic reset operations, we are able to mimic interaction with a large environment and generate thermal states of interacting many-body systems. We evaluate the algorithm by simulating preparation of thermal states of the transverse Ising model. Our algorithm can also be viewed as a quantum Markov chain Monte Carlo process that allows sampling of the Gibbs distribution of a multivariate model. To illustrate this we evaluate the accuracy of sampling Gibbs distributions of simple probabilistic graphical models using the algorithm.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗