Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical mechanics models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

pyTCR: A tropical cyclone rainfall model for python

pyTCR is a climatology software package developed in the Python programming language. It integrates the capabilities of several legacy physical models and increases computational efficiency to allow rapid estimation of tropical cyclone (TC) rainfall consistent with the large-scale environment. Specifically, pyTCR implements a horizontally distributed and vertically integrated model [Zhu et al., 2013] for simulating rainfall driven by TCs. Along storm tracks, rainfall is estimated by computing the cross-boundary-layer, upward water vapor transport caused by different mechanisms including frictional convergence, vortex stretching, large-scale baroclinic effect (i.e., wind shear), topographic forcing, and radiative cooling [Lu et al., 2018]. The package provides essential functionalities for modeling and interpreting spatio-temporal TC rainfall data. pyTCR requires a limited number of model input parameters, making it a convenient and useful tool for analyzing rainfall mechanisms driven by TCs. To sample rare (most intense) rainfall events that are often of great societal interest, pyTCR adapts and leverages outputs from a statistical-dynamical TC downscaling model [Lin et al., 2023] capable of rapidly generating a large number of synthetic TCs given a certain climate. As a result, pyTCR significantly reduces computational effort and improves the efficiency in capturing extreme TC rainfall events at the tail of the distributions from limited datasets. Furthermore, the TC downscaling model is forced entirely by large-scale environmental conditions from reanalysis data or coupled General Circulation Models (GCMs), simplifying the projection of TC-induced rainfall and wind speed under future climate using pyTCR. Finally, pyTCR can be coupled with hydrological and wind models to assess risks associated with independent and compound events (e.g., storm surges and freshwater flooding).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Multiphysics Multiscale Simulation Platform for Damage, Environmental Degradation, and Life Prediction of CMCs in Extreme Environments

This project successfully developed a multiphysics, multiscale computational framework to enhance the design and development of CMCs, with a focus on modeling highly nonlinear, time-dependent damage mechanisms and material degradation under extreme conditions, such as those experienced in turbine service environments. The project made significant advances in improving our understanding of progressive damage, oxidative degradation, and time-dependent inelastic deformation in CMCs, with particular attention to the role of uncertainties in predictions. Key outcomes include the integration of advanced material characterization, uncertainty quantification, and multiphysics constitutive models to predict the behavior of CMCs over their service life. A novel multiscale methodology was employed, which integrated microscale constituent behaviors with structural-scale responses, enabling the manufacturing defects in the microstructure that are prone to damage nucleation. Through the development of DL algorithms, the project advanced the prediction of damage initiation and crack propagation, taking into account the defect morphology and statistical variations across multiple scales. The framework was rigorously validated using thermomechanical experiments, which tested CMCs under various mechanical loadings at elevated temperatures, further enhancing the model's predictive capability. Overall, the research outcomes have provided a more accurate, reliable method for predicting CMC component life, significantly advancing material design, and improving component reliability in extreme environments. This work has strong implications for the optimization of turbine components and other high-performance applications where CMCs are used.

03 NATURAL GAS↗
DOE OSTI2023

The endpoint of partial deconfinement

We study the matrix quantum mechanics of two free hermitian N × N matrices subject to a singlet constraint in the microcanonical ensemble. This is the simplest example of a theory that at large N has a confinement/deconfinement transition. In the microcanonical ensemble, it also exhibits partial deconfinement with a Hagedorn density of states. We argue that the entropy of these configurations, based on a combinatorial counting of Young diagrams, are dominated by Young diagrams that have the VKLS shape. When the shape gets to the maximal depth allowed for a Young diagram of SU(N), namely N, we argue that the system stops exhibiting the Hagedorn behavior. The number of boxes (energy) at the transition is N 2 /4, independent of the charge of the state.

1/N Expansion↗
DOE OSTI2020

Learning to Predict Crystal Plasticity at the Nanoscale: Deep Residual Networks and Size Effects in Uniaxial Compression Discrete Dislocation Simulations

The density and configurational changes of crystal dislocations during plastic deformation influence the mechanical properties of materials. These influences have become clearest in nanoscale experiments, in terms of strength, hardness and work hardening size effects in small volumes. The mechanical characterization of a model crystal may be cast as an inverse problem of deducing the defect population characteristics (density, correlations) in small volumes from the mechanical behavior. In this work, we demonstrate how a deep residual network can be used to deduce the dislocation characteristics of a sample of interest using only its surface strain profiles at small deformations, and then statistically predict the mechanical response of size-affected samples at larger deformations. As a testbed of our approach, we utilize high-throughput discrete dislocation simulations for systems of widths that range from nano- to micro- meters. We show that the proposed deep learning model significantly outperforms a traditional machine learning model, as well as accurately produces statistical predictions of the size effects in samples of various widths. By visualizing the filters in convolutional layers and saliency maps, we find that the proposed model is able to learn the significant features of sample strain profiles.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Sterilization of N95 Respirators via Gamma Radiation: Comparison of Post-sterilization Efficacy

This study evaluated gamma irradiation for sterilization and reuse of two models of N95 respirators after gamma radiation sterilization as a method to increase availability of N95 respirators during a shortage. The Sandia National Laboratories Gamma Irradiation Facility was used to irradiate two different models of N95 filtering facepiece respirators at doses ranging from 0 kGy(tissue) to 50 kGy(tissue). The following tests were used to determine the efficacy of the respirator after irradiation sterilization: Ambient Aerosol Condensation Nuclei Counter Quantitative Fit Test, tensile test, strain cycling, oscillatory dynamic mechanical analysis, microscopic image analysis of fiber layers, and electrostatic field measurements. Both of the respirator models exhibited statistically significant changes after gamma irradiation as shown by the Quantitative Fit Test, electrostatic testing and the aerosol testing. The change in electrostatic capability of the filter reduced the efficiency of challenging particles near the 200 nm size by approximately 40-50%. Both tested respirators showed statistically significant changes associated with gamma sterilization. However, our results indicate that choices in materials and manufacturing methods to achieve N95 filtration lead to different magnitudes of damage when exposed to gamma radiation at sterilization relevant doses. This damage results in lower filtration performance. While our sample size (2 different types of respirators) was small, we did observe a change in electrostatic properties on a filter layer that coincided with the failure on the Quantitative Fit Test.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Sterilization of N95 Respirators via Gamma Radiation: Comparison of Post-sterilization Efficacy

This study evaluated gamma irradiation for sterilization and reuse of two models of N95 respirators after gamma radiation sterilization as a method to increase availability of N95 respirators during a shortage. The Sandia National Laboratories Gamma Irradiation Facility was used to irradiate two different models of N95 filtering facepiece respirators at doses ranging from 0 kGy(tissue) to 50 kGy(tissue). The following tests were used to determine the efficacy of the respirator after irradiation sterilization: Ambient Aerosol Condensation Nuclei Counter Quantitative Fit Test, tensile test, strain cycling, oscillatory dynamic mechanical analysis, microscopic image analysis of fiber layers, and electrostatic field measurements. Both of the respirator models exhibited statistically significant changes after gamma irradiation as shown by the Quantitative Fit Test, electrostatic testing and the aerosol testing. The change in electrostatic charge of the filter was correlated with a reduction in capturing particles near the 200 nm size by approximately 40-50%. Both tested respirators showed statistically significant changes associated with gamma sterilization. However, our results indicate that choices in materials and manufacturing methods to achieve N95 filtration lead to different magnitudes of damage when exposed to gamma radiation at sterilization relevant doses. This damage results in lower filtration performance. While our sample size (2 different types of respirators) was small, we did observe a change in electrostatic properties on a filter layer that coincided with the failure on the Quantitative Fit Test and reduction in aerosol filtering efficiency. Key Words: N95 respirators, respirators, airborne transmission, pandemic prevention, COVID-19, gamma sterilization

COVID-19↗
DOE OSTI2021

Quantifying High Temperature Corrosion

Alloys designed for high temperature service can form a variety of surface oxides or scales based on the alloy composition, service temperature, gas environment, surface deposits, etc. The most severe environments, including molten salts and liquid metals can result in significant metal loss, void formation and/or pitting. Standard and evolving practices are discussed for measuring these types of degradationincluding measuring oxide scale thickness, internal oxidation and/or metal loss, correlating them to mass change and reporting the results. Manual and automated methodologies to measure oxide thickness and metal loss are reviewed and compared in order to provide a summary of techniques and identify the best practices for quantifying a variety of materials and damage mechanisms and producing statistically meaningful results. Such techniques create datasets useful for improved understanding of corrosion degradation mechanisms, potentially better predictive models and enhanced data analytics.

Su, Yi Feng↗
DOE OSTI2020

Hadronic contributions to (g - 2) µ

The Muon g-2 Experiment at Fermilab, which recently started running, plans to reduce the uncertain- ties on the already very precisely measured anomalous magnetic moment of the muon by a factor of four. The goal of this effort is to probe the observed difference of more than three standard deviations between Standard-Model theory and experiment, one of the few persistent hints for physics beyond the Standard Model. The Fermilab experiment collected data from their first run last year with statistics comparable to BNL E821. They expect to release the measurement result in 2019. On the theoretical side, because the muon g - 2 arises from quantum- mechanical loop contributions in the Standard Model, it is sensitive to virtual effects of new particles, and places important constraints on Standard-Model extensions. To leverage the anticipated reduction in experimental errors, and determine unambiguously whether or not new-physics effects contribute to this quantity, the theoretical errors must be made more reliable and reduced to a commensurate precision. The Muon g-2 Theory Initiative was created to facilitate this development. The dominant sources of uncertainty in the Standard-Model prediction of the muon g -2 are from the hadronic contributions. The hadronic vacuum polarization (HVP) provides the leading correction followed by hadronic light-by-light (HLbL) scattering. There are a number of complementary theoretical efforts underway to better understand and quantify these contributions, including dispersive and data driven methods, lattice QCD, and effective field theories. Given the precision goals and the phenomenological importance, it is important to have more than one independent method for each of the two hadronic corrections, each with fully quantified uncertainties. Fostering the development of such methods is a prime goal of the initiative, as this will enable critical cross checks, and, upon combination, may yield gains in precision, to maximize the impact of E989. An important aspect of the Muon g-2 Theory Initiative’s activities are providing platforms that facilitate interactions between the different groups, as well as between the theoretical and experimental g - 2 communities. To this end, several workshops were organized in 2017 and 2018. The first meeting, held at Fermilab (June 3–6, 2017, St. Charles, IL, USA), served to kick-off the Initiative’s activities. Two meetings in early 2018 were focused respectively on the HVP and HLbL corrections. The HVP meeting was held at KEK (February 12–14, 2018, Tsukuba, Japan) and the HLbL meeting at the University of Connecticut (March 12–14, 2018, Storrs, CT, USA). The most recent workshop, which served as the second plenary meeting of the Initiative, was held at the University of Mainz (June 18–22, 2018, Mainz, Germany). An important outcome of these meetings are concrete plans for a first white paper, which is currently being written. We aim to post the white paper just prior to the release of the first E989 measurement, to present a clean theoretical prediction. The first white paper is focused on assessing and improving the reliability of the SM prediction. The purpose of the INT workshop in September 2019 is to start the next stage of the Initiative, focusing on the development of strategies to improve the theory uncertainties beyond the current level towards the E989 precision goal. We aim to accelerate theoretical developments on the hadronic contributions to the muon g - 2 so that the Standard-Model theory error can be brought to the needed precision, again in advance of the next release from E989. Hence the workshop will provide crucial theory support for a US experiment with broad impact. The Muon g-2 Theory Initiative relies on input from representatives of all the different communities that are engaged in this effort. It is therefore important that all these areas are properly represented. The funds from this grant will be used to enable more people, especially early-career scientists, to participate and make essential contributions to the workshop discussions. Since the workshop’s main goal is to kick-off the next stage of the theory initiative’s activities, support for this workshop from the DOE will help the theory community provide crucial theoretical support to the Fermilab Muon g-2 Experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

An analysis of fluff formation in metallic fuel via data analyzes from EBR-II experiments and BISON fuel code modeling

During the operation of EBR-II, it was found that a highly porous structure (over 40% area fraction) formed at the top of several fuel columns. Previous work has shown that this structure, designated fluff in this paper, contains a significant fraction of fuel elements (e.g., U and Pu) which could potentially impact neutronics. This work aims in analyzing the formation mechanism of this microstructure so its impact can be incorporated into future metallic fuel modeling codes and algorithms. This paper details a preliminary examination into the formation mechanisms of fluff by performing qualitative and statistical analysis of EBR-II experimental data. Additionally, the operating conditions that have the greatest impact on fluff formation were determined based on this data set. Also, BISON fuel code simulations were used to help postulate potential fluff formation mechanisms. From this analysis it was found that the largest contributors to fluff formation were fuel burnup and composition, with fluff formation exhibiting a roughly linear positive correlation with increasing burnup and a negative correlation with increasing Pu content. It was also found that higher pin operating temperatures decreased fluff formation but only for U-10Zr fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Differentially Private Map Matching (DPMM) v1.0

Human mobility trajectories provide valuable information for developing mobility applications, as they contain diverse and rich information about the users. User mobility data is valuable for various applications such as intelligent transportation systems (ITS), commercial business models, and disease-spread models. However, such spatio-temporal traces may pose a threat to user privacy. GPS trajectories in their raw form are not suitable for transportation studies, as they require matching locations with nearest road links — a process called map-matching. This software implements a differential privacy (DP)-based map-matching algorithm, called DPMM, that generates link-level location trajectories in a privacy-preserving manner to protect users' origin destinations (OD) and travel paths. OD privacy is achieved by injecting Planar Laplace noise to the user OD GPS points. Travel-path privacy is provided with randomized travel path construction using exponential DP mechanism. The injected noise level is selected adaptively, by considering the link density of the location and the functional category of the localized links. For path privacy, our mechanism samples waypoints and selects candidate paths between waypoints. DPMM provides privacy effectively with respect to link density instead of other trajectory samples in the database compared to other privacy mechanisms. Compared to the different baseline models our DP-based privacy model offers closer query responses to the raw data in terms of individual and aggregate trajectory-level statistics with an average at absolute deviation from the baseline for individual statistics on ϵ = 1.0. Beyond individual trajectory statistics, the DPMM outperforms the other benchmark DP-based mechanisms on different aggregate statistics with up to 8x improvement in utility.

Peisert, Sean [Lawrence Berkeley National Laborato↗
DOE OSTI2022

Atomic-scale mechanism of carbon nucleation from a deep crustal fluid by replica exchange reactive molecular dynamics simulation

Here we present a mechanistic model of carbon nucleation and growth from a fluid at elevated temperature (T) and pressure conditions, typical of those found in the shallow Earth’s lithosphere. Our model uses a replica exchange reactive molecular dynamics framework in which molecular configurations are swapped between adjacent T replica at regular intervals according to underlying statistical mechanics. This framework allows predicting complex molecular structures and thermodynamics while remaining computationally efficient. Here we simulate the reactivity of an unstable mixture of CO 2 and CH 4 at 1000 K and 1 GPa. We find that the path to thermodynamic equilibrium is initially entropy-driven, producing a diversity of short-lived species, including various alcohols with intermediate carbon oxidation states. Cyclic and polycyclic radicals that are sometimes resonance-stabilized form next and set the stage for carbon nucleation. The carbon exsolution process releases abundant water, is exothermic and starts with the nucleation of a large aggregate of hydrogenated graphene flakes from covalently bonded polycyclic units. The carbon backbone of this nucleus subsequently grows into a hydrogen-depleted fullerene-like structure, before evolving toward a partially bilayered graphene layer. Overall, our results show that the mechanism of graphitic C formation is certainly not bimolecular, and that it may involve a combination of key condensation and radical chain reactions. This will help understand the isotopic, and reactive characteristics of carbon-bearing fluids during their upward transit through the Earth’s mantle and crust. Moreover, the mechanistic insights outlined here present intriguing similarities with the process of soot and interstellar dust formation, which suggests that the widespread distribution of abiotic polyaromatic and graphitic material on Earth and beyond may reflect the prevalence of a fundamental chemical pathway.

58 GEOSCIENCES↗
DOE OSTI2021

First Principles Simulations of Electrified Interfaces in Electrochemistry

This chapter discusses some of the recent advances made in the first principles modeling of electrochemical catalysts. It also discusses the key development, namely the ability to explicitly treat the effects of surface electrification due to electrochemical processes and applied voltages in a computationally efficient manner. The chapter introduces the thermodynamics and statistical mechanics of electrified metal-solution interfaces. It then discusses the structure of the electrode-electrolyte interface and the effects of applied voltages, followed by an overview of the first-principles model and a motivating example. The chapter also provides a brief summary of classical thermodynamics and describes the basics of thermodynamic detour. It considers the thermodynamics of macroscopic systems and introduces several useful fundamental relations. Here, the chapter also considers a system to exist in a certain macrostate, which is a particular thermodynamic state specified by a set of fixed properties such as constant particle number, constant volume, and constant temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Bottom-Up Simulation, Reconstruction, and Quantification of Macromolecule Sequences from Experimental Polymerizations

Motivated by the canonical sequence–structure–function paradigm, tools to characterize chemical patterning in natural biomacromolecules, from proteins to nucleic acids, have grown exponentially in recent years. However, analogous strategies for synthetic macromolecules remain in nascent stages, complicated by sequence polydispersity and analytical limitations. To address this, we have developed a comprehensive and open-source Python package, PRISM (polymer rate insights and sequence modeling), an end-to-end workflow that provides a path from experimental kinetics measurements to quantitative and qualitative metrics for describing chemical patterning in stochastic polymers. First, a numerical integration strategy was constructed to simulate and fit experimental data from reversible addition–fragmentation chain transfer (RAFT) polymerization kinetics, enabling the facile estimation of relevant reactivity ratios. These ratios were then used in a mechanism-specific stochastic kinetic simulation strategy to simulate sequence ensembles corresponding to model systems spanning experimental copolymers, classes of statistical polymers (e.g., alternating, block, and gradient), and multiblock copolymers. Lastly, inspired by sequence homology metrics from bioinformatics, we introduce visualization strategies and quantitative metrics to facilitate comparisons of different sequence ensembles. As the sequence–structure–function paradigm becomes increasingly central in de novo design of synthetic macromolecules, this toolkit provides a first step toward accurate and representative sequence description and featurization.

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

Gauge-symmetrization method for energy-momentum tensors in high-order electromagnetic field theories

For electromagnetic field theories, canonical energy-momentum conservation laws can be derived from the underpinning spacetime translation symmetry according to the Noether procedure. However, the canonical energy-momentum tensors (EMTs) are neither symmetric nor gauge-symmetric (gauge invariant). The Belinfante-Rosenfeld (BR) method is a well-known procedure to symmetrize the EMTs, which also renders them gauge symmetric for first-order field theories. High-order electromagnetic field theories appear in the study of gyrokinetic systems for magnetized plasmas and the Podolsky system for the radiation reaction of classical charged particles. For these high-order field theories, gauge-symmetric EMTs are not necessarily symmetric and vice versa. In the present study, we develop a new gauge-symmetrization method for EMTs in high-order electromagnetic field theories. The Noether procedure is carried out using the Faraday tensor $F_{μν}$, instead of the 4-potential $A_{μ}$, to derive a canonical EMT $T^{μν}_{N}$. We show that the gauge-dependent part of $T^{μν}_{N}$ can be removed using the displacement-potential tensor $F^{σμν}$ ≡ $D^{σμ}A^{ν}/4π$, where $D^{σμ}$ is the antisymmetric electric displacement tensor. This method gauge-symmetrizes the EMT without necessarily making it symmetric, which is adequate for applications not involving general relativity. For first-order electromagnetic field theories, such as the standard Maxwell system, $F^{σμν}$ reduces to the familiar BR superpotential $S^{σμν}$, and the method developed can be used as a simpler procedure to calculate $S^{σμν}$ without employing the angular momentum tensor in 4D spacetime. When the electromagnetic system is coupled to classical charged particles, the gauge-symmetrization method for EMTs is shown to be effective as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Physics-Informed Gaussian Process Regression for States Estimation and Forecasting in Power Grids

Real-time state estimation and forecasting are critical for the efficient operation of power grids. In this paper, a physics-informed Gaussian process regression (PhI-GPR) method is presented and used for forecasting and estimating the phase angle, angular speed, and wind mechanical power of a three-generator power grid system using sparse measurements. In standard data-driven Gaussian process regression (GPR), parameterized models for the prior statistics are fit by maximizing the marginal likelihood of observed data. In the PhI-GPR method, we propose to compute the prior statistics offline by solving stochastic differential equations (SDEs) governing the power grid dynamics. The short-term forecast of a power grid system dominated by wind generation is complicated by the stochastic nature of the wind and the resulting uncertainty in wind mechanical power. Here, we assume that the power grid dynamics are governed by swing equations, with the wind mechanical power fluctuating randomly in time. We solve these equations for the mean and covariances of the power grid states using the Monte Carlo simulation method. We demonstrate that the proposed PhI-GPR method can accurately forecast and estimate observed and unobserved states. For the considered problem, PhI-GPR has computational advantages over the ensemble Kalman filter (EnKF) method: In PhI-GPR, ensembles are computed offline and independently of the data acquisition process, whereas for EnFK, ensembles are computed online with data acquisition, rendering real-time forecast more challenging. We also demonstrate that the PhI-GPR forecast is more accurate than the EnKF forecast when the random mechanical wind power is non-Markovian. In contrast, the two methods produce similar forecasts for the Markovian mechanical wind power. For observed states, we show that PhI-GPR provides a forecast comparable to the standard data-driven GPR; both forecasts are significantly more accurate than the autoregressive integrated moving average (ARIMA) forecast. We also show that the ARIMA forecast is more sensitive to observation frequency and measurement errors than the PhI-GPR forecast.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Harness the power of atomistic modeling and deep learning in biofuel separation

Biofuels offer a remarkable, sustainable energy source for a future of clean energy. The development of efficient biofuel separation plays a crucial role in achieving cost-effective utilization of biofuel. In this chapter, we provide an overview of the recent advancements in atomistic-level modeling and deep learning in the rational design of novel, efficient biofuel separation. The fundamental principles of quantum and statistical mechanics are covered in appropriate detail to highlight their underlying differences in theory. The methodologies of several molecular representations and deep learning algorithms applicable to biofuel separation are briefly demonstrated as well. The applications, successes, and risks of employing density functional theory, ab initio molecular dynamics, classical molecular dynamics, and deep learning are provided to showcase their recent accomplishments in biofuel separation as well as potential improvements in both methodology and application. Lastly, a vision for the future growth of these methods is illustrated.

deep learning, artificial intelligence, biofuels, ↗
DOE OSTI2024

Physics-informed machine learning of the correlation functions in bulk fluids

The Ornstein–Zernike (OZ) equation is the fundamental equation for pair correlation function computations in the modern integral equation theory for liquids. In this work, machine learning models, notably physics-informed neural networks and physics-informed neural operator networks, are explored to solve the OZ equation. The physics-informed machine learning models demonstrate great accuracy and high efficiency in solving the forward and inverse OZ problems of various bulk fluids. The results highlight the significant potential of physics-informed machine learning for applications in thermodynamic state theory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Trajectory Ensemble Methods Provide Single-Molecule Statistics for Quantum Dynamical Systems

The emergence of experiments capable of probing quantum dynamics at the single-molecule level requires the development of new theoretical tools capable of simulating and analyzing these dynamics beyond an ensemble-averaged description. In this article, we present an efficient method for sampling and simulating the dynamics of the individual quantum systems that make up an ensemble and apply it to study the nonequilibrium dynamics of the ubiquitous spin-boson model. We generate an ensemble of single-system trajectories, and we analyze this trajectory ensemble using tools from classical statistical mechanics. Our results demonstrate that the dynamics of quantum coherence is highly heterogeneous at the single-system level due to variations in the initial bath configuration, which significantly affects the transient exchange of coherence between the system and its bath. Here, we observe that single systems tend to retain coherence over time scales longer than that of the ensemble. We also compute a novel thermodynamic entanglement entropy that quantifies a thermodynamic driving force favoring system–bath entanglement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗