Next Article in Journal
The Infuence of Salicin on Rheological and Film-Forming Properties of Collagen
Next Article in Special Issue
Isomorph Invariance of Higher-Order Structural Measures in Four Lennard–Jones Systems
Previous Article in Journal
Obesity and Energy Substrate Transporters in Ovarian Cancer—Review
Previous Article in Special Issue
Phase-Transition Thermal Charging of a Channel-Shape Thermal Energy Storage Unit: Taguchi Optimization Approach and Copper Foam Inserts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sound Velocities of Generalized Lennard-Jones (n − 6) Fluids Near Freezing

1
Bauman Moscow State Technical University, 105005 Moscow, Russia
2
Joint Institute for High Temperatures, Russian Academy of Sciences, 125412 Moscow, Russia
Molecules 2021, 26(6), 1660; https://doi.org/10.3390/molecules26061660
Submission received: 15 February 2021 / Revised: 12 March 2021 / Accepted: 15 March 2021 / Published: 16 March 2021
(This article belongs to the Special Issue Thermodynamics and Transport Properties of Fluids)

Abstract

:
In a recent paper [S. Khrapak, Molecules 25, 3498 (2020)], the longitudinal and transverse sound velocities of a conventional Lennard–Jones system at the liquid–solid coexistence were calculated. It was shown that the sound velocities remain almost invariant along the liquid–solid coexistence boundary lines and that their magnitudes are comparable with those of repulsive soft-sphere and hard-sphere models at the fluid–solid phase transition. This implies that attraction does not considerably affect the magnitude of the sound velocities at the fluid–solid phase transition. This paper provides further evidence to this by examining the generalized Lennard–Jones (n − 6) fluids with n ranging from 12 to 7 and demonstrating that the steepness of the repulsive term has only a minor effect on the magnitude of the sound velocities. Nevertheless, these minor trends are identified and discussed.

1. Introduction

Sound velocities are important characteristics of a substance. They are directly related to long-wavelength excitations—phonons, which play a crucial role in condensed matter, materials science, and soft matter. Therefore, it is important to understand mechanisms that can affect and regulate sound velocities in various situations.
In simple atomic systems with smooth pairwise interaction potentials, the sound velocities can be expressed as sums over atoms involving the first and second derivatives of the interaction potential (see below). In isotropic situations (in gases and liquids), these sums can also be expressed as integrals involving the radial distribution function g ( r ) . The sound velocities are thus completely determined by the shape of the interaction potential and the atomic structure of the system (which are, of course, interrelated).
There have been indications that the shape of the interatomic interaction potential may not be an important factor for the sound velocities in special regions of the phase diagram. For example, the ratio of the sound to thermal velocity of many liquid metals and metalloids has about the same value 10 at the melting temperature [1,2,3]. A similar value (≃9.5) was reported from experiments with solid argon at the melting temperature [4]. Values close to ≃9 were reported for solid hydrogen and deuterium along the melting curve [5]. Rosenfeld pointed out that this “quasi-universal” property is also shared by the hard-sphere (HS) model [3]. More recently, it was demonstrated that this property is also exhibited by the purely repulsive soft inverse-power-law (IPL) model in a wide range of IPL exponents [6].
A generalization of the celebrated Lindemann’s criterion of melting of classical solids was proposed recently on the basis of statistical mechanics arguments [7]. With this generalized formulation, the expressions of the melting temperature are equivalent in three dimensions (3D) and two-dimensions (2D). Moreover, independently of dimensionality, the melting condition predicts that the ratio of the transverse sound velocity to the thermal velocity reaches a quasi-universal value along the melting curve (though the magnitudes can be different in 3D and 2D) [7,8].
Motivated by these observations, it was recently investigated in detail how the longitudinal and transverse sound velocities behave at the liquid–solid phase transition of 3D Lennard–Jones (LJ) systems [9]. It was observed that the sound velocities maintain quasi-universal values along the liquid–solid coexistence boundaries and that their magnitudes are comparable with those of repulsive soft-sphere and hard-sphere models. Hence, it was concluded that attraction does not affect the magnitude of sound velocities at the fluid–solid phase transition. This paper provides further evidence for this conclusion by examining the generalized Lennard–Jones (n − 6) fluids with n ranging from 12 to 7 to vary the relative strength of attraction. It will be documented that the steepness of the repulsive term has only a relatively weak effect on the magnitudes of the longitudinal and transverse sound velocities.
In principle, one could directly vary the attractive strength by considering a (12 − n) LJ system. However, since the long-range r 6 attraction has a physical justification (dispersion forces), while the r 12 repulsion is just a convenient mathematical approximation (among other possibilities) for a strong repulsion at short distances, the ( n 6 ) systems are much better investigated. Some of the available results will be used below in this paper.

2. Materials and Methods

In this work, the generalized (n − 6) LJ potential is considered:
ϕ ( r ) = C n ϵ ( σ / r ) n ( σ / r ) 6 ,
where
C n = n n 6 n 6 6 n 6 ,
so that for n = 12 , we get C 12 = 4 . Here, ϵ and σ are the energy and length scales and the exponent n determines the steepness of the repulsive term in the potential (1). The reduced pressure, energy, density, and temperature are expressed as P * = P σ 3 / ϵ , u * = U / N ϵ , ρ * = ρ σ 3 , and T * = T / ϵ , respectively.
For a spherically symmetric pairwise interaction potential ϕ ( r ) , the instantaneous longitudinal and transverse velocities can be expressed (in 3D) as follows [10,11,12]
m c l 2 = 1 30 j 2 r j ϕ ( r j ) + 3 r j 2 ϕ ( r j ) ,
m c t 2 = 1 30 j 4 r j ϕ ( r j ) + r j 2 ϕ ( r j ) .
Here, c l and c t are the instantaneous longitudinal and transverse sound velocities, m is the particle mass, the sums run over all neighbours of a test particle, and primes denote derivatives of the interaction potential with respect to the distance r. The averaging denotes that we perform summation over all N test particles in the system and then divide the sum by N. In isotropic fluids, such an averaging is often replaced by an integral involving the radial distribution function [10,13,14]: j ( ) 4 π ρ ( ) r 2 g ( r ) d r , where g ( r ) is the radial distribution function. For further details regarding the evaluation of sound velocities and their relation to the instantaneous elastic moduli, see Ref. [9].
For a pairwise interaction potential, the excess (over the ideal gas) contributions to the energy, u ex , and pressure, p ex , can be expressed via summations similar to those used above (we omit averaging to keep the notation compact):
u ex = 1 2 T j ϕ ( r j ) , p ex = 1 6 T j r j ϕ ( r j ) ,
where the reduced units are used: u ex = U ex / N T , p ex = P ex / ρ T (the temperature T is measured in energy units, k B = 1 ) . It is obvious then that for the ( n 6 ) LJ potential, we can express the reduced sound velocities in terms of reduced excess energy and pressure. After a simple and straightforward algebra, we obtain
c l 2 / v T 2 = A l ( n ) u ex + B l ( n ) p ex ,
c t 2 / v T 2 = A t ( n ) u ex + B t ( n ) p ex ,
with the n-dependent numerical coefficients
A l ( n ) = 6 n 5 , B l ( n ) = n ( 3 n + 1 ) 114 5 ( n 6 ) ,
A t ( n ) = 2 n 5 , B t ( n ) = n ( n 3 ) 18 5 ( n 6 ) .
The sound velocities in Equations (5) and (6) are expressed in units of thermal velocity, v T = T / m . Equations (5)–(8) can be considered as a generalization of the Zwanzig and Mountain result for the conventional (12-6) LJ potential. In this special case, the constants are A l ( 12 ) = 72 / 5 , B l ( 12 ) = 11 , A t ( 12 ) = 24 / 5 , and B t ( 12 ) = 3 [9,13].
The translation to LJ units is straightforward by recognizing that
P * P σ 3 ϵ = ρ T ( 1 + p ex ) σ 3 ϵ = ρ * T * ( 1 + p ex ) ,
u * U N ϵ = N T ( 3 / 2 + u ex ) N ϵ = T * ( 3 / 2 + u ex ) .
In the present paper, we focus on the fluid boundary of the fluid–solid coexistence (freezing line). In the crystalline state, although the isotropic treatment does not hold anymore, the angle-averaged sound velocities can be introduced in essentially the same way as in fluids. Calculated in this way, the sound velocities of the (12-6) LJ system at the solid coexistence boundary are numerically close to those at the fluid boundary, except for in the vicinity of the triple point (see Figure 4 from Ref. [9]).

3. Results and Discussion

The sound velocities of generalized (n − 6) LJ fluids under freezing conditions have been evaluated using the freezing parameters ( T * , P * , and ρ * ) tabulated by Sousa et al. [15] (simple analytical fits for the freezing and melting curves are also available [16], but are not used here). Note that the reduced energy tabulated there corresponds to the potential contribution only, u * = T * u ex . The results are presented in Figure 1.
In Figure 1a, the reduced sound velocities are plotted versus the reduced temperature. The main observations are as follows. For each repulsive exponent n, the longitudinal and transverse sound velocities remain practically constant along the freezing curve. The approximate constancy of the transverse sound velocity at the fluid–solid phase transition can be interpreted by employing the Lindemann’s melting criterion arguments [17]. In one of the formulations of this famous criterion, the melting condition can be expressed in terms of the instantaneous elastic longitudinal and shear moduli and thus in terms of the longitudinal and transverse sound velocities [7,18]. In particular, the condition contains a sum over inverse longitudinal and shear moduli, each multiplied by the number of the corresponding collective modes. Since the shear modulus is smaller than the longitudinal one and there are two shear and one longitudinal mode in 3D systems, the contribution from shear modes dominates (this is known as shear dominance), thus leading to the conclusion that c t / v T should be approximately constant at the fluid–solid phase transition [7]. The approximate constancy of the longitudinal sound velocity can be explained as follows. Although the effective softness of a given LJ ( n 6 ) potential somewhat varies along the freezing curve [19] (i.e., with T * ), the variation of c l / v T with effective softness is rather weak (provided the potential is not too soft and long-ranged), see Figure 1 from Ref. [9]. Thus, the reduced longitudinal sound velocity remains approximately constant for interactions that are not too soft.
Regarding the dependence of the reduced sound velocities on the repulsive exponent n at a fixed T * , a very weak, but systematic tendency of increasing c l / v T and c t / v T with n is observed. This tendency is consistent with a weak increase in the reduced sound velocities when the HS interaction limit is approached, which has been documented in Figure 1 from Ref. [9]. Increasing n leads to somewhat steeper repulsion and the reduced sound velocities increase accordingly, although very weakly.
To a reasonable approximation, the reduced sound velocities can be approximated by constant values of c l / v T 11 and c t / v T 6 . Deviations from actual values are within ∼10% for all exponents n in the range of T * considered. The sound velocity of purely repulsive soft inverse-power-law interaction potentials fall within the same range. For instance, the sound velocities of the IPL-12 system at freezing are c l / v T 11.7 and c t / v T 5.8 [9]. In this sense, it is legitimate to say that the presence of long-range attraction has a rather minor effect on the sound velocities of simple fluids near the fluid–solid phase transition.
The ratio of the longitudinal-to-transverse sound velocities is plotted in Figure 1b. It exhibits a slow monotonous increase with T * . In addition, at higher temperatures, there is an observable decline in c l / c t with increasing n. This tendency should be expected. For example, for the potential of soft Yukawa (screened Coulomb) interactions, the ratio c l / c t increases when the interaction becomes softer (screening length increases) and diverges on approaching the one-component plasma limit (see e.g., Figure 2 from Ref. [20]). At the same time, for all n considered here and the range of T * investigated, the ratio c l / c t remains in a relatively narrow range, between ≃1.7 and ≃2.0.

4. Conclusions

The instantaneous longitudinal and transverse sound velocities of the generalized ( n 6 ) Lennard–Jones fluids have been evaluated at freezing conditions for integers n between 7 and 12. For this purpose, the relation between the sound velocities (or elastic moduli) and excess energy and pressure derived by Zwanzig and Mountain [13] for the conventional (12-6) LJ potential has been generalized to the case of ( n 6 ) potential. The required freezing densities, pressures, and internal energies for temperatures in the range from T * = 1 to T * = 5 have been taken from Ref. [15].
The calculated ratios of sound velocities to the thermal velocity are practically constant along the freezing curves with the typical values c l / v T 11 and c t / v T 6 (note that these constant reduced values can result in quite diverse dimensional values when the temperature and the substance atomic mass are varied). These ratios are comparable with reduced sound velocities of repulsive soft sphere (IPL) and HS fluids at freezing. Thus, roughly speaking neither the presence of the long-range attraction nor the variation of the exponent in the repulsive term do not lead to considerable variations in the sound velocities. At the same time, minor systematic trends have been identified and discussed. These trends can be convincingly explained on the basis of the known (weak) dependence of the reduced sound velocities at freezing on the effective softness of the interaction potential.

Funding

Work on the effect of attraction on collective modes in fluids was supported by the Russian Science Foundation, Grant No. 20-12-00356.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

This study does not generate any new data.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Iida, T.; Guthrie, R. The Physical Properties of Liquid Metals; Oxford University Press: New York, NY, USA, 1988. [Google Scholar]
  2. Blairs, S. Sound velocity of liquid metals and metalloids at the melting temperature. Phys. Chem. Liq. 2007, 45, 399–407. [Google Scholar] [CrossRef]
  3. Rosenfeld, Y. Sound velocity in liquid metals and the hard-sphere model. J. Phys. Condens. Matter 1999, 11, L71–L74. [Google Scholar] [CrossRef]
  4. Ishizaki, K.; Spain, I.L.; Bolsaitis, P. Measurements of longitudinal and transverse ultrasonic wave velocities in compressed solidified argon and their relationship to melting theory. J. Chem. Phys. 1975, 63, 1401–1410. [Google Scholar] [CrossRef]
  5. Liebenberg, D.H.; Mills, R.L.; Bronson, J.C. Measurement of P, V, T, and sound velocity across the melting curve of n-H2 and n-D2 to 19 kbar. Phys. Rev. B 1978, 18, 4526–4532. [Google Scholar] [CrossRef]
  6. Khrapak, S.A. Note: Sound velocity of a soft sphere model near the fluid-solid phase transition. J. Chem. Phys. 2016, 144, 126101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Khrapak, S.A. Lindemann melting criterion in two dimensions. Phys. Rev. Res. 2020, 2, 012040. [Google Scholar] [CrossRef] [Green Version]
  8. Khrapak, S. Note: Melting criterion for soft particle systems in two dimensions. J. Chem. Phys. 2018, 148, 146101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Khrapak, S.A. Sound Velocities of Lennard-Jones Systems Near the Liquid-Solid Phase Transition. Molecules 2020, 25, 3498. [Google Scholar] [CrossRef] [PubMed]
  10. Schofield, P. Wavelength-dependent fluctuations in classical fluids: I. The long wavelength limit. Proc. Phys. Soc. 1966, 88, 149–170. [Google Scholar] [CrossRef]
  11. Balucani, U.; Zoppi, M. Dynamics of the Liquid State; Clarendon Press: Oxford, UK, 1994. [Google Scholar]
  12. Takeno, S.; Gôda, M. A Theory of Phonons in Amorphous Solids and Its Implications to Collective Motion in Simple Liquids. Prog. Theor. Phys. 1971, 45, 331–352. [Google Scholar] [CrossRef] [Green Version]
  13. Zwanzig, R.; Mountain, R.D. High-Frequency Elastic Moduli of Simple Fluids. J. Chem. Phys. 1965, 43, 4464–4471. [Google Scholar] [CrossRef]
  14. Hubbard, J.; Beeby, J. Collective motion in liquids. J. Phys. C 1969, 2, 556–571. [Google Scholar] [CrossRef]
  15. Sousa, J.M.G.; Ferreira, A.L.; Barroso, M.A. Determination of the solid-fluid coexistence of the n - 6 Lennard-Jones system from free energy calculations. J. Chem. Phys. 2012, 136, 174502. [Google Scholar] [CrossRef] [PubMed]
  16. Khrapak, S.A.; Ning, N. Freezing and melting equations for the n-6 Lennard-Jones systems. AIP Adv. 2016, 6, 055215. [Google Scholar] [CrossRef] [Green Version]
  17. Lindemann, F. The calculation of molecular vibration frequencies. Z. Phys. 1910, 11, 609. [Google Scholar]
  18. Buchenau, U.; Zorn, R.; Ramos, M.A. Probing cooperative liquid dynamics with the mean square displacement. Phys. Rev. E 2014, 90, 042312. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Khrapak, S.A.; Chaudhuri, M.; Morfill, G.E. Communication: Universality of the melting curves for a wide range of interaction potentials. J. Chem. Phys. 2011, 134, 241101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Khrapak, S.A. Unified description of sound velocities in strongly coupled Yukawa systems of different spatial dimensionality. Phys. Plasmas 2019, 26, 103703. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Longitudinal ( c l ) and transverse ( c t ) sound velocities of generalized LJ n 6 fluids at freezing vs the reduced temperature T * . In (a) the sound velocities are expressed in units of thermal velocity v T ; Upper curves are for the longitudinal mode, lower curves are for the transverse mode. In (b), the ratio of the longitudinal-to-transverse sound velocities, c l / c t , is plotted.
Figure 1. Longitudinal ( c l ) and transverse ( c t ) sound velocities of generalized LJ n 6 fluids at freezing vs the reduced temperature T * . In (a) the sound velocities are expressed in units of thermal velocity v T ; Upper curves are for the longitudinal mode, lower curves are for the transverse mode. In (b), the ratio of the longitudinal-to-transverse sound velocities, c l / c t , is plotted.
Molecules 26 01660 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khrapak, S. Sound Velocities of Generalized Lennard-Jones (n − 6) Fluids Near Freezing. Molecules 2021, 26, 1660. https://doi.org/10.3390/molecules26061660

AMA Style

Khrapak S. Sound Velocities of Generalized Lennard-Jones (n − 6) Fluids Near Freezing. Molecules. 2021; 26(6):1660. https://doi.org/10.3390/molecules26061660

Chicago/Turabian Style

Khrapak, Sergey. 2021. "Sound Velocities of Generalized Lennard-Jones (n − 6) Fluids Near Freezing" Molecules 26, no. 6: 1660. https://doi.org/10.3390/molecules26061660

APA Style

Khrapak, S. (2021). Sound Velocities of Generalized Lennard-Jones (n − 6) Fluids Near Freezing. Molecules, 26(6), 1660. https://doi.org/10.3390/molecules26061660

Article Metrics

Back to TopTop