Lubachevsky-Stillinger algorithm (LS algorithm, LSA,
or LS protocol) sometimes also called
Lubachevsky-Stillinger compression algorithm
or protocol) is
a numerical procedure that simulates or imitates
a physical process of compressing an assembly
of hard particles. As the LSA may need thousands of arithmetic operations even for a few particles,
it is usually carried out on a digital computer.
A real physical process of compression typically
involves a contracting container
boundary, such as a piston pressing against the particles. The LSA is able to simulate just
such a scenario, like in ???.
However,
in a more frequently used setting and in the
example reported in paper ???,
where the LSA was firstly
introduced,
the LSA compresses
the virtual particles by "swelling" or expanding them
in a fixed, final (but not necessarily bounded if a periodic
boundary condition is adopted as is done in ???)
virtual volume.
The absolute sizes of the particles are increasing but particle-to-particle relative sizes remain constant.
As a result, in a final, compressed, or "jammed" state,
some particles, the so-called "rattlers," turn out not to be jammed. Rattlers are mobile
within "cages" formed by their immobile, jammed neighbors
and the boundary, if any.
A substantial limitation of the original LS protocol
is that it was designed to practically work only
for spherical particles, though the spheres may be
of different sizes ??? ???. Any deviation from the spherical
(or circular in two dimensions) shape, even a simplest one, when spheres are replaced with ellipsoids (or ellipses in two dimensions) ??? ???, causes thus modified LSA to slow down dramatically.
But as long as the shape is spherical,
the LSA is able to handle particle ensembles
in tens to hundreds of thousands
on today's (2011) standard personal computers.
How useful the LSA is in dimensions higher than 3
is unknown.
References
- (LS 90): B.D. Lubachevsky and F.H. Stillinger,
Geometric Properties of Random Disk Packings, Journal of Statistical Physics 60, nos. 5/6 (1990) [1]