In this master thesis the criticality of the 3D XY model with cubic anisotropy has been examined by the means of Monte Carlo simulation. The 3D XY model is a special case of the
O(n) spin-vector model which describe numerous phenomena's. Here,
n is the number of vector dimensions and for the XY model
n=2. The
O(n)-model is invariant for global rotations. This means that it is arbitrary which direction in the symmetry breaking occurs. In real system one should always expect perturbations from a perfect system, here it is done by adding a cubic anisotropy (order
p=4) then the symmetry of the Hamiltonian is lowered and it will be advantageous for the spins to be ordered in directions
90o apart. This symmetry is called
C4, which is a part of the cyclic group, and it restricts the symmetry breaking to one of these four directions. \n\nThe criticality of this model is when the spins suddenly align as the temperature is lowered to a critical temperature. Sufficiently close to the critical temperature quantities like the magnetization, which is the sum of spins, will show a power law growth characterised by critical exponents specified by the class of universality of the transition. The critical exponents can be obtained by finite-size scaling. By applying renormalization theory on the critical region one can show that some of parameters in the Hamiltonian are: relevant and effect the scaling, irrelevant and does not effect the scaling and dangerously irrelevant meaning that they only affect some of the critical exponents. Previous Monte Carlo and renormalization calculus show that anisotropy which is relevant for
p<4 becomes dangerously irrelevant for
pgeq4. The question of interest is what kind of effect the dangerously irrelevant parameter in the Hamiltonian has on the scaling. This is interesting because the effect of dangerously irrelevant parameters is hardly studied and anisotropy arise in many crystal systems because of the lattice geometry. \n\nThe results are obtained by Monte Carlo simulations and the critical slowing down close to the phase transition is solved by the method of parallel tempering. This method was chosen because it is simple to implement in a parallel code.\n\nThe Binder cumulant for different quantities and the peak of the magnetic susceptibility was used support a critical temperature
Tc=2.2019textnormalk for no field applied and no anisotropy. The known critical exponents for the no field and no anisotropy was used. The scaling made a good collapse, meaning that the right
Tc and the correct critical exponents where used. The scaling did not take into account a magnetic field and therefore a good collapse of curves where not obtained in this case. When the cubic anisotropy was added the critical exponent for the magnetic susceptibility and the cubic order parameter changed, while the other exponents where unchanged. This supports the view that cubic anisotropy is dangerously irrelevant for the 3D XY model. The scaling of the magnetic susceptibility showed a possible new kind of scaling, a scaling which is different on the two sides of the transition. Although a more accurate determination of the critical exponent in the ordered phase require more extensive simulations. It appears that this exponent is larger than the exponent in the disordered phase, where one recovered the isotropic exponent. The difference in the critical exponents has been predicted to increase with
p and an investigation of the hexagonal anisotropy
p=6 has been started.
The anisotropic Heisenberg model has been investigated by interpolating anisotropy between the Ising and isotropic Heisenberg limits via general spin high temperature series (HTS) for the zero-field susceptibility of the ferromagnets for the square lattice. This is done by locating the critical temperature ( T c (2) ) and the susceptibility exponent (γ) from the analysis of the series by Pade approximant and ratio methods. The critical temperature is found to be spin dependent. It decreases as anisotropy changes from the Ising limit towards the isotropic limit and remains finite (lower than that of Stanley and Kaplan) for all the spins except for S =1/2 at the isotropic limit. Exponent γ is almost Ising like for the anisotropic system and anomalous jump in it occurs for the isotropic case. In addition, expressions for the magnetization are derived by using spin-wave theory. The limit of complete isotropy has been found not to be a special case for the susceptibility, in contrast to magnetization.
The renormalization group flows of the one-dimensional anisotropic XY model are obtained by different multiscale entanglement renormalization ansatz schemes. It is shown that the optimized disentangler removes the short-range entanglement by rotating the system in the parameter space spanned by the anisotropy and the magnetic field. It is understood from the study that the disentangler reduces the entanglement by mapping the system to another one in the same universality class but with smaller short-range entanglement. The phase boundary and corresponding critical exponents are calculated using different schemes with different block sizes, look-ahead steps and truncation dimensions. It is shown that a larger truncation dimension leads to more accurate results and that using a larger block size or look-ahead step improves the overall calculational consistency.
We study the entanglement near a quantum phase transition in a class of exactly solvable one‐dimensional quantum magnets, the anisotropic XY models. The critical behaviour (scaling) is reflected by the derivatives of the bipartite entanglent, quantified by the concurrence. However, we found a profound difference between classical correlation and non‐local quantum correlation: The range of the concurrence is found to behave in a non‐universal way, being extremely short‐ranged for the Ising model. This range tends to infinity towards the isotropic XY model.
We have calculated the spin-correlation functions (in the x, y, and z directions) in the linear spin-(1/2 XY system for both isotropic and anisotropic couplings. The study of zero-temperature correlation functions shows that a long-range order develops in the direction in which the coupling is stronger, but no such ordering exists for isotropic systems. We have studied the temperature variation of inverse correlation lengths in these directions and it has been found that at T=0 all the correlation lengths (${\ensuremath{\xi}}_{x}$, ${\ensuremath{\xi}}_{y}$, ${\ensuremath{\xi}}_{z}$) diverge in the isotropic case. On the other hand, in the anisotropic case, the correlation length diverges in the direction (either x or y) in which the coupling becomes stronger. The results are compared with the experimental data on
Cs2CoCl4. The critical exponents (\ensuremath{\eta}',\ensuremath{\nu}') of the correlation function and correlation length are also calculated near the critical temperature (T=0) for different anisotropies, and it is found that the system behaves like an Ising model when a little anisotropy is introduced.