Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal equilibrium”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Atmospheric-pressure ammonia synthesis on AuRu catalysts enabled by plasmon-controlled hydrogenation and nitrogen-species desorption

The Haber–Bosch process for ammonia synthesis contributes up to ~3% of global greenhouse gas emissions. Plasmonic catalysts strongly concentrate light and can alter the reaction intermediates via out-of-equilibrium processes, providing the potential for an alternative, less-energy-intensive pathway to synthesize ammonia. Here, in this study, we show that gold-ruthenium (AuRu) bimetallic nanoparticles can synthesize ammonia at room temperature and pressure using visible light. We create AuRu alloys with varying compositions and achieve ammonia production rates of ~60 μmol per gram of catalyst bed per hour. In situ infrared spectroscopy reveals that light accelerates the hydrogenation of nitrogen intermediates compared to conventional thermal catalysis. Through computational modelling, we demonstrate that photo-excited electrons enable associative hydrogenation pathways for nitrogen activation rather than direct nitrogen–nitrogen bond breaking. This light-assisted mechanism requires both hydrogen and light working together to overcome the nitrogen activation barrier, mimicking how biological enzymes produce ammonia naturally and providing fundamental insights for developing sustainable, energy-efficient chemical synthesis.

Yuan, Lin [Stanford Univ., CA (United States)] (OR↗
DOE OSTI2020

Electronic structure and optical properties of quantum crystals from first principles calculations in the Born–Oppenheimer approximation

In this work, we develop a formalism to accurately account for the renormalization of the electronic structure due to quantum and thermal nuclear motions within the Born–Oppenheimer approximation. We focus on the fundamental energy gap obtained from electronic addition and removal energies from quantum Monte Carlo calculations in either the canonical or grand-canonical ensembles. The formalism applies as well to effective single electron theories such as those based on density functional theory. We show that the electronic (Bloch) crystal momentum can be restored by marginalizing the total electron–ion wave function with respect to the nuclear equilibrium distribution, and we describe an explicit procedure to establish the band structure of electronic excitations for quantum crystals within the Born–Oppenheimer approximation. Based on the Kubo–Greenwood equation, we discuss the effects of nuclear motion on optical conductivity. Our methodology applies to the low temperature regime where nuclear motion is quantized and, in general, differs from the semi-classical approximation. We apply our method to study the electronic structure of C2/c-24 crystalline hydrogen at 200 K and 250 GPa and discuss the optical absorption profile of hydrogen crystals at 200 K and carbon diamond at 297 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Pions in hot dense matter and their astrophysical implications

Here we study the role of pions in hot dense matter encountered in astrophysics. We find that strong interactions enhance the number density of negatively charged pions and that this enhancement can be calculated reliably for a relevant range of density and temperature by using the virial expansion. We assess the influence of pions and muons on the equation of state (EOS) and weak-interaction rates in hot dense matter. We find that thermal pions increase the proton fraction and soften the EOS. We also find that charged current weak reactions involving pions and muons $\bar{ν}_ μ + μ^– → π$ and $ν_μ + π^– → μ^–$ make an important contribution to the opacity of muon neutrinos. This could influence the dynamics of core-collapse supernovae and neutron-star mergers. Finally, we note that pion-nucleon reactions can alter the evolution of the proton fraction when weak interactions are not in equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Observation of conformational dynamics in single light-harvesting proteins from cryptophyte algae

Photosynthetic organisms use pigment–protein complexes to capture the sunlight that powers most life on earth. Within these complexes, the position of the embedded pigments is all optimized for light harvesting. At the same time, the protein scaffold undergoes thermal fluctuations that vary the structure, and, thus, photophysics, of the complexes. While these variations are averaged out in ensemble measurements, single-molecule spectroscopy provides the ability to probe these conformational changes. We used single-molecule fluorescence spectroscopy to identify the photophysical substates reflective of distinct conformations and the associated conformational dynamics in phycoerythrin 545 (PE545), a pigment–protein complex from cryptophyte algae. Rapid switching between photophysical states was observed, indicating that ensemble measurements average over a conformational equilibrium. A highly quenched conformation was also identified, and its population increased under high light. This discovery establishes that PE545 has the characteristics to serve as a photoprotective site. Finally, unlike homologous proteins from the evolutionarily related cyanobacteria and red algae, quenching was not observed upon photobleaching, which may allow for robust photophysics without the need for rapid repair or replacement machinery. Collectively, these observations establish the presence of a rich and robust set of conformational states of PE545. Cryptophytes exhibit particularly diverse energetics owing to the variety of microenvironments in which they survive, and the conformational states and dynamics reported here may provide photophysical flexibility that contributes to their remarkable ability to flourish under diverse conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Characterization of compositional variation and solidification condition on rapid solidification behavior of beta stabilized titanium

Rapid solidification (RS) is a field of study that has increased in importance due to the rise of additive manufacturing (AM) techniques. The breakdown of equilibrium at the solid-liquid interface results in an inability to use traditional solidification models that have been applied to great success in casting. As such, there is a current need for a high-throughput testing method to screen alloys at rapid solidification conditions without the significant cost for powder atomization. In this work, twin hammer splat quenching (SQ) is utilized and theoretically developed as a technique uniquely placed to test alloy systems for AM applications. Using analytical thermal spray models and numerical heat transfer simulations, the cooling rate and solidification conditions of splat quenching are compared across several parameters to evaluate the key parameters for generating cooling rates relevant to AM process. Once a framework for the cooling rates and relevancy of SQ are developed, SQ was used to process alloy compositions in the Ti5553 alloy space (AMS 7026) to study the formation of detrimental secondary phases occurring upon rapid quenching from the melt. A martensitic phase is found to occur at the lower end of the standard, but the composition listed in the standard were found to be safe. Finally, the rapid solidification behavior of Ti-5553 is evaluated by comparing splat quenching, laser welding, and plasma atomization that offer a wide window of solidification conditions. The segregation and solidification structure are analyzed and compared to models to provide a deeper understanding of the solidification behavior of beta stabilized titanium alloys.

Williamson, Christopher J.↗
DOE OSTI2020

Microstructural modification of additively manufactured metals by electropulsing

Additive manufacturing (AM) promises rapid development cycles and fabrication of ready-to-use, geometrically-complex parts. The metallic parts produced by AM often contain highly non-equilibrium microstructures, e.g. chemical microsegregation and residual dislocation networks. While such microstructures can enhance some material properties, they are often undesirable. Many AM parts are thus heat-treated after fabrication, a process that significantly slows production. This study investigated if electropulsing, the process of sending high-current-density electrical pulses through a metallic part, could be used to modify the microstructures of AM 316 L stainless steel (SS) and AlSi10Mg parts fabricated by selective laser melting (SLM) more rapidly than thermal annealing. Electropulsing has shown promise as a rapid postprocessing method for materials fabricated using conventional methods, e.g. casting and rolling, but has never been applied to AM materials. For both the materials used in this study, as-fabricated SLM parts contained significant chemical heterogeneity, either chemical microsegregation (316 L SS) or a cellular interdendritic phase (AlSi10Mg). In both cases, annealing times on the order of hours at high homologous temperatures are necessary for homogenization. Using electropulsing, chemical microsegregation was eliminated in 316 L SS samples after 10, 16 ms electrical pulses. In AlSi10Mg parts, electropulsing produced spheroidized Si-rich particles after as few as 15, 16 ms electrical pulses with a corresponding increase in ductility. Finally, this study demonstrated that electropulsing can be used to modify the microstructures of AM metals.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Uncovering grain and subgrain microstructure at the scale of additive manufacturing melt tracks with a scalable cellular automaton solidification model

Metal additive manufacturing, characterized by rapid solidification, yields refined grains with a distinctive cellular subgrain microstructure that plays a pivotal role in determining material properties. Due to the significant computational expense demanded to simulate the required physics with submicron spatial resolution, their numerical simulations have been limited to proof-of-concept studies to either 2D or small subregions of a melt pool. In this study, an open-source, scalable, solidification code, muMatScale, based on the cellular automaton method, has been developed to predict the grain and the underlying subgrain microstructure over an entire melt pool. The model incorporates flexible parallelization schemes, utilizing MPI and OpenMP GPU Offloading, in addition to appropriate multi-physics specific to non-equilibrium rapid solidification in AM. The impact of nucleation parameters on grain microstructures was investigated with a focus on grain size variations and morphology transitions. With selected nucleation parameters, the simulation predicted the grain size, subgrain morphology, crystallographic orientation, and microsegregation aligned with experimental measurements. The model demonstrates that epitaxial grain growth is a dominant factor at the melt pool boundary, influencing grain size variation under different grain sizes in the build plate while maintaining consistent primary dendrite arm spacing under identical thermal conditions. Here, the highly efficient numerical model enables large-scale simulations with a spatial resolution of 100 nm or less, unveiling unprecedented insights into thermal and solutal diffusion driven grain growth, and the subgrains with microsegregation within grains in 3D across scales. muMatScale will enable the linking of submicron length-scale microstructure to part-level material behavior by investigating fundamental solidification problems at the intercellular scale in many-track and many-layer builds.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Comprehensive Characterization of Water Group Ion Composition and Distributions in Saturn's Magnetosphere With Cassini Plasma Spectrometer Data

Saturn's magnetosphere is continuously supplied with neutrals from the Enceladus plume and the icy rings, which undergo ionization and charge-exchange to form a complex water-group plasma environment. While the Cassini Plasma Spectrometer (CAPS) instrument has provided extensive compositional information, detailed separation of individual water-group ion species in time-of-flight (TOF) data has not previously been achieved. In this study, we perform forward modeling of CAPS-IMS energy-per-charge (E/Q) and TOF spectra obtained between 2004 and 2012 to resolve O + , OH + , H 2 O + , and H 3 O + and to characterize their plasma properties, including number density, temperature, and thermodynamic κ. Our results demonstrate that O + is the dominant thermal ion species throughout Saturn's magnetosphere, comprising up to ∼70% of the total ion population beyond ∼5 Saturn radii (R S ). In contrast, molecular ions such as OH + , H 2 O + , and H 3 O + dominate closer to Enceladus but rapidly dissociate into atomic ions between ∼5 and 10 R S . This radial region is also characterized by the steepest increase in plasma flow speed, which rises from ∼40% to ∼80% of rigid corotation. Simultaneously, ion velocity distributions approach Maxwell–Boltzmann equilibrium, as indicated by high kappa values. These findings provide new constraints on the ion–neutral chemistry that regulates the balance between molecular and atomic ions in Saturn's magnetosphere. They also emphasize the critical role of the 5–10 R S region as a transition zone for both plasma composition and dynamics. Our results refine previous CAPS-based studies and underscore the need to incorporate seasonal variability and ionospheric coupling into future global models of Saturn's plasma environment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Isothermal Microcalorimetric Analysis of Li/CF x Battery Discharge

Lithium/carbon monofluoride (Li/CF x ) batteries generate substantial amounts of heat during discharge in part due to the large deviation of the loaded voltage from the measured and thermodynamically predicted open circuit voltage. Here, in this study, we further analyze this system by estimating the equilibrium voltage (V eq ) and its temperature dependence (dV eq /dT) over the entire discharge range. Based on these results, the ohmic and entropic heat contributions to the overall heat flow are calculated from experimental data. The ohmic heat flow is consistent with an electrochemical impedance spectroscopy model and electrode density measurements, indicating decreasing porosity of the electrode during discharge. Entropic heating increases during discharge as the cell reaction becomes increasingly entropically unfavorable. The total energy dissipated by the cell (electrical and thermal) remains similar over the entire discharge, evidence for a two-phase reaction with a constant, and rate-independent, ΔG rxn of about -465 kJ mol -1 . However, near the end of discharge (>90% DOD), ΔG rxn changes significantly indicating possible secondary reactions.

25 ENERGY STORAGE↗
DOE OSTI2023

FES 2023 Theory & Simulation Performance Target

The Report for the 2023 Theory & Simulation Performance Target (TSPT) is presented herein. The 1st, 2nd and 3rd Quarter Reports are given in Sec. I, Sec. II and Sec. III. These are unchanged from those previously submitted, except that in this final document the equations, figures and citations have been automatically renumbered. The 4th Quarter Report is given in Sec. IV. The following material presents verification calculations of several widely used three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium codes in strongly shaped configurations; for both fixed and free-boundary equilibria; for the so-called vacuum case and for the case of high plasma pressure; and, for magnetic fields that have nested magnetic surfaces and for fields that have a magnetic island. The configurations studied includes the recently discovered “precise-QA” case [Phys. Rev. Lett., 128, 035001 (2002)] and other highly optimized configurations. Because of the strong (extreme) shaping, the latter of these pose challenges to the 3D MHD codes. A new method is introduced for constructing coordinates in arbitrarily shaped toroidal domains, and this can improve the robustness of 3D equilibrium calculations. Beyond the equilibrium calculation, this report demonstrates the capability of performing non-linear extended-MHD modelling of optimized stellarator equilibria. This capability allows the assessment not only of linear stability thresholds but also of the nonlinear evolution of unstable equilibria to new configurations that might be far from the original configuration, with significantly different transport properties and magnetic topology. They also enable the self-consistent determination of pressure profiles in non-integrable magnetic topology through the use of highly anisotropic thermal-transport equations. The calculations presented here represent the first demonstration of extended-MHD calculations in strongly-shaped, optimized quasi-symmetric configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Oxygenate Onboard Separation for Octane-on-Demand

This project further examines the use of Self-Assembled Monolayers on Mesoporous Supports ® SAMMS ® -based sorbent materials as a sorbent for alcohols from alcohol-gasoline blends in the context of an onboard separation approach for use in an octane-on-demand strategy. Several questions were posed by potential industry collaborators seeking to better understand how the SAMMS®-based materials would perform in a more realistic environment. Several conclusions can be made from the work conducted here, with the caveat that these experiments do not represent the results that would be obtained from continuous or long-term use of the sorbent, because of the short duration of the project. Vigorous extractions into warm gasoline did not reveal the presence of additional species in gas chromatographic analysis. Vibration testing for up to eight hours under aggressive conditions did not show particle attrition. Thermal desorption experiments showed that the SAMMS ® have a higher capacity for methanol, approximately 50 weight-percent of the sorbent, than for ethanol, approximately 20 weight-percent of the sorbent, and that the methanol is easier to extract. Testing of the A20 fuel blends was insufficient and requires a slightly more sophisticated approach than was attempted here. Additionally, further work would be needed to assess the rate at which the alcohol is absorbed into the sorbent. The testing conducted here suggests that equilibrium is reached in well under an hour. While this study provides additional insights into the use of SAMMS®-based sorbent materials for onboard alcohol separation, there is room for further work employing a benchtop testing apparatus similar to that described herein.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2020

Simple non-perturbative resummation schemes beyond mean-field II: Thermodynamics of scalar φ4 theory in 1 + 1 dimensions at arbitrary coupling

Recently, non-perturbative approximate solutions were presented that go beyond the well-known mean-field resummation. In this work, these non-perturbative approximations are used to calculate finite-temperature equilibrium properties for scalar [Formula: see text] theory in two dimensions such as the pressure, entropy density and speed of sound. Unlike traditional approaches, it is found that results are well-behaved for arbitrary temperature/coupling strength, are independent of the choice of the renormalization scale [Formula: see text], and are apparently converging as the resummation level is increased. Results also suggest the presence of a possible analytic cross-over from the high-temperature to the low-temperature regime based on the change in the thermal entropy density.

Astronomy & Astrophysics↗
DOE OSTI2024

Quantum Thermodynamics of Nonequilibrium Processes in Lattice Gauge Theories

A key objective in nuclear and high-energy physics is to describe nonequilibrium dynamics of matter, e.g., in the early Universe and in particle colliders, starting from the standard model of particle physics. Classical computing methods, via the framework of lattice gauge theory, have experienced limited success in this mission. Quantum simulation of lattice gauge theories holds promise for overcoming computational limitations. Because of local constraints (Gauss’s laws), lattice gauge theories have an intricate Hilbert-space structure. This structure complicates the definition of thermodynamic properties of systems coupled to reservoirs during equilibrium and nonequilibrium processes. We show how to define thermodynamic quantities such as work and heat using strong-coupling thermodynamics, a framework that has recently burgeoned within the field of quantum thermodynamics. Our definitions suit instantaneous quenches, simple nonequilibrium processes undertaken in quantum simulators. To illustrate our framework, we compute the work and heat exchanged during a quench in a Z 2 lattice gauge theory coupled to matter in 1+1 dimensions. Here, the thermodynamic quantities, as functions of the quench parameter, evidence a phase transition. For general thermal states, we derive a simple relation between a quantum many-body system’s entanglement Hamiltonian, measurable with quantum-information-processing tools, and the Hamiltonian of mean force, used to define strong-coupling thermodynamic quantities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Discrete stochastic model of point defect-dislocation interaction for simulating dislocation climb

Dislocation climb is an important high temperature process in metals plasticity, responsible for the phenomena such as creep, swelling, or hardening. Climb is defined by the ability of dislocations to leave their original glide plane by interacting with point defects. As such, dislocation climb is controlled by point defect diffusion/absorption/emission, all of which involve thermal activation. The existing thermodynamically consistent models for climb are generally formulated in a continuum framework, through the definition of effective defect fluxes and climb propensities in response to thermodynamic driving forces. However, the point-wise discrete nature of vacancies (and/or self-interstitials) confers a highly discrete nature to the climb dynamics, which is also strongly affected by elastic forces. The combination of discreteness, thermal activation, and elasticity makes this process too challenging for direct atomistic methods such as molecular dynamics. Here we develop a kinetic Monte Carlo model that captures vacancy generation and transport kinetics acting in conjuction with the evolving elastic fields provided by discrete dislocation dynamics simulations. The two models are coupled via the applied stresses and stress gradients generated by dislocation structures at vacancy locations. Our simulations reveal two surprising results. First that climb is dominated by vacancy emission even when the background vacancy concentration is much higher than the equilibrium one. And, second, that climb velocities might be much faster than otherwise believed when one uses the classical theories of climb. These effects are due to the locality of vacancy-dislocation processes, which are not captured in classical treatments that assume smooth vacancy fluxes and homogeneous concentrations. We apply the method to study elementary climb processes in body-centered cubic iron and furnish climb mobility functions to be used in parametric dislocation dynamics and/or crystal plasticity simulations. Here, we apply the technique to study non-conservative plastic bypass of spherical precipitates by edge dislocations and point out the differences between our discrete approach and existing continuum formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Applications of Thermochemical Modeling in Molten Salt Reactors

The extensively evaluated and consistent thermodynamic database, the Molten Salt Thermal Properties Database—Thermochemical (MSTDB-TC), was used along with additional thermodynamic values from other sources as examples of ways to examine molten salt reactor (MSR) fuel behavior. Relative stability with respect to halide potential and temperature for likely fuel and fission product components were mapped in Ellingham diagrams for the chloride and fluoride systems. The Ellingham diagrams provide a rich, visual means for identifying halide-forming components in proposed fuel/solvent salt systems. Thermochemical models and values from MSTDB-TC and ancillary sources were used in global equilibrium calculations to provide compositions for a close analysis of the behavior of a possible Molten Chloride Salt Fast Reactor and a Molten Salt Reactor Experiment-type system at high burnup (100 GWd/t). The results illustrated the oxidative nature of burnup in MSRs and provided information about redox behavior and possible control.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Asymptotic reversibility of thermal operations for interacting quantum spin systems via generalized quantum Stein’s lemma

Abstract For quantum spin systems in any spatial dimension with a local, translation-invariant Hamiltonian, we prove that asymptotic state convertibility from a quantum state to another one by a thermodynamically feasible class of quantum dynamics, called thermal operations, is completely characterized by the Kullback–Leibler (KL) divergence rate, if the state is translation-invariant and spatially ergodic. Our proof consists of two parts and is phrased in terms of a branch of the quantum information theory called the resource theory. First, we prove that any states, for which the min and max Rényi divergences collapse approximately to a single value, can be approximately reversibly converted into one another by thermal operations with the aid of a small source of quantum coherence. Second, we prove that these divergences collapse asymptotically to the KL divergence rate for any translation-invariant ergodic state. We show this via a generalization of the quantum Stein’s lemma for quantum hypothesis testing beyond independent and identically distributed situations. Our result implies that the KL divergence rate serves as a thermodynamic potential that provides a complete characterization of thermodynamic convertibility of ergodic states of quantum many-body systems in the thermodynamic limit, including out-of-equilibrium and fully quantum situations.

Physics↗
DOE OSTI2022

Femtosecond diffraction studies of the sodium chloride phase diagram under laser shock compression

The phase diagram of sodium chloride (NaCl) under laser shock compression has been studied at Linac Coherent Light Source (LCLS) at the x-ray free-electron laser facility. Both solid–solid (B1 → B2) and solid–liquid (B2 → liquid) transitions have been observed along the Hugoniot over nanosecond time scales. By combining structural measurements through in situ x-ray diffraction, pressure determination through velocimetry, and a thermal equation-of-state, the shock-compressed data are used to constrain the phase diagram of NaCl. Transformation into the B2 phase is found to occur at 28(2) GPa, and B2–liquid coexistence is observed between 54(4) and 66(6) GPa, with near full melt at 66(6) GPa. Late-time pressure release from an initial shocked B2-state results in a B2 → B1 back transformation. Our results show agreement with previous static compression data, suggesting that the time scale for melting is very rapid and that equilibrium states in NaCl are being accessed over nanosecond time scales. Finally, a multiphase equation-of-state description of NaCl incorporated into a one-dimensional hydrocode is used to interpret pressure and temperature evolution over these rapid time scales.

36 MATERIALS SCIENCE↗