Search NASA⌕ Search

SEARCH · Search NASA

Results for “STATISTICAL MECHANICS”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Direct estimation of the density of states for fermionic systems

Simulating time evolution is one of the most natural applications of quantum computers and is thus one of the most promising prospects for achieving practical quantum advantage. Here, we develop quantum algorithms to extract thermodynamic properties by estimating the density of states (DOS), which is a central object in quantum statistical mechanics. We introduce several key innovations that significantly improve the practicality and extend the generality of previous techniques. First, our approach allows one to estimate the DOS only for a specific subspace of the full Hilbert space. This is crucial for fermionic systems, since both canonical and grand canonical ensemble thermal equilibrium properties depend on subspaces of fixed number. Second, in our approach, by time evolving very simple, random initial states, such as randomly chosen computational basis states, we can exactly recover the DOS on average. Third, due to circuit-depth limitations, we only reconstruct the DOS up to a convolution with a Gaussian window—thus all imperfections that shift the energy levels by less than the width of the convolution window will not significantly affect the estimated DOS. For these reasons, we find the approach is a promising candidate for early quantum advantage as even short-time, noisy dynamics can yield a semiquantitative reconstruction of the DOS (convolution with a broad Gaussian window), while early fault-tolerant devices will likely enable higher-resolution DOS reconstruction through longer time evolutions. We demonstrate the practicality of our approach in representative Fermi-Hubbard and spin models and indeed find that our approach is highly robust against algorithmic errors in the time evolution and against gate noise. We further demonstrate that our approach is compatible with noisy intermediate-scale quantum (NISQ) computing NISQ-friendly variational techniques, introducing and leveraging a technique for variational time evolution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Effect of Entropic Constraints on the Thermodynamics of Molecular Adsorption in Nano‐Porous Materials

Abstract Gas separation is a critical industrial process that consumes a significant amount of energy due to the widely used techniques that are currently employed. Adsorptive materials—such as metal–organic frameworks (MOFs)—show promise as an energy‐efficient alternative. Of particular current interest are novel, temperature‐dependent separation processes in MOFs, such as the recently reported separation of ternary isomeric hydrocarbon mixtures within one and the same material. However, the mechanisms of these highly desirable separations remain poorly understood. Herein, through a combination of ab initio simulations and statistical mechanics, it is shown that the temperature dependence is the result of a constraint on the guest molecule's entropic degrees of freedom when loaded into the MOF, caused by the fortuitous tight fitting of the guest inside the pore. While the framework applies to all molecular adsorption in porous media, it is essential for the description of large molecules in small pores, which is demonstrated here using the separation of C6 isomers in Ca(H 2 tcpb) as a test case. The developed framework and analysis not only reveal the reason why separation occurs but also predict the temperatures at which it takes place, thus opening the door to newly designed MOFs with tailor‐made precision.

Chemistry↗
DOE OSTI2024

Theory of capillary tension and interfacial dynamics of motility-induced phases

The statistical mechanics of equilibrium interfaces has been well-established for over a half century. In the past decade, a wealth of observations have made increasingly clear that a new perspective is required to describe interfaces arbitrarily far from equilibrium. In this work, beginning from microscopic particle dynamics that break time-reversal symmetry, we derive the linear interfacial dynamics of coexisting motility-induced phases. Doing so allows us to identify the athermal energy scale that excites interfacial fluctuations and the nonequilibrium surface tension that resists these excitations. Our theory identifies that, in contrast to equilibrium fluids, this active surface tension contains contributions arising from nonconservative forces which act to suppress interfacial fluctuations and, crucially, is distinct from the mechanical surface tension of Kirkwood and Buff. Here we find that the interfacial stiffness scales linearly with the intrinsic persistence length of the constituent active particle trajectories, in agreement with simulation data. We demonstrate that at wavelengths much larger than the persistence length, the interface obeys surface-area minimizing Boltzmann statistics with our derived nonequilibrium interfacial stiffness playing a role identical to that of equilibrium systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Control-Affine Schrödinger Bridge and Generalized Bohm Potential

From a stochastic control perspective, the Schrödinger bridge is a density-valued continuous curve parameterized by time that connects a given pair of initial and terminal probability densities via minimum effort controlled Brownian motion. The control-affine Schrödinger bridge extends this idea to a generic control-affine Itô diffusion, possibly with an additive state cost. Here, in this letter, we recast the necessary conditions of optimality for the control-affine Schrödinger bridge problem as a two point boundary value problem for a quantum mechanical Schrödinger PDE with complex potential. This complex-valued potential is a generalization of the real-valued Bohm potential in quantum mechanics. Our derived potential is akin to the optical potential in nuclear physics where the real part of the potential encodes elastic scattering (transmission of wave function), and the imaginary part encodes inelastic scattering (absorption of wave function). The key takeaway is that the process noise that drives the evolution of probability densities induces an absorbing medium in the evolution of wave function. These results make new connections between control theory and non-equilibrium statistical mechanics through the lens of quantum mechanics.

Markov processes↗
DOE OSTI2026

Zentropy Theory for Transformative Functionalities of Magnetic and Superconducting Materials

The proposed research developed the zentropy theory through applications to complex magnetic materials and superconductors under the hypothesis that the emergent properties of complex magnetic materials and superconductors can be predicted by statistical mechanics of ergodic microstates with their partition functions computed from DFT-predicted free energies. The key objective is to develop approaches to systematically determine the types and number of microstates and the supercell size in DFT-based calculations through convergency of macroscopic functionalities, with the incorporation of our mixed-space approach accounting for the interactions between periodic supercells. In addition to use scientific intuitions to guide the design of important microstates, the key innovation of the proposed research is to integrate the domain knowledge and the material-property-descriptor database (MPDD) with 4 million microstates, which is supported by our deep neural network machine learning models (SIPFENN: structure-informed prediction of formation energy using neural networks) and integrated with our high throughput DFT Tool Kit (DFTTK). For complex magnetic materials, one of the objectives is to develop approaches to calculate short-range ordering from the statistical distribution of each microstate. For superconductors, the divergency of quasiparticle effective mass at a quantum critical point will be investigated, and the superconducting and non-superconducting microstates will be delineated through analysis of electronic band structure, density of states, charge density, and Fermi surface.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Entanglement Cost for Infinite-Dimensional Physical Systems

We prove that the entanglement cost equals the regularized entanglement of formation for any infinite-dimensional quantum state ρ ΑΒ with finite quantum entropy on at least one of the subsystems A or B. This generalizes a foundational result in quantum information theory that was previously formulated only for operations and states on finite-dimensional systems. The extension to infinite-dimensional systems is nontrivial because the conventional tools for establishing both the direct and converse bounds, i.e., strong typicality, monotonicity, and asymptotic continuity, are no longer directly applicable. To address this problem, we construct a new entanglement dilution protocol for infinite-dimensional states implementable by local operations and a finite amount of one-way classical communication (one-way LOCC), using weak and strong typicality multiple times. We also prove the optimality of this protocol among all protocols, even under infinite-dimensional separable operations, by developing an argument based on alternative forms of monotonicity and asymptotic continuity of the entanglement of formation for infinite-dimensional states. Along the way, we derive a new integral representation for the quantum entropy of infinite-dimensional states, which we believe to be of independent interest. Our results allow us to fully characterize an important operational entanglement measure—the entanglement cost—for all infinite-dimensional physical systems.

Complexity↗
DOE OSTI2025

Revisiting point defect thermodynamics in group IVB and VB transition metal carbides

We present a comprehensive re-examination of point defect thermodynamics in group IVB and VB transition metal carbides (TMCs) with the rocksalt structure using a combination of density functional theory (DFT) calculations and a statistical mechanical Wagner-Schottky model within the canonical ensemble. The most stable configurations of point defects were discovered using basin-hopping global optimization, driven by either a machine learning interatomic potential (MLIP) or DFT. A key finding is the identification of previously unreported dicarbon antisites—a C–C dimer occupying a metal site—as the structural (constitutional) defects on the carbon-rich side of stoichiometry in all group IVB and VB TMCs except TaC. Furthermore, dicarbon antisite-containing thermal defect complexes, such as quadruple and interbranch defects, can dominate in TMCs under specific stoichiometric and temperature conditions. In conclusion, by incorporating dicarbon antisites into the defect landscape, this work provides a revised understanding of the thermodynamics of point defects in TMCs.

Carbides↗
DOE OSTI2026

Lattice-scale variations in viscosity are correlated with solution structure at mineral-water interfaces

At solid-liquid interfaces, the viscosity increases markedly from the bulk due to the collective interactions of ions and water molecules, influencing phenomena relevant to nanofluidics, colloidal dynamics, and electrochemistry. Here, in this study, we investigated dissipative forces at the boehmite-water interface using 3D atomic force microscopy. We observed an increase in interfacial solution viscosity, η, by 10-100-fold as the nanoprobe approached the surface in normal direction, with up to four oscillatory features showing average peaks of η/η bulk = 44-71. Moreover, the viscosity showed sub-nanometer variations within 0.5 nm from the interface, templated by the underlying crystal lattice and correlated with the interfacial solution structure. Beyond a near-wall region of approximately 1.2 nm, the dissipative response was comparable to that in bulk solution. Molecular dynamics simulations, along with statistical mechanical analyses, provided details on hydrodynamic structures near the interface. Specifically, the lattice-dependent dissipative responses are correlated with extensive hydrogen bonding by interfacial water molecules, which increased friction, particularly along the [001] direction. These results demonstrate how solution viscosity at mineral-water interfaces is anisotropic and correlated with the local solution structure, providing insights into the dynamics of nanocrystal attachment.

Viscosity↗
DOE OSTI2025

Deriving the Landauer Principle From the Quantum Shannon Entropy

We derive an expression to determine the equilibrium probability distribution of a quantum state in contact with a noisy thermal environment that formally separates contributions from quantum and classical forms of probabilistic uncertainty. A statistical mechanical interpretation of this probability distribution enables us to derive an expression for the minimum free energy costs for arbitrary (reversible or irreversible) quantum state changes. In conclusion, based on this derivation, we demonstrate that–in contrast to classical systems–the free energy required to erase or reset a qubit depends sensitively on both the fidelity of the target state and on the physical properties of the environment, such as the number of quantum bath states, due primarily to the entropic effects of system-bath entanglement.

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

Molecular Dynamics Simulation and Theoretical Analysis of Structural Relaxation, Bond Exchange Dynamics, and Glass Transition in Vitrimers

Vitrimers are a class of polymer networks featuring dynamic covalent cross-links that can undergo associative bond exchange. There has been recent interest in these materials due to their promise as recyclable thermosets or self-healing polymers because of the ability of vitrimer networks to rearrange at the molecular level and undergo macroscopic flow. However, the practical use of these materials often occurs in the supercooled regime or glassy state, where the implications of dynamic bonds are complicated by the interplay between slow activated segmental dynamics, cross-link (i.e., bond-exchange) kinetics, and ultimately material properties. Here, in this paper, we combine coarse-grained molecular dynamics simulation and microscopic statistical mechanical theory to understand how cross-linking kinetics affect material dynamics and how this couples to segmental relaxation of the polymeric network strands across a spectrum of length and time scales, especially in the supercooled regime. We characterize the Kuhn segmental alpha relaxation time and bond exchange time for vitrimer systems across various cross-link densities, temperatures, and bond exchange rates. Simulation and theory both exhibit a bending-up behavior for bond exchange time upon cooling, suggesting a coupling between bond exchange dynamics and segmental relaxation that intensifies with faster bond exchange kinetics. We also found bond exchange dynamics have an impact on Kuhn segment alpha relaxation time, which is most significant at higher cross-link densities. Both these effects are most prominent when the bond exchange time is similar to the Kuhn segment alpha relaxation time, and the resulting coupling of these two relaxation processes is tied to both the probability of a free end to find a bonded pair and the time scale of the constraints imposed by the dynamic cross-links. This relationship is reflected by a cross-link dependence of a theoretical parameter which represents the quantitative degree of coupling between bond exchange and segmental dynamics. Overall, the combination of simulation and theory clarifies the intricate interaction between bond kinetics and segmental relaxation and demonstrates the ability to provide molecular-level insights into vitrimer dynamics over a wide temperature range.

dynamic relaxation↗
DOE OSTI2025

Segmental Dynamics and Vitrification in Associating Copolymer Melts: Role of Cluster Formation, Microdomains, and Cross-Linking

A microscopic statistical mechanical theory of the structure, self-assembly, and activated segmental relaxation is employed to study associating copolymer melts with high attractive sticker fractions, local clustering, and disordered microphase ordering. The stickers are dynamically pinned in a manner that does not affect equilibrium structure which mimics the much slower physical bond breaking process or postassembly cross-linking of sticky monomers. Local sticker clustering and microdomain spatial correlations significantly modify the activated relaxation of nonstickers and glass transition temperature, T g . A re-entrant glass-melting feature is predicted as sticker attraction strength is initially increased corresponding to a speed up of segmental relaxation, and hence reduction of T g relative to the cross-linked homopolymer network. A mechanistic analysis reveals three competing effects: a purely kinetic slowing down of nonstickers down due to cross-linking, disordering of the nonsticker local cage and weakening of effective forces they experience due to sticker physical clustering, and a longer range impact of microdomain scale correlations that results in nonmonotonic dynamical effects. At high enough attraction strength, a qualitative change emerges corresponding to a sticker fraction dependent elevation of T g , which eventually surpasses that of the cross-linked homopolymer network. Dynamically, the new physics arises from a complex evolution of the amplitude of the collective elastic field that dresses the large amplitude mobile segment hopping within a coupled local-nonlocal description of the alpha relaxation. Furthermore, the results are qualitatively consistent with recent experiments on associating PDMS and PPG telechelics of fixed sticker fraction but with chemically different end groups of variable attraction strengths. Finally, possible tests using simulation and the influence of material or model specific interaction potentials and other real world complications are discussed.

Copolymers↗
DOE OSTI2026

Tunable Nanostructures from Inverse Surfactants

Hierarchical materials in the natural world are often made through the self-assembly of amphiphilic molecules. Achieving similar structural complexity in synthetic materials requires understanding how various molecular parameters affect assembly behavior. In recent years, inverse surfactants─molecules with hydrophobic head groups and hydrophilic macromolecular tails─have been shown to self-assemble into supramolecular assemblies in aqueous solutions that show promise for a number of applications, including drug delivery. Here, we build an understanding of the morphological phase diagram of inverse surfactants using insights from scattering experiments, computer simulations, and statistical mechanics. The scattering and simulation results reveal that changing the headgroup size is an important molecular knob in controlling morphological transitions. The molecular size ratio of the hydrophobic group to the hydrophilic group emerges as a crucial dimensionless quantity in our theory and plays a determining role in setting the micelle structure and the transition from mesoscale to macroscale aggregates. Our minimal theory is able to qualitatively explain the key features of the morphological phase diagram, including the prevalence of fiber-like structures in comparison to spherical and planar micelles. Together, these findings provide a more complete picture of the molecular dependencies of assemblies of inverse surfactants, which we hope may aid in the de novo design of supramolecular structures.

Christakopoulos, Panagiotis [ORNL] (ORCID:00090004↗
DOE OSTI2026

An Atom-Precise Approach to Damp First-Order Phase Transitions and Its Implications for Neuromorphic Signal Processing

Neuromorphic computing inspired by mammalian intelligence aims to emulate the nonlinear dynamics of biological neurons and synapses to achieve fast, low-energy, and highly efficient information processing. Brain-inspired computing relies on the design and discovery of materials exhibiting nonlinear current–voltage profiles, frequently underpinned by electronic state transitions, to achieve spiking neurons and dynamically tunable synapses. A signature challenge in the design of artificial neurons is controlling the steepness of first-order transitions in active elements, as abrupt transitions are at risk of driving unstable voltage and temperature oscillations, which result in catastrophic device failure. A critical knowledge gap is the lack of structure–function correlations mapping the composition and atomistic structure of crystalline solids to nonlinear dynamical response characteristics. Here, we address the key question of how modification of atomistic structure correlates with alteration of neuron-like functionality. Constructing oscillator circuits from millimeter-scale single crystals enables high-resolution atomic structure solutions, which we use to demonstrate that the selective positioning of Pb cations modifies charge ordering along a one-dimensional CuxV2O5 framework even at low insertion stoichiometries, thereby providing an atom-precise design parameter for damping first-order transitions. We use temperature-variant X-ray diffraction and X-ray spectroscopy to elucidate the suppression of Cu-ion shuttling based on the precise positioning of Pb ions in seven-coordinated tunnel interstitial sites as the mechanistic basis for transition broadening, thus bridging a critical gap between statistical mechanics and quantum chemical descriptions of phase transitions. Such mechanistic understanding thus paves the way to site-selective modification strategies for modulating the sharpness of first-order transitions, with an exemplary demonstration here in tuning neuronal signal processing.

Crystal structure↗
DOE OSTI2025

Photon–photon chemical thermodynamics of frequency conversion processes in highly multimode systems

Abstract Frequency generation in highly multimode nonlinear optical systems is inherently a complex process, giving rise to an exceedingly convoluted landscape of evolution dynamics. While predicting and controlling the global conversion efficiencies in such nonlinear environments has long been considered impossible, here, we formally address this challenge even in scenarios involving a very large number of spatial modes. By utilizing fundamental notions from optical statistical mechanics, we develop a universal theoretical framework that effectively treats all frequency components as chemical reactants/products, capable of undergoing optical thermodynamic reactions facilitated by a variety of multi-wave mixing effects. These photon–photon reactions are governed by conservation laws that directly determine the optical temperatures and chemical potentials of the ensued chemical equilibria for each frequency species. In this context, we develop a comprehensive stoichiometric model and formally derive an expression that relates the chemical potentials to the optical stoichiometric coefficients, in a manner akin to atomic/molecular chemical reactions. This advancement unlocks new predictive capabilities that can facilitate the optimization of frequency generation in highly multimode photonic arrangements, surpassing the limitations of conventional schemes that rely exclusively on nonlinear optical dynamics. Notably, we identify a universal regime of Rayleigh–Jeans thermalization where an optical reaction at near-zero optical temperatures can promote the complete and entropically irreversible conversion of light to the fundamental mode at a target frequency. Our theoretical results are corroborated by numerical simulations in settings where second-harmonic generation, sum-frequency generation and four-wave mixing processes can manifest.

Optics↗
DOE OSTI2024

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

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

Transport coefficient approach for characterizing nonequilibrium dynamics in soft matter

Nonequilibrium states in soft condensed matter require a systematic approach to characterize and model materials, enhancing predictability and applications. Among the tools, X-ray photon correlation spectroscopy (XPCS) provides exceptional temporal and spatial resolution to extract dynamic insight into the properties of the material. However, existing models might overlook intricate details. We introduce an approach for extracting the transport coefficient, denoted as $J(t)$, from the XPCS studies. This coefficient is a fundamental parameter in nonequilibrium statistical mechanics and is crucial for characterizing transport processes within a system. Our method unifies the Green–Kubo formulas associated with various transport coefficients, including gradient flows, particle–particle interactions, friction matrices, and continuous noise. We achieve this by integrating the collective influence of random and systematic forces acting on the particles within the framework of a Markov chain. We initially validated this method using molecular dynamics simulations of a system subjected to changes in temperatures over time. Subsequently, we conducted further verification using experimental systems reported in the literature and known for their complex nonequilibrium characteristics. The results, including the derived $J(t)$ and other relevant physical parameters, align with the previous observations and reveal detailed dynamical information in nonequilibrium states. This approach represents an advancement in XPCS analysis, addressing the growing demand to extract intricate nonequilibrium dynamics. Further, the methods presented are agnostic to the nature of the material system and can be potentially expanded to hard condensed matter systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

A kinetic-based regularization method for data science applications

We propose a physics-based regularization technique for function learning, inspired by statistical mechanics. By drawing an analogy between optimizing the parameters of an interpolator and minimizing the energy of a system, we introduce corrections that impose constraints on the lower-order moments of the data distribution. This minimizes the discrepancy between the discrete and continuum representations of the data, in turn allowing to access more favorable energy landscapes, thus improving the accuracy of the interpolator. Our approach improves performance in both interpolation and regression tasks, even in high-dimensional spaces. Unlike traditional methods, it does not require empirical parameter tuning, making it particularly effective for handling noisy data. We also show that thanks to its local nature, the method offers computational and memory efficiency advantages over Radial Basis Function interpolators, especially for large datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Convex relaxation for Fokker–Planck equation

We propose an approach to directly estimate the moments or marginals for a high-dimensional equilibrium distribution in statistical mechanics by solving the high-dimensional Fokker–Planck equation in terms of low-order cluster moments or marginals. With this approach, we bypass the exponential complexity of estimating the full high-dimensional distribution and directly solve the simplified partial differential equations for low-order moments/marginals. Moreover, the proposed moment/marginal relaxation is fully convex and can be solved via off-the-shelf solvers. We further propose a time-dependent version of the convex programs to study non-equilibrium dynamics. In a specific setting, we show the proposed method can recover a mean-field-type equilibrium density. Numerical results are provided to demonstrate the performance of the proposed algorithm for high-dimensional systems.

Chen, Yian↗