Symmetric Mass Generation

The most well-known mechanism for fermions to acquire a mass is the Nambu-Goldstone-Anderson-Higgs mechanism, i.e. after a spontaneous symmetry breaking, a bosonic field that couples to the fermion mass term condenses, which grants a mass gap for the fermionic excitation. In the last few years, it was gradually understood that there is a new mechanism of mass generation for fermions without involving any symmetry breaking within an anomaly-free symmetry group. This new mechanism is generally referred to as the"Symmetric Mass Generation (SMG)."It is realized that the SMG has deep connections with interacting topological insulator/superconductors, symmetry-protected topological states, perturbative local and non-perturbative global anomaly cancellations, and deconfined quantum criticality. It has strong implications for the lattice regularization of chiral gauge theories. This article defines the SMG, summarizes current numerical results, introduces a novel unifying theoretical framework (including the parton-Higgs and the s-confinement mechanisms, as well as the symmetry-extension construction), and overviews various features and applications of SMG.

Global symmetry is a central concept in quantum field theories (QFTs). One of the most immediate implications of symmetry is that it restricts the terms that can appear in a field theory action in the path integral formulation, as the partition function Z must remain invariant under the symmetry transformation -the partition function invariance however can be up to an invertible complex phase factor Z → e iα Z known as a quantum anomaly [1,2]. The quantum anomaly associated with the global symmetry is also known as the 't Hooft anomaly [3], which has profound consequences. The invertible phase e iα that detects the anomaly also specifies a cobordism class of the partition function of one-higherdimensional invertible topological quantum field theory (invertible TQFT) [4] via the anomaly inflow [5,6].
The anomaly provides a concise and powerful organization principle to classify quantum field theories, and to dictate the influence of the ultraviolet (UV) kinematics on the infrared (IR) dynamics of a field theory, because the anomaly does not change under symmetric deformations of the theory, including the renormalization arXiv:2204.14271v1 [cond-mat.str-el] 29 Apr 2022 group (RG) flow. The preservation of the anomaly index from UV to IR is known as the anomaly matching. If a quantum field theory with a global symmetry G has a non-vanishing 't Hooft anomaly in G, its IR dynamics cannot be trivially gapped while preserving G. As a consequence, either the symmetry G is spontaneously broken in IR [7][8][9][10][11], or the anomaly must be saturated by anomalous gapless quantum criticality (e.g., the IR limit of the system corresponds to massless free theories or interacting conformal field theories [CFT]) or gapped topological order (e.g., the IR limit of the system is described by a topological quantum field theory [TQFT]) [12][13][14][15][16][17][18][19][20][21][22][23][24].
Many examples of anomalies involve massless fermions, such as the chiral fermion Adler-Bell-Jackiw anomaly [25,26]. The anomaly causes an obstruction to open a gap in the fermion spectrum (i.e. an obstruction to make the fermion correlation length finite) without breaking the symmetry. The symmetry-breaking mechanism to generate a mass gap is known as Nambu-Goldstone-Anderson-Higgs mechanism [27][28][29][30][31][32][33]. Here Nambu-Goldstone refers to the spontaneous symmetry breaking, while Anderson-Higgs refers to giving a mass by the elementary or composite boson condensation. However, there are instances where the anomaly vanishes for a collection of massless fermions, yet the symmetry is still restrictive enough to forbid any fermion bilinear mass term. In this case, although there is no obstruction towards gapping the fermions symmetrically, the mechanism to achieve the symmetric gapped state must go beyond the free-fermion (perturbatively free or weak-coupling) approach, which potentially leads to a non-perturbative strong-coupling approach in order to generate a finite excitation gap in the fermion many-body spectrum by non-trivial interaction effects. The strong-coupling [34] here refers to the coupling in the continuum field theory being nonperturbative, or the interaction energy being of the same order as the kinetic energy on the lattice scale (which may also be called the intermediate-strength interaction on a lattice). This phenomenon of gapping out massless fermions by interactions in an anomaly-free system without breaking the anomaly-free symmetry is now called the symmetric mass generation (SMG) [34,35].
The idea of gapping out massless fermions by interaction dates back to Eichten and Preskill [36] in an attempt to regularize chiral fermions on the lattice. The understanding of SMG is significantly deepened over the past few years, following the development in condensed matter theory regarding symmetry-protected topological (SPT) states [37][38][39]. SPT states are short-rangeentangled quantum many-body states respecting certain global symmetry G. The bulk of a SPT state is featureless (i.e. gapped, symmetric and non-degenerated). A non-trivial SPT state is most explicitly characterized by its non-trivial boundary features, which are endowed by the non-vanishing 't Hooft anomaly of G (or mixed G symmetry gauge-gravity anomaly in broader cases) in the boundary effective theory of the SPT state. The one-to-one correspondence between the bulk of a SPT state and its boundary anomaly [40][41][42][43][44] provides the basis to classify distinct SPT states by their distinct boundary anomalies.
The connection between SPT state and anomaly relates SMG to another topic: the interaction-reduced classification of fermionic SPT states. Fermionic SPT states [4,[45][46][47][48][49][50][51][52][53][54][55][56][57][58] are SPT states of fermionic systems with (at least) fermion parity symmetry Z F 2 (or more generally, invertible topological phases of fermions [4,57]). In the free-fermion limit [59][60][61][62][63], non-trivial fermionic SPT states are characterized by the symmetry-protected gapless fermion boundary modes, if the symmetry forbids any fermion bilinear mass on the boundary. However, the free-fermion analysis does not rule out the possibility of gapping out the boundary fermions by interaction though the SMG mechanism. Such situation can indeed happen when the boundary fermions are actually anomaly free, such that the bulk state should be classified as a trivial gapped state under symmetric interactions, even though it looks like a non-trivial SPT state in the non-interacting limit. This leads to the interaction-reduced classification of fermionic SPT states in the bulk, which is closely related to the SMG for gapless fermions on the boundary.
The first example of interaction-reduced classification was provided by Fidkowski and Kitaev [64,65] in (1+1)D fermion systems, where an explicit interaction was proposed to drive the SMG among fermion zero modes on the (0+1)D boundary. This motivated a sequence of works generalizing the discussion to (2+1)D [66][67][68][69], (3+1)D [15-17, 70, 71], and higher dimensions [72][73][74]. Studies along this direction reveals families of interactions that could potentially drive the SMG in different dimensions, paving ways for numerical verifications in concrete lattice models. It is also realized that the interaction must be carefully designed to drive the SMG: some symmetric interactions are helpful towards this goal, while other symmetric interactions are not [75].
Currently most known examples of interaction-reduced classifications all correspond to cancellations of nonperturbative global anomalies. However, there are also cases of SMG that involve cancellations of perturbative local anomalies. Important examples of such are the chiral fermions [76][77][78][79] in even dimensional spacetime. Regularizing chiral fermions on a lattice is a long-standing problem in the lattice gauge theory. The Nielsen-Ninomiya no-go theorem [80][81][82] states that it is not possible to gap the fermion doubler in a non-interacting local lattice model without breaking the chiral symmetry. One possibility to circumvent the no-go theorem is to consider interaction effects [36], which introduces fermion interaction to gap out the fermion doubler (the mirror fermion) in local lattice fermion models, leaving the normal fermion (the light fermion) untouched. However, the early attempts [79,[83][84][85][86][87][88][89][90][91][92] were not successful, either because certain anomalies were not carefully cancelled or because the appropriate gapping interaction has not been found. With a deeper understanding of the SMG mecha-nism, the problem was revisited by Wen [93] for Spin (10) chiral fermions in (3+1)D. The idea is further developed by subsequent works in the same dimension [94][95][96][97][98] as well as in lower-dimensional analogs [75,[99][100][101]. More recent numerical works have successfully shown that the SMG indeed provides a feasible solution to regularize chiral fermions [102][103][104][105].
The article is organized as follows. We will start by introducing some selective representative models of SMG (one in each spacetime dimension, in (0+1)D, (1+1)D, (2+1)D and (3+1)D respectively) in Sec. II. We conclude with Sec. II E by providing a general definition of SMG in all dimensions. We then review the numerical efforts in Sec. III, which is mainly focused on two tasks: (i) to establish the existence of SMG phases in Sec. III A and (ii) to investigate the nature of SMG transitions in Sec. III B. Based on these backgrounds, we then summarize the recent theoretical progress in Sec. IV, which aims to (i) understand the SMG phase by the fluctuating bilinear mass picture in Sec. IV A and (ii) describe the SMG transition by fermion fractionalization field theory in Sec. IV B. In particular, we unify two currently existing SMG mechanisms (namely the parton-Higgs mechanism and the s-confinement mechanism) under the same theoretical framework of fermion fractionalization. We also make connection to further understanding of SMG based on the symmetry extension construction [23,24] in Sec. IV D. In Sec. V, we discuss other aspects of SMG including Green's function zeros Sec. V A, the deconfined quantum criticality Sec. V B, and the Standard Model regularization Sec. V C. Finally, we summarize our review in Sec. VI.

A. (0+1)D SMG: Fidkowski-Kitaev Majorana Fermion Model
The simplest example of SMG happens in (0+1)D spacetime among a collection of Majorana fermion zero modes, as first shown by Fidkowski and Kitaev [64,65]. The model concerns a system of eight Majorana fermion modes, described by the Majorana fermion operators χ a (a = 1, 2, · · · , 8) satisfying {χ a , χ b } = 2δ ab . Consider an anti-unitary (time-reversal) symmetry Z T 2 : χ a → χ a , i → −i and the fermion parity symmetry Z F 2 : χ a → −χ a . Without involving interactions, the Hamiltonian must take a fermion bilinear form to preserve the fermion parity symmetry Z F 2 . But any Majorana fermion bilinear term iχ a χ b (one needs i to keep the operator Hermitian) will break the time-reversal symmetry Z T 2 , so the freefermion Hamiltonian has to vanish, e.g. H = 0, under both symmetry requirements, and these Majorana modes cannot be gapped on the free-fermion level. Eight Majorana fermion modes form four qubits, hence the system has a 2 4 = 16 fold degeneracy, and this degeneracy is protected by the Z T 2 × Z F 2 symmetry.
However, it is possible to create a many-body excitation gap by fermion interaction, leaving a unique ground state of this (0+1)D system. As introduced by Fidkowski and Kitaev [64], the following four-fermion interaction suffices to gap out all eight Majorana fermion modes without breaking the time-reversal symmetry Z T The coefficient is specified by V abcd = 4e|χ a χ b χ c χ d |4e , where |4e = (|0000 + |1111 )/ √ 2 is a many-body reference state written in the Fock state basis |n 1 n 2 n 3 n 4 , labeled by the fermion occupation numbers n i = (1 + iχ 2i−1 χ 2i )/2 (for i = 1, 2, 3, 4). Here the fermion number operator n i = c † i c i can also be expressed in terms of the complex fermion annihilation operator c i = (χ 2i−1 + iχ 2i )/2, which are constructed by pairing up the Majorana operators. The reference state |4e describes the quartet condensation of the complex fermions c i (where four fermions are created or annihilated together in the quantum superposition), which is also known as the charge-4e superconducting state [106][107][108][109][110][111][112] in condensed matter physics. It turns out that the state |4e respects the Z T 2 symmetry, and is the unique ground state of H FK with a finite excitation gap of 14 (energy units).
The fact that the system has eight Majorana zero modes is crucial for the SMG to occur. In contrast, if there are only four Majorana zero modes, the only Hamiltonian allowed by symmetry is where each pair of Majorana fermion operators defines a fermion parity operator P F i = iχ 2i−1 χ i = 2n i − 1 (associated with the ith complex fermion mode). This Hamiltonian always has a two-fold ground state degeneracy, regardless of the sign of g. When g < 0, the ground states |P F 1 = P F 2 = ±1 are bosonic (as the total fermion parity P F 1 P F 2 = +1 is even). When g > 0, the ground states |P F 1 = −P F 2 = ±1 are fermionic (as P F 1 P F 2 = −1). In either cases, the two-fold degenerated ground states transform into each other under Z T 2 , and form a Kramers doublet [113]. The ground state degeneracy implies that the Z T 2 symmetry is spontaneously broken. Hence even with interactions, four Majorana zero modes still cannot be symmetrically gapped. However, by doubling to eight Majorana fermions, it is then possible to couple the two Kramers doublets (originated from the first-and the lastfour Majorana fermions under interaction) together via a Heisenberg-type spin-spin interaction as [72,114] where the first two terms stabilizes two Kramers doublets, and S I/II stands for the effective spin-operator for each Kramers doublet. This will end up with a unique spin-singlet bosonic ground state with a finite gap to all excitations, which successfully gaps out all fermions without breaking the symmetry. Eq. (3) can be rewritten as H FK in Eq. (1) under some appropriate basis choice.
From the perspective of quantum anomaly, the (0+1)D Majorana fermions with a Z T 2 ×Z F 2 [115] internal symmetry (or the Pin − spacetime-internal symmetry [48,114], or the BDI symmetry class [116,117]) has a Z 8 class of non-perturbative global anomaly, whose anomaly index corresponds to the number of Majorana modes. With eight Majorana modes, the anomaly vanishes, meaning that the system can be trivially gapped without breaking symmetry. However, the Z T 2 symmetry is restrictive enough to rule out any fermion bilinear mass, making interaction a necessary ingredient in the fermion mass generation, which corresponds to the SMG mechanism.
The decomposition of the four-fermion interaction as a product of two fermion-bilinear operators Φ α Φ α in both Eq. (4) and Eq. (5) suggests a common physical picture to understand the SMG as fluctuating bilinear masses. Under the Hubbard-Stratonovich transformation, the SMG interactions in the above examples take the general form of Yukawa-Higgs interaction where φ α is a bosonic Yukawa field that couples to the fermion bilinear mass, and by integrating out φ α , the desired four fermion interaction is generated. If the Yukawa field condenses, i.e. φ α = 0, the fermions will be gapped (with a gap size proportional to the amplitude of φ α ) and the symmetry is also broken spontaneously. However, if it is possible to fluctuate the orientation of the Yukawa field smoothly in the spacetime (which is only the time here for the (0+1)D system) without bringing its local amplitude to zero, the fermion could potentially retain the excitation gap while restoring the required symmetry. This intuitive picture leads to fruitful understandings of SMG in higher dimensions [16,100,118], which will be further discussed in Sec. IV A. The SMG in (0+1)D is intimately related to the interaction-reduced classification of fermionic SPT states in (1+1)D. In the non-interacting limit, the (0+1)D Majorana zero modes can be viewed as the topological edge modes of a (1+1)D fermionic SPT state protected by the same Z T 2 × Z F 2 internal symmetry. A specific lattice model for such SPT root state is Kitaev's Majorana fermionic chain [119], which supports a dangling Majorana zero mode on each open end. So eight Majorana zero modes can be viewed as the boundary state of eight copies of Majorana chains. Without interaction, the Z T 2 × Z F 2 symmetric (BDI symmetry class) free fermionic SPT states are Z classified [59][60][61], where the index corresponds to the number of Majorana edge modes [120]. The fact that eight Majorana zero modes can be trivially gapped out by interaction without breaking the symmetry implies that eight copies of Majorana chain actually belong to the trivial SPT phase, as their interface with the vacuum state can be made featureless (i.e. gapped and non-degenerated) by the SMG. This indicates that the Z T 2 × Z F 2 (or Pin − ) symmetric interacting fermionic SPT states are Z 8 classified [45,64,65], which is consistent with the formal result Ω Pin − 2 (pt) = Z 8 by the cobordism classification [48,121]. The phenomenon that the fermion SPT classification is reduced from Z in the non-interacting limit to Z 8 under interaction is called the interaction-reduced classification.
The Z → Z 8 interaction-reduced classification implies that eight copies of the Majorana chain can be smoothly tuned (without closing the bulk gap) to the trivial phase under interaction. The bulk model is defined on a onedimensional lattice where each site i hosts eight Majorana modes χ i,a (a = 1, · · · , 8). In the free-fermion limit (g = 0), the u < v and u > v phases are separated by a quantum phase transition at which the single-particle band gap closes in the bulk. Ref. 64 shows that the quantum critical point can be circumvented by applying the Fidkowski-Kitaev interaction on every site with strength g (where V abcd in H int follows the same definition as that in Eq. (1)). This can be argued by inspecting the limiting cases (u = 0, v = 1) or (u = 1, v = 0) where Majorana fermions are fully dimerized along the lattice over even or odd bonds, see Fig. 1. The two dimerized states can both be smoothly tuned to the same strong coupling state i |4e i as g → ∞ without closing the many-body gap, which can be explicitly verified by diagonalizing the local Hamiltonian across a bond. This establishes a smooth deformation path between the u < v and u > v phases by going through the strong coupling regime as shown in Fig. 1, demonstrating that eight copies of the Majorana chain is indeed in the trivial phase. Moreover, the argument also implies that the gapless Majorana fermions (in the (1+1)D bulk) at the freefermion critical point (u = v) can be gapped out by turning on the interaction g (see Fig. 1), leading to another example of SMG [122]. The bulk SMG of the (1+1)D Majorana chain (at the critical point) is dynamically equivalent to the SMG for the boundary mode of a (2+1)D fermionic topological superconductor (TSC) made of eight copies of the p x ± ip y superconductor [123][124][125]. In general, each interaction-reduced classification of the fermion SPT state implies the existence of SMG mechanisms for the gapless fermions both on the boundary and at the bulk critical point between the SPT and the trivial state, where the same local interaction that gaps out the boundary mode can be used to gap out the bulk critical point as well. Using the connection of SMG phenomena between the boundary and the bulk, Ref. 72 obtains a series of interaction-reduced classification of fermionic SPT states in all spacetime dimensions systematically, hence extending the phenomenon of SMG to higher dimensions.

B. (1+1)D SMG: 3-4-5-0 Chiral Fermion Model
The 3-4-5-0 U(1) chiral fermion model by Wang-Wen [75,101] provides an instructive example of SMG in (1+1)D, which involves the cancellation of a perturbative anomaly and is not related to interaction-reduced fermionic SPT classifications. The (1+1)D chiral fermion theory is described by the Lagrangian where ψ a are complex fermion fields (for a = 1, 2, 3, 4) with different velocities v = (1, 1, −1, −1). The model can be assigned a U(1) symmetry, under which the fermions transform as U(1) : ψ a → e iθqa ψ a with q = (3, 4, 5, 0). The two left-moving Weyl fermions ψ 1 and ψ 2 are assigned with charge 3 and 4, and the two rightmoving Weyl fermions ψ 3 and ψ 4 are assigned with charge 5 and 0, hence the name of the 3-4-5-0 model. The seemingly peculiar charge assignment of the chiral fermions is designed to cancel the U(1) 't Hooft anomaly. The (1+1)D fermion with U(1) internal symmetry (or Spin c ≡ (Spin × U(1))/Z 2 spacetime-internal symmetry [48,114], or the A symmetry class [116,117]) has a Z perturbative anomaly, whose anomaly index is given by a v a q 2 a , which vanishes for the charge assignment of the 3-4-5-0 model. The (1+1)D chiral fermions can be also viewed as the chiral edge modes of a (2+1)D integer quantum Hall state, also known as the fermionic SPT state in symmetry class A. The integer quantum Hall states are Z classified by the quantized Hall conductance, and the classification does not further reduce under interaction. The bulk Hall conductance (in unit of e 2 /h) matches the boundary anomaly index.
The fact that the U(1) 't Hooft anomaly vanishes for the 3-4-5-0 model indicates that it should be possible to trivially gap out all fermions without breaking the U(1) symmetry. However, the U(1) symmetry is restrictive enough to prevent the gapping to happen on the freefermion level, because any fermion bilinear term that produces a gap must take the form of ψ † a ψ b or ψ a ψ b (with a ∈ {1, 2} and b ∈ {3, 4}) that mixes the left-and rightmoving fermions. Since the four flavors of fermions all carry distinct U(1) charges that do not add or subtract to zero, any flavor mixing fermion bilinear term will necessarily break the U(1) symmetry, which makes it impossible to symmetrically gap out these chiral fermions on the free-fermion level.
Nevertheless, the gapping can be achieved with fermion interaction, hence an example of SMG. The interaction to achieve the SMG would involve at least sixfermion terms, which can be derived using the null-vector condition [126][127][128][129][130] for quantum Hall edge states. By bosonization ψ a ∼ e iϕa , the (1+1)D chiral fermion system can be effectively described by a Luttinger liquid theory where K = diag(1, 1, −1, −1) and V = diag(1, 1, 1, 1) are diagonal matrices. In the bosonization language, backscattering fermion interactions can be introduced as each interaction labeled by a charge vector l α . It is possible to find two charge vectors l 1 = (1, −2, 1, 2) and l 2 = (2, 1, −2, 1) that satisfy the null-vector condition l α K −1 l β = 0 (for α, β = 1, 2). The null-vector condition ensures that the vertex operators O α = e il α ϕ are both self-boson and mutual-boson, which can be simultaneously condensed at large coupling g α . Once these vertex operators condense O α = 0, all fermion excitations will be gapped, since none of the fermion operator braid trivially with any of the vertex operator (as seen from l α K −1 = 0). Furthermore, the interaction term L int , as well as the condensate O α , preserves the U(1) symmetry since l α q = 0 for α = 1, 2. Therefore, L int provides a symmetric way to gap out all chiral fermions in the 3-4-5-0 model, realizing the SMG. In terms of the fermion field, the interaction can be translated from Eq. (10) to where the ∂ x operator (that generates infinitesimal translation) is inserted as a point-splitting regularization to avoid identical fermion operators appearing at the same spatial position. The multi-fermion interaction in Eq. (11) can be mediated by two independent Yukawa fields φ 1 , φ 2 via Integrating out the Yukawa fields [131] will generate the interaction in Eq. (11) to the leading order of g α ∼g 2 α (plus additional density-density interactions like ψ † 1 ψ 1 ψ † 3 ψ 3 or ψ † 2 ψ 2 ψ † 4 ψ 4 whose effect is only to renormalize the V matrix in the Luttinger liquid theory and can be safely ignored). The Yukawa fields φ 1 and φ 2 carry the U(1) charges −4 and −2 respectively. Directly condensing the Yukawa fields would provide Dirac/Majorana masses to all chiral fermions ψ a at the price of breaking the U(1) symmetry. Nevertheless, the SMG mechanism suggests an alternative scenario that these Yukawa fields are fluctuating in the disordered (a.k.a. strong-coupling symmetric) phase, such that the U(1) symmetry remains unbroken but the chiral fermions could still acquire a spectral gap via the Yukawa interaction in suitable parameter regimes. Similar idea was numerically explored in Ref. 92 without success, due to the incorrect design of the fermion interaction (which does not satisfy the nullvector condition). The correct design of fermion interaction in Eq. (11) or Eq. (12) includes only a restricted subset of symmetry-allowed interactions [75], which are helpful for the SMG. More recent numerical study in Ref. 105 has confirmed that the correct interaction indeed leads to the SMG phase.
The SMG transition happens when the interaction is beyond a finite critical strength, because the interaction is perturbatively irrelevant at the free-fermion fixed point, due to its high-order nature. However, strong enough interaction could lead to non-perturbative effects. Increasing the interaction strength generally tunes the Luttinger parameter (by renormalizing the V ma-trix) and alters operator scaling dimensions in the Luttinger liquid theory [105]. When the scaling dimension of the interaction term itself is tuned to marginal, the SMG transition will be triggered, which drives the system from the gapless phase to the featureless gapped phase [34,105]. In this case, the SMG transition will belong to the Berezinskii-Kosterlitz-Thouless (BKT) universality class.
1. Proof on the equivalence between the anomaly-free and SMG gapping conditions Although the earlier discussions focus on the SMG of (1+1)D 3-4-5-0 U(1) chiral fermion model, there is no obstacle to generalize to show any (1+1)D anomalyfree chiral fermion model with multiple U(1) symmetries can allow SMG, following [75]. Ref. 75 proves that the anomaly-free condition of the (1+1)D multiple U(1) chiral fermion theory is equivalent to the SMG gapping condition of the same theory. The proof can be achieved due to the exact bosonization-fermionization techniques in (1+1)D. Here we recall the proof to complete the discussion of (1+1)D SMG.
A generic (1+1)D anomaly-free multiple-U(1) chiral fermion theory have equal numbers of left and right moving Weyl fermions, N L = N R = N ∈ Z + , such that the total Weyl fermion number is 2N = N L + N R = 2Z, an even integer.
2. Its anomaly-free condition, on the other hand, requires: • Gravitational anomaly free (two-point one-loop Feynman diagram of grav 2 vertices vanish): The left and right chiral central charges c L = c R , which means N L = N R .
• Gauge anomaly free (two-point one-loop Feynman diagram of U(1) 2 vertices vanish): For each U(1) symmetry, with left-handed and right-handed Weyl fermion charge vector q L and q R respectively, the anomaly-free requires the square sum of each component q 2 L −q 2 R = 0. In terms of the symmetric bilinear form K for both bosonic (K b ) and fermionic (K f ) systems, the anomalyfree condition demands that again for α, β ∈ {1, 2, . . . , N } and a, b ∈ {1, 2, . . . , 2N }.
• The above two anomaly-free conditions are perturbative local anomalies. The (1+1)D nonperturbative global anomalies are classified by cobordism groups (denoted TP 3 or Ω 3 with the special orthogonal SO or Spin group and some internal U(1) symmetries, such SO × U(1), Spin × U(1) and Spin c ) which turn out to always vanish [132].
3. To explain that the SMG gapping condition holds implies the anomaly-free condition also holds, Ref. 75 shows that given the set of l α satisfying (15), one can find the set of q α simultaneously satisfying (16) and (17). This is true because that given l α , we can choose q α = Kl α = K −1 l α thanks to the property K = K −1 .
4. To explain that the anomaly-free condition holds implies the SMG gapping condition also holds, Ref. 75 shows that given the set of q α satisfying (17), one can find the set of l α simultaneously satisfying (16) and (15). This is true because that given q α , we can choose l α = Kq α = K −1 q α thanks to K = K −1 .
5. The above two remarks prove that the if and only if (sufficient and necessary) conditions to the equivalence of the anomaly-free condition and the SMG gapping condition. Once the set of q α and l α are found, they form 2N linear-independent integer-valued vectors spanning completely the 2N -dimensional vector space (known as the Narain lattice [133]). 6. What remains to be explained is why the SMG gapping condition defines a gapped boundary without any topological boundary degeneracy [128]. The idea is viewing the (1+1)D theory (14) as the boundary theory of (2+1)D invertible TQFT with a Chern-Simons ac- The A a is a multiplet 1-form gauge field. Hereafter all repeated indices are summed over. A stable boundary condition requires the variation of S bulk on the boundary 2-manifold (∂M) 2 vanished [127] under the boundary 1-form gauge field The differential δ of this variation is a symplectic form ω Sp = K ab 4π (∂M) 2 δA ∂,a ∧δA ∂,b on the space of boundary gauge fields. Consistent stable boundary conditions on (∂M) 2 defines a Lagrangian submanifold with respect to the symplectic form ω Sp in symplectic geometry.
• One consistent boundary condition sets one component of A ∂ vanished, such as K ab A b,t − V ab A b,x | ∂M = 0, which gives a gapless (1+1)D CFT (14).
• Another boundary condition sets the gauge degrees of freedom l α,a A a | ∂M = 0 vanish [53]. The boson modes ϕ a , originally related by the gauge transformation A a → A a + dλ a and ϕ a → ϕ a − λ a , now may condense on the boundary with nonzero vacuum expectation values exp i(l α,a ϕ a ) | ∂M = 0, more precisely, indeed exp i( lα,a | gcd(lα)| ϕ a ) | ∂M = 0, where gcd(l α ) ≡ gcd(l α,1 , l α,2 , . . . , l α,2N ) is the greatest common divisor (gcd) of the all components of l α . This condensation of ϕ can be triggered by the earlier sine-Gordon cosine term at a strong g coupling. The boundary vertex operator and bulk line operator are connected exp i(l α,a ϕ a )| ∂M + i l α,a A a | M . The gapped bulk and gapped boundary demand that the partition function Z evaluated on the 3-manifold M with the 2boundary ∂M has a finite value (in fact Z = 1 when the Z corresponds to counting the dimension of the Hilbert space for the invertible TQFT). This means that arbitrary link configuration of the bulk line operators should give a trivial braiding statistical phase to Z, so there are no unwanted quantum fluctuations destabilizing the gapped vacuum -namely, the mutual statistics exp(i2πl α K −1 l β ) = 1 and the self statistics exp(iπl α K −1 l α ) = 1 are trivial for all α, β ∈ {1, 2, . . . , N }. Hence we derive the correspondence between the N independent compatible SMG gapping terms and the N null-braiding statistics l α vectors [126][127][128][129]. This completes the proof [75]. The honeycomb lattice model is a simple lattice model for SMG in (2+1)D. The model is defined on a honeycomb lattice, with each lattice site i hosting four complex fermion modes, denoted as c ia (for a = 1, 2, 3, 4). The model is described by the Hamiltonian (18) where ij stands for the bond between nearest neighboring sites i and j on the honeycomb lattice, see Fig. 2(a). The interaction strength g is the only tuning parameter of this model. When g = 0, the free-fermion hopping on the honeycomb lattice gives rise to the graphene band structure [134], as shown in Fig. 2(b), that produces 4 × 2 = 8 gapless Dirac fermions ψ Qa at low energy, where 4 stands for the four internal flavors (a = 1, 2, 3, 4) and 2 comes of the fermion doubling (Q = K, K ) in the Brillouin zone. They can be described by a low-energy effective field theory Lagrangian where Q = K, K labels the valley (fermion doubling) freedom, γ µ = (σ 2 , σ 1 , σ 3 ) and ψ Qa are complex Grassmann spinors in (2+1)D (withψ Qa = ψ † Qa γ 0 ). They can also be viewed as 16 gapless Majorana fermions equivalently, by decomposing each complex Grassmann field into two real Grassmann fields.
The key physical symmetries that protect these gapless fermions in the non-interacting limit are the lattice translation symmetry and the anti-unitary sublattice symme- stands for a sign factor that takes ± on A/B sublattice of the honeycomb lattice) [59,61,63]. These two symmetries can be combined to create an emergent anti-unitary symmetry at low-energy Z T F 4 : ψ Qa → iγ 0 ψ † Qa , i → −i [135], whose generator T = T 3/4 R S consists of a sublattice symmetry Z S 2 generator S followed by a 3/4 fraction of the unit-cell translation T R . Although the 2D lattice translation symmetry Z 2 (generated by T R along two linearly independent Bravais lattice vectors R) cannot be fractionalized on the lattice, yet for the low-energy effective theory Eq. (19), the translation symmetry acting on the low-energy fermions ψ Q becomes an emergent valley U(1) symmetry T R : ψ K/K → e ±i2π/3 ψ K/K (± signs are associated with K and K valleys respectively), which can be fractionalized to T 3/4 R : ψ K/K → ±iψ K/K . The anomalous nature of the combined symmetry Z T F 4 is manifested by the fact that Z T F 4 is only an emergent symmetry at low-energy and becomes ill-defined on the lattice level.
The (2+1)D Majorana fermion with a Z T F 4 internal symmetry (or the Pin + spacetime-internal symmetry [48,114], or the DIII symmetry class [116,117]) has a Ω Pin + 4 (pt) = Z 16 non-perturbative global anomaly based on the cobordism, whose anomaly index corresponds to the number of gapless Majorana fermions. The honeycomb lattice model precisely has 16 Majorana fermions at low-energy, which is free of the Z 16 non-perturbative global anomaly. However, the Z T F 4 symmetry is still restrictive enough to rule out all possible fermion bilinear gapping terms from appearing in Eq. (19), which again calls for the SMG mechanism.
The four-fermion interaction g in Eq. (18) is one choice of the interaction that drives the desired SMG. On every site, this interaction is the same as the SU(4) symmetric interaction in Eq. (5) (up to unimportant density-density interactions, which do not affect the ground state but only to renormalize the gap size). The interaction g explicitly drives a four-fermion condensation, also known as the charge-4e superconducting order [106][107][108][109][110][111][112]. In the large g limit, the many-body ground state of the system is simply the product i |4e i of on-site ground states |4e i = (|0000 i +|1111 i )/ √ 2, which is symmetric, non-degenerated and gapped, realizing the SMG phase. Therefore, by tuning the strength g, one expects to drive an SMG transition at some intermediate g ∼ 1 (that is comparable with the bandwidth of the lattice fermion). Various numerical simulations of this model (and its variants) [118,[136][137][138][139] have suggested the existence of such a direct and continuous SMG transition between the gapless and the gapped phases in (2+1)D. A field theory description of the SMG quantum critical point was proposed in Ref. 35, which will be further reviewed in Sec. IV C 1.
The honeycomb lattice model also has an explicit SU(4) symmetry c ia → U ab c ib that rotates fermions among the four on-site flavors. Although the SU(4) symmetry does not affect the anomaly analysis in any essential way, it could help to remove all the SU(4)-breaking relevant perturbations (if there were any) at the critical point, which might help to promote a continuous SMG transition. Thus the SU(4) symmetry can also be included in the discussion, and the total internal symmetry becomes SU(4)× Z F 2 Z T F 4 (where SU(4) and Z T F 4 share the fermion parity Z F 2 subgroup).
The decomposition of the on-site SU(4)-symmetric interaction into a product of fermion bilinear operators Φ α i Φ α i as in Eq. (5) indicates that one can again introduce a Yukawa field φ α to mediate the fermion interaction, [35] whereΓ α = (σ 12 , σ 20 , σ 32 , iσ 21 , iσ 02 , iσ 23 ) are the flavorspace matrices for the flavor-sextet pairing that transform as representation 6 in SU (4). Importantly, the Yukawa field must transform as 4 , meaning that directly condensing the Yukawa field φ α = 0 will necessarily break the protecting symmetry Z T F 4 (as well as breaking the standing-by SU(4) symmetry). The SMG interaction thus provides a mechanism to allow the Yukawa field to condense locally without establishing long-range order, so as to maintain the fermion gap without breaking the symmetry.

D. (3+1)D SMG: Chiral Fermion Model
Lattice regularization of chiral fermions in (3+1)D has always been an important motivation to study SMG in lattice gauge theories. A rich class of SMG transitions has been proposed and analyzed in different 4d chiral fermion models recently [34,98]. Here we will review one simplest example from [34] to illustrate the essential features of these models. The example considers a collection of right-handed massless Weyl fermions ψ transforming under an internal SU(N ) symmetry in the ⊕ (N + 4)r epresentation, described by the Lagrangian where ψ is the collection of complex Grassmann spinor fields. More explicitly, the fermion field can be split into ψ = ψ λ ⊕ ψ ψ , such that ψ λ denotes the single Weyl fermion and ψ ψ denotes the (N + 4) multiples of¯Weyl fermions. The goal is to find a path of gapping these fermions without breaking the SU(N ) symmetry.
First of all, the symmetric gapping is possible because all fermions together cancel the SU(N ) 't Hooft anomaly, as can be verified by the following anomaly index calculation Secondly, the SU(N ) symmetry is restrictive enough to forbid all possible fermion bilinear masses. For chiral fermions, the only available fermion bilinear mass is the Majorana mass that pairs up fermions like (ψ iσ 2 ψ + h.c.). Since the Majorana mass is already antisymmetric in the spinor subspace, the flavor subspace must be symmetric. However, the fermion representation of the SU(N ) symmetry guarantees that its symmetric product does not contain a trivial representation 1, as So any bilinear mass term will transform non-trivially under SU(N ) and cannot be condensed in the presence of the SU(N ) symmetry. Therefore, one must rely on fermion interactions to achieve symmetric gapping, i.e. the SMG mechanism.
To design the appropriate interaction to realize the SMG transition, one might search for candidate fourfermion interactions by looking for trivial representations in ( ⊕ (N + 4)¯) ×4 . However, ( ⊕ (N + 4)¯) ×4 → 1 in general (unless for N = 2, 4), meaning that one might need to look into higher-order fermion interactions, which are even more irrelevant and less interesting to explore. Nevertheless, the three-fermion combination ( ⊕ (N + 4)¯) ×3 contains trivial representations of the form (ψ λ iσ 2 ψ ψ )ψ ψ , but this term is fermionic and cannot appear directly in the Lagrangian (otherwise the fermion parity symmetry would be broken). An interesting idea from [34,98] is to bring down (N +4)(N +3)/2 additional fermions ψ χ from the high-energy spectrum which transforms trivially under the SU(N ) symmetry, such that a four-fermion interaction can be constructed. In this interaction, the fermion flavors are contracted in the way such that ψ λ iσ 2 ψ ψ and ψ χ iσ 2 ψ ψ transform as and¯under SU(N ) respectively.  It is instructive to note that the chiral fermion model has an additional internal symmetry SU(N + 4), under which ψ λ , ψ ψ and ψ χ transform as 1,¯and respectively, as summarized in Tab. I. The representations are so assigned that the fermions are also free of the SU(N + 4) anomaly, as Moreover, there is no mixed anomaly between SU(N ) and SU(N + 4). The interaction L int in Eq. (23) also respects the larger SU(N )×SU(N +4) symmetry. Ref. 34 suggests that L int is a plausible SMG interaction that drives all fermions ψ = ψ λ ⊕ψ ψ ⊕ψ χ to the gapped phase without breaking the SU(N )×SU(N +4) symmetry. The mechanism will be further reviewed in Sec. IV C 2.
The interaction L int in Eq. (23) can be decomposed into the following Yukawa couplings where the Yukawa fields φ 1 and φ 2 are both in the (¯, ) representation of SU(N )×SU(N +4). The above Yukawa decomposition again provides an alternative way to understand the SMG transition by driving the Yukawa field into a strong-coupling symmetric phase.

E. SMG in General Dimensions
For SMG to happen in a fermion system with a choice of internal symmetry G in d-dimensional spacetime, the following two necessary conditions must be satisfied: • The anomaly index ν ∈ TP d+1 (G) of the system must vanish (ν = 0), where TP d+1 (G) denotes the classification of invertible topological phases (with low-energy invertible topological field theories) in (d + 1)-dimensional spacetime [4]. (We may also denote d + 1 ≡ D as the bulk dimension.) • The single fermion must be in a representation r G ψ of G, such that the antisymmetric product (denoted by × A ) representation r G ψ × A r G ψ does not contain the trivial representation 1 G in its direct sum decomposition.
If the first condition is violated (ν = 0), trivially and symmetrically gapping out the fermions is impossible with or without interaction, due to the anomaly obstruction. However, it may be possible to symmetrically gap out the fermions with a symmetric anomalous topological order described by its underlying low-energy topological quantum field theory (TQFT), see Sec. IV D. If the second condition is violated, a trivial and symmetric gap can already be achieved at the free-fermion level by a fermion bilinear condensation (through the trivial channel given by r G ψ × A r G ψ → 1 G ), such that fermion interaction is not necessary. When both conditions are satisfied, an interacting (strong-coupling) mechanism is needed to generate the fermion mass without breaking the G symmetry, i.e. the SMG mechanism.
These conditions are also sufficient for the existence of an interacting SMG phase (i.e. the featureless gapped phase). However, they do not necessarily imply a direct and continuous SMG transition of fermions from the gapless phase to the gapped phase. A general design of the four fermion interaction that leads to a single continuous transition between the gapless fermions and the SMG remains a challenging problem. Much of the progress relies on numerical simulations, as to be reviewed in the next section Sec. III.

A. Existence of SMG Phases
Because the SMG is generally a non-perturbative interaction effect (with respect to the gapless free-fermion fixed point), numerical simulations play important roles in the study of SMG. Given an anomaly-free gapless fermion system with a proposed gapping interaction, the numerical study has two major goals: • Establish the existence of the SMG phase in the strongly interacting limit.
• Investigate the nature of the SMG transition at the critical interaction strength.
The SMG phase refers to the strongly interacting trivial gapped state of fermions without symmetry breaking. It was also called the paramagnetic strong-coupling (PMS) phase [137,140], the strongly coupled symmetric gapped (SCSG) phase [94,95], or the featureless gapped phase [35,141] in the literature. Correspondingly, the gapless Dirac/Weyl/Majorana fermion phase in the weak interaction regime is sometimes also called the paramagnetic weak-coupling (PMW) phase [137,140], or the semi-metal (SM) phase [35,141]. The phase transition separating the SMG phase and the free-fermion phase is called the SMG transition.
where ψ ia denotes the fermion of the flavor a on the site i (or in the unit cell i). The free-fermion Hamiltonian H 0 generates gapless fermions at low-energy, and the interaction Hamiltonian H int applies properly designed interactions to drive the SMG. The interaction strength will be generally denoted as g in the following discussion. Different models are mainly distinct in the following aspects: • The low-energy fermions can be realized in the lattice model H 0 either as the gapless boundary modes of a fermionic SPT state in one higher spacetime dimension [105], or as the gapless bulk mode of a semi-metal state in the designated spacetime dimension (such as the honeycomb lattice fermion in (2+1)D [136,139] and the stagger fermion in general dimensions [140,143]).
• The interaction H int can either be explicit given by multi-fermion local interaction terms (denoted as ψ 4 for four-fermion interactions or ψ 6 for sixfermion interactions in Tab. II), or mediated by intermediate Yukawa-Higgs fields (denoted as YH) or non-Abelian gauge fields (denoted as QCD).
• In even spacetime dimensions, the interaction can be restricted to part of the fermions (forming the mirror sector) to study the SMG in chiral fermion systems [104,105].
The SMG phase has been successfully achieved in studies listed in Tab. II, which demonstrate the generality of the SMG mechanism in different spacetime dimensions with various symmetry assignments and under different forms of interaction. The SMG phase has the following defining features that can be checked in numerics: • Unique ground state with a gap to all excitations (including both fermionic and bosonic excitations), • Absence of spontaneous symmetry breaking (no fermion bilinear condensation), • Formation of the four-fermion (or higher multifermion) condensate that preserves the symmetry.
Directly computing the many-body excitation gap is challenging for most numerical approaches. One way to probe the excitation gap is to measure the correlation function in different channels An exponential decaying correlation function implies a finite excitation gap ∆ n ∼ 1/ξ n inversely proportional to the correlation length ξ n . In particular, the singlefermion correlation C (1) probes the single-particle gap ∆ 1 for fermionic excitations, and the fermion-bilinear correlation C (2) probes the gap ∆ 2 for collective bosonic excitations. Various studies [105,118,136,145] have used the correlation function approach to demonstrate the gap opening (mass generation) in the SMG phase, as illustrated in Fig. 3(a,b). Ref. 136 observed that the fermion single-particle gap ∆ 1 is a bit larger than the bosonic fluctuation gap ∆ 2 in the SMG phase, which seems to be consistent with the fluctuating bilinear mass picture (where the bosonic excitations are typically softer than the fermionic excitations), as to be discussed in Sec. IV A. More explicitly, Ref. 100 investigated the fermion single-particle spectrum of (1+1)D chiral gauge theory coupled to a fluctuating Higgs field with finite spacetime correlation. By choosing the Higgs coupling strength and the Higgs field correlation length appropriately, Ref. 100 was able to show that the fermion excitation gap remains open in the SMG phase, despite of the absence of Higgs condensation.
Many numerical works have also confirmed that there is no fermion bilinear condensation in the SMG phase, such that the mass generation is not due to spontaneous symmetry breaking. To check this statement, one can introduce a small source field m that couples to the fermion bilinear termψΓψ of interest. By showing that in the large system size limit (L → ∞) the induced fermion bilinear expectation value vanishes as the source field m is turned off one can rule out the fermion bilinear condensation, as illustrated in Fig. 3(c). Such check should be performed for all fermion bilinear terms that transform as different representations under the symmetry. This approach was adopted in various studies [104,139,140,143,145,146].
Another approach is to measure the correlation function of fermion bilinear operators. If the correlation function decays exponentially, then the bilinear operator is not long-ranged ordered (not condensed). This approach is used in Ref. 105 and 136.

B. Nature of the SMG Transitions
Much effort of the numerical study has been focused on understanding the nature of the SMG transition. Two questions can be asked: • Is the SMG transition a direct transition (i.e. without any intermediate phases setting in)? • Is the SMG transition continuous (i.e. not firstorder)?
If the answers to both questions were yes, the SMG transition would be a quantum critical point which admits a field-theory description in the continuum limit. Such case will be highly interesting from the theoretical perspective, because the phases on both sides cannot be distinguished by a symmetry-breaking order parameter, hence a direct and continuous SMG transition is necessarily an exotic quantum phase transition beyond the Landau-Ginzburg-Wilson paradigm. The nature of the SMG transition depends on the spacetime dimension and the symmetry. Tab. III summarizes numerical results in the literature. SMG transitions in (1+1)D are most well understood. As four-fermion interactions are marginal perturbations for (1+1)D massless fermions, the RG equation for the interaction strength g with respect to the RG scale = ln Λ takes the general form of dg/d = αg 2 to the leading order of g. If the coefficient α > 0, g > 0 is marginally relevant, which flows towards infinity for any finite g. This happens to any SMG transition driven by four-fermion interactions in (1+1)D: the system immediately enters the symmetric gapped phase as long as the interaction is turned on (with the correct sign), such that the SMG critical point is at g c = 0. The RG equation predicts that the excitation gap should open as (see Fig. 3 [146][147][148] which has been verified in Ref. 136 and 142 with different models. However, SMG can also be driven by higher-order fermion interactions. For example, the SMG in (1+1)D U(1)-symmetric chiral fermion 3-4-5-0 model is driven by six-fermion interactions, which are irrelevant at the massless free-fermion fixed point. As shown in Ref. 105, the SMG transition in this case happens at a finite interaction strength, and the transition is shown to be in the BKT universality class.
Numerical studies have also found evidence for direct and continuous SMG transitions in higher dimensions, as listed in Tab. III. A direct transition can be established by ruling out any intermediate phase between the gapless fermion phase and the SMG phase. The most probable intermediate phase is the spontaneous symmetry breaking (SSB) phase in which the fermion bilinear mass (the Yukawa field) condenses. Such intermediate SSB phase is often observed in (3+1)D systems. [146][147][148] However, there are also examples showing that it is possible to shrink the intermediate SSB phase by frustrating the Yukawa field (i.e. introducing local couplings of the Yukawa field in conflict with its natural ordering tendency), which could foster a direct SMG transition. [140,143,144] Given a direct SMG transition, one can further check whether the transition is continuous or first-order (discontinuous). The numerical evidence for a continuous transition include: • Continuous gap opening across the SMG transition.
• Universal scaling behavior of physical quantities near the transition.
The excitation gap (or the inverse correlation length) can be extracted from correlation functions (as previously explained around Eq. (27)). Ref. 118 and 136 have demonstrated that the gap opens smoothly across the SMG transition, as depicted in Fig. 3(b), in support of a continuous transition (i.e. a quantum critical point). As the correlation length diverges, physical quantities should exhibit universal scaling behaviors near the quantum critical point. The quantity that is often studied in numerics is the uniform static susceptibility of the Yukawa field, which is defined as where the integration is over the spacetime and L d stands for the spacetime volume. The following universal behavior is observed in Ref. 137 and 138 where f is a universal function, and η and ν are critical exponents. Such that the susceptibility generally diverges as χ ∼ |g −g c | −ν(2−η) near the SMG critical point, as illustrated in Fig. 3(d). The exponent η can also be determined from the power-law fitting of the correlation function at the critical point, This approach is used in Ref. 139. Currently, the scaling analysis has only been performed for (2+1)D systems, where the simulation can achieve a linear system size up to L = 60 [138], which enables a rather reliable estimate of critical exponents. As a comparison, for (3+1)D systems, the simulation can only reach L = 16 [146]. Numerical studies [137][138][139] of the (2+1)D su(4) symmetric SMG found η ∼ 1.05 and ν ∼ 1.30, which is close to the (large-N f limit) exponents (η = ν = 1) of the (2+1)D Gross-Neveu-Yukawa universality class [151,152]. However, the SMG mechanism is physically distinct from the symmetry-breaking mass generation described by the Gross-Neveu-Yukawa model. Their similar critical exponents motivate the idea to view the SMG transition as a hidden Gross-Neveu-Yukawa transition of fermionic partons, as to be discussed in Sec. IV B.

A. Fluctuating Bilinear Mass Picture
The fluctuating bilinear mass (Yukawa field) provides an intuitive physical picture for SMG. It suggests that the SMG can be generally understood in two steps: starting with gapless Weyl/Majorana fermions, first condense a Yukawa field that couples to fermion bilinear mass terms to gap out the fermions at the price of breaking the protecting symmetry, then fluctuate the phase (or direction) of the Yukawa field in the spacetime to restore the symmetry while maintaining the local amplitude of the Yukawa field finite to keep the fermion gap open. The picture can be described by the Yukawa-Higgs model with a symmetry G [100,145] where ψ a are the Weyl/Majorana fermions (written as real spinor Grassmann fields withψ a = ψ a γ 0 ) and φ α are the Yukawa bosons (written as real scalar fields) [153]. Both the fermion ψ and the boson φ fields are in nontrivial representations of the protecting symmetry G, denoted as r G ψ and r G φ respectively. The Yukawa coupling will be loosely denoted as φψΓψ for generality, where the vertex tensor Γ α ab is set by the Clebsch-Gordan coefficients of the fusion channel r G ψ × A r G ψ → r G φ . Precise forms of the Yukawa coupling can be found in Eqs. (6,12,20,25), which vary from model to model. S Y [ψ, φ] describes the gapless fermion ψ coupled to the Yukawa field φ, and S H [φ] describes the dynamics of the Yukawa field with V H (φ) being some Higgs potential.
Integrating out the Yukawa field φ leads to a pure fermion model of ψ. Assuming S H [φ] = 1 2g φ 2 takes the Gaussian form at the ultraviolet (UV) level, the fermion model will be It is assumed that the model is free of G-anomaly, but the symmetry G is still restrictive enough to forbid all possible bilinear mass terms, paving ways for SMG. The fluctuating Yukawa field effectively mediates the fermion interaction g that is responsible to drive the SMG.
The appearance of the SMG in the Yukawa-Higgs model Eq. (33) requires some delicate design of S H [φ] to achieve the appropriate infrared (IR) dynamics. Assuming the Higgs potential V H (φ) pins the Yukawa field φ to a finite amplitude (e.g. |φ| = 1) throughout the spacetime (as in a non-linear σ-model), while allows its orientation to fluctuate in its internal flavor space, numerical evidences [100] show that the fluctuating Yukawa field is possible to gap out the fermions ψ without spontaneously breaking the G symmetry, if the Yukawa field fluctuation has a finite but large correlation length ξ 1 (s.t. φ(x)φ(0) ∼ e −|x|/ξ ). The fluctuation must be smooth enough so as not to create sharp domain walls that trap gapless domain-wall fermions and close the fermion gap. But the fluctuation must also not be too smooth to establish long-range order of the Yukawa field and break the symmetry. The SMG should be achieved by balancing these two factors.
An alternative way to argue for the SMG is the topological defect condensation approach [16,17,94,95,154]. This approach also starts by condensing the Yukawa field with a finite amplitude, but then disordering the Yukawa field orientation by condensing topological point defects (e.g. vortices in 2D space or monopoles in 3D space) of the Yukawa field. Although the point defects could trap fermion zero modes and leads to gapless fermion excitations in the spectrum, these fermion zero modes can be gapped out by local interactions (such as the FK interaction) uniformly applied throughout the system. A non-trivial check in this approach is to show that for anomaly-free fermion systems (that admits SMG), the point defect always trap 8n Majorana zero modes that can be trivialized by interaction. The advantage of the topological defect condensation approach is that it does not rely on a delicate tuning of the smoothness of the Yukawa field fluctuation.
Further justification of the fluctuating bilinear mass picture comes from a more explicit trial wave function construction for the SMG state in (0+1)D [35]. Let |Ψ[φ] be the quantum many-body ground state of fermions ψ on a background configuration of the Yukawa field φ. The fluctuating bilinear mass picture Eq. (33) suggest that the SMG state (the featureless gapped state in the strong interacting limit) should be described by In the (0+1)D example, the ground state of the Yukawa )|0000 using notations introduced in Sec. II A. Assuming S H [φ] restricts the φ vector uniformly on a sphere, the spherical average results in |Ψ SMG ∝ (1+c † 1 c † 2 c † 3 c † 4 )|0000 which indeed matches the ground state |4e of the FK interaction. The construction also applies to su(4) symmetric SMGs in all higher dimensions, where the on-site interaction is the same as the (0+1)D case. In the limit that the Yukawa coupling is strong and the Yukawa field φ fluctuates independently on each site, the trial wave function |Ψ SMG in Eq. (35) approaches to product state i |4e i -the exact ground state in the strong coupling (g → ∞) limit. Away from the strong coupling limit, |Ψ SMG might still serve as a variational approximation of the SMG state, but this variational ansatz is less well-controlled beyond (0 + 1)D.
Although the fluctuating bilinear mass picture provides a simple understanding for the SMG phase/state, it breaks down near the SMG transition, where the Yukawa field amplitude is no longer large and stable. A major theoretical challenge is to develop field theory descriptions for the continuous SMG transition, which will be the topic of the following subsection Sec. IV B.

B. Fermion Fractionalization Field Theory
The fermion fractionalization [35,141] provides an overarching theoretical framework to understand various continuous SMG transitions in different dimensions. It unifies different theoretical pictures of SMG in literature.
Its key proposal is that the physical fermion ψ ∼ η n ψ fractionalizes into bosonic η and fermionic ψ partons that deconfine at and only at the transition, where the bosonic parton η may appear n multiple times in the fractionalization scheme. Note that upright greek letters will be used to denote parton fields in the parton theory.
The fermionic parton ψ is generally put in the same gapless phase as the physical fermion ψ. The physical free fermion phase corresponds to the condensed phase of the bosonic parton η, and the SMG phase corresponds to the disordered (gapped) phase of the bosonic parton η. Accompanied with the fractionalization, an additional gauge field a must be introduced to bind the partons together and to remove the redundant degrees of freedom introduced by fractionalization. Let K be the emergent gauge group (typically non-Abelian), then the partons are generally charged under an enlarged symmetry-gauge group G × K (where G is the symmetry group of the physical fermion), as summarized by Tab. IV. Representation assignments must be consistent with ψ ∼ η n ψ, φ ∼ψΓψ, and φ ∼ψΓψ.
The choice of representations should satisfy the following general rules: • The bosonic parton η is always in the (anti)fundamental representation¯K of K, such that its condensation can fully Higgs the gauge group K.
• The representation must be assigned in consistent with the fermion fractionalization ψ ∼ η n ψ, such that the following fusion channel must exist • The K-gauge field a is always in the trivial representation 1 G of G and the adjoint representation K of K.
Given the representation (r G ψ , 1 K ) of the physical fermion ψ, there is still the freedom to choose the gauge group K and the representation (r G ψ , r K ψ ) of the fermionic parton ψ. Once the choice is made, the multiplicity n and the representation (r G η ,¯K) of the bosonic parton η are fixed by Eq. (36). There could be multiple parton theories that describe the same SMG transition. Within the present framework, different parton theories are specified by different choices of K and (r G ψ , r K ψ ). The parton theory for the SMG transition generally takes the form of [35,141] where ψ and η are treated as real fields (withψ = ψ γ 0 ) to avoid unnecessary complications of complex conjugation, and D µ = ∂ µ − ia µ denotes the covariant derivative (assuming the gauge connection a is automatically represented in the representation of the matter field that it couples to). The parton theory describes the bosonic η and the fermionic ψ partons coupled together by the gauge field a, and is expected to provide an effective description of the low-energy physics around the SMG critical point. The bosonic parton mass r is treated as the only driving parameter in the parton theory, and is responsible for tuning the SMG transition. r < 0 corresponds to the gapless phase (the free-fermion fixed point), r > 0 corresponds to the trivially gapped phase (the strong-coupling fixed point), and r = 0 corresponds to the critical point where the SMG transition happens. The states of different fields in different cases are summarized in Tab. V. When r < 0, the bosonic parton η condenses (i.e. η = 0). As η carries the anti-fundamental representation of the gauge group K, the gauge structure is fully Higgs down by the bosonic parton condensation. At the same time, the fermionic parton ψ becomes the physical fermion ψ ∼ η n ψ effectively (as the fields ψ and ψ only differ by a constant η n which simply serves as a quasi-particle weight). The parton theory then reduces to the physical fermion theory S[ψ] = d d xψiγ µ ∂ µ ψ in the gapless free-fermion phase.
When r > 0, the bosonic parton η is gapped and decouples from the gauge field at low-energy. Below the gap of η field, the parton theory is effectively a quantum chromodynamics (QCD) theory S F [ψ, a] involving the fermionic parton ψ coupled to the gauge field a. Then it relies on an assumption that the QCD dynamics leads to a featureless gapped ground state (a.k.a. a trivial ground state), such that both the fermionic partons and the gauge bosons are gapped without breaking the G symmetry or generating topological orders. Of course, one necessary condition for the QCD theory to trivialize is that the theory must have vanishing G × K anomaly. However, even if the QCD theory is anomaly free, it still depends on the dynamical details to achieve a trivially gapped ground state. The major theoretical effort to understand SMG lies in designing appropriate mechanism to trivialize the QCD theory. Within the fermion fractionalization framework, two mechanisms has been proposed in the literature: (i) parton-Higgs (weak coupling) [34,35,141] and (ii) s-confinement (strong coupling) [34,98]. In some cases, the QCD trivialization may be consistently achieved by both mechanisms.

Parton-Higgs mechanism
The parton-Higgs mechanism introduces a collection of scalar fields φ (parton-Higgs fields) that couple to the fermionic parton ψ via Yukawa interactions, such that the QCD theory S F [ψ, a] is extended to a QCD-Yukawa-Higgs theory where the parton-Higgs field φ (as fermionic parton bilinear mass) is in the representation (r G φ , r K φ ) of G×K. The representation (r G φ , r K φ ) must be obtained from an antisymmetric product of fermionic parton representations for consistency The Yukawa vertex tensor Γ α is fixed by this fusion channel. The QCD-Yukawa-Higgs theory in Eq. (38) is reminiscent of the Yukawa-Higgs theory in Eq. (33) by promoting the physical fermion ψ and physical Yukawa field φ to their parton counterparts ψ and φ. However, the key difference lies in the different G × K representations of ψ, φ comparing to ψ, φ, as listed in Tab. IV. Unlike the physical Yukawa field φ, which is solely charged under G in a non-trivial representation r G φ , such that the condensation of φ inevitably breaks the G symmetry; the parton-Higgs field φ is in a joint representation (r G φ , r K φ ) of G×K, which could admit a G-symmetric condensation of φ. The sufficient and necessary condition for the existence of a condensed configuration φ that preserves the G symmetry is that there exists a subgroup G×K of G×K, such that when G×K is broken to G×K , the representation of φ admits a trivial branching channel: If the condition is met, a Higgs potential V H (φ) can be constructed to drive the condensation of φ in the above trivial branching channel. Such a condensation of the parton-Higgs field φ will (i) gap out all fermionic partons ψ (that couple to it), (ii) Higgs down the gauge group from H to its subgroup K ⊆ K, and (iii) preserving the symmetry group G. Then the only freedom remaining in the theory is the pure gauge fluctuation of K . If K remains a non-Abelian group, its gauge coupling can flow strong into the confine phase and gap out the remaining gauge bosons. In this way, all freedoms in the QCD-Yukawa-Higgs theory are symmetrically gapped and the system ends up in the SMG phase. The physical fermion ψ ∼ ηψ will also be gapped, as its partons are gapped (both η and ψ are gapped).

s-confinement mechanism
The s-confinement mechanism ("s" for "smooth") refers to the strong coupling dynamics of the QCD theory that confines fermionic partons and gauge bosons all together without breaking the G symmetry. This is possible if the 't Hooft anomaly of G is matched up in the theory [3]. One motivation for the s-confinement theory comes from the observation that the four-fermion SMG interaction ψ 4 (for physical fermions ψ) can be broken up into "3+1" as (ψ 3 )ψ. If there is some strong-coupling mechanism that binds the first three fermions into a composite fermion ψ comp ∼ ψ 3 at low-energy, the interaction can then be viewed as a fermion bilinear mass likē ψ free ψ comp (whereψ free ∼ ψ stands for the last dangling fermion), which provides a free-fermion understanding for the SMG. This picture is alternative to the "2+2" fluctuation bilinear mass picture, where the interaction is broken up into (ψ 2 )(ψ 2 ) with each ψ 2 being a composite boson (i.e. the Yukawa boson φ). However, the Yukawa boson is still interacting, which complicates the analysis of its low-energy dynamics. In contrast, the composite fermion flows to a free fermion critical point, which enables simple argument for gapping. Nevertheless, this still relies on a strong-coupling mechanism for the composite fermion to form in the first place.
The fermion fractionalization plays an important role in understanding the formation of composite fermions in the s-confinement picture. After fractionalization, physical fermions become fermionic partons that couple to an emergent gauge field, then the emergent gauge force can be employed to confine the fermions into the desired composite. Not all fermions are fractionalized in the sconfinement theory. The physical fermions ψ = ψ frac ⊕ ψ free is first divided into two sectors, where ψ frac ∼ η n ψ are part of the physical fermions that will be fractionalized to bosonic η and fermionic ψ partons, and ψ free are the remaining physical fermions that do not fractionalize. To drive the SMG transition, the first step is still to gap the bosonic parton η, leaving the low-energy freedom in a QCD theory with the fermionic parton ψ cou-pled to an emergent (non-Abelian) gauge field a. Under gauge confinement, the fermionic partons ψ are bound together in the IR to form fermion composites, denoted as ψ comp ∼ (ψ) 3 , just like fundamental quarks forming baryons. The necessary condition for the confinement to happen is the existence of a trivializing fusion channel in the gauge sector such that the gauge charges of ψ 3 can be neutralized to produce the physical composite fermion ψ comp . Then the fermion bilinear massψ free ψ comp can be introduced to gap out all fermions together, realizing the SMG. Moreover, the s-confinement mechanism admits the supersymmetric extension [34], where addition of supersymmetry provides extra control for the strong-coupling confinement dynamics at low-energy [155].

C. Examples of Fermion Fractionalization
In the following, we will briefly exemplify the fermion fractionalization theory in two SMG models introduced in Sec. II.

(2+1)D Honeycomb Lattice Model
The SMG in the (2+1)D honeycomb lattice model can be understood using the parton-Higgs mechanism within the framework of fermion fractionalization [35]. The internal symmetry in consideration is where Z T F 4 is the essential protecting symmetry. The physical fermion ψ is in the 4 1 representation of G. The proposed fermion fractionalization scheme ψ ∼ ηψ introduces an emergent gauge group K = SU(4), (43) and put the fermionic parton ψ in the (1 1 , 4) representation of G × K. This choice fully specifies the parton theory. The physical and fractionalized fields in the parton construction are summarized in Tab. VI, such that all consistency conditions are met among their symmetry representations. The SMG can be achieved through the parton-Higgs mechanism by gapping out the bosonic parton η and condensing the parton Yukawa field φ ∼ψΓψ. This breaks G × K to G × K , where the new gauge group is In particular, because the parton Yukawa field φ charges 2 under both the Z T F 4 symmetry group and the Z K 4 center of the gauge group K = SU(4), condensing φ necessarily locks the Z T F 4 and Z K 4 generators together, breaking which is consistent with the Z K 2 center of residual gauge group K = Sp(2). As G × K is broken to G × K , the representation of φ splits into which admits a trivial representation (1 0 , 1) in the decomposition, satisfying the trivializing condition in Eq. (40) for the parton-Higgs mechanism to work.
The SMG in this model can be understood by fermion fractionalization, where only part of the fermions ψ frac = ψ ψ ⊕ ψ χ are fractionalized: ψ ψ ∼ η † ψ ψ and ψ χ ∼ η 2 ψ χ , while ψ free = ψ λ remains untouched. The fractionalization introduces an emergent gauge group Tab. VII summarizes the representation of fields satisfying all consistency conditions. To drive the chiral fermions ψ to the SMG phase, the first step is to gap out the bosonic parton η. Then the theory contains a free fermion theory of ψ λ and a QCD theory of the fermionic partons ψ ψ ⊕ ψ χ coupled to the K gauge field a. The trivialization of the QCD theory can be achieved by either the parton-Higgs mechanism or the s-confinement mechanism (they are equivalent).
In the parton-Higgs mechanism, two parton-Higgs fields φ = φ 1 ⊕ φ 2 are introduced and coupled to the fermions as  N + 4))G × SU(N + 4)K , as well as the Euclidean spacetime symmetry Spin(4), which are relevant to the SMG in the (3+1)D chiral fermion model.

(SU(N ) × SU(N + 4))G SU(N + 4)K Spin
which is closely reminiscent of the Yukawa-Higgs decomposition Eq. (25) The fermion representation has transmuted, morphing from a chiral representation of G × K into a vector-like representation of the surviving G × K . The parton-Higgs field condensation φ = 0 will provide mass for almost all fermions via the Yukawa coupling in Eq. fermions can be simultaneously removed from the lowenergy spectrum once the K gauge fluctuation flows strong and confines. Hence no low-energy freedom is left over, and the system ends up in the SMG phase.
In the s-confinement mechanism, a gauge-neutral combination of fermionic partons must first exist, as required by Eq. (41). As the fermionic partons ψ ψ ⊕ ψ χ carry the gauge charges¯⊕ respectively (see Tab. VII), the gauge neural three-fermion combination can be found in the fusion channel So the fermionic partons can be confined into a composite fermion through this confinement channel as which transforms as (( , 1), 1) under G × K, precisely conjugate to the representation (( , 1), 1) of the remaining physical fermion ψ λ . Therefore, the composite fermion ψ comp can pair up with the remaining fermion ψ λ to produce a G-symmetric mass term ψ λ iσ 2 ψ comp , that gaps out all fermions from the low-energy spectrum, leading to the SMG state.

Symmetry Breaking vs Symmetry Extension vs SMG
From the symmetry perspective, the fermion mass generation generally falls into two categories: (i) the symmetry breaking mechanism and (ii) the symmetry extension mechanism [23,24], see Fig. 4. At a high level, the SMG is essentially an symmetry extension mechanism in the absence of the 't Hooft anomaly. (Strictly speaking, the SMG mechanism does not enlarge the physical symmetry, but rather extends it to a symmetry-gauge structure. Nevertheless, we will still call it by symmetry extension in a general sense.) by condensing a Yukawa-Higgs field that couples to a fermion bilinear mass term. More precisely, starting from a symmetry group G (specifically here an internal symmetry, global or gauged), G is broken down to an appropriate subgroup G ⊆ G to induce quadratic mass term for fermions. Mathematically, it can be described by an injective homomorphism ι:

Mass generation mechanisms
Here are some explicit examples: • Bardeen-Cooper-Schrieffer (BCS) type Z 2 -gauged superconductor with a low energy Z 2 TQFT, we have G = Z 2 and G = U(1) electromagnetic gauge group.
, which breaks the unitary internal G = U(1) symmetry of a Weyl fermion down to a fermion parity G = Z F 2 symmetry. 2. Symmetry/Gauge extension: In contrast, the symmetry extension [24] or gauge enhancement [157] provides a fermion mass generation mechanism that preserves the symmetry. Symmetry extension construction of gapped phases first appears in [23] based on the gauge bundle descriptions, then Ref. 24 refines the idea to the lattice models, group-cohomology cocycle or continuum field theory descriptions. It extends the symmetry group G (that can include the spacetime-internal symmetry, global or gauged) to a larger groupG by enlarging the Hilbert space with additional/redundant degrees of freedom, in order to trivialize the 't Hooft anomaly or to lift any other symmetry obstruction towards a gapped phase with a mass generation. Mathematically, it can be described by a surjective homomorphism r:G which can be understood as part of the group extension in an exact sequence (see Eq. (56)) [24], or can be further generalized as the fibrations of their classifying spaces and higher classifying spaces [158][159][160]. The symmetry extension can induce the topological mass in the presence of anomaly, and the symmetric mass in the absence of anomaly [114].
• Topological Mass Generation (TMG) -If the fermion system has a non-trivial 't Hooft anomaly in G, the anomaly will post an obstruction toward trivial gapping, which already rules out SMG, leaving the possibility to TMG. A non-vanishing perturbative local anomaly disallows any symmetric gapped phase (even with topological order), also it can never be trivialized by a symmetry extension. So, in order to implement the symmetry extension construction, the non-vanishing Ganomaly must be a nonperturbative global anomaly in G.
For simplicity, the discussion below focuses on a limited special case of (55). If the global anomaly in G can become anomaly-free inG, by pulling the G group back to the extendedG group via a short exact sequence and if the normal subgroup K is a discrete finite group, then the fermion can acquire a topological mass upon gauging K, which gives rise to a discrete K gauge TQFT in appropriate spacetime dimensions and under appropriate criteria detailed in Sec. IV D 2. This is also called the group extension of the original quotient group G extended by a normal subgroup K to the total groupG. The topological mass refers to the energy gap of a TQFT matching the 't Hooft anomaly of G. Finding the group extensionG often requires the essential use of algebraic topology criteria, such as the Lydon-Hochschild-Serre spectral sequence method [24]. Ref. 24,57,114,159,161, and 162 provide several explicit lattice Hamiltonian or lattice path integral constructions.
• Symmetric Mass Generation (SMG) -In the case of SMG, the fermion system is already anomaly free, but the physical symmetry G is too restrictive to allow any symmetric fermion-bilinear mass. However, with the symmetry extension described by the following short exact sequence the fermion-bilinear mass (i.e. the parton-Higgs field φ) can be charged under both G and K, which can possibly be condensed, breaking only K to its subgroup K without breaking G, as long as the condensed fermionbilinear mass transform trivially under G×K ⊆ G×K. Moreover, the extended normal subgroup K (and K ) can be (and will always be) gauged, such that the actual physical symmetry will not be enlarged by the symmetry extension. This is simply a rephrasing of the parton-Higgs mechanism in the fermion fractionalization framework discussed in Sec. IV B.

l a t e x i t s h a 1 _ b a s e 6 4 = " i / S u u v q G d D 5 g J / c V V 6 e h E P j N W v w = " > A A A N P H i c b d f L b t t G F A Z g J r 2 l j t o 4 6 b J A Q V Q p 0 E V i S G 7 R Z F M g s X y / S p Y l 3 y g Y w 8 O R T G h 4 w Z x h Y p v h M 3 T d P k O f o d 2 2 j 9 B 9 d 0 U X 3 X R Z d C j J m B G P t R q d 7 8 x Y Q + q 3 h n 4 q Q l S N x h / 3 7 r / 3 / g c f f v T g 4 4 W H t U 8 + f b T 4 + E k f k 0 w C 7 0 E i E n n i M + Q i j H l P h U r w k 1 R y F v m C H / v j V u n H b 7 j E M I m P 1 H X K B x E b
x e E w B K Z 0 6 W L x G 0 + w e C S 4 l 1 6 G n p w O F b 9 S u Z t I t 7 j F L D X 8 f e N i s d 5 Y a k x e L h 0 0 Z 4 P 6 q 9 o P n Z 9 / + u + L 9 s X j h / 9

Q Q J Z x G M F g i G e N x u p G u R M q h A E L x a 8 D H n K Y M x G / F w P Y x Z x H O S T 3 R X u V 7 o S u E P 9 e Y Z J r N x J 1 Z 6 R C 6 Y / M V 5 H c + v k L E J d 8 0 k x Y u q S F M u V c b 6 q w v H N c w j m i + V k l S Q C C 3 e u r q 8 Y k z J 5 W 6 1 n 6 U h y P p 5 f x I + e B Z O / W N m G r c 5 G G c Z o r H M N 3 3 M B O u S t z y z r l B K D k o c a 0 H D G S o L 5 0 L l 0 w y U P r + L n g x f w t J F L E y D 1 f s u K 8 O c h n 9 8 + t N 9 8 V 8 x 1 j r i Y d 7 + p N d 3 Z j C 7 L G p G l 5 k N v r 1 J d v 2 y v 9 Q V D k X n l 9 Z J Q H V Q x D g 2 E V s x S M Q l V 7 x n p V 6 1 r Y n a j k t h 8 Y P i B z L e x S T S 1 N q Y a x 7 f p d p W M 0 6 G + i v m o q j e 2 3 t D P n d m e V W f j H G N G 5 y d 2 Q 0 L m z z F S N 7 u b d t y x v 5 Y o 8 m k L M J G 3 C o 8 L Q W 7 6 7 T c i w L y p s f J 9 w d c z 9 v 3 8 N f k 2 Y r R F a o t o y 2 q q 0 Z X q a Z X a O 6 b n S d 6 o b R D a q b R j e p b h n d o r p t d J v q j t E d q r t G d 6 n u G d 2 j u m 9 0 n + q B U Z K R A N t G 2 1 Q 7 R j t U D 0 e U u 0 a 7 V I 9 M n p E t W e U / M c I s G + 0 T / X Y 6 D H V E 6 M n V E + N n l I 9 M 3 p W h q G S O b h N Q 5 k n E g e E F Y t J H h B a h Z 1 G w q s W k 0 Q g r F l M I o G w b j H J B M K G x S Q U C J s W k 1 Q g b F l M Y o G w b T H J B c K O x S Q Y C L s W k 2 Q g 7 F l M o o G w b z H J B s K B x S Q c C G 2 L S T o Q O h a T e C A c W k z y g d C 1 m A Q E c h i k h C E n s U k I g h 9 i 0 l G E I t J i F B O L G Y p A T h 1 G I S E Q z i 8 / I j 8 a R 1 J y k X D K V y P K Y m O t K d Y 1 R X G 0 q S + R H 8 p B X 2 w 5 7 d q K q l 2 6
U v m L b / S 5 b H K a 0 q e m K K / r 3 7 p K Q 5 z I a N Y g + l w q f Z D y Z D m o N u q d z w 5 u f K j O P Z + P w r g 8 d M r w q t B z u O 6 Z v f N k O L p U g + r e 7 1 x A V 4 W Y X 8 U u 3 b 1 U 6 z a Y k o t 8 + k y g b w 8 T 3 E + u c n 3 I f 1 H k T z 1 9 X V I v G n M Z P 2 8 2 o m z y / m l R 7 l 8 / J D S r j w R 0 0 F 9 e a n 6 3 9 G 1 H P y 0 8 c 6 a v B 8 7 n z p f O 1 0 7 T e e G 8 c j a d t t N z w P n R + d X 5 z f m 9 9 k v t z 9 p f t b + n r f f v z e Z 8 5 s y 9 a v / + D x u l H 0 0 = < / l a t e x i t > h i or h i = 0 condensed < l a t e x i t s h a 1 _ b a s e 6 4 = " b w S Y M l 7 7 m M n 9 W 1 V C 5 w t s C 9 q w 2 X I = "

> A A A N Q H i c b d f L b t t G F A Z g J r 2 l j t o 6 6 b J A Q V Q p 0 E V i S E Z R Z 5 l Y v t 8 t W Z I v F I z h 4 U g m N L x 0 z j C x z f A Z u i / Q N + h T t N v 2 A f o G 3 R X Z d l F 0 K M m Y E Y + 1 G p 3 v z F h D 6 r e G f i p C V I 3 G X w 8 e f v D h R x 9 / 8 u j T h c e 1 z z 7 / Y v H J 0 x 4 m m Q T e h U Q k 8 t R n y E U Y 8 6 4 K l e C n q e Q s 8 g X v + + N W 6 f 0 3 X G K Y x C f q J u W D i I 3 i c B g C U 7 p 0 u b j i C R a P B P f S q 9 C T 0 6 H i 1 y p 3 E + k W d 5 i l F s f 8 R 7 d x u V h v L D U m L 5 c O m r N B / V X t p + N f f / 7 v 6 6 P L J 4 / f e 0 E C W c R j B Y I h X j Q b q R r k T K o Q B C 8 W v A x 5 y m D M R v x C D 2 M W c R z k k x 0 W 7 r e 6 E r h D / Z m G S a z c S d W e k Q u m P z X e R H P r 5 C x C X f N J M W L q i h T L l X G + q s L x 7 Q s I 5 o v l Z J U k A g t 3 r q 6 v G p M y e V u t Z + l I c j 6 e X 8 S P n g e T v 1 j Z u B q + H O R h n G a K x z D d 9 z A T r k r c 8 u 6 5 Q S g 5 K H G j B w x k q C + d C 1 d M M l D 6 H i / o W / M W k i h i c Z B 7 v m T F R X O Q z + 6 h W 2 + + K + Y 7 x l x N O t 7 V m + 7 s 5 h Z k j U n T 8 i C 3 1 6 k v 3 7 V X + o O g y L 3 y + s g o D 6 o Y h g b D K m Y p G I W q d o 1 1 q 9 a x s D N R y W 0 / N H x I 5 l r Y o Z p a m l I N Y 9 v 1 u 0 r H a M Z D f R X z U V V v b b 2 l n z u z P a v O x j n G j M 5 P 7 I a E z J 9 j p G 5 2 N + 2 4 Z 3 8 t U e T T F m A i b x U e F 4 L c 9 L t v R I B 5 U 2 P l + 4 K v Z + z 7 + W v y b c J V o 6 t U W 0 Z b V N e M r l F d N 7 p O d c P o B t V N o 5 t U t 4 x u U d 0 2 u k 1 1 x + g O 1 V 2 j u 1 T 3 j O 5 R 3 T e 6 T / X A 6 A H V Q 6 M k I w E e G T 2 i e m z 0 m G r b a J t q x 2 i H 6 o n R E 6 p d o + Q / R o A 9 o z 2 q f a N 9 q q d G T 6 m e G T 2 j e m 7 0 v A x D J X N w l 4 Y y T y Q O C K s W k z w g t A o 7 j Y T X L C a J Q F i 3 m E Q C Y c N i k g m E T Y t J K B C 2 L C a p Q N i 2 m M Q C Y c d i k g u E X Y t J M B D 2 L C b J Q N i 3 m E Q D 4 c B i k g 2 E Q 4 t J O B C O L C b p Q D i 2 m M Q D o W 0 x y Q d C x 2 I S E I Q T i 0 l C E L o W k 4 g g 9 C w m G U H o W 0 x C g n B q M U k J w p n F J C Y I 5 x a f k x + N E 6 k 5 S b l k K p H l M T H X l e o a o 7 j a V J b I j 2 S b V 9 v a n H Z t R 9 U u X a n 8 x T f 6 X D Y 5 T e l T U 5 T X 9 W 9 d p S F O Z D R r E D 0 u l T 5 I e b I c V B v 1 z m c H N z 5 U F 5 7 P R 2 F c H j p l e F 3 o O V z 3 z N 5 5 M h x d q U F 1 7 / c u o K t C z K 9 i l + 5 f q n U X T M l F P n 0 u 0 L e H C e 4 n 1 7 k + 5 K 8 U + T N P X 5 f U i 8 Z c x i + a j S i b v H 9 W l P v X D w n N 6 i M B H f S W l 5 o / L H 1 / r J 8 W n j v T 1 y P n K + c b 5 z u n 6 a w 4 r 5 w t 5 8 j p O u D 8 4 v z u / O H 8 W f u t 9 n f t n 9 r 7 a e v D B 7 M 5 X z p z r 9 q / / w N O 4 S D 4 < / l a t e x i t >
h i or h i 6 = 0

< l a t e x i t s h a 1 _ b a s e 6 4 = " b u d R Z E F W 0 / O 2 k x C J 0 S G y B p T k V o I = " > A A A N I 3 i c b d f L b t t G F A Z g J r 2 l j t o 4 7 a 4 F C q J K g S 4 S Q z K K J M v E 8 v 0 q W x f f K B j D w 5 F M a H j p n G F i m y H Q N + i + + 3 Z V o H 2 P 7 o p u u u h D 9 A W K D i U Z M + K x A A P k + c 6 M P R z 9 1 s h P R Y i q 0 f j 7 3 v 3 3 3 v / g w 4 8 e f L z w s P b J p 4 8 W H 3 / W x y S T w H u Q i E S e + A y 5 C G P e U 6 E S / C S V n E W + 4 M f + u F X 6 8 R s u M U z i r r p O + S B i o z g c h s C U L l 0 s f u E J F o 8 E 9 7 K U K + b J 6 U 3 M v 3 c b F 4 v 1 x l J j 8 n L p R X N 2 U X 9 V + / H w 1 5 / + + 6 p 9 8 f j h v 1 6 Q Q B b x W I F g i O f N R q o G O Z M q B M G L B S 9 D n j I Y s x E / 1 5 c x i z g O 8 s k i C v c b X Q n c Y S L 1 T 6 z c S d U e k Q u m + B V e R 3 P z 5 C x C X f N J M W L q k h T L m X G + q s L x z T M I 5 o v l Y J U k A g t 3 r s 6 v F J M y e V u t Z + l I c j 6 e n 8 S P n g a T 3 1 h Z u B q + H O R h n G a K x z B d 9 z A T r k r c c o P c I J Q c l L j W F w x k q B + d C 5 d M M l B 6 G x f 0 1 r y F J I p Y H O S e L 1 l x 3 h z k s z 1 0 6 8 1 3 x X z H m K t J x 7 t 6 0 5 1 t b k H m m D Q t D 3 J 7 n v r y b X u l P w i K 3 C u f j 4 z y o I p h a D C s Y p a C U a h q z 1 i v a h 0 L O x O V 3 P Y D w w d k r I U d q q m l K d U w t l 3 f V T p G M x 7 q p 5 i P q n p j 6 w 3 9 u z P b s + p o n G P M 6 P j E b k j I + D l G 6 m Z 1 0 4 4 7 1 t c S R T 5 t A S b y V u F x I c i m 3 7 4 j A s y b G i v v F 3 w 9 Y 9 / P X 5 N 3 E 6 4 Y X a H a M t q i u m p 0 l e q a 0 T W q 6 0 b X q W 4 Y 3 a C 6 a X S T 6 p b R L a r b R r e p 7 h j d o b p r d J f q n t E 9 q v t G 9 6 k e G C U Z C b B t t E 3 1 0 O g h 1 S O j R 1 Q 7 R j t U u 0 a 7 V H t G y X + M A P t G + 1 S P j R 5 T P T F 6 Q v X U 6 C n V M 6 N n Z R g q m Y P b N J R 5 I n F A W L G Y 5 A G h V d h p J L x q M U k E w p r F J B I I 6 x a T T C B s W E x C g b B p M U k F w p b F J B Y I 2 x a T X C D s W E y C g b B r M U k G w p 7 F J B o I + x a T b C A c W E z C g d C 2 m K Q D 4 d B i E g + E I 4 t J P h A 6 F p O A I H Q t J g l B 6 F l M I o L Q t 5 h k B O H Y Y h I S h B O L S U o Q T i 0 m M U E 4 s / i M f G h 0 p e Y k 5 Z K p R J b H x F x X q n O M 4 m p T W S I f k k e 8 2 n b E a d d W V O 3 S l c p v f K P P Z Z P T l D 4 1 R X l d f 9 Z V G u J E R r M G 0 e d S 6
Y O U J 8 u L a q N e + e z g x o f q 3 P P 5 K I z L Q 6 c M r w o 9 h u u e 2 Z 0 n w 9 G l G l T X f u c E u i r E / C x 2 6 e 6 p W r f B l F z k n j 5 h q 3 J 7 m O B + c p X r Q / 6 L I n / i 6 e e S e t G Y y / h Z s x F l k / s n R b l + / S W h W f 1 K Q C / 6 y 0 v N 5 0 v f H e p v C 0 + d 6 e u B 8 6 X z t f O t 0 3 R e O K + c T a f t 9 B x w f n B + c X 5 z f q / 9 X P u j 9 m f t r 2 n r / X u z M Z 8 7 c 6 / a P / 8 D g 9 k U V A = = < / l a t e x i t >

d D i P 5 8 X m Q X l A y k = " > A A A N C X i c b d f L b t t G F A Z g J r 2 l j t o k 7 b J A Q V Q p 2 k V i S E b R d N f E 8 v 1 u y Z J v F I L h 4 U g m N L x g z j C x z X D Z V b t t 3 y O 7 o N s + R Z + h 3 f Q N O p R l z G h O C B i g z n d m p C H 5 W 6 M w F z G q V u v v O 3 c / + P C j j z + 5 9 + n C / c Z n n z 9 4 + O i L A W a F B N 6 H T G T y J G T I R Z z y v o q V 4 C e 5 5 C w J B T 8 O J 5 3 a j 1 9 x i X G W H q m r n A 8 T N k 7 j U Q x M 6 V J 3 / b u X D 5 u t x d b 0 8 O l J e 3 b S / P n t P 7 9 8 / e 7 w 3 4 O X j + 7 / F 0 Q Z F A l P F Q i G e N 5 u 5 W p Y M q l i E L x a C A r k O Y M J G / N z f Z q y h O O w n H 7 U y v 9 W V y J / l E n 9 l y p / W r V H l I I p f o l X y d w 8 J U t Q 1 0 J S T J i 6 I M V 6 Z p y v q n h y / R S i + W I 9 W G W Z w M q f q / N L x a T M X r v 1 I h 9
L z i f z k 4 T J k 2 j 6 j s 7 C 1 e i n Y R m n e a F 4 C j f r H h X C V 5 l f 3 w Y / i i U H J a 7 0 C Q M Z 6 0 v n w w W T D J S + W Q t B y l 9 D l i Q s j c o g l K w 6 b w / L Q L B 0 L L j f b L + p 5 j s m X E 0 7 3 j T b f i C n X R W Z Y 9 q 0 N C z t e Z p L t + 1 O f x R V Z V B f H 5 m U k Y t x b D B 2 s c j B K L j a N 9 Z 3 r W d h b 6 q S 2 7 5 v e J + M t b B H N b c 0 p x q n t u t X T s d 4 x i N 9 F c u x q 9 e 2 X t P P X d h e u K N x j r G g 4 z O 7 I S P j 5 x i p m 9 X d d L x n f R 1 R l T c t w E T Z q Q I u B L n p t 0 9 E h G V b o / O 8 4 I s Z h 2 H 5 g j x N u G x 0 m W r H a I f q i t E V q q t G V 6 m u G V 2 j u m 5 0 n e q G 0 Q 2 q m 0 Y 3 q W 4 Z 3 a K 6 b X S b 6 o 7 R H a q 7 R n e p 7 h n d o 7 p v l G Q k w g O j B 1 Q P j R 5 S 7 R r t U u 0 Z 7 V E 9 M n p E t W + U / M e I c G B 0 Q P X Y 6 D H V E 6 M n V E + N n l I 9   • The full gauge group K must be large enough to counteract any non-trivial action of the symmetry group G on the parton-Higgs field φ, i.e. G must acts projectively on φ, such that the G symmetry can remain unbroken under the condensation of φ. • To achieve SMG in this framework, the deformation path must pass through the strong coupling regime where neither the bosonic parton η nor the parton-Higgs field φ has an expectation value, and the full gauge group K is unbroken. • After the parton-Higgs field φ condenses, the remaining unbroken gauge group K would better be either trivial or non-Abelian, such that either there is no residual K gauge fluctuation or the residual K gauge fluctuation can be confined automatically. Otherwise, if K is Abelian, it becomes possible that the SMG critical point will expand into a critical phase, described by an Abelian K -gauge theory.

More on Symmetry Extension Construction
The symmetry extension constructionG r − → G in (55) based on the pullback G-symmetry to the extended G-symmetry can be interpreted in different languages for different communities [24]: • For condensed matter, a nontrivial SPT state in the G symmetry cannot be deformed to a trivial tensor product state via a finite-depth of local unitary transformations without breaking the G-symmetry. But the successful symmetry extension means that we can find an appropriateG such that the SPT state in the extended G symmetric Hilbert space can be deformed to a trivial tensor product state via a finite-depth of local unitary transformations still preserving theG-symmetry.
• For quantum field theory or high-energy physics, the successful symmetry extension means that the 't Hooft anomaly in G-symmetry becomes anomaly-free inGsymmetry.
• For mathematics, the successful symmetry extension means that a nontrivial class of cocycle, cohomology, or cobordism of the G-symmetry becomes a trivial class in theG-symmetry. Suppose the nontrivial class of Gsymmetry cocycle, cohomology, or cobordism in the D dimensions is labeled by ω G D , then the trivialization means that its pullback (namely r * ω G D ) becomes aG-symmetry coboundary in the D dimensions (namely ωG D = δβG D−1 ) which splits to theG cochain (namely βG D−1 ) in the D −1 dimensions [24]. In summary, given the ω G D of the Gsymmetry, the successful symmetry extension requires to find a solution of both the extendedG and βG D−1 to satisfy Here the bulk dimension is D = d + 1, while the theory (on the boundary) with 't Hooft anomaly has its dimension D − 1 = d.
Below some examples of the symmetry extension construction of gapped phases in various dimensions based on the simplest short exact sequence in (56) with a finite group K, are considered. In particular, two issues should be addressed: (1) Given a theory with 't Hooft anomaly in G symmetry (which we shall also call it a dD anomalous boundary theory of a (d + 1)D bulk SPT), canG be found to trivialize the anomaly? If so, what is the minimalG?
(2) If theG is found, there are two implications of the construction based on (56) [24]: •G-symmetric extended gapped phase (as ã G-symmetric gapped boundary of the bulk G-SPT). In this case, all K,G and G are not dynamically gauged.
• G-symmetric gapped dynamical K-gauge theory with a 't Hooft anomaly in G. When K is dynamically gauged, in some cases, the G is preserved at IR; in other cases, the G becomes spontaneous symmetry breaking (SSB) at IR. One should also be careful to distinguish the two different kinds of dynamics.
Ref. 24 proves that "given any unitary or anti-unitary finite group G, for d ≥ 1, there always exists a finite K-extension to aG-symmetric extended gapped phase to trivialize any nonperturbative global anomaly in G (see also the proof given later in [158])." It is physically meaningful to consider gauging K for d ≥ 2. Once K is dynamically gauged, Ref. 24 finds that G-symmetric gapped dynamical K-gauge theory with a 't Hooft anomaly in G can only be obtained for d ≥ 3, but also finds that the G-symmetry is spontaneously broken for d = 2.
The above result is established when G is an ordinary 0-symmetry that acts on the 0D point operator. The symmetry extension construction can also be generalized to higher global symmetries [163] such as the 1-symmetry that acts on the 1D line operator (see a recent review [164]). Here are some examples (Table VIII for a summary) including ordinary and higher symmetries: 1. d = 1, a (0+1)D anomalous theory and a (1+1)D bulk: • Two Kitaev's chains: In Sec. II A, the discussion was limited to gapping (0+1)D eight Majorana modes or (1+1)D eight Kitaev's chains by the SMG preserving the Z T 2 × Z F 2 , which is free from the Z 8 global anomaly. Beyond SMG, even for (0+1)D two Majorana modes or (1+1)D two Kitaev's chains with 2 mod 8 class of Z 8 global anomaly, it is still possible to gap the whole system, leaving a single ground state without breaking G = Z T 2 × Z F 2 , but instead by extending the symmetry to a dihedral group of order 8 asG = D F,T [114]. Namely, when T 2 = +1 is extended to a time-reversal symmetry fictionalization T 2 = −1 and T 4 = +1, theG = D F,T 8 symmetric interactions can lift up the degenerate Majorana zero modes. Crucially the time-reversal Z T 2 generator T does not commute with the fermion parity Z F 2 generator (−1) F , but T (−1) F T −1 = −(−1) F . This means that T switches a bosonic sector |B and a fermionic sector |F in the Hilbert space. This is called supersymmetry extension that trivializes this (0+1)D fermionic anomaly and also trivializes the (1+1)D fermionic SPT state [114]. Another related property is that the preserved symmetry demands the anomalous boundary theory of ±2 mod 8 Majorana zero modes must be N = 2 supersymmetric quantum mechanics with two supercharges [165]. • Related studies on the fractionalized symmetries on the boundary of the layers of Kitaev chains can also be found in [166][167][168][169][170]. In particular, Ref. 167 and 169 studies the pure Z 2 class global gravitational anomaly on the boundary of a single Kitaev chain (which is an invertible fermionic topological order beyond the SPT, known as the mathematical Arf invariant). The pure gravitational anomaly cannot be trivialized by any symmetry extension, but may be "trivialized" by coupling to a gravitational theory [167].
bdry of 2 Kitaev chains: bdry of 2 Kitaev chains: bdry of (2+1)D CZX model: bdry of (3+1)D kP(B2) higher-SPT: (2+1)D surface state SO × Z 2, [1] even k: [1] even k: The third column shows whether theG-symmetry-extended gapped phase exists. The fourth column shows how the boundary G-anomaly or bulk G-SPT classification is reduced inG. The last column shows the dynamics after gauging K: Even if G-symmetry-extended gapped phase exists, the K-gauged dynamics can induce either a G spontaneous symmetry-breaking (SSB) or a G-symmetric K gauge theory.
Here a (0+1)D theory has no parity P but only at most time-reversal T , and Majorana fermion has no charge conjugation C symmetry; so only the T fractionalization is found. In higher dimensions, the common theme along the direction of this phenomenon is the C-P -T fractionalization [171].
2. d = 2, a (1+1)D anomalous theory and a (2+1)D bulk: • The (1+1)D edge of a (2+1)D CZX model as a Z 2 -SPT state is known to allow a symmetric gapless or a symmetry-breaking gapped boundary [172]. However, the G = Z 2 can be extended to give aG = Z 4symmetry-extended gapped boundary. Unfortunately, gauging the normal subgroup K = Z 2 results in a (1+1)D discrete K-gauge theory with G = Z 2 spontaneous symmetry breaking, which is consistent with the standard lore that there is no (1+1)D non-invertible intrinsic topological order, at least in the bosonic systems.
Other applications of the symmetry-extension construction on the (1+1)D gauge theories and orbifolds can be found in a recent survey [173].
• The (2+1)D surface state of (3+1)D kP(B 2 ) higher-SPT state with k ∈ Z 4 : This higher-SPT state is protected by a 1-form symmetry (denoted as G = SO × Z 2, [1] ) which couples to a 2-form background field B 2 . The P(B 2 ) means the Pontryagin square of B 2 . The symmetry-extension construction can obtain a gapped phase for an even k via extending toG = Spin × Z 2, [1] (although gauging K = Z 2 results in G SSB), but the symmetry-extension trivialization is proven to not exist for an odd k [160]. Later Ref. 176 proves a no-go theorem that the symmetry-preserving TQFT also does not exist for an odd k. This means the (2+1)D surface state must be either symmetric gapless or symmetry-breaking for an odd k.

d = 4, a (3+1)D anomalous theory and a (4+1)D bulk:
• The (3+1)D boundary of a (4+1)D Z 2 -SPT state allows a symmetric gapless, symmetry-breaking, or symmetric gapped surface topological order boundary. The G = Z 2 can be extended to give aG = Z 4 -symmetryextended gapped boundary [24]. Gauging the normal subgroup K = Z 2 results in a (3+1)D discrete K-gauge theory with both electric and magnetic gauge charges carries a fractional G charge.
• The (3+1)D boundary of a (4+1)D w 3 (T M )B e 2 higher-SPT state is protected by a 1-form electric symmetry (denoted as G = Z e 2, [1] , coupled to a 2-form B e 2 field), while the w j (T M ) is the jth Stiefel-Whitney class of the tangent bundle T M of spacetime manifold M . The corresponding anomaly occurs as a part of the anomaly of (3+1)D SU(2) Yang-Mills theory coupled to two Weyl fermions in the adjoint representation of SU(2) (below called this theory as the adjoint QCD 4 , see [177][178][179] and [160]). The G = Z e 2, [1] can be extended to give aG = Z F 2 × Z e 2,[1] -symmetry-extended gapped boundary [160]. Gauging the normal subgroup K = Z 2 results in a (3+1)D discrete K-gauge theory such that its electric gauge charge has fermionic statistics.
• The (3+1)D boundary of a (4+1)D kAP(B e 2 ) higher-SPT state: Again this is part of the anomaly of the adjoint QCD 4 [178]. The higher-SPT is protected by a Z 8 -axial symmetry (coupled to a 1-form A field, with its fourth power of the symmetry generator equals to (−1) F ) and a 1-form electric symmetry Z e 2, [1] (coupled to a 2-form B e 2 field), which can be denoted as a spacetime-internal symmetry G = Spin × Z F 2 Z 8 × Z e 2, [1] . There is a k ∈ Z 4 class. The even k class can be trivialized by a Z 2 -extension toG = Spin × Z 8 × Z e 2, [1] . The odd k class cannot be trivialized by any symmetry extension [160]. Later Ref. 176 proves a no-go theorem that the symmetry-preserving TQFT also does not exist. The above results [160,176] turn out to rule out certain UV-IR duality proposal of the adjoint QCD 4 hypothesized in [179].
• The (3+1)D boundary of a (4+1)D 1 2w 1 (T M )P(B e 2 ) higher-SPT state: The boundary turns out to associate with the anomaly of the (3+1)D SU(2) Yang-Mills gauge theory with a θ = π F ∧ F topological term (denoted as SU(2) θ=π YM) [180][181][182], while thẽ w 1 (T M ) is the twisted first Stiefel-Whitney class of T M such that 1 2w 1 (T M )P(B e 2 ) is a mod 2 class. The SU(2) θ=π YM kinematically at UV has time-reversal Z T 2 and 1-form Z e 2, [1] symmetries. The G = Z T 2 × Z e 2, [1] can be extended toG = Z T 2 × Z e 4, [1] to trivialize the anomaly, thus theG-symmetric extended gapped phase can be constructed [181,182]. However, upon gauging K = Z 2, [1] , this induces the G spontaneous symmetry breaking (SSB) [182]. Later Ref. 176 proves a no-go theorem that the symmetry-preserving TQFT also does not exist. The above results together demand that the IR fate of SU(2) θ=π YM must be either symmetric gapless or symmetry-breaking only. = Z 16 , see [183,184] and [57]. It turns out that k = (−N f ) of such SPT state captures a global anomaly of (3+1)D 15N f -Weylfermion Standard Model (SM) where N f is the number of families of quarks and leptons [183,185]. If k is odd, Ref. 186 proves an obstruction, so the symmetry-gapped TQFT is not possible to saturate this odd k anomaly. If k is even, Ref. 184, 187, and 188 show that two layers of symmetry extensions can construct theG = Spin × Z 8symmetry extended gapped phase: the first layer 1 → Z 2 → Spin × Z 4 → Spin × Z F 2 Z 4,X → 1 and the second layer 1 → Z 2 → Spin × Z 8 → Spin × Z 4 → 1. These constructions may have applications to beyond the SM physics [187][188][189].

A. Green's Function Zeros
Given different mass generation mechanisms discussed above, a key question is how to diagnose the SMG mechanism: if a mass gap is observed in an interacting fermionic system, how do we know the mass gap is opened up through the SMG mechanism, rather than a more con-ventional symmetry-breaking mechanism, say by condensing a fermion-bilinear mass. Over the years it was gradually realized that one of the characteristic features of the SMG phase is that, the fermion Green's function G(iω) = − ψ(iω)ψ † (iω) in the Matsubara frequency space should have a zero at ω → 0, i.e. det G(iω → 0) = 0. This applies to the SMG phase either in the bulk or on the boundary of a fermionic system.
There are various ways to argue the necessary existence of this Green's function zero in a system with SMG. For the Green's function defined in the bulk, one general argument for the existence of zero is based on the topological number associated with the free fermionic SPT state [120,190,191]. For a fermionic system with translation symmetry, the topological number can be defined in the Matsubara frequency and momentum space using the Euclidean spacetime full fermion Green's function [192][193][194][195][196]: (59) where G is the matrix of the full fermion Green's function, and B is certain matrix in the flavor space. For a Chern insulator at even spatial dimension, B is an identity matrix, and n is equivalent to the Chern number. One can prove mathematically that the number n must be a quantized integer, and it can only change discontinuously when the Green's function has singularity. Schematic phase diagram between a trivial phase and a (fake) fermionic SPT phase that can be trivialized by interaction. Massless bulk fermions along the critical line in the weak-coupling limit undergoes the SMG as the interaction g exceeds a critical value gc. Although the two phases (technically one phase) are smoothly connected through the strong-coupling regime, the topological number n must still jump across the missing "phase boundary", along which the Green's function must have a zero.
The singularity of G happens at two types of "transitions". The first type of transition is a physical transition where det G −1 (iω = 0) = 0 vanishing to zero, i.e. the fermions become gapless. This corresponds to the quantum critical line between the fermionic SPT and trivial phases in the weak-coupling regime, as shown in Fig. 6. However, note that in the definition of the topological number n in Eq. (59), G −1 and G appear on an equal footing, hence n can also change when det G(iω = 0) = 0, i.e. when the Green's function has a zero. Hence when the free fermion SPT phase is trivialized by interactions, although there is no unavoidable phase transition between the SPT and the trivial phase, the topological number n still has to change discontinuously somewhere in the phase diagram, as shown in Fig. 6. Since there is no real physical transition, the number n has to change through the zero of the full interacting fermion's Green's function. This must hold throughout the SMG phase (line) regardless of the interaction strength g (as long as g > g c ).
Computing the full fermion's Green's function for a given interacting Hamiltonian is generally a challenging task. To argue that the Green's function zero is a general feature of the SMG state independent of microscopic details, one relies on the topological defect condensation construction for the SMG state, introduced in Sec. IV A. Starting with the SMG of eight Majorana zero modes in (0+1)D, the fermion Green's function can be explicitly computed [136,197]: where m is an effective mass proportional to the strength of the fermion interaction. Obviously, in (0+1)D, G ab (iω) approaches zero as iω → 0. Now one can construct a higher dimensional SMG state by decorating every topological point defect of the Yukawa field with eight interacting Majorana fermions in its SMG state, and then proliferating the point defects to put the Yukawa field in a disordered phase, following the strategy of the decorated domain wall construction [198]. The fermion Green's function can be evaluated in the spacetime by patching the (0+1)D Green's functions along the world line of the topological defect, and then path integrating all possible world line configurations. Following this approach, Ref. 199 was able to show that the fermion Green's function takes the general form of which universally exhibits a zero G(k µ = 0) = 0 in zero momentum-energy limit. Before the non-perturbative prove by Ref. 199 [35], where the physical fermions ψ fractionalize into deconfined bosonic η and fermionic ψ partons at and only at the critical point.
The concept of DQCP [200][201][202] was originally introduced to describe the direct continuous transition between the antiferromagnetic phase and the valence bond solid phase in (2+1)D quantum spin models. The two phases break distinct symmetries (spin-rotation and lattice-rotation symmetries respectively) and cannot be connected by a single quantum critical point without finetuning in the Landau-Ginzburg-Wilson paradigm. The DQCP provides an explanation for this exotic quantum critical point by fractionalizing the physical spin into deconfined spinons (partons) at and only at the critical point. The continuous SMG transition is also an exotic quantum critical point that requires the DQCP description. However, Unlike the conventional (bosonic) DQCP that fractionalizes bosonic degrees of freedoms, the continuous SMG transition is a fermionic version of DQCP, as it fractionalizes fermionic degrees of freedoms.
The SMG critical point can be viewed as the intersection among four phases: The four phases corresponds to the four different states of the bosonic parton η and the parton Higgs φ fields, as summarized in Fig. 5. Fig. 7(a) shows a schematic phase diagram. The DQCP emerges at the origin of the phase diagram, where both η and φ are critical. It seems that a direct transition between the massless fermion phase to the SMG phase requires fine tuning through the DQCP. However, if the fermionic parton QCD theory is dynamically unstable, the QCD phase will shrink to a single transition line between the massless fermion and the SMG phases, such that a direct continuous SMG transition can persist without fine tuning. Strictly speaking, the DQCP picture for the SMG transition is only valid and necessary in (2+1)D or higher dimensions. The SMG in lower dimensions (i.e. (0+1)D and (1+1)D) are exceptional and can be understood without involving fermion fractionalization or DQCP.
• In the (0+1)D Fidkowski-Kitaev model, the fourfermion interaction will immediately opens the gap, which can be understood by solving the quantum mechanical problem in Eq. (1) exactly. There is no notion of phase transition and quantum criticality in (0+1)D, not to mention DQCP.
• In (1+1)D, take the 3-4-5-0 model for example, the SMG can be understood within the Luttinger liquid framework as a BKT transition. The sixfermion SMG interaction in Eq. (11) has a bare scaling dimension ∆ int | g=0 = 5 in the free fermion limit, which is perturbatively irrelevant. However, with a non-perturbative (finite) interaction  marginal marginal marginal marginal marginal marginal marginal marginal marginal  marginal  marginal  marginal  marginal  marginal  marginal  marginal  marginal  marginal  marginal  marginal  marginal  marginal  marginal  marginal marginal marginal marginal marginal marginal marginal marginal marginal  strength, under the RG flow, the Luttinger parameters (as exact marginal parameters) will be modified by the interaction, leading to the decrease of the scaling dimension of the interaction term. When the scaling dimension drops below 2 (which is the spacetime dimension), the SMG interaction will become relevant, driving the system into the SMG phase. [34,105] So the SMG transition is triggered right at ∆ int | g=gc = 2, see Fig. 7(b).
• In higher dimensions ((2+1)D and above), interactions are always perturbatively irrelevant at the free fermion fixed point, such that an infinitesimal interaction g will not immediately drive the SMG transition. Therefore, the transition generally requires a finite critical interaction strength g c . The critical point (g = g c ) is expected to be an unstable fixed point under RG, which either flows to the free fermion fixed point (g = 0) or the SMG (gapped phase) fixed point (g = ∞), as illustrated in Fig. 7(c).
One implication of fermion fractionalization at the DQCP is that the physical fermion scaling dimension will generally be larger at the SMG transition compared to the free fermion fixed point [35]. Massless free fermion (Dirac/Weyl/Majorana) in d-dimensional spacetime has the scaling dimension ∆ ψ | g=0 = (d − 1)/2. At the SMG transition point (g = g c ), the fermion fractionalizes to the critical bosonic parton η and the gapless fermionic parton ψ. In the large-N limit, the scaling dimensions of partons are ∆ η = d/2 − 1 and ∆ ψ = (d − 1)/2, such that the leading order estimation of the physical fermion scaling dimension becomes ∆ ψ | g=gc = ∆ η + ∆ ψ = d − 3/2. For d > 2, the large-N estimation implies i.e. the fermion scaling dimension will jump to a higher value right at the SMG transition, as shown in Fig. 7(c). For d = 2, the above naive dimension counting seems to indicate ∆ ψ | g=gc = ∆ ψ | g=0 (which is not correct, unless g c = 0). More careful Luttinger liquid RG analysis [105] shows that the fermion scaling dimension increases continuously with the interaction until the transition happens, as shown in Fig. 7(b). Thus, the statement ∆ ψ | g=gc > ∆ ψ | g=0 still holds for the (1+1)D model.

C. Symmetric Mass Generation of the Standard Model
One important motivation driving the study of SMG is to seek for the lattice regularization of chiral gauge theories, in particular the Standard Model of particle physics. The goal is to gap the fermion doublers (mirror fermions) without affecting the original chiral fermions [36,41,93,95,97]. Possible routes to gap the Standard Model with either 15N f or 16N f Weyl fermions via the SMG are reviewed in the following, where N f stands for the family (or generation) number and can be taken to be N f = 3.

Symmetry Extension of the 15N
The (3+1)D Standard Model (SM) has a chiral internal symmetry group of the Lie algebra su(3) × su(2) × u(1) Y , which could correspond to either of the four versions of Lie groups for q = 1, 2, 3, 6: which are all compatible with all known particle representation data to date. The following discussion is applicable to any of the four versions of SM for q = 1, 2, 3, 6. For any specific version, the SM contains N f = 3 families (generations) of matter fermions. In each family, there can be either 15 or 16 left-handed Weyl fermions. The corresponding SM phases will be denoted as 15N f -SM and 16N f -SM respectively. The first 15 Weyl fermions transform in the following representations of su(3) × su(2) × u(1) Ȳ d R ⊕ l L ⊕ q L ⊕ū R ⊕ē R = (3, 1) 2 ⊕ (1, 2) −3 ⊕ (3, 2) 1 ⊕ (3, 1) −4 ⊕ (1, 1) 6 (64) in each family. Here u R and d R are up and down types of right-handed quarks. The q L is the su(2) doublet of up and down types of left-handed quarks. The e R is the right-handed electron. The l L is the su(2) doublet of neutrino and electron types of left-handed leptons. The 16th Weyl fermion corresponds to the sterile neutrinō ν R = (1, 1) 0 , which can be appended to Eq. (64) in any of the families. For simplicity, all SM fermions are assumed to be massless, without including the electroweak Higgs and its symmetry breaking. To apply SMG to 15N f -SM or 16N f -SM, one shall confirm that the two necessary conditions in Sec. II E are satisfied. For the second condition in Sec. II E, it is true in the SM. For the first condition, several recent works had checked the cobordism group classification of anomalies in the SM [183,185,187,203]. Given the spacetime-internal symmetry G = Spin × G SMq , the anomaly index ν ∈ TP 5 (G) belongs to the cobordism group TP 5 (Spin × G SMq ) = Z 5 × Z 2 , q = 1, 3, Z 5 , q = 2, 6.
As checked in Ref. 185,187, and 203, The perturbative local anomalies (Z classes) and nonperturbative global anomalies (Z n classes) all vanish (ν = 0) for both the 15N f -SM and 16N f -SM, such that the SMG is possible in either cases. However, if an additional continuous baryon minus lepton symmetry U(1) B−L is to be preserved, the spacetime-internal symmetry is enlarged to G = Spin× Z F 2 U(1) B−L ×G SMq = Spin c ×G SMq , then the anomaly index ν ∈ TP 5 (G) belongs to a different cobordism group TP 5 (Spin c × G SMq ) = Z 11 , q = 1, 2, 3, 6, which only vanishes (ν = 0) for 16N f Weyl fermions [187,203]. Therefore, the SMG preserving an additional U(1) B−L only works for the 16N f -SM. Razamat  Hereafter for the left-handed and right-handed notations, we always use the italic font L and R to denote that of spacetime symmetry (Spin group), while use the text font L and R for that of internal symmetry.
SMG phase (1) Fig. 8. The transition from the 27N f -LR phase to the 16N f -SM phase can be driven by condensing a scalar Higgs field h R = (1, 1, 2) 0,−3 that couples to the Weyl fermions by the following Higgs term (70) where R is an anti-symmetric tensor in the SU(2) R doublet subspace. The L and tensor for the SU(2) L and Lorentz su(2) subspace are omitted. All Lagrangian terms become a scalar in a trivial singlet representation of both spacetime and internal symmetries. For each family, the Higgs condensation h R = 0 (only its upper SU(2) R doublet component is nonzero and condenses) lifts 11 Weyl fermions with a mass gap (the lower half of both SU(2) R doublets ofd R and l L get mass, while the dimension-5 term (h † Rν R )(h † Rν R ) gives Majorana mass to only the upper half of the SU(2) R doubletν R ), leaves 16 Weyl fermions at low energy, and breaks the symmetry from G LRq,p down to G SMq × U(1) B−L .
To further lift the sterile neutrino, an additional scalar Higgs field h R = (1, 1, 1) 0,6 should be introduced (with U(1) B−L charge −1), such that the Higgs term becomes In each family, the Higgs condensation h R , h R = 0 leaves 15 Weyl fermions at low energy (the h R ( RνR ν R ) gives a Dirac mass to both upper and lower components of the SU(2) R doubletν R ), breaks the symmetry from G LRq,p down to G SMq with no U(1) B−L , and drives the transition from the 27N f -LR phase to the 15N f -SM phase.
Embedding the 15N f -SM into the 27N f -LR amounts to the following symmetry extension and embedding the 16N f -SM into the 27N f -LR amounts to the following symmetry extension where SU(2) R = SU(2)R×U(1)R Z gcd(p,2) ×U(1) B−L .

Symmetric Mass Generation in a 27N f -Weyl-Fermion Left-Right Model
The SMG in the 27N f -LR model can be achieved by both the parton-Higgs mechanism in Sec. IV B 1 and the s-confinement mechanism in Sec. IV B 2.
1. The s-confinement mechanism: According to Razamat and Tong [98], one first supersymmetrizes the fermions in Eq. (67) to their corresponding N = 1 chiral multiplets asd R ⊕ l L ⊕ q L ⊕ū R ⊕ē R ⊕ν R ⊕ d L ⊕l R → D ⊕ L ⊕ Q ⊕ U ⊕ E ⊕ N ⊕ D ⊕ L . Then gauge the SU(2) R symmetry by turning on a dynamical SU(2) R gauge field that couples to the su(2) R doublet: D, L, N. A dangerously irrelevant superpotential W UV at UV W UV = LLE + DDU + LDQ + LNL + DND + h.c. (74) becomes a G SMq -symmetric relevant deformation that pairs the 15 mesons (Ẽ,Ũ,Q,L,D) formed by D, L, N via s-confinement) with the remained 15 superfields (namely Q, U, E, D , L ) in a quadratic manner W IR at IR, consequently gapping out all fields as SMG: W IR =ẼE +ŨU +QQ +LL +DD + h.c. (75) Here the su(3) color and su(2) L flavor indices are suppressed, with the understanding that they should be contracted properly to make the Lagrangian a singlet. As W IR flows strong, all fermions are gapped from lowenergy, resulting in the SMG phase. When there are multiple families, independent SU(2) R gauge fields are introduced in each family, such that the total gauge group is SU(2) R1 × SU(2) R2 × · · · × SU(2) R N f . This guarantees that the s-confinement can induce the fully gapped SMG phase in each family independently.
2. The parton-Higgs mechanism: Tong [34] shows that N f families of 27 Weyl-fermion model can be fully gapped by preserving not only the SM internal symmetry group G SMq for q = 1, 2, 3, 6, but also an additional continuous baryon minus lepton symmetry U(1) B−L . The parton-Higgs mechanism introduces the scalar Higgs fields φ = (1, 2, 2) −3,3 . Ref. 34 suggests to fully gap (the 27 Weyl fermions per family) to achieve the SMG by adding φ(d R q L +ν Rl R +ē R l L ) + φ 2 d LūR + h.c. (76) when the condensation φ = 0 occurs. The internal symmetry breaking pattern of this SMG deformation is which leaves a continuous baryon minus lepton number symmetry preserved.

VI. SUMMARY AND OUTLOOK
We have reviewed the topic of symmetric mass generation (SMG) which has attracted considerable ongoing interests from both the condensed matter and high energy physics communities in recent years. The SMG transition is beyond the classic Landau-Ginzburg-Wilson paradigm, which potentially involves fractionalization of physical fermions at and only at the critical point. We have reviewed various aspects related to SMG, including interacting topological insulators (TI), topological superconductors (TSC) [204,205], SPT states [206,207], anomalies, lattice regularization of chiral gauge theories, and the current status of numerical efforts from both the lattice gauge theory and condensed matter communities. We have also discussed theoretical understanding related to SMG, such as the s-confinement mechanism, the connection between SMG and DQCP, etc.
Various numerical works have suggested that the SMG could be a continuous transition, and the critical exponents have been measured in some of the numerical works. These works pose a challenge to further analytical understanding of SMG. Unlike the standard Higgs-Yukawa types of theories where a large-N or small-expansion can be applied, theoretically we do not yet have a controlled theory where we can perform a reliable analytical calculation for the critical exponents for SMG, and compare with the numerical simulation. If the SMG is indeed a continuous transition and corresponds to a certain type of conformal field theories, then the conformal Bootstrap method can also provide very helpful insights into the nature of the transition. All these require further efforts from different disciplines of theoretical physics.
The problem of the nonperturbative regularization of chiral fermion and chiral gauge theories also potentially bring the high-energy lattice, mathematical physics, quantum information, numerical simulation, and condensed matter communities to work closer together. What else can stimulate different communities to work together other than tackling a profound mysterious problem? Looking back at science history, the phenomena of the anomaly-inflow [5,6] and domain wall fermion [76,78] (as the precursors of the TI/TSC/SPT states) had an old tradition rooted in the high-energy theory and lattice community in the mid 1980s and mid 1990s. Closely related phenomena like integer and fractional quantum Hall states are studied in the condensed matter community already in the early 1980s. Yet the acceleration came much later until the concrete materialization discovery of TI/TSC in 2005 (in both theories and experiments [204,205]) ignites the serious interest in the classification of the interacting many-body quantum systems of SPT states in the early 2010s [59-61, 206, 207]. This generates crossing-disincline interests in the classification of anomalies, cobordism classes, and the TQFTs in quantum matter, high-energy string theory, and mathematical physics [4,[41][42][43][44]48]. Look back and ponder this history, some curious minds might ask: Why were the discovery of SMG and classification of TI/TSC/SPT states not made in the high-energy lattice community even earlier, although the lattice gauge theory pioneers like Wilson and Kogut already had earlier achievements on the related lattice topics in the 1970s [208]? Perhaps the inputs of the condensed matter material experiments are crucial, also perhaps the crossing-boundary interdisciplinary ideas are much more welcomed in the recent condensed matter developments. As Preskill summarized: Combing more ideas together is better than each of the isolated ideas [209]. Hopefully combing the ideas of (1) domain wall fermions [76], (2) gapping the mirror fermions [36], and (3) the anomaly-free SMG [41,93,95,97], already gives the ample insights to fully solve the lattice definition of chiral fermion, chiral gauge theory, and chiral Standard Model problems.