We study the directional-ordering transition in the two-dimensional classical and quantum compass models on the square lattice by means of Monte Carlo simulations. An improved algorithm is presented which builds on the Wolff cluster algorithm in one-dimensional subspaces of the configuration space. This improvement allows us to study classical systems up to L=512. Based on this algorithm, we give evidence for the presence of strongly anomalous scaling for periodic boundary conditions which is much worse than anticipated before. We propose and study alternative boundary conditions for the compass model which do not make use of extended configuration spaces and show that they completely remove the problem with finite-size scaling. In the last part, we apply these boundary conditions to the quantum problem and present a considerably improved estimate for the critical temperature which should be of interest for future studies on the compass model. Our investigation identifies a strong one-dimensional magnetic ordering tendency with a large correlation length as the cause of the unusual scaling and moreover allows for a precise quantification of the anomalous length scale involved.
We study the directional-ordering transition in the two-dimensional classical and quantum compass models on the square lattice by means of Monte Carlo simulations. An improved algorithm is presented which builds on the Wolff cluster algorithm in one-dimensional subspaces of the configuration space. This improvement allows us to study classical systems up to L=512. Based on this algorithm, we give evidence for the presence of strongly anomalous scaling for periodic boundary conditions which is much worse than anticipated before. We propose and study alternative boundary conditions for the compass model which do not make use of extended configuration spaces and show that they completely remove the problem with finite-size scaling. In the last part, we apply these boundary conditions to the quantum problem and present a considerably improved estimate for the critical temperature which should be of interest for future studies on the compass model. Our investigation identifies a strong one-dimensional magnetic ordering tendency with a large correlation length as the cause of the unusual scaling and moreover allows for a precise quantification of the anomalous length scale involved.
A comprehensive study of the two-dimensional (2D) compass model on the square lattice is performed for classical and quantum spin degrees of freedom using Monte Carlo and quantum Monte Carlo methods. We employ state-of-the-art implementations using Metropolis, stochastic series expansion, and parallel tempering techniques to obtain the critical ordering temperatures and critical exponents. In a preinvestigation we reconsider the classical compass model where we study and contrast the finite-size scaling behavior of ordinary periodic boundary conditions against annealed boundary conditions. It is shown that periodic boundary conditions suffer from extreme finite-size effects which might be caused by closed-loop excitations on the torus. These excitations also appear to have severe effects on the Binder parameter. On this footing we report on a systematic Monte Carlo study of the quantum compass model. Our numerical results are at odds with recent literature on the subject which we trace back to neglecting the strong finite-size effects on periodic lattices. The critical temperatures are obtained as Tc=0.1464(2)J and Tc=0.055(1)J for the classical and quantum versions, respectively, and our data support a transition in the 2D Ising universality class for both cases.
We determine the finite-temperature phase diagram and critical behavior of the classical square-lattice Heisenberg-compass model using large-scale Monte Carlo simulations and finite-size scaling. Six symmetry distinct ordered phases are identified. The four phases that simultaneously break the spin-lattice C4 and in-plane spin-inversion symmetries undergo continuous transitions in the Ashkin-Teller universality class, with the associated critical lines terminating at four-state Potts points, beyond which the transitions become first order. In contrast, the two z-polarized phases display conventional two-dimensional Ising criticality. Our results reveal how the interplay between Heisenberg exchange and compass anisotropy organizes these distinct critical regimes, thereby completing the characterization of the model's thermal phase transitions.
Abstract. We obtained exact heat capacities of the quantum compass model on the square L×L clusters with L = 2, 3,4, 5 using Kernel Polynomial Method and compare them with heat capacity of a large compass ladder. Intersite correlations found in the ground state for these systems demonstrate that the quantum compass model differs from its classical version. Quantum compass model originating from the orbital (Kugel-Khomskii) superexchange in transition metal oxides has been recently studied both using analytical [1] and numerical [2, 3] approach. In spite of remarkable progress in numerical methods for two–dimensional (2D) Ising– like models [4, 5], exact solutions are necessary to investigate the structure of excited states. The one–dimensional (1D) compass model is exactly solvable by mapping to quantum Ising model [6] and exhibits interesting properties while approaching to the quantum critical point at zero temperature [7]. In addition, its ladder version, first considered by Douçot et al. [1] is solvable in a similar way and its partition function [8] can be obtained exactly in case of a large (but finite) system. There is no exact solution for the 2D compass model but the latest Monte Carlo data [3] prove that the model exhibits a phase transition at finite temperature both in quantum and classical version with symmetry breaking between x and z part of the Hamiltonian. This
We reconsider the tricritical Blume-Capel model on the square lattice with a magnetic field acting on the open boundaries in one direction. Periodic boundary conditions are applied in the other direction. We apply three types of Monte Carlo algorithms, local Metropolis updates, and cluster algorithms of the Wolff and geometric type, adapted to the symmetry properties of the model. Statistical analyses of the bulk magnetization, the bulk Binder ratio, the edge magnetization, and the connected product of the edge and bulk magnetizations lead to new results confirming the presence of a singular edge transition at H_{sc}≈0.68, as we reported earlier [Phys. Rev. E 71, 026109 (2005)PLEEE81539-375510.1103/PhysRevE.71.026109]. We provide a plausible answer concerning a discrepancy between the behavior of the edge Binder ratio reported in that work and our new results.
To demonstrate the application of computational physics in modern statistical physics, this paper investigates the critical behavior of a finite-size ferromagnetic system based on the two-dimensional square-lattice Ising model.A nearest-neighbor Ising Hamiltonian is constructed, and periodic boundary conditions are adopted to reduce finite-boundary effects.The Metropolis Monte Carlo algorithm is used to sample spin configurations at different temperatures, while energy, magnetization, specific heat, magnetic susceptibility and Binder cumulant are calculated to describe the transition from an ordered state to a disordered state.The simulation results show that the system exhibits clear long-range order at low temperatures, where the order parameter remains high.As temperature increases, thermal fluctuations become stronger, magnetization decreases rapidly, and the system gradually evolves into a disordered state.The specific heat curve shows a pronounced peak at approximately (T = 2.269), while the susceptibility peak shifts slightly toward the high-temperature side due to finite-size effects and the discrete temperature grid.Overall, the simulated critical region is consistent with the theoretical critical temperature of the two-dimensional square-lattice Ising model.Spin-configuration analysis further reveals enhanced fluctuations and domain reconstruction near the critical region.In addition, a thermalization criterion and an extended sampling test for L = 48 near the critical region are added to examine the influence of finite sampling and equilibration.These results indicate that Monte Carlo simulation of the two-dimensional Ising model can effectively reproduce continuous phase transition, critical fluctuation and finite-size effects, providing a basic computational case for understanding critical phenomena in complex systems and data-driven physical analysis.
We give an intuitive geometric explanation for the apparent breakdown of standard finite-size scaling in systems with periodic boundaries above the upper critical dimension. The Ising model and self-avoiding walk are simulated on five-dimensional hypercubic lattices with free and periodic boundary conditions, by using geometric representations and recently introduced Markov-chain Monte Carlo algorithms. We show that previously observed anomalous behavior for correlation functions, measured on the standard Euclidean scale, can be removed by defining correlation functions on a scale which correctly accounts for windings.