Predicting continuous and discontinuous phase decompositions using steepest-entropy-ascent quantum thermodynamics

Predicting Polymer Brush Behavior in Solvents Using the Steepest-Entropy-Ascent Quantum Thermodynamic Framework

The steepest-entropy-ascent quantum thermodynamic (SEAQT) framework is utilized to study the effects of temperature on polymer brushes. The brushes are represented by a discrete energy spectrum, and energy degeneracies obtained through the replica-exchange Wang-Landau algorithm. The SEAQT equation of motion is applied to the density of states to establish a unique kinetic path from an initial thermodynamic state to a stable equilibrium state. The kinetic path describes the brush's evolution in state space, as it interacts with a thermal reservoir. The predicted occupation probabilities along the kinetic path are used to determine the expected thermodynamic and structural properties. The polymer density profile of a polystyrene brush in cyclohexane solvent is predicted using the equation of motion, and it agrees qualitatively with the experimental density profiles. The Flory-Huggins parameter chosen to describe brush-solvent interactions affects the solvent distribution in the brush but has a minimal impact on the polymer density profile. Three types of nonequilibrium kinetic paths with differing amounts of entropy production are considered: a heating path, a cooling path, and a heating-cooling path. Properties such as tortuosity, radius of gyration, brush density, solvent density, and brush chain conformations are calculated for each path.

Predicting non-equilibrium folding behavior of polymer chains using the steepest-entropy-ascent quantum thermodynamic framework

The steepest-entropy-ascent quantum thermodynamic (SEAQT) framework is used to explore the influence of heating and cooling on polymer chain folding kinetics. The framework predicts how a chain moves from an initial non-equilibrium state to stable equilibrium along a unique thermodynamic path. The thermodynamic state is expressed by occupation probabilities corresponding to the levels of a discrete energy landscape. The landscape is generated using the Replica Exchange Wang-Landau method applied to a polymer chain represented by a sequence of hydrophobic and polar monomers with a simple hydrophobic-polar amino acid model. The chain conformation evolves as energy shifts among the levels of the energy landscape according to the principle of steepest entropy ascent. This principle is implemented via the SEAQT equation of motion. The SEAQT framework has the benefit of providing insight into structural properties under non-equilibrium conditions. Chain conformations during heating and cooling change continuously without sharp transitions in morphology. The changes are more drastic along non-equilibrium paths than along quasi-equilibrium paths. The SEAQT-predicted kinetics are fitted to rates associated with the experimental intensity profiles of cytochrome c protein folding with Rouse dynamics.

Predicting defect stability and annealing kinetics in two-dimensional PtSe2using steepest entropy ascent quantum thermodynamics

The steepest-entropy-ascent quantum thermodynamic (SEAQT) framework was used to calculate the stability of a collection of point defects in 2D PtSe₂and predict the kinetics with which defects rearrange during thermal annealing. The framework provides a non-equilibrium, ensemble-based framework with a self-consistent link between mechanics (both quantum and classical) and thermodynamics. It employs an equation of motion derived from the principle of steepest entropy ascent (maximum entropy production) to predict the time evolution of a set of occupation probabilities that define the states of a system undergoing a non-equilibrium process. The system is described by a degenerate energy landscape of eigenvalues, and the entropy is found from the occupation probabilities and the eigenlevel degeneracies. Scanning tunneling microscopy was used to identify the structure and distribution of point defects observed experimentally in a 2D PtSe₂film. A catalog of observed defects includes six unique point defects (vacancies and anti-site defects on Pt and Se sublattices) and twenty combinations of multiple point defects in close proximity. The defect energies were estimated with density functional theory, while the degeneracies, or density of states, for the 2D film with all possible combinations or arrangements of cataloged defects was constructed using a non-Markovian Monte-Carlo approach (i.e. the Replica-Exchange-Wang-Landau algorithm (Vogelet al2013Phys. Rev. Lett.110210603)) with a q-state Potts model. The energy landscape and associated degeneracies were determined for a 2D PtSe₂film two molecules thick and30×30unit cells in area (total of 5400 atoms). The SEAQT equation of motion was applied to the energy landscape to determine how an arbitrary density and arrangement of the six defect types evolve during annealing. Two annealing processes were modeled: heating from 77 K (-196 ∘C) to 523 K (250 ^∘C) and isothermal annealing at 523 K. The SEAQT framework predicted defect configurations, which were consistent with experimental STM images.

Kinetics of phase separation in thermally isolated critical binary fluids

Spinodal decomposition in a near-critical binary fluid is examined for experimental scenarios in which the liquid is quenched abruptly by changing the pressure and the subsequent phase separation occurs with no heat flow from the outside, i.e., adiabatically. Equations of motion for the system volume and effective temperature are derived. It is shown for this case that the nonequilibrium decomposition process is well approximated as one of constant entropy, i.e., as thermodynamically reversible. Quantitative comparison, with no adjustable parameters, is made with the experimental light scattering data of Bailey and Cannell [Phys. Rev. Lett. 70, 2110 (1993)PRLTAO0031-900710.1103/PhysRevLett.70.2110]. It is found that including these adiabatic effects accounts for most of the discrepancies between these experiments and previous isothermal theory. The equilibrium static critical properties of the isothermal theory are also examined, this discussion serving to justify some approximations in the current theory.

Low-temperature atomistic spin relaxation and non-equilibrium intensive properties using steepest-entropy-ascent quantum-inspired thermodynamics modeling

The magnetization of body-centered cubic iron at low temperatures is calculated with the steepest-entropy-ascent quantum thermodynamics (SEAQT) framework. This framework assumes that a thermodynamic property in an isolated system traces the path through state space with the greatest entropy production. Magnetization is calculated from the expected value of a thermodynamic ensemble of quantized spin waves based on the Heisenberg spin model applied to an ensemble of coupled harmonic oscillators. A realistic energy landscape is obtained from a magnon dispersion relation calculated using spin-density-functional-theory. The equilibrium magnetization as well as the evolution of magnetization from a non-equilibrium state to equilibrium are calculated from the path of steepest entropy ascent determined from the SEAQT equation of motion in state space. The framework makes it possible to model the temperature and time-dependence of magnetization without a detailed description of magnetic damping. The approach is also used to define intensive properties (temperature and magnetic field strength) that are fundamentally, i.e. canonically or grand canonically, valid for any non-equilibrium state. Given the assumed magnon dispersion relation, the SEAQT framework is used to calculate the equilibrium magnetization at different temperatures and external magnetic fields and the results are shown to closely agree with experiment for temperatures less than 500 K. The time-dependent evolution of magnetization from different initial states and interactions with a reservoir is also predicted.

Time-Energy and Time-Entropy Uncertainty Relations in Nonequilibrium Quantum Thermodynamics under Steepest-Entropy-Ascent Nonlinear Master Equations

In the domain of nondissipative unitary Hamiltonian dynamics, the well-known Mandelstam–Tamm–Messiah time–energy uncertainty relation τFΔH≥ℏ/2 provides a general lower bound to the characteristic time τF=ΔF/|d〈F〉/dt| with which the mean value of a generic quantum observable F can change with respect to the width ΔF of its uncertainty distribution (square root of F fluctuations). A useful practical consequence is that in unitary dynamics the states with longer lifetimes are those with smaller energy uncertainty ΔH (square root of energy fluctuations). Here we show that when unitary evolution is complemented with a steepest-entropy-ascent model of dissipation, the resulting nonlinear master equation entails that these lower bounds get modified and depend also on the entropy uncertainty ΔS (square root of entropy fluctuations). For example, we obtain the time–energy-and–time–entropy uncertainty relation (2τFΔH/ℏ)2+(τFΔS/kBτ)2≥1 where τ is a characteristic dissipation time functional that for each given state defines the strength of the nonunitary, steepest-entropy-ascent part of the assumed master equation. For purely dissipative dynamics this reduces to the time–entropy uncertainty relation τFΔS≥kBτ, meaning that the nonequilibrium dissipative states with longer lifetime are those with smaller entropy uncertainty ΔS.