Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSampling from a polytope and hard-disk Monte Carlo
143
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The hard-disk problem, the statics and the dynamics of equal two-dimensional hard spheres in a periodic box, has had a profound influence on statistical and computational physics. Markov-chain Monte Carlo and molecular dynamics were first discussed for this model. Here we reformulate hard-disk Monte Carlo algorithms in terms of another classic problem, namely the sampling from a polytope. Local Markov-chain Monte Carlo, as proposed by Metropolis et al. in 1953, appears as a sequence of random walks in high-dimensional polytopes, while the moves of the more powerful event-chain algorithm correspond to molecular dynamics evolution. We determine the convergence properties of Monte Carlo methods in a special invariant polytope associated with hard-disk configurations, and the implications for convergence of hard-disk sampling. Finally, we discuss parallelization strategies for event-chain Monte Carlo and present results for a multicore implementation.
rate research
Read More
The diagrammatic Monte Carlo (Diag-MC) method is a numerical technique which samples the entire diagrammatic series of the Greens function in quantum many-body systems. In this work, we incorporate the flat histogram principle in the diagrammatic Monte method and we term the improved version Flat Histogram Diagrammatic Monte Carlo method. We demonstrate the superiority of the method over the standard Diag-MC in extracting the long-imaginary-time behavior of the Greens function, without incorporating any a priori knowledge about this function, by applying the technique to the polaron problem
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 construct asymptotic arguments for the relative efficiency of rejection-free Monte Carlo (MC) methods compared to the standard MC method. We find that the efficiency is proportional to $exp{({const} beta)}$ in the Ising, $sqrt{beta}$ in the classical XY, and $beta$ in the classical Heisenberg spin systems with inverse temperature $beta$, regardless of the dimension. The efficiency in hard particle systems is also obtained, and found to be proportional to $(rhoc -rho)^{-d}$ with the closest packing density $rhoc$, density $rho$, and dimension $d$ of the systems. We construct and implement a rejection-free Monte Carlo method for the hard-disk system. The RFMC has a greater computational efficiency at high densities, and the density dependence of the efficiency is as predicted by our arguments.
The fidelity susceptibility is a general purpose probe of phase transitions. With its origin in quantum information and in the differential geometry perspective of quantum states, the fidelity susceptibility can indicate the presence of a phase transition without prior knowledge of the local order parameter, as well as reveal the universal properties of a critical point. The wide applicability of the fidelity susceptibility to quantum many-body systems is, however, hindered by the limited computational tools to evaluate it. We present a generic, efficient, and elegant approach to compute the fidelity susceptibility of correlated fermions, bosons, and quantum spin systems in a broad range of quantum Monte Carlo methods. It can be applied both to the ground-state and non-zero temperature cases. The Monte Carlo estimator has a simple yet universal form, which can be efficiently evaluated in simulations. We demonstrate the power of this approach with applications to the Bose-Hubbard model, the spin-$1/2$ XXZ model, and use it to examine the hypothetical intermediate spin-liquid phase in the Hubbard model on the honeycomb lattice.
We describe a Monte Carlo procedure which allows sampling of the disjoint configuration spaces associated with crystalline and fluid phases, within a single simulation. The method utilises biased sampling techniques to enhance the probabilities of gateway states (in each phase) which are such that a global switch (to the other phase) can be implemented. Equilibrium freezing-point parameters can be determined directly; statistical uncertainties prescribed transparently; and finite-size effects quantified systematically. The method is potentially quite general; we apply it to the freezing of hard spheres.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
