BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Wolfhard Janke appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We propose a simple method for uniform sampling of points on the surface of a hypersphere in arbitrarily many dimensions. By avoiding the evaluation of computationally expensive functions like logarithms, sines, cosines, or higher order roots the new method is faster than alternative techniques.
Population annealing (PA) is a population-based algorithm that is designed for equilibrium simulations of thermodynamic systems with a rough free energy landscape. It is known to be more efficient in doing so than standard Markov chain Monte Carlo alone. The algorithm has a number of parameters that can be fine-tuned to improve performance. While there is some theoretical and numerical work regarding most of these parameters, there appears to be a gap in the literature concerning the role of resampling in PA. Here, we present a numerical comparison of a number of resampling schemes for PA simulations of the 2D Ising model.
We study the nonequilibrium kinetics during the coil-globule transition of a flexible polymer chain with active beads after a quench from good to poor solvent condition using molecular dynamics simulation. Activity for each bead is introduced via the well-known Vicsek-like alignment rule due to which the velocity of a bead tries to align towards the average direction of its neighbors. We investigate the role of quenching temperature with varying activity during collapse of this polymer. We find that although for lower activities the kinetics remains qualitatively similar for different temperatures, for higher activity noticeable differences can be identified.
Simulations of systems with quenched disorder are extremely demanding, suffering from the combined effect of slow relaxation and the need of performing the disorder average. As a consequence, new algorithms, improved implementations, and alternative and even purpose-built hardware are often instrumental for conducting meaningful studies of such systems. The ensuing demands regarding hardware availability and code complexity are substantial and sometimes prohibitive. We demonstrate how with a moderate coding effort leaving the overall structure of the simulation code unaltered as compared to a CPU implementation, very significant speed-ups can be achieved from a parallel code on GPU by mainly exploiting the trivial parallelism of the disorder samples and the near-trivial parallelism of the parallel tempering replicas. A combination of this massively parallel implementation with a careful choice of the temperature protocol for parallel tempering as well as efficient cluster updates allows us to equilibrate comparatively large systems with moderate computational resources.
Met-enkephalin, one of the smallest opiate peptides and an important neurotransmitter, is a widely used benchmarking problem in the field of molecular simulation. Through its range of possible low-temperature conformations separated by free-energy barriers it was previously found to be hard to thermalize using straight canonical molecular dynamics simulations. Here, we demonstrate how one can use the recently proposed population annealing molecular dynamics scheme to overcome these difficulties. We show how the use of multi-histogram reweighting allows one to accurately estimate the density of states of the system and hence derive estimates such as the potential energy as quasi continuous functions of temperature. We further investigate the free-energy surface as a function of end-to-end distance and radius-of-gyration and observe two distinct basins of attraction.
Collapse of the polypeptide backbone is an integral part of protein folding. Using polyglycine as a probe, we explore the nonequilibrium pathways of protein collapse in water. We find that the collapse depends on the competition between hydration effects and intra-peptide interactions. Once intra-peptide van der Waal interactions dominate, the chain collapses along a nonequilibrium pathway characterized by formation of pearl-necklace-like local clusters as intermediates that eventually coagulate into a single globule. By describing this coarsening through the contact probability as a function of distance along the chain, we extract a time-dependent length scale that grows in linear fashion. The collapse dynamics is characterized by a dynamical critical exponent $z=0.5$ that is much smaller than the values of $z=1-2$ reported for non-biological polymers. This difference in the exponents is explained by the instantaneous formation of intra-chain hydrogen bonds and local ordering that may be correlated with the observed fast folding times in proteins.
Recent emerging interest in experiments of single-polymer dynamics urge computational physicists to revive their understandings, particularly in the nonequilibrium context. Here we briefly discuss the currently evolving approaches of investigating the evolution dynamics of homopolymer collapse in computer simulations. Primary focus of these approaches is to understand various dynamic scaling laws related to coarsening and aging during the collapse in space dimension $d=3$, using tools popular in nonequilibrium coarsening dynamics of particle or spin systems. In addition to providing an overview of those results, we also present new preliminary data for $d=2$.
From a coarse-grained perspective the motif of a self-activating species, activating a second species which acts as its own repressor, is widely found in biological systems, in particular in genetic systems with inherent oscillatory behavior. Here we consider a specific realization of this motif as a genetic circuit, in which genes are described as directly producing proteins, leaving out the intermediate step of mRNA production. We focus on the effect that inherent time scales on the underlying fine-grained scale can have on the bifurcation patterns on a coarser scale in time. Time scales are set by the binding and unbinding rates of the transcription factors to the promoter regions of the genes. Depending on the ratio of these rates to the decay times of the proteins, the appropriate averaging procedure for obtaining a coarse-grained description changes and leads to sets of deterministic equations, which differ in their bifurcation structure. In particular the desired intermediate range of regular limit cycles fades away when the binding rates of genes are of the same order or less than the decay time of at least one of the proteins. Our analysis illustrates that the common topology of the widely found motif alone does not necessarily imply universal features in the dynamics.
Driven diffusive systems such as the zero-range process (ZRP) and the pair-factorized steady states (PFSS) stochastic transport process are versatile tools that lend themselves to the study of transport phenomena on a generic level. While their mathematical structure is simple enough to allow significant analytical treatment, they offer a variety of interesting phenomena. With appropriate dynamics, the ZRP and PFSS models feature a condensation transition where for a supercritical density the translational symmetry breaks spontaneously and excess particles form a single-site or spatially extended condensate, respectively. In this paper we numerically study the typical time scales of the two stages of this condensation process: Nucleation and coarsening. Nucleation is the first stage of condensation where the bulk system relaxes to its stationary distribution and droplet nuclei form in the system. These droplets then gradually grow or evaporate in the coarsening regime to finally coalesce in a single condensate when the system finally relaxes to the stationary state. We use the ZRP condensation model to discuss the choice of the estimation method for the nucleation time scale and present scaling exponents for the ZRP and PFSS condensation models with respect to the choice of the typical droplet mass. We then proceed to present scaling exponents in the coarsening regime of the ZRP for partial-asymmetric dynamics and the PFSS model for symmetric and asymmetric dynamics.
We discuss the effects of open boundary conditions and boundary induced drift on condensation phenomena in the pair-factorized steady states transport process, a versatile model for stochastic transport with tunable nearest-neighbour interactions. Varying the specific type of the boundary implementation as well as the presence of a particle drift, we observe phase diagrams that are similar but richer compared to those of the simpler zero-range process with open boundary conditions. Tuning our model towards zero-range-process-like properties we are able to study boundary induced effects in the transition regime from zero-range interactions to short-range interactions. We discuss the emerging phase structure where spatially extended condensates can be observed at the boundaries as well as in the bulk system and compare it to the situation with periodic boundaries, where the dynamics leads to the formation of a single condensate in the bulk.
An anisotropic limit of the 3d plaquette Ising model, in which the plaquette couplings in one direction were set to zero, was solved for free boundary conditions by Suzuki (Phys. Rev. Lett. 28 (1972) 507), who later dubbed it the fuki-nuke, or "no-ceiling", model. Defining new spin variables as the product of nearest-neighbour spins transforms the Hamiltonian into that of a stack of (standard) 2d Ising models and reveals the planar nature of the magnetic order, which is also present in the fully isotropic 3d plaquette model. More recently, the solution of the fuki-nuke model was discussed for periodic boundary conditions, which require a different approach to defining the product spin transformation, by Castelnovo et al. (Phys. Rev. B 81 (2010) 184303). We clarify the exact relation between partition functions with free and periodic boundary conditions expressed in terms of original and product spin variables for the 2d plaquette and 3d fuki-nuke models, noting that the differences are already present in the 1d Ising model. In addition, we solve the 2d plaquette Ising model with helical boundary conditions. The various exactly solved examples illustrate how correlations can be induced in finite systems as a consequence of the choice of boundary conditions.
We propose a generic model of driven DNA under the influence of an oscillatory force of amplitude $F$ and frequency $ν$ and show the existence of a dynamical transition for a chain of finite length. We find that the area of the hysteresis loop, $A_{\rm loop}$, scales with the same exponents as observed in a recent study based on a much more detailed model. However, towards the true thermodynamic limit, the high-frequency scaling regime extends to lower frequencies for larger chain length $L$ and the system has only one scaling ($A_{\rm loop} \approx ν^{-1}F^2)$. Expansion of an analytical expression for $A_{\rm loop}$ obtained for the model system in the low-force regime revealed that there is a new scaling exponent associated with force ($A_{\rm loop} \approx ν^{-1}F^{2.5}$), which has been validated by high-precision numerical calculation. By a combination of analytical and numerical arguments, we also deduce that for large but finite $L$, the exponents are robust and independent of temperature and friction coefficient.
We present a simple non-equilibrium model of mass condensation with Lennard-Jones interactions between particles and the substrate. We show that when some number of particles is deposited onto the surface and the system is left to equilibrate, particles condense into an island if the density of particles becomes higher than some critical density. We illustrate this with numerically obtained phase diagrams for three-dimensional systems. We also solve a two-dimensional counterpart of this model analytically and show that not only the phase diagram but also the shape of the cross-sections of three-dimensional condensates qualitatively matches the two-dimensional predictions. Lastly, we show that when particles are being deposited with a constant rate, the system has two phases: a single condensate for low deposition rates, and multiple condensates for fast deposition. The behaviour of our model is thus similar to that of thin film growth processes, and in particular to Stranski-Krastanov growth.
We discuss the effects of particle exchange through open boundaries and the induced drive on the phase structure and condensation phenomena of a stochastic transport process with tunable short-range interactions featuring pair-factorized steady states (PFSS) in the closed system. In this model, the steady state of the particle hopping process can be tuned to yield properties from the zero-range process (ZRP) condensation model to those of models with spa- tially extended condensates. By varying the particle exchange rates as well as the presence of a global drift, we observe a phase transition from a free particle gas to a phase with condensates aggregated to the boundaries. While this transition is similar to previous results for the ZRP, we find that the mechanism is different as the presence of the boundary actually influences the interaction due to the non-zero interaction range.
We present a finite-size scaling analysis of the droplet condensation-evaporation transition of a lattice gas (in two and three dimensions) and a Lennard-Jones gas (in three dimensions) at fixed density. Parallel multicanonical simulations allow sampling of the required system sizes with precise equilibrium estimates. In the limit of large systems, we verify the theoretical leading-order scaling prediction for both the transition temperature and the finite-size rounding. In addition, we present an emerging intermediate scaling regime, consistent in all considered cases and with similar recent observations for polymer aggregation. While the intermediate regime locally may show a different effective scaling, we show that it is a gradual crossover to the large-system scaling behavior by including empirical higher-order corrections. This implies that care has to be taken when considering scaling ranges, possibly leading to completely wrong predictions for the thermodynamic limit. In this study, we consider a crossing of the phase boundary orthogonal to the usual fixed temperature studies. We show that this is an equivalent approach and, under certain conditions, may show smaller finite-size corrections.
We study the low-temperature properties of the classical three-dimensional compass or $t_{2g}$ orbital model on simple-cubic lattices by means of comprehensive large-scale Monte Carlo simulations. Our numerical results give evidence for a directionally ordered phase that is reached via a first-order transition at the temperature $T_0 = 0.098328(3) J / k_{\mathrm{B}}$. To obtain our results we employ local and cluster update algorithms, parallel tempering and multiple histogram reweighting as well as model-specific screw-periodic boundary conditions, which help counteract severe finite-size effects.
We study the aggregation transition of a finite theta-polymer system in dependence on the bending stiffness $κ$ with the help of parallel multicanonical simulations. In order to distinguish amorphous aggregates from polymer bundles we introduce an order parameter, measuring the correlation of the end-to-end vectors. With the help of this order parameter, we construct generic $T$-$κ$ phase diagrams for systems with $2$ and $8$ polymers and discuss the occurring phases from amorphous aggregates to bundle structures. For an intermediate stiffness range we find multiple aggregated phases which change with increasing number of polymers and discuss their nature with the help of microcanonical analyses. We show that the stiffness of semiflexible theta polymers is the distinguishing parameter for the emergent structural motifs.
Using a combination of the replica-exchange Monte Carlo algorithm and the multicanonical method, we investigate the influence of bending stiffness on the conformational phases of a bead-stick homopolymer model and present the pseudo-phase diagram for the complete range of semi-flexible polymers, from flexible to stiff. Although a simple model, we observe a rich variety of conformational phases, reminiscent of conformations observed for synthetic polymers or biopolymers. Changing the internal bending stiffness, the model exhibits different pseudo phases like bent, hairpin or toroidal. In particular, we find thermodynamically stable knots and transitions into these knotted phases with a clear phase coexistence, but almost no change in the mean total energy.
The purely plaquette 3d Ising Hamiltonian with the spins living at the vertices of a cubic lattice displays several interesting features. The symmetries of the model lead to a macroscopic degeneracy of the low-temperature phase and prevent the definition of a standard magnetic order parameter. Consideration of the strongly anisotropic limit of the model suggests that a layered, "fuki-nuke" order still exists and we confirm this with multicanonical simulations. The macroscopic degeneracy of the low-temperature phase also changes the finite-size scaling corrections at the first-order transition in the model and we see this must be taken into account when analysing our measurements.
We numerically survey predictions on the shapes and scaling laws of particle condensates that emerge as a result of spontaneous symmetry breaking in pair- factorized steady states of a stochastic transport process. The specific model consists of indistinguishable particles that stochastically hop between sites controlled by a tunable potential. We identify the different condensate shapes within their respective parameter regimes as well as determine precisely the condensate width scaling.
Using numerical methods we discuss the effects of open boundary conditions on condensation phenomena in the zero-range process (ZRP) and transport processes with pair-factorized steady states (PFSS), an extended model of the ZRP with nearest-neighbor interaction. For the zero-range process we compare to analytical results in the literature with respect to criticality and condensation. For the extended model we find a similar phase structure, but observe supercritical phases with droplet formation for strong boundary drives.
We study the thermodynamic properties of the three-dimensional Blume-Capel model on the simple cubic lattice by means of computer simulations. In particular, we implement a parallelized variant of the multicanonical approach and perform simulations by keeping a constant temperature and crossing the phase boundary along the crystal-field axis. We obtain numerical data for several temperatures in both the first- and second-order regime of the model. Finite-size scaling analyses provide us with transition points and the dimensional scaling behavior in the numerically demanding first-order regime, as well as a clear verification of the expected Ising universality in the respective second-order regime. Finally, we discuss the scaling behavior in the vicinity of the tricritical point.
In this paper we conduct a careful multicanonical simulation of the isotropic $3d$ plaquette ("gonihedric") Ising model and confirm that a planar, fuki-nuke type order characterises the low-temperature phase of the model. From consideration of the anisotropic limit of the model we define a class of order parameters which can distinguish the low- and high-temperature phases in both the anisotropic and isotropic cases. We also verify the recently voiced suspicion that the order parameter like behaviour of the standard magnetic susceptibility $χ_m$ seen in previous Metropolis simulations was an artefact of the algorithm failing to explore the phase space of the macroscopically degenerate low-temperature phase. $χ_m$ is therefore not a suitable order parameter for the model.
We construct the complete structural phase diagram of polymer adsorption at substrates with attractive stripe-like patterns in the parameter space spanned by the adsorption affinity of the stripes and temperature. Results were obtained by extensive generalized-ensemble Monte Carlo simulations of a generic model for the hybrid organic-inorganic system. By comparing with adhesion properties at homogeneous substrates, we find substantial differences in the formation of adsorbed polymer structures if translational invariance at the surface is broken by a regular pattern. Beside a more specific understanding of polymer adsorption processes, our results are potentially relevant for the design of macromolecular pattern recognition devices such as sensors.
We investigate the aggregation transition of theta polymers in spherical confinement with multicanonical simulations. This allows for a systematic study of the effect of density on the aggregation transition temperature for up to 24 monodisperse polymers. Our results for solutions in the dilute regime show that polymers can be considered isolated for all temperatures larger than the aggregation temperature, which is shown to be a function of the density. The resulting competition between single-polymer collapse and aggregation yields the lower temperature bound of the isolated chain approximation. We provide entropic and energetic arguments to describe the density dependence and finite-size effects of the aggregation transition for monodisperse solutions in finite systems. This allows us to estimate the aggregation transition temperature of dilute systems in a spherical cavity, using a few simulations of small, sufficiently dilute polymer systems.
We present the speedup from a novel parallel implementation of the multicanonical method on the example of a lattice gas in two and three dimensions. In this approach, all cores perform independent equilibrium runs with identical weights, collecting their sampled histograms after each iteration in order to estimate consecutive weights. The weights are then redistributed to all cores. These steps are repeated until the weights are converged. This procedure benefits from a minimum of communication while distributing the necessary amount of statistics efficiently. Using this method allows us to study a broad temperature range for a variety of large and complex systems. Here, a gas is modeled as particles on the lattice, which interact only with their nearest neighbors. For a fixed density this model is equivalent to the Ising model with fixed magnetization. We compare our results to an analytic prediction for equilibrium droplet formation, confirming that a single macroscopic droplet forms only above a critical density.
We study self-avoiding walks on three-dimensional critical percolation clusters using a new exact enumeration method. It overcomes the exponential increase in computation time by exploiting the clusters' fractal nature. We enumerate walks of over $10^4$ steps, far more than has ever been possible. The scaling exponent $ν$ for the end-to-end distance turns out to be smaller than previously thought and appears to be the same on the backbones as on full clusters. We find strong evidence against the widely assumed scaling law for the number of conformations and propose an alternative, which perfectly fits our data.
The three-dimensional purely plaquette gonihedric Ising model and its dual are investigated to resolve inconsistencies in the literature for the values of the inverse transition temperature of the very strong temperature-driven first-order phase transition that is apparent in the system. Multicanonical simulations of this model allow us to measure system configurations that are suppressed by more than 60 orders of magnitude compared to probable states. With the resulting high-precision data, we find excellent agreement with our recently proposed nonstandard finite-size scaling laws for models with a macroscopic degeneracy of the low-temperature phase by challenging the prefactors numerically. We find an overall consistent inverse transition temperature of 0.551334(8) from the simulations of the original model both with periodic and fixed boundary conditions, and the dual model with periodic boundary conditions. For the original model with periodic boundary conditions, we obtain the first reliable estimate of the interface tension, 0.12037(18), using the statistics of suppressed configurations.
We note that the standard inverse system volume scaling for finite-size corrections at a first-order phase transition (i.e., 1/L^3 for an L x L x L lattice in 3D) is transmuted to 1/L^2 scaling if there is an exponential low-temperature phase degeneracy. The gonihedric Ising model which has a four-spin interaction, plaquette Hamiltonian provides an exemplar of just such a system. We use multicanonical simulations of this model to generate high-precision data which provides strong confirmation of the non-standard finite-size scaling law. The dual to the gonihedric model, which is an anisotropically coupled Ashkin-Teller model, has a similar degeneracy and also displays the non-standard scaling.
The most conspicuous property of a semiflexible polymer is its persistence length, defined as the decay length of tangent correlations along its contour. Using an efficient stochastic growth algorithm to sample polymers embedded in a quenched two-dimensional hard-disk fluid, we find apparent wormlike chain statistics with a renormalized persistence length. We identify a universal form of the disorder renormalization that suggests itself as a quantitative measure of molecular crowding.
It is known that fixed boundary conditions modify the leading finite-size corrections for an L^3 lattice in 3d at a first-order phase transition from 1/L^3 to 1/L. We note that an exponential low-temperature phase degeneracy of the form 2^3L will lead to a different leading correction of order 1/L^2 . A 3d gonihedric Ising model with a four-spin interaction, plaquette Hamiltonian displays such a degeneracy and we confirm the modified scaling behaviour using high-precision multicanonical simulations. We remark that other models such as the Ising antiferromagnet on the FCC lattice, in which the number of "true" low-temperature phases is not macroscopically large but which possess an exponentially degenerate number of low lying states may display an effective version of the modified scaling law for the range of lattice sizes accessible in simulations.
We study the conformational properties of heteropolymers containing two types of monomers A and B, modeled as self-avoiding random walks on a regular lattice. Such a model can describe in particular the sequences of hydrophobic and hydrophilic residues in proteins (K.F. Lau and K.A. Dill, Macromolecules {\bf 22}, 3986 (1989)) and polyampholytes with oppositely charged groups (Y. Kantor and M. Kardar, Europhys. Lett.{\bf 28}, 169 (1994)). Treating the sequences of the two types of monomers as quenched random variables, we provide a systematic analysis of possible generalizations of this model. To this end we apply the pruned-enriched Rosenbluth chain-growth algorithm (PERM), which allows us to obtain the phase diagrams of extended and compact states coexistence as function of both the temperature and fraction of A and B monomers along the heteropolymer chain.
We implemented a parallel version of the multicanonical algorithm and applied it to a variety of systems with phase transitions of first and second order. The parallelization relies on independent equilibrium simulations that only communicate when the multicanonical weight function is updated. That way, the Markov chains efficiently sample the temporary distributions allowing for good estimations of consecutive weight functions. The systems investigated range from the well known Ising and Potts spin systems to bead-spring polymers. We estimate the speedup with increasing number of parallel processes. Overall, the parallelization is shown to scale quite well. In the case of multicanonical simulations of the $q$-state Potts model ($q\ge6$) and multimagnetic simulations of the Ising model, the optimal performance is limited due to emerging barriers.
The investigation of freezing transitions of single polymers is computationally demanding, since surface effects dominate the nucleation process. In recent studies we have systematically shown that the freezing properties of flexible, elastic polymers depend on the precise chain length. Performing multicanonical Monte Carlo simulations, we faced several computational challenges in connection with liquid-solid and solid-solid transitions. For this reason, we developed novel methods and update strategies to overcome the arising problems. We introduce novel Monte Carlo moves and two extensions to the multicanonical method.
We applied the simulated tempering and magnetizing (STM) method to the two-dimensional three-state Potts model in an external magnetic field in order to perform further investigations of the STM's applicability. The temperature as well as the external field are treated as dynamical variables updated during the STM simulations. After we obtained adequate information for several lattice sizes $L$ (up to $160\times 160$), we also performed a number of conventional canonical simulations of large lattices, especially in order to illustrate the crossover behavior of the Potts model in external field with increasing $L$. The temperature and external field for larger lattice size simulations were chosen by extrapolation of the detail information obtained by STM. We carefully analyzed the crossover scaling at the phase transitions with respect to the lattice size as well as the temperature and external field. The crossover behavior is clearly observed in the simulations in agreement with theoretical predictions.
We study the static properties of a semiflexible polymer exposed to a quenched random environment by means of computer simulations. The polymer is modeled as two-dimensional Heisenberg chain. For the random environment we consider hard disks arranged on a square lattice. We apply an off-lattice growth algorithm as well as the multicanonical Monte Carlo method to investigate the influence of both disorder occupation probability and polymer stiffness on the equilibrium properties of the polymer. We show that the additional length scale induced by the stiffness of the polymer extends the well-known phenomenology considerably. The polymer's response to the disorder is either contraction or extension depending on the ratio of polymer stiffness and void space extension. Additionally, the periodic structure of the lattice is reflected in the observables that characterize the polymer.
We study the adsorption of flexible polymer macromolecules on a percolation cluster, formed by a regular two-dimensional disordered lattice at critical concentration p_c of attractive sites. The percolation cluster is characterized by a fractal dimension d_s^{p_c}=91/49. The conformational properties of polymer chains grafted to such a fractal substrate are studied by means of the pruned-enriched Rosenbluth method (PERM). We find estimates for the surface crossover exponent governing the scaling of the adsorption energy in the vicinity of the transition point, ϕ_s^{p_c}=0.425\pm0.009, and for the adsorption transition temperature, T_A^{p_c}=2.64\pm0.02. As expected, the adsorption is diminished when the fractal dimension of the substrate is smaller than that of a plain Euclidean surface. The universal size and shape characteristics of a typical spatial conformation which attains a polymer chain in the adsorbed state are analyzed as well.
We perform equilibrium computer simulations of a two-dimensional pinned flexible polymer exposed to a quenched disorder potential consisting of hard disks. We are especially interested in the high-density regime of the disorder, where subtle structures such as cavities and channels play a central role. We apply an off-lattice growth algorithm proposed by Garel and Orland [J. Phys. A 23, L621 (1990)], where a distribution of polymers is constructed in parallel by growing each of them monomer by monomer. In addition we use a multicanonical Monte Carlo method in order to cross-check the results of the growth algorithm. We measure the end-to-end distribution and the tangent-tangent correlations. We also investigate the scaling behavior of the mean square end-to-end distance in dependence of the monomer number. While the influence of the potential in the low-density case is merely marginal, it dominates the configurational properties of the polymer for high densities.
We investigate the cooperative effects of a single finite chain of monomers near an attractive substrate by first constructing a conformational pseudo-phase diagram based on the thermal fluctuations of energetic and structural quantities. Then, the adsorption transition is analyzed in more detail. This is conveniently done by a microcanonical analysis of densities of states obtained by extensive multicanonical Monte Carlo simulations. For short chains and strong surface attraction, the microcanonical entropy turns out to be a convex function of energy in the transition regime. This is a characteristic physical effect and deserves a careful consideration in analyses of cooperative macrostate transitions in finite systems.
We compare the thermodynamic behavior of a finite single nongrafted polymer near an attractive substrate with that of a polymer grafted to that substrate. After we recently found first-order-like signatures in the microcanonical entropy at the adsorption transition in the nongrafted case, and given the fact that many studies on polymer adsorption in the past have been performed for grafted polymers, the question arises, to what extent and in what way does grafting change the nature of the adsorption transition? This question is tackled here using a coarse-grained off-lattice polymer model and covers not only the adsorption transition but also all other transitions a single polymer near an attractive substrate of varying strengths undergoes. Because of the impact of grafting especially on the translational but also on the conformational entropy of desorbed chains, the adsorption transition is affected the strongest. Our results are obtained by a combined canonical and microcanonical analysis of parallel tempering Monte Carlo data.
We discuss the hierarchy of subphase transitions in first-order-like nucleation processes for an exemplified aggregation transition of heteropolymers. We perform an analysis of the microcanonical entropy, i.e., the density of states is considered as the central statistical system quantity since it connects system-specific entropic and energetic information in a natural and unique way.
We discuss the advantages of parallelization by multithreading on graphics processing units (GPUs) for parallel tempering Monte Carlo computer simulations of an exemplified bead-spring model for homopolymers. Since the sampling of a large ensemble of conformations is a prerequisite for the precise estimation of statistical quantities such as typical indicators for conformational transitions like the peak structure of the specific heat, the advantage of a strong increase in performance of Monte Carlo simulations cannot be overestimated. Employing multithreading and utilizing the massive power of the large number of cores on GPUs, being available in modern but standard graphics cards, we find a rapid increase in efficiency when porting parts of the code from the central processing unit (CPU) to the GPU.
The design of hybrid peptide-solid interfaces for nanotechnological applications such as biomolecular nanoarrays requires a deep understanding of the basic mechanisms of peptide binding and assembly at solid substrates. Here we show by means of experimental and computational analyses that the adsorption properties of mutated synthetic peptides at semiconductors exhibit a clear sequence-dependent adhesion specificity. Our simulations of a novel hybrid peptide-substrate model reveal the correspondence between proline mutation and binding affinity to a clean silicon substrate. After synthesizing theoretically suggested amino-acid sequences with different binding behavior, we confirm the relevance of the selective mutations upon adhesion in our subsequent atomic force microscopy experiments.
In finite-size scaling analyses of Monte Carlo simulations of second-order phase transitions one often needs an extended temperature range around the critical point. By combining the parallel tempering algorithm with cluster updates and an adaptive routine to find the temperature window of interest, we introduce a flexible and powerful method for systematic investigations of critical phenomena. As a result, we gain one to two orders of magnitude in the performance for 2D and 3D Ising models in comparison with the recently proposed Wang-Landau recursion for cluster algorithms based on the multibondic algorithm, which is already a great improvement over the standard multicanonical variant.
In finite-size scaling analyses of Monte Carlo simulations of second-order phase transitions one often needs an extended temperature/energy range around the critical point. By combining the replica-exchange algorithm with cluster updates and an adaptive routine to find the range of interest, we introduce a new flexible and powerful method for systematic investigations of critical phenomena. As a result, we gain two further orders of magnitude in the performance for 2D and 3D Ising models in comparison with the recently proposed Wang-Landau recursion for cluster algorithms based on the multibondic algorithm, which is already a great improvement over the standard multicanonical variant.
By means of multicanonical computer simulations, we investigate thermodynamic properties of the aggregation of interacting semiflexible polymers. We analyze a mesoscopic bead-stick model, where nonbonded monomers interact via Lennard-Jones forces. Aggregation turns out to be a process, in which the constituents experience strong structural fluctuations, similar to peptides in coupled folding-binding cluster formation processes. In contrast to a recently studied related proteinlike hydrophobic-polar heteropolymer model, aggregation and crystallization are separate processes for a homopolymer with the same small bending rigidity. Rather stiff semiflexible polymers form a liquid-crystal-like phase, as expected. In analogy to the heteropolymer study, we find that the first-order-like aggregation transition of the complexes is accompanied by strong system-size dependent hierarchical surface effects. In consequence, the polymer aggregation is a phase-separation process with entropy reduction.
We analyze the structural behavior of a single polymer chain grafted to an attractive, flexible surface. Our model is composed of a coarse-grained bead-and-spring polymer and a tethered membrane. By means of extensive parallel tempering Monte Carlo simulations it is shown that the system exhibits a rich phase behavior ranging from highly ordered, compact to extended random coil structures and from desorbed to completely adsorbed or even partially embedded conformations. These findings are summarized in a pseudophase diagram indicating the predominant class of conformations as a function of the external parameters temperature and polymer-membrane interaction strength. By comparison with adsorption to a stiff membrane surface it is shown that the flexibility of the membrane gives rise to qualitatively new behavior such as stretching of adsorbed conformations.
A Monte Carlo simulation study of the critical and off-critical behavior of the Baxter-Wu model, which belongs to the universality class of the 4-state Potts model, was performed. We estimate the critical temperature window using known analytical results for the specific heat and magnetization. This helps us to extract reliable values of universal combinations of critical amplitudes with reasonable accuracy. Comparisons with approximate analytical predictions and other numerical results are discussed.
The loop gas approach to lattice field theory provides an alternative, geometrical description in terms of fluctuating loops. Statistical ensembles of random loops can be efficiently generated by Monte Carlo simulations using the worm update algorithm. In this paper, concepts from percolation theory and the theory of self-avoiding random walks are used to describe estimators of physical observables that utilize the nature of the worm algorithm. The fractal structure of the random loops as well as their scaling properties are studied. To support this approach, the O(1) loop model, or high-temperature series expansion of the Ising model, is simulated on a honeycomb lattice, with its known exact results providing valuable benchmarks.
We investigate solid-solid and solid-liquid transitions of elastic flexible off-lattice polymers with Lennard-Jones monomer-monomer interaction and anharmonic springs by means of sophisticated variants of multicanonical Monte Carlo methods. We find that the low-temperature behavior depends strongly and non-monotonically on the system size and exhibits broad similarities to unbound atomic clusters. Particular emphasis is dedicated to the classification of icosahedral and non-icosahedral low-energy polymer morphologies.
Besides the well-known effect of autocorrelations in time series of Monte Carlo simulation data resulting from the underlying Markov process, using the same data pool for computing various estimates entails additional cross correlations. This effect, if not properly taken into account, leads to systematically wrong error estimates for combined quantities. Using a straightforward recipe of data analysis employing the jackknife or similar resampling techniques, such problems can be avoided. In addition, a covariance analysis allows for the formulation of optimal estimators with often significantly reduced variance as compared to more conventional averages.
In this work we investigate structural properties of native states of a simple model for short flexible homopolymers, where the steric influence of monomeric side chains is effectively introduced by a thickness constraint. This geometric constraint is implemented through the concept of the global radius of curvature and affects the conformational topology of ground-state structures. A systematic analysis allows for a thickness-dependent classification of the dominant ground-state topologies. It turns out that helical structures, strands, rings, and coils are natural, intrinsic geometries of such tubelike objects.
Various two-dimensional O(N) models coupled to Euclidean quantum gravity, whose intrinsic dimension is four, are shown to belong to universality classes of nongravitating statistical models in a lower number of dimensions. It is speculated that the matching critical behaviors in the gravitating and dimensionally reduced models may be manifestations of the holographic principle.
We study the directional-ordering transition in the two-dimensional classical and quantum compass models on the square lattice by means of Monte Carlo simulations. An improved algorithm is presented which builds on the Wolff cluster algorithm in one-dimensional subspaces of the configuration space. This improvement allows us to study classical systems up to $L=512$. Based on the new algorithm we give evidence for the presence of strongly anomalous scaling for periodic boundary conditions which is much worse than anticipated before. We propose and study alternative boundary conditions for the compass model which do not make use of extended configuration spaces and show that they completely remove the problem with finite-size scaling. In the last part, we apply these boundary conditions to the quantum problem and present a considerably improved estimate for the critical temperature which should be of interest for future studies on the compass model. Our investigation identifies a strong one-dimensional magnetic ordering tendency with a large correlation length as the cause of the unusual scaling and moreover allows for a precise quantification of the anomalous length scale involved.
We study the influence of structural obstacles in a disordered environment on the size and shape characteristics of long flexible polymer macromolecules. We use the model of self-avoiding random walks on diluted regular lattices at the percolation threshold in space dimensions d=2, 3. Applying the Pruned-Enriched Rosenbluth Method (PERM), we numerically estimate rotationally invariant universal quantities such as the averaged asphericity A_d and prolateness S of polymer chain configurations. Our results quantitatively reveal the extent of anisotropy of macromolecules due to the presence of structural defects.
Investigating thermodynamic properties of liquid-solid transitions of flexible homopolymers with elastic bonds by means of multicanonical Monte Carlo simulations, we find crystalline conformations that resemble ground-state structures of Lennard-Jones clusters. This allows us to set up a structural classification scheme for finite-length flexible polymers and their freezing mechanism in analogy to atomic cluster formation. Crystals of polymers with "magic length" turn out to be perfectly icosahedral.
We study the peculiarities of stretching of globular polymer macromolecules in a disordered (crowded) environment, using the model of self-attracting self-avoiding walks on site-diluted percolative lattices in space dimensions d=3. Applying the pruned-enriched Rosenbluth chain-growth method (PERM), we construct the phase diagram of collapsed-extended state coexistence when varying temperature and stretching force. The change in shape characteristics of globular polymers under stretching is analyzed as well.
The fractal structure of high-temperature graphs of the three-dimensional Ising and XY models is investigated by simulating these graphs directly on a cubic lattice and analyzing them with the help of percolation observables. The Ising graphs are shown to percolate right at the Curie critical point. The diverging length scale relevant to the graphs in the vicinity of the percolation threshold is shown to be provided by the spin correlation length. The fractal dimension of the high-temperature graphs at criticality is estimated to be $D = 1.7349(65)$ for the Ising and $D = 1.7626(66)$ for the XY model.
We study conformational transitions of simple coarse-grained models for protein-like heteropolymers on the simple cubic lattice and off-lattice, respectively, by means of multicanonical sampling algorithms. The effective hydrophobic/polar models do not require the knowledge of the native topology for a given sequence of residues as input. Therefore these models are eligible to investigate general properties of the tertiary folding behaviour of such protein-like heteropolymers.
We have performed chain-growth simulations of minimalistic hybrid lattice models for polymers interacting with interfaces of attractive solid substrates in order to gain insights into the conformational transitions of the polymers in the adsorption process. Primarily focusing on the dependence of the conformational behavior on temperature and solubility we obtained pseudophase diagrams with a detailed structure of conformational subphases. In the study of hydrophobic-polar peptides in the vicinity of different types of substrates, we found a noticeable substrate specificity of the assembly of hydrophobic domains in the conformations dominating the adsorption subphases.
In these notes, the application of Feynman's sum-over-paths approach to thermal phase transitions is discussed. The paradigm of such a spacetime approach to critical phenomena is provided by the high-temperature expansion of spin models. This expansion, known as the hopping expansion in the context of lattice field theory, yields a geometric description of the phase transition in these models, with the thermal critical exponents being determined by the fractal structure of the high-temperature graphs. The graphs percolate at the thermal critical point and can be studied using purely geometrical observables known from percolation theory. Besides the phase transition in spin models and in the closely related $ϕ^4$ theory, other transitions discussed from this perspective include Bose-Einstein condensation, and the transitions in the Higgs model and the pure U(1) gauge theory.
We have performed exhaustive multicanonical Monte Carlo simulations of three 12-residue synthetic peptides in order to investigate the thermodynamic and structural properties as well as the characteristic helix-coil transitions. In these studies, we employ a realistic model where the interactions between all atoms are taken into account. Effects of solvation are also simulated by using an implicit solvent model.
These notes give examples of how suitably defined geometrical objects encode in their fractal structure thermal critical behavior. The emphasis is on the two-dimensional Potts model for which two types of spin clusters can be defined. Whereas the Fortuin-Kasteleyn clusters describe the standard critical behavior, the geometrical clusters describe the tricritical behavior that arises when including vacant sites in the pure Potts model. Other phase transitions that allow for a geometrical description discussed in these notes include the superfluid phase transition and Bose-Einstein condensation.
The three-dimensional lattice Higgs model with compact U(1) gauge symmetry and unit charge is investigated by means of Monte Carlo simulations. The full model with fluctuating Higgs amplitude is simulated, and both energy as well as topological observables are measured. The data show a Higgs and a confined phase separated by a well-defined phase boundary, which is argued to be caused by proliferating vortices. For fixed gauge coupling, the phase boundary consists of a line of first-order phase transitions at small Higgs self-coupling, ending at a critical point. The phase boundary then continues as a Kertesz line across which thermodynamic quantities are nonsingular. Symmetry arguments are given to support these findings.
The six-vertex F model on the square lattice constitutes the unique example of an exactly solved model exhibiting an infinite-order phase transition of the Kosterlitz-Thouless type. As one of the few non-trivial exactly solved models, it provides a welcome gauge for new numerical simulation methods and scaling techniques. In view of the notorious problems of clearly resolving the Kosterlitz-Thouless scenario in the two-dimensional XY model numerically, the F model in particular constitutes an instructive reference case for the simulational description of this type of phase transition. We present a loop-cluster update Monte Carlo study of the square-lattice F model, with a focus on the properties not exactly known such as the polarizability or the scaling dimensions in the critical phase. For the analysis of the simulation data, finite-size scaling is explicitly derived from the exact solution and plausible assumptions. Guided by the available exact results, the careful inclusion of correction terms in the scaling formulae allows for a reliable determination of the asymptotic behaviour.
We study the phase transition of the three-dimensional complex |psi|^4 theory by considering the geometrically defined vortex-loop network as well as the magnetic properties of the system in the vicinity of the critical point. Using high-precision Monte Carlo techniques we examine an alternative formulation of the geometrical excitations in relation to the global O(2)-symmetry breaking, and check if both of them exhibit the same critical behavior leading to the same critical exponents and therefore to a consistent description of the phase transition. Different percolation observables are taken into account and compared with each other. We find that different definitions of constructing the vortex-loop network lead to different results in the thermodynamic limit, and the percolation thresholds do not coincide with the thermodynamic phase transition point.
We study by extensive Monte Carlo simulations the effect of random bond dilution on the phase transition of the three-dimensional 4-state Potts model which is known to exhibit a strong first-order transition in the pure case. The phase diagram in the dilution-temperature plane is determined from the peaks of the susceptibility for sufficiently large system sizes. In the strongly disordered regime, numerical evidence for softening to a second-order transition induced by randomness is given. Here a large-scale finite-size scaling analysis, made difficult due to strong crossover effects presumably caused by the percolation fixed point, is performed.
The two-dimensional Ising model with Brascamp-Kunz boundary conditions has a partition function more amenable to analysis than its counterpart on a torus. This fact is exploited to exactly determine the full finite-size scaling behaviour of the Fisher zeroes of the model. Moreover, exact results are also determined for the scaling of the specific heat at criticality, for the specific-heat peak and for the pseudocritical points. All corrections to scaling are found to be analytic and the shift exponent $λ$ does not coincide with the inverse of the correlation length exponent $1/ν$.
We report on a new method to extract thermodynamic properties from the density of partition function zeroes on finite lattices. This allows direct determination of the order and strength of phase transitions numerically. Furthermore, it enables efficient distinguishing between first- and second-order transitions, elucidates crossover between them and illuminates the origins of finite-size scaling. The power of the method is illustrated in typical applications for both Fisher and Lee-Yang zeroes.
Using an elaborate set of simulational tools and statistically optimized methods of data analysis we investigate the scaling behavior of the correlation lengths of three-dimensional classical O($n$) spin models. Considering three-dimensional slabs $S^1\times S^1\times\mathbb{R}$, the results over a wide range of $n$ indicate the validity of special scaling relations involving universal amplitude ratios that are analogous to results of conformal field theory for two-dimensional systems. A striking mismatch of the $n\to\infty$ extrapolation of these simulations against analytical calculations is traced back to a breakdown of the identification of this limit with the spherical model.
Motivated by the results of two-dimensional conformal field theory (CFT) we investigate the finite-size scaling of the mass spectrum of an Ising model on three-dimensional lattices with a spherical cross section. Using a cluster-update Monte Carlo technique we find a linear relation between the masses and the corresponding scaling dimensions, in complete analogy to the situation in two dimensions. Amplitude ratios as well as the amplitudes themselves appear to be universal in this case.
We derive high-temperature series expansions for the free energy and susceptibility of the two-dimensional random-bond Ising model with a symmetric bimodal distribution of two positive coupling strengths J_1 and J_2 and study the influence of the quenched, random bond-disorder on the critical behavior of the model. By analysing the series expansions over a wide range of coupling ratios J_2/J_1, covering the crossover from weak to strong disorder, we obtain for the susceptibility with two different methods compelling evidence for a singularity of the form $χ\sim t^{-7/4} |\ln t|^{7/8}$, as predicted theoretically by Shalaev, Shankar, and Ludwig. For the specific heat our results are less convincing, but still compatible with the theoretically predicted log-log singularity.
Variational weak-coupling perturbation theory yields converging approximations, uniformly in the coupling strength. This allows us to calculate directly the coefficients of `strong-coupling' expansions. For the anharmonic oscillator we explain the physical origin of the empirically observed convergence behavior which is exponentially fast with superimposed oscillations.
We study 2D quantum gravity on spherical topologies using the Regge calculus approach with the $dl/l$ measure. Instead of a fixed non-regular triangulation which has been used before, we study for each system size four different random triangulations, which are obtained according to the standard Voronoi-Delaunay procedure. We compare both approaches quantitatively and show that the difference in the expectation value of $R^2$ between the fixed and the random triangulation depends on the lattice size and the surface area $A$. We also try again to measure the string susceptibility exponents through a finite-size scaling Ansatz in the expectation value of an added $R^2$ interaction term in an approach where $A$ is held fixed. The string susceptibility exponent $γ_{str}'$ is shown to agree with theoretical predictions for the sphere, whereas the estimate for $γ_{str}$ appears to be too negative.
We present a comparison of the performance of two non-local update algorithms for path integral Monte Carlo (PIMC) simulations, the multigrid Monte Carlo method and the staging algorithm. Looking at autocorrelation times for the internal energy we show that both refined algorithms beat the slowing down which is encountered for standard local update schemes in the continuum limit. We investigate the conditions under which the staging algorithm performs optimally and give a brief discussion of the mutual merits of the two algorithms.
We report a fairly detailed finite-size scaling analysis of the first-order phase transition in the three-dimensional 3-state Potts model on cubic lattices with emphasis on recently introduced quantities whose infinite-volume extrapolations are governed `only' by exponentially small terms. In these quantities no asymptotic power series in the inverse volume are involved which complicate the finite-size scaling behaviour of standard observables related to the specific-heat maxima or Binder-parameter minima. Introduced initially for strong first-order phase transitions in q-state Potts models with ``large enough'' q, the new techniques prove to be surprisingly accurate for a q value as small as 3. On the basis of the high-precision Monte Carlo data of Alves `et al.' [Phys. Rev. B43 (1991) 5846], this leads to a refined estimate of $β_t = 0.550,565(10)$ for the infinite-volume transition point.
We use Monte Carlo simulations to study pure 2D Euclidean quantum gravity with $R^2$-interaction on spherical topologies, employing Regge's formulation. We attempt to measure the string susceptibility exponent $γ_{\rm str}$ by using a finite-size scaling Ansatz in the expectation value of $R^2$, as has been done in a previous study by Bock and Vink ( hep-lat/9406018 ). By considerably extending the range and statistics of their study we find that this Ansatz is plagued by large systematic errors. The $R^2$ specific string susceptibility exponent $\GS'$ is found to agree with theoretical predictions, but its determination also is subject to large systematic errors and the presence of finite-size scaling corrections. To circumvent this obstacle we suggest a new scaling Ansatz which in principle should be able to predict both, $\GS$ and $\GS'$. First results indicate that this requires large system sizes to reduce the uncertainties in the finite-size scaling Ansätze. Nevertheless, our investigation shows that within the achievable accuracy the numerical estimates are still compatible with analytic predictions, contrary to the recent claim by Bock and Vink.
We report a high statistics simulation of Ising spins coupled to 2D quantum gravity in the Regge calculus approach using triangulated tori with up to $512^2$ vertices. For the constant area ensemble and the $dl/l$ functional measure we definitively can exclude the critical exponents of the Ising phase transition as predicted for dynamically triangulated surfaces. We rather find clear evidence that the critical exponents agree with the Onsager values for static regular lattices, independent of the coupling strength of an $R^2$ interaction term. For exploratory simulations using the lattice version of the Misner measure the situation is less clear.
Inspired by the multicanonical approach to simulations of first-order phase transitions we propose for $q$-state Potts models a combination of cluster updates with reweighting of the bond configurations in the Fortuin-Kastelein-Swendsen-Wang representation of this model. Numerical tests for the two-dimensional models with $q=7, 10$ and $20$ show that the autocorrelation times of this algorithm grow with the system size $V$ as $τ\propto V^α$, where the exponent takes the optimal random walk value of $α\approx 1$.
We report results of a Monte Carlo simulation of the $ϕ^4$ quantum chain. In order to enhance the efficiency of the simulation we combine multigrid simulation techniques with a refined discretization scheme. The resulting accuracy of our data allows for a significant test of an analytical approximation based on a variational ansatz. While the variational approximation is well reproduced for a large range of parameters we find significant deviations for low temperatures and large couplings.