Source author record

Gerard Barkema

Gerard Barkema appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.stat-mech cond-mat.soft

Catalog footprint

What is connected

3works

2topics

3close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2014arXiv

An introduction to Monte Carlo methods

Monte Carlo simulations are methods for simulating statistical systems. The aim is to generate a representative ensemble of configurations to access thermodynamical quantities without the need to solve the system analytically or to perform an exact enumeration. The main principles of Monte Carlo simulations are ergodicity and detailed balance. The Ising model is a lattice spin system with nearest neighbor interactions that is appropriate to illustrate different examples of Monte Carlo simulations. It displays a second order phase transition between a disordered (high temperature) and ordered (low temperature) phases, leading to different strategies of simulations. The Metropolis algorithm and the Glauber dynamics are efficient at high temperature. Close to the critical temperature, where the spins display long range correlations, cluster algorithms are more efficient. We introduce the rejection free (or continuous time) algorithm and describe in details an interesting alternative representation of the Ising model using graphs instead of spins with the Worm algorithm. We conclude with an important discussion of the dynamical effects such as thermalization and correlation time.

preprint2013arXiv

Unwinding relaxation dynamics of polymers

The relaxation dynamics of a polymer wound around a fixed obstacle constitutes a fundamental instance of polymer with twist and torque and it is of relevance also for DNA denaturation dynamics. We investigate it by simulations and Langevin equation analysis. The latter predicts a relaxation time scaling as a power of the polymer length times a logarithmic correction related to the equilibrium fluctuations of the winding angle. The numerical data support this result and show that at short times the winding angle decreases as a power-law. This is also in agreement with the Langevin equation provided a winding-dependent friction is used, suggesting that such reduced description of the system captures the basic features of the problem.

preprint2011arXiv

Equilibrium winding angle of a polymer around a bar

The winding angle probability distribution of a planar self-avoiding walk has been known exactly since a long time: it has a gaussian shape with a variance growing as $<θ^2>\sim \ln L$. For the three-dimensional case of a walk winding around a bar, the same scaling is suggested, based on a first-order epsilon-expansion. We tested this three-dimensional case by means of Monte Carlo simulations up to length $L\approx25\,000$ and using exact enumeration data for sizes $L\le20$. We find that the variance of the winding angle scales as $<θ^2>\sim (\ln L)^{2α}$, with $α=0.75(1)$. The ratio $γ= <θ^4>/<θ^2>^2=3.74(5)$ is incompatible with the gaussian value $γ=3$, but consistent with the observation that the tail of the probability distribution function $p(θ)$ is found to decrease slower than a gaussian function. These findings are at odds with the existing first-order $ε$-expansion results.

Institution

Affiliation not imported yet

This author record came from a source that does not expose affiliation metadata. Once the author claims the profile or we enrich the record from another provider, this section will link to the concrete institution.

Topic footprint