Thermodynamic Curvature of the Binary van der Waals Fluid

Abstract

The thermodynamic Ricci curvature scalar R has been applied in a number of contexts, mostly for systems characterized by 2D thermodynamic geometries. Calculations of R in thermodynamic geometries of dimension three or greater have been very few, especially in the fluid regime. In this paper, we calculate R for two examples involving binary fluid mixtures: a binary mixture of a van der Waals (vdW) fluid with only repulsive interactions, and a binary vdW mixture with attractive interactions added. In both of these examples, we evaluate R for full 3D thermodynamic geometries. Our finding is that basic physical patterns found for R in the pure fluid are reproduced to a large extent for the binary fluid.

1. Introduction

The thermodynamic Ricci curvature scalar R has yielded a number of interesting results in the study of fluids [1]. However, to this point, the great majority of the calculations of R have been made for pure fluids. Much less examined has been R for binary fluids. Pure fluids offer many research topics in a relatively simple thermodynamic geometric scenario. One element of this simplicity is that pure fluids may be represented by a two-dimensional (2D) thermodynamic Riemannian geometry, where just the scalar R gives the full curvature picture.

The dimension of the thermodynamic phase space grows by one for each added fluid component, and, as the dimension grows, the curvature rapidly becomes more complicated. For example, the binary fluid corresponds to a three-dimensional (3D) phase space, where there are six independent components of the full Riemannian curvature tensor [2]. Which of these components do we focus on for physical interpretation? In this paper, our results suggest that the curvature scalar R remains the fundamental physical quantity. (We offer no rigorous proof, but we do note that classical general relativity, based on four-dimensional Riemannian geometry, may be expressed as a variational principle based on the Hilbert action, with R as the Lagrangian [2]. This variational principle is certainly not in play in this paper, but we do hope that this significance for R translates to the thermodynamic scenario as well).

We calculate R for two thermodynamic examples represented by a 3D thermodynamic geometry: (1) a binary van der Waals (vdW) fluid with just repulsive interactions, and (2) a binary vdW fluid with attractive interactions added. To add conceptual structure to our presentation, we discuss to what extent the 3D results for R follow the same physical principles as the 2D ones. As we show, there is a great deal of correspondence. (Scenarios with a thermodynamic geometric dimension greater than two have also been considered for black hole thermodynamics; see Sahay for a brief review [3]. This theme is beyond the scope of this paper.)

A few cases have been worked out for R for the full 3D thermodynamic geometry. Ruppeiner and Davis [4] worked out R for the multicomponent ideal gas, with an arbitrary number of components. Kaviani and Dalafi-Rezaie [5] worked out R for the paramagnetic ideal gas, where, in addition to the temperature and the density, there is an external organizing magnetic field. Erdem [6] worked out R for the antiferromagnetic Ising model with a temperature and two ordering fields, with special attention to the critical phenomena near the Neel point.

One may also approach 3D thermodynamic geometries by working out R over 2D slices of the full 3D thermodynamic geometry. This has physical relevance if, for example, one of the independent variables is irrelevant to a phase transition. Ginoza [7] considered the binary fluid mixture in generality, and calculated R 2D along surfaces of constant T. Jaramillo-Gutiérrez et al. [8] calculated a constrained R 2D by fixing the sum of the mole numbers. These authors made several comments about how their constrained R relates to that of the pure fluid. Generally, the calculation is simplified in 2D, but, in an age of powerful mathematics software, considerations of simplicity should maybe no longer be so much of a driving force.

2. Expectations for 3D Outcomes for R

For guidance as to what 3D thermodynamic geometry might tell us, let us pose some questions motivated by findings in the pure fluid [9,10,11,12]: (a) In the ideal gas limit, is $\left|R\right|$ either zero or small? “Small” means $\left|R\right|\rho <1$, with $\rho$ the number density. In this event, the volume measured by $\left|R\right|$ is less than the average volume occupied by a single particle, and we are under the low $\left|R\right|$ limit [10,11]. Thermodynamics is challenged at such a size scale, and a precise physical interpretation for R is harder to come by (though researchers do try if the opportunity arises). (b) As a critical point is approached, does $\left|R\right|$ diverge in proportion to the correlation volume ${\xi }^3$? (c) Near a critical point, is the sign of R negative? A negative sign is consistent with effectively attractive interactions. (d) Are the values of R in two coexisting phases equal to each other near a critical point? (e) Are there interesting regimes of positive R? A positive sign is consistent with effectively repulsive interactions. (f) Does the binary repulsive fluid have instances of anomalous negative R, such as is present for the hard-sphere pure fluid [13]?

We find considerable physical correspondence between R for the pure fluid and for the binary fluid.

3. Thermodynamic Geometry of a Binary Fluid

For a binary fluid, the Helmholtz free energy per volume may be written in terms of its appropriate coordinates [4,14],

$$f=f(T,{\rho }_1,{\rho }_2),$$

(1)

where T is the temperature, and ${\rho }_1=n_1/V$ and ${\rho }_2=n_2/V$ are the molar densities of the two components 1 and 2. Here, $n_1$ and $n_2$ are the mole numbers of the components, and V is the volume.

The thermodynamic entropy metric originates from the thermodynamic fluctuation theory [1,15]. Consider a finite open subsystem ${\mathcal{A}}_V$, with fixed volume V, of an infinite closed thermodynamic fluid system $\mathcal{A}$; see Figure 1. The thermodynamic state $(T,{\rho }_1,{\rho }_2)$ of ${\mathcal{A}}_V$, fluctuates in time about the state of $\mathcal{A}$. The Gaussian fluctuation probability density is proportional to

$$exp\left[-\frac{V}{2}{\left(\mathsf{\Delta }\ell \right)}^2\right],$$

(2)

where the entropy metric is [1]

$${\left(\mathsf{\Delta }\ell \right)}^2=\frac{1}{k_{B}T}\left(\frac{\partial s}{\partial T}\right){\left(\mathsf{\Delta }T\right)}^2+\frac{1}{k_{B}T}\sum _{i,j=1}^2\left(\frac{\partial {\mu }_i}{\partial {\rho }_j}\right)\mathsf{\Delta }{\rho }_i\mathsf{\Delta }{\rho }_j.$$

(3)

Here, the entropy per volume of ${\mathcal{A}}_V$ is

$$s=-\frac{\partial f}{\partial T},$$

(4)

and the chemical potentials of the two fluid components are

$${\mu }_i=\frac{\partial f}{\partial {\rho }_i}.$$

(5)

$\mathsf{\Delta }T$, $\mathsf{\Delta }{\rho }_1$, and $\mathsf{\Delta }{\rho }_2$ are the differences in the temperature and density coordinates of ${\mathcal{A}}_V$ and $\mathcal{A}$. In addition, $k_B$ is Boltzmann’s constant. The thermodynamic metric is employed here in fluctuation theory. It has also been used in finite time thermodynamics [16] as a measure of dissipation.

Let us emphasize a point not always appreciated in the metric geometry of thermodynamic fluctuations. A major project in this geometry is to calculate the thermodynamic curvature R. We want R to tell us something about the intrinsic properties of the material comprising the system $\mathcal{A}$. This goal naturally requires the use of an open subsystem ${\mathcal{A}}_V$, so as to leave particles free to move in or out of ${\mathcal{A}}_V$, unimpeded by any surrounding artificial wall or membrane. Such a physical constraint would change the value of R in a way involving more than just the properties of the particles; this should be avoided. As emphasized by Callen [14], fluid thermodynamics requires us to set one of the thermodynamic variables aside as the fixed subsystem scale. For an open system, this scale is the volume V. We recommend that authors always work with open subsystems for calculating R.

In the coordinates $(T,{\rho }_1,{\rho }_2)$, the metric tensor g is composed of five nonzero elements, which can be read off from Equation (3):

$${\displaystyle \mathbf{g}=\frac{1}{k_{B}T}\left(\begin{array}{ccc}-\frac{{\partial }^{2}f}{\partial T^2}& 0& 0\\ 0& \frac{{\partial }^{2}f}{\partial {\rho }_1^2}& \frac{{\partial }^{2}f}{\partial {\rho }_1\partial {\rho }_2}\\ 0& \frac{{\partial }^{2}f}{\partial {\rho }_2\partial {\rho }_1}& \frac{{\partial }^{2}f}{\partial {\rho }_2^2}\end{array}\right).}$$

(6)

The fourth-rank curvature tensor has elements [2]

$$R_{\beta \gamma \delta }^{\alpha }={\mathsf{\Gamma }}_{\beta \gamma ,\delta }^{\alpha }-{\mathsf{\Gamma }}_{\beta \delta ,\gamma }^{\alpha }+{\mathsf{\Gamma }}_{\beta \gamma }^{\mu }{\mathsf{\Gamma }}_{\delta \mu }^{\alpha }-{\mathsf{\Gamma }}_{\beta \delta }^{\mu }{\mathsf{\Gamma }}_{\gamma \mu }^{\alpha },$$

(7)

where the Greek indices range from $1,2,3$, and denote the coordinates $(T,{\rho }_1,{\rho }_2)$, respectively. The Christoffel symbols are defined as

$${\mathsf{\Gamma }}_{\beta \gamma }^{\alpha }=\frac{1}{2}{g}^{\alpha \mu }(g_{\mu \gamma ,\beta }+g_{\mu \beta ,\gamma }-g_{\beta \gamma ,\mu }),$$

(8)

where $g^{\alpha \beta }$ denotes the elements of the inverse of the metric g with elements $g_{\alpha \beta }$. We use the comma notation (e.g., $,\beta )$ to denote the partial derivative with respect to some specific coordinate. Repeated indices are summed over. The second-rank Ricci tensor is given by

$$R_{\alpha \beta }=R_{\alpha \mu \beta }^{\mu },$$

(9)

and the Ricci curvature scalar is

$$R=g^{\mu \nu }{R}_{\mu \nu }.$$

(10)

It is straightforward to show that R has units of volume per particle for the binary fluid, the same units as for the pure fluid. These units alone label R as a measure of some sort of size scale within the system. The definition of the Riemannian curvature tensor is ambiguous with respect to a sign, and we use the sign convention of Weinberg [17], where the two spheres have negative R.

4. Thermodynamic Stability

Fluctuations must satisfy thermodynamic stability, requiring maximum entropy of $\mathcal{A}$ in equilibrium. This is obtained if the line element in Equation (3) is positive definite for all fluctuations. A necessary and sufficient condition for thermodynamic stability is that the metric coefficients in Equation (6) satisfy the three following conditions [1,18]:

$$g_{11}>0,$$

(11)

$$\left|\begin{array}{cc}{g}_{11}& g_{12}\\ g_{21}& g_{22}\end{array}\right|>0,$$

(12)

and

$$\left|\begin{array}{ccc}{g}_{11}& g_{12}& g_{13}\\ g_{21}& g_{22}& g_{23}\\ g_{31}& g_{32}& g_{33}\end{array}\right|>0.$$

(13)

We tested frequently for stability.

5. Helmholtz Free Energy for Binary van der Waals

We take the Helmholtz free energy per mole for binary vdW from Konynenburg and Scott [19]:

$$\begin{array}{c}{\displaystyle A_m=e\left(T\right)-\mathcal{R}Tlog\left(\frac{V_m-b}{V_{m0}}\right)}\hfill \\ {\displaystyle +\mathcal{R}T\left[(1-x)log(1-x)+xlog\left(x\right)\right]-\frac{a}{V_m},}\hfill \end{array}$$

(14)

where

$$e\left(T\right)=-\frac{3}{2}\mathcal{R}Tln\left(\frac{T}{T_0}\right)+ϵ$$

(15)

is the purely thermal part of the ideal gas contribution, $\mathcal{R}$ is the universal gas constant, $V_m=V/(n_1+n_2)$ is the molar volume, $V_{m0}$, $T_0$, and $ϵ$ are constants that do not appear in R, and x is the concentration

$$x=\frac{n_2}{n_1+n_2}.$$

(16)

The quantities a and b are functions of x reflecting the purely attractive and repulsive parts of the interparticle interactions, respectively:

$$a\left(x\right)={(1-x)}^2{a}_{11}+2x(1-x)a_{12}+x^2{a}_{22},$$

(17)

and

$$b\left(x\right)=(1-x)b_{11}+x{b}_{22},$$

(18)

with $a_{11}$, $a_{12}$, $a_{22}$, $b_{11}$, and $b_{22}$ as the five independent vdW coefficients.

The full (extensive) Helmholtz free energy is

$$A(T,n_1,n_2,V)=(n_1+n_2)A_m,$$

(19)

and the Helmholtz free energy per volume is

$$f(T,{\rho }_1,{\rho }_2)=\frac{A(T,n_1,n_2,V)}{V}.$$

(20)

For the pure fluid, vdW simplifies by using scaled units [15]. For binary vdW, scaling cannot be done in a natural way. One possibility is to scale in terms of the van der Waals coefficients of one of the fluid components [19], but this seems somewhat artificial. In this paper, we avoid scaled units, and simply adopt real units when needed.

6. Repulsive Binary van der Waals

First consider the case with only repulsive interactions, with zero a coefficients, $b_{11}=b_1$, and $b_{22}=b_2$. By Equation (20), we have

$$\begin{array}{c}\hfill f(T,{\rho }_1,{\rho }_2)=-k_{B}T({\rho }_1+{\rho }_2)log\left[\frac{1-b_1{\rho }_1-b_2{\rho }_2}{V_{m0}}\right]+k_{B}T{\rho }_{1}log\left({\rho }_1\right)\\ \hfill +k_{B}T{\rho }_{2}log\left({\rho }_2\right)-\frac{3}{2}{k}_{B}T({\rho }_1+{\rho }_2)log\left(\frac{T}{T_0}\right)+ϵ({\rho }_1+{\rho }_2),\end{array}$$

(21)

where we send $\mathcal{R}\to k_B$ to convert ${\rho }_1$ and ${\rho }_2$ to units of particles per volume.

By using Equation (21), and the process described in Section 3, we find:

$$\begin{array}{c}\hfill R=\frac{1}{2\left({\rho }_1+{\rho }_2\right)}\left[1-2b_2{\rho }_1-2b_1{\rho }_2-4b_1{\rho }_1-4b_2{\rho }_2+5b_1^2{\rho }_1{\rho }_2\right.\\ \hfill \left.+5b_2^2{\rho }_1{\rho }_2+5b_1^2{\rho }_1^2+5b_2^2{\rho }_2^2\right].\end{array}$$

(22)

Note, R for repulsive vdW has no dependence on T, nor does it depend on the constants in Equation (21), other than the b’s.

Since R has units of volume per particle, $R\times ({\rho }_1+{\rho }_2)$ is dimensionless. Setting $b_1=b_2=b$ yields

$$R\times ({\rho }_1+{\rho }_2)=\frac{1}{2}\left[1-b\left({\rho }_1+{\rho }_2\right)\right]\left[1-5b\left({\rho }_1+{\rho }_2\right)\right].$$

(23)

If $b\to 0$, this expression agrees with Ruppeiner and Davis [4] for the binary ideal gas with the thermal contribution in Equation (15). $R\times ({\rho }_1+{\rho }_2)$ is shown as a contour graph in Figure 2.

With only the b coefficients nonzero, we would naively expect $R>0$, consistent with repulsive interactions. Cases with $R<0$ clearly occur, mirroring the situation in the pure fluid where, for example, the gas of hard-spheres has negative R [13]. Such anomalous results might be dismissed as aberrations since we are below the low $\left|R\right|$ limit. However, it is interesting that the negative sign persists from the pure fluid into the binary fluid. This consistency contributes to the need for a more nuanced interpretation for the sign of R.

7. Attractive Binary van der Waals

Now turn on the attractive interactions in binary vdW. In this scenario, Konynenburg and Scott [19] classified nine distinct possibilities, depending on the values of the vdW a and b coefficients. We consider here only an instance of their Type I possibility, with a single curve of critical points continuously connecting the critical points of the pure fluid components $x=0$ and $x=1$. Our example has vdW coefficients $a_{11}=0.002$, $a_{22}=0.005$, $a_{12}=0.004$, and $b_{11}=b_{22}=0.00002$. (In MKS units, the units of the a’s are Joules ${\mathrm{meters}}^3/{\mathrm{mole}}^2$, and the units of the b’s are ${\mathrm{meters}}^3/\mathrm{mole}$. These vdW coefficients produce critical points very roughly in the zone of normal fluid Helium.)

To structure the discussion, note that if T is high, the attraction between the particles has little effect, and we expect no phase transitions. Now consider lowering T slowly, with the particle numbers and the volume held fixed. We might move along such an isochore in a laboratory $PVT$ experiment; see Figure 3. As T decreases, attractive interatomic interactions become more effective, and the binary fluid could eventually break into two phases.

These phase transitions have associated second-order curves of critical points [19]. For Type 1, every x has a single critical point with critical coordinates $T=T_c$, $P=P_c$, and $V_m=V_{mc}$. In their Appendix A, Konynenburg and Scott [19] describe an explicit procedure $x\to (T_c,P_c,V_{mc})$ for locating these critical points.

The critical curve is shown in Figure 4. The $det\left(\mathbf{g}\right)$ switches signs upon crossing this critical curve, with the thermodynamics stable above the curve and unstable below it, according to Equations (11)–(13). R diverges to negative infinity on the critical curve.

The theoretical expectations on approaching the critical curve are well-known. A general physical argument was given by Widom [20] in the context of the hyperscaling critical exponent relation

$$d\nu =2-\alpha ,$$

(24)

where d is the spatial dimension ($d=3$ here), $\alpha$ is the heat capacity critical exponent ($\alpha =0$ here), and $\nu$ is the correlation length critical exponent. On approaching the critical curve from above along a critical isochore with fixed x, fixed $V_m=V_{mc}$, and decreasing T, the correlation length is expected to diverge according to:

$${\xi }^d\propto {\left(\frac{T-T_c}{T_c}\right)}^{\alpha -2}.$$

(25)

A key finding in previous research is that $\left|R\right|\propto {\xi }^d$ near a critical point [1,21,22,23], (The constant of proportionality between R and ${\xi }^3$ calculated for both fluid and magnetic systems in distinct spatial dimensions appears to be $-1/2$: $-\frac{1}{2}R={\xi }^d$; see [23] for a brief review.) and so we expect

$$\left|R\right|\propto {\left(\frac{T-T_c}{T_c}\right)}^{-2}$$

(26)

along a critical isochore.

An interesting contrast to binary vdW is offered by pure fluid vdW, where R is computed with 2D thermodynamic geometry. It was shown [1] that the asymptotic divergence of R for the pure vdW fluid, with the thermal ideal gas contribution as in Equation (15), along the critical isochore is

$$R=-b{\left(\frac{T-T_c}{T_c}\right)}^{-2},$$

(27)

independent of the vdW constant a. Here, $b=b_{11}$ if $x=0$ and $b=b_{22}$ if $x=1$.

Physically, ${\xi }^3$ should be the same for a binary fluid in the limits $x\to 0$ or $x\to 1$ as for the corresponding pure fluids. Since R in the critical region is a measure of ${\xi }^3$, we might then physically expect the limiting R 3D’s to be at least approximately equal to the corresponding R 2D’s. Remarkably, this limiting correspondence holds very well, even though the 2D and 3D formulas for R are quite different.

We computed R along four critical isochores in 3D thermodynamic geometry, $x=0.2,0.4,0.6$, and $0.8$. Results near the critical point are shown in Figure 5. The four curves clearly overlap each other very closely, and show very little dependence on x. The asymptotic R is always negative, in accordance with expectations for effectively attractive interactions. The critical curves have asymptotic power law divergences, in accordance with Equation (26). The critical exponents of the four curves are all within better than $0.3\%$ of the expected value 2, with the power law fits done over the full range of the data shown in Figure 5. Erdem [6] got similar 3D power law behavior in the Ising antiferromagnet. Note as well that the smallest values of $\left|R\right|$ in Figure 5 are about 4 nm${}^3$/atom, so values near the particle level are on the asymptotic power law line. Such a large span of the critical regime was seen also in magnetic systems [23].

Additionally shown in Figure 5 is $-R$ for the pure fluids $x=0$ and 1 calculated with the 2D thermodynamic geometry. Asymptotically, we expect these curves to obey Equation (27), and we found that they do. The R 2D curves are in excellent agreement with the corresponding R 3D curves.

Next, we consider letting the binary fluid expand into the binary ideal gas state, which has the dimensionless quantity [4]

$$R\times ({\rho }_1+{\rho }_2)=\frac{1}{2}.$$

(28)

We approached the binary ideal gas limit for four fixed values of x by starting at $(T_c,V_{mc})$ for each x, and increasing $V_m$ at constant $T=T_c$. Results of this expansion are shown in Figure 6. Clearly, for each x, $R\times ({\rho }_1+{\rho }_2)\to 1/2$, in accord with Equation (28).

8. Conclusions

We have made a start in calculating the full three-dimensional thermodynamic Ricci curvature scalar R in an interacting binary fluid system. Our main finding is that the physical interpretation of R for the pure fluid extends very naturally to the binary fluid. The emerging physical picture for the thermodynamic R is thus quite robust. We calculated R for two scenarios involving the van der Waals model for the binary fluid. The first had exclusively repulsive interactions, and the second added attractive interactions, and critical phenomena.

Let us place our results in the context of the expectations in Section 2. (a) In the binary ideal gas limit, we expect $R\times ({\rho }_1+{\rho }_2)=1/2$ from an earlier exact calculation. Figure 6 clearly shows this known limit. This limit is also shown at the origin of Figure 2. Both these figures (in the ideal gas limit) have $\left|R\right|$ smaller than the low $\left|R\right|$ limit, and so the results here are not as strong as those in the critical regime, where $\left|R\right|$ is much larger. (b) Near critical points, $\left|R\right|$ was found to diverge with the expected critical exponent of 2, as shown in Figure 5. In the limits as the binary fluid goes to the pure fluid, we found excellent concordance with the pure fluid results calculated with the 2D thermodynamic geometry. These critical point results are the strongest that we present in this paper. (c) Near critical points, R was found to be negative, as shown in Figure 5. This negative sign is expected when interactions are effectively attractive. (d) Evaluating R along coexistence curves was beyond the scope of our paper. (e) Our regimes of positive R were all more or less expected, with values all under the low $\left|R\right|$ limit. A more systematic search of the thermodynamic phase space for more interesting cases (above the low $\left|R\right|$ limit) was beyond the scope of our paper. (f) We found that the repulsive vdW was negative at large densities, in accord with the anomaly found by Brańka et al., for the hard sphere [13]; see Figure 2. That this anomaly translates from the 2D to the 3D thermodynamic geometry is interesting.