1
|
McDonald J, von Spakovsky MR, Reynolds WT. Predicting Ion Sequestration in Charged Polymers with the Steepest-Entropy-Ascent Quantum Thermodynamic Framework. NANOMATERIALS (BASEL, SWITZERLAND) 2024; 14:458. [PMID: 38470788 DOI: 10.3390/nano14050458] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/10/2024] [Revised: 02/26/2024] [Accepted: 02/27/2024] [Indexed: 03/14/2024]
Abstract
The steepest-entropy-ascent quantum thermodynamic framework is used to investigate the effectiveness of multi-chain polyethyleneimine-methylenephosphonic acid in sequestering rare-earth ions (Eu3+) from aqueous solutions. The framework applies a thermodynamic equation of motion to a discrete energy eigenstructure to model the binding kinetics of europium ions to reactive sites of the polymer chains. The energy eigenstructure is generated using a non-Markovian Monte Carlo model that estimates energy level degeneracies. The equation of motion is used to determine the occupation probability of each energy level, describing the unique path through thermodynamic state space by which the polymer system sequesters rare-earth ions from solution. A second Monte Carlo simulation is conducted to relate the kinetic path in state space to physical descriptors associated with the polymer, including the radius of gyration, tortuosity, and Eu-neighbor distribution functions. These descriptors are used to visualize the evolution of the polymer during the sequestration process. The fraction of sequestered Eu3+ ions depends upon the total energy of the system, with lower energy resulting in greater sequestration. The kinetics of the overall sequestration are dependent on the steepest-entropy-ascent principle used by the equation of motion to generate a unique kinetic path from an initial non-equilibrium state.
Collapse
Affiliation(s)
- Jared McDonald
- Materials Science & Engineering Department, Virginia Tech, Blacksburg, VA 24061, USA
| | | | - William T Reynolds
- Materials Science & Engineering Department, Virginia Tech, Blacksburg, VA 24061, USA
| |
Collapse
|
2
|
Leamer JM, Dawson W, Bondar DI. Positivity preserving density matrix minimization at finite temperatures via square root. J Chem Phys 2024; 160:074107. [PMID: 38375902 DOI: 10.1063/5.0189864] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2023] [Accepted: 01/18/2024] [Indexed: 02/21/2024] Open
Abstract
We present a Wave Operator Minimization (WOM) method for calculating the Fermi-Dirac density matrix for electronic structure problems at finite temperature while preserving physicality by construction using the wave operator, i.e., the square root of the density matrix. WOM models cooling a state initially at infinite temperature down to the desired finite temperature. We consider both the grand canonical (constant chemical potential) and canonical (constant number of electrons) ensembles. Additionally, we show that the number of steps required for convergence is independent of the number of atoms in the system. We hope that the discussion and results presented in this article reinvigorate interest in density matrix minimization methods.
Collapse
Affiliation(s)
- Jacob M Leamer
- Department of Physics and Engineering Physics, Tulane University, 6823 St. Charles Ave., New Orleans, Louisiana 70118, USA
| | - William Dawson
- RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Denys I Bondar
- Department of Physics and Engineering Physics, Tulane University, 6823 St. Charles Ave., New Orleans, Louisiana 70118, USA
| |
Collapse
|
3
|
McDonald J, von Spakovsky MR, Reynolds WT. Predicting Polymer Brush Behavior in Solvents Using the Steepest-Entropy-Ascent Quantum Thermodynamic Framework. J Phys Chem B 2023; 127:10370-10391. [PMID: 38006350 DOI: 10.1021/acs.jpcb.3c02713] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2023]
Abstract
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.
Collapse
Affiliation(s)
- Jared McDonald
- Materials Science & Engineering Department, Virginia Tech, Blacksburg, Virginia 24061, United States
| | | | - William T Reynolds
- Materials Science & Engineering Department, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
4
|
McDonald J, von Spakovsky MR, Reynolds WT. Predicting non-equilibrium folding behavior of polymer chains using the steepest-entropy-ascent quantum thermodynamic framework. J Chem Phys 2023; 158:104904. [PMID: 36922120 DOI: 10.1063/5.0137444] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/16/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Jared McDonald
- Materials Science and Engineering Department, Virginia Tech, Blacksburg, Virginia 24061, USA
| | | | - William T Reynolds
- Materials Science and Engineering Department, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
5
|
Schuch D, Bonilla-Licea M. Uncertainty Relations in the Madelung Picture Including a Dissipative Environment. ENTROPY (BASEL, SWITZERLAND) 2023; 25:312. [PMID: 36832678 PMCID: PMC9955996 DOI: 10.3390/e25020312] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/24/2022] [Revised: 01/31/2023] [Accepted: 02/02/2023] [Indexed: 06/18/2023]
Abstract
In a recent paper, we have shown how in Madelung's hydrodynamic formulation of quantum mechanics, the uncertainties are related to the phase and amplitude of the complex wave function. Now we include a dissipative environment via a nonlinear modified Schrödinger equation. The effect of the environment is described by a complex logarithmic nonlinearity that vanishes on average. Nevertheless, there are various changes in the dynamics of the uncertainties originating from the nonlinear term. Again, this is illustrated explicitly using generalized coherent states as examples. With particular focus on the quantum mechanical contribution to the energy and the uncertainty product, connections can be made with the thermodynamic properties of the environment.
Collapse
Affiliation(s)
- Dieter Schuch
- Institute of Theoretical Physics, Goethe-University Frankfurt, Max-von-Laue-Str. 1, D-60438 Frankfurt am Main, Germany
| | | |
Collapse
|
6
|
Younis A, Baniasadi F, von Spakovsky MR, Reynolds WT. Predicting defect stability and annealing kinetics in two-dimensional PtSe 2using steepest entropy ascent quantum thermodynamics. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 35:075703. [PMID: 36395516 DOI: 10.1088/1361-648x/aca3f1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/20/2022] [Accepted: 11/17/2022] [Indexed: 06/16/2023]
Abstract
The steepest-entropy-ascent quantum thermodynamic (SEAQT) framework was used to calculate the stability of a collection of point defects in 2D PtSe2and 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 PtSe2film. 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 PtSe2film 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.
Collapse
Affiliation(s)
- Aimen Younis
- Mechanical Engineering Department, Virginia Tech, Blacksburg, VA 24061, United States of America
| | - Fazel Baniasadi
- Materials Science Engineering Department, Virginia Tech, Blacksburg, VA 24061, United States of America
| | - Michael R von Spakovsky
- Center for Energy Systems Research, Mechanical Engineering Department, Virginia Tech, Blacksburg, VA 24061, United States of America
| | - William T Reynolds
- Materials Science Engineering Department, Virginia Tech, Blacksburg, VA 24061, United States of America
| |
Collapse
|
7
|
Beretta GP. The fourth law of thermodynamics: steepest entropy ascent. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2020; 378:20190168. [PMID: 32223406 DOI: 10.1098/rsta.2019.0168] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 11/18/2019] [Indexed: 05/18/2023]
Abstract
When thermodynamics is understood as the science (or art) of constructing effective models of natural phenomena by choosing a minimal level of description capable of capturing the essential features of the physical reality of interest, the scientific community has identified a set of general rules that the model must incorporate if it aspires to be consistent with the body of known experimental evidence. Some of these rules are believed to be so general that we think of them as laws of Nature, such as the great conservation principles, whose 'greatness' derives from their generality, as masterfully explained by Feynman in one of his legendary lectures. The second law of thermodynamics is universally contemplated among the great laws of Nature. In this paper, we show that in the past four decades, an enormous body of scientific research devoted to modelling the essential features of non-equilibrium natural phenomena has converged from many different directions and frameworks towards the general recognition (albeit still expressed in different but equivalent forms and language) that another rule is also indispensable and reveals another great law of Nature that we propose to call the 'fourth law of thermodynamics'. We state it as follows: every non-equilibrium state of a system or local subsystem for which entropy is well defined must be equipped with a metric in state space with respect to which the irreversible component of its time evolution is in the direction of steepest entropy ascent compatible with the conservation constraints. To illustrate the power of the fourth law, we derive (nonlinear) extensions of Onsager reciprocity and fluctuation-dissipation relations to the far-non-equilibrium realm within the framework of the rate-controlled constrained-equilibrium approximation (also known as the quasi-equilibrium approximation). This article is part of the theme issue 'Fundamental aspects of nonequilibrium thermodynamics'.
Collapse
|
8
|
Yamada R, von Spakovsky MR, Reynolds WT. Low-temperature atomistic spin relaxation and non-equilibrium intensive properties using steepest-entropy-ascent quantum-inspired thermodynamics modeling. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2019; 31:505901. [PMID: 31470419 DOI: 10.1088/1361-648x/ab4014] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
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.
Collapse
Affiliation(s)
- Ryo Yamada
- Materials Science and Engineering Department, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, United States of America
| | | | | |
Collapse
|
9
|
Time-Energy and Time-Entropy Uncertainty Relations in Nonequilibrium Quantum Thermodynamics under Steepest-Entropy-Ascent Nonlinear Master Equations. ENTROPY 2019; 21:e21070679. [PMID: 33267393 PMCID: PMC7515176 DOI: 10.3390/e21070679] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/13/2019] [Revised: 07/07/2019] [Accepted: 07/08/2019] [Indexed: 11/29/2022]
Abstract
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.
Collapse
|
10
|
Yamada R, von Spakovsky MR, Reynolds WT. Predicting continuous and discontinuous phase decompositions using steepest-entropy-ascent quantum thermodynamics. Phys Rev E 2019; 99:052121. [PMID: 31212545 DOI: 10.1103/physreve.99.052121] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2018] [Indexed: 06/09/2023]
Abstract
The decomposition kinetics of a solid solution into separate phases are analyzed with an equation of motion initially developed to account for dissipative processes in quantum systems. This equation and the steepest-entropy-ascent quantum thermodynamic framework of which it is a part make it possible to track kinetic processes in systems at nonequilibrium while retaining the framework of classical equilibrium thermodynamics. The general equation of motion is particularized for the case of the decomposition of a binary alloy, and a solution model is used to build an approximate energy eigenstructure, or pseudoeigenstructure, for the alloy system. This equation is then solved with the pseudoeigenstructure to obtain a unique reaction path and the decomposition kinetics of the alloy. For a hypothetical solid solution with a miscibility gap at low temperatures, the conditions under which this framework predicts a continuous transformation path (spinodal decomposition) or a discontinuous one (nucleation and growth) are demonstrated.
Collapse
Affiliation(s)
- Ryo Yamada
- Materials Science and Engineering Department, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, USA
| | - Michael R von Spakovsky
- Center for Energy Systems Research, Mechanical Engineering Department, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, USA
| | - William T Reynolds
- Materials Science and Engineering Department, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, USA
| |
Collapse
|
11
|
CH₄ Adsorption Probability on GaN(0001) and (000-1) during Metalorganic Vapor Phase Epitaxy and Its Relationship to Carbon Contamination in the Films. MATERIALS 2019; 12:ma12060972. [PMID: 30909584 PMCID: PMC6470845 DOI: 10.3390/ma12060972] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/25/2019] [Revised: 03/19/2019] [Accepted: 03/21/2019] [Indexed: 11/16/2022]
Abstract
Suppression of carbon contamination in GaN films grown using metalorganic vapor phase epitaxy (MOVPE) is a crucial issue in its application to high power and high frequency electronic devices. To know how to reduce the C concentration in the films, a sequential analysis based on first principles calculations is performed. Thus, surface reconstruction and the adsorption of the CH₄ produced by the decomposition of the Ga source, Ga(CH₃)₃, and its incorporation into the GaN sub-surface layers are investigated. In this sequential analysis, the dataset of the adsorption probability of CH₄ on reconstructed surfaces is indispensable, as is the energy of the C impurity in the GaN sub-surface layers. The C adsorption probability is obtained based on steepest-entropy-ascent quantum thermodynamics (SEAQT). SEAQT is a thermodynamic ensemble-based, non-phenomenological framework that can predict the behavior of non-equilibrium processes, even those far from equilibrium. This framework is suitable especially when one studies the adsorption behavior of an impurity molecule because the conventional approach, the chemical potential control method, cannot be applied to a quantitative analysis for such a system. The proposed sequential model successfully explains the influence of the growth orientation, GaN(0001) and (000-1), on the incorporation of C into the film. This model can contribute to the suppression of the C contamination in GaN MOVPE.
Collapse
|
12
|
Yamada R, von Spakovsky MR, Reynolds WT. A method for predicting non-equilibrium thermal expansion using steepest-entropy-ascent quantum thermodynamics. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:325901. [PMID: 29964269 DOI: 10.1088/1361-648x/aad072] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Steepest-entropy-ascent quantum thermodynamics (SEAQT) is an intriguing approach that describes equilibrium and dynamic processes in a self-consistent way. To date, it has been applied primarily to gas phase systems because of the difficulty in generating the complex eigenstructures (eigenvalues and eigenfunctions) associated with solid or liquid phases. In this contribution, the SEAQT modeling is extended to the solid phase by constructing a so-called pseudo-eigenstructure, and its applicability is demonstrated by calculating the thermal expansion of metallic silver for three cases: (a) stable equilibrium, (b) along three irreversible paths from different initial non-equilibrium states to stable equilibrium, and (c) along an irreversible path between two stable equilibrium states. The SEAQT framework with an anharmonic pseudo-eigenstructure predicts reasonable values for the equilibrium thermal expansion. For the irreversible cases considered, the SEAQT approach makes it possible to predict the time-dependence of lattice relaxations from the initial state to the final equilibrium state.
Collapse
Affiliation(s)
- Ryo Yamada
- Materials Science and Engineering Department, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, United States of America
| | | | | |
Collapse
|
13
|
Militello B. Steepest entropy ascent for two-state systems with slowly varying Hamiltonians. Phys Rev E 2018; 97:052113. [PMID: 29906942 DOI: 10.1103/physreve.97.052113] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2018] [Indexed: 06/08/2023]
Abstract
The steepest entropy ascent approach is considered and applied to two-state systems. When the Hamiltonian of the system is time-dependent, the principle of maximum entropy production can still be exploited; arguments to support this fact are given. In the limit of slowly varying Hamiltonians, which allows for the adiabatic approximation for the unitary part of the dynamics, the system exhibits significant robustness to the thermalization process. Specific examples such as a spin in a rotating field and a generic two-state system undergoing an avoided crossing are considered.
Collapse
Affiliation(s)
- Benedetto Militello
- Dipartimento di Fisica e Chimica, Università degli Studi di Palermo, Via Archirafi 36, I-90123 Palermo, Italy and I.N.F.N. Sezione di Catania, Via Santa Sofia 64, I-95123 Catania, Italy
| |
Collapse
|
14
|
Fongang Achu G, Moukam Kakmeni FM, Dikande AM. Breathing pulses in the damped-soliton model for nerves. Phys Rev E 2018; 97:012211. [PMID: 29448340 DOI: 10.1103/physreve.97.012211] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2017] [Indexed: 06/08/2023]
Abstract
Unlike the Hodgkin-Huxley picture in which the nerve impulse results from ion exchanges across the cell membrane through ion-gate channels, in the so-called soliton model the impulse is seen as an electromechanical process related to thermodynamical phenomena accompanying the generation of the action potential. In this work, account is taken of the effects of damping on the nerve impulse propagation, within the framework of the soliton model. Applying the reductive perturbation expansion on the resulting KdV-Burgers equation, a damped nonlinear Schrödinger equation is derived and shown to admit breathing-type solitary wave solutions. Under specific constraints, these breathing pulse solitons become self-trapped structures in which the damping is balanced by nonlinearity such that the pulse amplitude remains unchanged even in the presence of damping.
Collapse
Affiliation(s)
- G Fongang Achu
- Complex Systems and Theoretical Biology Group (CoSTBiG), Laboratory of Research on Advanced Materials and Nonlinear Science (LaRAMaNS), Department of Physics, Faculty of Science, University of Buea, P.O. Box 63, Buea, Cameroon
| | - F M Moukam Kakmeni
- Complex Systems and Theoretical Biology Group (CoSTBiG), Laboratory of Research on Advanced Materials and Nonlinear Science (LaRAMaNS), Department of Physics, Faculty of Science, University of Buea, P.O. Box 63, Buea, Cameroon
| | - A M Dikande
- Complex Systems and Theoretical Biology Group (CoSTBiG), Laboratory of Research on Advanced Materials and Nonlinear Science (LaRAMaNS), Department of Physics, Faculty of Science, University of Buea, P.O. Box 63, Buea, Cameroon
| |
Collapse
|
15
|
Kusaba A, Li G, von Spakovsky MR, Kangawa Y, Kakimoto K. Modeling the Non-Equilibrium Process of the Chemical Adsorption of Ammonia on GaN(0001) Reconstructed Surfaces Based on Steepest-Entropy-Ascent Quantum Thermodynamics. MATERIALS (BASEL, SWITZERLAND) 2017; 10:E948. [PMID: 28809816 PMCID: PMC5578314 DOI: 10.3390/ma10080948] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/21/2017] [Revised: 08/11/2017] [Accepted: 08/11/2017] [Indexed: 11/17/2022]
Abstract
Clearly understanding elementary growth processes that depend on surface reconstruction is essential to controlling vapor-phase epitaxy more precisely. In this study, ammonia chemical adsorption on GaN(0001) reconstructed surfaces under metalorganic vapor phase epitaxy (MOVPE) conditions (3Ga-H and Nad-H + Ga-H on a 2 × 2 unit cell) is investigated using steepest-entropy-ascent quantum thermodynamics (SEAQT). SEAQT is a thermodynamic-ensemble based, first-principles framework that can predict the behavior of non-equilibrium processes, even those far from equilibrium where the state evolution is a combination of reversible and irreversible dynamics. SEAQT is an ideal choice to handle this problem on a first-principles basis since the chemical adsorption process starts from a highly non-equilibrium state. A result of the analysis shows that the probability of adsorption on 3Ga-H is significantly higher than that on Nad-H + Ga-H. Additionally, the growth temperature dependence of these adsorption probabilities and the temperature increase due to the heat of reaction is determined. The non-equilibrium thermodynamic modeling applied can lead to better control of the MOVPE process through the selection of preferable reconstructed surfaces. The modeling also demonstrates the efficacy of DFT-SEAQT coupling for determining detailed non-equilibrium process characteristics with a much smaller computational burden than would be entailed with mechanics-based, microscopic-mesoscopic approaches.
Collapse
Affiliation(s)
- Akira Kusaba
- Department of Aeronautics and Astronautics, Kyushu University, Fukuoka 819-0395, Japan.
| | - Guanchen Li
- Department of Engineering Science, University of Oxford, Parks Road, Oxford OX1 3PJ, UK.
- Center for Energy Systems Research (CESR), Mechanical Engineering Department, Virginia Tech, Blacksburg, VA 24061, USA.
| | - Michael R von Spakovsky
- Center for Energy Systems Research (CESR), Mechanical Engineering Department, Virginia Tech, Blacksburg, VA 24061, USA.
| | - Yoshihiro Kangawa
- Department of Aeronautics and Astronautics, Kyushu University, Fukuoka 819-0395, Japan.
- Research Institute for Applied Mechanics (RIAM), Kyushu University, Fukuoka 816-8580, Japan.
- Center for Integrated Research of Future Electronics (CIRFE), Institute of Materials and Systems for Sustainability (IMaSS), Nagoya University, Nagoya 464-8601, Japan.
| | - Koichi Kakimoto
- Department of Aeronautics and Astronautics, Kyushu University, Fukuoka 819-0395, Japan.
- Research Institute for Applied Mechanics (RIAM), Kyushu University, Fukuoka 816-8580, Japan.
| |
Collapse
|
16
|
Kim I, von Spakovsky MR. Ab initio relaxation times and time-dependent Hamiltonians within the steepest-entropy-ascent quantum thermodynamic framework. Phys Rev E 2017; 96:022129. [PMID: 28950609 DOI: 10.1103/physreve.96.022129] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2017] [Indexed: 11/07/2022]
Abstract
Quantum systems driven by time-dependent Hamiltonians are considered here within the framework of steepest-entropy-ascent quantum thermodynamics (SEAQT) and used to study the thermodynamic characteristics of such systems. In doing so, a generalization of the SEAQT framework valid for all such systems is provided, leading to the development of an ab initio physically relevant expression for the intrarelaxation time, an important element of this framework and one that had as of yet not been uniquely determined as an integral part of the theory. The resulting expression for the relaxation time is valid as well for time-independent Hamiltonians as a special case and makes the description provided by the SEAQT framework more robust at the fundamental level. In addition, the SEAQT framework is used to help resolve a fundamental issue of thermodynamics in the quantum domain, namely, that concerning the unique definition of process-dependent work and heat functions. The developments presented lead to the conclusion that this framework is not just an alternative approach to thermodynamics in the quantum domain but instead one that uniquely sheds new light on various fundamental but as of yet not completely resolved questions of thermodynamics.
Collapse
Affiliation(s)
- Ilki Kim
- Center for Energy Research and Technology, North Carolina A&T State University, Greensboro, North Carolina 27411, USA
| | | |
Collapse
|
17
|
Beretta GP, Al-Abbasi O, von Spakovsky MR. Steepest-entropy-ascent nonequilibrium quantum thermodynamic framework to model chemical reaction rates at an atomistic level. Phys Rev E 2017; 95:042139. [PMID: 28505826 DOI: 10.1103/physreve.95.042139] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2016] [Indexed: 11/07/2022]
Abstract
The steepest entropy ascent (SEA) dynamical principle provides a general framework for modeling the dynamics of nonequilibrium (NE) phenomena at any level of description, including the atomistic one. It has recently been shown to provide a precise implementation and meaning to the maximum entropy production principle and to encompass many well-established theories of nonequilibrium thermodynamics into a single unifying geometrical framework. Its original formulation in the framework of quantum thermodynamics (QT) assumes the simplest and most natural Fisher-Rao metric to geometrize from a dynamical standpoint the manifold of density operators, which represent the thermodynamic NE states of the system. This simplest SEAQT formulation is used here to develop a general mathematical framework for modeling the NE time evolution of the quantum state of a chemically reactive mixture at an atomistic level. The method is illustrated for a simple two-reaction kinetic scheme of the overall reaction F+H_{2}⇔HF+F in an isolated tank of fixed volume. However, the general formalism is developed for a reactive system subject to multiple reaction mechanisms. To explicitly implement the SEAQT nonlinear law of evolution for the density operator, both the energy and the particle number eigenvalue problems are set up and solved analytically under the dilute gas approximation. The system-level energy and particle number eigenvalues and eigenstates are used in the SEAQT equation of motion to determine the time evolution of the density operator, thus effectively describing the overall kinetics of the reacting system as it relaxes toward stable chemical equilibrium. The predicted time evolution in the near-equilibrium limit is compared to the reaction rates given by a standard detailed kinetic model so as to extract the single time constant needed by the present SEA model.
Collapse
Affiliation(s)
- G P Beretta
- Mechanical and Industrial Engineering Department, Università di Brescia, via Branze 38, 25123 Brescia, Italy
| | - Omar Al-Abbasi
- Center for Energy Systems Research, Mechanical Engineering Department, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - M R von Spakovsky
- Center for Energy Systems Research, Mechanical Engineering Department, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
18
|
Li G, von Spakovsky MR. Generalized thermodynamic relations for a system experiencing heat and mass diffusion in the far-from-equilibrium realm based on steepest entropy ascent. Phys Rev E 2016; 94:032117. [PMID: 27739710 DOI: 10.1103/physreve.94.032117] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2016] [Indexed: 06/06/2023]
Abstract
This paper presents a nonequilibrium thermodynamic model for the relaxation of a local, isolated system in nonequilibrium using the principle of steepest entropy ascent (SEA), which can be expressed as a variational principle in thermodynamic state space. The model is able to arrive at the Onsager relations for such a system. Since no assumption of local equilibrium is made, the conjugate fluxes and forces are intrinsic to the subspaces of the system's state space and are defined using the concepts of hypoequilibrium state and nonequilibrium intensive properties, which describe the nonmutual equilibrium status between subspaces of the thermodynamic state space. The Onsager relations are shown to be a thermodynamic kinematic feature of the system independent of the specific details of the micromechanical dynamics. Two kinds of relaxation processes are studied with different constraints (i.e., conservation laws) corresponding to heat and mass diffusion. Linear behavior in the near-equilibrium region as well as nonlinear behavior in the far-from-equilibrium region are discussed. Thermodynamic relations in the equilibrium and near-equilibrium realm, including the Gibbs relation, the Clausius inequality, and the Onsager relations, are generalized to the far-from-equilibrium realm. The variational principle in the space spanned by the intrinsic conjugate fluxes and forces is expressed via the quadratic dissipation potential. As an application, the model is applied to the heat and mass diffusion of a system represented by a single-particle ensemble, which can also be applied to a simple system of many particles. Phenomenological transport coefficients are also derived in the near-equilibrium realm.
Collapse
Affiliation(s)
- Guanchen Li
- Center for Energy Systems Research, Mechanical Engineering Department, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Michael R von Spakovsky
- Center for Energy Systems Research, Mechanical Engineering Department, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
19
|
Comparing the Models of Steepest Entropy Ascent Quantum Thermodynamics, Master Equation and the Difference Equation for a Simple Quantum System Interacting with Reservoirs. ENTROPY 2016. [DOI: 10.3390/e18050176] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
20
|
Li G, von Spakovsky MR. Steepest-entropy-ascent quantum thermodynamic modeling of the relaxation process of isolated chemically reactive systems using density of states and the concept of hypoequilibrium state. Phys Rev E 2016; 93:012137. [PMID: 26871054 DOI: 10.1103/physreve.93.012137] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2015] [Indexed: 06/05/2023]
Abstract
This paper presents a study of the nonequilibrium relaxation process of chemically reactive systems using steepest-entropy-ascent quantum thermodynamics (SEAQT). The trajectory of the chemical reaction, i.e., the accessible intermediate states, is predicted and discussed. The prediction is made using a thermodynamic-ensemble approach, which does not require detailed information about the particle mechanics involved (e.g., the collision of particles). Instead, modeling the kinetics and dynamics of the relaxation process is based on the principle of steepest-entropy ascent (SEA) or maximum-entropy production, which suggests a constrained gradient dynamics in state space. The SEAQT framework is based on general definitions for energy and entropy and at least theoretically enables the prediction of the nonequilibrium relaxation of system state at all temporal and spatial scales. However, to make this not just theoretically but computationally possible, the concept of density of states is introduced to simplify the application of the relaxation model, which in effect extends the application of the SEAQT framework even to infinite energy eigenlevel systems. The energy eigenstructure of the reactive system considered here consists of an extremely large number of such levels (on the order of 10^{130}) and yields to the quasicontinuous assumption. The principle of SEA results in a unique trajectory of system thermodynamic state evolution in Hilbert space in the nonequilibrium realm, even far from equilibrium. To describe this trajectory, the concepts of subsystem hypoequilibrium state and temperature are introduced and used to characterize each system-level, nonequilibrium state. This definition of temperature is fundamental rather than phenomenological and is a generalization of the temperature defined at stable equilibrium. In addition, to deal with the large number of energy eigenlevels, the equation of motion is formulated on the basis of the density of states and a set of associated degeneracies. Their significance for the nonequilibrium evolution of system state is discussed. For the application presented, the numerical method used is described and is based on the density of states, which is specifically developed to solve the SEAQT equation of motion. Results for different kinds of initial nonequilibrium conditions, i.e., those for gamma and Maxwellian distributions, are studied. The advantage of the concept of hypoequilibrium state in studying nonequilibrium trajectories is discussed.
Collapse
Affiliation(s)
- Guanchen Li
- Center for Energy Systems Research, Mechanical Engineering Department, Virginia Tech, Blacksburg, Virginia 24061
| | - Michael R von Spakovsky
- Center for Energy Systems Research, Mechanical Engineering Department, Virginia Tech, Blacksburg, Virginia 24061
| |
Collapse
|
21
|
Reina C, Zimmer J. Entropy production and the geometry of dissipative evolution equations. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:052117. [PMID: 26651657 DOI: 10.1103/physreve.92.052117] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/14/2015] [Indexed: 06/05/2023]
Abstract
Purely dissipative evolution equations are often cast as gradient flow structures, z ̇=K(z)DS(z), where the variable z of interest evolves towards the maximum of a functional S according to a metric defined by an operator K. While the functional often follows immediately from physical considerations (e.g., the thermodynamic entropy), the operator K and the associated geometry does not necessarily do so (e.g., Wasserstein geometry for diffusion). In this paper, we present a variational statement in the sense of maximum entropy production that directly delivers a relationship between the operator K and the constraints of the system. In particular, the Wasserstein metric naturally arises here from the conservation of mass or energy, and depends on the Onsager resistivity tensor, which, itself, may be understood as another metric, as in the steepest entropy ascent formalism. This variational principle is exemplified here for the simultaneous evolution of conserved and nonconserved quantities in open systems. It thus extends the classical Onsager flux-force relationships and the associated variational statement to variables that do not have a flux associated to them. We further show that the metric structure K is intimately linked to the celebrated Freidlin-Wentzell theory of stochastically perturbed gradient flows, and that the proposed variational principle encloses an infinite-dimensional fluctuation-dissipation statement.
Collapse
Affiliation(s)
- Celia Reina
- Department of Mechanical Engineering and Applied Mechanics, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Johannes Zimmer
- Department of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, United Kingdom
| |
Collapse
|
22
|
Montefusco A, Consonni F, Beretta GP. Essential equivalence of the general equation for the nonequilibrium reversible-irreversible coupling (GENERIC) and steepest-entropy-ascent models of dissipation for nonequilibrium thermodynamics. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:042138. [PMID: 25974469 DOI: 10.1103/physreve.91.042138] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/19/2014] [Indexed: 06/04/2023]
Abstract
By reformulating the steepest-entropy-ascent (SEA) dynamical model for nonequilibrium thermodynamics in the mathematical language of differential geometry, we compare it with the primitive formulation of the general equation for the nonequilibrium reversible-irreversible coupling (GENERIC) model and discuss the main technical differences of the two approaches. In both dynamical models the description of dissipation is of the "entropy-gradient" type. SEA focuses only on the dissipative, i.e., entropy generating, component of the time evolution, chooses a sub-Riemannian metric tensor as dissipative structure, and uses the local entropy density field as potential. GENERIC emphasizes the coupling between the dissipative and nondissipative components of the time evolution, chooses two compatible degenerate structures (Poisson and degenerate co-Riemannian), and uses the global energy and entropy functionals as potentials. As an illustration, we rewrite the known GENERIC formulation of the Boltzmann equation in terms of the square root of the distribution function adopted by the SEA formulation. We then provide a formal proof that in more general frameworks, whenever all degeneracies in the GENERIC framework are related to conservation laws, the SEA and GENERIC models of the dissipative component of the dynamics are essentially interchangeable, provided of course they assume the same kinematics. As part of the discussion, we note that equipping the dissipative structure of GENERIC with the Leibniz identity makes it automatically SEA on metric leaves.
Collapse
|
23
|
Li G, Al-Abbasi O, von Spakovsky MR. Atomistic-level non-equilibrium model for chemically reactive systems based on steepest-entropy-ascent quantum thermodynamics. ACTA ACUST UNITED AC 2014. [DOI: 10.1088/1742-6596/538/1/012013] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
24
|
|
25
|
Smith CE, Spakovsky MRV. Comparison of the non-equilibrium predictions of Intrinsic Quantum Thermodynamics at the atomistic level with experimental evidence. ACTA ACUST UNITED AC 2012. [DOI: 10.1088/1742-6596/380/1/012015] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
26
|
The Rate-Controlled Constrained-Equilibrium Approach to Far-From-Local-Equilibrium Thermodynamics. ENTROPY 2012. [DOI: 10.3390/e14020092] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
27
|
Niven RK. Steady state of a dissipative flow-controlled system and the maximum entropy production principle. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 80:021113. [PMID: 19792083 DOI: 10.1103/physreve.80.021113] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/10/2009] [Indexed: 05/28/2023]
Abstract
A theory to predict the steady-state position of a dissipative flow-controlled system, as defined by a control volume, is developed based on the maximum entropy principle of Jaynes, involving minimization of a generalized free-energy-like potential. The analysis provides a theoretical justification of a local, conditional form of the maximum entropy production principle, which successfully predicts the observable properties of many such systems. The analysis reveals a very different manifestation of the second law of thermodynamics in steady-state flow systems, which provides a driving force for the formation of complex systems, including life.
Collapse
Affiliation(s)
- Robert K Niven
- School of Engineering and Information Technology, The University of New South Wales at ADFA, Canberra, Australian Capital Territory 2600, Australia.
| |
Collapse
|
28
|
Narnhofer H, Posch HA, Thirring W. Emergence of order in quantum extensions of the classical quasispecies evolution. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2007; 76:041133. [PMID: 17994962 DOI: 10.1103/physreve.76.041133] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/12/2007] [Indexed: 05/25/2023]
Abstract
We study evolution equations which model selection and mutation within the framework of quantum mechanics. The main question is to what extent order is achieved for an ensemble of typical systems. As an indicator for mixing or purification, a quadratic entropy is used which assumes values between zero for pure states and (d-1)/d for fully mixed states. Here, d is the dimension. Whereas the classical counterpart, the quasispecies dynamics, has previously been found to be predominantly mixing, the quantum quasispecies (QS) evolution surprisingly is found to be strictly purifying for all dimensions. This is also typically true for an alternative formulation (AQS) of this quantum mechanical flow. We compare this also to analogous results for the Lindblad evolution. Although the latter may be viewed as a simple linear superposition of the purifying QS and AQS evolutions, it is found to be predominantly mixing. The reason for this behavior may be explained by the fact that the two subprocesses by themselves converge to different pure states, such that the combined process is mixing. These results also apply to high-dimensional systems.
Collapse
Affiliation(s)
- Heide Narnhofer
- Faculty of Physics, Universität Wien, Boltzmanngasse 5, A-1090 Wien, Austria.
| | | | | |
Collapse
|