Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Prior Distribution”

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 343 records · Page 19

Radiative dynamics and delayed emission in pure and doped InP/ZnSe/ZnS quantum dots

We have used time-correlated single photon counting to elucidate the radiative dynamics of InP/ZnSe/ZnS core/shell/shell quantum dots (QDs) that differ in the amount and distribution of excess indium. Stoichiometric QDs having an In:P atom ratio very near unity exhibit simple luminescence kinetics. The photoluminescence (PL) rises with the 40 ps instrument response function and exhibits a decay that is close to a single exponential with a time constant that decreases from 32 to 28 ns with increasing shell thickness. QDs having excess indium (In:P ratio of 1.15–1.63) show a significant component of a slower rise time assigned to transient population of indium-based hole traps in the ZnSe shell. They also have a slower PL decay, attributed to an equilibrium between these traps, which are optically dark, and the emissive valence-band state. This results in a radiative lifetime that increases from 32 to 48 ns with increasing shell thickness. Different treatments of the InP cores prior to shell deposition result in different core/shell interfaces as indicated by resonance Raman spectroscopy, as well as differences in the amplitude and timescale of the slow PL rise and the PL decay time. These are interpreted in terms of different radial distributions of the indium-based hole traps, which can be related to differences in the interfacial lattice strain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Convex Q-Learning in Continuous Time with Application to Dispatch of Distributed Energy Resources

Convex Q-learning is a recent approach to reinforcement learning, motivated by the possibility of a firmer theory for convergence, and the possibility of making use of greater a priori knowledge regarding policy or value function structure. This paper explores algorithm design in the continuous time domain, with a finite-horizon optimal control objective. The main contributions are (i) The new Q-ODE: a model-free characterization of the Hamilton-Jacobi-Bellman equation. (ii) A formulation of Convex Q-learning that avoids approximations appearing in prior work. The Bellman error used in the algorithm is defined by filtered measurements, which is necessary in the presence of measurement noise. (iii) Convex Q-learning with linear function approximation is a convex program. It is shown that the constraint region is bounded, subject to an exploration condition on the training input. (iv) The theory is illustrated in application to resource allocation for distributed energy resources, for which the theory is ideally suited.

Lu, Fan↗

Assessing the future of global energy-for-water

This study incorporates the energy demands of water abstraction, treatment, distribution, and post-use (wastewater) treatment into the Global Change Assessment Model (GCAM), an integrated human-Earth systems model, and analyzes a range of scenarios that estimate the future evolution of this demand of energy. The study builds on prior research on the historical demands of energy-for-water, fitting these demands of energy into existing nation-level energy balances, where demands for energy-for-water have heretofore been aggregated with other energy use in the commercial and public services, industrial, agricultural, and/or electric power sectors. This development allows for improved projections of future energy demands in general, and in this study, allows for assessment of the energy implications of different scenarios of improvements in water access and water quality that are consistent with the Sustainable Development Goals. In our baseline scenario, energy-for-water approximately doubles from 2010 to 2050, and in the scenario with the greatest level of water access, standards for treatment of wastewater, and irrigation expansion, the consequent energy-for-water demands increase three-fold from 2010 to 2050. Note also that these increases do not consider future scarcity-driven increases in energy requirements per unit of water abstraction, treatment, and distribution.

54 ENVIRONMENTAL SCIENCES↗

Early-time measure in eternal inflation

In a situation like eternal inflation, where our data is replicated at infinitely-many other space-time events, it is necessary to make a prior assumption about our location to extract predictions. The principle of mediocrity entails that we live at asymptotic late times, when the occupational probabilities of vacua has settled to a near-equilibrium distribution. In this paper we further develop the idea that we instead exist during the approach to equilibrium, much earlier than the exponentially-long mixing time. In this case we are most likely to reside in vacua that are easily accessed dynamically. Using first-passage statistics, we prove that vacua that maximize their space-time volume at early times have: 1. maximal ever-hitting probability; 2. minimal mean first-passage time; and 3. minimal decay rate. These requirements are succinctly captured by an early-time measure. The idea that we live at early times is a predictive guiding principle, with many phenomenological implications. First, our vacuum should lie deep in a funneled region, akin to folding energy landscapes of proteins. Second, optimal landscape regions are characterized by relatively short-lived vacua, with lifetime of order the de Sitter Page time. For our vacuum, this lifetime is ~ 10 130 years, which is consistent with the Standard Model estimate due to Higgs metastability. Third, the measure favors vacua with small, positive vacuum energy. This can address the cosmological constant problem, provided there are sufficiently many vacua in the entire ensemble of funnels. As a concrete example, we study the Bousso-Polchinski lattice of flux vacua, and find that the early-time measure favors lattices with the fewest number of flux dimensions. This favors compactifications with a large hierarchy between the lightest modulus and all other Kähler and complex structure moduli.

79 ASTRONOMY AND ASTROPHYSICS↗

Transport Estimations of Final State Interaction Effects on Short-range Correlation Studies Using the $(e,e'p)$ and $(e,e'pp)$ Reactions

Short range correlated (SRC) nucleon-nucleon pairs in nuclei are typically studied using measurements of electron-induced hard nucleon-knockout reactions (e.g. $(e,e'p)$ and $(e,e'pN)$), where the kinematics of the knocked-out nucleons are used to infer their initial state prior to the interaction. The validity of this inference relies on our understanding of the scattering reaction, most importantly how rescattering of the detected nucleons (final state interactions or FSI) distort their kinematical distributions. Recent SRC measurements on medium to heavy nuclei have been performed at high-$x_B$ (i.e., anti-parallel kinematics) where calculations of light nuclei indicate that such distortion effects are small. Here we study the impact of FSI on recent $^{12}$C$(e,e'p)$ and $^{12}$C$(e,e'pp)$ measurements using a transport approach. We find that while FSI can significantly distort the measured kinematical distributions of SRC breakup events, selecting high-$x_B$ anti-parallel events strongly suppresses such distortions. In addition, including the effects of FSI improves the agreement between Generalized Contact Formalism-based calculations and data and can help identify those observables that have minimal sensitivity to FSI effects. This result helps confirm the interpretation of experimental data in terms of initial-state momentum distributions and provides a new tool for the study of SRCs using lepton-scattering reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Underlying Roles of Polyol Additives in Promoting CO2 Capture in PEI/Silica Adsorbents

Solid-supported amines having low molecular weight branched poly(ethylenimine) (PEI) physically impregnated into porous solid supports are promising adsorbents for CO2 capture. Co-impregnating short-chain poly(ethylene glycol) (PEG) together with PEI alters the performance of the adsorbent, delivering improved amine efficiency (AE, mol CO2 sorbed / mol N) and faster CO2 uptake rates. To uncover the physical basis for this improved gas capture performance, we probed the distribution and mobility of the polymers in the pores via small angle neutron scattering (SANS), solid-state NMR, and molecular dynamic (MD) simulation studies. SANS and MD simulations reveal that PEG displaces wall-bound PEI, making amines more accessible for CO2 sorption. Solid-state NMR and MD simulation suggest intercalation of PEG into PEI domains, separating PEI domains and reducing amine-amine interactions, providing potential PEG-rich and amine-poor interfacial domains that bind CO2 weakly via physisorption while providing facile pathways for CO2 diffusion. Contrary to a prior literature hypothesis, no evidence is obtained for PEG facilitating PEI mobility in solid supports. Instead, the data suggest that PEG chains coordinate to PEI, form larger bodies with reduced mobility compared to PEI alone. We also demonstrate promising CO2 uptake and desorption kinetics at varied temperatures, given by favorable amine distribution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Real-time Implementation of Grid Code Compliant Grid Edge Energy Management System

Integrated distributed energy resources (DER) in a distribution system need to follow grid codes to avoid violations that result in DER/circuit segment disconnection. To comply with grid code requirements at the grid edge level, network constrained grid edge energy management system (EMS) can be deployed. The objective of grid edge EMS is to provide economic solution for active and reactive power DER setpoints at each dispatch interval and ensure voltage regulation to support secure interconnection of the grid edge segment to the distribution system with multiple inverter based DER units. In this work, real-time simulation of grid code compliant grid edge EMS is deployed in a realistic feeder circuit segment. For real-time simulation, communication between the grid edge EMS and DERs is done exploiting IEC 61850-7-420. It enables interoperability among different DERs and grid edge EMS. No prior art has deployed IEC 61850-7-420 GOOSE communication protocol for grid edge EMS. Conversion of IEC 61850 GOOSE messages to Modbus communication protocol is also performed to communicate with grid edge EMS in commodity-off the shelf embedded boards in this work. The real-time simulation in OPAL-RT real-time digital simulator shows the out-performance of grid edge EMS by reducing the voltage violation in the distribution circuit.

Energy management system↗

Computational studies of impurity migration during induction stirring of molten uranium

Understanding and controlling impurity behavior in actinide metal casting processes are foundational for efficient part production yet remain major challenges for researchers and industry. To help provide insight regarding impurity distribution during actinide metal casting, we have developed computational fluid dynamics (CFD) models for a laboratory-scale system using commercial and open-source code. Multiple simulation frameworks allow for improved confidence in the resulting outputs, while taking advantage of the maturity and convenience offered by commercial platforms and simultaneously maintaining the transparency and flexibility often provided by open-source software. Here, we describe multiple experiment-informed models designed to simulate a specific laboratory system in which uranium melt, containing a known starting concentration of carbon impurity, is electromagnetically stirred in an induction furnace. Here, the goal of the simulations is to predict the motion of uranium carbide microparticles in the velocity field of the melt. Prior to simulating the uranium-carbon system, numerical models were validated using a previously published nonradioactive experimental system. Effects of the size and shape of impurity particles in the models were investigated and agree with experimental findings. Simulation of smaller particles (< 50 µm) shows more homogenous distribution throughout the stirred melt. With increased particle size (100 µm), the body forces, which include the buoyancy force, dominate over the drag force, causing larger particles to move toward the crucible walls and upward in the system.

36 MATERIALS SCIENCE↗

Cesium Ion Exchange Testing Using Crystalline Silicotitanate with Hanford Tank Waste 241-AP-107

At the time of this testing, the Low-Activity Waste Pretreatment System (LAWPS) was to provide for the initial production of immobilized low-activity waste by feeding Hanford tank supernate from tank farms to the Hanford Tank Waste Treatment and Immobilization Plant (WTP) Low-Activity Waste (LAW) Facility for immobilization. Washington River Protection Solutions requested that Hanford tank waste collected from tank 241-AP-107 (hereafter called AP-107) be processed using conceived pretreatment steps (suspended solids removal by filtration, Cs removal by ion exchange) then vitrified. A small-scale test platform to demonstrate the solids filtration, Cs removal, and LAW vitrification was constructed and installed at Pacific Northwest National Laboratory. Bench-scale ion exchange testing with approximately 9 L of AP-107 supernate was conducted using crystalline silicotitanate (CST) ion exchange media. The IONSIV R9140-B CST was provided by Honeywell UOP, LLC in 2018 (Batch 2081000057). The ion exchange media was first tested with simulant and was previously described. This report describes the Cs ion exchange batch contact and column test results with the AP-107 tank waste. Batch contact testing helps to evaluate CST performance on tank waste supernate prior to processing it in the ion exchange columns. Batch contacts were performed with the waste at four Cs concentrations at a phase ratio of 200 (liquid volume to exchanger mass) with AP-107. The distribution coefficient (K d ) at the equilibrium condition of 8.57 µg Cs/mL (AP-107 feed condition) was determined to be 669 mL AP-107/g CST. With a CST bed density of 1.00 g/mL, this K d corresponded to a predicted 50% Cs breakthrough of 669 bed volumes (BVs). The Cs load capacity at the equilibrium feed condition was determined to be 7.5 mg Cs/g dry CST. The column testing was prototypic to the intended LAWPS operations in a lead-lag column format, although on a small-scale basis with 10-mL CST beds. The feed was processed downflow through the lead column and then through the lag column at ~2.2 BV/h. Loading continued until the lag column reached the WTP waste acceptance criteria (WAC) for receiving supernatant waste for vitrification (a function of the Na and 137 Cs concentrations). For AP-107, the WAC is 0.114% of the influent 137 Cs concentration; this required a Cs decontamination factor of 876. The Cs effluent from the lag column reached the WAC after processing ~410 BVs. To keep the subsequent product effluent below the WAC, a replacement lag column was prepared, the lead column was removed from service (after processing a total of 471 BVs), the lag column was put into the lead column position, and the replacement lag column was installed. Feed processing continued and after another ~290 BVs the Cs effluent from the lag column again exceeded the WAC. In both cases, the lead columns only reached 25% Cs breakthrough before removal. Although 50% Cs breakthrough was not reached, this value was estimated and averaged based on extrapolation of the loading curves (640 BVs) and agreed within 4% of the predicted 50% Cs breakthrough from batch contact test results (669 BVs). Table ES.1 summarizes the observed column performance and relevant Cs loading characteristics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Particle scale impacts on deconstruction energy of pine residues

The goal of this Case Study was to quantify the impacts of variable moisture and ash on hammer mill throughput and energy consumption and on generation of fines that are not able to be fed to conversion, as compared to a status quo Base Case system. Also considered was convertible carbon content (minimum carbon specification) and maximum ash content and the delivered feedstock cost impacts of not being able to feed residue not meeting both specifications to the conversion reactor. Laboratory data on the impacts of input particle size and moisture content on the exit particle size were received from FCIC Subtask 5.2 from their single particle impact population balance modeling study. Additional throughput and energy consumption data were obtained from FCIC Subtask 5.2 for the same grinder with a 6 mm screen in place. These data were utilized to develop the necessary response surface equations to perform throughput analysis using discrete event simulation. Because the ash contents in the separated fines had not been analyzed in the laboratory at the time of the model runs, we chose to assume that the ash distributed proportionally with total mass into the overs and unders in the disk screen following grinding. Key takeaways from this Case Study are that it is significantly more cost effective to hammer mill the residue prior to drying, even though the grinder throughput is lower and energy consumption is higher versus drying first before grinding. An effect of dry grinding versus high moisture grinding is the production of higher amounts of fines during dry grinding, leading to significantly more of the ground feedstock being rejected by conversion for being below a minimum particle size. With wet grinding the system is still able to produce more preprocessed feedstock meeting the minimum particle size specification even though the instantaneous throughput is lower than for the case of grinding dry feedstock. Additionally, even without the higher fines production from dry grinding, the status quo would still be more costly than wet grinding because the material is rejected after the drying energy has already been input for the dry grinding case. Finally, significant reductions in drying energy are obtained by drying after grinding, and those reductions are of far greater magnitude than the grinding energy increase.

09 BIOMASS FUELS↗

Influence of Powder Characteristics and Processing Methods on Creep Performance of Powder Metallurgy Hot Isostatic Pressed SS316

The U.S. nuclear energy expansion goals are driving the demand for manufacturing routes that can rapidly produce large, complex, near net shape components. Powder metallurgy hot isostatic pressing (PM HIP) is an advanced manufacturing technique that can be economically scaled-up, while alleviating the supply chain challenges that forging and casting face in terms of cost and lead time constraints. This makes PM-HIP a viable technology to aid and accelerate large-scale part manufacturing for nuclear applications. However, large-scale qualification and deployment of this technology require a thorough understanding of the influence of powder feedstock quality, powder handling history, hot isostatic pressing (HIP) parameters, and subsequent heat treatment on microstructural evolution and elevated temperature mechanical performance. The present work focusses on 316 austenitic stainless steel (SS316) which is most commonly used in high temperature environments for nuclear applications Results from this study show that PM HIPed SS316 meets ASME tensile requirements at room temperature and at elevated temperature. However, creep performance of PM-HIPed SS316 remains inferior to its wrought counterpart, demonstrating that tensile performance alone is not a reliable metric for long duration high temperature integrity. Further, this report delineates powder derived microstructural features that govern creep damage, with key evidences pointing to “microstructural inheritance” from gas atomized powder feedstocks. Commercial SS316 powders of varying chemical compositions and recycling histories were studies, and the results showed large differences in elemental segregation, oxide surface layers and secondary phase distributions. Multi-scale characterization revealed segregation of chromium, molybdenum, manganese and silicon at the boundaries and the precipitation of manganese-, silicon-, and molybdenum-oxides. During HIP consolidation, these surface oxides transform into decorated prior particle boundaries (PPBs) and grain boundary inclusions that persist through conventional post-HIP solution annealing treatment. The retained oxides in post-HIP microstructures were found to influence grain growth, precipitation behavior, and ultimately creep cavitation and fracture. Such post-HIP heat treatments are therefore limited by a complex trade-off between grain growth, and oxide coarsening which aggravate creep damage by acting as nucleation sites for cavities. The objective of this work is to establish an integrated processing–structure–property framework for PM-HIP 316 stainless steel by investigating the influence of powder feedstock characteristics in pre- and post-HIP processing as well as to understand the significance of post-HIP heat treatment on microstructural evolution and creep properties. The results from this report emphasize the significance of powder feedstock integrity in improving creep performance of PM-HIPed SS316, by highlight the effect of rapid solidification, elemental segregation, oxide formation and powder recycling on microstructural defect inheritance following HIP consolidation. Rather than considering HIP processing, solution annealing, and mechanical performance independently, this report treats powder production, HIP consolidation, post-HIP thermal processing, and creep deformation as interconnected stages within a continuous metallurgical process. The resulting framework will provide a scientific basis for developing feedstock engineering strategies capable of improving long-term reliability of PM-HIP stainless steels and accelerating their qualification for advanced nuclear applications.

Ajjarapu, Pavan [Oak Ridge National Laboratory (OR↗

Time-Sensitive Quantum Key Distribution (TSQKD)

This is the final report for the Time-Sensitive Quantum Key Distribution (TSQKD) project and summarizes the demonstration and longer-term field testing at an electric utility, EPB of Chattanooga, in Chattanooga, TN. The demonstration was held virtually (due to the COVID pandemic) and successfully executed the tests outlined in a prior report, “Task 10: Develop a QKD Time-Sensitive Networking Test Suite” submitted 9/30/2020. The portable test equipment explained in Task 10 and earlier task reports was tested at GE Research remotely connected with Qubitekk’s QKD equipment in San Diego, CA. Then the equipment was shipped to Qubitekk for physical integration with Qubitekk’s Quantum DataLoc™ Key Server 1570nm (http://qubitekk.com/products/qkd-for-industrial-control-systems-ics/). The portable test equipment was shipped from Qubitekk to EPB for demonstration and extended field testing. A summary of the demonstration and test results are described in this report. The demo and extended field test results are enabled by the NETCONF/YANG IEEE P1913 data management model developed as part of this project (Task 5). The NETCONF/YANG IEEE P1913 data management model demonstrated significant practical utility in providing a common means of monitoring and configuring the QKD system. All data in this report comes from the IEEE P1913 data management model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Widespread co‐occurrence of Sebacinales and arbuscular mycorrhizal fungi in switchgrass roots and soils has limited dependence on soil carbon or nutrients

Societal Impact Statement This work addresses a novel group of Sebacinales mycorrhizal fungi being considered for development as inoculants in managed ecosystems because of their potential access to nutrients in soil organic matter. By comparing the diversity, distributions, and abundances of the Sebacinales with the more common arbuscular mycorrhizal fungi in switchgrass, a biofuel crop, we demonstrate that current suggestions for a Sebacinales revolution in agriculture should be tempered by their ecology. In particular, Sebacinales are rare compared to arbuscular mycorrhizal fungi, and are only weakly associated with soil carbon, suggesting that ideas about improved soil carbon cycling associated with Sebacinales need to be sufficiently studied across a range of environmental conditions prior to their consideration for broad‐spectrum soil inoculants. Summary Arbuscular mycorrhizal (AM) fungi are widespread and important root symbionts, but recent work suggests that Sebacinales fungi may play an equally important role in both plant success and ecosystem carbon and nutrient cycling based on their worldwide occurrence and putative access to organic matter. However, the ecological impacts of Sebacinales will depend on their abundance and distribution relative to AM fungi and environmental soil carbon and nutrient patterns, which remain unexplored. We characterized Glomeromycota and Sebacinales fungi in switchgrass ( Panicum virgatum L.) roots and soils across 14 sites with diverse soil conditions. We examined group richness differences, co‐occurrence patterns, and how the relative abundance of these fungi related to soil carbon, nutrient stoichiometry, and host size. Sebacinales were widespread, but less diverse, common, and abundant than Glomeromycota. Moreover, co‐occurrences were predominantly random, suggesting relatively few interactions between these groups. Sebacinales increased relative to Glomeromycota in soils with more carbon, but explanatory power was limited. Based on our findings, we suggest that Sebacinales are likely complementary to AM fungi in roots. Expectations that Sebacinales have large effects on soil carbon and nutrient cycling may need to be reconsidered, at least based on their limited abundances relative to AM fungi in switchgrass. This is an important consideration as Sebacinales are candidates for use as inoculants in managed ecosystems.

Lee, Marissa R.↗

The PAU Survey: an improved photo- z sample in the COSMOS field

Here we present – and make publicly available – accurate and precise photometric redshifts in the ACS footprint from the COSMOS field for objects with i AB ≤ 23. The redshifts are computed using a combination of narrow-band photometry from PAUS, a survey with 40 narrow bands spaced at $100\,\mathring{\rm A}$ intervals covering the range from 4500 to $8500\,\mathring{\rm A}$, and 26 broad, intermediate, and narrow bands covering the UV, visible and near-infrared spectrum from the COSMOS2015 catalogue. We introduce a new method that models the spectral energy distributions as a linear combination of continuum and emission-line templates and computes its Bayes evidence, integrating over the linear combinations. The correlation between the UV luminosity and the $\mathrm{O\,{\small II}}$ line is measured using the 66 available bands with the zCOSMOS spectroscopic sample, and used as a prior which constrains the relative flux between continuum and emission-line templates. The flux ratios between the $\mathrm{O\,{\small II}}$ line and H α , H β and $\mathrm{O\,{\small III}}$ are similarly measured and used to generate the emission-line templates. Comparing to public spectroscopic surveys via the quantity Δ z ≡ (z photo – z spec )/(1 + z spec ), we find the photometric redshifts to be more precise than previous estimates, with σ 68 (Δ z ) ≈ (0.003, 0.009) for galaxies at magnitude i AB ~ 18 and i AB ~ 23, respectively, which is three times and 1.66 times tighter than COSMOS2015. Additionally, we find the redshifts to be very accurate on average, yielding a median of the Δ z distribution compatible with |median(Δ z )| ≤ 0.001 at all redshifts and magnitudes considered. Both the added PAUS data and new methodology contribute significantly to the improved results. The catalogue produced with the technique presented here is expected to provide a robust redshift calibration for current and future lensing surveys, and allows one to probe galaxy formation physics in an unexplored luminosity-redshift regime, thanks to its combination of depth, completeness, and excellent redshift precision and accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗

Powder Metallurgy HIP Process Study and Mechanical Property Evaluations for IN740H

The effects of hot isostatic pressing (HIP) parameters and powder size distribution (PSD) on the resulting microstructure and properties of powder metallurgy (P/M) IN740H have been investigated. Properties can be substantially tailored through the appropriate management of these parameters. Applying a thermal presoak to the input powder prior to the sintering stage of HIP can significantly increase average grain size and coarse grain fraction, enabling a trade for improved creep resistance at the expense of tensile strength. Here, powder thermal presoak was found to have a larger impact on microstructure and properties than the use of a coarser PSD alone. Additionally, substantial grain growth and grain boundary migration is observed in P/M IN740H at HIP temperatures above 1260°C. Given its inherently finer grain size, P/M HIP IN740H typically outperforms its cast and wrought counterparts in tensile response. Creep life at 700–800°C was inferior to wrought due to PPBs coincident with grain boundaries which served as cavity nucleation sites, leading to lower creep ductility. The ability to adjust the microstructure toward a more balanced creep and tensile behavior allows P/M HIP IN740H to be considered for high-temperature applications in advanced ultra-supercritical steam and supercritical carbon dioxide power cycles.

36 MATERIALS SCIENCE↗

Normalizing flows for likelihood-free inference with fusion simulations

Fluid-based scrape-off layer transport codes, such as UEDGE, are heavily utilized in tokamak analysis and design, but typically require user-specified anomalous transport coefficients to match experiments. Determining the uniqueness of these parameters and the uncertainties in them to match experiments can provide valuable insights to fusion scientists. Here, we leverage recent work in the area of likelihood-free inference (‘simulation-based inference’) to train a neural network, which enables accurate statistical inference of the anomalous transport coefficients given experimental plasma profile input. UEDGE is treated as a black-box simulator and runs multiple times with anomalous transport coefficients sampled from priors, and the neural network is trained on these simulations to emulate the posterior. The neural network is trained as a normalizing flow model for density estimation, allowing it to accurately represent complicated, high-dimensional distribution functions. With a fixed simulation budget, we compare a single-round procedure to a multi-round approach that guides the training simulations toward a specific target observation. Finally, we discuss the future possibilities for use of amortized models, which train on a wide range of simulations and enable fast statistical inference for results during experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantum algorithms from fluctuation theorems: Thermal-state preparation

Fluctuation theorems provide a correspondence between properties of quantum systems in thermal equilibrium and a work distribution arising in a non-equilibrium process that connects two quantum systems with Hamiltonians H 0 and H 1 = H 0 + V. Building upon these theorems, we present a quantum algorithm to prepare a purification of the thermal state of H 1 at inverse temperature β ≥ 0 starting from a purification of the thermal state of H 0 . The complexity of the quantum algorithm, given by the number of uses of certain unitaries, is $\mathcal{O}$ (e β(ΔA - w l )/2 ), where ΔA is the free-energy difference between H 1 and H 0 , and w l is a work cutoff that depends on the properties of the work distribution and the approximation error ϵ > 0. If the non-equilibrium process is trivial, this complexity is exponential i β∥V∥, where ∥V∥ is the spectral norm of V. This represents a significant improvement of prior quantum algorithms that have complexity exponential in β∥H 1 ∥ in the regime where ∥V∥$\ll$ ∥H 1 ∥. The dependence of the complexity in ϵ varies according to the structure of the quantum systems. It can be exponential in 1/ϵ in general, but we show it to be sublinear in 1/ϵ if H 0 and H 1 commute, or polynomial in 1/ϵ if H 0 and H 1 are local spin systems. The possibility of applying a unitary that drives the system out of equilibrium allows one to increase the value of w l and improve the complexity even further. To this end, we analyze the complexity for preparing the thermal state of the transverse field Ising model using different non-equilibrium unitary processes and see significant complexity improvements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗