English

Monte Carlo Algorithm for Simulating Reversible Aggregation of Multisite Particles

Quantitative Methods 2011-09-27 v2 Computational Physics

Abstract

We present an efficient and exact Monte Carlo algorithm to simulate reversible aggregation of particles with dedicated binding sites. This method introduces a novel data structure of dynamic bond tree to record clusters and sequences of bond formations. The algorithm achieves a constant time cost for processing cluster association and a cost between O(logM)\mathcal{O}(\log M) and O(M)\mathcal{O}(M) for processing bond dissociation in clusters with MM bonds. The algorithm is statistically exact and can reproduce results obtained by the standard method. We applied the method to simulate a trivalent ligand and a bivalent receptor clustering system and obtained an average scaling of O(M0.45)\mathcal{O}(M^{0.45}) for processing bond dissociation in acyclic aggregation, compared to a linear scaling with the cluster size in standard methods. The algorithm also demands substantially less memory than the conventional method.

Cite

@article{arxiv.1010.0339,
  title  = {Monte Carlo Algorithm for Simulating Reversible Aggregation of Multisite Particles},
  author = {Qiang Chang and Jin Yang},
  journal= {arXiv preprint arXiv:1010.0339},
  year   = {2011}
}

Comments

8 pages, 3 figures

R2 v1 2026-06-21T16:22:51.087Z