Next Article in Journal
A Simple Mechanism Causing Wealth Concentration
Next Article in Special Issue
On Energy–Information Balance in Automatic Control Systems Revisited
Previous Article in Journal
On Entropy of Probability Integral Transformed Time Series
Previous Article in Special Issue
Differential Invariants of Measurements, and Their Relation to Central Moments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Higher Order Structures in Thermodynamics

1
V.A. Trapeznikov Institute of Control Sciences of Russian Academy of Sciences, 65 Profsoyuznaya Str., 117997 Moscow, Russia
2
Faculty of Physics, Lomonosov Moscow State University, Leninskie Gory, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Entropy 2020, 22(10), 1147; https://doi.org/10.3390/e22101147
Submission received: 30 August 2020 / Revised: 7 October 2020 / Accepted: 10 October 2020 / Published: 12 October 2020
(This article belongs to the Special Issue Thermodynamics, Geometry and Control Theory)

Abstract

:
We present the development of the approach to thermodynamics based on measurement. First of all, we recall that considering classical thermodynamics as a theory of measurement of extensive variables one gets the description of thermodynamic states as Legendrian or Lagrangian manifolds representing the average of measurable quantities and extremal measures. Secondly, the variance of random vectors induces the Riemannian structures on the corresponding manifolds. Computing higher order central moments, one drives to the corresponding higher order structures, namely the cubic and the fourth order forms. The cubic form is responsible for the skewness of the extremal distribution. The condition for it to be zero gives us so-called symmetric processes. The positivity of the fourth order structure gives us an additional requirement to thermodynamic state.

1. Introduction

The geometrical interpretation of thermodynamic systems in equilibrium goes back already to the 19th century [1] and is reflected recently in Reference [2]. In modern terms, it is clear that thermodynamic states are Legendrian submanifolds, i.e., maximal integral manifolds of the structure contact form (for more details, see Reference [3,4,5]) in the contact space where the mentioned structure form is the first law of thermodynamics. Additional structures, such as Riemannian structures on these Legendrian manifolds, were studied in, for example, Reference [6]. Considering thermodynamics in the context of measurement of random vectors [7], one gets both structures as coming from the minimal information gain or Kullback-Leibler divergence principle [8]. Namely, Legendrian manifolds represent averages of measurable quantities (or extremal probability distributions) that are extensive thermodynamic variables, while the Riemannian structure is their variance, i.e., contact and Riemannian structures arise from the first two central moments of random vectors. It is worthy of note that both these structures were widely studied, but, at the same time, the higher order structures, corresponding to central moments of higher order, have not been addressed before.
In this paper, we develop the geometrical approach to thermodynamic states and extend it by considering the third and the fourth order moments and corresponding symmetric forms of the third and the fourth order on Legendrian manifolds of two types of gases, ideal and van der Waals. The third order symmetric form represents the skewness of the extremal probability distribution, as well as thermodynamic processes along which the skewness is equal to zero, we, therefore, call symmetric. We elaborate such processes for ideal gases explicitly, and, in the case of real gases represented by the van der Waals model, we show that there are domains on their Legendrian manifold where there are either three types of such processes, or one. It is natural to require that the fourth order symmetric form must be positive. It is known that the positivity of the variance leads to the notion of thermodynamic phases [7], and the positivity of the fourth order form gives us an additional separation of applicable phases.

2. Geometry, Measurement, Thermodynamics

In this section, we briefly recall how contact geometry naturally appears in the context of measurement, as well as how symmetric k-forms represent kth central moments. For details, we refer to Reference [4,7,9].
Let ( Ω , A , q ) be a probability space, where Ω is a sample space, A is a σ -algebra on Ω , and q is a probability measure. Then, a random vector X : ( Ω , A , q ) W , where W is a vector space, dim W = n is a measurement of x 0 W if E q X = x 0 . To measure another vector x W , one has to choose another measure d p = ρ d q , where ρ is a probability density, such that
E p X = Ω X ρ d q = x , Ω ρ d q = 1 .
To find ρ , we also use the principle of minimal information gain:
I ( ρ ) = Ω ρ ln ρ d q min .
This gives us the extremal probability distribution (see Reference [7]):
ρ = e λ , X Z ( λ ) ,
where λ W * , and Z ( λ ) = Ω e λ , X d q .
Introducing the Hamiltonian H ( λ ) = ln Z ( λ ) and using the first relation in (1), we obtain that the measurement belongs to a manifold
L = x = H λ W × W * ,
which is Lagrangian with respect to the symplectic form
ω = i = 1 n d λ i d x i ,
i.e., ω | L = 0 . Considering the relation x i = H λ i , i = 1 , , n as an equation for λ i , one gets (locally) λ i = λ i ( x ) . The information gain I ( x ) is a function on L which is related with Hamiltonian H ( λ ) as
I ( x ) = H ( λ ( x ) ) + λ ( x ) , x .
Therefore, we have I x i = λ i , and, if u is a coordinate on R , then the submanifold
L ^ = u = I ( x ) , λ i = I x i W × W * × R
is Legendrian with respect to the contact form θ = d u i = 1 n λ i d x i .
The kth moment of random vector X is defined by the well-known relation
m k ( X ) = Ω X k ρ d q .
In coordinates ( λ 1 , , λ n ) , it can be written as
m k ( X ) = Z λ i 1 λ i k Z d λ i 1 d λ i k .
Here, and further, we use the Einstein’s summation convention. Formula (4) can be found inductively using the definition of the partition function Z ( λ ) .
The kth central moment σ k is a symmetric k-from on W:
S k ( W ) σ k = Ω ( X m 1 ( X ) ) k ρ d q .
Applying the standard binomial theorem, one gets the following relation between moments m k and central moments σ k (also see Reference [9]):
σ k = i = 0 k ( 1 ) k i k i m i · m 1 ( k i ) ,
where · stands for the symmetric product.
Let us give a coordinate description of central moments of orders 2,3,4.
  • k = 2 .
    Theorem 1.
    On the Legendrian manifold L ^ the second central moment has the following form in coordinates ( λ 1 , , λ n )
    σ 2 = 2 H λ i 1 λ i 2 d λ i 1 d λ i 2 ,
    and in coordinates ( x 1 , , x n )
    σ 2 = 2 I x i 1 x i 2 d x i 1 d x i 2 ,
  • k = 3 .
    Theorem 2.
    On the Legendrian manifold L ^ the third central moment has the following form in coordinates ( λ 1 , , λ n )
    σ 3 = 3 H λ i 1 λ i 3 d λ i 1 d λ i 3 ,
    and in coordinates ( x 1 , , x n )
    σ 3 = 3 I x i 1 x i 3 d x i 1 d x i 3 ,
  • k = 4 .
    Theorem 3.
    On the Legendrian manifold L ^ the fourth central moment has the following form in coordinates ( λ 1 , , λ n )
    σ 4 = 4 H λ i 1 λ i 4 d λ i 1 d λ i 4 + 3 σ 2 · σ 2 .
Formulae (7), (9), (11) follow directly from substitution of (4)–(6) and using Z = exp ( H ) , while (8), (10) result from formula (2) relating the Hamiltonian H ( λ ) and information gain I ( x ) . Formula for σ 4 in terms of ( x 1 , , x n ) can be found, as well, but we do not provide it here because of its bulkiness.
Recall that, in thermodynamics of gases, x = ( e , v ) W , λ = T 1 , p T 1 W * , I ( x ) = S ( e , v ) , where e and v are specific inner energy and specific volume, respectively, T is temperature, p is pressure, and the contact structure is
θ = d s + T 1 d e + p T 1 d v ,
and the Legendrian manifold is given by
L ^ = s = S ( e , v ) , p = S v S e , T = 1 S e W × W * × R .
It is natural to require that all the even order symmetric forms must be positive on Legendrian manifolds. In case of k = 2 , this condition leads to the notion of phases [4,7]. Further, we will elaborate symmetric forms of orders k = 3 and k = 4 in more detail for ideal and van der Waals models of gases.
It is worth mentioning that central moments σ k are preserved by the affine group Aff ( W ) action. In Reference [9,10], the central moments are used to construct scalar differential invariants of Aff ( W ) .

3. Third Central Moment σ 3

The third central moment σ 3 is nothing but skewness of the probability distribution ρ . From the thermodynamic perspective, the symmetric 3-form σ 3 gives us a special type of thermodynamic processes, along which this form vanishes. From geometrical viewpoint, thermodynamic processes are considered to be contact vector fields preserving the Legendrian manifold L ^ . Assuming that these vector fields have no singular points, one may look for such processes in the form
X = e + q v .
Taking into account that I ( x ) = S ( e , v ) and substituting this relation to formula (10), the condition σ 3 ( X , X , X ) = 0 forces the following equation on the coefficient q:
S v v v q 3 + 3 S e v v q 2 + 3 S e e v q + S e e e = 0 .

3.1. Ideal Gas

For ideal gases, the entropy function is given by
S ( e , v ) = ln e n / 2 v ,
where n is the degree of freedom. Equation (12) in case of ideal gases takes the form
2 q 3 v 3 + n e 3 = 0 .
This equation has a unique real solution; therefore, in case of ideal gases, we have one type of symmetric processes shown in Figure 1. The vector field X takes the form
X = e + n 2 1 / 3 v e v .

3.2. van der Waals Gas

In case of van der Waals gas, we use the reduced thermodynamic variables, in terms of which the entropy function is
S ( e , v ) = ln e + 3 v 4 n / 3 ( 3 v 1 ) 8 / 3 .
Here, Equation (12) is
( 54 e 3 v 6 243 e 2 n v 5 + 243 e 2 n v 4 + 486 e 2 v 5 81 e 2 n v 3 729 e n v 4 + 9 e 2 n v 2 + 729 e n v 3 + 1458 e v 4 243 e n v 2 729 n v 3 + + 27 e n v + 729 n v 2 + 1458 v 3 243 n v + 27 n ) q 3 + + ( 243 e n v 6 + 243 e n v 5 81 e n v 4 + 9 e n v 3 ) q 2 + ( 243 n v 7 + 243 n v 6 81 n v 5 + 9 n v 4 ) q + + 27 n v 9 27 n v 8 + 9 n v 7 n v 6 = 0 .
In general, Equation (13) may have either three or one real roots, depending on the discriminant of the cubic (13). Therefore, we get domains where there are either three symmetric processes, or only one. They are shown in Figure 2. We use here coordinates ( T , v ) , where T is temperature, instead of ( e , v ) .
In this figure, the pink domain corresponds to that on the Legendrian manifold where there are three real roots of (13). The blue line separates the domains where the 2-form σ 2 is positive (above) corresponding to applicable domain (where the conditions of thermodynamic stability hold), and negative (under). This means that, in a pink part of an applicable domain of van der Waals gas, there is a family of 3 symmetric processes, while, in the white part above the blue line, only one of them remains. Thus, we can see that, in this case, there is only one transition from three symmetric processes to one.
In Figure 3, again, we are interested only in an applicable domain which is above the blue line. In this case, pink domains where there are three symmetric processes are separated from each other by white domain where there is only one symmetric process.

4. Fourth Central Moment σ 4

The analysis here is similar to that in the previous section. Consider the vector field
X = x 1 e + x 2 v .
The conditions σ 4 > 0 and σ 2 > 0 yield that the homogeneous polynomials P 1 ( x 1 , x 2 ) = σ 2 ( X , X ) and P 2 ( x 1 , x 2 ) = σ 4 ( X , X , X , X ) in x 1 and x 2 are positive.

4.1. Ideal Gas

In case of ideal gas, one gets
P 1 ( x 1 , x 2 ) = n x 1 2 2 e 2 + x 2 2 v 2 , P 2 ( x 1 , x 2 ) = 3 n ( n + 4 ) 4 e 4 x 1 4 + 3 n v 2 e 2 x 1 2 x 2 2 + 9 v 4 x 2 4 ,
and both conditions hold on the entire Legendrian manifold.

4.2. van der Waals Gas

In this case, the standard analysis shows that the condition σ 2 > 0 holds not everywhere as well as σ 4 > 0 does. The first condition gives us separation of the Legendrian manifold to liquid and gas phases, while the second one is an additional requirement to applicable with respect to σ 2 states. This is shown in Figure 4. In the pink domain, both conditions σ 2 > 0 and σ 4 > 0 hold. The red line separates domains where σ 2 has opposite signs.
We can see that the positivity of σ 4 gives us an additional clarification of the applicability of the van der Waals model.

5. Discussion

This paper presents some new results in equilibrium thermodynamics that arise from the measurement approach to thermodynamics. By measuring random vectors, components of which in thermodynamics are specific energy and volume, one gets a Legendrian manifold representing the average of random vector. The variance σ 2 of random vectors leads to the Riemannian structure on the corresponding Legendrian manifold, and condition of positivity of σ 2 forces a well-known condition of thermodynamic stability of a thermodynamic system in equilibrium. The next natural step is to elaborate higher order central moments, for instance, σ 3 , which is skewness, and σ 4 . They lead to the third and the fourth order symmetric forms on the Legendrian manifold, similarly to that of the second order induced by the second central moment σ 2 . It is clear that the even order forms should be positive on the thermodynamic Legendrian manifolds, which can be interpreted by the following way. Not all the points on the Legendrian manifolds correspond to real physical states, but only those of them where σ 2 , σ 4 , and so on are positive. Here, we have considered only σ 4 , but even this additional (comparing with traditional ones) structure has lead us to new applicable domains on van der Waals thermodynamic Legendrian manifold, while, for an ideal gas, σ 4 is positive everywhere. Namely, we can see from Figure 4 that the critical point ( T crit , v crit ) = ( 1 , 1 ) no longer belongs to the domain where both σ 2 and σ 4 are positive. Instead of this, we get something like a new critical point, which is the maximum of the curve separating domains where σ 4 > 0 and σ 4 < 0 . Note that this curve is of the similar form as that separating domains where σ 2 > 0 and σ 2 < 0 . This means that we get a more accurate applicability condition for van der Waals model, and it may be a considerable contribution to the theory of phase transitions, since, as we know, it is the sign changing of σ 2 that forces the first order phase transitions in van der Waals model [4,5,7].
We may say that applicable domain of real gases is an intersection of areas where all even central moments are positive. We could see that, by adding the fourth central moment, we shrink an applicable domain of the van der Waals model comparing with well-known one where σ 2 > 0 . The natural question is will there remain any domains where all even central moments are positive. Our hypothesis is that there will be such domains, but theoretical proof of this statement is an open question.

Author Contributions

Conceptualization, V.L.; Formal analysis, M.R.; Investigation, V.L. and M.R.; Writing—original draft, M.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by the Russian Foundation for Basic Research (project 18-29-10013) and by the Foundation for the Advancement of Theoretical Physics and Mathematics “BASIS” (project 19-7-1-13-3).

Acknowledgments

The authors would like to thank Eivind Schneider for his fruitful suggestions during the preparation of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gibbs, J.W. A Method of Geometrical Representation of the Thermodynamic Properties of Substances by Means of Surfaces. Trans. Conn. Acad. 1873, 1, 382–404. [Google Scholar]
  2. Mrugala, R. Geometrical formulation of equilibrium phenomenological thermodynamics. Rep. Math. Phys. 1978, 14, 419–427. [Google Scholar] [CrossRef]
  3. Arnold, V.I. Mathematical Methods of Classical Mechanics; Springer: New York, NY, USA, 1989. [Google Scholar]
  4. Kushner, A.; Lychagin, V.; Roop, M. Optimal Thermodynamic Processes for Gases. Entropy 2020, 22, 448. [Google Scholar] [CrossRef]
  5. Lychagin, V.; Roop, M. Critical phenomena in filtration processes of real gases. Lobachevskii J. Math. 2020, 41, 382–399. [Google Scholar] [CrossRef]
  6. Ruppeiner, G. Riemannian geometry in thermodynamic fluctuation theory. Rev. Mod. Phys. 1995, 67, 605–659. [Google Scholar] [CrossRef]
  7. Lychagin, V. Contact Geometry, Measurement, and Thermodynamics. In Nonlinear PDEs, Their Geometry and Applications; Kycia, R., Schneider, E., Ulan, M., Eds.; Birkhäuser: Cham, Switzerland, 2019; pp. 3–52. [Google Scholar]
  8. Kullback, S. Information Theory and Statistics; John Willey & Sons.: Chichester, UK, 1959. [Google Scholar]
  9. Schneider, E. Differential invariants of measurements, and their relation to central moments. Entropy 2020, 22, 1118. [Google Scholar] [CrossRef]
  10. Schneider, E. Differential Invariants in Thermodynamics. In Nonlinear PDEs, Their Geometry and Applications; Kycia, R., Schneider, E., Ulan, M., Eds.; Birkhäuser: Cham, Switzerland, 2019; pp. 223–232. [Google Scholar]
Figure 1. Symmetric process for ideal gas.
Figure 1. Symmetric process for ideal gas.
Entropy 22 01147 g001
Figure 2. Symmetric processes for van der Waals gas, n = 3 .
Figure 2. Symmetric processes for van der Waals gas, n = 3 .
Entropy 22 01147 g002
Figure 3. Symmetric processes for van der Waals gas, n = 13 .
Figure 3. Symmetric processes for van der Waals gas, n = 13 .
Entropy 22 01147 g003
Figure 4. Applicable domains for van der Waals gas, n = 3 .
Figure 4. Applicable domains for van der Waals gas, n = 3 .
Entropy 22 01147 g004

Share and Cite

MDPI and ACS Style

Lychagin, V.; Roop, M. On Higher Order Structures in Thermodynamics. Entropy 2020, 22, 1147. https://doi.org/10.3390/e22101147

AMA Style

Lychagin V, Roop M. On Higher Order Structures in Thermodynamics. Entropy. 2020; 22(10):1147. https://doi.org/10.3390/e22101147

Chicago/Turabian Style

Lychagin, Valentin, and Mikhail Roop. 2020. "On Higher Order Structures in Thermodynamics" Entropy 22, no. 10: 1147. https://doi.org/10.3390/e22101147

APA Style

Lychagin, V., & Roop, M. (2020). On Higher Order Structures in Thermodynamics. Entropy, 22(10), 1147. https://doi.org/10.3390/e22101147

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop