Researcher profile

Gerard T. Barkema

Gerard T. Barkema contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

7 published item(s)

preprint2022arXiv

A comparison of cluster algorithms for the bond-diluted Ising model

Monte Carlo cluster algorithms are popular for their efficiency in studying the Ising model near its critical temperature. We might expect that this efficiency extends to the bond-diluted Ising model. We show, however, that this is not always the case by comparing how the correlation times $τ_w$ and $τ_{\rm sw}$ of the Wolff and Swendsen-Wang cluster algorithms scale as a function of the system size $L$ when applied to the two-dimensional bond-diluted Ising model. We demonstrate that the Wolff algorithm suffers from a much longer correlation time than in the pure Ising model, caused by isolated (groups of) spins which are infrequently visited by the algorithm. With a simple argument we prove that these cause the correlation time $τ_w$ to be bounded from below by $L^{z_w}$ with a dynamical exponent $z_w=γ/ ν\approx 1.75$ for a bond concentration $p < 1$. Furthermore, we numerically show that this lower bound is actually taken for several values of $p$ in the range $0.5 < p < 1$. Moreover, we show that the Swendsen-Wang algorithm does not suffer from the same problem. Consequently, it has a much shorter correlation time, shorter than in the pure Ising model even. Numerically at $p = 0.6$, we find that its dynamical exponent is $z_{\rm sw} = 0.09(4)$.

preprint2022arXiv

Structural dynamics of polycrystalline graphene

The exceptional properties of the two-dimensional material graphene make it attractive for multiple functional applications, whose large-area samples are typically polycrystalline. Here, we study the mechanical properties of graphene in computer simulations and connect these to the experimentally relevant mechanical properties. In particular, we study the fluctuations in the lateral dimensions of the periodic simulation cell. We show that over short time scales, both the area A and the aspect ratio B of the rectangular periodic box show diffusive behavior under zero external field during dynamical evolution, with diffusion coefficients DA and DB that are related to each other. At longer times, fluctuations in A are bounded, while those in B are not. This makes the direct determination of DB much more accurate, from which DA can then be derived indirectly. We then show that the dynamic behavior of polycrystalline graphene under external forces can also be derived from DA and DB via the Nernst-Einstein relation. Additionally, we study how the diffusion coefficients depend on structural properties of the polycrystalline graphene, in particular, the density of defects.

preprint2020arXiv

Super slowing down in the bond-diluted Ising model

In models in statistical physics, the dynamics often slows down tremendously near the critical point. Usually, the correlation time $τ$ at the critical point increases with system size $L$ in power-law fashion: $τ\sim L^z$, which defines the critical dynamical exponent $z$. We show that this also holds for the 2D bond-diluted Ising model in the regime $p>p_c$, where $p$ is the parameter denoting the bond concentration, but with a dynamical critical exponent $z(p)$ which shows a strong $p$-dependence. Moreover, we show numerically that $z(p)$, as obtained from the autocorrelation of the total magnetisation, diverges when the percolation threshold $p_c=1/2$ is approached: $z(p)-z(1) \sim (p-p_c)^{-2}$. We refer to this observed extremely fast increase of the correlation time with size as {\it super slowing down}. Independent measurement data from the mean-square deviation of the total magnetisation, which exhibits anomalous diffusion at the critical point, supports this result.

preprint2013arXiv

Dynamical Eigenmodes of Star and Tadpole Polymers

The dynamics of phantom bead-spring chains with the topology of a symmetric star with $f$ arms and tadpoles ($f=3$, a special case) is studied, in the overdamped limit. In the simplified case where the hydrodynamic radius of the central monomer is $f$ times as heavy as the other beads, we determine their dynamical eigenmodes exactly, along the lines of the Rouse modes for linear bead-spring chains. These eigenmodes allow full analytical calculations of virtually any dynamical quantity. As examples we determine the radius of gyration, the mean square displacement of a tagged monomer, and, for star polymers, the autocorrelation function of the vector that spans from the center of the star to a bead on one of the arms.

preprint2010arXiv

Elastic Lattice Polymers

We study a model of &#34;elastic&#34; lattice polymer in which a fixed number of monomers $m$ is hosted by a self-avoiding walk with fluctuating length $l$. We show that the stored length density $ρ_m = 1 - <l>/m$ scales asymptotically for large $m$ as $ρ_m=ρ_\infty(1-θ/m + ...)$, where $θ$ is the polymer entropic exponent, so that $θ$ can be determined from the analysis of $ρ_m$. We perform simulations for elastic lattice polymer loops with various sizes and knots, in which we measure $ρ_m$. The resulting estimates support the hypothesis that the exponent $θ$ is determined only by the number of prime knots and not by their type. However, if knots are present, we observe strong corrections to scaling, which help to understand how an entropic competition between knots is affected by the finite length of the chain.

preprint2010arXiv

Simulations of Two-Dimensional Unbiased Polymer Translocation Using the Bond Fluctuation Model

We use the Bond Fluctuation Model (BFM) to study the pore-blockade times of a translocating polymer of length $N$ in two dimensions, in the absence of external forces on the polymer (i.e., unbiased translocation) and hydrodynamic interactions (i.e., the polymer is a Rouse polymer), through a narrow pore. Earlier studies using the BFM concluded that the pore-blockade time scales with polymer length as $τ_d \sim N^β$, with $β=1+2ν$, whereas some recent studies with different polymer models produce results consistent with $β=2+ν$, originally predicted by us. Here $ν$ is the Flory exponent of the polymer; $ν=0.75$ in 2D. In this paper we show that for the BFM if the simulations are extended to longer polymers, the purported scaling $τ_d \sim N^{1+2ν}$ ceases to hold. We characterize the finite-size effects, and study the mobility of individual monomers in the BFM. In particular, we find that in the BFM, in the vicinity of the pore the individual monomeric mobilities are heavily suppressed in the direction perpendicular to the membrane. After a modification of the BFM which counters this suppression (but possibly introduces other artifacts in the dynamics), the apparent exponent $β$ increases significantly. Our conclusion is that BFM simulations do not rule out our theoretical prediction for unbiased translocation, namely $β=2+ν$.