Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “expansion techniques”

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 217 records · Page 12

Thermionic microwave gun for terahertz and synchrotron light sources

Conventional thermionic microwave and radiofrequency (RF) guns can offer high average beam current, which is important for synchrotron light and terahertz (THz) radiation source facilities, as well as for industrial applications. For example, the Advanced Photon Source at Argonne National Laboratory is a national synchrotron radiation light source research facility that utilizes thermionic RF guns. Yet, these existing thermionic guns are bulky, difficult to handle and install, easily detuned, very sensitive to thermal expansion, and are due for a major upgrade and replacement. Here, we present the design of a new, more stable and reliable gun with optimized electromagnetic performance, improved thermal engineering and a more robust cathode mounting technique, which is a critical step to improve the performance of existing and future light sources, industrial accelerators, and electron beam-driven terahertz sources. We will also present a fabricated gun prototype and show results of high-power and beam tests.

47 OTHER INSTRUMENTATION↗

Measuring the density, viscosity, and surface tension of molten titanates using electrostatic levitation in microgravity

Rare earth and barium titanates are useful as ferroelectric, dielectric, and optical materials. Measurements of their thermophysical properties in the liquid state can help guide melt processing technologies for their manufacture and advance understanding of fragile liquids' behavior and glass formation. Here, we report the density, thermal expansion, viscosity, and surface tension of molten BaTi 2 O 5 , BaTi 4 O 9 , and 83TiO 2 -17RE 2 O 3 (RE = La or Nd). Measurements were made using electrostatic levitation and droplet oscillation techniques in microgravity, which provide access to quiescent liquid droplets and deep supercooling of 510–815 K below the equilibrium melting points. Densities were measured over 900–2400 K. Viscosities were similar for all four compositions, increasing from ~10 mPa s near 2100 K to ~30 mPa s near 1750 K. Surface tensions were 450–490 dyn cm -1 for the rare earth titanates and 383–395 dyn cm -1 for the barium titanates; surface tensions of all compositions had small or negligible temperature dependence over 1700–2200 K. In conclusion, for solids recovered after melt quenching, x-ray microtomography revealed the fracture mechanics in crystalline products and minimal internal porosity in glass products, likely arising from entrapped gas bubbles. Internal microstructures were generally similar for products processed either in microgravity or in a terrestrial aerodynamic levitator.

42 ENGINEERING↗

Novel One-Step Production of Carbon-Coated Sn Nanoparticles for High-Capacity Anodes in Lithium-Ion Batteries

Lithium-ion batteries offer the highest energy density of any currently available portable energy storage technology. By using different anode materials, these batteries could have an even greater energy density. One material, tin, has a theoretical lithium capacity (994 mAh/g) over three-times higher than commercial carbon anode materials. Unfortunately, to achieve this high capacity, bulk tin undergoes a large volume expansion, and the material pulverizes during cycling, giving a rapid capacity fade. To mitigate this issue, tin must be scaled down to the nano-level to take advantage of unique micromechanics at the nanoscale. Synthesis techniques for Sn nanoparticle anodes are costly and overly complicated for commercial production. A novel one-step process for producing carbon-coated Sn nanoparticles via spark plasma erosion (SPE) shows great promise as a simple, inexpensive production method. The SPE method, characterization of the resulting particles, and their high-capacity reversible electrochemical performance as anodes are described. With only a 10% addition of these novel SPE carbon-coated Sn particles, one anode composition demonstrated a reversible capacity of ~460 mAh/g, achieving the theoretical capacity of that particular electrode formulation. These SPE carbon-coated Sn nanoparticles are drop-in ready for present commercial lithium-ion anode processing and would provide a ~10% increase in the total capacity of current commercial lithium-ion cells.

25 ENERGY STORAGE↗

Rare Earth Element Detection and Quantification in Coal and Rock Mineral Matrices

As global demand for rare earth elements (REEs) increases, maintaining the production and supply chain is critical. Technologies capable of being used in the field and in situ in the subsurface for rapid REE detection and quantification facilitates the efficient mining of known resources and exploration of new and unconventional resources. Laser-induced breakdown spectroscopy (LIBS) is a promising technique for rapid elemental analysis both in the laboratory and in the field. Multiple articles have been published evaluating LIBS for detection and quantification of REEs; however, REEs in their natural deposits have not been adequately studied. In this work, detection and quantification of two REEs, La and Nd, have been studied in both synthetic and natural mineral matrices at concentrations relevant to REE extraction. Measurements were performed on REE-containing rock and coal samples (natural and synthetic) utilizing different LIBS instruments and techniques, specifically a commercial benchtop instrument, a custom benchtop instrument (single- and double-pulse modes), and a custom LIBS probe currently being developed for in situ, subsurface, borehole wall detection and quantification of REEs. Plasma expansion, emission intensity, detection limits, and double-pulse signal enhancement were studied. The limits of detection (LOD) were found to be 10/14 ppm for La and 15/25 ppm for Nd in simulated coal/rock matrices in single-pulse mode. Signal enhancement of 3.5 to 6-fold was obtained with double-pulse mode as compared to single-pulse operation.

La detection↗

A DATA EFFICIENT SPARSE MODELING FRAMEWORK FOR POWER ESTIMATION IN WATER TREATMENT SENSING OPERATIONS

With increasing freshwater scarcity, advanced process design mechanisms such as Closed-Circuit Reverse Osmosis (CCRO) and Digital/Physical Twin systems are gaining traction in water treatment and reuse operations. While digital and physical twin models enable improved system insight and control, their development is often expensive and computationally intensive, requiring large volumes of synthetic or experimental data to characterize underlying process dynamics. This work introduces a sparse surrogate modeling framework to estimate power consumption from measured flow and pressure variables, along with their nonlinear polynomial and interaction expansions. To ensure model reliability and reduce overfitting, a two-stage pipeline is proposed. First, a dynamic data filtering algorithm is employed to remove uninformative observations and transient operational states. Second, a sparse penalized regression technique is applied to select a minimal set of parsimonious features. The proposed model achieves high sparsity, retaining only 7 out of 34 candidate features (≈79.41% sparsity) while delivering a root mean square error (RMSE) of 0.072 on the test dataset.

Mukherjee, Subrata [ORNL] (ORCID:0000000309930338)↗

Gold–Thiolate Nanocluster Dynamics and Intercluster Reactions Enabled by a Machine Learned Interatomic Potential

Monolayer protected metal clusters comprise a rich class of molecular systems and are promising candidate materials for a variety of applications. While a growing number of protected nanoclusters have been synthesized and characterized in crystalline forms, their dynamical behavior in solution, including prenucleation cluster formation, is not well understood due to limitations both in characterization and first-principles modeling techniques. Recent advancements in machine-learned interatomic potentials are rapidly enabling the study of complex interactions such as dynamical behavior and reactivity on the nanoscale. Here, we develop an Au-S-C-H atomic cluster expansion (ACE) interatomic potential for efficient and accurate molecular dynamics simulations of thiolate-protected gold nanoclusters (Au n (SCH 3 ) m ). Trained on more than 30,000 density functional theory calculations of gold nanoclusters, the interatomic potential exhibits ab initio level accuracy in energies and forces and replicates nanocluster dynamics including thermal vibration and chiral inversion. Long dynamics simulations (up to 0.1 μs time scale) reveal a mechanism explaining the thermal instability of neutral Au 25 (SR) 18 clusters. Specifically, we observe multiple stages of isomerization of the Au 25 (SR) 18 cluster, including a chiral isomer. Additionally, we simulate coalescence of two Au 25 (SR) 18 clusters and observe series of clusters where the formation mechanisms are critically mediated by ligand exchange in the form of [Au-S] n rings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effective bias expansion for 21-cm cosmology in redshift space

A near-future detection of the 21-cm signal from the epoch of reionization will provide unique opportunities to probe the underlying cosmology, provided that such cosmological information can be extracted with precision. To this end, we further develop effective field theory (EFT) inspired techniques for the 21-cm brightness temperature field during the epoch of reionization, incorporating renormalized bias and a treatment of redshift space distortions. Notably, we confirm that in redshift space, measures of the 21-cm brightness, e.g., the power spectrum, should have irreducible contributions that lack a bias coefficient and therefore contain direct, astrophysics-free information about the cosmological density field; in this work, we study this effect beyond linear order. To validate our theoretical treatment, we fit the predicted EFT Fourier-space shapes to the thesan suite of hydrodynamical simulations of reionization at the field level, where the considerable number of modes prevents overfitting. We find agreement at the level of a few percent between the 21-cm power spectrum from the EFT fits and simulations over the wave number range k ≲ 0.8 h/Mpc and neutral fraction x HI ≳ 0.4, which is imminently measurable by the Hydrogen Epoch of Reionization Array and future experiments. Furthermore, the ability of the EFT to describe the 21-cm signal extends to simulations that have different astrophysical prescriptions for reionization as well as simulations with interacting dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗

Mechanically-, Thermally, and Chemically-Robust High-Temperature Ceramic Composites

In an effort to lower the levelized cost of electricity generated by Concentrated Solar Power (CSP) plants, the heat-to-electricity conversion efficiency may be increased by operating such plants with higher turbine inlet temperatures using closed Brayton cycles with high-pressure, supercritical carbon dioxide (sCO 2 ) as the working fluid. For example, if the turbine inlet temperature is raised to ≥750°C, in combination with sCO2 closed Brayton cycles, the relative heat-to-electricity efficiency may be raised by more than 20% (compared to subcritical steam Rankine cycles operating with turbine inlet temperatures ≤550°C). The associated reduction in the cost of dispatchable electricity from such CSP plants (utilizing high-temperature thermal energy storage) would be an important breakthrough on the path towards direct competition with fossil-fuel-based plants. A key barrier to achieve such high-temperature efficient operation has been the limited thermomechanical performance of metal alloys used in compact, primary printed circuit- type heat exchangers (PCHEXs) for heat transfer to high-pressure sCO 2 . The maximum allowed stresses for the use of conventional stainless steels and nickel-based superalloys at high sCO 2 pressures (≥20 MPa) decline rapidly at temperatures ≥550°C. This project has been focused on demonstrating the desired high-temperature properties and attractive manufacturing characteristics of mechanically-, thermally-, and chemically-robust ceramic/metal composites (cermets), in order to allow for the use of such composites in primary PCHEXs for heat transfer to sCO 2 at ≥750°C and ≥20 MPa in CSP plants. Prior work with high-temperature, co-continuous carbide/refractory metal cermets, such as ZrC/W composites, has already demonstrated: i) the attractive combination of high-temperature properties exhibited by such materials (e.g., high values of stiffness, failure strength, and thermal conductivity at 800°C), and ii) cost-effective processes for manufacturing PCHEX plates comprised of such cermets with tailorable channel patterns. However, ZrC/W cermets were not found to be inherently corrosion resistant in sCO 2 at 750°C and 20 MPa. While the application of a Cu coating to ZrC/W surfaces, along with a modest addition (50 ppm) of CO to CO 2 , rendered such composites resistant to corrosion for 1000 h in such sCO 2 -based fluids at 750°C, such additional steps would increase the design complexity of such cermet-based components. An alternative approach (the focus of this project) is to develop mechanically-robust, cost- effective cermets that are inherently resistant to high-temperature oxidation in air and CO 2 . The overall goals of this project have been to demonstrate: i) the oxidation resistance of at least one cermet at 750°C in air and in CO 2 (i.e., with a projected annual corrosion ≤30 μm), ii) the mechanically robust nature of at least one cermet at 750°C (i.e., with an average failure strength ≥200 MPa), and iii) that the selected cermet can be manufactured via low-cost forming techniques.

14 SOLAR ENERGY↗

Second-order renormalized Hamiltonian of Yukawa theory

Using the renormalization group procedure for effective particles we calculate the effective Hamiltonians in the theory of a fermion field coupled to a scalar field via the Yukawa interaction. The theory is renormalized by the addition of counterterms. Necessary counterterms are determined by computing matrix elements of the effective Hamiltonian. All calculations are performed up to the second order in the expansion in powers of the coupling constant. Renormalized effective Hamiltonians are well-defined symmetric forms acting in the Fock space as opposed to the renormalized bare Hamiltonian, which is not well defined without regularization. We introduce computational techniques that should streamline higher-order calculations and may be of independent interest.

Ab initio calculations↗

Dipolar Spin Ice Regime Proximate to an All-In-All-Out Néel Ground State in the Dipolar-Octupolar Pyrochlore Ce 2 Sn 2 O 7

The dipolar-octupolar (DO) pyrochlores, R 2 M 2 O 7 ( R = Ce , Sm , Nd ), are key players in the search for realizable novel quantum spin liquid (QSL) states as a large parameter space within the DO pyrochlore phase diagram is theorized to host QSL states of both dipolar and octupolar nature. New single crystals and powders of Ce 2 Sn 2 O 7 , synthesized by hydrothermal techniques, present an opportunity for a new characterization of the exchange parameters in Ce 2 Sn 2 O 7 using the near-neighbor X Y Z model Hamiltonian associated with DO pyrochlores. Utilizing quantum numerical linked cluster expansion fits to heat capacity and magnetic susceptibility measurements, and classical Monte Carlo calculations to the diffuse neutron diffraction of the new hydrothermally grown Ce 2 Sn 2 O 7 samples, we place Ce 2 Sn 2 O 7 ’s ground state within the ordered dipolar all-in-all-out (AIAO) Néel phase, with quantum Monte Carlo calculations showing a transition to long-range order at temperatures below those accessed experimentally. Indeed, our new neutron diffraction measurements on the hydrothermally grown Ce 2 Sn 2 O 7 powders show a broad signal at low scattering wave vectors, reminiscent of a spin ice, in striking contrast from previous powder neutron diffraction on samples grown from solid-state synthesis, which found diffuse scattering at high scattering wave vectors associated with magnetic and suggested an octupolar quantum spin ice state. We conclude that new hydrothermally grown Ce 2 Sn 2 O 7 samples host a finite-temperature proximate dipolar spin ice phase, above the expected transition to AIAO Néel order. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hierarchical Gaussian Random Field Sampling for Multilevel Markov Chain Monte Carlo: Coupling Stochastic Partial Differential Equation and the Karhunen–Loève Decomposition

This work introduces structure preserving hierarchical decompositions for sampling Gaussian random fields (GRFs) within the context of multilevel Bayesian inference in high-dimensional space. Existing scalable hierarchical sampling methods, such as those based on stochastic partial differential equations (SPDEs), often reduce the dimensionality of the sample space at the cost of accuracy of inference. Other approaches, such that those based on Karhunen-Loève (KL) expansions, offer sample space dimensionality reduction but sacrifice GRF representation accuracy and ergodicity of the Markov chain Monte Carlo (MCMC) sampler and are computationally expensive for high-dimensional problems. The proposed method integrates the dimensionality reduction capabilities of KL expansions with the scalability of SPDE-based sampling, thereby providing a robust, unified framework for high-dimensional uncertainty quantification (UQ) that is scalable and accurate, preserves ergodicity, and offers dimensionality reduction of the sample space. The hierarchy in our multilevel algorithm is derived from the geometric multigrid hierarchy. By constructing a hierarchical decomposition that maintains the covariance structure across the levels in the hierarchy, the approach enables efficient coarse-to-fine sampling while ensuring that all samples are drawn from the desired distribution. The effectiveness of the proposed method is demonstrated on a benchmark subsurface flow problem, demonstrating its effectiveness in improving computational efficiency and statistical accuracy. Furthermore, our proposed technique is more efficient and accurate and displays better convergence properties than existing methods for high-dimensional Bayesian inference problems.

Gaussian random fields↗

Powder X-ray diffraction of altrenogest, C 21 H 26 O 2

The crystal structure of altrenogest has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional techniques. Altrenogest crystallizes in space group P2 1 2 1 2 1 (#19) with a = 7.286 916(16), b = 10.580 333(19), c = 22.266 08(7) Å, V = 1716.671(6) Å 3 , and Z = 4 at 295 K. Thermal expansion between 113 and 295 K is anisotropic. An O–H···O hydrogen bond links the molecules into chains along the c-axis. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗

A new method for solving the linearized 1D Vlasov–Poisson system yielding a new class of solutions

We describe a new method for solving the linearized 1D Vlasov–Poisson system by using properties of Cauchy-type integrals. Our method remedies critical flaws of the two standard methods, reveals a previously unrecognized Gaussian-in-time-like decay, and can also account for an externally applied electric field. The Landau approximation involves deforming the Bromwich contour around the poles closest to the real axis due to the analytically continued dielectric function, finding the long-time behavior for a stable system: Landau damping. Jackson's generalization encircles all poles while sending the contour to infinity, assuming its contribution vanishes, which is not true in general. This gives incorrect solutions for physically reasonable configurations and can exhibit pathological behavior, of which we show examples. The van Kampen method expresses the solution for a stable equilibrium as a continuous superposition of waves, resulting in an opaque integral. Case's generalization includes unstable systems and predicts a decaying discrete mode for each growing discrete mode, an apparent contradiction to both the Jackson solution and ours. We show, without imposing additional constraints, that the decaying modes are never present in the time evolution due to an exact cancellation with part of the continuum. Our solution is free of integral expressions, is obtained using algebra and Laurent series expansions, does not rely on analytic continuations, and results in a correct asymptotically convergent form in the case of infinite sums. The analysis used can be readily applied in higher-dimensional, electromagnetic systems and also provides a new technique for evaluating certain inverse Laplace transforms.

Physics↗

Considerations for in situ, real time measurement of plasma-material interactions using Digital Holographic imaging

Digital Holographic (DH) imaging is a laser-based measurement technique, which can be used to monitor a material surface as it is being exposed to fusion-relevant plasmas. Both single-laser and dual-laser DH measurements have been demonstrated ex situ at Oak Ridge National Laboratory (ORNL). A DH diagnostic system is being assessed for deployment at ORNL to make in situ, real-time measurements of plasma erosion/redeposition in the Prototype Material Plasma Exposure eXperiment (Proto-MPEX). In situ, real time measurements pose unique challenges to diagnostic systems, which are not encountered by ex situ analysis techniques. The absolute surface height (tracked at the nm to μm position) during plasma exposure can be modified by: 1) surface movement due to vibration or thrust exerted by the plasma device, 2) surface growth due to the thermal expansion of the material under steady MW/m2 or transient GW/m2 plasma heat loads, and 3) surface modification due to plasma erosion/redeposition of the material substrate. To assess these effects, 1) the vibration spectrum of the diagnostic and the pulsed plasma device have been measured, 2) the thermal growth of the surface has been measured for an applied heat flux, and 3) the surface modification has been measured ex situ post plasma exposure, and ex situ post laser ablation as a proxy for the plasma. This paper will provide an overview of the results that have been achieved in the development of a DH diagnostic for in situ real-time measurements of plasma exposed surfaces in the Proto-MPEX device.

Biewer, Theodore↗

Galaxy bispectrum in the spherical Fourier-Bessel basis

The bispectrum, the three-point correlation in Fourier space, is a crucial statistic for studying many effects targeted by the next-generation galaxy surveys, such as primordial non-Gaussianity (PNG) and general relativistic (GR) effects on large scales. In this work we develop a formalism for the bispectrum in the spherical Fourier-Bessel (SFB) basis—a natural basis for computing correlation functions on the curved sky, as it diagonalizes the Laplacian operator in spherical coordinates. Working in the SFB basis allows for line-of-sight effects such as redshift space distortions and GR to be accounted for exactly, i.e., without having to resort to perturbative expansions to go beyond the plane-parallel approximation. Only analytic results for the SFB bispectrum exist in the literature given the intensive computations needed. We numerically calculate the SFB bispectrum for the first time, enabled by a few techniques: We implement a template decomposition of the redshift-space kernel Z 2 into Legendre polynomials, and separately treat the PNG and velocity-divergence terms. We derive an identity to integrate a product of three spherical harmonics connected by a Dirac delta function as a simple sum and use it to investigate the limit of a homogeneous and isotropic Universe. Furthermore, we present a formalism for convolving the signal with separable window functions and use a toy spherically symmetric window to demonstrate the computation and give insights into the properties of the observed bispectrum signal. While our implementation remains computationally challenging, it is a step toward a feasible full extraction of information on large scales via a SFB bispectrum analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Advanced Models and Fast Simulation Algorithms for MARS Plants

With the successful commercialization of modular multilevel converters (MMC) and rapid increase in photovoltaic (PV) penetration and energy storage systems (ESS), one of the research areas being explored dynamically is integration of PV systems and ESS systems to MMCs. On this front, an integrated development of PV plants, ESS, and alternating- and direct-current (ac/dc) systems through Multi-Port Autonomous Reconfigurable Solar (MARS) plant is proposed and being researched upon. However, simulation of MARS systems on traditional existing software takes very long time (order of days) due to the presence of a large number of states and non-linear devices. Therefore, development of fast high fidelity models is of paramount importance to test the hardware design, control design and planning of power system expansions. In this context, two research objectives are addressed in this paper: firstly, an ultra fast single- or multi- CPU simulation algorithm to simulate the MARS system based on state-space models, hybrid discretization algorithm with a relaxation technique that reduces the imposed computational burden and a multi-rate method is presented. Secondly, a reduced order model of MARS is also presented. Both the developed techniques are validated with respect to reference PSCAD/EMTDC model.

Marthi, Phani Ratna Vanamali↗

Application of electron backscatter diffraction techniques to quantify effects of aging on sub-grain and spatial heterogeneity in NMC cathodes

Identification and evaluation of structural heterogeneity in spent cathode materials is crucial to developing appropriate remediation strategies for novel recycling processes. Native heterogeneities may be exacerbated during the cell's operational lifetime, as sub-particle-scale variations induce anisotropic expansion and contraction upon cycling. Structural transformations resulting from repeated cycling and calendar aging predominantly occur at the secondary particle surface and at the grain boundaries (GBs) between primary particles. However, the diffusion and stress build up around and across GBs are poorly understood. In this study, electron backscatter diffraction (EBSD) is employed to track sub-grain lattice structure across a statistically relevant number of Li(Ni0.33Mn0.33Co0.33)O2 (NMC-111) particles. Specifically, differences in lattice misorientation – measured as the deviation from the grain's average orientation – are tracked as a function of position within the electrode (near current collector, middle, and near separator) and as a function of electrochemical cycling. Further, a novel method of structural analysis is developed, offering insight into sub-grain diffusion behavior by comparing lattice misorientation near the grain boundary versus in the grain bulk. The present results suggest that electrode-scale spatial heterogeneity in lattice structure is induced by initial manufacturing conditions, and that radial gradients in lattice misorientation evolve at the primary particle scale with repeated electrochemical cycling.

25 ENERGY STORAGE↗

A Practical Framework for Simulating Time-Resolved Spectroscopy Based on a Real-Time Dyson Expansion

Time-resolved spectroscopy is a powerful tool for probing electron dynamics in molecules and solids, revealing transient phenomena on subfemtosecond time scales. The interpretation of experimental results is often enhanced by parallel numerical studies, which can provide insight and validation for experimental hypotheses. However, developing a theoretical framework for simulating time-resolved spectra remains a significant challenge. The most suitable approach involves the many-body nonequilibrium Green's function formalism, which accounts for crucial dynamical many-body correlations during time evolution. While these dynamical correlations are essential for observing emergent behavior in time-resolved spectra, they also render the formalism prohibitively expensive for large-scale simulations. Substantial effort has been devoted to reducing this computational cost─through approximations and numerical techniques─while preserving the key dynamical correlations. The ultimate goal is to enable first-principles simulations of time-dependent systems ranging from small molecules to large, periodic, multidimensional solids. Here, in this perspective, we outline key challenges in developing practical simulations for time-resolved spectroscopy, with a particular focus on Green's function methodologies. We highlight a recent advancement toward a scalable framework: the real-time Dyson expansion (RT-DE) [Phys. Rev. Lett. 2024, 133, 226902]. We introduce the theoretical foundation of RT-DE and discuss strategies for improving scalability, which have already enabled simulations of system sizes beyond the reach of previous fully dynamical approaches. We conclude with an outlook on future directions for extending RT-DE to first-principles studies of dynamically correlated, nonequilibrium systems.

Reeves, Cian C. [Univ. of California, Santa Barbar↗