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 73 records · Page 4

Data-driven learning of Mori–Zwanzig operators for isotropic turbulence

Developing reduced-order models for turbulent flows, which contain dynamics over a wide range of scales, is an extremely challenging problem. In statistical mechanics, the Mori–Zwanzig (MZ) formalism provides a mathematically exact procedure for constructing reduced-order representations of high-dimensional dynamical systems, where the effects due to the unresolved dynamics are captured in the memory kernel and orthogonal dynamics. Turbulence models based on MZ formalism have been scarce due to the limited knowledge of the MZ operators, which originates from the difficulty in deriving MZ kernels for complex nonlinear dynamical systems. In this work, we apply a recently developed data-driven learning algorithm, which is based on Koopman's description of dynamical systems and Mori's linear projection operator, on a set of fully resolved isotropic turbulence datasets to extract the Mori–Zwanzig operators. With data augmentation using known turbulence symmetries, the extracted Markov term, memory kernel, and orthogonal dynamics are statistically converged and the generalized fluctuation–dissipation relation can be verified. The properties of the memory kernel and orthogonal dynamics, and their dependence on the choices of observables are investigated to address the modeling assumptions that are commonly used in MZ-based models. A series of numerical experiments are then constructed using the extracted kernels to evaluate the memory effects on prediction. The results show that the prediction errors are strongly affected by the choice of observables and can be further reduced by including the past history of the observables in the memory kernel.

97 MATHEMATICS AND COMPUTING↗

Linear and nonlinear viscoelasticity of concentrated thermoresponsive microgel suspensions

We present an integrated experimental and theoretical study of the dynamics and rheology of self-crosslinked, slightly charged, temperature responsive soft poly(N-isopropylacrylamide) (pNIPAM) microgels over a wide range of concentration and temperature spanning the sharp change in particle size and intermolecular interactions across the lower critical solution temperature (LCST). Dramatic, non-monotonic changes in viscoelasticity are observed as a function of temperature, with distinct concentration dependence in the dense fluid, glassy, and soft-jammed regimes. Motivated by our experimental observations, we formulate a minimalistic model for the size dependence of a single microgel particle and the change of the interparticle interaction from purely repulsive to attractive upon heating. Using microscopic equilibrium and time-dependent statistical mechanical theories, theoretical predictions are quantitatively compared with experimental measurements of the shear modulus. Good agreement is found for the nonmonotonic temperature behavior that originates as a consequence of the competition between reduced microgel packing fraction and increasing interparticle attractions. Furthermore, testable predictions are made for nonlinear rheological properties such as the yield stress and strain. To our knowledge, this is the first attempt to quantitatively understand in a unified manner the viscoelasticity of dense, temperature-responsive microgel suspensions spanning a wide range of temperatures and concentrations.

42 ENGINEERING↗

Remembering the Work of Phillip L. Geissler: A Coda to His Scientific Trajectory

Phillip L. Geissler made important contributions to the statistical mechanics of biological polymers, heterogeneous materials, and chemical dynamics in aqueous environments. He devised analytical and computational methods that revealed the underlying organization of complex systems at the frontiers of biology, chemistry, and materials science. In this retrospective we celebrate his work at these frontiers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular-based analysis of nanoparticle solvation: Classical density functional approach

Proper statistical mechanics understanding of nanoparticle solvation processes requires an accurate description of the molecular structure of the solvent. Achieving this goal with standard molecular dynamics (MD) simulation methods is challenging due to large length scales. An alternative approach to this problem can be formulated using classical density functional theory (cDFT), where a full configurational description of the positions of all the atoms is replaced by collective atomic site densities in the molecule. Using an example of the negatively charged silica-like system in an aqueous polar environment represented by a two-site water model, we demonstrate here that cDFT can reproduce MD data at a fraction of the computational cost. An important implication of this result is the ability to understand how the solvent molecular features may affect the system’s properties at the macroscopic scale. A concrete example highlighted in this work is the analysis of nanoparticle interactions with sizes of up to 100 nm in diameter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamically consistent incorporation of the Langmuir adsorption model into compressible fluctuating hydrodynamics

For a gas–solid interfacial system where chemical species undergo reversible adsorption, we develop a mesoscopic stochastic modeling method that simulates both gas-phase hydrodynamics and surface coverage dynamics by coupling the Langmuir adsorption model with compressible fluctuating hydrodynamics. To this end, we derive a thermodynamically consistent mass–energy update scheme that accounts for how the mass and energy variables in the gas and surface subsystems should be updated according to the changes in the number of molecules of each species in each subsystem due to adsorption and desorption events. By performing a stochastic analysis for the ideal Langmuir model and the full hydrodynamic system, we analytically confirm that our mass–energy update scheme captures thermodynamic equilibrium predicted by equilibrium statistical mechanics. We find that an internal energy correction term is needed, which is attributed to the difference in the mean kinetic energy of gas molecules colliding with the surface from that computed from the Maxwell–Boltzmann distribution. By performing an equilibrium simulation study for an ideal gas mixture of CO and Ar, with CO undergoing reversible adsorption, we validate our overall simulation method and implementation.

Adsorption↗

Breaking the curse of dimensionality: Solving configurational integrals for crystalline solids by tensor networks

Accurately evaluating configurational integrals for dense solids remains a central and difficult challenge in the statistical mechanics of condensed systems. Here, we present a tensor network approach that reformulates the high-dimensional configurational integral for identical-particle crystals into a sequence of computationally efficient summations. We represent the integrand as a high-dimensional tensor and apply tensor-train (TT) decomposition together with a custom TT-cross interpolation. This approach circumvents the need to explicitly construct the full tensor. We introduce tailored rank-1 and rank-2 schemes optimized for sharply peaked Boltzmann probability densities, typical for identical-particle crystals. When applied to the calculation of internal energy and pressure-temperature curves for crystalline Cu and Ar at high (GPa) pressures, as well as the alpha-to-beta phase transition diagram of Sn, our method accurately reproduces molecular dynamics simulation results using tight-binding, machine learning, hierarchical interacting particle–neural network, and modified embedded atom method potentials,all within seconds of computation time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Self-consistent hopping theory of activated relaxation and diffusion of dilute penetrants in dense crosslinked polymer networks

We generalize a microscopic statistical mechanical theory of the activated dynamics of dilute spherical penetrants in glass-forming liquids to study the influence of crosslinking in polymer networks on the penetrant relaxation time and diffusivity over a wide range of temperature and crosslink fraction (f n ). Our calculations are relevant to recent experimental studies of a nm-sized molecule diffusing in poly-(n-butyl methacrylate) networks. The theory predicts the penetrant relaxation time increases exponentially with the glass transition temperature, T g (f n ), which grows roughly linearly with the square root of fn due to the coupling of local hopping to longer-range collective elasticity. Moreover, T g is also found to be proportional to a geometric confinement parameter defined as the ratio of the penetrant diameter to the mean network mesh size. The decoupling ratio of the penetrant and Kuhn segment alpha times displays a complex non-monotonic dependence on f n and temperature that is well collapsed based on the variable T g (f n )/T. A model for the penetrant diffusion constant that combines activated relaxation and entropic mesh confinement is proposed, which results in a significantly stronger suppression of mass transport with degree of effective supercooling than predicted for the penetrant alpha time. Furthermore, this behavior corresponds to a new network-based type of “decoupling” of diffusion and relaxation. In contrast to the diffusion of larger nanoparticles in high temperature rubbery networks, our analysis in the supercooled regime suggests that for the penetrants studied the mesh confinement effects are of secondary importance relative to the consequences of crosslink-induced slowing down of activated hopping of glassy physics origin.

Chemistry↗

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↗

Reversible to irreversible transitions for cyclically driven particles on periodic obstacle arrays

Here, we examine the collective dynamics of disks moving through a square array of obstacles under cyclic square wave driving. Below a critical density, we find that the system organizes into a reversible state in which the disks return to the same positions at the end of every drive cycle. Above this density, the dynamics are irreversible and the disks do not return to the same positions after each cycle. The critical density depends strongly on the angle θ between the driving direction and a symmetry axis of the obstacle array, with the highest critical densities appearing at commensurate angles such as θ = 0° and θ = 45° and the lowest critical densities falling at θ = arctan (0.618), the inverse of the golden ratio, where the flow is the most degenerate. As the density increases, the number of cycles required to reach a reversible state grows as a power law with an exponent near ν = 1.36, similar to what is found in periodically driven colloidal and superconducting vortex systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Assessment of software methods for estimating protein-protein relative binding affinities

A growing number of computational tools have been developed to accurately and rapidly predict the impact of amino acid mutations on protein-protein relative binding affinities. Such tools have many applications, for example, designing new drugs and studying evolutionary mechanisms. In the search for accuracy, many of these methods employ expensive yet rigorous molecular dynamics simulations. By contrast, non-rigorous methods use less exhaustive statistical mechanics, allowing for more efficient calculations. However, it is unclear if such methods retain enough accuracy to replace rigorous methods in binding affinity calculations. This trade-off between accuracy and computational expense makes it difficult to determine the best method for a particular system or study. Here, eight non-rigorous computational methods were assessed using eight antibody-antigen and eight non-antibody-antigen complexes for their ability to accurately predict relative binding affinities (ΔΔG) for 654 single mutations. In addition to assessing accuracy, we analyzed the CPU cost and performance for each method using a variety of physico-chemical structural features. This allowed us to posit scenarios in which each method may be best utilized. Most methods performed worse when applied to antibody-antigen complexes compared to non-antibody-antigen complexes. Rosetta-based JayZ and EasyE methods classified mutations as destabilizing (ΔΔG < -0.5 kcal/mol) with high (83–98%) accuracy and a relatively low computational cost for non-antibody-antigen complexes. Some of the most accurate results for antibody-antigen systems came from combining molecular dynamics with FoldX with a correlation coefficient (r) of 0.46, but this was also the most computationally expensive method. Overall, our results suggest these methods can be used to quickly and accurately predict stabilizing versus destabilizing mutations but are less accurate at predicting actual binding affinities. This study highlights the need for continued development of reliable, accessible, and reproducible methods for predicting binding affinities in antibody-antigen proteins and provides a recipe for using current methods.

59 BASIC BIOLOGICAL SCIENCES↗

The response field and the saddle points of quantum mechanical path integrals

Highlights: • Moyal quantum mechanics and Marinov’s path integral. • Classical and semiclassical limits of Marinov’s path integral. • Oscillating functional integrals. • Instantons of the Marinov’s path integral. In quantum statistical mechanics, Moyal’s equation governs the time evolution of Wigner functions and of more general Weyl symbols that represent the density matrix of arbitrary mixed states. A formal solution to Moyal’s equation is given by Marinov’s path integral. In this paper we demonstrate that this path integral can be regarded as the natural link between several conceptual, geometric, and dynamical issues in quantum mechanics. A unifying perspective is achieved by highlighting the pivotal role which the response field, one of the integration variables in Marinov’s integral, plays for pure states even. The discussion focuses on how the integral’s semiclassical approximation relates to its strictly classical limit; unlike for Feynman type path integrals, the latter is well defined in the Marinov case. The topics covered include a random force representation of Marinov’s integral based upon the concept of “Airy averaging”, a related discussion of positivity-violating Wigner functions describing tunneling processes, and the role of the response field in maintaining quantum coherence and enabling interference phenomena. The double slit experiment for electrons and the Bohm–Aharonov effect are analyzed as illustrative examples. Furthermore, a surprising relationship between the instantons of the Marinov path integral over an analytically continued (“Wick rotated”) response field, and the complex instantons of Feynman-type integrals is found. The latter play a prominent role in recent work towards a Picard–Lefschetz theory applicable to oscillatory path integrals and the resurgence program.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Geometric and topological entropies of sphere packing

Here, we present a statistical mechanical description of randomly packed spherical particles, where the average coordination number is treated as a macroscopic thermodynamic variable. The overall packing entropy is shown to have two contributions: geometric, reflecting statistical weights of individual configurations, and topological, which corresponds to the number of topologically distinct states. Both of them are computed in the thermodynamic limit for isostatic and weakly underconstrained packings in 2D and 3D. The theory generalizes concepts of granular and glassy configurational entropies for the case of nonjammed systems. It is directly applicable to sticky colloids and predicts an asymptotic phase behavior of sticky spheres in the limit of strong binding.

36 MATERIALS SCIENCE↗

Thermodynamics of order and randomness in dopant distributions inferred from atomically resolved imaging

Abstract Exploration of structure-property relationships as a function of dopant concentration is commonly based on mean field theories for solid solutions. However, such theories that work well for semiconductors tend to fail in materials with strong correlations, either in electronic behavior or chemical segregation. In these cases, the details of atomic arrangements are generally not explored and analyzed. The knowledge of the generative physics and chemistry of the material can obviate this problem, since defect configuration libraries as stochastic representation of atomic level structures can be generated, or parameters of mesoscopic thermodynamic models can be derived. To obtain such information for improved predictions, we use data from atomically resolved microscopic images that visualize complex structural correlations within the system and translate them into statistical mechanical models of structure formation. Given the significant uncertainties about the microscopic aspects of the material’s processing history along with the limited number of available images, we combine model optimization techniques with the principles of statistical hypothesis testing. We demonstrate the approach on data from a series of atomically-resolved scanning transmission electron microscopy images of Mo x Re 1- x S 2 at varying ratios of Mo/Re stoichiometries, for which we propose an effective interaction model that is then used to generate atomic configurations and make testable predictions at a range of concentrations and formation temperatures.

25 ENERGY STORAGE↗

Linking low and high temperature NO oxidation mechanisms over Brønsted acidic chabazite to dynamic changes of the active site

The oxidation of NO over Brønsted acid sites in chabazite (CHA) zeolites shows an atypical temperature dependence; at low temperature the apparent activation energy is negative, but it becomes positive as the temperature exceeds a transition temperature. To explain this behavior we used density functional theory and statistical mechanics to investigate high and low temperature mechanisms for this reaction and propose a dynamic active site change in response to temperature variation. Our simulations show that the apparent activation barrier in the low temperature regime is more negative over Brønsted acidic CHA as compared to the siliceous zeolite framework. This effect is attributed to further enthalpic stabilization of the transition states by physical interaction with the H-CHA Brønsted acid sites. At elevated temperature, our calculations support both the existence and the significant catalytic role of NO+ in providing a modified active site. The temperature dependent transformation of the active site from H-CHA to NO-CHA sites may occur via two plausible ion-exchange mechanisms that define a transition temperature for the reaction. This transition temperature can be tuned by incorporating different trivalent metal atoms (B, Al, Ga or In) within the CHA framework. We found the lowest transition temperature for H-[In] CHA and H-[Ga]CHA. Here, the ability to control the dynamic response of the active site and the associated switch between low and high temperature mechanism with negative and positive apparent activation energy, respectively, is of fundamental interest for the design of zeolite catalysts operating in the presence of NO.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Collective Nanoparticle Dynamics Associated with Bridging Network Formation in Model Polymer Nanocomposites

The addition of nanoparticles (NPs) to polymers is a powerful method to improve the mechanical and other properties of macromolecular materials. Such hybrid polymer–particle systems are also rich in fundamental soft matter physics. Among several factors contributing to mechanical reinforcement, a polymer-mediated NP network is considered to be the most important in polymer nanocomposites (PNCs). Here, we present an integrated experimental–theoretical study of the collective NP dynamics in model PNCs using X-ray photon correlation spectroscopy and microscopic statistical mechanics theory. Silica NPs dispersed in unentangled or entangled poly(2-vinylpyridine) matrices over a range of NP loadings are used. Static collective structure factors of the NP subsystems at temperatures above the bulk glass transition temperature reveal the formation of a network-like microstructure via polymer-mediated bridges at high NP loadings above the percolation threshold. The NP collective relaxation times are up to 3 orders of magnitude longer than the self-diffusion limit of isolated NPs and display a rich dependence with observation wavevector and NP loading. A mode-coupling theory dynamical analysis that incorporates the static polymer-mediated bridging structure and collective motions of NPs is performed. It captures well both the observed scattering wavevector and NP loading dependences of the collective NP dynamics in the unentangled polymer matrix, with modest quantitative deviations emerging for the entangled PNC samples. Here, we identify an unusual and weak temperature dependence of collective NP dynamics, in qualitative contrast with the mechanical response. Hence, the present study has revealed key aspects of the collective motions of NPs connected by polymer bridges in contact with a viscous adsorbing polymer medium and identifies some outstanding remaining challenges for the theoretical understanding of these complex soft materials.

36 MATERIALS SCIENCE↗

Building materials genome from ground‐state configuration to engineering advance

Individual phases are commonly considered as the building blocks of materials. However, the accurate theoretical prediction of properties of individual phases remains elusive. The top-down approach by decoding genomic building blocks of individual phases from experimental observations is nonunique. The density functional theory (DFT), as a state-of-the-art solution of quantum mechanics, prescribes the existence of a ground-state configuration at 0 K for a given system. It is self-evident that the ground-state configuration alone is insufficient to describe a phase at finite temperatures as symmetry-breaking non-ground-state configurations are excited statistically at temperatures above 0 K. Our multiscale entropy approach (recently terms as Zentropy theory) postulates that the entropy of a phase is composed of the sum of the entropy of each configuration weighted by its probability plus the configurational entropy among all configurations. Consequently, the partition function of each configuration in statistical mechanics needs to be evaluated by its free energy rather than total energy. The combination of the ground-state and symmetry-breaking non-ground-state configurations represents the building blocks of materials and can be used to quantitatively predict free energy of individual phases with the free energy of each configuration predicted from DFT as well as all properties derived from free energy of individual phases.

CALPHAD↗

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↗