We study the cubic fixed point for
N=3 and 4 by using finite-size scaling applied to data obtained from Monte Carlo simulations of the
N-component ${\ensuremath{\phi}}^{4}$ model on the simple cubic lattice. We generalize the idea of improved models to a two-parameter family of models. The two-parameter space is scanned for the point, where the amplitudes of the two leading corrections to scaling vanish. To this end, a dimensionless quantity is introduced that monitors the breaking of the
O(N) invariance. For
N=4, we determine the correction exponents ${\ensuremath{\omega}}_{1}=0.763(24)$ and ${\ensuremath{\omega}}_{2}=0.082(5)$. In the case of
N=3, we obtain
Y4=0.0142(6) for the renormalization group exponent of the cubic perturbation at the
O(3)-invariant fixed point, while the correction exponent ${\ensuremath{\omega}}_{2}=0.0133(8)$ at the cubic fixed point. Simulations close to the improved point result in the estimates $\ensuremath{\nu}=0.7202(7)$ and $\ensuremath{\eta}=0.0371(2)$ of the critical exponents of the cubic fixed point for
N=4. For
N=3, at the cubic fixed point, the
O(3) symmetry is only mildly broken and the critical exponents differ only by little from those of the
O(3)-invariant fixed point. We find $\ensuremath{-}0.000\phantom{\rule{0.16em}{0ex}}01⪅{\ensuremath{\eta}}_{\text{cubic}}\ensuremath{-}{\ensuremath{\eta}}_{O(3)}⪅0.000\phantom{\rule{0.16em}{0ex}}07$ and ${\ensuremath{\nu}}_{\text{cubic}}\ensuremath{-}{\ensuremath{\nu}}_{O(3)}=\ensuremath{-}0.000\phantom{\rule{0.16em}{0ex}}61(10)$.
We use extensive Monte Carlo transfer-matrix calculations on infinite strips of widths L up to 30 lattice spacing and a finite-size scaling analysis to obtain critical exponents and conformal anomaly number c for the two-dimensional XY Ising model. This model is expected to describe the critical behavior of a class of systems with simultaneous U(1) and
Z2 symmetries of which the fully frustrated XY model is a special case. The effective values obtained for c show a significant decrease with L at different points along the line where the transition to the ordered phase takes place in a single transition. Extrapolations based on power-law corrections give values consistent with c=3/2 although larger values cannot be ruled out. Critical exponents are obtained more accurately and are consistent with previous Monte Carlo simulations suggesting critical behavior and with recent calculations for the frustrated XY model.
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.
We have performed large-scale numerical simulations for the ± J Ising model on the simple cubic and the square lattices with asymmetric weight of ferromagnetic and antiferromagnetic bonds. The Monte Carlo renormalization group method is used to estimate T c , γ/ν, β/ν and ν at various points along the phase boundary between ferromagnetic and paramagnetic phases including the tricritical point. The obtained values of T c and the resulting phase diagrams of the simple cubic and the square lattices are much more accurate than those estimated so far. In the simple cubic lattice, we find evidence for the possibility that universality of the ferromagnetic critical exponents holds along the critical line except at the tricritical point, while weak universality holds even at the tricritical point. Weak universality seems to hold also in the square lattice.
We study the three-component ${\ensuremath{\phi}}^{4}$ model on the simple cubic lattice in the presence of a cubic perturbation. To this end, we perform Monte Carlo simulations in conjunction with a finite-size scaling analysis of the data. The analysis of the renormalization group (RG) flow of a dimensionless quantity provides us with the accurate estimate ${Y}_{4}\ensuremath{-}{\ensuremath{\omega}}_{2}=0.000\phantom{\rule{0.16em}{0ex}}81(7)$ for the difference of the RG eigenvalue
Y4 at the
O(3)-symmetric fixed point and the correction exponent ${\ensuremath{\omega}}_{2}$ at the cubic fixed point. We determine an effective exponent ${\ensuremath{\nu}}_{\text{eff}}$ of the correlation length that depends on the strength of the breaking of the
O(3) symmetry. Field theory predicts that depending on the sign of the cubic perturbation, the RG flow is attracted by the cubic fixed point, or runs to an ever increasing amplitude, indicating a fluctuation-induced first-order phase transition. We demonstrate directly the first-order nature of the phase transition for a sufficiently strong breaking of the O(3) symmetry. We obtain accurate results for the latent heat, the correlation length in the disordered phase at the transition temperature, and the interface tension for interfaces between one of the ordered phases and the disordered phase. We study how these quantities scale with the RG flow, allowing quantitative predictions for weaker breaking of the
O(3) symmetry.