skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Efficient reactive Brownian dynamics

 [1];  [2];  [3]
  1. Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA
  2. Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA, Department of Mathematics, UC Berkeley, Berkeley, California 94720, USA
  3. Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
Publication Date:
Sponsoring Org.:
OSTI Identifier:
Grant/Contract Number:
AC02-05CH11231; SC0008271
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Journal of Chemical Physics
Additional Journal Information:
Journal Volume: 148; Journal Issue: 3; Related Information: CHORUS Timestamp: 2018-02-14 18:32:46; Journal ID: ISSN 0021-9606
American Institute of Physics
Country of Publication:
United States

Citation Formats

Donev, Aleksandar, Yang, Chiao-Yu, and Kim, Changho. Efficient reactive Brownian dynamics. United States: N. p., 2018. Web. doi:10.1063/1.5009464.
Donev, Aleksandar, Yang, Chiao-Yu, & Kim, Changho. Efficient reactive Brownian dynamics. United States. doi:10.1063/1.5009464.
Donev, Aleksandar, Yang, Chiao-Yu, and Kim, Changho. 2018. "Efficient reactive Brownian dynamics". United States. doi:10.1063/1.5009464.
title = {Efficient reactive Brownian dynamics},
author = {Donev, Aleksandar and Yang, Chiao-Yu and Kim, Changho},
abstractNote = {},
doi = {10.1063/1.5009464},
journal = {Journal of Chemical Physics},
number = 3,
volume = 148,
place = {United States},
year = 2018,
month = 1

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on January 18, 2019
Publisher's Accepted Manuscript

Save / Share:
  • The [FeFe] hydrogenase from the green alga Chlamydomonas reinhardtii can catalyze the reduction of protons to hydrogen gas using electrons supplied from photosystem I and transferred via ferredoxin. To better understand the association of the hydrogenase and the ferredoxin, we have simulated the process over multiple timescales. A Brownian dynamics simulation method gave an initial thorough sampling of the rigid-body translational and rotational phase spaces, and the resulting trajectories were used to compute the occupancy and free-energy landscapes. Several important hydrogenase-ferredoxin encounter complexes were identified from this analysis, which were then individually simulated using atomistic molecular dynamics to provide moremore » details of the hydrogenase and ferredoxin interaction. The ferredoxin appeared to form reasonable complexes with the hydrogenase in multiple orientations, some of which were good candidates for inclusion in a transition state ensemble of configurations for electron transfer.« less
  • Highlights: {yields} Molecular Dynamics codes implemented on GPUs have achieved two-order of magnitude computational accelerations. {yields} Brownian Dynamics and Dissipative Particle Dynamics simulations require a large number of random numbers per time step. {yields} We introduce a method for generating small batches of pseudorandom numbers distributed over many threads of calculations. {yields} With this method, Dissipative Particle Dynamics is implemented on a GPU device without requiring thread-to-thread communication. - Abstract: Brownian Dynamics (BD), also known as Langevin Dynamics, and Dissipative Particle Dynamics (DPD) are implicit solvent methods commonly used in models of soft matter and biomolecular systems. The interaction ofmore » the numerous solvent particles with larger particles is coarse-grained as a Langevin thermostat is applied to individual particles or to particle pairs. The Langevin thermostat requires a pseudo-random number generator (PRNG) to generate the stochastic force applied to each particle or pair of neighboring particles during each time step in the integration of Newton's equations of motion. In a Single-Instruction-Multiple-Thread (SIMT) GPU parallel computing environment, small batches of random numbers must be generated over thousands of threads and millions of kernel calls. In this communication we introduce a one-PRNG-per-kernel-call-per-thread scheme, in which a micro-stream of pseudorandom numbers is generated in each thread and kernel call. These high quality, statistically robust micro-streams require no global memory for state storage, are more computationally efficient than other PRNG schemes in memory-bound kernels, and uniquely enable the DPD simulation method without requiring communication between threads.« less
  • The autocorrelation function of the velocity of an infinitely dilute Na/sup +/ ion in aqueous solution, and the autocorrelation function of the force exerted on a stationary Na/sup +/ under the same conditions are evaluated by molecular dynamics calculations. The results are used to test the accuracy of Brownian motion assumptions which are basic to hydrodynamic models of ion dynamics in solution. The self-diffusion coefficient of the Na/sup +/ ion predicted by Brownian motion theory is (0.65 +- 0.1) x 10/sup -5/cm/sup 2//s. This value is about 60% greater than the one obtained for the proper dynamics of the finitemore » mass ion, (0.4 +- 0.1) x 10/sup -5/cm/sup 2//s. The numerically correct velocity autocorrelation function is nonexponential, and the autocorrelation of the force on the stationary ion does not decay faster than the ion velocity autocorrelation function. Motivated by previous hydrodynamic modeling of friction kernels, we examine the approximation in which the memory function for the velocity autocorrelation function is identified with the autocorrelation function of the force on the stationary ion. The overall agreement between this approximation for the velocity autocorrelation function and the numerically correct answer is quite good.« less
  • A novel method is presented for introducing fluctuations in one-body dynamics. It consists of employing a Brownian force in the kinetic equations. For nuclear matter within the spinodal zone, the magnitude of the Brownian force can be determined by demanding correspondence with the growth of the most unstable mode, as given by Boltzmann-Langevin simulations. The method is illustrated and tested for idealized two-dimensional matter and promises to provide a practical means for addressing catastrophic nuclear processes.