Related papers: Large-scale Monte Carlo simulation of two-dimensio…
For a second-order phase transition the critical energy range of interest is larger than the energy range covered by a canonical Monte Carlo simulation at the critical temperature. Such an extended energy range can be covered by performing…
We study the dynamics of a classical, two-component plasma in two dimensions, in the vicinity of the Kosterlitz-Thouless (KT) transition where the system passes from a dielectric low-temperature phase (consisting of bound pairs) to a…
We introduce a method that ensures efficient computation of one-dimensional quantum systems with long-range interactions across all temperatures. Our algorithm operates within a quasi-polynomial runtime for inverse temperatures up to…
The three-dimensional anisotropic classical XY ferromagnet has been investigated by extensive Monte Carlo simulation using the Metropolis single spin flip algorithm. The magnetization ($M$) and the susceptibility ($\chi$) are measured and…
A model based on microcanonical Monte Carlo method is used to study the application of the temperature gradient along a two-dimensional (2D) Ising system. We estimate the system size effects on thermal conductivity, $K$, for a nano-scale…
An extensive equilibrium Monte Carlo simulation is performed on the 3D isotropic Heisenberg SG model with the random nearest-neighbor Gaussian coupling, with particular interest in its chiral-glass (CG) and spin-glass (SG) orderings. For…
Curie temperature and exponents are studied for the three-dimensional double-exchange model. Applying the O(N) Monte Carlo algorithm, we perform systematic finite-size scaling analyses on the data up to $20^3$ sites. The obtained values of…
Understanding how interacting particles approach thermal equilibrium is a major challenge of quantum simulators. Unlocking the full potential of such systems toward this goal requires flexible initial state preparation, precise time…
It is by now well-known that ground states of gapped one-dimensional (1d) quantum-many body systems with short-range interactions can be studied efficiently using classical computers and matrix product state techniques. A corresponding…
We develop a quantum Monte Carlo procedure to compute the Renyi mutual information of an interacting quantum many-body system at non-zero temperature. Performing simulations on a spin-1/2 XXZ model, we observe that for a subregion of fixed…
We study the antiferromagnetic {\it XY} model on a triangular lattice by extensive Monte Carlo simulations, focusing on its ordering and critical properties. Our result clearly shows that two separate transitions occur at two distinct…
Berezinskii-Kosterlitz-Thouless transition of the classical XY model is re-investigated, combining the Tensor Network Renormalization (TNR) and the Level Spectroscopy method based on the finite-size scaling of the Conformal Field Theory. By…
Recently, Velazquez and Curilef have proposed a methodology to extend Monte Carlo algorithms based on canonical ensemble, which is aimed to overcome slow sampling problems associated with temperature-driven discontinuous phase transitions.…
A modern graphics processing unit (GPU) is able to perform massively parallel scientific computations at low cost. We extend our implementation of the checkerboard algorithm for the two dimensional Ising model [T. Preis et al., J. Comp.…
We introduce an algorithm to systematically improve the efficiency of parallel tempering Monte Carlo simulations by optimizing the simulated temperature set. Our approach is closely related to a recently introduced adaptive algorithm that…
We consider quenches in non-conserved two-dimensional XY systems between any two temperatures below the Kosterlitz-Thouless transition. The evolving systems are defect free at coarse-grained scales, and can be exactly treated. Correlations…
Using Monte Carlo simulations, we systematically investigate the non-equilibrium dynamics of the chiral degree of freedom in the two-dimensional fully frustrated XY model. The critical initial increase of the staggered chiral magnetization…
Using Monte Carlo methods, we compute the finite-size scaling function of the helicity modulus $\Upsilon$ of the two-dimensional O(3) model and compare it to the low temperature expansion prediction. From this, we estimate the range of…
We develop a trap-size scaling theory for trapped particle systems at quantum transitions. As a theoretical laboratory, we consider a quantum XY chain in an external transverse field acting as a trap for the spinless fermions of its…
We adopt the neural network flow (NN flow) method to study the Berezinskii-Kosterlitz-Thouless (BKT) phase transitions of the 2-dimensional q-state clock model with $q\ge 4$. The NN flow consists of a sequence of the same units to proceed…