Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userRadial distribution function of Lennard-Jones fluids in shear flows from intermediate asymptotics
75
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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.
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.
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.
It is well known from the quantum theory of strongly correlated systems that poles (or more subtle singularities) of dynamic correlation functions in complex plane usually correspond to the collective or localized modes. Here we address singularities of velocity autocorrelation function $Z$ in complex $omega$-plain for the one-component particle system with isotropic pair potential. We have found that naive few poles picture fails to describe analytical structure of $Z(omega)$ of Lennard-Jones particle system in complex plain. Instead of few isolated poles we see the singularity manifold of $Z(omega)$ forming branch cuts that suggests Lennard-Jones velocity autocorrelation function is a multiple-valued function of complex frequency. The brunch cuts are separated from the real axis by the well-defined gap. The gap edges extend approximately parallel to the real frequency axis. The singularity structure is very stable under increase of the temperature; we have found its trace at temperatures even several orders of magnitude higher than the melting point. Our working hypothesis that the branch cut origin is related to the interference in $Z$ of one-particle kinetics and collective hydrodynamic motion.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
