We use the Monte Carlo method to study an antiferromagnetical Ising spin system on a centred honeycomb lattice, which is composed of two kinds of 1/2 spin particles A and B. There exist two different bond energies JA–A and JA–B in this lattice. Our study is focused on how the ratio of JA–B to JA–A influences the critical behaviour of this system by analysing the physical quantities, such as the energy, the order parameter, the specific heat, susceptibility, etc each as a function of temperature for a given ratio of JA–B to JA–A. Using these results together with the finite-size scaling method, we obtain a phase diagram for the ratio JA–B/JA–A. This work is helpful for studying the phase transition problem of crystals composed of compounds.
We use the Monte Carlo method to study an antiferromagnetical Ising spin system on a centred honeycomb lattice, which is composed of two kinds of 1/2 spin particles A and B. There exist two different bond energies JA–A and JA–B in this lattice. Our study is focused on how the ratio of JA–B to JA–A influences the critical behaviour of this system by analysing the physical quantities, such as the energy, the order parameter, the specific heat, susceptibility, etc each as a function of temperature for a given ratio of JA–B to JA–A. Using these results together with the finite-size scaling method, we obtain a phase diagram for the ratio JA–B/JA–A. This work is helpful for studying the phase transition problem of crystals composed of compounds.
Abstract. We use improved Monte-Carlo algorithms to study the antiferromagnetic 2D-Ising model with competing interactions J1 on nearest neighbour and J2 on next-nearest neighbour bonds. The finite-temperature phase diagram is divided by a critical point at J2 = J1/2 where the groundstate is highly degenerate. To analyse the phase boundaries we look at the specific heat and the energy distribution for various ratios of J2/J1. We find a first order transition for small J2> J1/2 and the transition temperature suppressed to TC = 0 at the critical point. 1.
We use improved Monte-Carlo algorithms to study the antiferromagnetic 2D-Ising model with competing interactions J 1 on nearest neighbour and J 2 on next-nearest neighbour bonds. The finite-temperature phase diagram is divided by a critical point at J 2 = J 1 /2 where the groundstate is highly degenerate. To analyse the phase boundaries we look at the speciflc heat and the energy distribution for various ratios of J 2 / 1 . We find a first order transition for small J 2 > J 1 /2 and the transition temperature suppressed to T C = 0 at the critical point.
The quantum phase transition, scaling behaviors, and thermodynamics in the spin-1/2 quantum Heisenberg model with antiferromagnetic coupling J>0 in the armchair direction and ferromagnetic interaction J^{'}<0 in the zigzag direction on a honeycomb lattice are systematically studied using the continuous-time quantum Monte Carlo method. By calculating the Binder ratio Q_{2} and spin stiffness ρ in two directions for various coupling ratios α=J^{'}/J under different lattice sizes, we found that a quantum phase transition from the dimerized phase to the stripe phase occurs at the quantum critical point α_{c}=-0.93. Through the finite-size scaling analysis on Q_{2}, ρ_{x}, and ρ_{y}, we determined the critical exponent related to the correlation length ν to be 0.7212(8), implying that this transition falls into a classical Heisenberg O(3) universality. A zero magnetization plateau is observed in the dimerized phase, whose width decreases with increasing α. A phase diagram in the coupling ratio α-magnetic field h plane is obtained, where four phases, including dimerized, stripe, canted stripe, and polarized, are identified. It is also unveiled that the temperature dependence of the specific heat C(T) for different α's intersects precisely at one point, similar to that of liquid ^{3}He under different pressures and several magnetic compounds under various magnetic fields. The scaling behaviors of Q_{2}, ρ, and C(T) are carefully analyzed. The susceptibility is compared with the experimental data to give the magnetic parameters of both compounds.
Monte Carlo simulations of the Ising model on the triangular lattice with the antiferromagnetic nn interaction J and the ferromagnetic nnn interaction J ' are carried out. Several physical quantities (magnetizations, energies, susceptibilities, specific heats of sublattices and of the total system) are calculated. A phase which resembles the partially disordered phase proposed by Mekata is observed in the finite temperature region. But it is shown in this region that the sublattice magnetizations vanish at the thermodynamic limit and the Kosterlitz-Thouless phase transition is predicted by the finite size analysis.
Frustrated spin systems exhibit a wide range of phenomena such as novel phases, residual entropy and glassy dynamics. In contrast to the usual assumption of universality, microscopic details such as the underlying lattice structure often play a crucial role in these systems. In this dissertation the frustrated J1-J2 Ising model on the honeycomb lattice is considered. Earlier work shows that the coupling ratio R = J2 /|J1 | = −1/4 separates regions of two different ground states. Above this threshold the ground state is (anti-)ferromagnetically ordered for (anti-)ferromagnetic J1 , whereas below it infinitely many states of minimal energy exist. A recent Monte Carlo (MC) study by Žukovič finds a continuous transition with Ising critical exponents for R ≥ −0.2, leaving open the question of criticality for −1/4 < R < −0.2, for which the study was obfuscated by the presence of strong metastability. The regime with the highly degenerate ground states (R < −1/4), on the other hand, is characterized by feature-rich specific heat curves, and slow low-temperature dynamics reminiscent of spin glasses, although there is no disorder in this system. Prior to this work, no order parameter has been identified, and there is no clear understanding of its potentially critical behavior. Here, the question of criticality for −1/4 < R < −0.2 is addressed, and Ising universality is found for lower values of R, namely, at least for R ≥ −0.23. This result is achieved by carrying out finite-size scaling analyses of response functions like the magnetic susceptibility, and of the complex partition function zeros. In both cases, the underlying numerical data was generated using population annealing (PA) MC, a framework devised for spin glasses but generally applicable to systems with metastability. Furthermore, the region of degenerate ground states is studied using different boundary conditions, which are found to strongly affect the thermodynamic observations. A broken symmetry at zero temperature is identified, which can be quantified using a complex nematic order parameter. The MC data are consistent with a continuous symmetry-breaking transition. However, due to the limited system sizes that could be studied, and the presence of strong finite-size effects, one cannot rule out a higher-order transition or a cross-over without any singularities. In the process of studying this system, also the understanding of the PA framework was further developed. Specifically, the role of the resampling step in the algorithm was investigated. While it is crucial for efficient equilibration, no benefit is found in noisier resampling schemes such as the multinomial resampling method. Instead, two resampling schemes, which have the lowest sampling variance, are identified as preferential.:1 Introduction I Background and Methodology 2 Background 2.1 Fundamentals of Statistical Mechanics 2.1.1 Equilibrium ensembles 2.1.2 Phase transitions and critical behavior 2.1.3 Common spin models 2.2 Frustrated Spin Models 2.2.1 Frustrated J1-J2 Ising model on the square lattice 2.2.2 Gonihedric model in two dimensions 2.3 Monte Carlo Simulations 2.3.1 Fundamentals 2.3.2 Autocorrelation and error bars in MCMC time series 2.3.3 Possible implementations 2.3.4 Generalized ensemble methods 3 Methodology 3.1 General Algorithmic Design 3.2 Population Annealing 3.2.1 Algorithm 3.2.2 Assessing population annealing performance 3.2.3 Adaptive temperature steps 3.2.4 Weighted averages 3.2.5 Population annealing and first-order transitions 3.3 Estimating Critical Exponents from Finite Systems 3.3.1 Standard finite size scaling 3.3.2 Complex partition function zeros II Results 4 Algorithmic Improvements and Extensions to the Population Annealing Framework 4.1 Resampling Step in Population Annealing 4.1.1 Conceptual overview 4.1.2 Possible resampling methods 4.1.3 Benchmarking different resampling methods 4.1.4 Simulation results 4.1.5 PA simulations with finite populations in the limit of perfect equilibration 4.1.6 PA simulations in the limit of infinite populations and perfect equilibration 4.1.7 Conclusion 4.2 Using Population Annealing on Many Computing Nodes 4.2.1 Literature overview and different rebalancing schemes 4.2.2 Simulations of simulations 4.2.3 Conclusion and outlook 4.3 Adaptive Number of Spin Updates in Population Annealing 4.4 Microcanonical Population Annealing 5 The Frustrated J1-J2 Ising Model on the Honeycomb Lattice 5.1 Overview 5.2 Regime of Ising Universality Class 5.2.1 Finite-size scaling analysis 5.2.2 Energy barriers 5.2.3 Exploring the complex partition function zeros 5.3 Finite Ground State Entropy per Spin for J2/|J1| = −1/4 5.4 Highly Degenerate Striped Phase for J2/|J1| < −1/4 5.4.1 Boundary conditions 5.4.2 Energy minimization and ground states 5.4.3 Topological defects 5.4.4 Observables 5.4.5 Simulation details 5.4.6 Simulation results for R = −0.5 5.4.7 Simulation results for R = −1 5.4.8 Discussion 5.4.9 Conclusion and outlook 6 Parallel tempering and its infinite-swap limit 6.1 Formalism 6.2 Details on Implementation 6.3 Numerical Results from INS-PT Simulations 6.3.1 One-dimensional Ising model 6.3.2 Two-dimensional Ising model 6.3.3 q-state Potts model 6.3.4 Lennard-Jones clusters 6.4 Discussion 7 Conclusions and outlook Bibliography Appendices A Algorithmic Improvements and Extensions to the Population Annealing framework: Further details A.1 Generating correlated energy time series B J1-J2 honeycomb Ising model: Further plots and details B.1 GPU domain decomposition B.2 Finite-size scaling using Metropolis population annealing B.3 The decoupled antiferromagnet limit B.4 Zero temperature value of the nematic order parameter B.5 Percolation properties for J2 = −0.5 B.6 Binder parameter of the one-dimensional Ising model C Application for Ethical Approval
In this work we study a disordered binary Ising model on the square lattice. The model system consists of two different particles with spin-1/2 and spin-1, which are randomly distributed on the lattice. It has been considered only spin nearest-neighbor exchange interactions with J>0. This system can represent a disordered magnetic binary alloy AxB1−x, obtained from the high temperature quenching of a liquid mixture. The results were obtained by the use of Monte Carlo simulations for several lattice sizes L, temperature T and concentration x of ions A with spin-1/2. We found its critical temperature, through the reduced fourth-order Binder cumulant for the several values of the concentration x of particles (spin-1/2, spin-1), and also the magnetization, the susceptibility and the specific heat as a function of temperature T.