arxiv: v1 [cond-mat.stat-mech] 5 Oct 2018

Transcription

1 On the canonical distributions of a thermal particle in the weakly confining potential of special type Tatsuaki Wada Region of Electrical and Electronic Systems Engineering, Ibaraki Uniersity, Nakanarusawa-cho, Hitachi , Japan Antonio M. Scarfone Istituto dei Sistemi Complessi, Consiglio Nazionale delle Ricerche (ISC-CNR), c/o Politecnico di Torino, 1019 Torino, Italy arxi: [cond-mat.stat-mech] 5 Oct 018 Hiroshi Matsuzoe Department of Computer Science and Engineering, Nagoya Institute of Technology, Gokiso-cho, Showa-ku, Nagoya, , Japan (Dated: October 8, 018) We consider a thermal particle in a weakly confining potential characterized by the inerse hyperbolic sine function arsinh of the particle elocity and the control parameter c. The stationary state of the Fokker-Planck equation is shown to be a canonical probability distribution. Furthermore an appropriate re-parametrization relates this stationary state with the κ-deformed Gaussian. I. INTRODUCTION When one describes basic equations for physical systems, some linear constitutie relations are often used. For example, Hooke s law F e x for the elastic force F e of a spring against the change in its length x, Stokes law F S f for the frictional force F S exerted on an object in a fluid with the flow elocity f, and so on. Howeer, as a real spring deiates from Hooke s law, we know that any linear constitutie relation describes an idealized situation, and it is merely approximated, or linearized one for describing real physical properties. In other words, non-linearity plays a crucial role to describe more realistic physical systems. In the standard linear Fokker-Plank (FP) equation, which is one of the most fundamental equations in statistical physics, a parabolic potential / is such an approximated potential function in the elocity space for describing an idealized situation of a thermal system. ThedriftforceF d ()asthegradientforceoftheharmonic potential obeys a linear constitutie relation F d (). It is well known that for a thermal particle diffusing in a parabolic potential, the steady state solution of the corresponding linear FP equation is a Gaussian probability distribution function p() exp( β /). Since any particle with high thermal energy can be captured in a region around the minimum of this parobolic potential /, we call it a strongly confining potential. Instead of such a strongly confining potential, in this contribution we focus on a thermally diffusing particle in a weakly confining potential defined by U(; c ) := c arsinh tatsuaki.wada.to@c.ibaraki.ac.jp ( c ), (1) with a control parameter c. In the limit of c, this potential reduces to a parabolic potential /. Accordingly U(; c ) can be considered as a deformation of the parabolic potential /, or as a deformation of the kinetic energy of a thermal particle with unit mass m = 1. For a finite alue of c this potential is weak (or shallow) in order to capture a thermal particle in a region, say R c, around the minimum point (the origin = 0 in elocity space) of this potential U(; c ). Consequently, when a particle has an enough amount of thermal energy, it cannot be captured in the region R c of this potential U(; c ). We hence call it a weakly confining potential. The corresponding FP equation is still linear in a probability density function p but has a nonlinear drift force caused by this nonlinear potential U(; c ) [1]. Introducing such a nonlinear drift force is not a new idea. In fact Lutz [] showed that the anomalous transport of an atom in optical lattice was well described by the linear FP equation with another special type of the nonlinear drift force, and that its stationary state is a q-deformed Gaussian. This work is originally motiated by our recent studies on the information geometric structures of a κ-deformed exponential family [3 6]. Information geometry (IG) [7] proides us a powerful framework for studying a family of probability distributions by identifying the space of probability distributions with a differentiable manifold endowed with a Riemannian metric and an affine connection. An exponential family of probability distributions is mostfamiliar in IG, andit playedanimportant roleespecially in early deelopments of IG. We had studied some IG structures on the κ-deformed exponential families of probability distributions, which are non-gaussians and with heay-tails. For the κ-deformed exponential families, we constructed the suitable statistical manifolds and showed some information geometric structures such as κ-generalized Fisher metrics, θ- and η-potentials, duallyflat structures, κ-generalized diergence functions, and

2 so on [3 5]. In recent years much attention has been paid for studying the statistical physics concerning on some deformed exponential families of probability distributions. Some of them are Tsallis q-deformed exponential [8] and Kaniadakis κ-deformed exponential families [9 11] in nonextensie statistical mechanics. The κ-deformed exponential function is defined by exp κ (x) (κx+ κ x ) 1 κ, () for a real deformed parameter κ. The κ-deformed exponential function and its inerse function, i.e., κ-deformed logarithmic function, are important ingredients of the generalized statistical physics based on κ-entropy [9 11]. The statistical physics based on the κ-deformed functions has been deeloped oer a decade. The reiew article[1] summarizes the theoretical foundations and mathematical formalism generated by the κ-deformed functions, and proides plentiful references including many fields of applications. Recently, the κ-deformed Fourier series[13] and Fourier transform [14] are deeloped, and the open question on the composition law of κ-entropy for statistically independent systems was soled [15]. Despite such efforts for clarifying the physical situations described by the κ-deformed functions, the physical meaning of the deformed parameter κ is still missing. In this contribution we consider the thermal probability distributions for the weakly confining potential U(; c ) of Eq. (1) in the basic framework of statistical physics. In contrast to the well-known standard case of Gaussian distribution for strongly confining potential /, it is found that the quasi-equilibriumthermalprobability distribution for this weakly confining potential U(; c ) is non-gaussian with heay-tails. The corresponding FP equation of the thermal probability distribution for this weakly confining potential describes an anomalous diffusion (or anomalous transport) in a parameter region in which the second moment may dierge [1]. In addition, we relate the canonical distribution of a thermal particle in this weakly confining potential to the κ-exponential distribution (κ-deformed Gaussian), p() exp κ β, (3) for the strongly confining potential / by introducing a suitable re-parameterization of the control parameter c. The paper is organized as follows. The next section proides a brief summary of the κ-deformed functions. In section III we consider a thermal particle in a confining potential and show that for a weakly confining potential of special type(1), its constitutie relation has a elocitydependent nonlinearity. In addition, the stationary state is shown to be a κ-deformed Gaussian. In the region κ > /3of the deformed parameterκ, our model exhibits anomalous transport. Although the mean kinetic energy / dierges, the proposed specific aerage does not dierge, and it satisfies a generalized equipartition relation. Final section is deoted our conclusion and perspectie. II. SOME κ-deformed FUNCTIONS We here briefly summarize some κ-deformed functions and the associated useful relations which are based on the κ-entropy S κ described by S κ dp()ln κ p() = ln κ (p), (4) where is the elocity ofathermal particle, and is the expectation w.r.t p(). Note that S κ is a κ-generalization of the Gibbs-Shannon entropy by replacing the standard logarithm with the κ-logarithm. Here the κ-logarithm ln κ (x) [10, 11] is a deformed function of the standard logarithm ln(x) for a real ariable x > 0, and κ is a real deformed parameter. ln κ (x) xκ x κ κ = 1 κ sinh[ κln(x) ]. (5) Its inerse function is gien by Eq. (), which is also expressed as [ 1 exp κ (x) = exp κ arsinh( κx ) ]. (6) Intheκ 0limit, theκ-exponentialandtheκ-logarithm reduce to the standard exponential exp(x) and logarithm ln(x), respectiely. We next introduce another κ-deformed function: u κ (x) xκ +x κ = cosh [ κln(x) ], (7) which is the conjugate (or co-function) of ln κ x, as similar as that cos(x) is the co-function of sin(x). In the κ 0 limit, this κ-deformed function reduces to the unit constant function u 0 (x) = 1. By using this u κ (x), the deriatie of the κ-exponential is expressed as d dx exp κ(x) = exp κ (x) κ x = exp κ (x) u κ [exp κ (x)], (8) and the deriatie of κ-logarithm is expressed as d dx ln κ(x) = u κ(x) x, (9) respectiely. Since all κ-deformed functions are symmetric under the sign change of the deformed parameter κ to κ, through out this paper we assume that the deformed parameter κ takes a positie real alue.

3 III. THERMAL PDF FOR A WEAKLY CONFINING POTENTIAL We consider a thermal particle in a confining potential U() for the elocity of the particle. It is described by the Langein equation for the oerdamped Brownian dynamics: dx(t) = (t), dt d(t) = α d U()+ζ(t), (10) dt d where α denotes a drift coefficient and ζ(t) is a random force due to random density fluctuations of the enironment around the particle. For the sake of simplicity the mass of the thermal particle is set to m = 1. The randomforceζ(t)canbecharacterizedbyitsfirstandsecond moments ζ(t) = 0, ζ(t)ζ(t ) = αk B Tδ(t t ). (11) It is well known that when U() is a parabolic potential U s () = /, the drift force obeysthe linear constitutie relationf d () = α d d U s() = α, andthatagaussian pdf p() = 1 ] [ β Z(β) exp, (1) is the stationary state of a thermal particle in this parabolic potential. In contrast, we consider the weakly confining potential (1) with a controlling parameter c. Figure 1 shows this weakly confining potential U(; c ) with c = 1. For a / ln U(; 1) FIG. 1. The weakly confining potential U(; c) (1) with the controlling parameter c = 1 (red cure). For the sake of comparison, a parabolic (green cure) and logarithmic (blue cure) potentials are also shown. small much less than c, the potential U(; c ) is well approximated with /, and for a large much lager than c, the potential U(; c ) behaes as c ln, i.e., { U(; c ) / ( c ) c ln ( (13) c) The constitutie relation for this weakly confining potential becomes K(; c ) := α d d U(; α c) =, (14) 4 c wherek(; c ) isthe nonlineardrift forceassociatedwith the potential U(; c ). From the behaior of U(; c ) described in Eq. (13) we see that for a small c the drift force K(; c ) is well approximated with the linear constitutie relation K(; c ) = α, while for a larger c K(; c ) 1/, which means that the particle cannot be captured in the region R c around the minimum point of this potential U(; c ) if the particle elocity is much larger than c. Now, we consider the canonical pdf p() of a thermal particle in the weakly confining potential U(; c ), p() exp[ βu(; c )] = exp [ β c arsinh ( c )], (15) where β is the coldness, or inerse temperature 1/(k B T). At this stage let us re-parametrize the controlling parameter c as c = 1 κβ, (16) where κ is a new real parameter. With this reparametrization, the parameter c is a function of κ and β, i.e., c = c (κ,β). We then find that [ )] ] 1 p() exp ( κβ κ arsinh = exp κ [ β, (17) where we used the relation (6) in the last step. We thus see that the canonical pdf of a thermal particle in the weakly confining potential is the κ-deformed Gaussian. Next we consider the following linear FP equation [1] t p(,t) = ( K(; c )p(,t) )+D p(,t), (18) where K(; c ) denotes the nonlinear drift coefficient in Eq. (14), and D a constant diffusion coefficient. The stationary condition t p s() = 0 leads to lnp s() = K(; c) D This relation becomes = α c D lnp s () = α c D arsinh α D = ( ( arsinh ( c c ) 4 c ). (19) ) +const. (0)

4 We thus obtain the stationary pdf of the FP equation (18) as [ )] α p s () exp c ( D arsinh c. (1) Next we introduce the parameters as κ = D αc, β = α D, () which are consistent with the relation (16). Rewriting Eq. (1) by using these parameters, we hae p s () = 1 [ )] 1 ( κβ Z(β) exp κ arsinh = 1 ) Z(β) exp κ ( β, (3) which is the κ-deformed Gaussian. In this analysis it is remarkable that: i) as shown in Eq. (), the deformed parameter κ is expressed in terms of the microscopic physical parameters α, D and the controlling parameter c which determines the nonlinear constitutie relation (14). In the limit of c the deformed parameter κ reduces to 0; ii) the parameter β is expressed as the ration of the friction coefficient α to the diffusion coefficient D, in analogy with the fluctuation-dissipation relation; iii) for the region κ > of the deformed parameter κ, the integral dp s () dierges, and consequently p s () is no longerapdf since it cannotbe normalized. By using the aboe parameters (), this condition is equialently expressed as c < k B T = 1/β = D/α, which physically means the captured energy c / of the weakly confining potential U(; c ) at = c is smaller than the mean thermal energy k B T/. Then a thermal particle which has an enough high thermal energy cannot be captured in the region R c around the minimum of U(; c ). In other words, for a fixed alue of c, only the limited alue β > 1/c of β is allowed for conergence of the pdf p s (); i) in the deformed parameter region /3 < κ, the second moment = d p s (), (4) of the κ-gaussian becomes infinite, in which the mean kinetic energy / dierges, and it is a hall-mark of anomalous transport. We will discuss this point later. By introducing the Lyapuno functional F[p] as with F[p] U[p] D α SBG [p], (5) S BG [p] = d p(,t)lnp(,t), (6) U[p] d U(; c )p(,t), (7) we show [1] that df dt = d [ U(; c )p+ D ] p p α plnp t = d [U(; c )+ Dα ] (lnp+1) [ K(; c )p+d p ] = d p α [ K(; c )+D p lnp ] 0, (8) which states the time eolution of F[p] is non-increasing, i.e., F[p] is the Lyapuno functional w.r.t. the FP equation (18). Consequently the F[p] is minimized by the stationary solution, which is a κ-generalized Gaussian. Haing described the FP analysis of a thermal particle in the weakly confining potential U(; c ), we next point out the iolation of the standard formula U() = lnz(β), (9) β on the aerage energy for a canonical pdf. It is well known that the aerage energy U() for a canonical pdf p() = 1 Z(β) exp[ βu()], (30) is obtained by the famous relation (9). In the deriation of this formula we assume that a potential U() is independent of β. Howeer this assumption is no longer alid for the weakly confining potential U(; c ), since c = c (κ,β) is a function of κ and β. This weakly confining potential U(; c ) depends not only on but also on β; in other words, U(; c ) is temperature dependent. Consequently the aboe formula (9) is no longer alid. Taking into account of the κ- and β-dependence of c, we can obtain d p s () = lnz(β), (31) 4 β c by the straightforward calculations. The r.h.s. can be easily calculated by using the expression (3) and the analytical expression [10] for the κ-deformed Gaussian integral, as follows Then Z(β) = = π βκ d exp κ ( β ) Γ( 1 κ 1 4 ) ( κ )Γ( 1 κ ) κ 0 π β. (3) β lnz(β) = β ln β = 1 β. (33)

5 The l.h.s. of Eq. (31) is the half of the following aerage energy d d U(; c) = d p s (). (34) 4 c By using the relation 1 β d d exp κ ) ( β = exp κ β, (35) 4 c and integrating by part, the r.h.s. of Eq. (34) becomes ) βz(β) exp κ ( β + 1 dp s (). (36) β Finally, let us consider the case in the parameter region /3 < κ <, in which the second moment of the κ-deformed Gaussian dierges as stated before. In this parameter region, the first term in (36) becomes zero and the stationary pdf p s () is normalized, we then obtain d d U(; c) = 1 β = D α, (37) which reminds us a generalization of equipartition theorem [16], p p H = k B T, (38) where H is the Hamiltonian of a system in thermal equilibrium with the temperature T, and p is a generalized momentum conjugate w.r.t a generalized position q. Thus, while the mean kinetic energy dierges, the aerage energy (34) remains finite and is characterized by the ratio D/α of the diffusion coefficient D to the drift coefficient α. We note that the aerage energy (34) is characterized with the so called un-normalized κ-escort expectation, or κ-canonical expectation [17], which is written in our model by E κ,p [f()] := d f() d dx exp κ(x) x= β exp κ β = d f() κ β ( /) exp κ β = d f(), (39) 4 c where we used relation (8). Then, the κ-canonical expectation of the kinetic energy / is expressed as [ ] E κ,p = d Z(β)p s () = Z(β), (40) 4 β c where the relation (31) was used in the last step. In the limit of κ 0, the relation (40) reduces to the standard relation Z(β) = Z(β). (41) β IV. CONCLUSION AND PERSPECTIVE We hae considered a thermal particle in a weakly confining potential U(; c ) with a control parameter c. We hae shown that its stationary state is a canonical probability distribution which is the κ-deformed Gaussian. In this model the parameter c characterizes the -dependency of the weakly confining potential U(; c ), and consequently determines the nonlinear constitutie relation (14). Then, the deformed parameter κ is determined by c andβ asshownin relations(). In thisway, we proide a possible physical meaning of the deformed parameter κ in this model. Recently, Chaudhuri [18] discussed the stochastic dynamics of the actie Brownian particles which are modeled in terms of nonlinear elocity dependent force. Using the associated FP equation, he deried the expression for the total entropy production. Since our model also described by the nonlinear elocity dependent force, it is worth to further study our model from the iew of the entropy production and associated fluctuation theorem. It is interesting to further study the relation (40) from the iew points of both statistical physics and Information geometry. Another possible further study is to generalize our model by utilizing Naudts φ-exponential function[19]. It is a unification of some deformed-exponential functions. It includes the q- and κ-exponential functions as special cases. The inerse function of the φ-exponential function is called φ-logarithmic function, which is defined by ln φ (x) x 1 ds φ(s), (4) for a positie increasing function φ(s). The κ-logarithm is a special case of Naudt s φ-logarithm with φ(s) = s s κ +s κ. (43) ACKNOWLEDGMENTS The first named author (T.W.) is supported by Japan Society for the Promotion of Science (JSPS) Grants-in- Aid for Scientific Research (KAKENHI) Grant Number JP17K The third named author (H. M.) is partially supported by the JSPS Grants-in-Aid for Scientific Research (KAKENHI) Grant Number JP and JP15K0484.

6 [1] T. Wada, Eur. Phys. J. B (010). [] E. Lutz, F. Renzoni, Nature Physics, 9, (013). [3] A.M. Scarfone, T. Wada, J. Phys. A: Math. Theor. 47, 7500 (014). [4] T. Wada, and A.M. Scarfone, Entropy 17, (015). [5] T. Wada, H. Matsuzoe, and A.M. Scarfone, Entropy 17, (015). [6] A.M. Scarfone, H. Matsuzoe, and T. Wada Entropy 0, 436 (018). [7] S-I. Amari, Information Geometry and Its Applications; Springer: Tokyo, Japan, 016. [8] C. Tsallis, Introduction to Nonextensie Statistical Mechanics: Approaching a Complex World Springer: NY, USA, 009. [9] G. Kaniadakis, A.M Scarfone, Physica A, 305, (00). [10] G. Kaniadakis, Phys. Re. E 66, (00). [11] G. Kaniadakis, Phys. Re. E 7, (005). [1] G. Kaniadakis, Entropy, 15, (013). [13] A.M. Scarfone, Entropy, 17, (015). [14] A.M. Scarfone, Physica A, 480, (017). [15] G. Kaniadakis, A. M. Scarfone, A. Sparaigna, and T. Wada, Phys. Re. E 95, 0511 (017). [16] R. C. Tolman, Phys. Re. 11, 6175, (1918). [17] H. Matsuzoe, and T. Wada, Entropy 17, (015). [18] D. Chaudhuri, Phys. Re. E 94, (016). [19] J. Naudts, Generalized Thermostatistics Springer: Berlin, Germany, 011.