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

Peter Virnau

Peter Virnau 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-mechcond-mat.softcond-mat.mtrl-sciBiomoleculesBiological Physicsphysics.comp-phcond-mat.mes-hallGraphicsmath-phmath.MP

Catalog footprint

What is connected

26works
10topics
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

26 published item(s)

preprint2022arXiv

AlphaFold predicts the most complex protein knot and composite protein knots

The computer artificial intelligence system AlphaFold has recently predicted previously unknown three-dimensional structures of thousands of proteins. Focusing on the subset with high-confidence scores, we algorithmically analyze these predictions for cases where the protein backbone exhibits rare topological complexity, i.e. knotting. Amongst others, we discovered a $7_1$-knot, the most topologically complex knot ever found in a protein, as well several 6-crossing composite knots comprised of two methyltransferase or carbonic anhydrase domains, each containing a simple trefoil knot. These deeply embedded composite knots occur evidently by gene duplication and interconnection of knotted dimers. Finally, we report two new five-crossing knots including the first $5_1$-knot. Our list of analyzed structures forms the basis for future experimental studies to confirm these novel knotted topologies and to explore their complex folding mechanisms.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Constructing coarse-grained skyrmion potentials from experimental data with Iterative Boltzmann Inversion

In an effort to understand skyrmion behavior on a coarse-grained level, skyrmions are often described as 2D quasi particles evolving according to the Thiele equation. Interaction potentials are the key missing parameters for predictive modeling of experiments. We apply the Iterative Boltzmann Inversion technique commonly used in soft matter simulations to construct potentials for skyrmion-skyrmion and skyrmion-magnetic material boundary interactions from a single experimental measurement without any prior assumptions of the potential form. We find that the two interactions are purely repulsive and can be described by an exponential function for experimentally relevant skyrmions. This captures the physics on experimental time and length scales that are of interest for most skyrmion applications and typically inaccessible to atomistic or micromagnetic simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Skyrmion pinning energetics in thin film systems

A key issue for skyrmion dynamics and devices are pinning effects present in real systems. While posing a challenge for the realization of conventional skyrmionics devices, exploiting pinning effects can enable non-conventional computing approaches if the details of the pinning in real samples are quantified and understood. We demonstrate that using thermal skyrmion dynamics, we can characterize the pinning of a sample and we ascertain the spatially resolved energy landscape. To understand the mechanism of the pinning, we probe the strong skyrmion size and shape dependence of the pinning. Magnetic microscopy imaging demonstrates that in contrast to findings in previous investigations, for large skyrmions the pinning originates at the skyrmion boundary and not at its core. The boundary pinning is strongly influenced by the very complex pinning energy landscape that goes beyond the conventional rigid quasi-particle description. This gives rise to complex skyrmion shape distortions and allows for dynamic switching of pinning sites and flexible tuning of the pinning.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Commensurability between element symmetry and the number of skyrmions governing skyrmion diffusion in confined geometries

Magnetic skyrmions are topological magnetic structures, which exhibit quasi-particle properties and can show enhanced stability against perturbation from thermal noise. Recently, thermal Brownian diffusion of these quasi-particles has been found in continuous films and applications in unconventional computing have received significant attention, which however require structured elements. Thus, as the next necessary step, we here study skyrmion diffusion in confined geometries and find it to be qualitatively different: The diffusion is governed by the interplay between the total number of skyrmions and the structure geometry. In particular, we ascertain the effect of circular and triangular geometrical confinement and find that for triangular geometries the behavior is drastically different for the cases when the number of skyrmions in the element is either commensurate or incommensurate with a symmetric filling of the element. This influence of commensurability is corroborated by simulations of a quasi-particle model.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Critical behaviour in active lattice models of motility-induced phase separation

Lattice models allow for a computationally efficient investigation of motility-induced phase separation (MIPS) compared to off-lattice systems. Simulations are less demanding and thus bigger systems can be accessed with higher accuracy and better statistics. In equilibrium, lattice and off-lattice models with comparable interactions belong to the same universality class. Whether concepts of universality also hold for active particles is still a controversial and open question. Here, we examine two recently proposed active lattice systems that undergo MIPS and investigate numerically their critical behaviour. In particular, we examine the claim that these systems and MIPS in general belong to the Ising universality class. We also take a more detailed look on the influence and role of rotational diffusion and active velocity in these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Semiflexible Polymers Under Good Solvent Conditions Interacting With Repulsive Walls

Solutions of semiflexible polymers confined by repulsive planar walls are studied by density functional theory and Molecular Dynamics simulations, to clarify the competition between the chain alignment favored by the wall and the depletion caused by the monomer-wall repulsion. A coarse-grained bead-spring model with bond bending potential is studied, varying both the contour length and the persistence length of the polymers, as well as the monomer concentration in the solution (good solvent conditions are assumed throughout, and solvent molecules are not included explicitly). The profiles of monomer density and pressure tensor components near the wall are studied, and the surface tension of the solution is obtained. While the surface tension slightly decreases with chain length for flexible polymers, it clearly increases with chain length for stiff polymers. Thus, at fixed density and fixed chain length the surface tension also increases with increasing persistence length. Chain ends always are enriched near the wall, but this effect is much larger for stiff polymers than for flexible ones. Also the profiles of the mean square gyration radius components near the wall and the nematic order parameter are studied to clarify the conditions where wall-induced nematic order occurs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Crystal nuclei in melts: A Monte Carlo simulation of a model for attractive colloids

As a model for a suspension of hard-sphere like colloidal particles where small nonadsorbing dissolved polymers create a depletion attraction, we introduce an effective colloid-colloid potential closely related to the Asakura-Oosawa model but that does not have any discontinuities. In simulations, this model straightforwardly allows the calculation of the pressure from the Virial formula, and the phase transition in the bulk from the liquid to crystalline solid can be accurately located from a study where a stable coexistence of a crystalline slab with a surrounding liquid phase occurs. For this model, crystalline nuclei surrounded by fluid are studied both by identifying the crystal-fluid interface on the particle level (using suitable bond orientational order parameters to distinguish the phases) and by "thermodynamic" means. I.e., the latter method amounts to compute the enhancement of chemical potential and pressure relative to their coexistence values. We show that the chemical potential can be obtained from simulating thick films, where one wall with a rather long range repulsion is present, since near this wall the Widom particle insertion method works, exploiting the fact that the chemical potential in the system is homogeneous. Finally, the surface excess free energy of the nucleus is obtained, for a wide range of nuclei volumes. From this method, it is established that classical nucleation theory works, showing that for the present model the anisotropy of the interface excess free energy of crystals and their resulting nonspherical shape has only a very small effect on the barrier.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Finite size effects on liquid-solid phase coexistence and the estimation of crystal nucleation barriers

A fluid in equilibrium in a finite volume $V$ with particle number $N$ at a density $ρ= N/V$ exceeding the onset density $ρ_f $ of freezing may exhibit phase coexistence between a crystalline nucleus and surrounding fluid. Using a method suitable for the estimation of the chemical potential of dense fluids we obtain the excess free energy due to the surface of the crystalline nucleus. There is neither a need to precisely locate the interface nor to compute the (anisotropic) interfacial tension. As a test case, a soft version of the Asakura-Oosawa model for colloid polymer-mixtures is treated. While our analysis is appropriate for crystal nuclei of arbitrary shape, we find the nucleation barrier to be compatible with a spherical shape, and consistent with classical nucleation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Sequence determines degree of knottedness in a coarse-grained protein model

Knots are abundant in globular homopolymers but rare in globular proteins. To shed new light on this long-standing conundrum, we study the influence of sequence on the formation of knots in proteins under native conditions within the framework of the hydrophobic-polar (HP) lattice protein model. By employing large scale Wang-Landau simulations combined with suitable Monte Carlo trial moves we show that, even though knots are still abundant on average, sequence introduces large variability in the degree of self-entanglements. Moreover, we are able to design sequences which are either almost always or almost never knotted. Our findings serve as proof of concept that the introduction of just one additional degree of freedom per monomer (in our case sequence) facilitates evolution towards a protein universe in which knots are rare.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Anisotropic interfacial tension, contact angles, and line tensions: A graphics-processing-unit-based Monte Carlo study of the Ising model

As a generic example for crystals where the crystal-fluid interface tension depends on the orientation of the interface relative to the crystal lattice axes, the nearest neighbor Ising model on the simple cubic lattice is studied over a wide temperature range, both above and below the roughening transition temperature. Using a thin film geometry $L_x \times L_y \times L_z$ with periodic boundary conditions along the z-axis and two free $L_x \times L_y$ surfaces at which opposing surface fields $\pm H_{1}$ act, under conditions of partial wetting, a single planar interface inclined under a contact angle $θ< π/2$ relative to the yz-plane is stabilized. In the y-direction, a generalization of the antiperiodic boundary condition is used that maintains the translational invariance in y-direction despite the inhomogeneity of the magnetization distribution in this system. This geometry allows a simultaneous study of the angle-dependent interface tension, the contact angle, and the line tension (which depends on the contact angle, and on temperature). All these quantities are extracted from suitable thermodynamic integration procedures. In order to keep finite size effects as well as statistical errors small enough, rather large lattice sizes (of the order of 46 million sites) are found necessary, availability of very efficient code implementation of graphics processing units (GPUs) was crucial for the feasibility of this study.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Determination of the origin and magnitude of logarithmic finite-size effects on interfacial tension: Role of interfacial fluctuations and domain breathing

The ensemble-switch method for computing wall excess free energies of condensed matter is extended to estimate the interface free energies between coexisting phases very accurately. By this method, system geometries with linear dimensions $L$ parallel and $L_z$ perpendicular to the interface with various boundary conditions in the canonical or grandcanonical ensemble can be studied. Using two- and three-dimensional Ising models, the nature of the occurring logarithmic finite size corrections is studied. It is found crucial to include interfacial fluctuations due to "domain breathing".

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Logarithmic Finite-Size Effects on Interfacial Free Energies: Phenomenological Theory and Monte Carlo Studies

The computation of interfacial free energies between coexisting phases (e.g.~saturated vapor and liquid) by computer simulation methods is still a challenging problem due to the difficulty of an atomistic identification of an interface, and due to interfacial fluctuations on all length scales. The approach to estimate the interfacial tension from the free energy excess of a system with interfaces relative to corresponding single-phase systems does not suffer from the first problem but still suffers from the latter. Considering $d$-dimensional systems with interfacial area $L^{d-1}$ and linear dimension $L_z$ in the direction perpendicular to the interface, it is argued that the interfacial fluctuations cause logarithmic finite-size effects of order $\ln (L) / L^{d-1}$ and order $\ln (L_z)/L ^{d-1}$, in addition to regular corrections (with leading order $\text{const}/L^{d-1}$). A phenomenological theory predicts that the prefactors of the logarithmic terms are universal (but depend on the applied boundary conditions and the considered statistical ensemble). The physical origin of these corrections are the translational entropy of the interface as a whole, "domain breathing" (coupling of interfacial fluctuations to the bulk order parameter fluctuations of the coexisting domains), and capillary waves. Using a new variant of the ensemble switch method, interfacial tensions are found from Monte Carlo simulations of $d=2$ and $d=3$ Ising models and a Lennard Jones fluid. The simulation results are fully consistent with the theoretical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase Behavior of Active Swimmers in Depletants: Molecular Dynamics and Integral Equation Theory

We study the structure and phase behavior of a binary mixture where one of the components is self-propelling in nature. The inter-particle interactions in the system were taken from the Asakura-Oosawa model, for colloid-polymer mixtures, for which the phase diagram is known. In the current model version the colloid particles were made active using the Vicsek model for self-propelling particles. The resultant active system was studied by molecular dynamics methods and integral equation theory. Both methods produce results consistent with each other and demonstrate that the Vicsek model based activity facilitates phase separation, thus broadening the coexistence region.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scaling behavior of topologically constrained polymer rings in a melt

Large scale molecular dynamics simulations on graphic processing units (GPUs) are employed to study the scaling behavior of ring polymers with various topological constraints in melts. Typical sizes of rings containing $3_1$, $5_1$ knots and catenanes made up of two unknotted rings scale like $N^{1/3}$ in the limit of large ring sizes $N$. This is consistent with the crumpled globule model and similar findings for unknotted rings. For small ring lengths knots occupy a significant fraction of the ring. The scaling of typical ring sizes for small $N$ thus depends on the particular knot type and the exponent is generally larger than 0.4.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Monte Carlo Tests of Nucleation Concepts in the Lattice Gas Model

The conventional theory of homogeneous and heterogeneous nucleation in a supersaturated vapor is tested by Monte Carlo simulations of the lattice gas (Ising) model with nearest-neighbor attractive interactions on the simple cubic lattice. The theory considers the nucleation process as a slow (quasi-static) cluster (droplet) growth over a free energy barrier $ΔF^*$, constructed in terms of a balance of surface and bulk term of a "critical droplet" of radius $R^*$, implying that the rates of droplet growth and shrinking essentially balance each other for droplet radius $R=R^*$. For heterogeneous nucleation at surfaces, the barrier is reduced by a factor depending on the contact angle. Using the definition of "physical" clusters based on the Fortuin-Kasteleyn mapping, the time-dependence of the cluster size distribution is studied for "quenching experiments" in the kinetic Ising model, and the cluster size $\ell ^*$ where the cluster growth rate changes sign is estimated. These studies of nucleation kinetics are compared to studies where the relation between cluster size and supersaturation is estimated from equilibrium simulations of phase coexistence between droplet and vapor in the canonical ensemble. The chemical potential is estimated from a lattice version of the Widom particle insertion method. For large droplets it is shown that the "physical clusters" have a volume consistent with the estimates from the lever rule. "Geometrical clusters" (defined such that each site belonging to the cluster is occupied and has at least one occupied neighbor site) yield valid results only for temperatures less than 60% of the critical temperature, where the cluster shape is non-spherical. We show how the chemical potential can be used to numerically estimate $ΔF^*$ also for non-spherical cluster shapes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulation of fluid-solid coexistence in finite volumes: A method to study the properties of wall-attached crystalline nuclei

The Asakura-Oosawa model for colloid-polymer mixtures is studied by Monte Carlo simulations at densities inside the two-phase coexistence region of fluid and solid. Choosing a geometry where the system is confined between two flat walls, and a wall-colloid potential that leads to incomplete wetting of the crystal at the wall, conditions can be created where a single nanoscopic wall-attached crystalline cluster coexists with fluid in the remainder of the simulation box. Following related ideas that have been useful to study heterogeneous nucleation of liquid droplets at the vapor-liquid coexistence, we estimate the contact angles from observations of the crystalline clusters in thermal equilibrium. We find fair agreement with a prediction based on Young's equation, using estimates of interface and wall tension from the study of flat surfaces. It is shown that the pressure versus density curve of the finite system exhibits a loop, but the pressure maximum signifies the "droplet evaporation-condensation" transition and thus has nothing in common with a van der Waals-like loop. Preparing systems where the packing fraction is deep inside the two-phase coexistence region, the system spontaneously forms a "slab state", with two wall-attached crystalline domains separated by (flat) interfaces from liquid in full equilibrium with the crystal in between; analysis of such states allows a precise estimation of the bulk equilibrium properties at phase coexistence.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Anomalous Structure and Scaling of Ring Polymer Brushes

A comparative simulation study of polymer brushes formed by grafting at a planar surface either flexible linear polymers (chain length $N_L$) or (non-catenated) ring polymers (chain length $N_R=2 N_L$) is presented. Two distinct off-lattice models are studied, one by Monte Carlo methods, the other by Molecular Dynamics, using a fast implementation on graphics processing units (GPUs). It is shown that the monomer density profiles $ρ(z)$ in the $z$-direction perpendicular to the surface for rings and linear chains are practically identical, $ρ_R(2 N_L, z)=ρ_L(N_L, z)$. The same applies to the pressure, exerted on a piston at hight z, as well. While the gyration radii components of rings and chains in $z$-direction coincide, too, and increase linearly with $N_L$, the transverse components differ, even with respect to their scaling properties: $R_{gxy}^{(L)} \propto N_L^{1/2}$, $R_{gxy}^{(R)} \propto N_L^{0.4}$. These properties are interpreted in terms of the statistical properties known for ring polymers in dense melts.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hard sphere fluids confined between soft repulsive walls: A comparative study using Monte Carlo and density functional methods

Hard-sphere fluids confined between parallel plates a distance $D$ apart are studied for a wide range of packing fractions, including also the onset of crystallization, applying Monte Carlo simulation techniques and density functional theory. The walls repel the hard spheres (of diameter $σ$) with a Weeks-Chandler-Andersen (WCA) potential $V_{WCA}(z) = 4 ε[(σ_w/z)^{12}-(σ_w/z)^6 + 1/4]$, with range $σ_w = σ/2$. We vary the strength $ε$ over a wide range and the case of simple hard walls is also treated for comparison. By the variation of $ε$ one can change both the surface excess packing fraction and the wall-fluid $(γ_{wf})$ and wall-crystal $(γ_{wc})$ surface free energies. Several different methods to extract $γ_{wf}$ and $γ_{wc}$ from Monte Carlo (MC) simulations are implemented, and their accuracy and efficiency is comparatively discussed. The density functional theory (DFT) using Fundamental Measure functionals is found to be quantitatively accurate over a wide range of packing fractions; small deviations between DFT and MC near the fluid to crystal transition need to be studied further. Our results on density profiles near soft walls could be useful to interpret corresponding experiments with suitable colloidal dispersions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Methods to Compute Pressure and Wall Tension in Fluids containing Hard Particles

Colloidal systems are often modelled as fluids of hard particles (possibly with an additional soft attraction, e.g. caused by polymers also contained in the suspension). in simulations of such systems, the virial theorem cannot be straightforwardly applied to obtain the components of the pressure tensor. In systems confined by walls, it is hence also not straightforward to extract the excess energy due to the wall (the "wall tension") from the pressure tensor anisotropy. A comparative evaluation of several methods to circumvent this problem is presented, using as examples fluids of hard spheres and the Asakura-Oosawa model of colloid-polymer mixtures with a size ratio $q=0.15$ (for which the effect of the polymers can be integrated out to yield an effective attractive potential between the colloids). Factors limiting the accuracy of the various methods are carefully discussed, and controlling these factors very good mutual agreement between the various methods is found.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Monte Carlo Methods for Estimating Interfacial Free Energies and Line Tensions

Excess contributions to the free energy due to interfaces occur for many problems encountered in the statistical physics of condensed matter when coexistence between different phases is possible (e.g. wetting phenomena, nucleation, crystal growth, etc.). This article reviews two methods to estimate both interfacial free energies and line tensions by Monte Carlo simulations of simple models, (e.g. the Ising model, a symmetrical binary Lennard-Jones fluid exhibiting a miscibility gap, and a simple Lennard-Jones fluid). One method is based on thermodynamic integration. This method is useful to study flat and inclined interfaces for Ising lattices, allowing also the estimation of line tensions of three-phase contact lines, when the interfaces meet walls (where "surface fields" may act). A generalization to off-lattice systems is described as well. The second method is based on the sampling of the order parameter distribution of the system throughout the two-phase coexistence region of the model. Both the interface free energies of flat interfaces and of (spherical or cylindrical) droplets (or bubbles) can be estimated, including also systems with walls, where sphere-cap shaped wall-attached droplets occur. The curvature-dependence of the interfacial free energy is discussed, and estimates for the line tensions are compared to results from the thermodynamic integration method. Basic limitations of all these methods are critically discussed, and an outlook on other approaches is given.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Capillary Condensation in Cylindrical Pores: Monte Carlo Study of the Interplay of Surface and Finite Size Effects

When a fluid that undergoes a vapor to liquid transition in the bulk is confined to a long cylindrical pore, the phase transition is shifted (mostly due to surface effects at the walls of the pore) and rounded (due to finite size effects). The nature of the phase coexistence at the transition depends on the length of the pore: For very long pores the system is axially homogeneous at low temperatures. At the chemical potential where the transition takes place fluctuations occur between vapor-like and liquid-like states of the cylinder as a whole. At somewhat higher temperatures (but still far below bulk criticality) the system at phase coexistence is in an axially inhomogeneous multi-domain state, where long cylindrical liquid-like and vapor-like domains alternate. Using Monte Carlo simulations for the Ising/lattice gas model and the Asakura-Oosawa model of colloid-polymer mixtures the transition between these two different scenarios is characterized. It is shown that the density distribution changes gradually from a double-peak structure to a triple-peak shape, and the correlation length in axial direction (measuring the equilibrium domain length) becomes much smaller than the cylinder length. The (rounded) transition to the disordered phase of the fluid occurs when the axial correlation length has decreased to a value comparable to the cylinder diameter. It is also suggested that adsorption hysteresis vanishes when the transition from the simple domain state to the multi-domain state of the cylindrical pore occurs. We predict that the difference between the pore critical temperature and the hysteresis critical emperature should increase logarithmically with the length of the pore.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Monte Carlo Test of the Classical Theory for Heterogeneous Nucleation Barriers

Flat walls facilitate the condensation of a supersaturated vapor: Classical theory of heterogeneous nucleation predicts that the free energy barrier $ΔF_{\rm het}^*$ which needs to be overcome for the formation of sphere-cap shaped nucleation seeds is smaller than the barrier $ΔF^*_{\rm hom}$ for spherical droplets in the bulk by a factor $0<f(θ)<1$, which only depends on the contact angle $θ$. In this letter we compute both $ΔF^*_{\rm hom}$ and $ΔF^*_{\rm het}$ from Monte Carlo simulations and test the theory for the lattice gas model (for which $θ$ can be readily controlled). Even though the theory is only based on macroscopic arguments, it is shown to hold for experimentally relevant nanoscopic nucleation seeds ($20\leqΔF^*_{\rm hom}/k_BT\leq 200)$ if (independently estimated) line tension effects are considered.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multi-GPU Accelerated Multi-Spin Monte Carlo Simulations of the 2D Ising Model

A modern graphics processing unit (GPU) is able to perform massively parallel scientific computations at low cost. We extend our implementation of the checkerboard algorithm for the two dimensional Ising model [T. Preis et al., J. Comp. Phys. 228, 4468 (2009)] in order to overcome the memory limitations of a single GPU which enables us to simulate significantly larger systems. Using multi-spin coding techniques, we are able to accelerate simulations on a single GPU by factors up to 35 compared to an optimized single Central Processor Unit (CPU) core implementation which employs multi-spin coding. By combining the Compute Unified Device Architecture (CUDA) with the Message Parsing Interface (MPI) on the CPU level, a single Ising lattice can be updated by a cluster of GPUs in parallel. For large systems, the computation time scales nearly linearly with the number of GPUs used. As proof of concept we reproduce the critical temperature of the 2D Ising model using finite size scaling techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Rounding of Phase Transitions in Cylindrical Pores

Phase transitions of systems confined in long cylindrical pores (capillary condensation, wetting, crystallization, etc.) are intrinsically not sharply defined but rounded. The finite size of the cross section causes destruction of long range order along the pore axis by spontaneous nucleation of domain walls. This rounding is analyzed for two models (Ising/lattice gas and Asakura-Oosawa model for colloid-polymer mixtures) by Monte Carlo simulations and interpreted by a phenomenological theory. We show that characteristic differences between the behavior of pores of finite length and infinitely long pores occur. In pores of finite length a rounded transition occurs first, from phase coexistence between two states towards a multi-domain configuration. A second transition to the axially homogeneous phase follows near pore criticality.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Implementation and performance analysis of bridging Monte Carlo moves for off-lattice single chain polymers in globular states

Bridging algorithms are global Monte Carlo moves which allow for an efficient sampling of single polymer chains. In this manuscript we discuss the adaptation of three bridging algorithms from lattice to continuum models, and give details on the corrections to the acceptance rules which are required to fulfill detailed balance. For the first time we are able to compare the efficiency of the moves by analyzing the occurrence of knots in globular states. For a flexible homopolymer chain of length N=1000, independent configurations can be generated up to two orders of magnitude faster than with slithering snake moves.

0 citations0 reviewsNo code linkedOpen access