Dynamics Brownian

The conceptual forerunner to mesoscale dynamics is Brownian dynamics. Brownian simulations used equations of motion modified by a random force 

The interested reader is referred to the rich polymeric materials literature for more details on the conformation and dynamics of polymers. 

Robert Brown s observation that small pollen grains are in perpetual irregular motion on the surface of water laid dormant for nearly eight decades, before Einstein explained it as the result of particle collisions with solvent molecules (Fig. 11.6). Of course, even Einstein was not entirely certain that what he was modeling was what Brown had observed. In 1905, Einstein wrote It is possible that the motions to be discussed here are identical with the so-called Brownian molecular motion however, the data available to me on the latter are so imprecise that I could not form a judgment on the question. It was actually Jean Perrin who accurately measured Brownian motion and connected Einstein s theory to Brown s observations in 1909. 

Einstein s theory of Brownian motion was ground-breaking because it substantiated the atomic hypothesis, or the Atoms Doctrine, as Perrin calls it in his wonderful book. First enunciated by the ancient Greeks, the idea that all matter consists of minute, indivisible entities called atoms became a sound scientific hypothesis with the work of Dalton, Boltzmann, and Mendeleyev. But it was not until Einstein worked on the 

Einstein determined, using osmotic pressure arguments, that the mean-squared displacement of a particle from an original position is proportional to time, 

It was three years later, in 1908, that Paul Langevin devised an alternative theory that gives the same correct result as Einstein s. Because Langevin s theory is applicable to a large class of random processes, it is the one predominantly discussed to explain Brownian motion. Here we describe its important elements. Langevin devised the following equation that describes the Brownian motion of a particle of mass M  

A valence bond method which uses different orbitals for different valence bond configurations. 

The theorem that states that there is no nonvanishing configuration interaction matrix elements between the ground state determinantal wavefunction and those determinants resulting from the excitation of one electron into an empty orbital of the initial SCF calculation. 

BSA = bovine serum albumin El = electrostatic interactions HI = hydrodynamic interactions PB = Poisson-Boltz-mann TIM = triose phosphate isomerase WEBDS = weighted ensemble Brownian dynamics simulations. 

In 1827, Robert Brown, a British botanist, observed pollen grains moving erratically in water. It was not until 1905 that Albert Einstein and Marian von Smoluchowski, independently, explained this phenomenon. What Brown had observed was the effect of the water molecules randomly colliding with the large pollen particle. These random collisions lead to an erratic displacement, i.e., diffusion, of the pollen particle in water. This chaotic motion is conunonly referred to as Brownian motion. Therefore Brownian motion is really a diffusion process. Using the concept of a one-dimensional random walk, Einstein found that the average displacement of a particle in 

The dynamics of the diffusion process can be simulated using Newton s second law of motion. The form of Newton s equation of motion used in simulating Brownian motion is 

---

Madura J D, Davis M E, Gilson, M K, Wade R C, Luty B A and McCammon J A 1994 Biological applications of electrostatic calculations and Brownian dynamics simulations Rev. Comput. Chem. 5 229-67... 

Northrup S H and Erickson H P 1992 Kinetics of protein-protein association explained by Brownian dynamics computer simulation Proc. Natl Acad. Sci. USA 89 3338-42... 

The first requirement is the definition of a low-dimensional space of reaction coordinates that still captures the essential dynamics of the processes we consider. Motions in the perpendicular null space should have irrelevant detail and equilibrate fast, preferably on a time scale that is separated from the time scale of the essential motions. Motions in the two spaces are separated much like is done in the Born-Oppenheimer approximation. The average influence of the fast motions on the essential degrees of freedom must be taken into account this concerns (i) correlations with positions expressed in a potential of mean force, (ii) correlations with velocities expressed in frictional terms, and iit) an uncorrelated remainder that can be modeled by stochastic terms. Of course, this scheme is the general idea behind the well-known Langevin and Brownian dynamics. 

Madura et al. 1995] Madura, J.D., Briggs, J.M., Wade, R.C., Davis, M.E., Luty, B.A., Ilin, A., Antosiewicz, J., Gilson, M.K., Bagheri, B., Scott, L.R., McCammon, J.A. Electrostatics and Diffusion of Molecules in Solution Simulations with the University of Houston Brownian Dynamics Program. Comp. Phys. Comm. 91 (1995) 57-95... 

J. A. Electrostatics and disffusion of molecules in solution simulations with the university of houston brownian dynamics program. Comp. Phys. Commun. 91 (1996) 57-95. 

R. C. Wade, M. E. Davis, B. A. Luty, J. D. Madura, and J. A. McCammon. Gating of the active site of triose phosphate isomerase Brownian dynamics simulations of flexible peptide loops in the enzyme. Biophys. J., 64 9-15, 1993. 

We further discuss how quantities typically measured in the experiment (such as a rate constant) can be computed with the new formalism. The computations are based on stochastic path integral formulation [6]. Two different sources for stochasticity are considered. The first (A) is randomness that is part of the mathematical modeling and is built into the differential equations of motion (e.g. the Langevin equation, or Brownian dynamics). The second (B) is the uncertainty in the approximate numerical solution of the exact equations of motion. 

We further note that the Langevin equation (which will not be discussed in detail here) is an intermediate between the Newton s equations and the Brownian dynamics. It includes in addition to an inertial part also a friction and a random force term ... 

Gunsteren W F and H J C Berendsen 1982. Algorithms for Brownian Dynamics. Molecular Physics 45 637-547. 

Gunsteren W F, H J C Berendsen and J A C Rullmann 1981. Stochastic Dynamics for Molecules with Constraints. Brownian Dynamics of n-Alkanes. Molecular Physics 44 69-95. 

Rossky P J, J D Doll and H L Friedman 1978. Brownian Dynamics as Smart Monte Carlo Simulatio Journal of Chemical Physics 69 4628-4633. 

Dissipative particle dynamics (DPD) is a technique for simulating the motion of mesoscale beads. The technique is superficially similar to a Brownian dynamics simulation in that it incorporates equations of motion, a dissipative (random) force, and a viscous drag between moving beads. However, the simulation uses a modified velocity Verlet algorithm to ensure that total momentum and force symmetries are conserved. This results in a simulation that obeys the Navier-Stokes equations and can thus predict flow. In order to set up these equations, there must be parameters to describe the interaction between beads, dissipative force, and drag. 

ME Davis, JD Madura, BA Luty, JA McCammon. Electrostatics and diffusion of molecules m solution Simulations with the University of Houston Brownian dynamics program. Comput Phys Commun 62 187-197, 1991. 

Analysis of neutron data in terms of models that include lipid center-of-mass diffusion in a cylinder has led to estimates of the amplitudes of the lateral and out-of-plane motion and their corresponding diffusion constants. It is important to keep in mind that these diffusion constants are not derived from a Brownian dynamics model and are therefore not comparable to diffusion constants computed from simulations via the Einstein relation. Our comparison in the previous section of the Lorentzian line widths from simulation and neutron data has provided a direct, model-independent assessment of the integrity of the time scales of the dynamic processes predicted by the simulation. We estimate the amplimdes within the cylindrical diffusion model, i.e., the length (twice the out-of-plane amplitude) L and the radius (in-plane amplitude) R of the cylinder, respectively, as follows ... 

But a computer simulation is more than a few clever data structures. We need algorithms to manipulate our system. In some way, we have to invent ways to let the big computer in our hands do things with the model that is useful for our needs. There are a number of ways for such a time evolution of the system the most prominent is the Monte Carlo procedure that follows an appropriate random path through configuration space in order to investigate equilibrium properties. Then there is molecular dynamics, which follows classical mechanical trajectories. There is a variety of dissipative dynamical methods, such as Brownian dynamics. All these techniques operate on the fundamental degrees of freedom of what we define to be our model. This is the common feature of computer simulations as opposed to other numerical approaches. 

Colloidal particles experience kicks from the surrounding atoms or molecules of the solvent. This leads to Brownian dynamics in colloidal suspensions (Fig. 14). The study of dynamics is challenging as, of course, first the equilibrium of the system has to be understood. One often knows the short-time dynamics that govern the system and is interested in long-time properties. 

Short-time Brownian motion was simulated and compared with experiments [108]. The structural evolution and dynamics [109] and the translational and bond-orientational order [110] were simulated with Brownian dynamics (BD) for dense binary colloidal mixtures. The short-time dynamics was investigated through the velocity autocorrelation function [111] and an algebraic decay of velocity fluctuation in a confined liquid was found [112]. Dissipative particle dynamics [113] is an attempt to bridge the gap between atomistic and mesoscopic simulation. Colloidal adsorption was simulated with BD [114]. The hydrodynamic forces, usually friction forces, are found to be able to enhance the self-diffusion of colloidal particles [115]. A novel MC approach to the dynamics of fluids was proposed in Ref. 116. Spinodal decomposition [117] in binary fluids was simulated. BD simulations for hard spherocylinders in the isotropic [118] and in the nematic phase [119] were done. A two-site Yukawa system [120] was studied with... 

FIG. 14 Brownian dynamics. The arrows indicate that the particle trajectories show diffusive behavior. 

FIG. 1 Results of a Brownian dynamic simulation for a two-dimensional coulombic system with specific interactions [40]. 

Subsequent work by Johansson and Lofroth [183] compared this result with those obtained from Brownian dynamics simulation of hard-sphere diffusion in polymer networks of wormlike chains. They concluded that their theory gave excellent agreement for small particles. For larger particles, the theory predicted a faster diffusion than was observed. They have also compared the diffusion coefficients from Eq. (73) to the experimental values [182] for diffusion of poly(ethylene glycol) in k-carrageenan gels and solutions. It was found that their theory can successfully predict the diffusion of solutes in both flexible and stiff polymer systems. Equation (73) is an example of the so-called stretched exponential function discussed further later. 

Johansson, L Lofroth, J-E, Diffusion and Interaction in Gels and Solutions. 4 Hard Sphere Brownian Dynamics Simulations, Journal of Chemical Physics 98, 7471, 1993. 

---