Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constitutive 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.

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

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Plug-and-Play Methods for Integrating Physical and Learned Models in Computational Imaging: Theory, algorithms, and applications

Plug-and-play (PnP) priors constitute one of the most widely used frameworks for solving computational imaging problems through the integration of physical models and learned models. PnP leverages high-fidelity physical sensor models and powerful machine learning methods for prior modeling of data to provide state-of-the-art reconstruction algorithms. PnP algorithms alternate between minimizing a data fidelity term to promote data consistency and imposing a learned regularizer in the form of an image denoiser. Recent highly successful applications of PnP algorithms include biomicroscopy, computerized tomography (CT), magnetic resonance imaging (MRI), and joint ptychotomography. This article presents a unified and principled review of PnP by tracing its roots, describing its major variations, summarizing main results, and discussing applications in computational imaging. Additionally, we also point the way toward further developments by discussing recent results on equilibrium equations that formulate the problem associated with PnP algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Redshift inference from the combination of galaxy colours and clustering in a hierarchical Bayesian model – Application to realistic N -body simulations

ABSTRACT Photometric galaxy surveys constitute a powerful cosmological probe but rely on the accurate characterization of their redshift distributions using only broad-band imaging, and can be very sensitive to incomplete or biased priors used for redshift calibration. A hierarchical Bayesian model has recently been developed to estimate those from the robust combination of prior information, photometry of single galaxies, and the information contained in the galaxy clustering against a well-characterized tracer population. In this work, we extend the method so that it can be applied to real data, developing some necessary new extensions to it, especially in the treatment of galaxy clustering information, and we test it on realistic simulations. After marginalizing over the mapping between the clustering estimator and the actual density distribution of the sample galaxies, and using prior information from a small patch of the survey, we find the incorporation of clustering information with photo-z’s tightens the redshift posteriors and overcomes biases in the prior that mimic those happening in spectroscopic samples. The method presented here uses all the information at hand to reduce prior biases and incompleteness. Even in cases where we artificially bias the spectroscopic sample to induce a shift in mean redshift of $\Delta \bar{z} \approx 0.05,$ the final biases in the posterior are $\Delta \bar{z} \lesssim 0.003.$ This robustness to flaws in the redshift prior or training samples would constitute a milestone for the control of redshift systematic uncertainties in future weak lensing analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Variational multiscale reinforcement learning for discovering reduced order closure models of nonlinear spatiotemporal transport systems

Abstract A central challenge in the computational modeling and simulation of a multitude of science applications is to achieve robust and accurate closures for their coarse-grained representations due to underlying highly nonlinear multiscale interactions. These closure models are common in many nonlinear spatiotemporal systems to account for losses due to reduced order representations, including many transport phenomena in fluids. Previous data-driven closure modeling efforts have mostly focused on supervised learning approaches using high fidelity simulation data. On the other hand, reinforcement learning (RL) is a powerful yet relatively uncharted method in spatiotemporally extended systems. In this study, we put forth a modular dynamic closure modeling and discovery framework to stabilize the Galerkin projection based reduced order models that may arise in many nonlinear spatiotemporal dynamical systems with quadratic nonlinearity. However, a key element in creating a robust RL agent is to introduce a feasible reward function, which can be constituted of any difference metrics between the RL model and high fidelity simulation data. First, we introduce a multi-modal RL to discover mode-dependant closure policies that utilize the high fidelity data in rewarding our RL agent. We then formulate a variational multiscale RL (VMRL) approach to discover closure models without requiring access to the high fidelity data in designing the reward function. Specifically, our chief innovation is to leverage variational multiscale formalism to quantify the difference between modal interactions in Galerkin systems. Our results in simulating the viscous Burgers equation indicate that the proposed VMRL method leads to robust and accurate closure parameterizations, and it may potentially be used to discover scale-aware closure models for complex dynamical systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Rates of Sea‐Level Rise Are Highly Sensitive to Ice Viscosity Parameters in Model Benchmarks

Glacier flow plays a major role in current and future rates of globally averaged sea-level rise. The viscosity of glacial ice, controlling the rate of flow, decreases as stress increases and is highly sensitive to the value of the stress exponent, $n$, in the constitutive equation for viscous flow. Glaciologists and climate modelers almost exclusively assume $n=3$ when modeling ice flow and projecting sea-level rise through forward modeling. However, recent work suggests that $n\approx 4$ better fits observations, prompting the question: How sensitive are projections of sea-level rise to the value of $n$? We use an established community ice flow model and standard benchmark experiments designed as an idealized representation of Pine Island Glacier, West Antarctica. While initializing an $n=3$ model to match observations of an $n=4$ ice sheet is possible, we find that incorrectly assuming $n=3$ when in fact $n=4$ dramatically underestimates rates of sea-level rise. The scale of this error grows nonlinearly with the magnitude of the climate forcing, acting to increase projection uncertainties. Additionally, we find that models often account for this stress-dependent rheology mismatch during model initialization in a way that masks this rheological effect in the short term while leaving model outputs vulnerable to larger biases in longer-term projections. Initializations to observations of Pine Island Glacier display similar rheology-mismatch fingerprints to our idealized example.

climate sensitivity↗
DOE OSTI2022

SAM-ML: Integrating data-driven closure with nuclear system code SAM for improved modeling capability

Advanced reactors often involve complicated thermal-fluid (T-F) phenomena. Modeling such phenomena with the traditional one-dimensional (1-D) system code is a challenging task. The System Analysis Module (SAM), a modern nuclear system code, has developed a coarse mesh multi-dimensional (multi-D) flow model to capture the spatial effect of T-F phenomena in advanced reactors. As a coarse mesh solver, constitutive relations are required for SAM's multi-D model for unresolved fine-scale physics, such as turbulence. Here this work presents a novel approach that integrates neural networks as data-driven closure for SAM's multi-D flow model. The data-driven closure is trained with fine-resolution data to ensure its accuracy while maintaining a coarse mesh setup to ensure its efficiency and consistency with SAM. We demonstrate the applicability of this SAM-ML capability in an open volume thermal stratification problem, where a neural network model serves as the eddy viscosity closure. A customized interface between the neural network and SAM is developed to ensure flexible and efficient data exchange. The SAM-ML results demonstrate superior performance compared to SAM's built-in zero-equation eddy viscosity closure. The case study shows that although the generalization capability of the data-driven closure still needs to be improved for different transient case or different geometric setup, SAM -ML demonstrates good potential for challenging simulation problems with improved accuracy and computational efficiency.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

An Improved Fracture Mechanics-Informed Multiscale Thermomechanical Damage Model for Ceramic Matrix Composites

This paper extends recent work done by the authors in modeling length scale-dependent damage behavior of ceramic matrix composites (CMCs) to include effects of local anisotropy introduced by matrix cracking. This model captures scale-dependent damage initiation and propagation behavior of the brittle matrix by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress/strain constitutive relationships at each length scale. The damage ISV captures the effects of matrix cracking and growth by using fracture mechanics and the self-consistent scheme to determine the reduced stiffness of the cracked matrix. Matrix cracks, which activate when stress intensity factors near manufacturing induced cavities exceed the fracture toughness of the material, are assumed to be transversely isotropic in the plane of the crack, and matrix anisotropy occurs when the damaged stiffness tensor is rotated from the crack plane to the global axes. The crack progression and temporal evolution of the damage ISV are governed by fracture mechanics and crack growth kinetics. The model effectively captures first matrix cracking, which is the first significant deviation from linear elasticity. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide (SiC) and a 2D woven five-harness satin (5HS) carbon fiber SiC matrix (C/SiC) CMC.

Skinner, Travis↗
DOE OSTI2023

Establishment of a genome editing tool using CRISPR-Cas9 ribonucleoprotein complexes in the non-model plant pathogen Sphaerulina musiva

CRISPR-Cas9 is a versatile genome editing system widely used since 2013 to introduce site-specific modifications into the genomes of model and non-model species. This technology is used in various applications, from gene knock-outs, knock-ins, and over-expressions to more precise changes, such as the introduction of nucleotides at a targeted locus. CRISPR-Cas9 has been demonstrated to be easy to establish in new species and highly efficient and specific compared to previous gene editing strategies such as Zinc finger nucleases and transcription activator-like effector nucleases. Grand challenges for emerging CRISPR-Cas9 tools in filamentous fungi are developing efficient transformation methods for non-model organisms. In this paper, we have leveraged the establishment of CRISPR-Cas9 genome editing tool that relies on Cas9/sgRNA ribonucleoprotein complexes (RNPs) in the model species Trichoderma reesei and developed the first protocol to efficiently transform the non-model species, Sphaerulina musiva. This fungal pathogen constitutes a real threat to the genus Populus, a foundational bioenergy crop used for biofuel production. Herein, we highlight the general considerations to design sgRNAs and their computational validation. We also describe the use of isolated protoplasts to deliver the CRISPR-Cas9 RNP components in both species and the screening for targeted genome editing events. The development of engineering tools in S. musiva can be used for studying genes involved in diverse processes such as secondary metabolism, establishment, and pathogenicity, among many others, but also for developing genetic mitigation approaches. The approach described here provides guidance for potential development of transformation systems in other non-model spore-bearing ascomycetes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Projecting Antarctica's contribution to future sea level rise from basal ice shelf melt using linear response functions of 16 ice sheet models (LARMIP-2)

Abstract. The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4 mm of global sea level rise, with a standard deviation of 3.7 mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 ∘C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles. We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Machine-learning based model reduction for partial differential equations

We develop a novel synergistic approach between model reduction and machine learning. The specific goal of this project is to aid in the construction of reduced order models for basis functions that are custom-made to represent the solution of partial differential equations. Partial differential equations (PDEs) are one of the main mathematical tools for describing physical phenomena. However, due to either efficiency or necessity, for many real-world problems, we are interested in constructing reduced order models (ROMs) which focus only on the explicit computation of subsets of the active spatio-temporal scales in the problem, while treating the interaction with the rest of the scales approximately. The task of accurate representation of such interactions (usually called memory terms) constitutes a vast area of research known as model reduction. PI Stinis has significant expertise in the construction of ROMs for complex systems. In addition, in recent work with the project key participant Qadeer, they have utilized machine learning to acquire custom-made basis functions (CBFs) to expand the solutions of PDEs. In the proposed work, we will merge the two concepts by constructing ROMs for subsets of the CBFs needed to represent the solution of a PDE. Specifically, we will use the Mori-Zwanzig model reduction formalism to construct ROMs for subsets of CBFs for nonlinear PDEs of various complexity, as well as investigate the usage of CBFs in the spectral vanishing viscosity method for problems that can form shocks in finite time. The outcome of the research is aimed to be proof-of-concept about a novel synergistic approach between model reduction and machine learning, thus advancing the field of scientific machine learning. Such a capability will benefit the efficient modeling of physical systems appearing in various areas of interest to the DOE.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A radioisotope - enabled reactive transport model for deep vadose zone carbon

In mountainous regions, which constitute the principle source of recharge to major rivers and regional aquifers, infiltration occurs through fractured, partially saturated, weathered bedrock that acts as a boundary layer between saturated aquifers and surface soil. Commonly this deep vadose zone (DVZ) is many meters thick, and yet its role in regulating the generation, retention and mobility of reactive solutes, including nutrients, contaminants and weathering products, is largely unknown. In particular, many of the key reactions that drive the formation of the weathered DVZ and the quality of water moving through it are redox processes, regulated by the availability of organic carbon and oxygen below the soil layer. The role of the DVZ is thus also poorly constrained in the context of carbon stocks and mobility, particularly in lithologies that are naturally high in organic carbon, such as shales. The overarching hypothesis of this study is that upland regions developed in geologic settings with abundant petrogenic carbon store and actively cycle carbon in the weathered DVZ below the soil and above the water table at rates that are significant and currently unconstrained. In order to quantify this cycling, the current study combines novel instrumentation techniques allowing new direct sampling of DVZ systems with advanced numerical reactive transport simulations of carbon transport and transformation. Critically, these simulations will explicitly treat the three isotopes of carbon (the abundant 12C, the stable rare 13C and the radioactive 14C) in a unified framework, thus clearly parsing between the contributions of modern surface derived carbon and lithologic carbon sources in integrated measurements of fluid and gas phase fluxes. This novel model capability will be applied to test the role of DVZ carbon cycling as a regulator of water quality and geological weathering in two complementary field sites both located in organic carbon rich shale lithologies. The first is the Eel River Critical Zone Observatory (ERCZO) in Mendocino County, California, and the second is the Lawrence Berkeley National Laboratory Watershed Function Scientific Focus Area (SFA) in the East River watershed, near Crested Butte, Colorado. At the ERCZO, a novel vadose zone monitoring system has been installed in a 20 m thick, partially saturated, weathered shale hillslope, and preliminary data already indicate substantial CO2 flux generated many meters below the soil surface. At the SFA field site, an instrumented hillslope transect indicates a more complex multi-dimensional fluid and solute transport regime, which will serve as a key test of the calibrated models. Collectively, this project will advance understanding of the cycling of carbon belowground and in relation to transport pathways across the poorly constrained DVZ characteristic of primary water recharge areas. The key product of this work will be enhanced isotope simulation capabilities that are robust and publicly available for application across a broad diversity of systems.

58 GEOSCIENCES↗
DOE OSTI2025

Cohesive instability in elastomers: insights from a crosslinked Van der Waals fluid model

Abstract The resistance to volumetric deformations displayed by polymer networks is largely due to secondary and tertiary interactions between neighboring polymer chains. These interactions are both entropic and enthalpic in nature but are fundamentally different from the entropic forces that resist shearing in these networks. In this paper, we introduce a new depiction of elastomers as a crosslinked Van der Waals fluid. Starting from first principles, we develop constitutive equations that are implemented in a continuum model as well as a discrete network model. Our models predict that the failure of polymer networks may be driven by an instability in the underlying polymer bulk ‘fluid’ or by the breaking of polymer chains, depending on the loading path taken. The results of this study indicate that material failure in elastomers exposed to a purely triaxial state, such as in a poker chip experiment, may be driven by an entirely different mode of instability than those deformed in pure shear, such as in a uniaxial tension experiment.

Lamont, Samuel C.↗
DOE OSTI2024

Machine Learning Aided Modeling of Granular Materials: A Review

Artificial intelligence (AI) has become a buzzy word since Google’s AlphaGo beat a world champion in 2017. In the past five years, machine learning as a subset of the broader category of AI has obtained considerable attention in the research community of granular materials. This work offers a detailed review of the recent advances in machine learning-aided studies of granular materials from the particle-particle interaction at the grain level to the macroscopic simulations of granular flow. This work will start with the application of machine learning in the microscopic particle-particle interaction and associated contact models. Then, different neural networks for learning the constitutive behaviour of granular materials will be reviewed and compared. Finally, the macroscopic simulations of practical engineering or boundary value problems based on the combination of neural networks and numerical methods are discussed. We hope readers will have a clear idea of the development of machine learning-aided modelling of granular materials via this comprehensive review work.

42 ENGINEERING↗
DOE OSTI2025

Coupled THM modeling of bentonite heating and hydration in tank tests with a new temperature-dependent water retention model

This study presents a coupled thermo-hydro-mechanical (THM) model for simulating the heating and hydration behavior of bentonite, a buffer material in deep geological repositories (DGRs). The model incorporates a new temperature-dependent soil water retention curve which captures the thermal-induced shift in water retention behavior. It also distinguishes between liquid and gas permeability, modeling intrinsic gas permeability as a function of accessible porosity to improve vapor transport and desaturation predictions. The model was validated against two large-scale tank tests, demonstrating good agreement with measured temperature, relative humidity, and water inflow data. It revealed a complex porosity evolution driven by thermal expansion, vapor movement, vapor condensation, and hydration-induced swelling during heating and hydration processes. The simulation results also suggest that the permeability of the hydration layer plays a critical role in controlling water intake. Clogging of this layer can significantly reduce the volume of water inflow during the hydration phase. Furthermore, while the model effectively captures key THM behavior, further development of the mechanical constitutive law is required to account for possible thermo-elasto-plastic volume changes and microstructural effects. Overall, the model provides a robust tool for evaluating the evolution of bentonite-based barrier material in DGRs.

Guo, Guanlong [Lawrence Berkeley National Laborato↗
DOE OSTI2024

ChatGPT and Other Large Language Models for Cybersecurity of Smart Grid Applications

Cybersecurity breaches targeting electrical substations constitute a significant threat to the integrity of the power grid, necessitating comprehensive defense and mitigation strategies. Any anomaly in information and communication technology (ICT) should be detected for secure communications between devices in digital substations. This paper proposes large language models (LLMs), e.g., ChatGPT, for the cybersecurity of IEC 61850-based communications. Multi-cast messages such as generic object oriented system events (GOOSE) and sampled values (SV) are used for case studies. The proposed LLM-based cybersecurity framework includes, for the first time, data pre-processing of communication systems and human-in-the-loop (HITL) training (considering the cybersecurity guidelines recommended by humans). The results show a comparative analysis of detected anomaly data carried out based on the performance evaluation metrics for different LLMs. A hardware-in-the-loop (HIL) testbed is used to generate and extract a dataset of IEC 61850 communications.

ChatGPT↗
DOE OSTI2021

Explore Spatio‐Temporal Learning of Large Sample Hydrology Using Graph Neural Networks

Abstract Streamflow forecasting over gauged and ungauged basins play a vital role in water resources planning, especially under the changing climate. Increased availability of large sample hydrology data sets, together with recent advances in deep learning techniques, has presented new opportunities to explore temporal and spatial patterns in hydrological signatures for improving streamflow forecasting. The purpose of this study is to adapt and benchmark several state‐of‐the‐art graph neural network (GNN) architectures, including ChebNet, Graph Convolutional Network (GCN), and GraphWaveNet, for end‐to‐end graph learning. We explicitly represent river basins as nodes in a graph, learn the spatiotemporal nodal dependencies, and then use the learned relations to predict streamflow simultaneously across all nodes in the graph. The efficacy of the developed GNN models is investigated using the Catchment Attributes and MEteorology for Large‐sample Studies (CAMELS) data set under two settings, fixed graph topology (transductive learning), and variable graph topology (inductive learning), with the latter applicable to prediction in ungauged basins (PUB). Results indicate that GNNs are generally robust and computationally efficient, achieving similar or better performance than a baseline model trained using the long short‐term memory (LSTM) network. Further analyses are conducted to interpret the graph learning process at the edge and node levels and to investigate the effect of different model configurations. We conclude that graph learning constitutes a viable machine learning‐based method for aggregating spatiotemporal information from a multitude of sources for streamflow forecasting

Sun, Alexander Y.↗
DOE OSTI2026

Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas

Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

Electrical conductivity↗
DOE OSTI2024

On how structures convey non-diffusive turbulence spreading

Abstract We report on comprehensive experimental studies of turbulence spreading in edge plasmas. These studies demonstrate the relation of turbulence spreading and entrainment to intermittent convective density fluctuation events or bursts (i.e. blobs and holes). The non-diffusive character of turbulence spreading is thus elucidated. The turbulence spreading velocity (or mean jet velocity) manifests a linear correlation with the skewness of density fluctuations, and increases with the auto-correlation time of density fluctuations. Turbulence spreading by positive density fluctuations is outward, while spreading by negative density fluctuations is inward. The degree of symmetry breaking between outward propagating blobs and inward propagating holes increases with the amplitude of density fluctuations. Thus, blob-hole asymmetry emerges as crucial to turbulence spreading. These results highlight the important role of intermittent convective events in conveying the spreading of turbulence, and constitute a fundamental challenge to existing diffusive models of spreading.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Probabilistic simulation of quantum circuits using a deep-learning architecture

The fundamental question of how to best simulate quantum systems using conventional computational resources lies at the forefront of condensed matter and quantum computation. It impacts both our understanding of quantum materials and our ability to emulate quantum circuits. Here we present an exact formulation of quantum dynamics via factorized generalized measurements which maps quantum states to probability distributions with the advantage that local unitary dynamics and quantum channels map to local quasistochastic matrices. This representation provides a general framework for using state-of-the-art probabilistic models in machine learning for the simulation of quantum many-body dynamics. Using this framework, we have developed a practical algorithm to simulate quantum circuits using an attention network based on a powerful neural network ansatz responsible for the most recent breakthroughs in natural language processing. We demonstrate our approach by simulating circuits that build Greenberger-Horne-Zeilinger and linear graph states of up to 60 qubits, as well as a variational quantum eigensolver circuit for preparing the ground state of the transverse field Ising model on several system sizes. Our methodology constitutes a modern machine learning approach to the simulation of quantum physics with applicability both to quantum circuits as well as other quantum many-body systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗