Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userClock Monte Carlo methods
82
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We propose the clock Monte Carlo technique for sampling each successive chain step in constant time. It is built on a recently proposed factorized transition filter and its core features include its O(1) computational complexity and its generality. We elaborate how it leads to the clock factorized Metropolis (clock FMet) method, and discuss its application in other update schemes. By grouping interaction terms into boxes of tunable sizes, we further formulate a variant of the clock FMet algorithm, with the limiting case of a single box reducing to the standard Metropolis method. A theoretical analysis shows that an overall acceleration of ${rm O}(N^kappa)$ ($0 ! leq ! kappa ! leq ! 1$) can be achieved compared to the Metropolis method, where $N$ is the system size and the $kappa$ value depends on the nature of the energy extensivity. As a systematic test, we simulate long-range O$(n)$ spin models in a wide parameter regime: for $n ! = ! 1,2,3$, with disordered algebraically decaying or oscillatory Ruderman-Kittel-Kasuya-Yoshida-type interactions and with and without external fields, and in spatial dimensions from $d ! = ! 1, 2, 3$ to mean-field. The O(1) computational complexity is demonstrated, and the expected acceleration is confirmed. Its flexibility and its independence from the interaction range guarantee that the clock method would find decisive applications in systems with many interaction terms.
rate research
Read More
In this work, we discuss the implications of a recently obtained equilibrium fluctuation-dissipation relation on the extension of the available Monte Carlo methods based on the consideration of the Gibbs canonical ensemble to account for the existence of an anomalous regime with negative heat capacities $C<0$. The resulting framework appears as a suitable generalization of the methodology associated with the so-called textit{dynamical ensemble}, which is applied to the extension of two well-known Monte Carlo methods: the Metropolis importance sample and the Swendsen-Wang clusters algorithm. These Monte Carlo algorithms are employed to study the anomalous thermodynamic behavior of the Potts models with many spin states $q$ defined on a $d$-dimensional hypercubic lattice with periodic boundary conditions, which successfully reduce the exponential divergence of decorrelation time $tau$ with the increase of the system size $N$ to a weak power-law divergence $taupropto N^{alpha}$ with $alphaapprox0.2$ for the particular case of the 2D 10-state Potts model.
Previously, we have presented a methodology to extend canonical Monte Carlo methods inspired on a suitable extension of the canonical fluctuation relation $C=beta^{2}<delta E^{2}>$ compatible with negative heat capacities $C<0$. Now, we improve this methodology by introducing a better treatment of finite size effects affecting the precision of a direct determination of the microcanonical caloric curve $beta (E) =partial S(E) /partial E$, as well as a better implementation of MC schemes. We shall show that despite the modifications considered, the extended canonical MC methods possibility an impressive overcome of the so-called textit{super-critical slowing down} observed close to the region of a temperature driven first-order phase transition. In this case, the dependence of the decorrelation time $tau$ with the system size $N$ is reduced from an exponential growth to a weak power-law behavior $tau(N)propto N^{alpha}$, which is shown in the particular case of the 2D seven-state Potts model where the exponent $alpha=0.14-0.18$.
The Widom-Rowlinson model of a fluid mixture is studied using a new cluster algorithm that is a generalization of the invaded cluster algorithm previously applied to Potts models. Our estimate of the critical exponents for the two-component fluid are consistent with the Ising universality class in two and three dimensions. We also present results for the three-component fluid.
We present a Monte Carlo method for the direct evaluation of the difference between the free energies of two crystal structures. The method is built on a lattice-switch transformation that maps a configuration of one structure onto a candidate configuration of the other by `switching one set of lattice vectors for the other, while keeping the displacements with respect to the lattice sites constant. The sampling of the displacement configurations is biased, multicanonically, to favor paths leading to `gateway arrangements for which the Monte Carlo switch to the candidate configuration will be accepted. The configurations of both structures can then be efficiently sampled in a single process, and the difference between their free energies evaluated from their measured probabilities. We explore and exploit the method in the context of extensive studies of systems of hard spheres. We show that the efficiency of the method is controlled by the extent to which the switch conserves correlated microstructure. We also show how, microscopically, the procedure works: the system finds gateway arrangements which fulfill the sampling bias intelligently. We establish, with high precision, the differences between the free energies of the two close packed structures (fcc and hcp) in both the constant density and the constant pressure ensembles.
Fractional derivatives are nonlocal differential operators of real order that often appear in models of anomalous diffusion and a variety of nonlocal phenomena. Recently, a version of the Schrodinger Equation containing a fractional Laplacian has been proposed. In this work, we develop a Fractional Path Integral Monte Carlo algorithm that can be used to study the finite temperature behavior of the time-independent Fractional Schrodinger Equation for a variety of potentials. In so doing, we derive an analytic form for the finite temperature fractional free particle density matrix and demonstrate how it can be sampled to acquire new sets of particle positions. We employ this algorithm to simulate both the free particle and $^{4}$He (Aziz) Hamiltonians. We find that the fractional Laplacian strongly encourages particle delocalization, even in the presence of interactions, suggesting that fractional Hamiltonians may manifest atypical forms of condensation. Our work opens the door to studying fractional Hamiltonians with arbitrarily complex potentials that escape analytical solutions.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
