No Arabic abstract
Longitudinal and transverse sound velocities of Lennard-Jones systems are calculated at the liquid-solid coexistence using the additivity principle. The results are shown to agree well with the ``exact values obtained from their relations to excess energy and pressure. Some consequences, in particular, in the context of the Lindemanns melting rule and Stokes-Einstein relation between the self-diffusion and viscosity coefficients are discussed. Comparison with available experimental data on the sound velocities of solid argon at melting conditions is provided.
In recent years lines along which structure and dynamics are invariant to a good approximation, so-called isomorphs, have been identified in the thermodynamic phase diagrams of several model liquids and solids. This paper reports computer simulations of the transverse and longitudinal collective dynamics at different length scales along an isomorph of the Lennard-Jones system. Our findings are compared to corresponding results along an isotherm and an isochore. Confirming the theoretical prediction, the reduced-unit dynamics of the transverse momentum density is invariant to a good approximation along the isomorph at all time and length scales. Likewise, the wave-vector dependent shear-stress autocorrelation function is found to be isomorph invariant. A similar invariance is not seen along the isotherm or the isochore. Using a spatially non-local hydrodynamic model for the transverse momentum-density time-autocorrelation function, the macroscopic shear viscosity and its wave dependence are determined, demonstrating that the shear viscosity is isomorph invariant on all length scales studied. This analysis implies the existence of a novel length scale which characterizes each isomorph. The transverse sound-wave velocity, the Maxwell relaxation time, and the rigidity shear modulus are also isomorph invariant. In contrast, the reduced-unit dynamics of the mass density is not invariant at length scales longer than the inter-particle distance. By fitting to a generalized hydrodynamic model, we extract values for the wave-vector-dependent thermal diffusion coefficient, sound attenuation coefficient, and adiabatic sound velocity. The isomorph variation of these quantities in reduced units at long length scales can be eliminated by scaling with $gamma$, a fundamental quantity in the isomorph theory framework, an empirical observation that remains to be explained theoretically.
Determining the microstructure of colloidal suspensions under shear flows has been a challenge for theoretical and computational methods due to the singularly-perturbed boundary-layer nature of the problem. Previous approaches have been limited to the case of hard-sphere systems and suffer from various limitations in their applicability. We present a new analytic scheme based on intermediate asymptotics which solves the Smoluchowski diffusion-convection equation including both intermolecular and hydrodynamic interactions. The method is able to recover previous results for the hard-sphere fluid and, for the first time, to predict the radial distribution function (rdf) of attractive fluids such as the Lennard-Jones (LJ) fluid. In particular, a new depletion effect is predicted in the rdf of the LJ fluid under shear. This method can be used for the theoretical modelling and understanding of real fluids subjected to flow, with applications ranging from chemical systems to colloids, rheology, plasmas, and atmospherical science.
We calculate the density of states of a binary Lennard-Jones glass using a recently proposed Monte Carlo algorithm. Unlike traditional molecular simulation approaches, the algorithm samples distinct configurations according to self-consistent estimates of the density of states, thereby giving rise to uniform internal-energy histograms. The method is applied to simulate the equilibrium, low-temperature thermodynamic properties of a widely studied glass former consisting of a binary mixture of Lennard-Jones particles. We show how a density-of-states algorithm can be combined with particle identity swaps and configurational bias techniques to study that system. Results are presented for the energy and entropy below the mode coupling temperature.
Liquids displaying strong virial-potential energy correlations conform to an approximate density scaling of their structural and dynamical observables. This scaling property does not extend to the entire phase diagram, in general. The validity of the scaling can be quantified by a correlation coefficient. In this work a simple scheme to predict the correlation coefficient and the density-scaling exponent is presented. Although this scheme is exact only in the dilute gas regime or in high dimension d, a comparison with results from molecular dynamics simulations in d = 1 to 4 shows that it reproduces well the behavior of generalized Lennard-Jones systems in a large portion of the fluid phase.
Combining the recent Piskulich-Thompson approach [Z. A. Piskulich and W. H. Thompson, {it J. Chem. Phys.} {bf 152}, 011102 (2020)] with isomorph theory, from a single simulation, the structure of a single-component Lennard-Jones (LJ) system is obtained at an arbitrary state point in almost the whole liquid region of the temperature-density phase diagram. The LJ system exhibits two temperature range where the vant Hoffs assumption that energetic and entropic forces are temperature independent is valid. A method to evaluate the structure at an arbitrary state point along an isochore from the knowledge of structures at two temperatures on the isochore is also discussed. We argue that, in general, the structure of any hidden scale-invariant system obeying the vant Hoffs assumption in the whole range of temperatures can be determined in the whole liquid region of the phase diagram from only a single simulation.