Lattice gas with nearest-neighbor exclusion in a shear-like field

There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.

Abstract

We present Monte Carlo simulations of the lattice gas with nearest-neighbor exclusion and Kawasaki (hopping) dynamics (hard square lattice gas), under the influence of a nonuniform drive, on the square lattice. The drive, which favors motion along the +x axis and inhibits motion in the opposite direction, varies linearly in the y direction. Our lattice has rigid walls at the end points in the y direction and periodic boundaries along the drive. We find that this model has transition to a sublattice-ordered phase at a density of about 0.298, lower than in equilibrium (rhoc <FONT FACE=Symbol>@</FONT> 0.37), but somewhat higher than in the uniformly driven case at maximal bias (rhoc <FONT FACE=Symbol>@</FONT> 0.272). For smaller global densities (r < 0.33), the ordering occurs with particle accumulation in the low-drive region. Above this density we observe a surprising reversal in the density profile, with particles migrating to the high-drive region.

Related collections

Most cited references 4

Record: found
Abstract: found
Article: found

Is Open Access

Why do ultrasoft repulsive particles cluster and crystallize? Analytical results from density functional theory

Bianca M Mladek, Dieter Gottwald, Gerhard Kahl … (2007)

We demonstrate the accuracy of the hypernetted chain closure and of the mean-field approximation for the calculation of the fluid-state properties of systems interacting by means of bounded and positive-definite pair potentials with oscillating Fourier transforms. Subsequently, we prove the validity of a bilinear, random-phase density functional for arbitrary inhomogeneous phases of the same systems. On the basis of this functional, we calculate analytically the freezing parameters of the latter. We demonstrate explicitly that the stable crystals feature a lattice constant that is independent of density and whose value is dictated by the position of the negative minimum of the Fourier transform of the pair potential. This property is equivalent with the existence of clusters, whose population scales proportionally to the density. We establish that regardless of the form of the interaction potential and of the location on the freezing line, all cluster crystals have a universal Lindemann ratio L = 0.189 at freezing. We further make an explicit link between the aforementioned density functional and the harmonic theory of crystals. This allows us to establish an equivalence between the emergence of clusters and the existence of negative Fourier components of the interaction potential. Finally, we make a connection between the class of models at hand and the system of infinite-dimensional hard spheres, when the limits of interaction steepness and space dimension are both taken to infinity in a particularly described fashion.

0 comments Cited 430 times – based on 0 reviews

Preprint

     Review now

Bookmark

Record: found
Abstract: found
Article: found

Is Open Access

Path Integral Molecular Dynamics within the Grand Canonical-like Adaptive Resolution Technique: Simulation of Liquid Water

, (2015)

Quantum effects due to the spatial delocalization of light atoms are treated in molecular simulation via the path integral technique. Among several methods, Path Integral (PI) Molecular Dynamics (MD) is nowadays a powerful tool to investigate properties induced by spatial delocalization of atoms; however computationally this technique is very demanding. The abovementioned limitation implies the restriction of PIMD applications to relatively small systems and short time scales. One possible solution to overcome size and time limitation is to introduce PIMD algorithms into the Adaptive Resolution Simulation Scheme (AdResS). AdResS requires a relatively small region treated at path integral level and embeds it into a large molecular reservoir consisting of generic spherical coarse grained molecules. It was previously shown that the realization of the idea above, at a simple level, produced reasonable results for toy systems or simple/test systems like liquid parahydrogen. Encouraged by previous results, in this paper we show the simulation of liquid water at room conditions where AdResS, in its latest and more accurate Grand-Canonical-like version (GC-AdResS), is merged with two of the most relevant PIMD techniques available in literature. The comparison of our results with those reported in literature and/or with those obtained from full PIMD simulations shows a highly satisfactory agreement.