Related papers: Strengths and Weaknesses of Parallel Tempering
We present a new Monte Carlo algorithm that allows the simultaneous determination of a few extremal eigenpairs of a very large matrix. It extends the power method and uses a new sampling method, the sewing method, that does a large state…
We have developed a new simulation algorithm for free-energy calculations. The method is a multidimensional extension of the replica-exchange method. While pairs of replicas with different temperatures are exchanged during the simulation in…
Planar thermal equilibration is studied using direct numerical simulations of ultracold two-dimensional (2D) ion crystals in a Penning trap with a rotating wall. The large magnetic field of the trap splits the modes that describe in-plane…
The replica exchange method is a powerful tool for overcoming slow relaxation in molecular simulations, but its efficiency depends strongly on the choice of the number and interval of replicas and their exchange probabilities. Here, we…
When auditing a redistricting plan, a persuasive method is to compare the plan with an ensemble of neutrally drawn redistricting plans. Ensembles are generated via algorithms that sample distributions on balanced graph partitions. To audit…
We consider the nonequilibrium dynamics of a a real scalar field in a degenerate double-well potential. The system is prepared in the lowest free energy state in one of the wells and the dynamics is driven by the coupling of the field to a…
Combinatorial optimization problems play crucial roles in real-world applications, and many studies from a physics perspective have contributed to specialized hardware for high-speed computation. However, some combinatorial optimization…
Umbrella sampling is an efficient method for the calculation of free energy changes of a system along well-defined reaction coordinates. However, when multiple parallel channels along the reaction coordinate or hidden barriers in directions…
By extending the concept of diffusion to the potential energy landscapes (PELs), we introduce the mean-squared energy difference (MSED) as a novel quantity to investigate the intrinsic properties of glass. MSED can provide a clear…
It is analyzed whether the potential energy landscape of a glass-forming system can be effectively mapped on a random model which is described in statistical terms. For this purpose we generalize the simple trap model of Bouchaud and…
We present a quantitative comparison between extensive Monte Carlo simulations and self-consistent field calculations on the phase diagram and wetting behavior of a symmetric, binary (AB) polymer blend confined into a film. The flat walls…
Monte Carlo simulation provides a powerful tool for understanding and exploring thermodynamic phase equilibria in many-particle interacting systems. Among the most physically intuitive simulation methods is Gibbs ensemble Monte Carlo…
We study a system of $y=2$ coupled copies of a well-known constraint satisfaction problem (random hypergraph bicoloring) to examine how the ferromagnetic coupling between the copies affects the properties of the solution space. We solve the…
A rich literature has been produced on the quantum states of atoms and molecules confined into infinite potential wells with a specified symmetry. Apart from their interest as basic quantum systems, confined atoms and molecules are useful…
We represent the slow, glassy equilibrium dynamics of a line in a two-dimensional random potential landscape as driven by an array of asymptotically independent two-state systems, or loops, fluctuating on all length scales. The assumption…
Monte Carlo sampling methods often suffer from long correlation times. Consequently, these methods must be run for many steps to generate an independent sample. In this paper a method is proposed to overcome this difficulty. The method…
We propose an algorithm for molecular dynamics or Monte Carlo simulations that uses an interpolation procedure to estimate potential energy values from energies and gradients evaluated previously at points of a simplicial mesh. We chose an…
When approaching the continuum limit in lattice QCD or other theories in a setup with topological sectors, conventional update algorithms experience a particularly severe form of critical slowing down that is caused by high action barriers…
A new Monte Carlo algorithm for 2-dimensional spin glasses is presented. The use of clusters makes possible global updates and leads to a gain in speed of several orders of magnitude. As an example, we study the 2-dimensional +/-J…
We present a massive equilibrium simulation of the three-dimensional Ising spin glass at low temperatures. The Janus special-purpose computer has allowed us to equilibrate, using parallel tempering, L=32 lattices down to T=0.64 Tc. We…