Multi-particle collision dynamics

Last updated

Multi-particle collision dynamics (MPC), also known as stochastic rotation dynamics (SRD), [1] is a particle-based mesoscale simulation technique for complex fluids which fully incorporates thermal fluctuations and hydrodynamic interactions. [2] Coupling of embedded particles to the coarse-grained solvent is achieved through molecular dynamics. [3]

Contents

Method of simulation

The solvent is modelled as a set of point particles of mass with continuous coordinates and velocities . The simulation consists of streaming and collision steps.

During the streaming step, the coordinates of the particles are updated according to

where is a chosen simulation time step which is typically much larger than a molecular dynamics time step.

After the streaming step, interactions between the solvent particles are modelled in the collision step. The particles are sorted into collision cells with a lateral size . Particle velocities within each cell are updated according to the collision rule

where is the centre of mass velocity of the particles in the collision cell and is a rotation matrix. In two dimensions, performs a rotation by an angle or with probability . In three dimensions, the rotation is performed by an angle around a random rotation axis. The same rotation is applied for all particles within a given collision cell, but the direction (axis) of rotation is statistically independent both between all cells and for a given cell in time.

If the structure of the collision grid defined by the positions of the collision cells is fixed, Galilean invariance is violated. It is restored with the introduction of a random shift of the collision grid. [4]

Explicit expressions for the diffusion coefficient and viscosity derived based on Green-Kubo relations are in excellent agreement with simulations. [5] [6]

Simulation parameters

The set of parameters for the simulation of the solvent are:

The simulation parameters define the solvent properties, [1] such as

where is the dimensionality of the system.

A typical choice for normalisation is . To reproduce fluid-like behaviour, the remaining parameters may be fixed as . [7]

Applications

MPC has become a notable tool in the simulations of many soft-matter systems, including

Related Research Articles

<span class="mw-page-title-main">Hierarchy problem</span> Unsolved problem in physics

In theoretical physics, the hierarchy problem is the problem concerning the large discrepancy between aspects of the weak force and gravity. There is no scientific consensus on why, for example, the weak force is 1024 times stronger than gravity.

In physics, a Feshbach resonance can occur upon collision of two slow atoms, when they temporarily stick together forming an unstable compound with short lifetime. It is a feature of many-body systems in which a bound state is achieved if the coupling(s) between at least one internal degree of freedom and the reaction coordinates, which lead to dissociation, vanish. The opposite situation, when a bound state is not formed, is a shape resonance. It is named after Herman Feshbach, a physicist at MIT.

In non-equilibrium thermodynamics, GENERIC is an acronym for General Equation for Non-Equilibrium Reversible-Irreversible Coupling. It is the general form of dynamic equation for a system with both reversible and irreversible dynamics. GENERIC formalism is the theory built around the GENERIC equation, which has been proposed in its final form in 1997 by Miroslav Grmela and Hans Christian Öttinger.

In physical cosmology, cosmological perturbation theory is the theory by which the evolution of structure is understood in the Big Bang model. Cosmological perturbation theory may be broken into two categories: Newtonian or general relativistic. Each case uses its governing equations to compute gravitational and pressure forces which cause small perturbations to grow and eventually seed the formation of stars, quasars, galaxies and clusters. Both cases apply only to situations where the universe is predominantly homogeneous, such as during cosmic inflation and large parts of the Big Bang. The universe is believed to still be homogeneous enough that the theory is a good approximation on the largest scales, but on smaller scales more involved techniques, such as N-body simulations, must be used. When deciding whether to use general relativity for perturbation theory, note that Newtonian physics is only applicable in some cases such as for scales smaller than the Hubble horizon, where spacetime is sufficiently flat, and for which speeds are non-relativistic.

The Bose–Hubbard model gives a description of the physics of interacting spinless bosons on a lattice. It is closely related to the Hubbard model that originated in solid-state physics as an approximate description of superconducting systems and the motion of electrons between the atoms of a crystalline solid. The model was introduced by Gersch and Knollman in 1963 in the context of granular superconductors. The model rose to prominence in the 1980s after it was found to capture the essence of the superfluid-insulator transition in a way that was much more mathematically tractable than fermionic metal-insulator models.

<span class="mw-page-title-main">Quantum vortex</span> Quantized flux circulation of some physical quantity

In physics, a quantum vortex represents a quantized flux circulation of some physical quantity. In most cases, quantum vortices are a type of topological defect exhibited in superfluids and superconductors. The existence of quantum vortices was first predicted by Lars Onsager in 1949 in connection with superfluid helium. Onsager reasoned that quantisation of vorticity is a direct consequence of the existence of a superfluid order parameter as a spatially continuous wavefunction. Onsager also pointed out that quantum vortices describe the circulation of superfluid and conjectured that their excitations are responsible for superfluid phase transitions. These ideas of Onsager were further developed by Richard Feynman in 1955 and in 1957 were applied to describe the magnetic phase diagram of type-II superconductors by Alexei Alexeyevich Abrikosov. In 1935 Fritz London published a very closely related work on magnetic flux quantization in superconductors. London's fluxoid can also be viewed as a quantum vortex.

The Gross–Pitaevskii equation describes the ground state of a quantum system of identical bosons using the Hartree–Fock approximation and the pseudopotential interaction model.

Stokesian dynamics is a solution technique for the Langevin equation, which is the relevant form of Newton's 2nd law for a Brownian particle. The method treats the suspended particles in a discrete sense while the continuum approximation remains valid for the surrounding fluid, i.e., the suspended particles are generally assumed to be significantly larger than the molecules of the solvent. The particles then interact through hydrodynamic forces transmitted via the continuum fluid, and when the particle Reynolds number is small, these forces are determined through the linear Stokes equations. In addition, the method can also resolve non-hydrodynamic forces, such as Brownian forces, arising from the fluctuating motion of the fluid, and interparticle or external forces. Stokesian Dynamics can thus be applied to a variety of problems, including sedimentation, diffusion and rheology, and it aims to provide the same level of understanding for multiphase particulate systems as molecular dynamics does for statistical properties of matter. For rigid particles of radius suspended in an incompressible Newtonian fluid of viscosity and density , the motion of the fluid is governed by the Navier–Stokes equations, while the motion of the particles is described by the coupled equation of motion:

<span class="mw-page-title-main">Light front quantization</span> Technique in computational quantum field theory

The light-front quantization of quantum field theories provides a useful alternative to ordinary equal-time quantization. In particular, it can lead to a relativistic description of bound systems in terms of quantum-mechanical wave functions. The quantization is based on the choice of light-front coordinates, where plays the role of time and the corresponding spatial coordinate is . Here, is the ordinary time, is one Cartesian coordinate, and is the speed of light. The other two Cartesian coordinates, and , are untouched and often called transverse or perpendicular, denoted by symbols of the type . The choice of the frame of reference where the time and -axis are defined can be left unspecified in an exactly soluble relativistic theory, but in practical calculations some choices may be more suitable than others.

Microrheology is a technique used to measure the rheological properties of a medium, such as microviscosity, via the measurement of the trajectory of a flow tracer. It is a new way of doing rheology, traditionally done using a rheometer. There are two types of microrheology: passive microrheology and active microrheology. Passive microrheology uses inherent thermal energy to move the tracers, whereas active microrheology uses externally applied forces, such as from a magnetic field or an optical tweezer, to do so. Microrheology can be further differentiated into 1- and 2-particle methods.

In the context of the physical and mathematical theory of percolation, a percolation transition is characterized by a set of universal critical exponents, which describe the fractal properties of the percolating medium at large scales and sufficiently close to the transition. The exponents are universal in the sense that they only depend on the type of percolation model and on the space dimension. They are expected to not depend on microscopic details such as the lattice structure, or whether site or bond percolation is considered. This article deals with the critical exponents of random percolation.

<span class="mw-page-title-main">Kicked rotator</span>

The kicked rotator, also spelled as kicked rotor, is a paradigmatic model for both Hamiltonian chaos and quantum chaos. It describes a free rotating stick in an inhomogeneous "gravitation like" field that is periodically switched on in short pulses. The model is described by the Hamiltonian

<span class="mw-page-title-main">Metadynamics</span> Scientific computer simulation method

Metadynamics is a computer simulation method in computational physics, chemistry and biology. It is used to estimate the free energy and other state functions of a system, where ergodicity is hindered by the form of the system's energy landscape. It was first suggested by Alessandro Laio and Michele Parrinello in 2002 and is usually applied within molecular dynamics simulations. MTD closely resembles a number of newer methods such as adaptively biased molecular dynamics, adaptive reaction coordinate forces and local elevation umbrella sampling. More recently, both the original and well-tempered metadynamics were derived in the context of importance sampling and shown to be a special case of the adaptive biasing potential setting. MTD is related to the Wang–Landau sampling.

The Vicsek model is a mathematical model used to describe active matter. One motivation of the study of active matter by physicists is the rich phenomenology associated to this field. Collective motion and swarming are among the most studied phenomena. Within the huge number of models that have been developed to catch such behavior from a microscopic description, the most famous is the model introduced by Tamás Vicsek et al. in 1995.

<span class="mw-page-title-main">Two-body Dirac equations</span> Quantum field theory equations

In quantum field theory, and in the significant subfields of quantum electrodynamics (QED) and quantum chromodynamics (QCD), the two-body Dirac equations (TBDE) of constraint dynamics provide a three-dimensional yet manifestly covariant reformulation of the Bethe–Salpeter equation for two spin-1/2 particles. Such a reformulation is necessary since without it, as shown by Nakanishi, the Bethe–Salpeter equation possesses negative-norm solutions arising from the presence of an essentially relativistic degree of freedom, the relative time. These "ghost" states have spoiled the naive interpretation of the Bethe–Salpeter equation as a quantum mechanical wave equation. The two-body Dirac equations of constraint dynamics rectify this flaw. The forms of these equations can not only be derived from quantum field theory they can also be derived purely in the context of Dirac's constraint dynamics and relativistic mechanics and quantum mechanics. Their structures, unlike the more familiar two-body Dirac equation of Breit, which is a single equation, are that of two simultaneous quantum relativistic wave equations. A single two-body Dirac equation similar to the Breit equation can be derived from the TBDE. Unlike the Breit equation, it is manifestly covariant and free from the types of singularities that prevent a strictly nonperturbative treatment of the Breit equation.

<span class="mw-page-title-main">Interatomic potential</span> Functions for calculating potential energy

Interatomic potentials are mathematical functions to calculate the potential energy of a system of atoms with given positions in space. Interatomic potentials are widely used as the physical basis of molecular mechanics and molecular dynamics simulations in computational chemistry, computational physics and computational materials science to explain and predict materials properties. Examples of quantitative properties and qualitative phenomena that are explored with interatomic potentials include lattice parameters, surface energies, interfacial energies, adsorption, cohesion, thermal expansion, and elastic and plastic material behavior, as well as chemical reactions.

In solid state physics, the Luttinger–Ward functional, proposed by Joaquin Mazdak Luttinger and John Clive Ward in 1960, is a scalar functional of the bare electron-electron interaction and the renormalized one-particle propagator. In terms of Feynman diagrams, the Luttinger–Ward functional is the sum of all closed, bold, two-particle irreducible diagrams, i.e., all diagrams without particles going in or out that do not fall apart if one removes two propagator lines. It is usually written as or , where is the one-particle Green's function and is the bare interaction.

Spin squeezing is a quantum process that decreases the variance of one of the angular momentum components in an ensemble of particles with a spin. The quantum states obtained are called spin squeezed states. Such states have been proposed for quantum metrology, to allow a better precision for estimating a rotation angle than classical interferometers. Recently, it was shown that these states cannot provide a better precision.

A fracton is an emergent topological quasiparticle excitation which is immobile when in isolation. Many theoretical systems have been proposed in which fractons exist as elementary excitations. Such systems are known as fracton models. Fractons have been identified in various CSS codes as well as in symmetric tensor gauge theories.

References

  1. 1 2 Gompper, G.; Ihle, T.; Kroll, D. M.; Winkler, R. G. (2009). Multi-Particle Collision Dynamics: A Particle-Based Mesoscale Simulation Approach to the Hydrodynamics of Complex Fluids. pp. 1–87. arXiv: 0808.2157 . Bibcode:2009acsa.book....1G. doi:10.1007/978-3-540-87706-6_1. ISBN   978-3-540-87705-9. S2CID   8433369.{{cite book}}: |journal= ignored (help)
  2. Malevanets, Anatoly; Kapral, Raymond (1999). "Mesoscopic model for solvent dynamics". The Journal of Chemical Physics. 110 (17): 8605–8613. Bibcode:1999JChPh.110.8605M. doi:10.1063/1.478857.
  3. 1 2 Malevanets, Anatoly; Kapral, Raymond (2000). "Solute molecular dynamics in a mesoscale solvent". The Journal of Chemical Physics. 112 (16): 7260–7269. Bibcode:2000JChPh.112.7260M. doi: 10.1063/1.481289 . S2CID   73679245.
  4. Ihle, T.; Kroll, D. M. (2003). "Stochastic rotation dynamics. I. Formalism, Galilean invariance, and Green-Kubo relations". Physical Review E. 67 (6): 066705. Bibcode:2003PhRvE..67f6705I. doi:10.1103/PhysRevE.67.066705. PMID   16241378.
  5. Ihle, T.; Tüzel, E.; Kroll, D. M. (2004). "Resummed Green-Kubo relations for a fluctuating fluid-particle model". Physical Review E. 70 (3): 035701. arXiv: cond-mat/0404305 . Bibcode:2004PhRvE..70c5701I. doi:10.1103/PhysRevE.70.035701. PMID   15524580. S2CID   11272882.
  6. Ihle, T.; Tüzel, E.; Kroll, D. M. (2005). "Equilibrium calculation of transport coefficients for a fluid-particle model". Physical Review E. 72 (4): 046707. arXiv: cond-mat/0505434 . Bibcode:2005PhRvE..72d6707I. doi:10.1103/PhysRevE.72.046707. PMID   16383567. S2CID   14413944.
  7. 1 2 J. Elgeti "Sperm and Cilia Dynamics" PhD thesis, Universität zu Köln (2006)
  8. Padding, J. T.; Louis, A. A. (2004). "Hydrodynamic and Brownian Fluctuations in Sedimenting Suspensions". Physical Review Letters. 93 (22): 220601. arXiv: cond-mat/0409133 . Bibcode:2004PhRvL..93v0601P. doi:10.1103/PhysRevLett.93.220601. PMID   15601076. S2CID   119504730.
  9. Hecht, Martin; Harting, Jens; Bier, Markus; Reinshagen, Jörg; Herrmann, Hans J. (2006). "Shear viscosity of claylike colloids in computer simulations and experiments". Physical Review E. 74 (2): 021403. arXiv: cond-mat/0601413 . Bibcode:2006PhRvE..74b1403H. doi:10.1103/PhysRevE.74.021403. PMID   17025421. S2CID   19998245.
  10. Mussawisade, K.; Ripoll, M.; Winkler, R. G.; Gompper, G. (2005). "Dynamics of polymers in a particle-based mesoscopic solvent" (PDF). The Journal of Chemical Physics. 123 (14): 144905. Bibcode:2005JChPh.123n4905M. doi:10.1063/1.2041527. PMID   16238422.
  11. Ripoll, M.; Winkler, R. G.; Gompper, G. (2007). "Hydrodynamic screening of star polymers in shear flow". The European Physical Journal E. 23 (4): 349–354. Bibcode:2007EPJE...23..349R. doi:10.1140/epje/i2006-10220-0. PMID   17712520. S2CID   36780360.
  12. Noguchi, Hiroshi; Gompper, Gerhard (2005). "Dynamics of fluid vesicles in shear flow: Effect of membrane viscosity and thermal fluctuations" (PDF). Physical Review E. 72 (1): 011901. Bibcode:2005PhRvE..72a1901N. doi:10.1103/PhysRevE.72.011901. PMID   16089995.
  13. K.-W. Lee and Marco G. Mazza (2015). "Stochastic rotation dynamics for nematic liquid crystals". Journal of Chemical Physics. 142 (16): 164110. arXiv: 1502.03293 . Bibcode:2015JChPh.142p4110L. doi:10.1063/1.4919310. PMID   25933755. S2CID   36839435.