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.
Systems of particles in a confining potential exhibit a spatially dependent density which fundamentally alters the nature of phase transitions that occur. A specific instance of this situation, which is being extensively explored currently, concerns the properties of ultra-cold, optically trapped atoms. Of interest is how the superfluid-insulator transition is modified by the inhomogeneity, and, indeed, the extent to which a sharp transition survives at all. This paper explores a classical analog of these systems, the Blume-Capel model with a spatially varying single ion anisotropy and/or temperature gradient. We present results both for the nature of the critical properties and for the validity of the local density approximation which is often used to model the inhomogeneous case. We compare situations when the underlying uniform transition is first and second order.
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.
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.
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.