A variety of mesoscale simulation techniques has been developed over the last decades. The central idea of the methods like dissipative particle dynamics, smooth particle hydrodynamics, lattice Boltzmann, or multiparticle collision dynamics (stochastic rotation dynamics) is the replacement of complete set of molecules in a complex fluid by a smaller number of particles possessing the same properties as the initial set. In addition to reducing the number of degrees of freedom, this replacement allows one to avoid direct modelling of fast motions and, therefore, the need for small time steps in simulations. While the systematic coarse-graining methods based on rigorous statistical mechanics exist for static properties, the mesoscale methods for complex fluids often lose the connection to the original microscopic description. Different approaches exists to set up the parameters of coarse-grained system based on matching (i) equilibrium properties of the original fluid (density, pressure, size of hydration layer, surface tension, contact angle), (ii) dynamic properties (viscosity, friction coefficients, molecular relaxation times). The practice shows, however, that fixing the static properties does not guarantee the accuracy of the resulting dynamic behaviour.

We believe that there is a need of methodological advances in order to compute out of atomistic simulations the different reversible (effective potentials, equations of state, surface tension) and irreversible parts (frictions, viscosities, relaxation times) of the coarse-grained dynamic equations. The main obstacle in achieving this goal is the need of dealing with functions of many variables. Pair-wise assumptions seem to be the easiest solution, but may be there are other possibilities.

Attempts to bridge the timescales to microsecond and longer regimes have been made using various mesoscale methods, such as lattice-Boltzmann methods [1-4], particle based off-lattice methods such as dissipative particle dynamics [5-8], multi-particle-collision (or stochastic rotational) dynamics (MPCD) [9-12], and Brownian dynamics of stokeslets with hydrodynamics based on the Green function formalism [13-14]. On the basis of these techniques a few hybrid schemes were developed and used to study phenomenology of non-equilibrium processes, dynamic scaling, etc. The DPD approach has been successfully used for modelling surfactant self-assembly, including aggregation of phospholipids. The MPCD approach was applied for modelling dynamics of long polymer chains near obstacles, shear driven dynamics of star polymers, vesicles, and red blood cells, colloidal sedimentation [15]. The lattice Boltzmann method was applied to model flow in narrow channels, past rough surfaces, and surfaces with complex boundary conditions including partial slippage [16]. A few quantitative studies addressed dynamics of disperse systems: colloidal dispersions or polymer solutions where, however, only one or two types of particles with a single interaction potential were involved [15,17,18]. A recent themed issue in Phys.Chem.Chem.Phys  [19] offers a general view of the state of the art on Coarse-grained modeling of soft condensed matter. 

From a theoretical point of view, coarse-graining has a solid foundation in what is also known as non-equilibrium statistical mechanics with a clear treatment of both static and dynamic properties, the latter in the form of Green-Kubo expressions [20,21]. While the framework is well established, its practical implementation has only been successfully addressed for static properties. Just recently dynamic properties are being considered in a systematic way [22,23]