Open access peer-reviewed chapter

Differential Geometry and Macroscopic Descriptions in Nonequilibrium Process

Written By

Claudia B. Ruscitti, Laura B. Langoni and Augusto A. Melgarejo

Submitted: 12 December 2019 Reviewed: 26 March 2020 Published: 13 May 2020

DOI: 10.5772/intechopen.92274

From the Edited Volume

Advances on Tensor Analysis and their Applications

Edited by Francisco Bulnes

Chapter metrics overview

667 Chapter Downloads

View Full Metrics

Abstract

The method of Riemannian geometry is fruitful in equilibrium thermodynamics. From the theory of fluctuations it has been possible to construct a metric for the space of thermodynamic equilibrium states. Inspired by these geometric elements, we will discuss the geometric-differential approach of nonequilibrium systems. In particular we will study the geometric aspects from the knowledge of the macroscopic potential associated with the Uhlenbeck-Ornstein (UO) nonequilibrium process. Assuming the geodesic curve as an optimal path and using the affine connection, known as α-connection, we will study the conditions under which a diffusive process can be considered optimal. We will also analyze the impact of this behavior on the entropy of the system, relating these results with studies of instabilities in diffusive processes.

Keywords

  • nonequilibrium processes
  • Uhlenbeck-Ornstein process
  • statistical manifold
  • α-connections
  • macroscopic potential

1. Introduction

The use of Riemannian geometry associated with the space of thermodynamic equilibrium states has been successful. In this framework the geometric elements are constructed from the knowledge of the thermodynamic potential. In this sense, one of the most interesting ideas is associated with the study of the phase transitions visualized by means of the singularities of the scalar curvature [1, 2]. As a geometric-differential approach of nonequilibrium systems, we consider in our study the geometric properties of a statistical manifold associated with trajectory-dependent entropy [3]. In statistics mechanics it is known that any statistical system has an associated metric-affine manifold having a special affine connection whether it is in equilibrium or not. The affine connection, called the α-connection [4], is a generalization of the Levi-Civita connection in Riemannian geometry, and in the case α=0 the metric-affine manifold reduces to the so-called Riemannian manifold. A first element that appears for nonequilibrium systems is the visualization of the phase transitions through the curvature tensor [5]. A second issue is the study of temporal evolution of a statistical system in order to study the optimal evolution. In other words, it is the analysis of the dynamic behavior of the system and the study of the conditions that are optimal.

In the previous context, this chapter will focus on the analysis of the optimal evolution. In particular we consider the Uhlenbeck-Ornstein (UO) nonequilibrium process described by the probability density function (PDF) solution of the Fokker-Planck Equation [6]. From this probability density function, we build a two-dimensional metric-affine manifold in the coordinates μσ, where μ is the mean and σ the standard deviation. In this coordinates and for the connections α=0 and α=1, the system evolves on a geodesic curve [7].

However, due to simplicity in geometric construction, we are interested in studying the behavior of the system in coordinates θ1θ2, where θ1=1/2σ2 and θ2=μ/σ2. In these coordinates the probability density function belongs to the exponential family, and using this formal expression and analogously with the equilibrium probability density function, we build a macroscopic potential ψθ1θ2 for the UO process. From the geometry constructed using the ψθ1θ2 potential, we show that for α=3 and α=2, the system evolves on a geodesic curve. In particular for α=3, we show that the manifold is flat and for the steady state the macroscopic potential and entropy have the same functional dependence. Thinking the geodesic curve as an optimal trajectory, our results allow us to conjecture that the entropy describes the steady state of an optimal evolution.

In the second section of this chapter, we summarize the most relevant aspects of the theory of the statistical manifold. The geometric development associated with the fundamental solution of the Fokker-Planck equation of UO process is found in the third section. The fourth section is devoted to the construction of the potential. In the fifth section, we analyze the geometric relationship between macroscopic potential and entropy. In the sixth section, we present our conclusions and perspectives.

Advertisement

2. Elements of statistical manifold

In this section we briefly review the information of geometrical theory [4] that is used to analyze geometrically a family of probability density functions (PDF) and its application to thermodynamics. Let pxθ be a PDF described by a random variable x and parameters θ=θ1θ2θm that characterize a system. The set of PDFs

M=pxθ:θΩmE1

becomes an m-dimensional statistical manifold having θi coordinates. According to the information geometrical theory, we can make a metric tensor gij

gijθ=Eilxθjl(xθ)=Eijlxθ,E2

where lxθ=lnpxθ and E. means the expectation operation with respect to pxθ. The last expression is obtained by the use of the normalization condition Eilxθ=0. This metric tensor is the Fisher information matrix in information theory.

In the statistical manifold M, we can introduce a natural derivative of the vector field B toward the tangent vector A, denoted by AαB. It is obtained through the covariant coefficients [7]:

Γijkαθ=Eijlxθ+1α2il(xθ)jl(xθ)klxθ.E3

We now restrict our attention to a special family of probability density function, called an exponential family, which is described by [8].

pxθ=expCx+imθiFixψθ,E4

where Cx and Fix are arbitrary functions of x and ψθ is a function of θi coordinates.

Particularly for a probability density function belonging to the exponential family, from Eqs. (2) and (3) the covariant coefficients and metric tensor are written as [9, 10].

Γijkαθ=1α2ijkψθ,E5
gijθ=2ψθθiθj.E6

These coefficients, for α, determine a one-parameter family of affine connections, and each element of this family is called an α-connection. An affine connection allows one to compare vectors in nearby tangent spaces [9]. Moreover, the covariant coefficients satisfy the following relation:

Γijkα=gkmΓijmα,E7

where gkm is the inverse matrix of the metric and gij and Γijkα are the contravariant coefficients.

In the case α=0, the coefficients reduce to the Levi-Civita’s connection:

Γijk0=12igjk+jgikkgij.E8

From Eq. (3) the curvature tensor for an α-connection is written as

Rijkmα=iΓjksαjΓiksαgsm+ΓirmαΓjkrαΓjrmαΓikrα.E9

Since Γijkαθ vanishes identically for α=1, any exponential family of PDF constitutes an uncurved space when the α=1 connection is used. The one connection is therefore called the exponential connection as we mentioned above [11]. The relationships (Eqs. (5) and (6)) lead to the simplification of the geometric construction associated with the probability density function (Eq. (4)).

For systems in thermodynamic equilibrium, the parameters θi may include inverse temperature, chemical potential, pressure, magnetic field, and so on. Also, ψθ represents the thermodynamic potential of the system [8]. It is worth noting that, if we consider the equilibrium density functions and α=0, this formalism reproduces the geometric structure found by Ruppeiner for the space of equilibrium states. In the Ruppeiner formalism, the metric is constructed using the theory of fluctuations, and the geometric elements are obtained as the second derivatives of the corresponding thermodynamic potential [2]. In this paper we have chosen the approach of statistical manifold because it allows us to address nonequilibrium problems.

Advertisement

3. Fokker-Planck equation and macroscopic potential

A diffusion process can be thought of as a process of Uhlenbeck-Ornstein. The Uhlenbeck-Ornstein process is a stochastic process that, roughly speaking, describes the velocity of a massive Brownian particle under the influence of friction. The probability density function Pxt of the Uhlenbeck-Ornstein process satisfies the Fokker-Planck equation [6]:

Pt=1τxxx0P+D2Px2.E10

The fundamental solution of this linear parabolic partial differential equation, and the initial condition consisting of a unit point mass at location y, is:

Pxt=12πDτ1e2t/τexp12xx0yx0et/τ21e2t/τE11

which is the Gaussian density function with mean

μ=x0+yx0et/τE12

and variance

σ2=1e2t/τE13

where x0 represents the average length of the displacement, D is the diffusion coefficient, and τ is a characteristic time. Without loss of generality, in the rest of the work, we consider y=0.

Considering Eqs. (12) and (13), we think the function (Eq. (11)) as a probability density function dependent on two parameters μ and σ, formally:

pxμσ=12πσ2exμ2/2σ2.E14

Considering these parameters as the coordinates μσ of the manifold M and taking into account Eqs. (2) and (4), it can be seen that for α=0 and α=1, the system evolves on a geodesic curve [7].

Inspired by the simplicity of relations (Eqs. (5) and (6)), we use an alternative description of the UO process through the coordinates:

θ1θ2=12σ2μσ2.E15

In these coordinates the probability density function (Eq. (14)) belongs to the exponential family and is written as

pxθ1θ2=expθ1x2θ2xψθ1θ2E16

where

ψθ1θ2=θ224θ1lnθ1π.E17

By analogy with the probability density functions for systems in equilibrium, we will call the relationship (Eq. (17)) “nonequilibrium potential” [8].

Using the potential (Eq. (17)) we calculate the coefficients (Eq. (5))

Γ111α=α121θ13+3θ222θ14,Γ222α=0Γ112α=Γ121α=Γ211α=α12θ22θ13Γ122α=Γ221α=Γ212α=α14θ12E18

Using Eqs. (18) and (19), we calculate the curvature tensor component in the coordinates θ1θ2

R1212αθ1θ2=α3α18θ13.E19

From Eq. (19) we observe that there are two α values for which the manifold is flat, α=3 and α=1. The case α=1 is a direct consequence of Eq. (5).

In a manifold with a connection, we can generalize the straight line of Euclidean geometry. The generalized straight line is called geodesic, and it is defined by the characteristic that its tangent vector does not change its direction. It satisfies the Eq. (11)

d2θidu2+Γjkidθkdudθjdu=d2u/ds2du/ds2dθiduE20

with an arbitrary parameter u. The special parameter s is called the affine parameter. We note that if u is a linear transformation of s, the right-hand side vanishes. Let us now prove that the curve defined by Eq. (15) of the UO process is a geodesic. We choose the Newtonian time t for the arbitrary parameter u. From Eq. (20) for the study of geodetic curves, we consider that the temporal dependence of the coordinates θ1θ2 is given by

θ1t=121e2t/τ,θ2t=x01et/τ1e2t/τ.E21

The tangent vector coordinate dθi/dt and the acceleration coordinate d2θi/dt2 are

dθ1dt=cosh2t/τ4Dτ2,d2θ1dt2=cotht/τcsch2t/τ2Dτ2E22

and

dθ2dt=x0sech2t/2τ4Dτ2,d2θ2dt2=2x0sinh4t/2τcsch3t/τDτ3.E23

Substituting Eqs. (22) and (23) in Eq. (20) for the coordinate θ1, we have

12Dτ22α1α3cotht/τα1x02=d2t/ds2dt/ds2,E24

and for the coordinate θ2

cscht/τ2Dτ221α2cosht/τ+α1sinht/τ+cscht/τ2Dτ21+et/τα1x02=d2t/ds2dt/ds2.E25

The results (Eqs. (24) and (25)) are compatible if we choose α as

α=1+2x02.E26

Taking as reference the geometric study in the coordinates μσ [7], we consider two possible relationships between the characteristic time τ and the other two parameters, τ=x02/D or τ=x02/2D. The first case leads to a flat geometry, that is, if τ=x02/D, we have that α=3 and for Eq. (19)R12123θ1θ2=0. This case is equivalent to the connection α=1 in the coordinates μσ, where the manifold is also flat and the affine time s and Newtonian time t are related by the differential equation

1τ=d2t/ds2dt/ds2E27

whose solution has the form

t=τlns++b.E28

On the other hand if τ=x02/2D, we have that α=2 and the affine and Newtonian time are now related by

cotht/ττ=d2t/ds2dt/ds2.E29

Due to its nature, the equation (Eq. (29)) has been solved numerically. For each τ we find solutions that make sense in the context of our problem. In Figure 1 we show two examples of these solutions.

Figure 1.

t as a function of s for two values of τ. τ=10 in the solid line and τ=1 in the dashed line.

In the coordinates μσ, Eqs. (27) and (29) correspond to the choices α=1 and α=0, respectively [7]. On the other hand in the context of geometric construction from the potential (Eq. (17)), the main interpretation of the equation (Eq. (26)) is that the temporal evolution described by Eq. (21) coincides with a geodesic curve of space θ1θ2. In this sense we find two relationships τ=x02/D and τ=x02/2D between the parameters x0, D, and τ for which the evolution is optimized.

It is interesting to note that if we assume that v=x0/τ, the relationship τ=x02/D leads us to D=v2τ. This is a well-known and a widely used relationship, for example, in the understanding of bacterial mobility through a diffusive model. In this case the characteristic time τ represents the tumbling time, and ν represents the average speed of bacteria [12].

Advertisement

4. Optimal trajectory and entropy

Thinking the geodesic curve as an optimal trajectory, from Eq. (21) we will analyze the optimal evolution of the system in terms of the behavior of the macroscopic potential (Eq. (17)) and the entropy function. To begin with this discussion, we first notice that the common definition of a nonequilibrium Gibbs entropy

Hθ=pxθlnpxθdx=hxθE30

suggests to define a trajectory-dependent entropy for the particle (or “system”) hxθ=lnpxθ [3]. The trajectory-dependent entropy is related to the function lxθ used in the construction of the metric of the statistical manifold M by means of the relation hxθ=lxθ. In our analysis we are interested in the entropy averaged over trajectories Hθ.

Moreover, given that the probability density function (Eq. (16)) belongs to the exponential family, we can calculate the entropy by means of the relation [13].

Hθ=ψθ1ψθ1θ2ψθ2.E31

In any case, using Eqs. (30) or (31), we obtain

Hθ=12lnθ1π.E32

In order to study the temporal evolution of potential ψ and entropy H, we use the relations Eqs. (12) and (13) in Eq. (21). In particular, we focus on investigating the asymptotic behavior, that is, the behavior for t. In this limit and using Eq. (21), we have that θ1=1/2 and θ2=x0/. Therefore, the asymptotic behavior of the potential and entropy is written as

ψ=x022+12ln2πDτE33

and

H=12+12ln2πDτ.E34

If we think that in the coordinates θ1θ2, with α=3, the system evolves on a geodesic of a flat manifold (R12123θ1θ2=0), we have τ=x02/D, and both the potential ψ and the entropy H tend to the same asymptotic value. In other words the entropy describes the asymptotic behavior of an optimal evolution:

ψ=H=12+ln2πDτ.E35

In this last result, it can be seen that the asymptotic behavior depends fundamentally on the transport coefficient D. In this sense, the relevant variable in the description of the process is the diffusion coefficient D. It is interesting to note that the relations (Eq. (21)) are not invariant when we change t by t. In this context, the evolution of the system is irreversible.

While in the coordinates θ1θ2 are two α-connections which leads to the system on a geodesic curve (α=3 and α=2), only α=3 makes that the macroscopic potential and entropy have the same functional dependence. It is interesting to note that for this choice of α, the manifold in the coordinates θ1θ2 is flat.

In the study of processes, the speed at which they occur is important. In this sense we will study the rate of change of potential:

ψ̇=x024Dτ2secht2τ+12τcothtτ1==12τsecht2τ+12τcothtτ10,E36

and the entropy

Ḣ=12τcothtτ10.E37

The equal sign in Eqs. (36) and (37) corresponds to the steady state [14]. In this regard, we associate the asymptotic behavior with the steady state, and we can indicate that the entropy describes the steady state of an optimal evolution.

From the perspective of the probability density function (PDF), when the diffusion coefficient takes small values, the distribution (Eq. (14)) diffuses from a uniform state to a sharp concentrated state, that is, from a stable state to an unstable state [15]. Although in this chapter we have not studied the behavior of the curvature tensor in the μσ coordinates, these instabilities observed from the perspective of the PDF can be expressed in the singularity of the curvature tensor R12120μσ. An example of this behavior is the formation of density patterns in the populations of self-propelled bacteria whose mobility can be investigated in terms of the diffusion coefficient D [5]. We also observe from Eq. (35) that for D0, the entropy and the macroscopic potential have a singular behavior. The behavior that we observe both the tensor R12120μσ and the potential ψ allows us to conjecture that the instabilities of the system can also be observed in the singularities or discontinuities of the macroscopic potential.

Advertisement

5. Discussion and perspectives

In this chapter we have studied the PDF (Eq. (14)) that is the fundamental solution of the Fokker-Planck equation associated with the UO process. Our main interest was relating the geometric aspects of the process with the steady-state behavior. In our analysis we used the theoretical framework of the statistical manifold M with α-connections for two different coordinates μσ and θ1θ2 [4, 7]. In the first case, there exist two interesting values of α, namely, α=1 and α=0, for which the process evolves on a geodesic of space μσ with different values of τ. However, in the search for a simpler geometric construction, we find that for the coordinates θ1θ2, we can define a macroscopic potential ψθ1θ2, and the values of α that lead to the system evolving on a geodesic curve are α=3 and α=2. In our study we show that the connection α=1 for the coordinates μσ corresponds to the connection α=3 for the coordinates θ1θ2 in the sense that both connections lead to a flat curvature and the same relationship between the parameters of the system. An important consequence of this behavior is that when the system evolves over geodesics, the macroscopic potential ψ and entropy H have the same functional dependence in the steady state. If we think of the geodesic curve as an optimal trajectory, our results allow us to conjecture that the entropy describes the steady state of an optimal evolution in a flat manifold.

Additionally and poorly developed in this chapter is the use of geometric aspects in the study of instabilities in nonequilibrium system. In equilibrium thermodynamics this information is contained in the scalar curvature of the manifold of equilibrium states [2]. For nonequilibrium problems, the instabilities are associated with the singularities or discontinuities of the curvature tensor [5, 15]. In the case of diffusive problems, the instabilities can be found by studying the singularities of R12120μσ. In terms of diffusivity, we see that as it decreases, the system has different macroscopic behavior. An example of this behavior is the formation of density patterns in the populations of self-propelled bacteria whose mobility can be investigated in terms of the diffusion coefficient D [12]. From the perspective of the macroscopic potential ψ, the instabilities can be associated to singularities or discontinuities of ψ. In this sense, from a wider point of view, we consider that the potential ψ represents an alternative way to study the phase transitions in nonequilibrium systems.

Advertisement

Acknowledgments

This work was supported by the Universidad Nacional de La Plata, Argentina (UNLP). AM, LL, and CR are professors at the UNLP.

References

  1. 1. Ruppeiner G. Thermodynamics: A Riemannian geometric model. Physical Review A. 1979;20:1608-1613
  2. 2. Ruppeiner G. Riemannian geometric approach to critical points: General theory. Physical Review E. 1998;57:5135-5145
  3. 3. Seifert U. Entropy production along a stochastic trajectory and an integral fluctuation theorem. Physical Review Letters. 2005;95:040602
  4. 4. Amari S. Differential-geometrical methods in statistics. In: Lecture Notes in Statistics. Vol. 28. Heindelberg: Springer; 1985
  5. 5. Melgarejo A, Langoni L, Ruscitti C. Instability in bacterial populations and the curvature tensor. Physica A: Statistical Mechanics and its Applications. 2016;458:189-193
  6. 6. Uhlenbeck G, Ornstein L. On the theory of the Brownian motion. Physics Review. 1930;36:823-841
  7. 7. Ruscitti C, Langoni L, Melgarejo A. Description of instabilities in Uhlenbeck–Ornstein process. Physica Scripta. 2019;94:115010
  8. 8. Brody D, Hook D. Information geometry in vapour-liquid equilibrium. Journal of Physics A: Mathematical and Theoretical. 2009;42:023001
  9. 9. Amari S. Differential geometry of curved exponential families-curvatures and information loss. Annals of Statistics. 1982;10:357-385
  10. 10. Oshima H, Obata T, Hara H. Riemann scalar curvature of ideal quantum gases obeying Gentile’s statistics. Journal of Physics A: Mathematical and General. 1999;32:6373-6383
  11. 11. Zhang J. On monotone embedding in information geometry. Entropy. 2015;17:4485-4499
  12. 12. Cates M, Marenduzzo D, Pagonabarraga I, Tailleur J. Arrested phase separation in reproducing bacteria creates a generic route to pattern formation. Proceedings of the National Academy of Sciences. 2010;107:11715-11720
  13. 13. Jaynes ET. Information theory and statistical mechanics. Physics Review. 1957;106:620-630
  14. 14. Kondepudi D, Prigogine I. Modern Thermodynamics from Heat Engines to Dissipative Structures. 2nd ed. United Kingdom: John Wiley & Sons; 2015
  15. 15. Obata T, Hara H, Endo K. Differential geometry of non-equilibrium processes. Physical Review A. 1992;45:6997-7001

Written By

Claudia B. Ruscitti, Laura B. Langoni and Augusto A. Melgarejo

Submitted: 12 December 2019 Reviewed: 26 March 2020 Published: 13 May 2020