Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDirected-loop Monte Carlo simulations of vertex models
55
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We show how the directed-loop Monte Carlo algorithm can be applied to study vertex models. The algorithm is employed to calculate the arrow polarization in the six-vertex model with the domain wall boundary conditions (DWBC). The model exhibits spatially separated ordered and ``disordered regions. We show how the boundary between these regions depends on parameters of the model. We give some predictions on the behavior of the polarization in the thermodynamic limit and discuss the relation to the Arctic Circle theorem.
rate research
Read More
Population annealing is a recent addition to the arsenal of the practitioner in computer simulations in statistical physics and beyond that is found to deal well with systems with complex free-energy landscapes. Above all else, it promises to deliver unrivaled parallel scaling qualities, being suitable for parallel machines of the biggest calibre. Here we study population annealing using as the main example the two-dimensional Ising model which allows for particularly clean comparisons due to the available exact results and the wealth of published simulational studies employing other approaches. We analyze in depth the accuracy and precision of the method, highlighting its relation to older techniques such as simulated annealing and thermodynamic integration. We introduce intrinsic approaches for the analysis of statistical and systematic errors, and provide a detailed picture of the dependence of such errors on the simulation parameters. The results are benchmarked against canonical and parallel tempering simulations.
We use Monte Carlo simulations to study ${rm Ni Fe_2O_4}$ nanoparticles. Finite size and surface effects differentiate them from their bulk counterparts. A continuous version of the Wang-Landau algorithm is used to calculate the joint density of states $g(M_z, E)$ efficiently. From $g(M_z, E)$, we obtain the Bragg-Williams free energy of the particle, and other physical quantities. The hysteresis is observed when the nanoparticles have both surface disorder and surface anisotropy. We found that the finite coercivity is the result of interplay between surface disorder and surface anisotropy. If the surface disorder is absent or the surface anisotropy is relatively weak, the nanoparticles often exhibit superparamagnetism.
The unconstrained ensemble describes completely open systems whose control parameters are chemical potential, pressure, and temperature. For macroscopic systems with short-range interactions, thermodynamics prevents the simultaneous use of these intensive variables as control parameters, because they are not independent and cannot account for the system size. When the range of the interactions is comparable with the size of the system, however, these variables are not truly intensive and may become independent, so equilibrium states defined by the values of these parameters may exist. Here, we derive a Monte Carlo algorithm for the unconstrained ensemble and show that simulations can be performed using chemical potential, pressure, and temperature as control parameters. We illustrate the algorithm by applying it to physical systems where either the system has long-range interactions or is confined by external conditions. The method opens up a new avenue for the simulation of completely open systems exchanging heat, work, and matter with the environment.
We propose a new generalized-ensemble algorithm, which we refer to as the multibaric-multithermal Monte Carlo method. The multibaric-multithermal Monte Carlo simulations perform random walks widely both in volume space and in potential energy space. From only one simulation run, one can calculate isobaric-isothermal-ensemble averages at any pressure and any temperature. We test the effectiveness of this algorithm by applying it to the Lennard-Jones 12-6 potential system with 500 particles. It is found that a single simulation of the new method indeed gives accurate average quantities in isobaric-isothermal ensemble for a wide range of pressure and temperature.
We present Quantum Monte Carlo simulations of a generalization of the Feynman-Kikuchi model which includes the possibility of vacancies and interactions between the particles undergoing exchange. By measuring the winding number (superfluid density) and density structure factor, we determine the phase diagram, and show that it exhibits regions which possess both superfluid and charge ordering.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
