Revised Enskog theory and molecular dynamics simulations of the viscosities and thermal conductivity of the hard-sphere fluid and crystal.

Hard-sphere (HS) shear, longitudinal, cross, and bulk viscosities and the thermal conductivity are obtained by molecular dynamics (MD) simulations, covering the entire density range from the dilute fluid to the solid crystal near close-packing. The transport coefficient data for the HS crystal are largely new and display, unlike for the fluid, a surprisingly simple behavior in that they can be represented well by a simple function of the density compressibility factor. In contrast to the other four transport coefficients (which diverge), the bulk viscosity in the solid is quite small and decreases rapidly with increasing density, tending to zero in the close-packed limit. The so-called cross viscosity exhibits a different behavior to the other viscosities, in being negative over the entire solid range, and changes sign from negative to positive on increasing the density in the fluid phase. The extent to which the viscosity tensor and thermal conductivity of the HS crystal can be represented by revised Enskog theory (RET) is investigated. The RET expressions are sums of an instantaneous (I), a kinetic (K), and a so-called α part. The I part of the transport coefficients evaluated directly by MD are statistically indistinguishable from those of the corresponding kinetic theory (Enskog and RET) expressions. For the K part the integral over the spatial two-particle distribution function at contact was determined and the α part was estimated using the direct correlation function and density functional theory approximations. All three parts were determined in this work which allowed the accuracy of RET for solid systems to be assessed rigorously. It is found that in the case of the thermal conductivity the predictions of RET are in excellent agreement with the MD results. Also, for the shear viscosity the agreement over the entire solid phase is quite good but is considerably worse for the three remaining viscosities in the solid phase.

Determining the Kerr constant in optically isotropic liquid crystals.

A new scheme is investigated for evaluating the temperature dependence and dispersion relation of the Kerr constant (K) of an optically isotropic medium in isotropic and blue phases (BPs) liquid crystals. The scheme employs the measurement of the component of the transmitted light intensity of double modulated frequency using the modified in-plane-switching cell geometry (based on metallic film electrodes). It overcomes to a large extent the problem of a nonuniform electric field, employs relatively small driving voltages, and allows K to be measured directly. It is shown that the dispersion relation based on the single-band birefringence model describes well both blue and isotropic liquid crystal phases. It is found that the experimental data indicate that the temperature-dependent coefficients in this relation have a simple linear form in the isotropic phase, which allows a general model for the temperature and wavelength dependence of the Kerr constant in the isotropic liquid crystal phase to be formulated. In the BPs the temperature dependence of the experimental data deviate from the simple linear trend, but follow well an inverse exponential form.

Structural properties of additive binary hard-sphere mixtures. III. Direct correlation functions.

An analysis of the direct correlation functions c_{ij}(r) of binary additive hard-sphere mixtures of diameters σ_{s} and σ_{b} (where the subscripts s and b refer to the "small" and "big" spheres, respectively), as obtained with the rational-function approximation method and the WM scheme introduced in previous work [S. Pieprzyk et al., Phys. Rev. E 101, 012117 (2020)2470-004510.1103/PhysRevE.101.012117], is performed. The results indicate that the functions c_{ss}(r<σ_{s}) and c_{bb}(r<σ_{b}) in both approaches are monotonic and can be well represented by a low-order polynomial, while the function c_{sb}(r<1/2(σ_{b}+σ_{s})) is not monotonic and exhibits a well-defined minimum near r=1/2(σ_{b}-σ_{s}), whose properties are studied in detail. Additionally, we show that the second derivative c_{sb}^{''}(r) presents a jump discontinuity at r=1/2(σ_{b}-σ_{s}) whose magnitude satisfies the same relationship with the contact values of the radial distribution function as in the Percus-Yevick theory.

Structural properties of additive binary hard-sphere mixtures. II. Asymptotic behavior and structural crossovers.

The structural properties of additive binary hard-sphere mixtures are addressed as a follow-up of a previous paper [S. Pieprzyk et al., Phys. Rev. E 101, 012117 (2020)]2470-004510.1103/PhysRevE.101.012117. The so-called rational-function approximation method and an approach combining accurate molecular dynamics simulation data, the pole structure representation of the total correlation functions, and the Ornstein-Zernike equation are considered. The density, composition, and size-ratio dependencies of the leading poles of the Fourier transforms of the total correlation functions h_{ij}(r) of such mixtures are presented, those poles accounting for the asymptotic decay of h_{ij}(r) for large r. Structural crossovers, in which the asymptotic wavelength of the oscillations of the total correlation functions changes discontinuously, are investigated. The behavior of the structural crossover lines as the size ratio and densities of the two species are changed is also discussed.

A comprehensive study of the thermal conductivity of the hard sphere fluid and solid by molecular dynamics simulation.

This work reports a new set of hard sphere (HS) thermal conductivity coefficient, λ, data obtained by Molecular Dynamics (MD) computer simulation, over a density range covering the dilute fluid to near the close-packed solid, and for a large number of particles (up to N = 13 1072) and long simulation times. The N-dependence of the thermal conductivity is shown to be proportional to N-2/3 to a good approximation over a wide range of system sizes, which enabled λ values in the thermodynamic limit to be predicted accurately. The fluid and solid λ can be represented well by the Enskog theory (ET) formula, λE, times a density-dependent correction term, which is close to unity for the fluid and practically constant for the solid. The convergence of the MD λ data back towards ET in the metastable fluid starts just above the freezing density. For the HS solid and dense fluid it was found that the thermal conductivity is nearly linear in pressure, as has been observed experimentally for a number of solids. Simple excess entropy scaling over the higher density fluid phase region was found, and Rosenfeld's exponential relationship can be fitted to the simulation data for the solid to a high degree of accuracy. The simulation analysis has revealed a number of new trends in the behaviour of the HS thermal conductivity which could be useful in building more accurate models for heat conduction in experimental systems.

Thermodynamic and dynamical properties of the hard sphere system revisited by molecular dynamics simulation.

Revised thermodynamic and dynamical properties of the hard sphere (HS) system are obtained from extensive molecular dynamics calculations carried out with large system sizes (number of particles, N) and long times. Accurate formulas for the compressibility factor of the HS solid and fluid branches are proposed, which represent the metastable region and take into account its divergence at close packing. Some basic second-order thermodynamic properties are obtained and a maximum in some of their derivatives in the metastable fluid region is found. The thermodynamic parameters associated with the melting-freezing transition have been determined to four digit accuracy, which generates accurate new values for the coexistence properties of the HS system. For the self-diffusion coefficient, D, it is shown that relatively large systems (N > 104) are required to achieve an accurate linear extrapolation of D to the infinite size limit with a D vs. N-1/3 plot. Moreover, it is found that there is a density dependence of the value of the slope in the linear regime. The density dependent correction becomes practically insignificant at higher densities and the hydrodynamic formula found in the literature is still accurate. However, with decreasing density the density dependence of the size correction cannot be neglected, which indicates that other sources of N-dependence, apart from those derived on purely hydrodynamic grounds, may also be important (and as yet unaccounted for). A detailed analytic representation of the density dependence of the HS self-diffusion coefficient and the HS viscosity, η, is given. It is shown that the HS viscosity near freezing and in the metastable region can be described well by the Krieger-Dougherty equation. Both D and η start to scale at high densities and in the metastable region in such a way that Dηp = const, where p ≃ 0.97, and D → 0 and η → ∞ at a packing fraction of 0.58, which coincides with some previous predictions of the HS glass transition density.

First derivative of the hard-sphere radial distribution function at contact.

Molecular dynamics simulations have been carried out of the radial distribution function of the hard sphere fluid for a range of densities in the equilibrium fluid and just into the metastable region. The first derivative of the hard-sphere radial distribution function at contact was computed and its density dependence fitted to a simple analytic form. Comparisons were made with semi-empirical formulae from the literature, and of these the formula proposed by Tao et al (1992 Phys. Rev. A 46 8007) was found to be in best agreement with the simulation data, although it slightly underestimates the derivative at the higher packing fractions in excess of about 0.45. Close to contact, within a few per cent of the particle diameter, the radial distribution function can be represented well by a second order polynomial. An exponential function, which has some useful analytic features, can also be applied in this region.