How would you integrate the equations of motion in dissipative particle dynamics simulations?

In this work we assess the quality and performance of several novel dissipative particle dynamics integration schemes that have not previously been tested independently. Based on a thorough comparison we identify the respective methods of Lowe and Shardlow as particularly promising candidates for future studies of large-scale properties of soft matter systems

The Chiral Phase Transition in Dissipative Dynamics

Numerical simulations of the chiral phase transition in the (3+1)dimensional O(4)-model are presented. The evolutions of the chiral field follow purely dissipative dynamics, starting from random chirally symmetric initial configurations down to the true vacuum with spontaneously broken symmetry. The model stabilizes topological textures which are formed together with domains of disoriented chiral condensate (DCC) during the roll-down phase. The classically evolving field acts as source for the emission of pions and $\sigma$ mesons. The exponents of power laws for the growth of angular correlations and for emission rates are extracted. Fluctuations in the abundance ratios for neutral and charged pions are compared with those for uncorrelated sources as potential signature for the chiral phase transition after heavy-ion collisions. It is found that the presence of stabilizing textures (baryons and antibaryons) prevents sufficiently rapid growth of DCC-domain size, so observability of anomalous tails in the abundance ratios is unlikely. However, the transient formation of growing DCC domains causes sizable broadening of the distributions as compared to the statistical widths of generic sources.Comment: 28 pages, 8 figure