Arxiv:2005.01990V2 [Cond-Mat.Str-El] 14 Aug 2020
Total Page:16
File Type:pdf, Size:1020Kb
2 where hi; ji is a pair of the nearest neibghbor sites and the III. AWAY FROM QUANTUM CRITICAL iπyi energy unit is t = 1. The hopping is tij = te exp(iAij) POINT along the x-direction on the y = yi bond and tij = t exp(iAij) along the y-direction. The vector potential is Before discussing quantum criticality, we investigate 57 0 given in the string gauge with the period Lx ×Ly where the magnetic catalysis when the system is away from the 0 Lx is the superlattice unit period used in the iDMRG cal- critical point. Firstly, we consider a weak interaction culations for the system size Lx×Ly = 1×Ly. Typically, V = 0:50t < Vc = 1:30t for which the system at B = 0 0 we use Lx = 20 for Ly = 6 and Lx = 10 for Ly = 10. 0 is a Dirac semimetal renormalized by the interaction. See also Appendix A. Aij = 0 corresponds to the con- As exemplified in Fig. 1, dependence of M(B) on the ventional π-flux square lattice without an applied field, bond dimension χ used in the calculation is negligibly while Aij 6= 0 describes an applied magnetic field for a small for L = 6, and it can be safely extrapolated to P y plaquette p, Bp = hiji2p Aij. The magnetic field is χ ! 1 even for Ly = 10. Standard deviations of the spatially uniform and an integer multiple of a unit value extrapolations are less than 1% and within the symbols. 0 B = n × δB (n = 1; 2; ··· ;LxLy) allowed by the super- Such an extrapolation can be done also for other values 0 lattice size, where δB = 2π=LxLy. The lattice constant of V as mentioned before, and all results shown in this a as a length unit and the electric charge e have been set study are extrapolated ones. as a = 1; e = 1, and the magnetic field is measured in the unit of B0 = 2π. 0.1 The V -term is a repulsive nearest neighbor interaction 0.08 leading to the CDW state and the quantum phase transi- 0.06 tion with B = 0 takes place at V = Vc ' 1:30t according M 0.04 to the quantum Monte Carlo calculations for the bulk two dimensional system, where the criticality belongs to 0.02 the (2 + 1)-dimensional chiral Ising universality class3{8. 0 On the other hand, our cylinder system is anisotropic and 0 0.5 1 1.5 2 the CDW quantum phase transition at B = 0 is simply 1/χ (x 10-3) (1+1)-dimensional Ising transition and critical interac- tion strength depends on the system size Ly, which may 58 FIG. 1. Extrapolation of the CDW order parameter M(B) be regarded as a finite size effect . However, the system for the χ ! 1 limit at V = 0:50t. The blue (red) symbols are can be essentially two-dimensional in space underp a mag- for Ly = 6(Ly = 10) and the curves are power law fittings. netic field when the magnetic length lB = 1= B becomes Each curve corresponds to a magnetic field in the range 0 ≤ shorter than the system size Ly. We will use this prop- B ≤ 0:06B0. erty to discuss the (2 + 1)-dimensional criticality. Note that the critical interaction strength Vc ' 1:30t will be In Fig. 2 (a), we show the CDW order parameter M confirmed later within the present framework. extrapolated to χ ! 1 for the system sizes Ly = 6 and In the following, we focus on the CDW order parame- Ly = 10 at V = 0:50t. The calculated results almost con- ter, verge for Ly = 6; 10 and are independent of the system size, except for B = 0 where there is a finite size effect due to l = 1, although there is some accidental devia- B 1 X tion around B ' 0:1B0. Therefore, these results give the M = (−1)jijn ; (2) L0 L i CDW order parameter essentially in the thermodynamic x y i limit Ly ! 1. In order to understand impacts of quan- tum fluctuations, we also performed a mean field calcu- where the summation runs over the superlattice unit cell. lation for a comparison59. The critical interaction within In the iDMRG calculation, we introduce a finite bond the mean field approximation is found to be Vc = 0:78t dimension χ up to χ = 1600 as a cut-off to approximate and the interaction V = 0:30t corresponds to the same the ground state wavefunction in the form of a matrix coupling strength in terms of the normalized interaction product state, and we can obtain the true ground state g = (V − Vc)=Vc = −0:62. The iDMRG reuslts of M by a careful extrapolation to χ ! 1 from the finite χ (blue symbols) are larger than the corresponding mean 51{56 results (see also Appendix A). Generally, iDMRG field results (red symbols), M > MMF, which suggests with finite χ gives a good approximation especially when that quantum fluctuations enhance the order parameter the system considered is gapful. As we will show, an even for the present weak V . It is noted that the or- extrapolation to χ ! 1 works well, because our system der parameter behaves roughly as M(B) ∼ B as seen has a gap in presence of a non-zero B due to the magnetic for small B, which is consistent with the large Nf field 24{31 catalysis of the broken discrete symmetry Z2 where there theory . is no gapless Nambu-Goldstone mode. For a comparison, Similarly to the weak interaction case, the CDW order we also discuss a two-leg ladder system with Ly = 2 in parameter M calculated by iDMRG (blue symbols) is en- Appendix B. hanced at a strong interaction V = 2:0t > Vc compared 3 (a) iDMRG A. Scaling argument The enhancement of M(B) by the quantum fluctua- tions can be even more pronounced near the quantum critical point. Figure 3 shows the CDW order parameter mean field (a) (b) mean field iDMRG iDMRG mean field (b) iDMRG FIG. 2. (a) The CDW order parameter M at a weak cou- pling. The blue symbols are the iDMRG results at V = mean field 0:50t < Vc for Ly = 6 (squares) and Ly = 10 (circles), while the red symbols are the mean field results (V = 0:30t) for the same system sizes. (b) M at a strong coupling V = 2:0t > Vc calculated by iDMRG (blue) and V = 1:20t by the mean field approximation (red). The interactions for iDMRG and the mean field approximation correspond to the same value of FIG. 3. (a) The CDW order parameter M at the quantum the normalized coupling constant g. critical point V = Vc = 1:30t calculated by iDMRG together with the mean field result at V = 0:78t corresponding to g = 0. Definitions of the symbols are the same as in Fig. 2. (b) M in the log-log plot. The black solid line is the power to the mean field result MMF (red symbols) at the corre- law fitting M ∼ B0:355, while the black dashed line is the sponding interaction V = 1:20t (or equivalently g = 0:54) p large Nf result M ∼ B shown for the eyes. as shown in Fig. 2 (b). However, this is governed by the B = 0 values and increase of M(B) by the magnetic field at V = Vc = 1:30t (blue symbols) together with the mean is roughly comparable to that of MMF(B). The result field result for V = 0:78t (red symbols), corresponding to that M > MMF already at B = 0 is because they behave as M(V; B = 0) ∼ gβ with β ' 0:5 ∼ 0:6 < 13{8 while g = 0. Clearly, the iDMRG result is significantly larger βMF than the mean field result, and the enhancement is much MMF(V; B = 0) ∼ g with βMF = 1 near the quantum critical point, and these critical behaviors essentially de- stronger than that in the weak interaction case. There termine magnitudes of the CDW order parameters away are some deviations between the results for Ly = 6 and from the critical points. For B 6= 0, the order parameter Ly = 10 for small magnetic fields, B . 0:01B0, due to a behaves roughly as M(B) − M(0) ∼ B2 in agreement long magnetic length lB, and the CDW order gets more 24{31 strongly stabilized when the system size Ly increases with the large Nf field theory . from Ly = 6 to Ly = 10. This should be a general ten- dency since the CDW phase at B = 0 extends to a smaller interaction region when the system size increases58. From this observation, we can discuss scaling behaviors of the IV. NEAR QUANTUM CRITICAL POINT CDW order parameter in the thermodynamic limit as a function of B near the quantum critical point. Indeed, as shown in Fig. 3 (b), the calculated M except for the In this section, we discuss quantum criticality of the smallest values of B converge for different system sizes magnetic catalysis based on a variant of finite size scal- Ly = 6; 10, and M(B) exhibits a power law behavior for ing ansatzes. Then, we establish a global phase diagram 0:02B0 . B . 0:1B0. The finite size effects are neg- around the quantum critical point in the interaction- ligible in this range of the magnetic field, and further- magnetic field plane, in close analogy with the well- more the scaling behavior would hold for smaller mag- known finite temperature phase diagram near a quantum netic fields down to B = 0 in a thermodynamic system critical point.