Disorder is an inevitable part of any condensed matter system and therefore its study has always been of great importance. The effect of quenched randomness on a system that exhibits a continuous phase transition in the absence of any impurity has been studied in the past and the results are relatively well understood. However, the effect of quenched randomness on \\emph{first-order} phase transitions is still not well understood. In this dissertation, we study the effect of quenched bond-randomness on the classical and quantum first-order phase transitions.In Chapter 2, we study the effect of the disordered three-color Ashkin-Teller model, whose pure version undergoes a first-order phase transition. We show that the rounding of the first-order transition of the pure model due to the impurities is manifested as a critical point. We conclusively rule out that the model belongs to the universality class of the two-dimensional Ising model. Furthermore, we find that the exponents
beta and
nu vary with disorder and the four-spin coupling strength.In Chapter 3, we extend our study of the disordered three-color Ashkin-Teller model. Utilizing extensive cluster Monte Carlo simulations on large lattice sizes of up to
128times128 spins, each of which is represented by three colors taking values
pm1, we show that the rounding of the first-order phase transition is an emergent criticality.We find that the critical exponents,
nu and
beta, change as the strength of disorder or the four-spin coupling varies, and we show that the correlation length critical exponent violates the lower bound
2/Dlenu, where
D is the dimension of the system.In Chapter 4, we study the \\emph{quantum} three-color Ashkin--Teller model and show that the quantum critical point in
(1+1) dimension is an unusual one, with activated scaling at the critical point and Griffiths-McCoy phase away from it. We find that the behavior is similar to the transverse random field Ising model, even though the pure system has a first-order transition in this case.
A cluster Monte Carlo algorithm for the Ashkin-Teller (AT) model is constructed according to the guidelines of a general scheme for such algorithms. Its dynamical behavior is tested for the square-lattice AT model. We perform simulations on the line of critical points along which the exponents vary continuously, and find that critical slowing down is significantly reduced. We find continuous variation of the dynamical exponent z along the line, following the variation of the ratio of specific-heat and correlation-length exponents \ensuremath{\alpha}/\ensuremath{\nu}, in a manner which satisfies the Li-Sokal bound
zcluster\ensuremath{\ge}\ensuremath{\alpha}/\ensuremath{\nu}, that was so far proved only for Potts models.
An intriguing result of statistical mechanics is that a first-order phase transition can be rounded by disorder coupled to energylike variables. In fact, even more intriguing is that the rounding may manifest itself as a critical point, quantum or classical. In general, it is not known, however, what universality classes, if any, such criticalities belong to. In order to shed light on this question we examine in detail the disordered three-color Ashkin-Teller model by Monte Carlo methods. Extensive analyses indicate that the critical exponents define a new universality class. We show that the rounding of the first-order transition of the pure model due to the impurities is manifested as criticality. However, the magnetization critical exponent, β, and the correlation length critical exponent, ν, are found to vary with disorder and the four-spin coupling strength, and we conclusively rule out that the model belongs to the universality class of the two-dimensional Ising model.
We consider the sample to sample fluctuations that occur in the value of a thermodynamic quantity
P in an ensemble of finite systems with quenched disorder, at equilibrium. The variance of
P,
VP, which characterizes these fluctuations is calculated as a function of the systems' linear size
l, focusing on the behavior at the critical point. The specific model considered is the bond-disordered Ashkin-Teller model on a square lattice [Phys. Rev. 64, 178 (1943)]. Using extensive Monte Carlo simulations, several bond-disordered Ashkin-Teller models were examined, including the bond-disordered Ising model and the bond-disordered four-state Potts model. It was found that far from criticality all thermodynamic quantities which were examined (energy, magnetization, specific heat, susceptibility) are strongly self-averaging, that is ${V}_{P}\ensuremath{\sim}{l}^{\ensuremath{-}d}$ (where
d=2 is the dimension). At criticality though, the results indicate that the magnetization
M and the susceptibility $\ensuremath{\chi}$ are non-self-averaging, i.e., $\frac{{V}_{\ensuremath{\chi}}}{{\ensuremath{\chi}}^{2}}, \frac{{V}_{M}}{{M}^{2}}\ensuremath{\nrightarrow}0$. The energy
E at criticality is clearly weakly self-averaging, that is ${V}_{E}\ensuremath{\sim}{l}^{\ensuremath{-}{y}_{v}}$ with
0<yv<d. Less conclusively, and possibly only as a transient behavior, the specific heat too is found to be weakly self-averaging. A phenomenological theory of finite size scaling for disordered systems is developed, based on physical considerations similar to those leading to the Harris criterion. Its main prediction is that when the specific heat exponent $\ensuremath{\alpha}<0$ ($\ensuremath{\alpha}$ of the disordered model) then, for a quantity
P which scales as ${l}^{\ensuremath{\rho}}$ at criticality, its variance
VP will scale asymptotically as ${l}^{2\ensuremath{\rho}+\frac{\ensuremath{\alpha}}{\ensuremath{\nu}}}$. The theory is not applicable in the asymptotic limit ($l\ensuremath{\rightarrow}\ensuremath{\infty}$) to the bond-disordered Ashkin-Teller model where $\frac{\ensuremath{\alpha}}{\ensuremath{\nu}}={0}_{+}$. Nonetheless in the accessible range of lattice sizes we found very good agreement between the theory and the data for ${V}_{\ensuremath{\chi}}$ and
VE. The theory may also be compatible with the data for the variance of the magnetization
VM and the variance of the specific heat
VC, but evidence for this is less convincing.
We evaluate by Monte Carlo simulations various singular thermodynamic quantities
X for ensembles of quenched random Ising and Ashkin-Teller models. The measurements are taken at
Tc and we study how the distributions
P(X) (and, in particular, their relative squared width,
RX) over the ensemble depend on the system size
l. The Ashkin-Teller model was studied in the regime where bond randomness is irrelevant and we found weak self-averaging; ${R}_{X}\ensuremath{\sim}{l}^{\ensuremath{\alpha}/\ensuremath{\nu}}\ensuremath{\rightarrow}0,$ where $\ensuremath{\alpha}<0$ and $\ensuremath{\nu}$ are the exponents (of the pure model fixed point) governing the transition. For the site-dilute Ising model on a cubic lattice, known to be governed by a random fixed point, we find that
RX tends to a constant, as predicted by Aharony and Harris. We tested whether this constant is universal. However, this constant is different for canonical and grand canonical disorder. We identify the pseudocritical temperature of each sample
i, Tc(i,l), as the temperature at which the susceptibility reaches its maximal value ${\ensuremath{\chi}}^{\mathrm{max}}.$ The distribution of these sample dependent
Tc(i,l) was investigated; we found that its variance scales as $[\ensuremath{\delta}{T}_{c}(l){]}^{2}\ensuremath{\sim}{l}^{\ensuremath{-}2/\ensuremath{\nu}}.$ Our previously proposed finite size scaling ansatz for disordered systems was tested and found to hold. We did observe deviations from a single function, which imply that sample dependent scaling functions are needed. These deviations are, however, relatively small and hence to obtain a fixed statistical error it may be more computationally efficient to measure ${\ensuremath{\chi}}^{\mathrm{max}}$ than the commonly used $\ensuremath{\chi}{(T}_{c}^{\ensuremath{\infty}}).$
The results of the large-scale Monte Carlo simulations of the continuous phase transitions in the three-dimensional Ashkin-Teller model on a cubic lattice are presented in the two-parameter space. The invariance of the ratio of the square of the second moment of the order parameter to its fourth moment in the critical region and the Metropolis algorithm are exploited. Using the universality hypothesis and the finite-size-scaling analysis, the accurate positions of the continuous phase-transition lines are calculated and their Ising character is confirmed. Approximate locations of the tricritical points are also estimated.