Quick search Find article
Quick search
Find article

A fast ergodic algorithm for generating ensembles of equilateral random polygons

R Varela1, K Hinson2, J Arsuaga3,4 and Y Diao2,5

Show affiliations


Knotted structures are commonly found in circular DNA and along the backbone of certain proteins. In order to properly estimate properties of these three-dimensional structures it is often necessary to generate large ensembles of simulated closed chains (i.e. polygons) of equal edge lengths (such polygons are called equilateral random polygons). However finding efficient algorithms that properly sample the space of equilateral random polygons is a difficult problem. Currently there are no proven algorithms that generate equilateral random polygons with its theoretical distribution. In this paper we propose a method that generates equilateral random polygons in a 'step-wise uniform' way. We prove that this method is ergodic in the sense that any given equilateral random polygon can be generated by this method and we show that the time needed to generate an equilateral random polygon of length n is linear in terms of n. These two properties make this algorithm a big improvement over the existing generating methods. Detailed numerical comparisons of our algorithm with other widely used algorithms are provided.


PACS

87.14.G- Nucleic acids

36.20.Fz Constitution (chains and sequences)

87.15.B- Structure of biomolecules

05.20.-y Classical statistical mechanics

MSC

92D20 Protein sequences, DNA sequences

92C40 Biochemistry, molecular biology

37A30 Ergodic theorems, spectral theory, Markov operators (For operator ergodic theory, see mainly 47A35)

Subjects

Soft matter, liquids and polymers

Atomic and molecular physics

Biological physics

Statistical physics and nonlinear systems

Dates

Issue 9 (6 March 2009)

Received 10 November 2008, in final form 7 January 2009

Published 4 February 2009



View by subject




Export








Please login to access our web services, or create an account if you don't yet have one.

You must have cookies enabled in your web browser to be able to login.

Username
Password

Forgotten your password? Get a new one here.