BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
Wolfhard Janke contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph 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$.