Thermodynamic Properties of the Parabolic-Well Fluid

Abstract

The thermodynamic properties of the parabolic-well fluid are considered. The intermolecular interaction potential of this model, which belongs to the class of the so-called van Hove potentials, shares with the square-well and the triangular well potentials the inclusion of a hard-core and an attractive well of relatively short range. The analytic second virial coefficient for this fluid is computed explicitly and an equation of state is derived with the aid of the second-order thermodynamic perturbation theory in the macroscopic compressibility approximation and taking the hard-sphere fluid as the reference system. For this latter, the fully analytical expression of the radial distribution function, consistent with the Carnahan-Starling equation of state as derived within the rational function approximation method, is employed. The results for the reduced pressure of the parabolic-well fluid as a function of the packing fraction and two values of the range of the parabolic-well potential at different temperatures are compared with Monte Carlo and Event‐driven molecular dynamics simulation data. Estimates of the values of the critical temperature are also provided.

1 Introduction

The issue of Frontiers in Physics this paper belongs to is devoted to commemorating the celebration of fifty consecutive annual Winter Meetings in Statistical Physics in Mexico. Therefore, we have chosen to write on a subject that has been present in these meetings from the beginning; namely, the thermodynamic properties of fluids that we are persuaded can still offer some interesting results.

We begin by recalling that, in an attempt to prove the validity of the thermodynamic limit of classical statistical mechanics, van Hove [1] introduced in 1949 a potential consisting of a hard core of radius and a finite-range attractive tail. The actual form of this so-called van Hove potential iswhere r is the distance, b corresponds to the range, and corresponds to the lower bound of the attractive tail, whose form is rather arbitrary. It should be pointed out that two popular models of intermolecular potentials used in liquid state physics, namely, the triangle-well potential and the square-well potential, fulfill the condition of being van Hove potentials and their thermodynamic properties have been thoroughly studied (see, for instance, Refs. [2–11] for the former model and Refs. [12–24] for the latter and references therein). Surprisingly, as far as we know, the parabolic-well potential, which is also a van Hove potential, has not been used for that purpose. The main aim of this paper is to contribute to partly remedying this situation.

We consider a parabolic-well fluid whose molecules interact with a potential of the formwhere is the reduced distance (r being the distance), σ is the diameter of the hard core, is the well depth, and is the potential range. As it occurs with other relatively simple models, the main asset of this model potential is probably that, despite being an idealized representation, it nevertheless contains the main features of true molecular interactions in fluids, namely, a repulsive hard-core and an attractive interaction that continuously goes to zero as the intermolecular distance increases. In this regard, it is interesting to recall what Widom [25] pointed out in the case of the square-well fluid: “Where I speak of the necessity to treat accurately the effects of the attractive or repulsive forces, I do not mean that it is important to know the corresponding part of with quantitative accuracy. Indeed, even if the of Figure 1 were idealized as a square-well potential, as in Figure 3, but the statistical mechanical consequences of such a potential were then determined without further approximation, there would undoubtedly result in an essentially correct description of all the macroscopic properties of matter throughout a vast region of the p and T plane, including the neighborhoods of the triple and critical points. Thus, what matters is not the quantitative accuracy of the assumed , but rather the qualitative accuracy of the resulting spatial correlations of molecular positions; the triple and critical points are distinguished by having the relevant qualitative features of this correlation, and the nature of its propagation through the fluid, determined primarily by the short-range repulsive forces between molecules, or by the longer ranged attractive forces, respectively.” Something similar may be said about the parabolic-well potential. In fact, an interesting asset of this model is that its thermodynamic properties are readily amenable for treatment within the second-order thermodynamic perturbation theory of Barker and Henderson [26]. Within this approach, in order to derive the Helmholtz free energy of the parabolic-well fluid, two ingredients are required: on the one hand, one needs the Helmholtz free energy of the reference hard-sphere fluid of diameter σ. On the other hand, one also requires an expression for the radial distribution function of the hard-sphere fluid. In this work, we will profit from the availability of a method [27], the so-called rational function approximation (RFA) method, to (analytically) obtain an approximate which is thermodynamically consistent with the equation of state of the hard-sphere fluid and make use of this fact to derive the equation of state of the parabolic-well fluid.

FIGURE 1

The paper is organized as follows. In the next section, we recall the main aspects of the RFA method for the computation of and provide the explicit expression for this quantity in the first coordination shell. This is followed in Section 3 with the completely analytic derivation of the equation of state of the parabolic-well fluid within the second-order Barker-Henderson thermodynamic perturbation theory in the macroscopic compressibility approximation and taking the hard-sphere fluid as the reference system. Section 4 contains some illustrative results for the reduced pressure of the parabolic-well fluid and a comparison with our own Monte Carlo and Event-driven Molecular Dynamics simulation data. We close the paper in the final section with further discussion and some concluding remarks.

2 The Rfa Method For The Computation of The Radial Distribution Function of The Hard-Sphere Fluid

In this section, we provide the analytic result for the radial distribution function (rdf) of the hard-sphere fluid , as derived with the RFA method [27], and its explicit expression in the range . We begin by recalling two important relationships between the thermodynamic and structural properties of the hard-sphere fluid derived from statistical mechanics. On the one hand, the compressibility factor (where p is the pressure, ρ the number density, the Boltzmann constant, and T the absolute temperature) of the hard-sphere fluid is related to the contact value of the rdf throughwhere is the packing fraction. On the other hand, the hard-sphere isothermal susceptibility is related to the rdf through

In the RFA method [27], the Laplace transform of is taken to be given bywhereand the six coefficients , , , , , and (which depend on the packing fraction) may be evaluated in an algebraic form by imposing the following requirements: (i) must be finite and hence the first two integral moments of the total correlation function , i.e., with , must be well defined; (ii) the approximation must be thermodynamically consistent with a prescribed equation of state; i.e., the thermodynamic relationship must be satisfied. Using the first requirement, one finds that , , , and are linear functions of and . Imposing the requirement (ii) leads to explicit expressions for and in terms of and [27]. Finally, the expressions for all the coefficients are as follows:

Here, and are the compressibility factor and isothermal susceptibility arising in the Percus-Yevick theory. To close the problem, one has to give an expression for , so all the procedure is a function of this choice. For a given , the radial distribution function is given bywith being the Heaviside step function and

Explicitly, using the residues theorem,wherewith being the four roots of ; namely,where

As we will indicate below, once has been chosen, Eqs. 5–26 are all that is needed to evaluate the first- and second-order perturbation terms for the free energy of the parabolic-well fluid within the Barker-Henderson thermodynamic perturbation theory taking the hard-sphere fluid as the reference system. To close this section and for later use, we now write the explicit expression for the radial distribution function up to the first coordination shell which readswhere

3 Thermodynamic Perturbation Theory and the Equation of State of the Parabolic-Well Fluid

Perturbation approaches for the computation of thermodynamic properties of fluids are well established theoretical tools [28, 29]. In the Barker-Henderson perturbation theory [26], one splits the potential into a hard-sphere part and a perturbation part; namely, , whereand

Once this separation has been made, the Helmholtz free energy per particle of the parabolic-well fluid is expressed as a power series in the inverse of the reduced temperature , which up to second order reads

Here, N is the number of particles and stands for the Helmholtz free energy of the reference hard-sphere fluid while and (this latter in the so-called macroscopic compressibility approximation) are given, respectively, byand

Note that we have made use of the fact that vanishes for and of the expression for given in Eq. 30, respectively, to set the lower and upper limits of the integrals in Eqs. 32, 33. In turn, the equation of state of the parabolic-well fluid in this approximation is given byAnd, the chemical potential may be readily obtained aswhere is given in Eq. 31, together with Eqs. 32, 33, and Z is given in Eq. 34. So, provided we choose , which of course also determines and , and take to be the one computed with the RFA approach and such compressibility factor, the completely analytic formulation of the second-order Barker-Henderson thermodynamic perturbation theory in the macroscopic compressibility approximation for the parabolic-well fluid taking the hard-sphere fluid as the reference system has been derived. In our subsequent calculations, we will be restricted to relatively narrow wells () so that Eq. 27 for will be used. Furthermore, and will be chosen to be those corresponding to the Carnahan-Starling (CS) equation of state [30]; namely,and

Further, from the CS equation of state, it also follows that

The availability of the completely analytic (albeit approximate) forms of the Helmholtz free energy and the equation of state of the system (which are themselves not very illuminating and therefore will not be explicitly written down [31]) allows us in principle to compute, for a given value of λ, the compressibility factor using Eq. 34, the vapor-liquid coexistence curve from the equality of pressures, and chemical potentials of the two phases and also to obtain the critical point in the usual way.

Preliminary results for the isotherms will be presented in the following section, together with a comparison with our simulation data. But before presenting such results, we will take advantage of the simple form of the intermolecular potential of this fluid to compute its second virial coefficient. This is given by

Equation 39, which to the best of our knowledge has not been reported before, allows us to obtain the Boyle temperature of the parabolic-well fluid as a function of λ by equating to zero and solving numerically for . In Table 1, we show some particular values and, for comparison, we also include the values corresponding to triangle-well and square-well fluids with the same range λ.

To close this section, we will also take advantage of the knowledge of the second virial coefficient, to obtain estimates of the critical temperature according to the Vliegenthart and Lekkerkerker criterion [32], namely, from equating this coefficient with , where is the volume of the spherical core. The results for given values of the range are given in Table 2, where we have also included such estimates for the cases of the triangle-well and square-well fluids with the same range λ.

Note that, for all three model fluids, the values of both the reduced Boyle temperatures and the estimates of the reduced critical temperatures increase as the range λ is increased. Also note that the geometrical form of the well influences such values as reflected in the fact that, for the same value of the range, the ones corresponding to the triangle-well fluid are smaller than those of the parabolic-well fluid which, in turn, are smaller than those of the square-well fluid.

4 Illustrative Results

Now we return to our main aim. In order to assess the value of the thermodynamic perturbation theory approach presented in the previous section, we have carried out NVT Monte Carlo (MC) simulations to compute the pressure of parabolic-well fluids for various values of the range and supercritical temperatures for later comparison with our theoretical results. The details of such simulations are as follows. The number of particles in our simulations is and we have considered a cubic box, of length L and with periodic boundary conditions. Reduced units are used, so that lengths are expressed in units of σ (, with l a pure number), the reduced temperature is , the packing fraction is (where is the volume), and the reduced pressure is .

For the sake of illustration, we report here the results of the simulations for and 1.75, along various isotherms. For each isotherm, eight different packing fraction values to 0.5 with were simulated in order to compute the reduced pressure . Each run was carried out using Monte Carlo steps (MCS) discarding the first MCS for equilibration, and the properties were measured every 20 MCS and averaged every 1000 MCS; furthermore, for each packing fraction, the values of were averaged over 20 parallel simulations to obtain better statistics.

Finally, the pressure was calculated using the expressionHere, , is the radial distribution function of the parabolic-well fluid (computed in the usual way [33] with the subscript PW standing for parabolic well), and the second term on the right-hand side of Eq. 40, obtained following a similar procedure to the one used by Rotenberg [12] in the case of the square-well fluid, accounts for the hard-core contribution to the parabolic-well potential.

In Figures 1‐3, we show the comparison between the results of the isotherms obtained with the thermodynamic perturbation theory and from simulation. Note the good agreement between theoretical and simulation results for all the values of above the critical temperature that we considered.

On the other hand, the subcritical isotherms were obtained by means of Molecular Dynamics (MD) Event-driven simulations. We have performed event-driven simulations of elastic smooth spheres carried out with the DynamO software package [34]. The spheres interact by a stepped parabolic-well type potential [35–37], a discretized version consisting of a sequence of 15 steps of widths , steps more than reasonable in most instances [37]. We have used σ, , m, and as units of length, time, mass, and temperature, respectively. The MD event-driven simulations were performed for one value of the range of the potential, namely, , along the isotherms and 0.7. In the first stage, we have performed NVT simulations with an Andersen thermostat during collisions, and after the equilibration, the second stage of NVE simulations with a duration of collisions was performed [38]. The full pressure tensor for the system was determined by means of the expressionwhere the kinetic pressure is given bywith being a dyadic product which yields a matrix result and the masses of the particles are set as . The contribution to the pressure due to interactions is given bywhere the summation is over each two-particle event interaction, is the total simulation time, is the momentum impulse on particle i, and is the separation vector between the interacting particles. Finally, the hydrostatic pressure, which in this instance coincides with the reduced pressure, was computed from the trace of the tensor

In Figure 2, we show the comparison of the subcritical isotherms obtained with the thermodynamic perturbation theory and from simulation for a range . The typical van der Waals loop is clearly seen for the theoretical isotherm with (which grossly underestimates the simulation data) and is still present in the isotherms with and , respectively. On the other hand, we have checked both through Monte Carlo and Event-driven MD simulations and also through the outcome of the thermodynamic perturbation theory that the isotherm with (not shown), which according to the Vliegenthart and Lekkerkerker criterion should be the critical one, is a supercritical isotherm. In fact, the simulation data indicate that the real critical isotherm for this value of the range lies above but close to the one corresponding to the theoretical curve for .

FIGURE 2

While it is clear from Figures 1–3 that the qualitative trends observed in all the simulation results are correctly accounted for by the theory, a better perspective of its performance may be gained by looking at the quantitative differences. Therefore, in Table 3, we display the actual numerical values for a couple of isotherms. In both cases, it is clear that the good qualitative agreement seen in Figures 1, 2, respectively, is not accompanied by quantitative agreement. In fact, the first theoretical isotherm ( and ), which is a supercritical isotherm, yields an underestimation of the reduced pressure when compared to the simulation values. On the other hand, for the second isotherm ( and ), which is subcritical, the general overall trend is that the theoretical curve overestimates the value of the reduced pressure. As one would expect, in the case of the supercritical isotherms, the quantitative agreement is improved as the reduced temperature is increased.

FIGURE 3

5 Concluding Remarks

In this paper, we have addressed the study of the thermodynamic properties of a fluid whose molecules interact through a parabolic-well potential. For this model, we obtained the exact second virial coefficient which in turn allowed us to compute the Boyle temperature and to estimate the critical temperature for arbitrary values of the potential range λ. The parabolic-well potential is in the same family as the triangle-well potential and the square-well potential, being in some sense intermediate between the other two. A reflection of this is the behavior of both the Boyle temperatures and the estimates of the critical temperatures in which, for a fixed range, the values corresponding to the parabolic-well potential lie between the ones corresponding to the other two. Whether this points out to a deeper relationship between the geometrical shape of the well and the location of the critical point in van Hove fluids is not clear to us at this stage but might be worth considering in the future.

In order to obtain further analytic results, we considered a thermodynamic perturbation theory approach for this fluid within the Barker-Henderson second-order macroscopic compressibility approximation and taking the hard-sphere fluid as the reference fluid. Restricting ourselves to values of the range in the interval and evaluating the radial distribution function of the hard-sphere fluid according to the RFA method with the CS equation of state, we were able to derive (albeit approximate) fully analytic expressions for the Helmholtz free energy, the equation of state, and the chemical potential of the parabolic-well fluid. With such expressions, we were able to compute theoretically various isotherms for a given potential range. These were subsequently compared to our own Monte Carlo NVT and Event-driven MD simulation results. It must be emphasized that these simulation data are to our knowledge the only ones available in the literature for this system.

It should be clear that the calculations that we have presented in the previous section are still preliminary but we want to stress that further work on this subject is currently being carried out. Nevertheless, at this stage, a few additional comments are in order. We begin by pointing out that the qualitative agreement between the results for the isotherms above the critical one obtained from thermodynamic perturbation theory and those stemming out of NVT Monte Carlo simulations, as well as the improvement of the quantitative agreement as the reduced temperature is increased, although clearly rewarding, are not very surprising in view of the fact that our theoretical approximation relies on the convergence of the perturbation expansion for high temperatures. Also rewarding is the fact that the results of the Event-driven MD simulation for the isotherm in the case in which the range is , which is a subcritical isotherm, are also well accounted for by the curve obtained using thermodynamic perturbation theory. The same happens with the isotherm with . On the other hand, the gross underestimation of the theoretical curve for the subcritical isotherm with and the same value of the range indicates that the convergence of the perturbation series is very poor for this reduced temperature. In any case, it is fair to say that the present theoretical approach provides a good starting point for the study of the thermodynamic properties of parabolic-well fluids. Future work with the same approach contemplates the computation of the critical point and the liquid-vapor coexistence curve of such fluids. Finally, since we are persuaded that the parabolic-well fluid may still offer some other insights on the thermodynamic behavior of fluids, we hope that the results of the present paper may also motivate others to conduct more studies using this model.

References

• 1.

  van HoveL. Quelques propriétés générales de L'intégrale de configuration D'un système de particules avec interaction. Physica (1949) 15:951–61. 10.1016/0031-8914(49)90059-2

  • CrossRef
  • Google Scholar

• 2.

  LargoJSolanaJR. A simplified perturbation theory for equilibrium properties of triangular-well fluids. Phys Stat Mech Appl (2000) 284:68–78. 10.1016/s0378-4371(00)00232-6

  • CrossRef
  • Google Scholar

• 3.

  Betancourt-CárdenasFFGalicia-LunaLASandlerSI. Thermodynamic properties for the triangular-well fluid. Mol Phys (2007) 105:2987–98. 10.1080/00268970701725013

  • CrossRef
  • Google Scholar

• 4.

  Betancourt-CárdenasFFGalicia-LunaLABenavidesALRamírezJASchöll-PaschingerE. Thermodynamics of a long-range triangle-well fluid. Mol Phys (2008) 106:113–26. 10.1080/00268970701832397

  • CrossRef
  • Google Scholar

• 5.

  ZhouS. Thermodynamics and phase behavior of a triangle-well model and density-dependent variety. J Chem Phys (2009) 130:014502–12. 10.1063/1.3049399

  • CrossRef
  • Google Scholar

• 6.

  KoyuncuM. Equation of state of a long-range triangular-well fluid. Mol Phys (2011) 109:565–73. 10.1080/00268976.2010.538738

  • CrossRef
  • Google Scholar

• 7.

  GuérinH. Improved analytical thermodynamic properties of the triangular-well fluid from perturbation theory. J Mol Liq (2012) 170:37–40. 10.1016/j.molliq.2012.03.014

  • CrossRef
  • Google Scholar

• 8.

  RiveraLDRoblesMLópez de HaroM. Equation of state and liquid–vapour equilibrium in a triangle-well fluid. Mol Phys (2012) 110:1327–33. 10.1080/00268976.2012.655338

  • CrossRef
  • Google Scholar

• 9.

  BárcenasMOdriozolaGOreaP. Coexistence and interfacial properties of triangle-well fluids. Mol Phys (2014) 112:2114–21. 10.1080/00268976.2014.887801

  • CrossRef
  • Google Scholar

• 10.

  TrejosVMMartínezAValadez-PérezNE. Statistical fluid theory for systems of variable range interacting via triangular-well pair potential. J Mol Liq (2018) 265:337–46. 10.1016/j.molliq.2018.05.116

  • CrossRef
  • Google Scholar

• 11.

  BenavidesALCervantesLATorres-ArenasJ. Analytical equations of state for triangle-well and triangle-shoulder potentials. J Mol Liq (2018) 271:670–6. 10.1016/j.molliq.2018.08.110

  • CrossRef
  • Google Scholar

• 12.

  RotenbergA. Monte Carlo equation of state for hard spheres in an attractive square well. J Chem Phys (1965) 43:1198–201. 10.1063/1.1696904

  • CrossRef
  • Google Scholar

• 13.

  BarkerJAHendersonD. Perturbation theory and equation of state for fluids: the square‐well potential. J Chem Phys (1967) 47:2856–61. 10.1063/1.1712308

  • CrossRef
  • Google Scholar

• 14.

  LuksKDKozakJJ. Adv Chem Phys (1978) 37:139–201.

  • Google Scholar

• 15.

  CarleyDD. Thermodynamic properties of a square‐well fluid in the liquid and vapor regions. J Chem Phys (1983) 78:5776–81. 10.1063/1.445462

  • CrossRef
  • Google Scholar

• 16.

  del RíoFLiraL. Properties of the square-well fluid of variable width. Mol Phys (1987) 61:275–92. 10.1080/00268978700101141

  • CrossRef
  • Google Scholar

• 17.

  del RíoFLiraL. Properties of the square‐well fluid of variable width. II. The mean field term. J Chem Phys (1987) 87:7179–83. 10.1063/1.453361

  • CrossRef
  • Google Scholar

• 18

  BenavidesALdel RíoF. Properties of the square-well fluid of variable width. Mol Phys (1989) 68:983–1000. 10.1080/00268978900102691

  • CrossRef
  • Google Scholar

• 19.

  López-RendónRReyesYOreaP. Thermodynamic properties of short-range square well fluid. J Chem Phys (2006) 125:084508–5. 10.1063/1.2338307

  • CrossRef
  • Google Scholar

• 20.

  Rivera-TorresSdel RíoFEspíndola-HerediaRKolafaJMalijevskýA. Molecular dynamics simulation of the free-energy expansion of the square-well fluid of short ranges. J Mol Liq (2013) 185:44–9. 10.1016/j.molliq.2012.12.005

  • CrossRef
  • Google Scholar

• 21.

  ElliotJRSchultzAJKofkeDA. Combined temperature and density series for fluid‐phase properties. I. Square-well spheres. J Chem Phys (2015) 147:1141101–12. 10.1063/1.4930268

  • CrossRef
  • Google Scholar

• 22.

  PadillaLBenavidesAL. The constant force continuous molecular dynamics for potentials with multiple discontinuities. J Chem Phys (2017) 147: 034502–6. 10.1063/1.4993436

  • CrossRef
  • Google Scholar

• 23.

  SastreFMoreno-HilarioESotelo-SernaMGGil-VillegasA. Microcanonical-ensemble computer simulation of the high-temperature expansion coefficients of the Helmholtz free energy of a square-well fluid. Mol Phys (2018) 116:351–60. 10.1080/00268976.2017.1392051

  • CrossRef
  • Google Scholar

• 24.

  RíoFd.GuzmánOMartínezFO. Global square-well free-energy model via singular value decomposition. Mol Phys (2018) 116:2070–82. 10.1080/00268976.2018.1461943

  • CrossRef
  • Google Scholar

• 25.

  WidomB. Intermolecular Forces and the Nature of the Liquid State: liquids reflect in their bulk properties the attractions and repulsions of their constituent molecules. Science (1967) 157:375–82. 10.1126/science.157.3787.375

  • CrossRef
  • Google Scholar

• 26.

  BarkerJAHendersonD. What is “liquid”? Understanding the states of matter. Rev Mod Phys (1976) 48:587–671. 10.1103/revmodphys.48.587

  • CrossRef
  • Google Scholar

• 27.

  López de HaroMYusteSBSantosA. Alternative approaches to the equilibrium properties of hard-sphere liquids. In: MuleroA, editor. Theory and simulation of hard-sphere fluids and related systems, lecture notes in physics 753. Berlin, Germany: Springer (2008). p. 183–245.

  • Google Scholar

• 28.

  ZhouSSolanaJR. Progress in the perturbation approach in fluid and fluid-related theories. Chem Rev (2009) 109:2829–58. 10.1021/cr900094p

  • CrossRef
  • Google Scholar

• 29.

  SolanaJR. Perturbation theories for the thermodynamic properties of fluids and solids. Boca Raton, Florida: CRC Press (2013).

  • Google Scholar

• 30.

  CarnahanNFStarlingKE. Equation of state for nonattracting rigid spheres. J Chem Phys (1969) 51:635–6. 10.1063/1.1672048

  • CrossRef
  • Google Scholar

• 31.

  These expressions are available, however, in a Mathematica code that we have employed and which we are willing to supply if requested.

  • Google Scholar

• 32.

  VliegenthartGALekkerkerkerHNW. Predicting the gas-liquid critical point from the second virial coefficient. J Chem Phys (2000) 112:5364–9. 10.1063/1.481106

  • CrossRef
  • Google Scholar

• 33.

  FrenkelDSmitB. Understanding molecular simulation: from algorithms and applications. San Francisco, CA: Academic press (2002).

  • Google Scholar

• 34.

  BannermanMNSargantRLueL. DynamO: a free ${\cal O}$(N) general event-driven molecular dynamics simulator. J Comput Chem (2011) 32:3329–38. 10.1002/jcc.21915

  • CrossRef
  • Google Scholar

• 35.

  ChapelaGAScrivenLEDavisHT. Molecular dynamics for discontinuous potential. IV. Lennard‐Jonesium. J Chem Phys (1989) 91:4307–13. 10.1063/1.456811

  • CrossRef
  • Google Scholar

• 36.

  ThomsonCLueLBannermanMN. Mapping continuous potentials to discrete forms. J Chem Phys (2014) 140:034105–13. 10.1063/1.4861669

  • CrossRef
  • Google Scholar

• 37.

  López de HaroMRodríguez-RivasAYusteSBSantosA. Structural properties of the jagla fluid. Phys Rev E (2018) 98:012138–48. 10.1103/physreve.98.012138

  • CrossRef
  • Google Scholar

• 38.

  BannermanMNLueLWoodcockLV. J Chem Phys (2010) 132:084607–12. 10.1063/1.3328823

  • CrossRef
  • Google Scholar