Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userTethered Monte Carlo: computing the effective potential without critical slowing down
162
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present Tethered Monte Carlo, a simple, general purpose method of computing the effective potential of the order parameter (Helmholtz free energy). This formalism is based on a new statistical ensemble, closely related to the micromagnetic one, but with an extended configuration space (through Creutz-like demons). Canonical averages for arbitrary values of the external magnetic field are computed without additional simulations. The method is put to work in the two dimensional Ising model, where the existence of exact results enables us to perform high precision checks. A rather peculiar feature of our implementation, which employs a local Metropolis algorithm, is the total absence, within errors, of critical slowing down for magnetic observables. Indeed, high accuracy results are presented for lattices as large as L=1024.
rate research
Read More
We study the performance of Monte Carlo simulations that sample a broad histogram in energy by determining the mean first-passage time to span the entire energy space of d-dimensional ferromagnetic Ising/Potts models. We first show that flat-histogram Monte Carlo methods with single-spin flip updates such as the Wang-Landau algorithm or the multicanonical method perform sub-optimally in comparison to an unbiased Markovian random walk in energy space. For the d=1,2,3 Ising model, the mean first-passage time tau scales with the number of spins N=L^d as tau propto N^2L^z. The critical exponent z is found to decrease as the dimensionality d is increased. In the mean-field limit of infinite dimensions we find that z vanishes up to logarithmic corrections. We then demonstrate how the slowdown characterized by z>0 for finite d can be overcome by two complementary approaches - cluster dynamics in connection with Wang-Landau sampling and the recently developed ensemble optimization technique. Both approaches are found to improve the random walk in energy space so that tau propto N^2 up to logarithmic corrections for the d=1 and d=2 Ising model.
We study the three-dimensional Ising model at the critical point in the fixed-magnetization ensemble, by means of the recently developed geometric cluster Monte Carlo algorithm. We define a magnetic-field-like quantity in terms of microscopic spin-up and spin-down probabilities in a given configuration of neighbors. In the thermodynamic limit, the relation between this field and the magnetization reduces to the canonical relation M(h). However, for finite systems, the relation is different. We establish a close connection between this relation and the probability distribution of the magnetization of a finite-size system in the canonical ensemble.
The confinement of critical fluctuations in soft media induces critical Casimir forces acting on the confining surfaces. The temperature and geometry dependences of such forces are characterized by universal scaling functions. A novel approach is presented to determine them for films via Monte Carlo simulations of lattice models. The method is based on an integration scheme of free energy differences. Our results for the Ising and the XY universality class compare favourably with corresponding experimental results for wetting layers of classical binary liquid mixtures and of 4He, respectively.
We investigate the interaction potential of superconducting vortices at the full quantum level. We formulate the interaction potential in a constrained path integral and calculate it by the quantum Monte Carlo simulation. The vortex-vortex potential is attractive (type-I), repulsive (type-II), and flat (critical) depending on a coupling constant. The vortex-antivortex potential also depends on the coupling constant at long range but is always attractive at short range.
A path integral Monte Carlo method based on the worm algorithm has been developed to compute the chemical potential of interacting bosonic quantum fluids. By applying it to finite-sized systems of helium-4 atoms, we have confirmed that the chemical potential scales inversely with the number of particles to lowest order. The introduction of a simple scaling form allows for the extrapolation of the chemical potential to the thermodynamic limit, where we observe excellent agreement with known experimental results for helium-4 at saturated vapor pressure. We speculate on future applications of the proposed technique, including its use in studies of confined quantum fluids.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
