Chiral Tricritical Point: A New Universality Class in Dirac Systems Shuai Yin,1 Shao-Kai Jian,1, ∗ and Hong Yao1, 2, 3, y 1Institute for Advanced Study, Tsinghua University, Beijing 100084, China 2State Key Laboratory of Low Dimensional Quantum Physics, Tsinghua University, Beijing 100084, China 3Collaborative Innovation Center of Quantum Matter, Beijing 100084, China (Dated: May 28, 2018) Tricriticality, as a sister of criticality, is a fundamental and absorbing issue in condensed-matter physics. It has been verified that the bosonic Wilson-Fisher universality class can be changed by gapless fermionic modes at criticality. However, the counterpart phenomena at tricriticality have rarely been explored. In this Letter, we study a model in which a tricritical Ising model is coupled to massless Dirac fermions. We find that the massless Dirac fermions result in the emergence of a new tricritical point, which we refer to as the chiral tricritical point (CTP), at the phase boundary between the Dirac semimetal and the charge-density wave insulator. From functional renormalization group analysis of the effective action, we obtain the critical behaviors of the CTP, which are qualitatively distinct from both the tricritical Ising universality and the chiral Ising universality. We further extend the calculations of the chiral tricritical behaviors of Ising spins to the case of Heisenberg spins. The experimental relevance of the CTP in two-dimensional Dirac semimetals is also discussed. Introduction.| Discovering and understanding uni- H⊥ versality classes is among the central issues in modern 0.6 ◼ condensed-matter physics as the fundamental concept 1st order of universality class is crucial to classifying the infrared 0.4 ◼ Dirac Semimetal CTP 2 physics in many-body systems [1, 2]. The universality nd 0.2 ◼ CDW order class is dictated by the renormalization group (RG) flow near the corresponding fixed point [3, 4]. According to ◼ ◼ ◼ ◼ ◼ ◼ ◼ 0 0.2 0.4 0.6 0.8 1.21.0 the number of relevant directions of the RG flow, the H singular behavior near the fixed point can be classified into criticality, tricriticality, and so on [5]. In contrast FIG. 1. The quantum phase diagram as a function of Zee- to the criticality controlled by only one relevant direc- man fields H? and Hk with h = 0:2; t = 1;J1 = 2J2 = 1=6. The blue dashed and blue solid lines refer to first- and tion, tricriticality often shows more intriguing properties, second-order phase transitions between the charge-density since it has two relevant directions [5, 6]. Thus, tricrit- wave (CDW) insulator and the Dirac semimetal. The red icality has aroused enormous theoretical as well as ex- point at (Hk;H?) = (0:7; 0:375) denotes the chiral tricritial perimental studies in condensed-matter systems [7{17]. point (CTP). The gray dashed and gray solid lines refer to For instance, the tricritical Ising fixed point in 1 + 1 first- and second-order phase transitions of the Ising field dimensions enjoys superconformal symmetry [8]. The without coupling to fermions (h = 0). quantum tricritical point (TP) connecting the deconfined quantum critical point [9] and first-order transition in two-dimensional spin systems exhibits O(4) symmetry different universality class, usually called the chiral Ising and self-duality [10]. Furthermore, even supersymmetric universality class (or the Gross-Neveu-Yukawa universal- quantum electrodynamics can emerge at the pair-density- ity class) if the Ising field couples to the Dirac fermions wave quantum TP on the surface of a three-dimensional as a chiral mass [18{28]. Recently, it was proposed in topological insulator [11]. Among them, some theoretical Ref. [29] by two of us and collaborators that the pres- predictions [12, 13] have been verified in experiments [14{ ence of gapless fermionic modes can even change the type 17]. of phase transition from first-order to continuous, realiz- ing the so-called fermion-induced quantum critical points arXiv:1711.10473v3 [cond-mat.str-el] 25 May 2018 On the other hand, quantum phase transitions in sys- (also see Refs. [30{32]). However, for the more colorful tems of interacting Dirac fermions have been investigated tricriticality, the effects induced by the massless Dirac intensively in recent years [3, 4]. Besides their fundamen- fermions have rarely been studied. Some questions then tal relevance in relativistic quantum field theories, such arise: To what extent is the tricriticality affected by the fermion systems also emerge in the low-energy sector of Dirac fermion? Will a new fixed point emerge? If so, various condensed-matter systems [1, 2]. It was realized how are we to characterize the critical properties in this that gapless fermionic excitations at quantum criticality new universality class? play an essential role in determining the critical behav- To answer these questions, in this Letter we report iors. For instance, coupling massless Dirac fermions to a first theoretical study of the chiral tricritical point the Ising quantum critical point results in a qualitatively (CTP). As suggested by its name, the CTP refers to a tri- 2 critical point coupled to massless Dirac fermions. More plicity, because m decouples from the MF Hamiltonian specifically, we first construct a microscopic model on the for J− = 0. Such a simplification will not qualitatively honeycomb lattice, where the Ising field couples to the change the physics discussed here, because the global Z2 Q x itinerant fermions as a staggered on-site potential. Us- symmetry|namely, the symmetry generated by i σi | ing mean-field (MF) calculations, we find that this mi- is already broken by the external longitudinal field H?. croscopic model features a zero-temperature transition As a function of H? and Hk, the Hamiltonian Hb ex- between a charge-density wave insulator (CDW) and a hibits a phase transition from hφi = 0 to hφi 6= 0 (break- Dirac semimetal as shown in Fig. 1, where the transition ing the lattice inversion or Z2 symmetry that exchanges can be either discontinuous or continuous depending on two sublattices), which is continuous or discontinuous two tuning parameters. The CTP emerges at the location depending on the external fields (see the Supplemental where first- and second-order transitions meet. Material for details) and a quantum tricritical point sep- To understand the universal behaviors of the CTP, arating these two types of transitions, as shown by the we carry out functional renormalization group (FRG) gray line in Fig. 1. The Landau free energy up to sixth calculations of the effective theory near the transition order near the transitions is given by (see the Supple- point. By solving the RG equations, we find four nontriv- P3 00 2n mental Material for details) fb = n=0 λ2nφ , where ial fixed points, corresponding to the Ising, chiral Ising, 00 2 2 3=2 2 00 2 2 00 λ2 / (H? + Hk ) − Hk J+, λ4 / Hk − 4H?, and λ6 tricritical Ising, and chiral tricritical Ising universality is a positive constant. The tricritical point is located classes, distinguished by the relevant directions. The at λ2 = λ4 = 0, which can be accessed by tuning two former three universal classes have been extensively ex- parameters, H? and Hk. Since the critical behaviors of plored [5, 6, 19, 21, 33], while the last one| i.e., the this quantum tricritical point are governed by a φ6 the- chiral tricritical Ising universality class|is new and is ory whose upper critical dimension is 3, the MF theory our main focus in the present Letter. Various critical should be reliable. Various critical exponents of this tri- exponents characterizing this novel universality class are critical Ising universality class are listed in Table I. obtained. We further study tricritical points of Heisen- To investigate a CTP, we consider massless Dirac berg transitions in Dirac semimetals [19, 21, 34{37] by fermions on the honeycomb lattice which couple to the extending the Ising fields to Heisenberg spins. The rel- Ising variables as follows: evance of these universality classes to real materials like X y X z y graphene or graphene-like systems will also be discussed. Hf = −t (ciscjs + H:c:) + h σi ciscis; (3) Lattice model and phase diagram.| It is known that hiji;s i;s the (2+1)-dimensional Ising model under both transver- y sal and longitudinal Zeeman fields can realize a tricritical where cis labels the operator creating a fermion on site i point at the phase boundary between first- and second- with spin polarization s, t is the hopping amplitude, and order transitions [7]. The Hamiltonian on the honeycomb h refers to the coupling strength between fermions and lattice is given by Ising variables. Note that in Eq. (3) there is an implicit P y chemical potential term −µ c cis. The condensa- X z z X z z X z X x i;s is Hb =J1 σi σj −J2 σi σj −H? σi −Hk σi ; (1) tion of the Ising field φ gives rise to staggered on-site po- hiji hhijii i i tentials for the fermions on the honeycomb lattice which where σ denotes the Ising spins defined at each site, and can gap out the Dirac fermions, as shown by the dis- h i and hh ii refer to nearest neighbors (NN) and next- persions of the fermions in the Ising symmetry-breaking ~ 2 2 2 nearest neighbors (NNN), respectively. J and J refer to phase: "(k) = ±{h φ +t [3+2(cos k1 +cos k2 +cos(k1 − 1 2 1 k ))]g 2 , where the lattice constant is set to be 1 for sim- the NN antiferromagnetic and NNN ferromagnetic cou- 2 p 1 3 ~ plings, respectively. H? and Hk refer to the longitudinal plicity, ande ^1 = (1; 0),e ^2 = ( 2 ; 2 ), and ki =e ^i · k.