Related papers: Monte Carlo simulation of ice models
A $\theta$ term, which couples to topological charge, is added to the two-dimensional lattice CP^3 model and U(1) gauge theory. Monte Carlo simulations are performed and compared to strong-coupling character expansions. In certain…
The cactus approximation in the cluster variation method is applied to the spin ice system with nearest neighbor ferromagnetic coupling. The temperature dependences of the entropy and the specific heat show qualitatively good agreement with…
We introduce a Monte Carlo method, as a modification of existing cluster algorithms, which allows simulations directly on systems of infinite size, and for quantum models also at beta=infinity. All two-point functions can be obtained,…
The critical behavior of the disordered ferromagnetic Ising model is studied numerically by the Monte Carlo method in a wide range of variation of concentration of nonmagnetic impurity atoms. The temperature dependences of correlation…
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…
High-order virtual excitations play an important role in microscopic models of nuclear reactions at intermediate energies. However, the factorial growth of their complexity has prevented their consistent inclusion in ab initio many-body…
We created an efficient algorithm suitable for graphics processing units (GPUs) to perform Monte Carlo simulations of a subset of reaction-diffusion models. The algorithm uses techniques that are specific to GPU programming, and combines…
A simple site-diluted classical XY model is proposed to study the magnetic properties of the ferromagnetic polycrystalline Fe[Se$_2$CN(C$_2$H$_5$)$_2$]$_2$Cl diluted with diamagnetic Zn[S$_2$CN(C$_2$H$_5$)$_2$]$_2$. An extra superexchange…
Taking the two-dimensional Ising model for example, short-time behavior of critical dynamics with a conserved order parameter is investigated by Monte Carlo simulations. Scaling behavior is observed, but the dynamic exponent $z$ is updating…
In finite-size scaling analyses of Monte Carlo simulations of second-order phase transitions one often needs an extended temperature/energy range around the critical point. By combining the replica-exchange algorithm with cluster updates…
The computational cost of a Monte Carlo algorithm can only be meaningfully discussed when taking into account the magnitude of the resulting statistical error. Aiming for a fixed error per particle, we study the scaling behavior of the…
Monte Carlo methods play important part in modern statistical physics. The application of these methods suffer from two main difficulties.The first is caused by the relatively small number of particles that can participate in any numerical…
Ultracold neutrons (UCN) with kinetic energies up to 300 neV can be stored in material or magnetic confinements for hundreds of seconds. This makes them a very useful tool for probing fundamental symmetries of nature, by searching for…
We simulate TIP3P water using a constrained Monte Carlo algorithm to generate electrostatic interactions eliminating the need to sum over long ranged Coulomb interactions. We study discretization errors when interpolating charges using…
We propose an efficient Monte Carlo algorithm for the off-lattice simulation of dense hard sphere polymer melts using cluster moves, called event chains, which allow for a rejection-free treatment of the excluded volume. Event chains also…
Quantum Monte Carlo (QMC) simulations constitute nowadays one of the most powerful methods to study strongly correlated quantum systems, provided that no "sign problem" arises. However, many systems of interest, including highly frustrated…
We put forward a simple procedure for extracting dynamical information from Monte Carlo simulations, by appropriate matching of the short-time diffusion tensor with its infinite-dilution limit counterpart, which is supposed to be known.…
If a stochastic system during some periods of its evolution can be divided into non-interacting parts, the kinetics of each part can be simulated independently. We show that this can be used in the development of efficient Monte Carlo…
On the base of a Feynman-Kac--type formula involving Poisson stochastic processes, recently a Monte Carlo algorithm has been introduced, which describes exactly the real- or imaginary-time evolution of many-body lattice quantum systems. We…
A Monte Carlo algorithm to simulate the isothermal recrystallization process of snow is presented. The snow metamorphism is approximated by two mass redistribution processes, surface diffusion and sublimation-deposition. The algorithm is…