We introduce new Langevin-type equations describing the rotational and translational motion of rigid bodies interacting through conservative and non-conservative forces and hydrodynamic coupling. In the absence of non-conservative forces, the Langevin-type equations sample from the canonical ensemble. The rotational degrees of freedom are described using quaternions, the lengths of which are exactly preserved by the stochastic dynamics. For the proposed Langevin-type equations, we construct a weak 2nd order geometric integrator that preserves the main geometric features of the continuous dynamics. The integrator uses Verlet-type splitting for the deterministic part of Langevin equations appropriately combined with an exactly integrated Ornstein-Uhlenbeck process. Numerical experiments are presented to illustrate both the new Langevin model and the numerical method for it, as well as to demonstrate how inertia and the coupling of rotational and translational motion can introduce qualitatively distinct behaviours.
REFERENCES
1.
I.
Snook
, The Langevin and Generalised Langevin Approach to the Dynamics of Atomic, Polymeric and Colloidal Systems
(Elsevier
, 2006
).2.
T. E.
Ouldridge
, “The importance of thermodynamics for molecular systems, and the importance of molecular systems for thermodynamics
,” Nat. Comput.
(published online).3.
D.
Frenkel
and B.
Smit
, Understanding Molecular Simulation
, 2nd ed. (Academic Press
, New York
, 2002
).4.
G. N.
Milstein
and M. V.
Tretyakov
, Phys. D
229
, 81
(2007
).5.
N. G.
Van Kampen
, Stochastic Processes in Physics and Chemistry
(Elsevier
, Amsterdam
, 2007
).6.
B. H.
Zimm
, J. Chem. Phys.
24
, 269
(1956
).7.
J. D.
Park
, J. S.
Myung
, and K. H.
Ahn
, Korean J. Chem. Eng.
33
, 3069
(2016
).8.
G. K.
Batchelor
, J. Fluid Mech.
52
, 245
(1971
).9.
M.
Długosz
and J.
Trylska
, BMC Biophys.
4
, 3
(2011
).10.
J.
Skolnick
, J. Chem. Phys.
145
, 100901
(2016
).11.
A.
Najafi
and R.
Golestanian
, Phys. Rev. E
69
, 062901
(2004
).12.
C. M.
Polley
, G. P.
Alexander
, and J. M.
Yeomans
, Phys. Rev. Lett.
99
, 228103
(2007
).13.
D.
Ermak
and J.
McCammon
, J. Chem. Phys.
69
, 1352
(1978
).14.
E.
Dickinson
, S. A.
Allison
, and J. A.
McCammon
, J. Chem. Soc., Faraday Trans.
81
, 591
(1985
).15.
J. F.
Brady
and G.
Bossis
, Annu. Rev. Fluid Mech.
20
, 111
(1988
).16.
D. J.
Evans
, Mol. Phys.
34
, 317
(1977
).17.
R. L.
Davidchack
, R.
Handel
, and M. V.
Tretyakov
, J. Chem. Phys.
130
, 234101
(2009
).18.
I. M.
Ilie
, W. J.
Briels
, and W. K.
den Otter
, J. Chem. Phys.
142
, 114103
(2015
).19.
R. L.
Davidchack
, T. E.
Ouldridge
, and M. V.
Tretyakov
, J. Chem. Phys.
142
, 144114
(2015
).20.
T. F.
Miller
III, M.
Eleftheriou
, P.
Pattnaik
, A.
Ndirango
, D.
Newns
, and G. J.
Martyna
, J. Chem. Phys.
116
, 8649
(2002
).21.
B.
Ai
and F.
Li
, Soft Matter
13
, 2536
(2017
).22.
R.
Kutteh
, J. Chem. Phys.
132
, 174107
(2010
).23.
R. Z.
Hasminskii
, Stochastic Stability of Differential Equations
(Sijthoff & Noordhoff
, 1980
).24.
C.
Soize
, The Fokker-Planck Equation for Stochastic Dynamical Systems and Its Explicit Steady State Solutions
(World Scientific
, Singapore
, 1994
).25.
J.
Dunkel
and I. M.
Zaid
, Phys. Rev. E
80
, 021903
(2009
).26.
M.
Reichert
, “Hydrodynamic interactions in colloidal and biological systems
,” Ph.D. thesis, University of Constance
, Germany
, 2006
.27.
R.
Joubad
, G. A.
Pavliotis
, and G.
Stoltz
, J. Stat. Phys.
158
, 1
(2015
).28.
G. N.
Milstein
and M. V.
Tretyakov
, Stochastic Numerics for Mathematical Physics
(Springer
, Berlin
, 2004
).29.
G. N.
Milstein
and M. V.
Tretyakov
, IMA J. Numer. Anal.
23
, 593
(2003
).30.
G. N.
Milstein
and M. V.
Tretyakov
, SIAM J. Numer. Anal.
34
, 2142
(1997
).31.
A.
Dullweber
, B.
Leimkuhler
, and R.
McLachlan
, J. Chem. Phys.
107
, 5840
(1997
).32.
E.
Hairer
, C.
Lubich
, and G.
Wanner
, Geometric Numerical Integration
, Volume 31 of Springer Series in Computational Mathematics (Springer-Verlag
, Berlin
, 2006
).33.
J. H.
Mentink
, M. V.
Tretyakov
, A.
Fasolino
, M. I.
Katsnelson
, and T.
Rasing
, J. Phys.: Condens. Matter
22
, 176001
(2010
).34.
B.
Leimkuhler
and C.
Matthews
, Appl. Math. Res. Express
2013
, 34
.35.
B.
Leimkuhler
and C.
Matthews
, Proc. R. Soc. A
472
, 20160138
(2016
).36.
G. N.
Milstein
and M. V.
Tretyakov
, SIAM J. Sci. Comput.
18
, 1067
(1997
).37.
M.
Reichert
and H.
Stark
, Phys. Rev. E
69
, 031407
(2004
).38.
S.
Martin
, M.
Reichert
, H.
Stark
, and T.
Gisler
, Phys. Rev. Lett.
97
, 248301
(2006
).39.
K.
Hinsen
, Comput. Phys. Commun.
88
, 327
(1995
), the Fortran source code is available at http://dirac.cnrs-orleans.fr/HYDROLIB.40.
B.
Cichocki
, B. U.
Felderhof
, K.
Hinsen
, E.
Wajnryb
, and J.
Bławzdziewicz
, J. Chem. Phys.
100
, 3780
(1994
).41.
J.-C.
Meiners
and S. R.
Quake
, Phys. Rev. Lett.
82
, 2211
(1999
).42.
P.
Bartlett
, S. I.
Henderson
, and S. J.
Mitchell
, Philos. Trans. R. Soc., A
359
, 883
(2001
).43.
J. P.
Eckmann
, L.
Thomas
, and P.
Wittwer
, J. Phys. A: Math. Gen.
14
, 3153
(1981
).44.
S. C.
Glotzer
and M. J.
Solomon
, Nat. Mater.
6
, 557
(2007
).45.
A. B.
Pawar
and I.
Kretzschmar
, Macromol. Rapid Commun.
31
, 150
(2010
).46.
Q.
Chen
, S. C.
Bae
, and S.
Granick
, Nature
469
, 381
(2011
).47.
A. C.
Newton
, J.
Groenewold
, W. K.
Kegel
, and P. G.
Bolhuis
, Proc. Nat. Acad. Sci. U. S. A.
112
, 15308
(2015
).48.
F.
Takagi
, N.
Koga
, and S.
Takada
, Proc. Nat. Acad. Sci. U. S. A.
100
, 11367
(2003
).49.
T. E.
Ouldridge
, R. L.
Hoare
, A. A.
Louis
, J. P. K.
Doye
, J.
Bath
, and A. J.
Turberfield
, ACS Nano
7
, 2479
(2013
).50.
R. R. F.
Machinek
, T. E.
Ouldridge
, N. E. C.
Haley
, J.
Bath
, and A. J.
Turberfield
, Nat. Commun.
5
, 5324
(2014
).51.
D.
Hinckley
, J. P.
Lequieu
, and J. J.
de Pablo
, J. Chem. Phys.
141
, 035102
(2014
).52.
A. J. C.
Ladd
, Phys. Rev. Lett.
70
, 1339
(1993
).53.
R. D.
Groot
and P. B.
Warren
, J. Chem. Phys.
107
, 4423
(1997
).54.
A.
Malavanets
and R.
Kapral
, J. Chem. Phys.
110
, 8605
(1999
).55.
E.
Lauga
, Phys. Rev. Lett.
106
, 178101
(2011
).56.
A.
Vijaykumar
, P. G.
Bolhuis
, and P. R.
ten Wolde
, Faraday Discuss.
195
, 421
(2016
).© 2017 Author(s).
2017
Author(s)
You do not currently have access to this content.
