Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMonte Carlo simulations for Ising spins with spin greater than 1/2 applied to the square and triangular lattices with antiferromagnetic interactions and comparing results using Kawasaki and Glauber dynamics
439
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
This paper has a pedagogical introduction. We describe the correct method for performing Monte Carlo simulations of Ising model systems with spin greater than one half. Correct and incorrect procedures are clearly outlined and the consequences of using the incorrect procedure are shown. The difference between Kawasaki and Glauber dynamics is then outlined and both methods are applied to the antiferromagnetic square and triangular lattices for S =1.
rate research
Read More
We study the thermodynamics of Ising spins on the triangular kagome lattice (TKL) using exact analytic methods as well as Monte Carlo simulations. We present the free energy, internal energy, specific heat, entropy, sublattice magnetizations, and susceptibility. We describe the rich phase diagram of the model as a function of coupling constants, temperature, and applied magnetic field. For frustrated interactions in the absence of applied field, the ground state is a spin liquid phase with integer residual entropy per spin $s_0/k_B={1/9} ln 72approx 0.4752...$. In weak applied field, the system maps to the dimer model on a honeycomb lattice, with irrational residual entropy 0.0359 per spin and quasi-long-range order with power-law spin-spin correlations that should be detectable by neutron scattering. The power-law correlations become exponential at finite temperatures, but the correlation length may still be long.
The Ising model S=1/2 and the S=1 model are studied by efficient Monte Carlo schemes on the (3,4,6,4) and the (3,3,3,3,6) Archimedean lattices. The algorithms used, a hybrid Metropolis-Wolff algorithm and a parallel tempering protocol, are briefly described and compared with the simple Metropolis algorithm. Accurate Monte Carlo data are produced at the exact critical temperatures of the Ising model for these lattices. Their finite-size analysis provide, with high accuracy, all critical exponents which, as expected, are the same with the well known 2d Ising model exact values. A detailed finite-size scaling analysis of our Monte Carlo data for the S=1 model on the same lattices provides very clear evidence that this model obeys, also very well, the 2d Ising model critical exponents. As a result, we find that recent Monte Carlo simulations and attempts to define effective dimensionality for the S=1 model on these lattices are misleading. Accurate estimates are obtained for the critical amplitudes of the logarithmic expansions of the specific heat for both models on the two Archimedean lattices.
The ground state and thermodynamics of a generalized spin-1/2 Ising-Heisenberg diamond chain with the second-neighbor interaction between nodal spins are calculated exactly using the mapping method based on the decoration-iteration transformation. Rigorous results for the magnetization, susceptibility, and heat capacity are investigated in dependence on temperature and magnetic field for the frustrated diamond spin chain with the antiferromagnetic Ising and Heisenberg interactions. It is demonstrated that the second-neighbor interaction between nodal spins gives rise to a greater diversity of low-temperature magnetization curves, which may include an intermediate plateau at two-third of the saturation magnetization related to the classical ferrimagnetic (up-up-up-down-up-up-...) ground state with translationally broken symmetry besides an intermediate one-third magnetization plateau reflecting the translationally invariant quantum ferrimagnetic (monomer-dimer) spin arrangement.
We implement a new and accurate numerical entropic scheme to investigate the first-order transition features of the triangular Ising model with nearest-neighbor ($J_{nn}$) and next-nearest-neighbor ($J_{nnn}$) antiferromagnetic interactions in ratio $R=J_{nn}/J_{nnn}=1$. Important aspects of the existing theories of first-order transitions are briefly reviewed, tested on this model, and compared with previous work on the Potts model. Using lattices with linear sizes $L=30,40,...,100,120,140,160,200,240,360$ and 480 we estimate the thermal characteristics of the present weak first-order transition. Our results improve the original estimates of Rastelli et al. and verify all the generally accepted predictions of the finite-size scaling theory of first-order transitions, including transition point shifts, thermal, and magnetic anomalies. However, two of our findings are not compatible with current phenomenological expectations. The behavior of transition points, derived from the number-of-phases parameter, is not in accordance with the theoretically conjectured exponentially small shift behavior and the well-known double Gaussian approximation does not correctly describe higher correction terms of the energy cumulants. It is argued that this discrepancy has its origin in the commonly neglected contributions from domain wall corrections.
The recently fabricated two-dimensional magnetic materials Cu9X2(cpa)6.xH2O (cpa=2-carboxypentonic acid; X=F,Cl,Br) have copper sites which form a triangular kagome lattice (TKL), formed by introducing small triangles (``a-trimers) inside of each kagome triangle (``b-trimer). We show that in the limit where spins residing on b-trimers have Ising character, quantum fluctuations of XXZ spins residing on the a-trimers can be exactly accounted for in the absence of applied field. This is accomplished through a mapping to the kagome Ising model, for which exact analytic solutions exist. We derive the complete finite temperature phase diagram for this XXZ-Ising model, including the residual zero temperature entropies of the seven ground state phases. Whereas the disordered (spin liquid) ground state of the pure Ising TKL model has macroscopic residual entropy ln72=4.2767... per unit cell, the introduction of transverse(quantum) couplings between neighboring $a$-spins reduces this entropy to 2.5258... per unit cell. In the presence of applied magnetic field, we map the TKL XXZ-Ising model to the kagome Ising model with three-spin interactions, and derive the ground state phase diagram. A small (or even infinitesimal) field leads to a new phase that corresponds to a non-intersecting loop gas on the kagome lattice, with entropy 1.4053... per unit cell and a mean magnetization for the b-spins of 0.12(1) per site. In addition, we find that for moderate applied field, there is a critical spin liquid phase which maps to close-packed dimers on the honeycomb lattice, which survives even when the a-spins are in the Heisenberg limit.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
