  A new method is presented for mesoscopic simulations of particle dispersions
in liquid crystal solvents. It allows efficient first-principle simulations of
the dispersions involving many particles with many-body interactions mediated
by the solvents. Demonstrations have been performed for the aggregation of
colloid dispersions in two-dimensional nematic and smectic-C* solvents
neglecting hydrodynamic effects, which will be taken into account in the near
future.
