Pub Date : 2024-06-21DOI: 10.1007/s10955-024-03293-w
Miroslav Grmela
Approach of mesoscopic state variables to time independent equilibrium sates (zero law of thermodynamics) gives birth to the classical equilibrium thermodynamics. Approach of fluxes and forces to fixed points (equilibrium fluxes and forces) that drive reduced mesoscopic dynamics gives birth to the rate thermodynamics that is applicable to driven systems. We formulate the rate thermodynamics and dynamics, investigate its relation to the classical thermodynamics, to extensions involving more details, to the hierarchy reformulations of dynamical theories, and to the Onsager variational principle. We also compare thermodynamic and dynamic critical behavior observed in closed and open systems. Dynamics and thermodynamics of the van der Waals gas provides an illustration.
{"title":"Thermodynamics and Rate Thermodynamics","authors":"Miroslav Grmela","doi":"10.1007/s10955-024-03293-w","DOIUrl":"https://doi.org/10.1007/s10955-024-03293-w","url":null,"abstract":"<p>Approach of mesoscopic state variables to time independent equilibrium sates (zero law of thermodynamics) gives birth to the classical equilibrium thermodynamics. Approach of fluxes and forces to fixed points (equilibrium fluxes and forces) that drive reduced mesoscopic dynamics gives birth to the rate thermodynamics that is applicable to driven systems. We formulate the rate thermodynamics and dynamics, investigate its relation to the classical thermodynamics, to extensions involving more details, to the hierarchy reformulations of dynamical theories, and to the Onsager variational principle. We also compare thermodynamic and dynamic critical behavior observed in closed and open systems. Dynamics and thermodynamics of the van der Waals gas provides an illustration.</p>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":null,"pages":null},"PeriodicalIF":1.6,"publicationDate":"2024-06-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141514561","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-06-21DOI: 10.1007/s10955-024-03290-z
Alexander D. Kolesnik
We consider the symmetric Markov random flight, also called the persistent random walk, performed by a particle that moves at constant finite speed in the Euclidean space (mathbb {R}^m, ; mge 2,) and changes its direction at Poisson-distributed time instants by taking it at random according to the uniform distribution on the surface of the unit ((m-1))-dimensional sphere. Such stochastic motion has become a very popular object of modern statistical physics because it can serve as an appropriate model for describing the isotropic finite-velocity transport in multidimensional Euclidean spaces. In recent decade this approach was also developed in the framework of the run-and-tumble theory. In this article we study one of the most important characteristics of the multidimensional symmetric Markov random flight, namely, its characteristic function. We derive two series representations of the characteristic function of the process with respect to Bessel functions with variable indices and with respect to the powers of time variable. The coefficients of these series are given by recurrent relations, as well as in the form of special determinants. As an application of these results, an asymptotic formula for the second moment function (mu _{(2,2,2)}(t), ; t>0,) of the three-dimensional Markov random flight, is presented. The moment function (mu _{(2,0,0)}(t), ; t>0,) is obtained in an explicit form.
{"title":"Series Representations for the Characteristic Function of the Multidimensional Markov Random Flight","authors":"Alexander D. Kolesnik","doi":"10.1007/s10955-024-03290-z","DOIUrl":"https://doi.org/10.1007/s10955-024-03290-z","url":null,"abstract":"<p>We consider the symmetric Markov random flight, also called the persistent random walk, performed by a particle that moves at constant finite speed in the Euclidean space <span>(mathbb {R}^m, ; mge 2,)</span> and changes its direction at Poisson-distributed time instants by taking it at random according to the uniform distribution on the surface of the unit <span>((m-1))</span>-dimensional sphere. Such stochastic motion has become a very popular object of modern statistical physics because it can serve as an appropriate model for describing the isotropic finite-velocity transport in multidimensional Euclidean spaces. In recent decade this approach was also developed in the framework of the run-and-tumble theory. In this article we study one of the most important characteristics of the multidimensional symmetric Markov random flight, namely, its characteristic function. We derive two series representations of the characteristic function of the process with respect to Bessel functions with variable indices and with respect to the powers of time variable. The coefficients of these series are given by recurrent relations, as well as in the form of special determinants. As an application of these results, an asymptotic formula for the second moment function <span>(mu _{(2,2,2)}(t), ; t>0,)</span> of the three-dimensional Markov random flight, is presented. The moment function <span>(mu _{(2,0,0)}(t), ; t>0,)</span> is obtained in an explicit form.</p>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":null,"pages":null},"PeriodicalIF":1.6,"publicationDate":"2024-06-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141502807","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-06-20DOI: 10.1007/s10955-024-03291-y
Zhimeng Ouyang
The justification of hydrodynamic limits in non-convex domains has long been an open problem due to the singularity at the grazing set. In this paper, we investigate the unsteady neutron transport equation in a general bounded domain with the in-flow, diffuse-reflection, or specular-reflection boundary condition. Using a novel kernel estimate, we demonstrate the optimal (L^2) diffusive limit in the presence of both initial and boundary layers. Previously, this result was only proved for convex domains when the time variable is involved. Our approach is highly robust, making it applicable to all basic types of physical boundary conditions.
{"title":"Diffusive Limit of the Unsteady Neutron Transport Equation in Bounded Domains","authors":"Zhimeng Ouyang","doi":"10.1007/s10955-024-03291-y","DOIUrl":"https://doi.org/10.1007/s10955-024-03291-y","url":null,"abstract":"<p>The justification of hydrodynamic limits in non-convex domains has long been an open problem due to the singularity at the grazing set. In this paper, we investigate the unsteady neutron transport equation in a general bounded domain with the in-flow, diffuse-reflection, or specular-reflection boundary condition. Using a novel kernel estimate, we demonstrate the optimal <span>(L^2)</span> diffusive limit in the presence of both initial and boundary layers. Previously, this result was only proved for convex domains when the time variable is involved. Our approach is highly robust, making it applicable to all basic types of physical boundary conditions.</p>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":null,"pages":null},"PeriodicalIF":1.6,"publicationDate":"2024-06-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141502808","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-06-18DOI: 10.1007/s10955-024-03279-8
Caio Alves, Rodrigo Ribeiro, Rémy Sanchis
In this paper, we investigate the global clustering coefficient (a.k.a transitivity) and clique number of graphs generated by a preferential attachment random graph model with an additional feature of allowing edge connections between existing vertices. Specifically, at each time step t, either a new vertex is added with probability f(t), or an edge is added between two existing vertices with probability (1-f(t)). We establish concentration inequalities for the global clustering and clique number of the resulting graphs under the assumption that f(t) is a regularly varying function at infinity with index of regular variation (-gamma ), where (gamma in [0,1)). We also demonstrate an inverse relation between these two statistics: the clique number is essentially the reciprocal of the global clustering coefficient.
{"title":"Clustering and Cliques in Preferential Attachment Random Graphs with Edge Insertion","authors":"Caio Alves, Rodrigo Ribeiro, Rémy Sanchis","doi":"10.1007/s10955-024-03279-8","DOIUrl":"https://doi.org/10.1007/s10955-024-03279-8","url":null,"abstract":"<p>In this paper, we investigate the global clustering coefficient (a.k.a transitivity) and clique number of graphs generated by a preferential attachment random graph model with an additional feature of allowing edge connections between existing vertices. Specifically, at each time step <i>t</i>, either a new vertex is added with probability <i>f</i>(<i>t</i>), or an edge is added between two existing vertices with probability <span>(1-f(t))</span>. We establish concentration inequalities for the global clustering and clique number of the resulting graphs under the assumption that <i>f</i>(<i>t</i>) is a regularly varying function at infinity with index of regular variation <span>(-gamma )</span>, where <span>(gamma in [0,1))</span>. We also demonstrate an inverse relation between these two statistics: the clique number is essentially the reciprocal of the global clustering coefficient.</p>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":null,"pages":null},"PeriodicalIF":1.6,"publicationDate":"2024-06-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141502809","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-06-14DOI: 10.1007/s10955-024-03287-8
Tiecheng Xu
For (alpha ge 1), let (g:{mathbb {N}}rightarrow {mathbb {R}}_+) be given by (g(0)=0), (g(1)=1), (g(k)=(k/k-1)^alpha ), (kge 2). Consider the homogeneous zero range process on a discrete set in which a particle jumps from a site x, occupied by k particles, to site y with rate (g(k)p(y-x)) for some fixed probability (p:{mathbb {Z}}rightarrow [0,1]). Armendáriz and Loulakis (Probab Theory Relat Fields 145:175–188, 2009, https://doi.org/10.1007/s00440-008-0165-7) proved a strong form of the equivalence of ensembles for the invariant measure of the supercritical zero range process with (alpha >2). We generalize their result to all (alpha ge 1).
For (α ge 1), let (g:{mathbb {N}}rightarrow {mathbb {R}}_+) be given by (g(0)=0), (g(1)=1), (g(k)=(k/k-1)^α ), (kge 2).考虑离散集合上的同质零范围过程,在这个过程中,一个粒子以某种固定概率(p:{/mathbb {Z}}rightarrow [0,1])从一个被k个粒子占据的位置x跳到位置y,速率为(g(k)p(y-x))。Armendáriz和Loulakis(Probab Theory Relat Fields 145:175-188,2009,https://doi.org/10.1007/s00440-008-0165-7)为具有(alpha >2)的超临界零范围过程的不变度量证明了集合等价的强形式。我们将他们的结果推广到所有的(alpha ge 1 )。
{"title":"Condensation of the Invariant Measures of the Supercritical Zero Range Processes","authors":"Tiecheng Xu","doi":"10.1007/s10955-024-03287-8","DOIUrl":"https://doi.org/10.1007/s10955-024-03287-8","url":null,"abstract":"<p>For <span>(alpha ge 1)</span>, let <span>(g:{mathbb {N}}rightarrow {mathbb {R}}_+)</span> be given by <span>(g(0)=0)</span>, <span>(g(1)=1)</span>, <span>(g(k)=(k/k-1)^alpha )</span>, <span>(kge 2)</span>. Consider the homogeneous zero range process on a discrete set in which a particle jumps from a site <i>x</i>, occupied by <i>k</i> particles, to site <i>y</i> with rate <span>(g(k)p(y-x))</span> for some fixed probability <span>(p:{mathbb {Z}}rightarrow [0,1])</span>. Armendáriz and Loulakis (Probab Theory Relat Fields 145:175–188, 2009, https://doi.org/10.1007/s00440-008-0165-7) proved a strong form of the equivalence of ensembles for the invariant measure of the supercritical zero range process with <span>(alpha >2)</span>. We generalize their result to all <span>(alpha ge 1)</span>.</p>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":null,"pages":null},"PeriodicalIF":1.6,"publicationDate":"2024-06-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141514562","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-06-04DOI: 10.1007/s10955-024-03267-y
Kohei Fukai
We rigorously prove that the local conserved quantities in the one-dimensional Hubbard model are uniquely determined for each locality up to the freedom to add lower-order ones. From this, we can conclude that the local conserved quantities are exhausted by those obtained from the expansion of the transfer matrix.
{"title":"Proof of Completeness of the Local Conserved Quantities in the One-Dimensional Hubbard Model","authors":"Kohei Fukai","doi":"10.1007/s10955-024-03267-y","DOIUrl":"https://doi.org/10.1007/s10955-024-03267-y","url":null,"abstract":"<p>We rigorously prove that the local conserved quantities in the one-dimensional Hubbard model are uniquely determined for each locality up to the freedom to add lower-order ones. From this, we can conclude that the local conserved quantities are exhausted by those obtained from the expansion of the transfer matrix.</p>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":null,"pages":null},"PeriodicalIF":1.6,"publicationDate":"2024-06-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141253531","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-05-30DOI: 10.1007/s10955-024-03281-0
Akira Onuki
We present a statistical mechanical theory of multi-component fluids, where we consider the correlation functions of the number densities and the energy density in the grand canonical ensemble. In terms of their space integrals we express the partial volumes ({{bar{v}}}_i), the partial enthalpies ({{bar{H}}}_i), and other thermodynamic derivatives. These ({{bar{v}}}_i) and ({{bar{H}}}_i) assume simple forms for binary mixtures and for ternary mixtures with a dilute solute. They are then related to the space-dependent thermal fluctuations of the temperature and the pressure. The space averages of these fluctuations are those introduced by Landau and Lifshits in the isothermal-isobaric (T-p) ensemble. We also give expressions for the long-range (nonlocal) correlations in the canonical and T-p ensembles, which are inversely proportional to the system volume. For a mixture solvent, we examine the solvent-induced solute–solute attraction and the osmotic enthalpy changes due to the solute doping using the correlation function integrals.
{"title":"Extension of Kirkwood–Buff Theory: Partial Enthalpies, Fluctuations of Energy Density, Temperature, and Pressure, and Solute-Induced Effects in a Mixture Solvent","authors":"Akira Onuki","doi":"10.1007/s10955-024-03281-0","DOIUrl":"https://doi.org/10.1007/s10955-024-03281-0","url":null,"abstract":"<p>We present a statistical mechanical theory of multi-component fluids, where we consider the correlation functions of the number densities and the energy density in the grand canonical ensemble. In terms of their space integrals we express the partial volumes <span>({{bar{v}}}_i)</span>, the partial enthalpies <span>({{bar{H}}}_i)</span>, and other thermodynamic derivatives. These <span>({{bar{v}}}_i)</span> and <span>({{bar{H}}}_i)</span> assume simple forms for binary mixtures and for ternary mixtures with a dilute solute. They are then related to the space-dependent thermal fluctuations of the temperature and the pressure. The space averages of these fluctuations are those introduced by Landau and Lifshits in the isothermal-isobaric (<i>T</i>-<i>p</i>) ensemble. We also give expressions for the long-range (nonlocal) correlations in the canonical and <i>T</i>-<i>p</i> ensembles, which are inversely proportional to the system volume. For a mixture solvent, we examine the solvent-induced solute–solute attraction and the osmotic enthalpy changes due to the solute doping using the correlation function integrals.</p>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":null,"pages":null},"PeriodicalIF":1.6,"publicationDate":"2024-05-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141192461","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-05-30DOI: 10.1007/s10955-024-03283-y
Stéphane Brull, Doheon Kim, Myeong-Su Lee, Seok-Bae Yun
Unlike the single species gases, the transport coefficients such as Fick, Soret, Dufour coefficients arise in the hydrodynamic limit of multi-species gas mixtures. To the best of the authors’ knowledge, no multi-component relaxational models is reported that produces all these values correctly. In this paper, we establish the existence of unique stationary mild solutions to the BGK models for gas mixtures which produces the correct Fick coefficients in the Navier–Stokes limit for inert gases (Brull in Eur J Mech B 33:74–86, 2012), and for reactive gases (Brull and Schneider in Commun Math Sci 12(7):1199–1223, 2014) in a unified manner.
与单种气体不同,菲克系数、索雷特系数和杜富尔系数等传输系数产生于多种类气体混合物的流体力学极限。据作者所知,目前还没有任何多组分弛豫模型能正确地产生所有这些数值。在本文中,我们为气体混合物的 BGK 模型建立了唯一的静态温和解,这些解在纳维-斯托克斯极限中为惰性气体(Brull 发表于 Eur J Mech B 33:74-86, 2012)和反应性气体(Brull 和 Schneider 发表于 Commun Math Sci 12(7):1199-1223, 2014)统一生成了正确的菲克系数。
{"title":"Stationary Mixture BGK Models with the Correct Fick Coefficients","authors":"Stéphane Brull, Doheon Kim, Myeong-Su Lee, Seok-Bae Yun","doi":"10.1007/s10955-024-03283-y","DOIUrl":"https://doi.org/10.1007/s10955-024-03283-y","url":null,"abstract":"<p>Unlike the single species gases, the transport coefficients such as Fick, Soret, Dufour coefficients arise in the hydrodynamic limit of multi-species gas mixtures. To the best of the authors’ knowledge, no multi-component relaxational models is reported that produces all these values correctly. In this paper, we establish the existence of unique stationary mild solutions to the BGK models for gas mixtures which produces the correct Fick coefficients in the Navier–Stokes limit for inert gases (Brull in Eur J Mech B 33:74–86, 2012), and for reactive gases (Brull and Schneider in Commun Math Sci 12(7):1199–1223, 2014) in a unified manner.</p>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":null,"pages":null},"PeriodicalIF":1.6,"publicationDate":"2024-05-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141192405","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-05-30DOI: 10.1007/s10955-024-03284-x
Henryk Gzyl
In this work, we develop a mathematical framework to model a quantum system whose time evolution may depend on the state of a randomly changing environment that evolves according to a Markovian process. When the environment changes its state, three possible things may occur: the quantum system starts evolving according to a new Hamiltonian, it may suffer an instantaneous perturbation that changes its state or both things may happen simultaneously. We consider the case of quantum systems with finite dimensional Hilbert state space, in which case the observables are described by Hermitian matrices. We show how to average over the environment to predict the expected value of the density matrix with which one can compute the expected values of the observables of interest.
{"title":"Quantum Random Evolutions","authors":"Henryk Gzyl","doi":"10.1007/s10955-024-03284-x","DOIUrl":"https://doi.org/10.1007/s10955-024-03284-x","url":null,"abstract":"<p>In this work, we develop a mathematical framework to model a quantum system whose time evolution may depend on the state of a randomly changing environment that evolves according to a Markovian process. When the environment changes its state, three possible things may occur: the quantum system starts evolving according to a new Hamiltonian, it may suffer an instantaneous perturbation that changes its state or both things may happen simultaneously. We consider the case of quantum systems with finite dimensional Hilbert state space, in which case the observables are described by Hermitian matrices. We show how to average over the environment to predict the expected value of the density matrix with which one can compute the expected values of the observables of interest.</p>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":null,"pages":null},"PeriodicalIF":1.6,"publicationDate":"2024-05-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141192406","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-05-30DOI: 10.1007/s10955-024-03282-z
Sahil Kumar Singh, Abhishek Dhar, Herbert Spohn, Anupam Kundu
We consider the relaxation of an initial non-equilibrium state in a one-dimensional fluid of hard rods. Since it is an interacting integrable system, we expect it to reach the Generalized Gibbs Ensemble (GGE) at long times for generic initial conditions. Here we show that there exist initial conditions for which the system does not reach GGE even at very long times and in the thermodynamic limit. In particular, we consider an initial condition of uniformly distributed hard-rods in a box with the left half having particles with a singular velocity distribution (all moving with unit velocity) and the right half particles in thermal equilibrium. We find that the density profile for the singular component does not spread to the full extent of the box and keeps moving with a fixed effective speed at long times. We show that such density profiles can be well described by the solution of the Euler equations almost everywhere except at the location of the shocks, where we observe slight discrepancies due to dissipation arising from the initial fluctuations of the thermal background. To demonstrate this effect of dissipation analytically, we consider a second initial condition with a single particle at the origin with unit velocity in a thermal background. We find that the probability distribution of the position of the unit velocity quasi-particle has diffusive spreading which can be understood from the solution of the Navier–Stokes (NS) equation of the hard rods. Finally, we consider an initial condition with a spread in velocity distribution for which we show convergence to GGE. Our conclusions are based on molecular dynamics simulations supported by analytical arguments.
{"title":"Thermalization and Hydrodynamics in an Interacting Integrable System: The Case of Hard Rods","authors":"Sahil Kumar Singh, Abhishek Dhar, Herbert Spohn, Anupam Kundu","doi":"10.1007/s10955-024-03282-z","DOIUrl":"https://doi.org/10.1007/s10955-024-03282-z","url":null,"abstract":"<p>We consider the relaxation of an initial non-equilibrium state in a one-dimensional fluid of hard rods. Since it is an interacting integrable system, we expect it to reach the Generalized Gibbs Ensemble (GGE) at long times for generic initial conditions. Here we show that there exist initial conditions for which the system does not reach GGE even at very long times and in the thermodynamic limit. In particular, we consider an initial condition of uniformly distributed hard-rods in a box with the left half having particles with a singular velocity distribution (all moving with unit velocity) and the right half particles in thermal equilibrium. We find that the density profile for the singular component does not spread to the full extent of the box and keeps moving with a fixed effective speed at long times. We show that such density profiles can be well described by the solution of the Euler equations almost everywhere except at the location of the shocks, where we observe slight discrepancies due to dissipation arising from the initial fluctuations of the thermal background. To demonstrate this effect of dissipation analytically, we consider a second initial condition with a single particle at the origin with unit velocity in a thermal background. We find that the probability distribution of the position of the unit velocity quasi-particle has diffusive spreading which can be understood from the solution of the Navier–Stokes (NS) equation of the hard rods. Finally, we consider an initial condition with a spread in velocity distribution for which we show convergence to GGE. Our conclusions are based on molecular dynamics simulations supported by analytical arguments.</p>","PeriodicalId":667,"journal":{"name":"Journal of Statistical Physics","volume":null,"pages":null},"PeriodicalIF":1.6,"publicationDate":"2024-05-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141192355","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}