We study the Néel-paramagnetic quantum phase transition in two-dimensional dimerized S=1/2 Heisenberg antiferromagnets using finite-size scaling of quantum Monte Carlo data. We resolve the long-standing issue of the role of cubic interactions arising in the bond-operator representation when the dimer pattern lacks a certain symmetry. We find nonmonotonic (monotonic) size dependence in the staggered (columnar) dimerized model, where cubic interactions are (are not) present. We conclude that there is a new irrelevant field in the staggered model, but, at variance with previous claims, it is not the leading irrelevant field. The new exponent is ω_{2}≈1.25 and the prefactor of the correction L^{-ω_{2}} is large and comes with a different sign from that of the conventional correction with ω_{1}≈0.78. Our study highlights competing scaling corrections at quantum critical points.
Scaling, universality, and renormalization are three pillars of modern critical phenomena. According to the hypothesis of universality, continuous phase transitions fall into classes mainly determined by spatial dimensionality and symmetry of the order parameter. The latter is usually reflected by the degeneracy of the ground state of the Hamiltonian. However, for certain systems at criticality, a higher symmetry emerges in the order parameter, and the associated critical behavior may become very rich. Such emergent symmetry has been found in spin ice systems, deconfined quantum critical points, high-
Tc superconductors, and so forth, and are often accompanied by very interesting critical phenomena. New results presented here are based on Monte Carlo simulations and finite-size scaling of a series of
q-state Potts models on the simple cubic lattice with ferromagnetic interactions in one lattice direction and antiferromagnetic interactions in the other two directions. The staggered magnetization appears to display an emergent continuous O(
n) symmetry with
n=
q-1, as illustrated by two-dimensional intersections of the distribution functions. Also the estimated critical exponents are consistent with the O(
n=
q-1) universality classes.
We investigate the quantum phase transition in an
S=1/2 dimerized Heisenberg antiferromagnet in three spatial dimensions. By performing large-scale quantum Monte Carlo simulations and detailed finite-size scaling analyses, we obtain high-precision results for the quantum critical properties at the transition from the magnetically disordered dimer-singlet phase to the antiferromagnetically ordered N\'eel phase. This transition breaks
O(N) symmetry with
N=3 in
D=3+1 dimensions. This is the upper critical dimension, where multiplicative logarithmic corrections to the leading mean-field critical properties are expected; we extract these corrections, establishing their precise forms for both the zero-temperature staggered magnetization
ms and the N\'eel temperature
TN. We present a scaling ansatz for
TN, including logarithmic corrections, which agrees with our data and indicates exact linearity with
ms, implying a complete decoupling of quantum and thermal fluctuation effects even arbitrarily close to the quantum critical point. We also demonstrate the predicted
N-independent leading and subleading logarithmic corrections in the size dependence of the staggered magnetic susceptibility. These logarithmic scaling forms have not previously been identified or verified by unbiased numerical methods, and we discuss their relevance to experimental studies of dimerized quantum antiferromagnets such as
TlCuCl3.
In the study of strongly interacting condensed-matter systems controlled microscopic theories hold a key position. Spin-wave theory, large-N expansion, and
epsilon-expansion are some of the few successful cornerstones. In this doctoral thesis work, we have developed a novel large-
d expansion method,
d being the spatial dimension, to study model Hamiltonians hosting a quantum phase transition between a paramagnet and a magnetically ordered phase. A highlight of this technique is that it can consistently describe the entire phase diagram of the above mentioned models, including the quantum critical point. Note that most analytical techniques either efficiently describe only one of the phases or suffer from divergences near the critical point. The idea of large-
d formalism is that in this limit, non-local fluctuations become unimportant and that a suitable product state delivers exact expectation values for local observables, with corrections being suppressed in powers of
1/d. It turns out that, due to momentum summation properties of the interaction structure factor, all diagrams are suppressed in powers of
1/d leading to an analytic expansion. We have demonstrated this method in two important systems namely, the coupled-dimer magnets and the transverse-field Ising model. Coupled-dimer magnets are Heisenberg spin systems with two spins, coupled by intra-dimer antiferromagnetic interaction, per crystallographic unit cell (dimer). In turn, spins from neighboring dimers interact via some inter-dimer interaction. A quantum paramagnet is realized for a dominant intra-dimer interaction, while a magnetically ordered phase exists for a dominant (or of the same order as intra-dimer interaction) inter-dimer interaction. These two phases are connected by a quantum phase transition, which is in the Heisenberg O(3) universality class. Microscopic analytical theories to study such systems have been restricted to either only one of the phases or involve uncontrolled approximations. Using a non-linear bond-operator theory for spins with S=
1/2, we have calculated the
1/d expansion of static and dynamic observables for coupled dimers on a hypercubic lattice at zero temperature. Analyticity of the
1/d expansion, even at the critical point, is ensured by correctly identifying suitable observables using the mean-field critical exponents. This method yields gapless excitation modes in the continuous symmetry broken phase, as required by Goldstone\'s theorem. In appropriate limits, our results match with perturbation expansion in small ratio of inter-dimer and intra-dimer coupling, performed using continuous unitary transformations, as well as the spin-wave theory for spin-
1/2 in arbitrary dimensions. We also discuss the Brueckner approach, which relies on small quasiparticle density, and derive the same
1/d expansion for the dispersion relation in the disordered phase. Another success of our work is in describing the amplitude (Higgs) mode in coupled-dimer magnets. Our novel method establishes the popular bond-operator theory as a controlled approach. In
d=2, the results from our calculations are in qualitative agreement with the quantum Monte Carlo study of the square-lattice bilayer Heisenberg AF spin-
1/2 model. In particular, our results are useful to identify the amplitude (Higgs) mode in the QMC data. The ideas of large-
d are also successfully applied to the transverse-field Ising model on a hypercubic lattice. Similar to bond operators, we have introduced auxiliary Bosonsic operators to set up our method in this case. We have also discussed briefly the bilayer Kitaev model, constructed by antiferromagnetically coupling two layers of the Kitaev model on a honeycomb lattice. In this case, we investigate the dimer quantum paramagnetic phase, realized in the strong inter-layer coupling limit. Using bond-operator theory, we calculate the mode dispersion in this phase, within the harmonic approximation. We also conjecture a zero-temperature phase diagram for this model.
We study the critical properties in cubic systems of antiferromagnetically coupled spin dimers near magnetic-field-induced quantum phase transitions. The quantum critical points in the zero-temperature phase diagram are determined from quantum Monte Carlo simulations. Furthermore, scaling properties of the uniform magnetization and the staggered transverse magnetization across the quantum phase transition in magnetic fields are calculated. Excellent agreement is observed between the quantum Monte Carlo simulations and previous analytical and experimental results near the magnetic field-induced transition. At the zero-field quantum critical point, logarithmic corrections to mean-field behavior are observed, which match well recent experiments for the pressure-induced quantum phase transition.