Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Isotope effects in supercooled H 2 O and D 2 O and a corresponding-states-like rescaling of the temperature and pressure

Water shows anomalous properties that are enhanced upon supercooling. The unusual behavior is observed in both H 2 O and D 2 O, however, with different temperature dependences for the two isotopes. It is often noted that comparing the properties of the isotopes at two different temperatures (i.e., a temperature shift) approximately accounts for many of the observations—with a temperature shift of 7.2 K in the temperature of maximum density being the most well-known example. However, the physical justification for such a shift is unclear. Motivated by recent work demonstrating a “corresponding-states-like” rescaling for water properties in three classical water models that all exhibit a liquid–liquid transition and critical point, the applicability of this approach for reconciling the differences in the temperature- and pressure-dependent thermodynamic properties of H 2 O and D 2 O is investigated here. Utilizing previously published data and equations-of-state for H 2 O and D 2 O, we show that the available data and models for these isotopes are consistent with such a low temperature correspondence. These observations provide support for the hypothesis that a liquid–liquid critical point, which is predicted to occur at low temperatures and high pressures, is the origin of many of water’s anomalies.

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

Unusual dynamics of tetrahedral liquids caused by the competition between dynamic heterogeneity and structural heterogeneity

Tetrahedral liquids exhibit intriguing thermodynamic and transport properties because of the various ways tetrahedra can be packed and connected. Recently, an unusual temperature dependence of the stretching exponent β in a model tetrahedral liquid ZnCl 2 from T m + 85 K to T m + 35 K has been reported using neutron-spin echo spectroscopy. This discovery stands in sharp contrast to other glass-forming liquids. In this study, we conducted neural network force field driven molecular dynamic simulations of ZnCl 2 . We found a non-monotonic temperature dependence of β from liquid to supercooled liquid temperatures. Further structural decomposition and dynamic analysis suggest that this unusual dynamic behavior is a result of the competition between the decrease in the diversity of tetrahedra motifs (structural heterogeneity) and the increase in glassy dynamic heterogeneity. Furthermore, this result may contribute to new understandings of the structural relaxation of other network liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evidence of short chains in liquid sulfur

High energy x-ray pair distribution function measurements show the average coordination number of the first shell in liquid sulfur is 1.86 ± 0.04 across the λ-transition, not precisely 2.0 as widely accepted. This indicates that upon melting, liquid sulfur does not comprise solely of S 8 rings but also possesses a significant number of short chains. Intensities of the pre-peak and first diffraction peak of the x-ray structure factor and third peak height of the pair distribution function all show deviations at the λ-transition temperature T λ , associated with the break-up of S 8 rings and the start of oligomer polymerization. A significant number of non-bonded or loosely bonded “interstitial atoms,” with an average coordination number of 0.20 ± 0.005, are also observed in the so-called “forbidden zone” between the first and second shells upon melting. The number of interstitial atoms is found to decrease to a minimum at the λ-transition, but the majority persist into the high temperature polymerized liquid. Furthermore, the existence of short chains and nearby interstitial atoms represent the two main factors required to initiate the S 8 -ring to chain transition, as proposed by recent molecular dynamics simulations.

Chemical bonding↗
DOE OSTI2025

Medium-range order and compositional correlation in metallic glasses

The compositional atomic ordering in metallic glasses was studied by simulation focusing on the medium-range order (MRO). Many metallic alloy liquids and glasses show MRO characterized by the oscillations in the atomic pair-distribution function (PDF) beyond the first peak, which decay exponentially with distance. To study the effects of the local chemical order on MRO, we examine the compositionally resolved PDF and its MRO for models of various binary metallic alloy glasses. We show that compositional ordering is limited mostly to the nearest-neighbor atoms and the MRO is largely independent of the compositional order. For some elements that strongly repel each other in the alloy, a second MRO periodicity is observed owing to the distinct correlations among them. These results are discussed in light of the idea that the MRO oscillations in the PDF describe the correlations in the atomic density fluctuations, rather than the detailed local atomic structure.

Atomic structure↗
DOE OSTI2025

Foundations of magnetohydrodynamics

In this tutorial, a derivation of magnetohydrodynamics (MHD) valid beyond the usual ideal gas approximation is presented. Non-equilibrium thermodynamics is used to obtain conservation equations and linear constitutive relations. When coupled with Maxwell's equations, this provides closed fluid equations in terms of material properties of the plasma, described by the equation of state and transport coefficients. These properties are connected to microscopic dynamics using the Irving–Kirkwood procedure and Green–Kubo relations. Symmetry arguments and the Onsager–Casimir relations allow one to vastly simplify the number of independent coefficients. Importantly, expressions for current density, heat flux, and stress (conventionally Ohm's law, Fourier's law, and Newton's law) take different forms in systems with a non-ideal equation of state. The traditional form of the MHD equations, which is usually obtained from a Chapman–Enskog solution of the Boltzmann equation, corresponds to the ideal gas limit of the general equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Charge regulation effects on colloidal mixture nanoparticles

Changes in pH within a system containing dissociable sites affect the protonation and deprotonation of these groups, thereby influencing their physical properties. In response, the system modifies their surface charge, affecting electrostatic interactions, aggregation, stability, and structural behavior. Although the pH can be tuned in experiments, it is difficult to model this phenomenon using simulations or theoretical approaches. Here, we perform hybrid Monte Carlo-molecular dynamics simulations to model charge regulation effects in an equimolar colloidal charged system. We compare charge regulation effects with those of a system in which the charges of colloidal nanoparticles are not dissociable. The comparison between the two cases modifies the phase diagram, and it changes the volume fraction where a percolation network of nanoparticles is found. Charge regulation is found to destroy network formation, as the charge in the nanoparticles is modified because of the cooperativity dependency of the degree of charge dissociation sites among the nanoparticles favoring cluster formation. Furthermore, our work suggests that the ionic and/or electronic conductivity in functionalized nanoparticles can be modified by changing pH values. It also guides the experimental design of oppositely charged nanoparticles as inks for 3D printing processes.

Classical statistical mechanics↗
DOE OSTI2025

Minimal implicit-solvent coarse-grained simulation of Pluronic block copolymers with ionic liquids

Pluronic block copolymers, composed of poly(ethylene oxide) (PEO) and poly(propylene oxide) (PPO) in a triblock structure (PEO–PPO–PEO), are well known for their amphiphilic character and ability to self‐assemble into micelles in aqueous solution. The addition of ionic liquids (ILs) can further modulate the core–shell structures of these copolymers, influencing their stability, critical micellization temperature, and size. However, fully atomistic simulations often become prohibitively expensive due to the size and complexity of these systems. In this work, coarse‐grained simulations using a minimal implicit‐solvent model were performed to examine how two classes of ILs, namely, 1‐alkyl‐3‐methylimidazolium ([C n C 1 im]) and 1‐alkyl‐3‐methylpyrrolidinium ([C n C 1 pyrr]), change the micellization of Pluronic block copolymers in aqueous solution. The effects of IL concentration and alkyl group length were investigated, and the model greatly improved the efficiency of simulating large‐scale micelle systems. Furthermore, the numerical simulations are qualitatively compared with experimental investigations. Our results show that adding ILs expands the micelle core by embedding IL tails among the PPO blocks, thereby increasing overall micelle size. Less polar ILs generally induce more pronounced micellar growth. However, the effect of IL tail length on conformation and micellar packing is non‐monotonic. Up to moderate chain lengths (around C8–C10), the IL tails can extend sufficiently to increase local separation within the micelle; at longer tail lengths, enhanced hydrophobic clustering and steric hindrance cause the tails to bend or fold, capping further expansion. In addition, although block copolymer chains tend to pack more closely in the presence of longer‐tailed ILs, the random coil size of an individual polymer chain does not necessarily shrink. Meanwhile, these insights provide a deeper understanding of how Pluronic/IL systems interact, informing applications in drug delivery, cosmetics, food, and environmental engineering. Finally, our minimal implicit‐solvent model can be applied to larger systems and longer timescales, substantially reducing computational cost while reproducing key structural trends observed experimentally.

Atomistic simulations↗
DOE OSTI2026

Optimizing time integration for accurate recovery of shockwave interface location in radiography

We present simulations and experiments of time integrated radiographic imaging of a moving 1D shock wave front and a quantitative method for determining the statistical error in locating the shock front as a function of integration time and noise in the radiograph. We discuss the trade-off between increasing motion blur, which leads to decreased shock front location certainty, and increasing signal-to-noise, which leads to improved image quality with increasing integration time. We find an optimum integration time between a short integration time, where noise limits the error, and a long integration time, where motion blurring limits the error. This methodology can be used to tune experimental configurations to obtain the highest quality radiograph for a given experimental configuration.

Bremsstrahlung↗
DOE OSTI2026

Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas

Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

Electrical conductivity↗
DOE OSTI2026

Plasma hydrodynamics from mean force kinetic theory

Mean force kinetic theory is used to evaluate the electrical conductivity, thermal conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity of a two-component plasma. The results are compared with molecular dynamics simulations. These simulations are made possible by assuming a repulsive Coulomb force for all interactions. Good agreement is found for all coefficients up to a Coulomb coupling strength of Γ ≈ 20. This is over 100-times larger than the coupling strength at which traditional theories break down. It is concluded that mean force kinetic theory is a promising framework for extending hydrodynamics to dense plasmas.

Electrical conductivity↗
DOE OSTI2026

An interactive machine learning platform for analyzing multi-particle coincidence data from cold target recoil ion momentum spectroscopy

We present SCULPT (Supervised Clustering and Uncovering Latent Patterns with Training), a comprehensive software platform for analyzing tabulated high-dimensional multi-particle coincidence data from Cold Target Recoil Ion Momentum Spectroscopy (COLTRIMS) experiments. The software addresses critical challenges in modern momentum spectroscopy by integrating advanced machine learning techniques with physics-informed analysis in an interactive web-based environment. SCULPT implements uniform manifold approximation and projection for non-linear dimensionality reduction to reveal correlations in high-dimensional data. We also discuss potential extensions to deep autoencoders for feature learning and genetic programming for automated discovery of physically meaningful observables. A novel adaptive confidence scoring system provides quantitative reliability assessments by evaluating user-selected clustering quality metrics with predefined weights that reflect each metric’s robustness. The platform features configurable molecular profiles for different experimental systems, interactive visualization with selection tools, and comprehensive data filtering capabilities. Utilizing a subset of SCULPT’s capabilities, we analyze photo-double-ionization data measured using the COLTRIMS method for three-body dissociation of the D 2 O molecule, revealing distinct fragmentation channels and their correlations with physics parameters. The software’s modular architecture and web-based implementation make it accessible to the broader atomic and molecular physics community, significantly reducing the time required for complex multi-dimensional analyses. This opens the door to finding and isolating rare events exhibiting non-linear correlations on the fly during experimental measurements, which can help steer exploration and improve the efficiency of experiments.

Artificial neural networks↗
DOE OSTI2022

Nonequilibrium statistical thermodynamics of multicomponent interfaces

Nonequilibrium interfacial thermodynamics has important implications for crucial biological, physical, and industrial-scale transport processes. Here, we discuss a theory of local equilibrium for multiphase multicomponent interfaces that builds upon the “sharp” interface concept first introduced by Gibbs, allowing for a description of nonequilibrium interfacial processes such as those arising in evaporation, condensation, adsorption, etc. By requiring that the thermodynamics be insensitive to the precise location of the dividing surface, one can identify conditions for local equilibrium and develop methods for measuring the values of intensive variables at the interface. We then use extensive, high-precision nonequilibrium molecular dynamics (NEMD) simulations to verify the theory and establish the validity of the local equilibrium hypothesis. In particular, we demonstrate that equilibrium equations of state are also valid out of equilibrium, and can be used to determine interfacial temperature and chemical potential(s) that are consistent with nonequilibrium generalizations of the Clapeyron and Gibbs adsorption equations. We also show, for example, that, far from equilibrium, temperature or chemical potential differences need not be uniform across an interface and may instead exhibit pronounced discontinuities. However, even in these circumstances, we demonstrate that the local equilibrium hypothesis and its implications remain valid. These results provide a thermodynamic foundation and computational tools for studying or revisiting a wide variety of interfacial transport phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Adaptive nonequilibrium design of actin-based metamaterials: Fundamental and practical limits of control

The adaptive and surprising emergent properties of biological materials self-assembled in far-from-equilibrium environments serve as an inspiration for efforts to design nanomaterials. In particular, controlling the conditions of self-assembly can modulate material properties, but there is no systematic understanding of either how to parameterize external control or how controllable a given material can be. Here, we demonstrate that branched actin networks can be encoded with metamaterial properties by dynamically controlling the applied force under which they grow and that the protocols can be selected using multi-task reinforcement learning. These actin networks have tunable responses over a large dynamic range depending on the chosen external protocol, providing a pathway to encoding “memory” within these structures. Interestingly, we obtain a bound that relates the dissipation rate and the rate of “encoding” that gives insight into the constraints on control—both physical and information theoretical. Taken together, these results emphasize the utility and necessity of nonequilibrium control for designing self-assembled nanostructures.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Arrested coarsening and large density fluctuations in driven particle mixtures in two dimensions

Abstract Using molecular dynamics simulations, we study a driven, nonadditive binary mixture of spherical particles confined to move in two dimensions and immersed in an explicit solvent consisting of point particles with purely repulsive interactions. We show that, without a drive, the mixture of spherical particles phase separates and coarsens with kinetics consistent with an Ising-like conserved dynamics. Conversely, when the drive is applied, the coarsening is arrested and the system develops large density fluctuations. We show that the drive creates domains of a characteristic size which decreases with an increasing force. Furthermore, we find that these domains are anisotropic and can be oriented either parallel or perpendicular to the drive direction. Finally, we connect our findings to existing theories of strongly-driven systems, pointing out the importance of introducing the explicit solvent particles to break the Galilean invariance of the system.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Generalized Wigner theorem for noninvertible symmetries

In this article, we establish the conditions under which a conservation law associated with a non-invertible operator may be realized as a symmetry in quantum physics. As established by Wigner, all quantum symmetries must be represented by either unitary or antiunitary transformations. Relinquishing an implicit assumption of invertibility, we demonstrate that the fundamental invariance of quantum transition probabilities under the application of symmetries mandates that all non-invertible symmetries may only correspond to projective unitary or antiunitary transformations, i.e., partial isometries. This extends the notion of physical states beyond conventional rays in Hilbert space to equivalence classes in an extended, gauged Hilbert space, thereby broadening the traditional understanding of symmetry transformations in quantum theory. Our generalized theorem applies irrespective of the origin of the (non)invertible symmetry, holds in arbitrary spatial dimensions, and is independent of the Hamiltonian or action. We explore its physical consequences and, using simple model systems, illustrate how the distinction between invertible and non-invertible symmetries can sometimes be tied to the choice of boundary conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Strong Hilbert space fragmentation and fractons from subsystem and higher-form symmetries

Here, we introduce a route to Hilbert space fragmentation in high dimensions leveraging the group-word formalism. We show that taking strongly fragmented models in one dimension and “lifting” to higher dimensions using subsystem symmetries can yield strongly fragmented dynamics in higher dimensions, with subdimensional (e.g., lineonic) excitations. This provides a route to higher-dimensional strong fragmentation, and also a route to fractonic behavior. Meanwhile, lifting one-dimensional fragmented models to higher dimensions using higher-form symmetries yields models with topologically robust fragmentation. In three or more spatial dimensions, one can also “mix and match” subsystem and higher-form symmetries, leading to canonical fracton models such as X cube. We speculate that this approach could also yield a route to non-Abelian fractons. These constructions unify a number of phenomena that have been discussed in the literature, as well as furnishing models with unique properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Generalized geometric speed limits for quantum observables

Leveraging quantum information geometry, we derive generalized quantum speed limits on the rate of change of the expectation values of observables. These bounds subsume and, for Hilbert space dimension ≥3, tighten existing bounds—in some cases by an arbitrarily large multiplicative constant. Our theoretical results are supported by illustrative examples and an experimental demonstration using a superconducting qutrit. We also derive two upper bounds on the generalized quantum Fisher information in terms of the condition number of the density matrix. One of these bounds applies only to coherent dynamics and depends also on the variance of the Hamiltonian. The other bound depends also on the so-called Wigner-Yanase skew information. These bounds generalize well-known bounds on the symmetric logarithmic derivative quantum Fisher information and are tighter than the existing bounds for sufficiently mixed states (e.g., for sufficiently high temperature thermal states).

open quantum systems & decoherence↗
DOE OSTI2026

Nonequilibrium universality of the nonreciprocally coupled 𝑂⁡(𝑛 1 ) × 𝑂⁡(𝑛 2 ) model

Nonequilibrium dynamics play an important role in all contexts of physics, both classical and quantum as well as living and nonliving, so it is crucial to develop a foundational understanding of nonequilibrium phase transitions. In this work we investigate an important class of nonequilibrium dynamics in the form of nonreciprocal interactions. In particular we study how nonreciprocal coupling between two 𝑂⁡(𝑛𝑖) order parameters (with 𝑖 = 1,2) affects the universality at a multicritical point, extending the analysis of J. T. Young et al. [Phys. Rev. X 10, 011039 (2020)], which considered the case 𝑛 1 = 𝑛 2 = 1, i.e., a ℤ 2 × ℤ 2 model. We show that nonequilibrium fixed points (NEFPs) emerge for a broad range of 𝑛 1 ,𝑛 2 and exhibit intrinsically nonequilibrium critical phenomena, namely a violation of fluctuation-dissipation relations at all scales and underdamped oscillations near criticality in contrast to the overdamped relaxational dynamics of the corresponding equilibrium models. Furthermore, the NEFPs exhibit an emergent discrete scale invariance in certain physically relevant regimes of 𝑛 1 ,𝑛 2 , but not others, depending on whether the critical exponent 𝜈 is real or complex. The boundary between these two regions is described by an exceptional point in the renormalization group (RG) flow, leading to distinctive features in correlation functions and the phase diagram. Another contrast with the previous work is the number and stability of the NEFPs as well as the underlying topology of the RG flow. Lastly, we investigate an extreme form of nonreciprocity where one order parameter is independent of the other order parameter but not vice versa. Unlike the ℤ 2 × ℤ 2 model, which becomes nonperturbative in this case, we identify a distinct nonequilibrium universality class whose dependent field similarly violates fluctuation-dissipation relations but does not exhibit discrete scale invariance or underdamped oscillations near criticality.

Critical phenomena↗