BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

A. C. Maggs

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

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechphysics.comp-phphysics.chem-phphysics.data-anquant-phBiological PhysicsComputational Complexitycond-mat.dis-nncond-mat.mes-hallmath-phmath.MPq-bio

Catalog footprint

What is connected

24works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

24 published item(s)

preprint2021arXiv

Hard-disk dipoles and non-reversible Markov chains

We benchmark event-chain Monte Carlo (ECMC) algorithms for tethered hard-disk dipoles in two dimensions in view of application of ECMC to water models in molecular simulation. We characterize the rotation dynamics of dipoles through the integrated autocorrelation times of the polarization. The non-reversible straight, reflective, forward, and Newtonian ECMC algorithms are all event-driven, and they differ only in their update rules at event times. They realize considerable speedups with respect to the local reversible Metropolis algorithm. We also find significant speed differences among the ECMC variants. Newtonian ECMC appears particularly well-suited for overcoming the dynamical arrest that has plagued straight ECMC for three-dimensional dipolar models with Coulomb interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Large-scale dynamics of event-chain Monte Carlo

Event-chain Monte Carlo (ECMC) accelerates the sampling of hard-sphere systems, and has been generalized to the potentials used in classical molecular simulation. Rather than imposing detailed balance on transition probabilities, the method enforces a weaker global-balance condition in order to guarantee convergence to equilibrium. In this paper we generalize the factor-field variant of ECMC to higher space dimensions. In the two-dimensional fluid phase, factor-field ECMC saturates the lower bound $z=0$ for the dynamical scaling exponent for local dynamics, whereas molecular dynamics is characterized by $z=1$ and local Metropolis Monte Carlo by $z=2$. In the presence of hexatic order, factor fields are not found to speed up the convergence. We indicate applications to the physics of glasses, and note that generalizations of factor fields could couple to orientational order.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Sparse hard-disk packings and local Markov chains

We propose locally stable sparse hard-disk packings, as introduced by Böröczky, as a model for the analysis and benchmarking of Markov-chain Monte Carlo (MCMC) algorithms. We first generate such packings in a square box with periodic boundary conditions and analyze their properties. We then study how local MCMC algorithms, namely the Metropolis algorithm and several versions of event-chain Monte Carlo (ECMC), escape from configurations that are obtained by slightly reducing all disk radii by a relaxation parameter. A scaling analysis is confirmed by simulation results. We obtain two classes of ECMC, one in which the escape time varies algebraically with the relaxation parameter (as for the local Metropolis algorithm) and another in which the escape time scales as the logarithm of the relaxation parameter. We discuss the connectivity of the hard-disk sample space, the ergodicity of local MCMC algorithms, as well as the meaning of packings in the context of the NPT ensemble. Our work is accompanied by open-source, arbitrary-precision software for Böröczky packings (in Python) and for straight, reflective, forward, and Newtonian ECMC (in Go).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multithreaded event-chain Monte Carlo with local times

We present a multithreaded event-chain Monte Carlo algorithm (ECMC) for hard spheres. Threads synchronize at infrequent breakpoints and otherwise scan for local horizon violations. Using a mapping onto absorbing Markov chains, we rigorously prove the correctness of a sequential-consistency implementation for small test suites. On x86 and ARM processors, a C++ (OpenMP) implementation that uses compare-and-swap primitives for data access achieves considerable speed-up with respect to single-threaded code. The generalized birthday problem suggests that for the number of threads scaling as the square root of the number of spheres, the horizon-violation probability remains small for a fixed simulation time. We provide C++ and Python open-source code that reproduces all our results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Phase Separation by Entanglement of Active Polymerlike Worms

We investigate the aggregation and phase separation of thin, living T. tubifex worms that behave as active polymers. Randomly dispersed active worms spontaneously aggregate to form compact, highly entangled blobs, a process similar to polymer phase separation, and for which we observe power-law growth kinetics. We find that the phase separation of active polymerlike worms does not occur through Ostwald ripening, but through active motion and coalescence of the phase domains. Interestingly, the growth mechanism differs from conventional growth by droplet coalescence: the diffusion constant characterizing the random motion of a worm blob is independent of its size, a phenomenon that can be explained from the fact that the active random motion arises from the worms at the surface of the blob. This leads to a fundamentally different phase-separation mechanism that may be unique to active polymers.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nonequilibrium Tuning of the Thermal Casimir Effect

In net-neutral systems correlations between charge fluctuations generate strong attractive thermal Casimir forces and engineering these forces to optimize nanodevice performance is an important challenge. We show how the normal and lateral thermal Casimir forces between two plates containing Brownian charges can be modulated by decorrelating the system through the application of an electric field, which generates a nonequilibrium steady state with a constant current in one or both plates, reducing the ensuing fluctuation-generated normal force while at the same time generating a lateral drag force. This hypothesis is confirmed by detailed numerical simulations as well as an analytical approach based on stochastic density functional theory.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The range and nature of effective interactions in hard-sphere solids

Colloidal systems observed in video microscopy are often analysed using the displacements correlation matrix of particle positions. In non-thermal systems, the inverse of this matrix can be interpreted as a pair-interaction potential between particles. If the system is thermally agitated, however, only an effective interaction is accessible from the correlation matrix. We show how this effective interaction differs from the non-thermal case by comparing with high-statistics numerical data from hard-sphere crystals.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cavity averages for hard spheres in the presence of polydispersity and incomplete data

We develop a cavity-based method which allows to extract thermodynamic properties from position information in hard-sphere/disk systems. So far, there are 'available-volume' and 'free-volume' methods. We add a third one, which we call 'available-volume-after-takeout', and which is shown to be mathematically equivalent to the others. In applications, where data sets are finite, all three methods show limitations, and they do this in different parameter ranges. We illustrate the principal equivalence and the limitations on data from molecular dynamics -- In particular, we test robustness against missing data. We have in mind experimental limitations where there is a small polydispersity, say 4% in the particle radii, but individual radii cannot be determined. We observe that, depending on the used method, the errors in such a situation are easily 100% for the pressure and 10kT for the chemical potentials. Our work is meant as guideline to the experimentalist for choosing the right one of the three methods, in order to keep the outcome of experimental data analysis meaningful.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Convexity and stiffness in energy functions for electrostatic simulations

We study the properties of convex functionals which have been proposed for the simulation of charged molecular systems within the Poisson-Boltzmann approximation. We consider the extent to which the functionals reproduce the true fluctuations of electrolytes and thus the one-loop correction to mean field theory -- including the Deby-Hückel correction to the free energy of ionic solutions. We also compare the functionals for use in numerical optimization of a mean field model of a charged polymer and show that different functionals have very different stiffnesses leading to substantial differences in accuracy and speed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

General theory of asymmetric steric interactions in electrostatic double layers

We study the Poisson-Boltzmann equation in the context of dense charged fluids where steric effects become important. We generalise the lattice gas theory by introducing a Flory-Huggins entropy for ions of differing volumes and then compare the effective free energy density to other approximations, valid for more realistic equations of state, such as the Carnahan-Starling approximation and find strong differences in the shapes of the free energy functions. We solve the Carnahan-Starling model in the high density limit, and demonstrate a slow, power-law convergence at high potentials. We elucidate how equivalent convex free energy functions can be constructed that describe steric effects in a manner which is more convenient for numerical minimisation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Electrostatic interactions in the presence of surface charge regulation: exact results

We study the problem of charge regulation and its effects on electrostatic interactions between dissociable charge groups immersed in a univalent electrolyte, within a family of one dimensional exactly solvable models. We consider the case of both charge regulated plates, but also the interaction of pairs of finite size dielectric "particles". Using the transfer matrix formalism we are able to determine the disjoining pressure as well as the correlations between the charge and the dipole moments of the objects as a function of their separation and electrolyte concentration

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Adding an energy-like conservation law to the leapfrog integrator

The leapfrog integrator is widely used because of its excellent stability in molecular dynamics simulation. This is recognized as being due to the existence of a discrete variational structure of the equations. We introduce a modified leapfrog method which includes an additional energy-like conservation law by embedding a molecular dynamics simulation within a larger dynamical system.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Phonons in two-dimensional soft colloidal crystals

The vibrational modes of pristine and polycrystalline monolayer colloidal crystals composed of thermosensitive microgel particles are measured using video microscopy and covariance matrix analysis. At low frequencies, the Debye relation for two dimensional harmonic crystals is observed in both crystal types; at higher frequencies, evidence for van Hove singularities in the phonon density of states is significantly smeared out by experimental noise and measurement statistics. The effects of these errors are analyzed using numerical simulations. We introduce methods to correct for these limitations, which can be applied to disordered systems as well as crystalline ones, and we show that application of the error correction procedure to the experimental data leads to more pronounced van Hove singularities in the pristine crystal. Finally, quasi-localized low-frequency modes in polycrystalline two-dimensional colloidal crystals are identified and demonstrated to correlate with structural defects such as dislocations, suggesting that quasi-localized low-frequency phonon modes may be used to identify local regions vulnerable to rearrangements in crystalline as well as amorphous solids.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Solving Fluctuation-Enhanced Poisson-Boltzmann Equations

Electrostatic correlations and fluctuations in ionic systems can be described within an extended Poisson-Boltzmann theory using a Gaussian variational form. The resulting equations are challenging to solve because they require the solution of a non-linear partial differential equation for the pair correlation function. This has limited existing studies to simple approximations or to one-dimensional geometries. In this paper we show that the numerical solution of the equations is greatly simplified by the use of selective inversion of a finite difference operator which occurs in the theory. This selective inversion preserves the sparse structure of the problem and leads to substantial savings in computer effort. In one and two dimensions further simplifications are made by using a mixture of selective inversion and Fourier techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A minimizing principle for the Poisson-Boltzmann equation

The Poisson-Boltzmann equation is often presented via a variational formulation based on the electrostatic potential. However, the functional has the defect of being non-convex. It can not be used as a local minimization principle while coupled to other dynamic degrees of freedom. We formulate a convex dual functional which is numerically equivalent at its minimum and which is more suited to local optimization.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Elastic fluctuations as observed in a confocal slice

Recent confocal experiments on colloidal solids motivate a fuller study of the projection of three-dimensional fluctuations onto a two-dimensional confocal slice. We show that the effective theory of a projected crystal displays several exceptional features, such as non-standard exponents in the dispersion relations. We provide analytic expressions for the effective two-dimensional elastic properties which allow one to work back from sliced experimental observations to three-dimensional elastic constants.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Legendre transforms for electrostatic energies

We review the use of Legendre transforms in the formulation of electrostatic energies in condensed matter. We show how to render standard functionals expressed in terms of the electrostatic potential, Phi, convex - at the cost of expressing them in terms of the vector field D. This leads to great simplification in the formulation of numerical minimisation of electrostatic energies coupled to other physical degrees of freedom. We also demonstrate the equivalence of recent functionals for dielectrics derived using field theory methods to classical formulations in terms of the electric polarisation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Anisotropic elasticity in confocal studies of colloidal crystals

We consider the theory of fluctuations of a colloidal solid observed in a confocal slice. For a cubic crystal we study the evolution of the projected elastic properties as a function of the anisotropy of the crystal using numerical methods based on the fast Fourier transform. In certain situations of high symmetry we find exact analytic results for the projected fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Truncated correlations in video microscopy of colloidal solids

Studies by video microscopy on fluctuating colloids measure the real-space cross-correlations in particle motion. This set of correlations is then treated as a matrix, in order to study the spectrum and mode structure. We show that in general the modes are modified by the truncation of the full real-space correlations. We perform a theoretical analysis of the truncation, find the boundary conditions imposed by the truncation, and propose practical windowing strategies to eliminate artefacts. We study the problem from various perspectives, to compile a survey for experimentalists.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Energy gaps in quantum first-order mean-field-like transitions: The problems that quantum annealing cannot solve

We study first-order quantum phase transitions in models where the mean-field traitment is exact, and the exponentially fast closure of the energy gap with the system size at the transition. We consider exactly solvable ferromagnetic models, and show that they reduce to the Grover problem in a particular limit. We compute the coefficient in the exponential closure of the gap using an instantonic approach, and discuss the (dire) consequences for quantum annealing.

0 citations0 reviewsNo code linkedOpen access