Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDensity scaling of generalized Lennard-Jones fluids in different dimensions
95
0
0.0
(
0
)
Added by
Lorenzo Costigliola Dr
Publication date
2019
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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.
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.
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.
This paper investigates the relation between the density-scaling exponent $gamma$ and the virial potential-energy correlation coefficient $R$ at several thermodynamic state points in three dimensions for the generalized $(2n,n)$ Lennard-Jones (LJ) system for $n=4, 9, 12, 18$, as well as for the standard $n=6$ LJ system in two, three, and four dimensions. The state points studied include many low-density states at which the virial potential-energy correlations are not strong. For these state points we find the roughly linear relation $gammacong 3nR/d$ in $d$ dimensions. This result is discussed in light of the approximate extended inverse power law description of generalized LJ potentials [N. P. Bailey et al., J. Chem. Phys. 129, 184508 (2008)]. In the plot of $gamma$ versus $R$ there is in all cases a transition around $Rapprox 0.9$, above which $gamma$ starts to decrease as $R$ approaches unity. This is consistent with the fact that $gammarightarrow 2n/d$ for $Rrightarrow 1$, a limit that is approached at high densities and/or temperatures at which the repulsive $r^{-2n}$ term dominates the physics.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
