Prepared for submission to JHEP
Introduction
Topological defect lines are by now a standard tool in two-dimensional conformal field theory: they encode the symmetries of a theory, including the non-invertible ones, they act on bulk and boundary fields, and they constrain renormalisation-group flows [1 , 2 , 3 , 4 , 5 ] . In a diagonal minimal model M ( p , p + 1 ) M(p,p+1) M ( p , p + 1 ) the topological defects are the Verlinde lines D ( r , s ) D_{(r,s)} D ( r , s ) , labelled like the primary fields, and they form a fusion category. Kormos, Runkel and Watts (KRW) [6 ] asked what happens when such a line is perturbed by one of its own relevant defect fields: the defect is no longer topological, and it flows, in the manner of a boundary RG flow, to a new conformal interface, which may be topological again or may be a genuinely transmitting, non-topological defect. Perturbed topological defects of minimal models have been studied from several angles—chiral perturbations, which preserve pure transmission while generally breaking conformal invariance away from fixed points [7 , 8 ] ; RG interfaces between the UV and IR theories of a bulk flow [9 , 10 ] ; the defect flows of the tricritical Ising model [11 ] —but the non-topological endpoints of the KRW flows are known only through perturbation theory in 1 / p 1/p 1/ p , truncated conformal space approach [6 , 12 , 13 ] and through level-truncated string field theory [14 ] . A first-principles, non-perturbative lattice determination of these interfaces has been missing.
In this paper, we study the flow of the ( 2 , 1 ) (2,1) ( 2 , 1 ) defect D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) driven by the non-chiral, spinful defect field ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) \phi_{(1,3)} \otimes\bar{\phi}_{(3,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) , at p = 4 p=4 p = 4 (tricritical Ising, c = 7 / 10 c=7/10 c = 7/10 ), at p = 5 p=5 p = 5 (M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) , c = 4 / 5 c=4/5 c = 4/5 ), and at continuous p p p in between. In the two-parameter space of KRW, D ( 1 , 2 ) D_{(1,2)} D ( 1 , 2 ) perturbed by λ l ϕ ( 1 , 3 ) ⊗ 1 + λ r 1 ⊗ ϕ ˉ ( 1 , 3 ) \lambda_l\phi_{(1,3)}\otimes1+\lambda_r 1\otimes\bar{\phi}_{(1,3)} λ l ϕ ( 1 , 3 ) ⊗ 1 + λ r 1 ⊗ ϕ ˉ ( 1 , 3 ) , the defect D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) is the endpoint of the purely chiral flow and ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) \phi_{(1,3)} \otimes\bar{\phi}_{(3,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) is the first-order image of the anti-chiral field along it; perturbing D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) by it with one sign leads to the new conformal defect C (the endpoint of the diagonal λ l = λ r > 0 \lambda_l=\lambda_r>0 λ l = λ r > 0 flow from D ( 1 , 2 ) D_{(1,2)} D ( 1 , 2 ) ), with the other sign to the identity defect. Our main results are the following. At p = 5 p=5 p = 5 the flow ends on the genuinely conformal interface C 5 C_5 C 5 , whose g g g -function, transmission, operator content and leading irrelevant coupling we determine; this is the first non-perturbative lattice determination of a KRW defect. At p = 4 p=4 p = 4 , where KRW’s truncated-conformal-space analysis could not locate the fixed point, we find that the flow ends on the invertible topological line η = D ( 3 , 1 ) \eta=D_{(3,1)} η = D ( 3 , 1 ) , an endpoint that had not been proposed before, so that the tricritical Ising model has no non-topological KRW defect at all. Following the family C p C_p C p at continuous p p p then shows how the two cases are connected: C p C_p C p moves smoothly from C 5 C_5 C 5 to η \eta η as p → 4 p\to4 p → 4 . Notice the η = D ( 1 , 4 ) \eta=D_{(1,4)} η = D ( 1 , 4 ) at p = 4 p=4 p = 4 . The D ( 1 , 4 ) D_{(1,4)} D ( 1 , 4 ) interface exists at any p p p . When p → 4 p\to4 p → 4 , the C p C_p C p interface merges with D ( 1 , 4 ) D_{(1,4)} D ( 1 , 4 ) because the relevant operator driving the RG flow between them becomes a null-state.
Two facts make a lattice construction possible, and they are the two themes of the paper. The first is a property of the flow itself: it preserves the ( r , 1 ) (r,1) ( r , 1 ) Verlinde lines of M ( p , p + 1 ) M(p,p+1) M ( p , p + 1 ) , which form a fusion category C p \mathcal{C}_p C p with the fusion rules of S U ( 2 ) p − 2 SU(2)_{p-2} S U ( 2 ) p − 2 [4 , 5 ] .1 { 1 , σ ′ = D ( 2 , 1 ) , η = D ( 3 , 1 ) } \{\mathbf{1},\sigma'=D_{(2,1)},\eta=D_{(3,1)}\} { 1 , σ ′ = D ( 2 , 1 ) , η = D ( 3 , 1 ) } , the Ising category, at p = 4 p=4 p = 4 ; Fib⊠ Z 2 \boxtimes\mathbb{Z}_2 ⊠ Z 2 with D ( 2 , 1 ) = w = τ F i b ⊗ η D_{(2,1)}=\mathbf{w}=\tau_{\mathrm{Fib}}\otimes\eta D ( 2 , 1 ) = w = τ Fib ⊗ η of quantum dimension φ = 1 + 5 2 \varphi= \frac{1+\sqrt{5}}{2} φ = 2 1 + 5 at p = 5 p = 5 p = 5 ). We make this statement precise in § 2 : a line x x x wrapping the defect has to be attached to it by a junction, and since x ⊗ D ( 2 , 1 ) x \otimes D_{(2,1)} x ⊗ D ( 2 , 1 ) has two fusion channels there are two inequivalent ways of doing so, the over- and under-crossing “lassos”. The deformation selects one of them—the lasso that reads the holomorphic label of the defect fields, under which ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) \phi_{(1,3)} \otimes\bar{\phi}_{(3,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) is neutral—and we explain, in the three-dimensional Reshetikhin–Turaev description of the theory on Σ × [ − 1 , 1 ] \Sigma\times[-1,1] Σ × [ − 1 , 1 ] , why an over-crossing reads the holomorphic label and an under-crossing the antiholomorphic one. This is a point that does not arise in KRW’s original setting, where the flow starts from D ( 1 , 2 ) D_{(1,2)} D ( 1 , 2 ) and the junction space is one-dimensional. The second fact is a property of the lattice: the dilute A p − 1 A_{p-1} A p − 1 RSOS chain [15 , 16 ] is precisely an anyon chain [17 , 18 ] of C p C_p C p , for every integer p p p considered here. Its height variables are S U ( 2 ) p − 2 SU(2)_{p-2} S U ( 2 ) p − 2 labels, its domain walls are w w w anyons2 , and its Hamiltonian is built from the fusion data. A ( r , 1 ) (r,1) ( r , 1 ) defect line is therefore just one more anyon of the same category threaded through the chain, and it can be made exactly topological on the lattice by coupling it to the chain with the same F F F - and R R R -moves that define the category, in the spirit of the lattice topological defects of [19 , 20 ] and of the Temperley–Lieb constructions of [21 , 22 , 23 ] . (For the special case p = 4 p = 4 p = 4 , where C 4 C_4 C 4 is the Ising fusion category, a closely related one-dimensional model built directly on the fusion spaces of the Ising category, with an explicit vacancy label and the categorical symmetry realised exactly, has recently been studied by Roychowdhury and Wang [24 ] within the fusion-category lattice models of [25 ] .) Because the symmetry preserved along the flow is the one the lattice realises exactly, the whole flow can be followed with the symmetry imposed, sector by sector; and the perturbing field has a transparent meaning in this language, as a quantum-dimension-weighted combination of the projectors onto the two fusion channels of the defect anyon and its neighbouring domain wall.
The paper is organised as follows. § 2 reviews the KRW flows and establishes, in the continuum, which symmetry the flow D ( 2 , 1 ) → C D_{(2,1)} \to C D ( 2 , 1 ) → C preserves and why. § 3 builds the lattice model: the dilute height chain, the topological seam, the lattice symmetry operators, the clean deformation operator, the symmetry-resolved numerics and the observables used to characterise an interface. § 4 and § 5 report the results at p = 4 p = 4 p = 4 and p = 5 p = 5 p = 5 . § 6 reformulates the construction as a loop model valid at continuous p p p , and § 7 follows the family C p C_p C p for 4 ≤ p ≤ 5 4 \leq p \leq5 4 ≤ p ≤ 5 and discusses its endpoint at p = 4 p = 4 p = 4 . § 8 summarises and compares with conformal perturbation theory, integrable defect flows and string field theory.
The KRW flow and the symmetry it preserves
Perturbed topological defects: the work of Kormos, Runkel and Watts
A topological defect D D D of a conformal field theory is a line that commutes with both copies of the Virasoro algebra, [ L n , D ] = [ L ˉ n , D ] = 0 [L_n,D] = [\bar{L}_n,D] = 0 [ L n , D ] = [ L ˉ n , D ] = 0 ; it can be deformed freely as long as it does not cross field insertions, and its action on the bulk fields is encoded in a defect operator on the Hilbert space of the theory. In a diagonal minimal model the elementary topological defects are the Verlinde lines D ( r , s ) D_{(r,s)} D ( r , s ) [1 , 2 ] , one for each Kac label, which act on the primary field ϕ ( r ′ , s ′ ) \phi_{(r',s')} ϕ ( r ′ , s ′ ) by the Verlinde ratio S ( r , s ) , ( r ′ , s ′ ) / S ( 1 , 1 ) , ( r ′ , s ′ ) S_{(r,s),(r',s')}/S_{(1,1),(r',s')} S ( r , s ) , ( r ′ , s ′ ) / S ( 1 , 1 ) , ( r ′ , s ′ ) and fuse according to the minimal-model fusion rules. Their quantum dimension d ( r , s ) = S ( r , s ) , ( 1 , 1 ) / S ( 1 , 1 ) , ( 1 , 1 ) d_{(r,s)} = S_{(r,s),(1,1)}/S_{(1,1),(1,1)} d ( r , s ) = S ( r , s ) , ( 1 , 1 ) / S ( 1 , 1 ) , ( 1 , 1 ) is the g g g -function of the line [26 ] , and only the ( 1 , 1 ) (1,1) ( 1 , 1 ) line and the invertible line η = D ( p − 1 , 1 ) ≅ D ( 1 , p ) \eta= D_{(p-1,1)} \cong D_{(1,p)} η = D ( p − 1 , 1 ) ≅ D ( 1 , p ) have g = 1 g = 1 g = 1 . A defect D D D also carries its own fields, the defect fields, which live on the line; their space is the twisted sector of D ⊗ D ˉ D \otimes\bar{D} D ⊗ D ˉ (§ 2.2 ), and it contains in general chiral fields of weight ( h , 0 ) (h,0) ( h , 0 ) and ( 0 , h ) (0,h) ( 0 , h ) as well as non-chiral ones.
Kormos, Runkel and Watts [6 ] considered the perturbation of a topological defect by its relevant defect fields and asked where the resulting flow ends. Their starting point is the defect D ( 1 , 2 ) D_{(1,2)} D ( 1 , 2 ) of M ( p , p + 1 ) M(p,p+1) M ( p , p + 1 ) , whose relevant defect fields are the two chiral fields ϕ = ϕ ( 1 , 3 ) ⊗ 1 \phi= \phi_{(1,3)} \otimes1 ϕ = ϕ ( 1 , 3 ) ⊗ 1 and ϕ ˉ = 1 ⊗ ϕ ˉ ( 1 , 3 ) \bar{\phi} = 1 \otimes\bar{\phi}_{(1,3)} ϕ ˉ = 1 ⊗ ϕ ˉ ( 1 , 3 ) of weights ( h 1 , 3 , 0 ) (h_{1,3},0) ( h 1 , 3 , 0 ) and ( 0 , h 1 , 3 ) (0,h_{1,3}) ( 0 , h 1 , 3 ) , h 1 , 3 = ( p − 1 ) / ( p + 1 ) h_{1,3} = (p-1)/(p+1) h 1 , 3 = ( p − 1 ) / ( p + 1 ) , so that the perturbed defect is
D ( 1 , 2 ) ( λ l , λ r ) = D ( 1 , 2 ) exp ( λ l ∫ ϕ + λ r ∫ ϕ ˉ ) (1) D_{(1,2)}(\lambda_l,\lambda_r) = D_{(1,2)} \exp\left(\lambda_l \int\phi+ \lambda_r \int\bar{\phi}\right) \tag*{(1)} D ( 1 , 2 ) ( λ l , λ r ) = D ( 1 , 2 ) exp ( λ l ∫ ϕ + λ r ∫ ϕ ˉ ) (1)
with y = 1 − h 1 , 3 = 2 / ( p + 1 ) y = 1 - h_{1,3} = 2/(p+1) y = 1 − h 1 , 3 = 2/ ( p + 1 ) the RG eigenvalue of both couplings. Their analysis of this two-parameter space of flows combines three tools: exact arguments for the chiral flows, conformal perturbation theory in y y y , and the truncated conformal space approach for the non-chiral directions. The resulting flow diagram, their figure 1, is reproduced schematically in Figure 1 , and three of its features organise everything that follows.
Figure 1. The space of flows of the (1, 2) defect perturbed by λ l ϕ + λ r ϕ ˉ \lambda_l\phi+\lambda_r\bar{\phi} λ l ϕ + λ r ϕ ˉ , after figure 1 of Kormos, Runkel and Watts [6 ] (schematic; p > 3 p > 3 p > 3 ). The chiral flows along the axes end on the topological defects D ′ = D ( 2 , 1 ) D' = D_{(2,1)} D ′ = D ( 2 , 1 ) (positive couplings) and I I I (negative couplings). All flows in the first quadrant end on the new conformal defect C C C , those in the second and fourth quadrants on I I I , and the negative diagonal on the factorising defect F F F . The dashed lines are the boundaries of the basins: flows from D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) to C C C , driven by the spinful field ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) \phi_{(1,3)} \otimes\bar{\phi}_{(3,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) or its mirror image, and flows from D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) to I I I . The flow D ( 2 , 1 ) → C D_{(2,1)} \to C D ( 2 , 1 ) → C along the dashed boundary of the first quadrant is the one studied in this paper.
First, a purely chiral perturbation (λ r = 0 \lambda_r = 0 λ r = 0 or λ l = 0 \lambda_l = 0 λ l = 0 ) preserves pure transmission and the chiral algebra on the unperturbed side [7 , 8 ] , but need not preserve conformal invariance along the flow. A conformal endpoint of such a flow is topological; matching the action on the boundary state ∥ 1 , 1 ⟩ ⟩ \|1,1\rangle\rangle ∥1 , 1 ⟩⟩ fixes the endpoints ([6 ] , §3.1): the positive chiral flows of D ( 1 , 2 ) D_{(1,2)} D ( 1 , 2 ) end on D ′ = D ( 2 , 1 ) D' = D_{(2,1)} D ′ = D ( 2 , 1 ) and the negative ones on the identity defect I I I . The chiral flows have recently been put on a systematic footing in [27 , 28 ] . In the present paper the perturbation is non-chiral and the endpoint is a genuinely conformal defect (or, at p = 4 p = 4 p = 4 , a topological line reached by a non-chiral flow), a situation outside the scope of that analysis. Second, a non-chiral perturbation, λ l λ r ≠ 0 \lambda_l\lambda_r \ne0 λ l λ r = 0 , destroys the commutation with both Virasoro algebras and produces a genuinely conformal, non-topological interface. KRW find that every flow in the second and fourth quadrants ends on I I I , that the negative diagonal λ l = λ r < 0 \lambda_l = \lambda_r < 0 λ l = λ r < 0 ends on a factorising defect F F F , the superposition of the ( r , 1 ) (r,1) ( r , 1 ) boundary conditions, and that the whole first quadrant flows to a new conformal defect C C C , the subject of this paper; near c = 1 c = 1 c = 1 they thus find six fixed points, five of which are known defects and one, C C C , a new one [6 ] , abstract and §6). This picture is established for large p p p and confirmed by TCSA down to M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) ; for M ( 4 , 5 ) M(4,5) M ( 4 , 5 ) they could not locate the non-chiral fixed point and left its nature open ([6 ] , §6.2.4). For C C C they computed the g g g -function and the reflection coefficient in conformal perturbation theory in powers of y y y , i.e. at large p p p : ln g C − ln g D ( 2 , 1 ) = − π 2 4 y 3 + O ( y 4 ) \ln g_C - \ln g_{D_{(2,1)}} = -\frac{\pi^2}{4}y^3 + O(y^4) ln g C − ln g D ( 2 , 1 ) = − 4 π 2 y 3 + O ( y 4 ) [6 , 14 ] and R C = 1 − T = 9 8 π 2 y 4 + O ( y 5 ) R_C = 1 - \mathcal{T} = \frac{9}{8}\pi^2 y^4 + O(y^5) R C = 1 − T = 8 9 π 2 y 4 + O ( y 5 ) [12 ] . Third, and most relevant here, the boundaries of these basins of attraction are themselves flows out of D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) : the boundary of the first quadrant consists of flows from D ′ D' D ′ to C C C , generated by the spinful defect field ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) \phi_{(1,3)} \otimes\bar{\phi}_{(3,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) on the side λ r > λ l \lambda_r > \lambda_l λ r > λ l and by its mirror image ϕ ( 3 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) \phi_{(3,3)} \otimes\bar{\phi}_{(1,3)} ϕ ( 3 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) on the side λ l > λ r \lambda_l > \lambda_r λ l > λ r , while the boundary of the second and fourth quadrants is a flow from D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) to I I I ([6], §6.2.3). In other words, once D ( 1 , 2 ) D_{(1,2)} D ( 1 , 2 ) has reached D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) along a chiral axis, the remaining perturbation is the first-order image of the other chiral field, the spinful field ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) \phi_{(1,3)} \otimes\bar{\phi}_{(3,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) of dimension Δ = h 1 , 3 + h 3 , 3 = 1 − 2 ( p − 1 ) / ( p ( p + 1 ) ) \Delta= h_{1,3} + h_{3,3} = 1 - 2(p - 1)/(p(p + 1)) Δ = h 1 , 3 + h 3 , 3 = 1 − 2 ( p − 1 ) / ( p ( p + 1 )) , and D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) perturbed by it flows to C C C for one sign of the coupling and to I I I for the other. It is this flow, D ( 2 , 1 ) → C D_{(2,1)} \to C D ( 2 , 1 ) → C , that we realise on the lattice.
KRW also observed that the perturbed defects (1) commute with the topological lines D ( r , 1 ) D_{(r,1)} D ( r , 1 ) for all values of the couplings (their §2.4.2), and used this commutation to constrain the possible infrared fixed points (their §3.2). For the flow out of D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) the corresponding statement is both true and subtle, and since it is what allows the lattice construction we devote the rest of this section to it: the preserved symmetry is not “the line D ( r , 1 ) D_{(r,1)} D ( r , 1 ) ” but a specific way of attaching that line to the defect, which the deformation selects.
Defect operators of D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 )
In radial quantisation about a point of the ring, the states of a ring crossed by a line D D D at two antipodal points are the operators sitting on a straight D D D line through the origin, i.e. the defect operators of D D D . For a topological D D D this space is [2 , 4 ]
H D × D = ⨁ k N D D k H k , H k = ⨁ i , j N k i j V i ⊗ V ‾ j , (2) \mathcal{H}_{D\times D}=\bigoplus_k N^k_{DD}\mathcal{H}_k,\qquad\mathcal{H}_k=\bigoplus_{i,j}N^j_{ki}V_i\otimes\overline{V}_j, \tag*{(2)} H D × D = k ⨁ N D D k H k , H k = i , j ⨁ N k i j V i ⊗ V j , (2)
where H k \mathcal{H}_k H k is the k k k -twisted sector of the bulk theory (H 1 \mathcal{H}_1 H 1 is the diagonal spectrum), i = ( r ′ , s ) i=(r',s) i = ( r ′ , s ) , j = ( r ′ ′ , s ) j=(r'',s) j = ( r ′′ , s ) label Virasoro representations and L L L are the fusion multiplicities. For D = D ( 2 , 1 ) D=D_{(2,1)} D = D ( 2 , 1 ) , D × D = D ( 1 , 1 ) + D ( 3 , 1 ) D\times D=D_{(1,1)}+D_{(3,1)} D × D = D ( 1 , 1 ) + D ( 3 , 1 ) , so the ring carries H 1 ⊕ H ( 3 , 1 ) \mathcal{H}_1\oplus\mathcal{H}_{(3,1)} H 1 ⊕ H ( 3 , 1 ) : the diagonal fields V i ⊗ V ‾ i V_i\otimes\overline{V}_i V i ⊗ V i and the fields V ( r ′ , s ) ⊗ V ‾ ( r ′ ′ , s ) V_{(r',s)}\otimes\overline{V}_{(r'',s)} V ( r ′ , s ) ⊗ V ( r ′′ , s ) with r ′ ′ ∈ r ′ ⊗ 3 r''\in r'\otimes3 r ′′ ∈ r ′ ⊗ 3 . A ring crossed only once by D D D carries instead the twisted sector H D = ⨁ i , j N D i j V i ⊗ V ‾ j \mathcal{H}_D=\bigoplus_{i,j}N^j_{Di}V_i\otimes\overline{V}_j H D = ⨁ i , j N D i j V i ⊗ V j , the space of operators at which the line ends; this is what the single-strand computation of § 4 uses.
Lassos: a crossing is a choice
A wrapped ( r , 1 ) (r,1) ( r , 1 ) line x x x on the ring must cross the defect line at every strand, and a crossing is not canonical. It is a junction u ∈ Hom ( x ⊗ D , D ⊗ x ) u\in\operatorname{Hom}(x\otimes D,D\otimes x) u ∈ Hom ( x ⊗ D , D ⊗ x ) , a space whose dimension is the number of fusion channels of x ⊗ D x\otimes D x ⊗ D : it is spanned by the fuse-and-split junctions u D ′ u_{D'} u D ′ , in which x x x and D D D fuse into an intermediate line D ′ ∈ x ⊗ D D'\in x\otimes D D ′ ∈ x ⊗ D and split again. For the KRW starting point D ( 1 , 2 ) D_{(1,2)} D ( 1 , 2 ) and x = ( r , 1 ) x=(r,1) x = ( r , 1 ) the fusion ( r , 1 ) ⊗ ( 1 , 2 ) = ( r , 2 ) (r,1)\otimes(1,2)=(r,2) ( r , 1 ) ⊗ ( 1 , 2 ) = ( r , 2 ) has a single channel, so the junction is unique up to normalisation. For our starting point D ( 2 , 1 ) = w D_{(2,1)}=\mathbf{w} D ( 2 , 1 ) = w and x = w x=\mathbf{w} x = w one has w ⊗ w = 1 ⊕ adj \mathbf{w}\otimes\mathbf{w}=\mathbf{1}\oplus\operatorname{adj} w ⊗ w = 1 ⊕ adj , so the junction space is two-dimensional (for x = η x=\eta x = η it is again one-dimensional, so the crossing is unique up to a scalar; the over- and under-crossings can still differ by a monodromy phase). A loop of x x x with a chosen junction at each crossing is a lasso , an element of the tube algebra of D D D [5 , 29 ] . As long as D D D is topological every lasso commutes with the Hamiltonian, but different lassos are different operators. Two of them are distinguished: the over- and under-crossings, in which x x x passes in front of or behind D D D at the crossing. They are the two half-braidings of x x x with D D D —the junctions compatible with moving x x x past the fusion vertices of the defect—and this is what makes them special among all junctions: a lasso built from either one acts diagonally on the defect fields, with eigenspaces the objects ( i , j ˉ ) (i,\bar{j}) ( i , j ˉ ) of the Drinfeld centre C ⊠ C ‾ \mathcal{C}\boxtimes\overline{\mathcal{C}} C ⊠ C , the irreducible representations of the tube algebra [29 , 30 ] , whereas a generic junction need not act diagonally within an irreducible tube-algebra representation or on the corresponding defect multiplicity spaces. It does not mix distinct irreducible representation blocks. We denote the two half-braidings by B B B and B − 1 B^{-1} B − 1 ; on the lattice they are the braid matrix B = F R F − 1 B=FRF^{-1} B = F R F − 1 of § 3 and its inverse, and their explicit expansion in the channel basis is given in § 3.3 , eq. (19) . On a defect field V i ⊗ V ‾ j V_i\otimes\overline{V}_j V i ⊗ V j the over-crossing lasso acts by the Verlinde eigenvalue of the left label, S x i / S 0 i S_{xi}/S_{0i} S x i / S 0 i , and the under-crossing lasso by that of the right label, S x j / S 0 j S_{xj}/S_{0j} S x j / S 0 j .
The reason a crossing reads one chiral label is most transparent in the three-dimensional description of the diagonal theory [3 , 4 ] , Figure 2 . There the CFT on a surface Σ \Sigma Σ is the Reshetikhin–Turaev topological field theory on the slab Σ × [ − 1 , 1 ] \Sigma\times[-1,1] Σ × [ − 1 , 1 ] : the holomorphic degrees of freedom live on the top sheet Σ × { + 1 } \Sigma\times\{+1\} Σ × { + 1 } , the antiholomorphic ones on the bottom sheet Σ × { − 1 } \Sigma\times\{-1\} Σ × { − 1 } , and all defect lines, the symmetry lines x x x as well as the defect D D D , are Wilson lines in the middle layer Σ × { 0 } \Sigma\times\{0\} Σ × { 0 } . A field with labels ( i , j ˉ ) (i,\bar{j}) ( i , j ˉ ) is the meeting point of a Wilson line i i i coming down from the top sheet and a Wilson line j j j coming up from the bottom sheet; for a bulk field of the diagonal theory i = j i=j i = j , for a defect field of D D D the two labels are independent and the lines meet on D D D . A loop of x x x in the middle layer around the insertion point crosses D D D twice, and “over” and “under” at a crossing are literally the two sides of the middle layer in the third direction. If x x x passes over D D D at both crossings the loop can be lifted off the middle layer into the holomorphic half of the slab, where it links only the i i i -line and contracts to S x i / S 0 i S_{xi}/S_{0i} S x i / S 0 i ; if it passes under at both crossings it can be pushed into the antiholomorphic half, where it links only the j j j -line and gives S x j / S 0 j S_{xj}/S_{0j} S x j / S 0 j . A loop with one over- and one under-crossing can be moved to neither side, links both lines at once, and need not act as a scalar within a tube-algebra representation; this does not imply mixing between distinct irreducible representation blocks. For a bulk field i = j i=j i = j and the two lassos agree, which is why the action of a Verlinde line on bulk fields never requires a junction to be specified; for a defect field with i ≠ j i \ne j i = j they differ. On the lattice, moving the wrapped line through the strand vertex with the bare F F F -move keeps it in front of the strand in the fusion tree, which is the over-crossing, while passing behind the strand requires the ribbon twists θ a \theta_a θ a of the labels (§ 3.3 ); which of the two reads the holomorphic label is a matter of orientation conventions, fixed here empirically by the undeformed-ring check of Table 1.
Figure 2. Why a lasso reads one chiral label. In the Reshetikhin–Turaev description on Σ × [ − 1 , 1 ] \Sigma\times[-1,1] Σ × [ − 1 , 1 ] the holomorphic data live on the top sheet, the antiholomorphic data on the bottom sheet, and all defect lines in the middle layer; a defect field ( i , j ˉ ) (i,\bar{j}) ( i , j ˉ ) at P P P is the meeting point of a Wilson line i i i from above and j ˉ \bar{j} j ˉ from below. (a) A loop of x x x passing above D D D at both crossings can be lifted into the holomorphic half, where it links only i i i and gives S x i / S 0 i S_{xi}/S_{0i} S x i / S 0 i . (b) Passing below at both crossings it can be pushed into the antiholomorphic half, where it links only j ˉ \bar{j} j ˉ and gives S x j ˉ / S 0 j ˉ S_{x\bar{j}}/S_{0\bar{j}} S x j ˉ / S 0 j ˉ .
The deformation fixes the lasso
Once D D D is deformed by a defect field ϕ \phi ϕ , a lasso stays a symmetry only if the line x x x can be moved through the insertions of ϕ \phi ϕ along D D D with the chosen junction,
u ∘ ( i d x ⊗ ϕ ) = ( ϕ ⊗ i d x ) ∘ u , (3) u \circ(\mathrm{id}_x \otimes\phi) = (\phi\otimes\mathrm{id}_x) \circ u, \tag*{(3)} u ∘ ( id x ⊗ ϕ ) = ( ϕ ⊗ id x ) ∘ u , (3)
read with time running along D D D : the insertion of ϕ \phi ϕ on D D D is a point operator ϕ : D → D \phi: D \to D ϕ : D → D , and i d x ⊗ ϕ \mathrm{id}_x \otimes\phi id x ⊗ ϕ (resp. ϕ ⊗ i d x \phi\otimes\mathrm{id}_x ϕ ⊗ id x ) denotes this insertion while x x x runs to the left (resp. right) of D D D , i.e. before (resp. after) the crossing u : x ⊗ D → D ⊗ x u : x \otimes D \to D \otimes x u : x ⊗ D → D ⊗ x . Both sides are maps x ⊗ D → D ⊗ x x \otimes D \to D \otimes x x ⊗ D → D ⊗ x , and the equation states that the crossing can be slid along D D D past the insertion point of ϕ \phi ϕ . Closing x x x into a small loop around the insertion shows that this holds exactly when ϕ \phi ϕ is neutral under that lasso, Y x ϕ = d x ϕ Y_x\phi= d_x\phi Y x ϕ = d x ϕ : the two crossings of the loop can then be brought together and cancel, leaving a free loop worth d x d_x d x (Figure 3 ). Since the over-crossing reads the left label of ϕ \phi ϕ and the under-crossing its right label, the rule is simply that a deformation selects the lasso that reads its neutral label. For the driving field ϕ = ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) \phi= \phi_{(1,3)} \otimes\bar{\phi}_{(3,3)} ϕ = ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) on D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) (holomorphic label first) the criterion is decisive: the left label ( 1 , 3 ) (1,3) ( 1 , 3 ) has r ′ = 1 r' = 1 r ′ = 1 and hence the same eigenvalue d x d_x d x as the identity under every ( r , 1 ) (r,1) ( r , 1 ) line, so ϕ \phi ϕ is neutral under the over-crossing lasso; the right label ( 3 , 3 ) (3,3) ( 3 , 3 ) has Y ( 2 , 1 ) = − ϕ − 1 ≠ ϕ Y_{(2,1)} = -\phi^{-1} \ne\phi Y ( 2 , 1 ) = − ϕ − 1 = ϕ at p = 5 p=5 p = 5 , where ( 3 , 3 ) = ( 2 , 3 ) (3,3) = (2,3) ( 3 , 3 ) = ( 2 , 3 ) , and Y ( 2 , 1 ) = − 2 ≠ 2 Y_{(2,1)} = -\sqrt{2} \ne\sqrt{2} Y ( 2 , 1 ) = − 2 = 2 at p = 4 p=4 p = 4 , where ( 3 , 3 ) = ( 1 , 2 ) (3,3) = (1,2) ( 3 , 3 ) = ( 1 , 2 ) , so the under-crossing lasso is broken. The symmetry preserved along the flow D ( 2 , 1 ) → C D_{(2,1)} \to C D ( 2 , 1 ) → C is therefore not “the line x x x ”, which can be attached to the defect in two ways, but the specific lasso that crosses over the defect at every intersection. The mirror field ϕ ( 3 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) \phi_{(3,3)} \otimes\bar{\phi}_{(1,3)} ϕ ( 3 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) , charged on the left and neutral on the right, fixes the under-crossing lasso and breaks the over-crossing one; the two deformations are exchanged by the swap of the two sides of the defect, which is why the flow D ( 2 , 1 ) → C D_{(2,1)} \to C D ( 2 , 1 ) → C is not swap-symmetric while both its endpoints are. The two fields are the two members of the driving doublet at Δ = 11 / 15 \Delta= 11/15 Δ = 11/15 in Table 1, which the over-crossing lasso places in different sectors (a + a+ a + and c + c+ c + ): the lasso preserved by the deformation is the one that distinguishes the deformation from its mirror image, and a deformation by the sum of the two fields would preserve neither lasso.
Figure 3. Why the deformation selects a lasso. (a) The junction u ∈ Hom ( x ⊗ D , D ⊗ x ) u\in\operatorname{Hom}(x\otimes D,D\otimes x) u ∈ Hom ( x ⊗ D , D ⊗ x ) : unique for the KRW starting point D ( 1 , 2 ) D_{(1,2)} D ( 1 , 2 ) , where ( r , 1 ) ⊗ ( 1 , 2 ) = ( r , 2 ) (r,1)\otimes(1,2)=(r,2) ( r , 1 ) ⊗ ( 1 , 2 ) = ( r , 2 ) has one channel, and two-dimensional for D ( 2 , 1 ) = w D_{(2,1)}=w D ( 2 , 1 ) = w , where w ⊗ w = 1 ⊕ adj w\otimes w=1\oplus\operatorname{adj} w ⊗ w = 1 ⊕ adj , with the over- and under-crossings B B B , B − 1 B^{-1} B − 1 as the distinguished choices. (b) A lasso survives the deformation only if x x x can be moved through every insertion of ϕ \phi ϕ with its junction, eq. (3) , i.e. only if ϕ \phi ϕ is neutral under it. For ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) \phi_{(1,3)}\otimes\bar{\phi}_{(3,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) only the over-crossing lasso, which reads the neutral left label, survives; for the mirror field only the under-crossing one.
With two strands deformed by the same field, the preserved operator is the lasso with the over-crossing at both strands, and with this choice (2) is diagonalised by the left labels:
Y ( r , 1 ) ∣ V ( r ′ , s ) ⊗ V ‾ ( r ′ ′ , s ) = S ( r , 1 ) , ( r ′ , s ) S ( 1 , 1 ) , ( r ′ , s ) = ( − 1 ) ( r − 1 ) s sin ( π r r ′ p + 1 p ) sin ( π r ′ r p + 1 p ) (4) Y_{(r,1)}\big|_{V_{(r',s)}\otimes\overline{V}_{(r'',s)}}=\frac{S_{(r,1),(r',s)}}{S_{(1,1),(r',s)}}=(-1)^{(r-1)s}\frac{\sin\left(\pi rr'\frac{p+1}{p}\right)}{\sin\left(\pi r'r\frac{p+1}{p}\right)} \tag*{(4)} Y ( r , 1 ) V ( r ′ , s ) ⊗ V ( r ′′ , s ) = S ( 1 , 1 ) , ( r ′ , s ) S ( r , 1 ) , ( r ′ , s ) = ( − 1 ) ( r − 1 ) s sin ( π r ′ r p p + 1 ) sin ( π r r ′ p p + 1 ) (4)
independently of r ′ ′ r'' r ′′ and of the sector H k \mathcal{H}_k H k the field comes from. The eigenvalues are the quantum dimensions of S U ( 2 ) p − 2 SU(2)_{p-2} S U ( 2 ) p − 2 up to signs: at p = 4 p=4 p = 4 , Y ( 2 , 1 ) ∈ { + 2 , 0 , − 2 } Y_{(2,1)}\in\{+\sqrt{2},0,-\sqrt{2}\} Y ( 2 , 1 ) ∈ { + 2 , 0 , − 2 } and η = Y ( 3 , 1 ) = ± 1 \eta=Y_{(3,1)}=\pm1 η = Y ( 3 , 1 ) = ± 1 ; at p = 5 p=5 p = 5 , Y ( 2 , 1 ) = ± φ , ± φ − 1 Y_{(2,1)}=\pm\varphi,\pm\varphi^{-1} Y ( 2 , 1 ) = ± φ , ± φ − 1 and η = Y ( 4 , 1 ) = ± 1 \eta=Y_{(4,1)}=\pm1 η = Y ( 4 , 1 ) = ± 1 , with Y ( 2 , 1 ) = Y τ F i b η Y_{(2,1)}=Y_{\mathrm{\tau Fib}}\eta Y ( 2 , 1 ) = Y τ Fib η .
Since the preserved lasso commutes with the deformed Hamiltonian, these eigenvalues are conserved along the whole flow: every operator of the infrared interface inherits a ( Y ( 2 , 1 ) , η ) (Y_{(2,1)},\eta) ( Y ( 2 , 1 ) , η ) label from the ultraviolet defect operator it descends from, so that the spectrum of C C C can be organised, and its finite-size extrapolation performed, sector by sector (§ 4 , § 5 ); at the endpoint η \eta η , x ⊗ η x\otimes\eta x ⊗ η is simple, so the crossing junction is unique up to normalization. Any scalar difference between the over- and under-crossings must still be tracked; in the two-crossing defect-operator geometry the corresponding monodromy factors cancel.
Kormos, Runkel and Watts prove, by a matrix-element argument ([6 ] , §2.4.2), that the perturbed defect (1) commutes with all the lines D ( r , 1 ) D_{(r,1)} D ( r , 1 ) for every value of the couplings, and use this to constrain the infrared fixed points ([6 ] , §3.2). In their setting the statement is unambiguous for two reasons that the preceding analysis makes explicit: the junction space of x = ( r , 1 ) x=(r,1) x = ( r , 1 ) with D ( 1 , 2 ) D_{(1,2)} D ( 1 , 2 ) is one-dimensional, so there is only one lasso to speak of, and the chiral perturbing fields carry the labels ( 1 , 3 ) ⊗ ( 1 , 1 ) (1,3)\otimes(1,1) ( 1 , 3 ) ⊗ ( 1 , 1 ) and ( 1 , 1 ) ⊗ ( 1 , 3 ) (1,1)\otimes(1,3) ( 1 , 1 ) ⊗ ( 1 , 3 ) , both of ( 1 , s ) (1,s) ( 1 , s ) type and hence neutral under D ( r , 1 ) D_{(r,1)} D ( r , 1 ) whichever side of the defect the line reads, so that (3) holds automatically. The junction never has to be discussed. Our flow starts from D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) , a member of the same family as the symmetry lines, so that the junction space is two-dimensional, and it is driven by a spinful field with one charged label, so that the two lassos are inequivalent on the deformed defect. The requirement that the deformation commute with the symmetry therefore does more than it does in KRW: it selects one element of the junction space, and this is the precise sense in which “the ( r , 1 ) (r,1) ( r , 1 ) symmetry is preserved along the flow D ( 2 , 1 ) → C D_{(2,1)}\to C D ( 2 , 1 ) → C ”.
This is also what makes a lattice realisation possible. The preserved symmetry is a line of the fusion category C p \mathcal{C}_p C p attached to the defect by a definite junction, so it can be realised exactly on a lattice whose degrees of freedom are the fusion spaces of that category, with the junction transcribed into the F F F -moves that carry the wrapped line across the defect. The dilute A p − 1 A_{p-1} A p − 1 RSOS chain is such a lattice, and the next section constructs on it the defect, the wrapped lines with the over-crossing junction, and the deformation; the conserved charges (4) are then imposed exactly, sector by sector, along the whole flow.
The lattice model
The dilute RSOS chain
Degrees of freedom. We use the following convention throughout the paper: the height variables live on the bonds of the chain, and the anyons (the domain walls, or their absence) live on the sites between consecutive bonds. Consider a ring of L L L bonds. Each bond j j j carries an integer height
a j ∈ { 1 , … , p − 1 } (three values for tricritical Ising, four for M ( 5 , 6 ) ) , (5) a_j \in\{1,\ldots,p-1\} \qquad\text{(three values for tricritical Ising, four for }M(5,6)\text{)}, \tag*{(5)} a j ∈ { 1 , … , p − 1 } (three values for tricritical Ising, four for M ( 5 , 6 ) ) , (5)
subject to the admissibility rule that neighbouring heights differ by at most one,
∣ a j − a j + 1 ∣ ≤ 1 at every site . (6) \lvert a_j-a_{j+1}\rvert\le1 \qquad\text{at every site}. \tag*{(6)} ∣ a j − a j + 1 ∣ ≤ 1 at every site . (6)
Between two consecutive bonds there is a site . A site carries no variable of its own; we call it empty if the two heights around it are equal and occupied if they differ,
site between bonds a , b : a = b ⇒ empty , a ≠ b ⇒ occupied . (7) \text{site between bonds }a,b:\qquad a=b\Rightarrow\text{ empty},\qquad a\ne b\Rightarrow\text{ occupied}. \tag*{(7)} site between bonds a , b : a = b ⇒ empty , a = b ⇒ occupied . (7)
An occupied site is a domain wall between two height domains. The word “dilute” refers to the fact that empty sites are allowed (in the ordinary, dense A p − 1 A_{p-1} A p − 1 RSOS model every site would be a wall) see Figure 4 . The dense A p − 1 A_{p-1} A p − 1 models are those of Andrews, Baxter and Forrester [31 ] , whose critical points realise the unitary minimal models [32 , 33 ] ; the dilute A p − 1 A_{p-1} A p − 1 models were constructed by Warnaar, Nienhuis and Seaton [15 ] from the O ( n ) O(n) O ( n ) loop model of [34 , 35 ] .
Figure 4. The dilute A p − 1 A_{p-1} A p − 1 RSOS chain (here p = 5 p=5 p = 5 , four heights). Bottom: each bond j j j carries a height a j ∈ { 1 , … , p − 1 } a_j \in\{1,\ldots,p-1\} a j ∈ { 1 , … , p − 1 } (boxes), neighbours differing by at most one; the site between two bonds is occupied by a w (filled red dot) when the heights differ and empty (open circle) when they are equal. Top: the same configuration drawn as a height profile—the heights are domains and every domain wall carries a w. The Hamiltonian (11) hops these walls and creates/annihilates them in pairs.
Local Hilbert space. The heights are the only degrees of freedom, so the local Hilbert space of a bond is
H j = span { ∣ a j ⟩ : a j = 1 , … , p − 1 } , dim H j = p − 1 , (8) \mathcal{H}_j = \operatorname{span}\{ \lvert a_j \rangle: a_j = 1,\ldots,p-1 \}, \qquad\dim\mathcal{H}_j = p-1, \tag*{(8)} H j = span {∣ a j ⟩ : a j = 1 , … , p − 1 } , dim H j = p − 1 , (8)
a three-state space at p = 4 p=4 p = 4 and a four-state space at p = 5 p=5 p = 5 . The many-body space is not the full tensor product ⨂ j H j \bigotimes_j \mathcal{H}_j ⨂ j H j but its admissible subspace, spanned by the height configurations obeying (6) at every site. Every operator introduced below (the Hamiltonian, the seam, the wrapped lines, the deformation) preserves this subspace. In exact diagonalisation we enumerate the admissible configurations and work in that basis; in DMRG we keep the full ( p − 1 ) (p-1) ( p − 1 ) -dimensional space on each bond and let the Hamiltonian enforce the constraint (forbidden height jumps have no matrix elements, so an admissible state never leaves the admissible subspace).
Anyon-chain reading. The same data have a second, equivalent description that will be essential for the defect. Identify a height a a a with an object of C p C_p C p , using the S U ( 2 ) p − 2 SU(2)_{p-2} S U ( 2 ) p − 2 spin label j = a − 1 2 j=\frac{a-1}{2} j = 2 a − 1 , and let w denote the spin-1 2 \frac{1}{2} 2 1 object (category label 1), of quantum dimension
d w = 2 cos π p = S 2 S 1 , S a ≡ sin π a p (9) d_{\mathrm{w}} = 2\cos\frac{\pi}{p} = \frac{S_2}{S_1}, \qquad S_a \equiv\sin\frac{\pi a}{p} \tag*{(9)} d w = 2 cos p π = S 1 S 2 , S a ≡ sin p π a (9)
(d w = 2 d_{\mathrm{w}}=\sqrt{2} d w = 2 at p = 4 p=4 p = 4 , d w = φ = 1.618 d_{\mathrm{w}}=\varphi=1.618 d w = φ = 1.618 at p = 5 p=5 p = 5 ; S a S_a S a is the vector of quantum dimensions, which will appear in all the weights below). Then the admissibility rule (6) says exactly that a j + 1 ∈ a j ⊗ ( 1 ⊕ w ) a_{j+1}\in a_j\otimes(1\oplus\mathrm{w}) a j + 1 ∈ a j ⊗ ( 1 ⊕ w ) : an empty site fuses the trivial object into the chain, an occupied site fuses a w. In other words, the height configuration { a j } \{a_j\} { a j } is the fusion-tree basis of a chain of w anyons sitting on the occupied sites, with the heights recording the intermediate total charge as one moves along the ring. “A w lives on every domain wall” is the precise statement, not just a picture. The two kinds of labels play different categorical roles, and it is worth stating them precisely, since this is the structure the defect construction will exploit. The site labels 1 and w w w (and later the defect labels x x x ) are objects of the fusion category C = C p \mathcal{C} = \mathcal{C}_p C = C p : they are the things that fuse, braid and carry quantum dimensions. The heights, on the other hand, are objects of a module category M \mathcal{M} M over C \mathcal{C} C : they are acted upon by the objects of C \mathcal{C} C (a ↦ a ⊗ w a \mapsto a \otimes w a ↦ a ⊗ w ), and the admissibility graph of the height model is the fusion graph of this action. For the A p − 1 A_{p-1} A p − 1 models M \mathcal{M} M is C \mathcal{C} C itself, acting on itself by fusion (the regular module category), which is why the heights can be identified with S U ( 2 ) p − 2 SU(2)_{p-2} S U ( 2 ) p − 2 labels at all; in the D D D - and E E E -type height models of Pasquier [33 ] , the heights instead label vertices of other ADE graphs, illustrating why the action on heights must be specified separately from the fusion rules, and the module-category structure is the general statement that an anyon chain is defined by a pair ( C , M ) (\mathcal{C},\mathcal{M}) ( C , M ) [17 , 18 ] . This is the lattice counterpart of the correspondence between module categories and conformal boundary conditions, or equivalently between module categories and the ways a chiral algebra can be completed to a full CFT [36 , 37 ] ; on the lattice it is the organising principle of the fusion-category lattice models of [20 , 25 ] and of the general framework of dualities and non-invertible symmetries in one-dimensional quantum lattice models developed by Lootens, Delcamp, Ortiz and Verstraete [38 , 39 ] , which is reviewed in Delcamp’s IHES lectures [40 ] .
The Hamiltonian. The critical dynamics is a nearest-neighbour Hamiltonian
H 0 = − ∑ j = 1 L h j , (10) H_0 = -\sum_{j=1}^{L} h_j, \tag*{(10)} H 0 = − j = 1 ∑ L h j , (10)
where h j h_j h j acts on three consecutive bonds—a left spectator l = a j − 1 l = a_{j-1} l = a j − 1 , the middle bond a = a j a = a_j a = a j and a right spectator r = a j + 1 r = a_{j+1} r = a j + 1 —and only ever changes the middle height, a → a ′ a \to a' a → a ′ , keeping the result admissible. With P l = ∣ l ⟩ ⟨ l ∣ P_l = |l\rangle\langle l| P l = ∣ l ⟩ ⟨ l ∣ the projector on a spectator height,
h j = ∑ l , r ∑ a , a ′ W a → a ′ ( l , r ) P l ∣ a ′ ⟩ ⟨ a ∣ P r , (11) h_j = \sum_{l,r}\sum_{a,a'} W_{a\to a'}^{(l,r)} P_l \lvert a'\rangle\langle a\rvert P_r, \tag*{(11)} h j = l , r ∑ a , a ′ ∑ W a → a ′ ( l , r ) P l ∣ a ′ ⟩ ⟨ a ∣ P r , (11)
so h j h_j h j is a ( p − 1 ) 3 × ( p − 1 ) 3 (p-1)^3 \times(p-1)^3 ( p − 1 ) 3 × ( p − 1 ) 3 matrix specified by the numbers W a → a ′ ( l , r ) W_{a\to a'}^{(l,r)} W a → a ′ ( l , r ) . These are the Hamiltonian limit of the dilute A p − 1 A_{p-1} A p − 1 face weights of Nienhuis [34 ] and of Warnaar, Nienhuis and Seaton [15 , 16 ] , taken at the critical branch with crossing parameter
λ = π ( p − 1 ) 4 p ( λ = 3 π 16 at p = 4 , λ = π 5 at p = 5 ) . (12) \lambda= \frac{\pi(p-1)}{4p} \qquad\left(\lambda= \frac{3\pi}{16}\text{ at }p=4,\quad\lambda= \frac{\pi}{5}\text{ at }p=5\right). \tag*{(12)} λ = 4 p π ( p − 1 ) ( λ = 16 3 π at p = 4 , λ = 5 π at p = 5 ) . (12)
They depend only on the shape of the move, i.e. on which of the two internal sites are domain walls before (o L = [ l ≠ a ] o_L = [l \ne a] o L = [ l = a ] , o R = [ a ≠ r ] o_R = [a \ne r] o R = [ a = r ] ) and after (o L ′ = [ l ≠ a ′ ] o'_L = [l \ne a'] o L ′ = [ l = a ′ ] , o R ′ = [ a ′ ≠ r ] o'_R = [a' \ne r] o R ′ = [ a ′ = r ] ) the move, and they take one of six values:
empty ( a = a ′ , o L = o R = 0 ) : w ∅ = sin 3 λ N 0 , one wall ( a = a ′ , o L ≠ o R ) : w ∥ = − sin 2 λ cos 3 λ N 0 , two walls ( a = a ′ , o L = o R = 1 ) : w = = − sin 5 λ N 0 , hop ( o L o R = 0 , wall moves side ) : w ↷ = sin 3 λ N 0 , create/annihilate ( l = r , one bump ) : w ∪ ∩ = sin 2 λ N 0 S m S l , TL ( l = r , both walls, any a ′ ) : w e = − sin λ N 0 S a S a ′ S l . (13) \begin{aligned}
\text{empty }(a=a',\ o_L=o_R=0):\quad&w_{\varnothing}=\frac{\sin3\lambda}{\mathcal{N}_0},\\
\text{one wall }(a=a',\ o_L\ne o_R):\quad&w_{\parallel}=-\frac{\sin2\lambda\cos3\lambda}{\mathcal{N}_0},\\
\text{two walls }(a=a',\ o_L=o_R=1):\quad&w_{=}=-\frac{\sin5\lambda}{\mathcal{N}_0},\\
\text{hop }(o_Lo_R=0,\ \text{wall moves side}):\quad&w_{\curvearrowright}=\frac{\sin3\lambda}{\mathcal{N}_0},\\
\text{create/annihilate }(l=r,\ \text{one bump}):\quad&w_{\cup\cap}=\frac{\sin2\lambda}{\mathcal{N}_0}\sqrt{\frac{S_m}{S_l}},\\
\text{TL }(l=r,\ \text{both walls, any }a'):\quad&w_e=-\frac{\sin\lambda}{\mathcal{N}_0}\frac{\sqrt{S_aS_{a'}}}{S_l}. \tag*{(13)}
\end{aligned} empty ( a = a ′ , o L = o R = 0 ) : one wall ( a = a ′ , o L = o R ) : two walls ( a = a ′ , o L = o R = 1 ) : hop ( o L o R = 0 , wall moves side ) : create/annihilate ( l = r , one bump ) : TL ( l = r , both walls, any a ′ ) : w ∅ = N 0 sin 3 λ , w ∥ = − N 0 sin 2 λ cos 3 λ , w = = − N 0 sin 5 λ , w ↷ = N 0 sin 3 λ , w ∪∩ = N 0 sin 2 λ S l S m , w e = − N 0 sin λ S l S a S a ′ . (13)
with N 0 = sin 2 λ sin 3 λ \mathcal{N}_0=\sin2\lambda\sin3\lambda N 0 = sin 2 λ sin 3 λ a convenient overall normalisation. The six moves are drawn in Figure 5 . Two points of reading are worth making explicit, since they decide the matrix. In the creation/annihilation weight m m m is the bump height, the middle height of the configuration in which both sites are walls, whichever side of the move it is on: W l → m ( l , l ) = W m → l ( l , l ) = sin 2 λ N 0 S m / S l W^{(l,l)}_{l\to m}=W^{(l,l)}_{m\to l}=\frac{\sin2\lambda}{\mathcal{N}_0}\sqrt{S_m/S_l} W l → m ( l , l ) = W m → l ( l , l ) = N 0 s i n 2 λ S m / S l for ∣ m − l ∣ = 1 |m-l|=1 ∣ m − l ∣ = 1 , which makes the Hermiticity manifest. The Temperley–Lieb weight applies to every pair of configurations with l = r l=r l = r and both sites occupied, including a ′ = a a'=a a ′ = a , so that on such a configuration it adds to the diagonal: W a → a ( l , l ) = w e = − sin λ N 0 S a S l W^{(l,l)}_{a\to a}=w_e=-\frac{\sin\lambda}{\mathcal{N}_0}\frac{S_a}{S_l} W a → a ( l , l ) = w e = − N 0 s i n λ S l S a . (These are the Hamiltonian limits of the dilute face weights, and the resulting matrix is the one used in all our numerics.) In the anyon-chain reading h j h_j h j is simply a sum of projectors on the fusion channels of the two site anyons (the dilute generalisation of the golden-chain Hamiltonian H = − ∑ j P j , j + 1 1 H=-\sum_j P^1_{j,j+1} H = − ∑ j P j , j + 1 1 of [17 ] , in which every site carries a w w w ); the channel projectors are introduced in § 3.4 , where they are needed to build the deformation. The first three are diagonal (the middle height does not change): they are the energies of an empty site, of an isolated wall, and of two adjacent walls. The last three are off-diagonal: a wall hops across the middle bond; a pair of walls is created from, or annihilated into, an empty stretch (this needs l = r l=r l = r , since the two new walls must carry opposite height steps); and the Temperley–Lieb move (with the diagonal part just mentioned), in which two adjacent walls with l = r l=r l = r are replaced by another pair with the same end heights. The square-root factors in the last two lines are what makes h j h_j h j Hermitian in the height basis and are the only place where the quantum dimensions S a S_a S a enter; they are inherited from the loop-model origin of the weights, where a closed domain wall (a loop) carries a fugacity d w d_w d w . With periodic boundary conditions the low-energy spectrum of H 0 H_0 H 0 is that of M ( p , p + 1 ) M(p,p+1) M ( p , p + 1 ) , with c = 7 / 10 c=7/10 c = 7/10 , v = 1.778 v=1.778 v = 1.778 at p = 4 p=4 p = 4 and c = 4 / 5 c=4/5 c = 4/5 , v = 5 / 3 v=5/3 v = 5/3 at p = 5 p=5 p = 5 [16 , 41 ] ; we use these sound velocities to convert all gaps into scaling dimensions, Δ = ( E − E 0 ) L / ( 2 π v ) \Delta=(E-E_0)L/(2\pi v) Δ = ( E − E 0 ) L / ( 2 π v ) [42 , 43 , 44 ] .
Figure 5. The six face-weight moves of (13) , each drawn as a “before → after” configuration of the three bonds ( l , a , r ) → ( l , a ′ , r ) (l,a,r) \to(l,a',r) ( l , a , r ) → ( l , a ′ , r ) . Boxes are heights, the red box is the middle height after the move, and the dots between boxes show whether the site is a domain wall (w, filled) or empty (open). From left to right and top to bottom: no walls (empty); one wall (it stays); two walls (they stay); a wall hops across the middle; a wall pair is created/annihilated (l = r l=r l = r ); and the Temperley–Lieb generator e e e (l = r l=r l = r , both sites occupied), which is also the operator used to build the deformation of Sec. 3.4.
The topological ( r , 1 ) (r,1) ( r , 1 ) interfaces as braided seams
What we want. A defect line is topological if it can be deformed freely without changing any correlation function [2 , 3 , 4 ] . In radial quantisation on the ring, the defect is a static worldline that pierces the spatial circle at one point, between two sites; on the lattice
“topological” means that this insertion commutes with the Hamiltonian in the sense that the strand can be moved from one site to the next without changing the spectrum, and more strongly that the wrapped version of the same line (the operator that transports it once around the ring) commutes with H H H exactly. We now build the ( r , 1 ) (r,1) ( r , 1 ) lines with this property.
The line as an anyon. In the anyon-chain reading of § 3.1 , the ( r , 1 ) (r,1) ( r , 1 ) Verlinde line is the object of C p \mathcal{C}_p C p with label x = r − 1 x=r-1 x = r − 1 (so x = 1 x=1 x = 1 is w = D ( 2 , 1 ) \mathrm{w}=D_{(2,1)} w = D ( 2 , 1 ) and x = p − 2 x=p-2 x = p − 2 is the invertible line η = D ( p − 1 , 1 ) \eta=D_{(p-1,1)} η = D ( p − 1 , 1 ) ). Inserting the defect at a point of the ring means threading a static anyon x x x through the chain there. In the fusion-tree basis this is done by adding one extra bond c j c_j c j between a j a_j a j and a j + 1 a_{j+1} a j + 1 , whose height records the total charge after fusing x x x into the chain,
c j ∈ a j ⊗ x ( for x = w : ∣ c j − a j ∣ = 1 ) . (14) c_j \in a_j \otimes x \qquad(\text{for }x=\mathrm{w}: \lvert c_j-a_j\rvert=1). \tag*{(14)} c j ∈ a j ⊗ x ( for x = w : ∣ c j − a j ∣ = 1 ) . (14)
So the strand of anyon x x x is one more height variable, subject to a stricter adjacency rule than an ordinary bond (for x = w x=\mathrm{w} x = w the site between a j a_j a j and c j c_j c j is always a domain wall), and the next ordinary site, between c j c_j c j and a j + 1 a_{j+1} a j + 1 , obeys the usual dilute rule (6) . Away from the strand nothing changes.
Face operators across the strand: braiding. The only question is how the Hamiltonian acts near the strand (Figure 6 ). The face operator h j h_j h j of (11) needs the two physical neighbours of the middle bond a j a_j a j , namely a j − 1 a_{j-1} a j − 1 and a j + 1 a_{j+1} a j + 1 , but the bond immediately to the right of a j a_j a j is now the strand c j c_j c j , which is not a physical bond, and h j h_j h j knows nothing about anyon x x x . The face operator must therefore skip c j c_j c j and act on the next physical bond a j + 1 a_{j+1} a j + 1 , and it must do so in a way that keeps the strand consistent with the fusion rules. The resolution is to move the strand out of the way before acting and to put it back afterwards.
Figure 6. How the Hamiltonian acts across the strand. (1) The face operator h j h_j h j needs the physical neighbours ( a j − 1 , a j , a j + 1 ) (a_{j-1},a_j,a_{j+1}) ( a j − 1 , a j , a j + 1 ) of the middle bond a j a_j a j , but the strand bond c j c_j c j sits between a j a_j a j and a j + 1 a_{j+1} a j + 1 . (2) The strand x x x is braided past the anyon y y y on site j + 1 j+1 j + 1 with the unitary B B B of (15) ; the label between a j a_j a j and the strand is now b ′ ∈ a j ⊗ y b' \in a_j \otimes y b ′ ∈ a j ⊗ y , an ordinary-looking bond, and the ordinary face operator acts on the triple ( a j − 1 , a j , b ′ ) (a_{j-1},a_j,b') ( a j − 1 , a j , b ′ ) . (3) The strand is braided back with B − 1 B^{-1} B − 1 . The net effect, h j s e a m = B − 1 h j f a c e B h_j^{\mathrm{seam}}=B^{-1}h_j^{\mathrm{face}}B h j seam = B − 1 h j face B , changes a j a_j a j and the strand height c j c_j c j but leaves a j ± 1 a_{j\pm1} a j ± 1 untouched, exactly as an ordinary h j h_j h j would leave its spectators. Right: the same three steps drawn as world-lines in the skein picture of Figure 21 (b), time upwards: the thick line is the ( 2 , 1 ) (2,1) ( 2 , 1 ) strand, the thin lines are the domain walls of sites j j j and j + 1 j+1 j + 1 (present when the site is occupied); the strand x x x crosses in front of the y y y line (B B B ), the face operator acts on the now adjacent walls of sites j j j and j + 1 j+1 j + 1 (grey box), and the strand crosses back, again in front of y y y (B − 1 B^{-1} B − 1 , the inverse braid), so that the two crossings cancel and the world-lines are not linked. The over/under choice is the orientation of the braid discussed in § 2.3 .
Moving an anyon x x x past the anyon y y y on the neighbouring site (here y = 1 y=1 y = 1 if site j + 1 j+1 j + 1 is empty and y = w y=w y = w if it is occupied) is a braiding, which in the fusion-tree basis is the unitary
B = F a y x R y x ( F a x y ) − 1 , (15) B=F^{a y x}R^{y x}(F^{a x y})^{-1}, \tag*{(15)} B = F a y x R y x ( F a x y ) − 1 , (15)
where F F F is the associator (6 j 6j 6 j -symbol) of the category, R y x R^{y x} R y x the diagonal braiding phase [30 , 45 ] , and the matrices act on the intermediate label—i.e. on the strand height c j c_j c j , which after the braid becomes an ordinary-looking label b ′ ∈ a j ⊗ y b'\in a_j\otimes y b ′ ∈ a j ⊗ y . The seam face operator is then
h j s e a m = B − 1 ( h j f a c e ) B , (16) h_j^{\mathrm{seam}}=B^{-1}\left(h_j^{\mathrm{face}}\right)B, \tag*{(16)} h j seam = B − 1 ( h j face ) B , (16)
where h j f a c e h_j^{\mathrm{face}} h j face is the ordinary operator (11) acting on the triple ( a j − 1 , a j , b ′ ) (a_{j-1},a_j,b') ( a j − 1 , a j , b ′ ) . Because B B B is unitary, h j s e a m h_j^{\mathrm{seam}} h j seam is Hermitian and has the same spectrum as h j f a c e h_j^{\mathrm{face}} h j face ; the strand bond is never the middle bond of any face operator, so the Hamiltonian still has exactly L L L terms, one of which is braided. For p = 4 p=4 p = 4 the ingredients are completely explicit: with x = x ′ = w = σ ′ x=x'=w=\sigma' x = x ′ = w = σ ′ and η \eta η the invertible line,
F w w w w = 1 2 ( 1 1 1 − 1 ) , R 1 w w = e − 7 i π / 8 , R η w w = e + 5 i π / 8 , F^{www}_{w}=\frac{1}{\sqrt{2}}\begin{pmatrix}1&1\\1&-1\end{pmatrix},\qquad R^{ww}_{1}=e^{-7i\pi/8},\qquad R^{ww}_{\eta}=e^{+5i\pi/8}, F w w w w = 2 1 ( 1 1 1 − 1 ) , R 1 w w = e − 7 iπ /8 , R η w w = e + 5 iπ /8 ,
so (15) is a 2 × 2 2 \times2 2 × 2 unitary on the two fusion channels w × w = 1 ⊕ η w \times w = 1 \oplus\eta w × w = 1 ⊕ η . For p = 5 p = 5 p = 5 the same formulas use d w = φ d_w = \varphi d w = φ and R R R -phases that are tenth roots of unity, R 1 w w = e − i 4 π / 5 R^{ww}_{1} = e^{-i4\pi/5} R 1 w w = e − i 4 π /5 , R a d j w w = e + 3 i π / 5 R^{ww}_{\mathrm{adj}} = e^{+3i\pi/5} R adj w w = e + 3 iπ /5 . Nothing else is needed: no extra couplings, no fine tuning.
Which category. Fusion rules do not determine a category, and it is worth stating precisely which data the construction uses, since the seam operator is sensitive to them. The ( r , 1 ) (r,1) ( r , 1 ) lines have the fusion rules of S U ( 2 ) p − 2 SU(2)_{p-2} S U ( 2 ) p − 2 , label a = r − 1 = 2 j a = r - 1 = 2j a = r − 1 = 2 j , but they are the Wilson lines of the chiral modular category of M ( p , p + 1 ) M(p,p+1) M ( p , p + 1 ) restricted to its ( r , 1 ) (r,1) ( r , 1 ) sector, whose twists are θ ( r , 1 ) = e 2 π i h ( r , 1 ) \theta_{(r,1)} = e^{2\pi i h_{(r,1)}} θ ( r , 1 ) = e 2 π i h ( r , 1 ) rather than e 2 π i j ( j + 1 ) / p e^{2\pi i j(j+1)/p} e 2 π ij ( j + 1 ) / p : the two differ by the factor i [ a ] i^{[a]} i [ a ] , [ a ] = a m o d 2 [a] = a \bmod2 [ a ] = a mod 2 , since h ( 2 , 1 ) = 3 4 p + 1 4 h_{(2,1)} = \frac{3}{4p} + \frac{1}{4} h ( 2 , 1 ) = 4 p 3 + 4 1 . Accordingly the F F F -symbols are the quantum 6 j 6j 6 j -symbols of S U ( 2 ) p − 2 SU(2)_{p-2} S U ( 2 ) p − 2 multiplied by the Z 2 \mathbb{Z}_2 Z 2 3-cocycle ( − 1 ) [ a ] [ b ] [ c ] [ d ] (-1)^{[a][b][c][d]} ( − 1 ) [ a ] [ b ] [ c ] [ d ] (which flips the Frobenius–Schur indicator of the half-integer-spin objects: at p = 4 p = 4 p = 4 this is the Ising fusion category, with F w w w w = 1 2 ( 1 1 1 − 1 ) F^{www}_{w} = \frac{1}{\sqrt{2}}\begin{pmatrix}1 & 1\\1 & -1\end{pmatrix} F w w w w = 2 1 ( 1 1 1 − 1 ) up to vertex gauge, and not S U ( 2 ) 2 SU(2)_2 S U ( 2 ) 2 ), expressed in a real orthogonal vertex gauge compatible with the positive height-basis projectors of (24) ; commutation of the wrapped lines (21) with H 0 H_0 H 0 is a consistency check on this common choice, rather than a definition of the gauge; and the R R R -symbols are
R c a b = i [ a ] [ b ] ( − 1 ) ( a + b − c ) / 2 q j c ( j c + 1 ) − j a ( j a + 1 ) − j b ( j b + 1 ) , q = e i π / p , (17) R^{ab}_{c} = i^{[a][b]}(-1)^{(a+b-c)/2}q^{j_c(j_c+1)-j_a(j_a+1)-j_b(j_b+1)}, \qquad q = e^{i\pi/p}, \tag*{(17)} R c ab = i [ a ] [ b ] ( − 1 ) ( a + b − c ) /2 q j c ( j c + 1 ) − j a ( j a + 1 ) − j b ( j b + 1 ) , q = e iπ / p , (17)
which solve the hexagon equations with these F F F -symbols and reproduce the minimal-model twists (θ w = 1 d w ∑ c d c R c w w = e 2 π i h ( 2 , 1 ) \theta_w = \frac{1}{d_w}\sum_c d_c R^{ww}_c = e^{2\pi i h_{(2,1)}} θ w = d w 1 ∑ c d c R c w w = e 2 π i h ( 2 , 1 ) ). The prefactor i [ a ] [ b ] i^{[a][b]} i [ a ] [ b ] is not a convention that can be dropped. The braid (15) is 1 1 1 when one of the two slots is empty and R c w w R^{ww}_c R c w w on the channels when both are occupied, and since the dilute Hamiltonian creates and annihilates pairs of walls the relative phase between these two blocks enters h s e a m = B − 1 h B h^{\mathrm{seam}} = B^{-1}hB h seam = B − 1 h B : with the S U ( 2 ) p − 2 SU(2)_{p-2} S U ( 2 ) p − 2 phases − q − 3 / 2 -q^{-3/2} − q − 3/2 , q 1 / 2 q^{1/2} q 1/2 in place of (17) the one-strand ring at p = 4 p = 4 p = 4 , L = 6 L = 6 L = 6 , has lowest scaled levels 0.298 0.298 0.298 , 0.354 0.354 0.354 , 0.757 0.757 0.757 instead of the D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) values 0.1375 0.1375 0.1375 , 0.4375 0.4375 0.4375 , 0.6375 0.6375 0.6375 , which (17) reproduces to 10 − 3 10^{-3} 1 0 − 3 at this size. In the loop-model language of Section 6 the same data appear as the Kauffman braid A 1 + A − 1 e A1 + A^{-1}e A 1 + A − 1 e with A = i q 1 / 2 A = iq^{1/2} A = i q 1/2 , eq. (47) , which agrees with (17) channel by channel (A = i q 1 / 2 = R a d j w w A = iq^{1/2} = R^{ww}_{\mathrm{adj}} A = i q 1/2 = R adj w w , A + A − 1 d w = − i q − 3 / 2 = R 1 w w A + A^{-1}d_w = -iq^{-3/2} = R^{ww}_{1} A + A − 1 d w = − i q − 3/2 = R 1 w w ). All numerical results of the paper, in both formulations, use this one set of data; explicit evaluations of the associators, braid operators and channel projectors must use the same vertex gauge. The braid phases and twists are given above, while (24) fixes the height-basis normalization of the identity-channel projector.
Why this is topological: the wrapped lines. The construction inherits the consistency of the category. The operator Y x Y_x Y x that transports anyon x x x once around the ring is assembled from the same F F F -symbols, one per slot,
( Y x ) { b } { a } = ∏ k F b k b k − 1 a k x a k − 1 ℓ k , ℓ k = anyon on slot k ( 1 , w , or x at the strand ) , (18) (Y_x)_{\{b\}\{a\}} = \prod_k F^{x a_{k-1} \ell_k}_{b_k b_{k-1} a_k}, \qquad\ell_k = \text{anyon on slot } k \ (1,\ w,\ \text{or } x \text{ at the strand}), \tag*{(18)} ( Y x ) { b } { a } = k ∏ F b k b k − 1 a k x a k − 1 ℓ k , ℓ k = anyon on slot k ( 1 , w , or x at the strand ) , (18)
and it commutes with H 0 H_0 H 0 by the pentagon identity and with every braided term (16) by the hexagon identity [30 , 45 ] ; this is the same mechanism that makes the wrapped lines of [19 , 20 , 46 ] exact lattice symmetries. There is a duality here with the dense RSOS models: there it is the ( 1 , s ) (1,s) ( 1 , s ) lines that can be made exactly topological on the lattice, as centralisers of the affine Temperley–Lieb algebra [21 , 23 ] , while the ( r , 1 ) (r,1) ( r , 1 ) lines are only topological in the scaling limit; in the dilute algebra the roles are exchanged, as noted in [23 ] , and the ( r , 1 ) (r,1) ( r , 1 ) lines are the exact ones. So the wrapped ( r , 1 ) (r,1) ( r , 1 ) lines are exact symmetries of the lattice Hamiltonian with a strand, which is the lattice statement that the strand is a topological defect of the same category. We verify this numerically on the ring used throughout (two ( 2 , 1 ) (2,1) ( 2 , 1 ) strands, Figure 7 ):
Figure 7. Two-interface ring. Heights live on the bonds (black dots mark the sites between them); two ( 2 , 1 ) (2,1) ( 2 , 1 ) strands (red), built with the same braiding, are inserted at antipodal points, and the dashed line is a wrapped Verlinde line Y ( r , 1 ) Y_{(r,1)} Y ( r , 1 ) used to resolve the fusion and Z 2 \mathbb{Z}_2 Z 2 charges. Radial quantisation about a strand turns the ring eigenstates into the defect operators of the interface.
∥ [ Y ( 2 , 1 ) , H ] ∥ ∞ , ∥ [ Y ( p − 1 , 1 ) , H ] ∥ ∞ < 10 − 14 . \lVert[Y_{(2,1)},H]\rVert_{\infty},\ \lVert[Y_{(p-1,1)},H]\rVert_{\infty}<10^{-14}. ∥[ Y ( 2 , 1 ) , H ] ∥ ∞ , ∥[ Y ( p − 1 , 1 ) , H ] ∥ ∞ < 1 0 − 14 .
The commuting lines are also conserved charges : the simultaneous eigenvalues of Y ( 2 , 1 ) Y_{(2,1)} Y ( 2 , 1 ) and Y ( p − 1 , 1 ) = η Y_{(p-1,1)}=\eta Y ( p − 1 , 1 ) = η label every eigenstate by its fusion sector and its Z 2 \mathbb{Z}_2 Z 2 parity, and this labelling survives the deformation of § 3.4 , which commutes with the same lines.
Check against the defect spectrum. At zero deformation the ring with two strands is the lattice version of M ( p , p + 1 ) M(p,p+1) M ( p , p + 1 ) on a cylinder crossed twice by D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) , so its low spectrum must be the space of operators living on the line [2 ] , H D × D = H 1 ⊕ H ( 3 , 1 ) \mathcal{H}_{D\times D}=\mathcal{H}_1\oplus\mathcal{H}_{(3,1)} H D × D = H 1 ⊕ H ( 3 , 1 ) (the untwisted spectrum plus the ( 3 , 1 ) (3,1) ( 3 , 1 ) -twisted sector, since D ( 2 , 1 ) × D ( 2 , 1 ) = D ( 1 , 1 ) + D ( 3 , 1 ) D_{(2,1)}\times D_{(2,1)}=D_{(1,1)}+D_{(3,1)} D ( 2 , 1 ) × D ( 2 , 1 ) = D ( 1 , 1 ) + D ( 3 , 1 ) ). This is what we find; the comparison, level by level and charge by charge, is made in § 2.3 (Table 1). At p = 4 p=4 p = 4 the g g g -function of the undeformed strand also reproduces g D ( 2 , 1 ) = 2 g_{D_{(2,1)}}=\sqrt{2} g D ( 2 , 1 ) = 2 (§ 4 ).
The two-strand ring. All computations use a ring of L L L physical sites with two D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) strands at antipodal points (Figure 7 ), both built with the same braiding B B B . Two strands are needed for the total fusion charge of the ring to contain the identity, so that the vacuum sector is present (D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) is self-conjugate, and D ( 2 , 1 ) × D ( 2 , 1 ) ∋ 1 D_{(2,1)}\times D_{(2,1)}\ni1 D ( 2 , 1 ) × D ( 2 , 1 ) ∋ 1 ). Radial quantisation about one strand then turns the ring eigenstates into the defect operators of the interface, and the two-strand geometry is also what allows the transmission coefficient to be measured from a current–current correlator across a single strand (§ 5.2 ).
Antiunitary reflection symmetry. The reflection P P P of the ring about the midpoint between the two strands exchanges the two sides of each strand and maps the braiding B B B to its complex conjugate, so P H P − 1 = H ∗ PHP^{-1}=H^* P H P − 1 = H ∗ . Writing complex conjugation in the height basis as K K K , the combination P K PK P K is therefore an antiunitary symmetry of H H H . This symmetry constrains the eigenstates but does not provide an additional invariant ± \pm ± quantum number: rephasing a state changes its apparent antiunitary eigenphase. The sector labels used below are the eigenvalues of the unitary or Hermitian symmetry operators, in particular ( Y ( 2 , 1 ) , η ) (Y_{(2,1)},\eta) ( Y ( 2 , 1 ) , η ) . In the real F F F -gauge, Y x Y_x Y x is real and [ P K , Y x ] = [ P , Y x ] K [PK,Y_x]=[P,Y_x]K [ P K , Y x ] = [ P , Y x ] K ; consequently complex conjugation alone cannot turn a failure of P P P to commute with Y x Y_x Y x into commutation. The geometric exchange of holomorphic and antiholomorphic labels must be accompanied by the corresponding change of lasso junction. (A variant of the ring with B ∗ B^* B ∗ at the second seam makes P P P itself a unitary symmetry, with the corresponding change in the wrapped-line junctions; we used it only as a check.)
The wrapped lines on the lattice
The junction compatible with the seam. Section 2.3 showed that a line x x x crossing the strand needs a junction u ∈ Hom ( x ⊗ x ′ , x ′ ⊗ x ) u\in\operatorname{Hom}(x\otimes x',x'\otimes x) u ∈ Hom ( x ⊗ x ′ , x ′ ⊗ x ) , with x ′ x' x ′ the strand label, and that for x = x ′ = w x=x'=\mathrm{w} x = x ′ = w there are two independent ones, the over- and under-crossings B B B and B − 1 B^{-1} B − 1 , of which the deformation preserves the over-crossing (§ 2.4 , Figure 3 ). We now explain which junction the lattice realises and why it is the right one. The wrapped line Y x Y_x Y x is built as follows (Figure 8 ). The ring is a fusion tree: a chain of trivalent vertices, one per slot, at which the site anyon ℓ k \ell_k ℓ k (1 or w \mathrm{w} w , and x ′ x' x ′ at the strand slot) is fused into the running label, the heights a k a_k a k being the labels of the segments between the vertices and the chain closing on itself. The line x x x is a closed loop running parallel to the ring; it is fused into the tree at one point of one segment, which splits that segment into a part still carrying the old label and a part carrying the new label b ∈ a ⊗ x b\in a\otimes x b ∈ a ⊗ x (Figure 8 (a), top: x x x attached to segment 5, with a 5 a_5 a 5 on one side of the attachment point and b 5 b_5 b 5 on the other). The attachment point is then moved along the tree one vertex at a time. Moving it across the vertex of slot k k k is the associator of the category, Figure 8 (b): the fusion tree in which x x x is fused into the label a k a_k a k after the vertex is re-expanded in the basis in which x x x is fused into the label a k − 1 a_{k-1} a k − 1 before it, with amplitude F b k − 1 a k − 1 ℓ k x a k b k F^{x a_k b_k}_{b_{k-1}a_{k-1}\ell_k} F b k − 1 a k − 1 ℓ k x a k b k , so that the segment k k k now carries b k b_k b k throughout and the segment k − 1 k-1 k − 1 is split into a k − 1 a_{k-1} a k − 1 and the new label b k − 1 ∈ a k − 1 ⊗ x b_{k-1}\in a_{k-1}\otimes x b k − 1 ∈ a k − 1 ⊗ x (Figure 8 (a), bottom). In the real gauge used here the F F F -symbols are orthogonal matrices, so the same symbol is the amplitude for the move in either direction. After the attachment point has gone once around the ring every segment carries a new label b k b_k b k , the two ends of the x x x line meet again on the segment they started from, and the loop closes; the matrix element of Y x Y_x Y x between the configurations { a } \{a\} { a } and { b } \{b\} { b } is the product of the F F F -symbols of all the moves, with the indices read cyclically, and the closure of the loop is automatic because the sum over the new labels in (21) below includes the label of the starting segment. The strand is simply one more slot of the tree, with site anyon ℓ k = x ′ \ell_k=x' ℓ k = x ′ and outgoing label c j c_j c j , so the x x x line is moved across it by the same F F F -move as across an ordinary site, the strand height changing as c j → c j ′ ∈ c j ⊗ x c_j\to c'_j\in c_j\otimes x c j → c j ′ ∈ c j ⊗ x ; no R R R -matrix appears anywhere in Y x Y_x Y x . This pure-F F F crossing is nevertheless a definite junction and not a third option. In the channel basis of x ⊗ x ′ x \otimes x' x ⊗ x ′ (the channel projectors P D ′ P_{D'} P D ′ of § 3.4 are the junctions u D ′ u_{D'} u D ′ of § 2.3 ) the two half-braidings expand with the R R R -matrix eigenvalues as coefficients, B = ∑ D ′ R D ′ x x ′ u D ′ B = \sum_{D'} R^{xx'}_{D'} u_{D'} B = ∑ D ′ R D ′ x x ′ u D ′ and B − 1 = ∑ D ′ ( R D ′ x x ′ ) − 1 u D ′ B^{-1} = \sum_{D'} (R^{xx'}_{D'})^{-1}u_{D'} B − 1 = ∑ D ′ ( R D ′ x x ′ ) − 1 u D ′ ; for x = x ′ = w x = x' = w x = x ′ = w , with the R R R -symbols (17) and q = e i π / p q = e^{i\pi/p} q = e iπ / p ,
Figure 8. The wrapped line on the lattice. (a) The ring as a fusion tree (eight slots, the strand at slot 3 with its anyon x ′ x' x ′ , red, and strand height c c c ; the chain closes on itself, slot 7 being followed by slot 0), with the symmetry line x x x (blue) fused into one segment: top, x x x attached to segment 5, which carries a 5 a_5 a 5 on one side of the attachment point and b 5 ∈ a 5 ⊗ x b_5 \in a_5 \otimes x b 5 ∈ a 5 ⊗ x on the other, the segments already passed (6, 7) carrying their new labels; bottom, the attachment point moved one vertex to the left, across the vertex of slot 5. (b) The move itself, the associator of the category: the tree with x x x fused into the label after the vertex equals the sum over b k − 1 b_{k-1} b k − 1 of the trees with x x x fused into the label before it, with amplitude F b k b k − 1 a k x a k − 1 ℓ k F^{x a_{k-1} \ell_k}_{b_k b_{k-1} a_k} F b k b k − 1 a k x a k − 1 ℓ k ; at the strand slot ℓ k = x ′ \ell_k = x' ℓ k = x ′ and the label that changes is the strand height, and the move is the same. Repeating the move around the whole ring changes every segment label, a k → b k a_k \to b_k a k → b k (c → c ′ c \to c' c → c ′ at the strand), closes the loop, and gives the cyclic product of F F F -symbols of eq. (21) . (c) The crossing of x x x with the strand decomposed in the fusion channels of x ⊗ x ′ x \otimes x' x ⊗ x ′ , eq. (19) : the pure-F F F move of (b) is one of the two junctions, B B B in the convention of the text, and the other, B − 1 B^{-1} B − 1 , has the conjugate coefficients and is realised by the twist-dressed line (20) .
B = − i q − 3 / 2 u 1 + i q 1 / 2 u a d j , B − 1 = i q 3 / 2 u 1 − i q − 1 / 2 u a d j , (19) B = -i q^{-3/2} u_1 + i q^{1/2} u_{\mathrm{adj}}, \qquad B^{-1} = i q^{3/2} u_1 - i q^{-1/2} u_{\mathrm{adj}}, \tag*{(19)} B = − i q − 3/2 u 1 + i q 1/2 u adj , B − 1 = i q 3/2 u 1 − i q − 1/2 u adj , (19)
so the two junctions differ by the relative phase q ± 2 = e ± 2 π i / p q^{\pm2} = e^{\pm2\pi i/p} q ± 2 = e ± 2 π i / p between the two channels, Figure 8 (c), which is exactly the phase that distinguishes passing in front of the strand from passing behind it. Why does a line built from F F F -moves alone commute with the braided seam term (16) , and which of the two junctions is it? Both questions are answered by noticing what the planar F F F -move does and does not record. Place the trunk of the fusion tree and the line x x x in the plane of Figure 8 ; the strand leaf x ′ x' x ′ joins the trunk at its vertex V V V , and in three dimensions it may come down onto the plane from above or from below. The F F F -move that slides the attachment point of x x x past V V V is a planar operation on the trunk and on x x x , and it is the same matrix in both cases. The two cases differ in the side of x ′ x' x ′ on which x x x passes—under x ′ x' x ′ if the leaf comes from above, over x ′ x' x ′ if it comes from below—and they are related by one full passage of x ′ x' x ′ through the plane around x x x . Figure 9 shows the two situations, before and after the move, flattened into the plane: to compare the move with a junction in Hom ( x ⊗ x ′ , x ′ ⊗ x ) \operatorname{Hom}(x \otimes x',x' \otimes x) Hom ( x ⊗ x ′ , x ′ ⊗ x ) the leaf is brought down to the side of x x x , and it then crosses the trunk segment a k − 1 a_{k-1} a k − 1 either over it (the leaf comes from above the plane) or under it (from below). In the fusion-tree basis the passage from one picture to the other is a number: the two differ by the double braiding of x ′ x' x ′ with the trunk label it is fused to, which is the same as one full turn of the leaf x ′ x' x ′ about its vertex V V V . On the lattice this full turn is the diagonal operator U U U with eigenvalue θ c / ( θ a θ x ′ ) \theta_c/(\theta_a\theta_{x'}) θ c / ( θ a θ x ′ ) on a configuration whose strand vertex carries ( a , x ′ → c ) (a,x' \to c) ( a , x ′ → c ) , with θ i = e 2 π i h i \theta_i = e^{2\pi i h_i} θ i = e 2 π i h i the twist of the ( r , 1 ) (r,1) ( r , 1 ) label i i i . Since Y x Y_x Y x is an operator from the configuration before the move to the one after it, passing to the other reading means applying the full turn on both sides, Y x ′ = U Y x U − 1 Y'_x = UY_xU^{-1} Y x ′ = U Y x U − 1 . For the move of Figure 8 (b), in which x x x passes V V V from right to left, the configuration before the move has x x x to the right of V V V and the trunk at V V V carries ( a k − 1 , x ′ → a k ) (a_{k-1},x' \to a_k) ( a k − 1 , x ′ → a k ) , so U − 1 U^{-1} U − 1 contributes the inverse of θ a k / ( θ a k − 1 θ x ′ ) \theta_{a_k}/(\theta_{a_{k-1}}\theta_{x'}) θ a k / ( θ a k − 1 θ x ′ ) ; the configuration after it has x x x to the left of V V V , the trunk at V V V carries ( b k − 1 , x ′ → b k ) (b_{k-1},x' \to b_k) ( b k − 1 , x ′ → b k ) with b k − 1 ∈ a k − 1 ⊗ x b_{k-1} \in a_{k-1} \otimes x b k − 1 ∈ a k − 1 ⊗ x , and U U U contributes θ b k / ( θ b k − 1 θ x ′ ) \theta_{b_k}/(\theta_{b_{k-1}}\theta_{x'}) θ b k / ( θ b k − 1 θ x ′ ) . The θ x ′ \theta_{x'} θ x ′ cancels and what remains is the ratio θ a k − 1 θ b k / ( θ a k θ b k − 1 ) \theta_{a_{k-1}}\theta_{b_k}/(\theta_{a_k}\theta_{b_{k-1}}) θ a k − 1 θ b k / ( θ a k θ b k − 1 ) , the full turn of x ′ x' x ′ around a k − 1 ⊗ x a_{k-1}\otimes x a k − 1 ⊗ x divided by the full turn around a k − 1 a_{k-1} a k − 1 alone, i.e. the monodromy of x ′ x' x ′ with x x x ; it is precisely what distinguishes the two junctions B B B and B − 1 = B − 2 B B^{-1}=B^{-2}B B − 1 = B − 2 B of (19) . ( U U U itself does not commute with H H H , so the fact that U Y x U − 1 UY_xU^{-1} U Y x U − 1 is again a symmetry is not automatic; it is the lattice statement that both lassos are symmetries of a topological defect.) So the one planar operator corresponds to two three-dimensional situations, and the two lassos of § 2.3 are realised on the lattice by the pure-F F F line Y x Y_x Y x of (21) and by the same line with the F F F -symbol at the strand vertex multiplied by that ratio,
Figure 9. The two junctions as one planar move. The trunk and the line x x x lie in the plane; the strand leaf x ′ x^{\prime} x ′ , attached at the vertex V V V , comes down onto the plane from above (left column) or from below (right column). Top row: the configuration before the move of Figure 8 (b), with x x x attached to the right of V V V and the trunk at V V V carrying a k − 1 a_{k-1} a k − 1 and a k a_k a k ; bottom row: after the move, with x x x attached to the left of V V V and the trunk at V V V carrying b k − 1 ∈ a k − 1 ⊗ x b_{k-1} \in a_{k-1} \otimes x b k − 1 ∈ a k − 1 ⊗ x and b k b_k b k . The planar F F F -move connecting the two rows is the same in both columns; the columns differ in the side of x ′ x^{\prime} x ′ on which x x x passes. Flattened into the plane (the leaf brought down to the side of x x x so that the move can be compared with a junction in Hom ( x ⊗ x ′ , x ′ ⊗ x ) \operatorname{Hom}(x \otimes x^{\prime},x^{\prime} \otimes x) Hom ( x ⊗ x ′ , x ′ ⊗ x ) ), the leaf crosses the trunk over it (left) or under it (right), and in each row the two columns differ by one full turn of x ′ x^{\prime} x ′ about V V V , the double braiding of x ′ x^{\prime} x ′ with the trunk label it is fused to: θ a k / ( θ a k − 1 θ x ′ ) \theta_{a_k}/(\theta_{a_{k-1}}\theta_{x^{\prime}}) θ a k / ( θ a k − 1 θ x ′ ) before the move and θ b k / ( θ b k − 1 θ x ′ ) \theta_{b_k}/(\theta_{b_{k-1}}\theta_{x^{\prime}}) θ b k / ( θ b k − 1 θ x ′ ) after it. The left column is the pure-F F F line Y x Y_x Y x , the right column the dressed line Y x ′ = U Y x U − 1 Y_x^{\prime}=U Y_x U^{-1} Y x ′ = U Y x U − 1 of (20) , whose matrix element carries the ratio of the two factors.
Y x ′ = U Y x U − 1 : F b k b k − 1 a k x a k − 1 x ′ ⟶ F b k b k − 1 a k x a k − 1 x ′ θ a k − 1 θ b k θ a k θ b k − 1 at the strand slot (20) Y'_x = UY_xU^{-1}: \qquad F^{x\,a_{k-1}\,x'}_{b_k\,b_{k-1}\,a_k} \longrightarrow F^{x\,a_{k-1}\,x'}_{b_k\,b_{k-1}\,a_k}\frac{\theta_{a_{k-1}}\theta_{b_k}}{\theta_{a_k}\theta_{b_{k-1}}}\quad\text{at the strand slot} \tag*{(20)} Y x ′ = U Y x U − 1 : F b k b k − 1 a k x a k − 1 x ′ ⟶ F b k b k − 1 a k x a k − 1 x ′ θ a k θ b k − 1 θ a k − 1 θ b k at the strand slot (20)
(with the twists θ \theta θ of the category whose R R R -matrix builds the seam (15) ; for a seam built from B ∗ B^* B ∗ the conjugate twists are used). The planar line commutes with the seam term for either orientation of the braid, since a line that passes on one fixed side of two leaves is indifferent to how those two leaves pass each other (the naturality of the braiding, which in F F F - and R R R -symbols is the hexagon identity); the dressed line commutes with the seam whose orientation matches its dressing. We checked all of this on the undeformed single-strand rings at p = 4 p=4 p = 4 and p = 5 p=5 p = 5 (L = 6 L=6 L = 6 ), to machine precision: Y x Y_x Y x commutes with the seam built from B B B , with the one built from B ∗ B^* B ∗ and with a ring carrying one of each; Y x ′ Y'_x Y x ′ of (20) commutes with H 0 H_0 H 0 and with Y x Y_x Y x , is normal with exactly the Verlinde spectrum ({ ± φ , ± φ − 1 } \{\pm\varphi,\pm\varphi^{-1}\} { ± φ , ± φ − 1 } at p = 5 p=5 p = 5 ), and reads the right label where Y x Y_x Y x reads the left one: on the mirror pair ( 1 , 2 ) ⊗ ( 2 , 2 ) (1,2)\otimes(2,2) ( 1 , 2 ) ⊗ ( 2 , 2 ) , ( 2 , 2 ) ⊗ ( 1 , 2 ) (2,2)\otimes(1,2) ( 2 , 2 ) ⊗ ( 1 , 2 ) at Δ = 3 20 \Delta=\frac{3}{20} Δ = 20 3 , Y ( 2 , 1 ) Y_{(2,1)} Y ( 2 , 1 ) takes the values ( − φ , + φ − 1 ) (-\varphi,+\varphi^{-1}) ( − φ , + φ − 1 ) and Y ( 2 , 1 ) ′ Y'_{(2,1)} Y ( 2 , 1 ) ′ the values ( + φ − 1 , − φ ) (+\varphi^{-1},-\varphi) ( + φ − 1 , − φ ) (for the invertible line η \eta η the two coincide at p = 5 p=5 p = 5 ; at p = 4 p=4 p = 4 , its monodromy with a single σ ′ \sigma' σ ′ strand is θ η ⊗ σ ′ / ( θ η θ σ ′ ) = − 1 \theta_{\eta\otimes\sigma'}/(\theta_\eta\theta_{\sigma'})=-1 θ η ⊗ σ ′ / ( θ η θ σ ′ ) = − 1 , giving Y η ′ = − Y η Y'_\eta=-Y_\eta Y η ′ = − Y η on the single-strand ring and Y η ′ = Y η Y'_\eta=Y_\eta Y η ′ = Y η on the two-strand ring). The lattice thus reproduces § 2.4 literally: on the undeformed ring both lassos are exact symmetries, because the strand is topological, and the deformation decides between them.
Which lasso the deformation preserves, and why. The clean field Φ c l \Phi_{\mathrm{cl}} Φ cl of § 3.4 is built from the channel projectors P 1 P^1 P 1 , P a d j P^{\mathrm{adj}} P adj of the pair formed by the strand and its neighbouring domain wall, eq. (24) . These projectors are planar objects of the fusion tree: they are written in the same basis as Y x Y_x Y x , with the strand leaf treated exactly like a site leaf, i.e. attached to the trunk from the same side of the plane. The closed wrapped categorical line can be moved past these compatible local fusion morphisms by naturality and the pentagon identity, so [ Y x , Φ c l ] = 0 [Y_x,\Phi_{\mathrm{cl}}]=0 [ Y x , Φ cl ] = 0 holds exactly, in the same way as [ Y x , H 0 ] = 0 [Y_x,H_0]=0 [ Y x , H 0 ] = 0 . The dressed line Y x ′ = U Y x U − 1 Y'_x=UY_xU^{-1} Y x ′ = U Y x U − 1 passes the strand on the other side of the plane, and Φ c l \Phi_{\mathrm{cl}} Φ cl ties the strand to its neighbour: [ Y x ′ , Φ c l ] = 0 [Y'_x,\Phi_{\mathrm{cl}}]=0 [ Y x ′ , Φ cl ] = 0 would require [ Y x , U − 1 Φ c l U ] = 0 [Y_x,U^{-1}\Phi_{\mathrm{cl}}U]=0 [ Y x , U − 1 Φ cl U ] = 0 , but the full turn U U U of the strand leaf does not commute with a projector that fuses the strand with another leaf (U − 1 Φ c l U U^{-1}\Phi_{\mathrm{cl}}U U − 1 Φ cl U is not planar), and indeed ∥ [ Y x ′ , Φ c l ] ∥ = O ( 1 ) \lVert[Y'_x,\Phi_{\mathrm{cl}}]\rVert=O(1) ∥[ Y x ′ , Φ cl ]∥ = O ( 1 ) while ∥ [ Y x , Φ c l ] ∥ < 10 − 14 \lVert[Y_x,\Phi_{\mathrm{cl}}]\rVert<10^{-14} ∥[ Y x , Φ cl ]∥ < 1 0 − 14 . In the language of § 2.4 this is the statement that the field created by Φ c l \Phi_{\mathrm{cl}} Φ cl is neutral under one lasso and charged under the other: on the two-strand ring at p = 5 p=5 p = 5 (L = 6 L=6 L = 6 and 8 8 8 ) the state Φ c l ∣ v a c ⟩ \Phi_{\mathrm{cl}}|\mathrm{vac}\rangle Φ cl ∣ vac ⟩ (vacuum component removed) is an eigenvector of both lassos to numerical precision at the tested sizes, ∥ ( Y ( 2 , 1 ) − φ ) Φ c l ∣ v a c ⟩ ∥ < 10 − 12 \lVert(Y_{(2,1)}-\varphi)\Phi_{\mathrm{cl}}|\mathrm{vac}\rangle\rVert<10^{-12} ∥( Y ( 2 , 1 ) − φ ) Φ cl ∣ vac ⟩∥ < 1 0 − 12 and ∥ ( Y ( 2 , 1 ) ′ + φ − 1 ) Φ c l ∣ v a c ⟩ ∥ < 10 − 12 \lVert(Y'_{(2,1)}+\varphi^{-1})\Phi_{\mathrm{cl}}|\mathrm{vac}\rangle\rVert<10^{-12} ∥( Y ( 2 , 1 ) ′ + φ − 1 ) Φ cl ∣ vac ⟩∥ < 1 0 − 12 , with eigenvalue + φ +\varphi + φ under Y ( 2 , 1 ) Y_{(2,1)} Y ( 2 , 1 ) , the neutral value, and − φ − 1 -\varphi^{-1} − φ − 1 under Y ( 2 , 1 ) ′ Y'_{(2,1)} Y ( 2 , 1 ) ′ , the value of the right label ( 3 , 3 ) ≡ ( 2 , 3 ) (3,3)\equiv(2,3) ( 3 , 3 ) ≡ ( 2 , 3 ) , i.e. it is the driving doublet member ( 1 , 3 ) ⊗ ( 3 , 3 ) (1,3)\otimes(3,3) ( 1 , 3 ) ⊗ ( 3 , 3 ) of Table 1 read by its left label under Y x Y_x Y x and by its right label under Y x ′ Y'_x Y x ′ . The deformation therefore preserves Y x Y_x Y x and breaks Y x ′ Y'_x Y x ′ , exactly as (3) requires of the junction that reads the deformation-neutral label. Which of the two sides of the plane is called the holomorphic one is a convention, the same convention that names the lattice deformation ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) \phi_{(1,3)}\otimes\bar{\phi}_{(3,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) rather than its mirror image; with it, the pure-F F F line is the over-crossing lasso of § 2.4 , the one that survives the deformation, which is why [ Y x , H ] = 0 [Y_x,H]=0 [ Y x , H ] = 0 holds exactly (numerically to 10 − 14 10^{-14} 1 0 − 14 ) at every coupling, for Φ c l \Phi_{\mathrm{cl}} Φ cl as well as for H 0 H_0 H 0 , and why the charges it assigns are those of the left label (Table 1 below).
The lattice symmetry operators. Explicitly, the wrapped line x x x is the operator on height configurations
Y x ∣ a 0 , a 1 , … ⟩ = ∑ { b } ∏ k F b k b k − 1 a k x a k − 1 ℓ k ∣ b 0 , b 1 , … ⟩ , b k ∈ a k ⊗ x , ℓ k = { 1 site k empty w site k occupied x ′ strand slot (21) Y_x |a_0,a_1,\ldots\rangle= \sum_{\{b\}} \prod_k F^{x\ a_{k-1}\ \ell_k}_{b_k\ b_{k-1}\ a_k} |b_0,b_1,\ldots\rangle,\qquad b_k \in a_k \otimes x,\qquad\ell_k =
\begin{cases}
1 & \text{site } k \text{ empty}\\
w & \text{site } k \text{ occupied}\\
x^{\prime} & \text{strand slot}
\end{cases} \tag*{(21)} Y x ∣ a 0 , a 1 , … ⟩ = { b } ∑ k ∏ F b k b k − 1 a k x a k − 1 ℓ k ∣ b 0 , b 1 , … ⟩ , b k ∈ a k ⊗ x , ℓ k = ⎩ ⎨ ⎧ 1 w x ′ site k empty site k occupied strand slot (21)
the product running over all slots of the ring, strand slots included; it is obtained by placing the x x x line parallel to the ring and moving it across the tree one vertex at a time, each vertex contributing the associator F F F . For the invertible line η = Y ( p − 1 , 1 ) \eta=Y_{(p-1,1)} η = Y ( p − 1 , 1 ) the fusion a ⊗ η a \otimes\eta a ⊗ η has a single channel, b = p − a b=p-a b = p − a , and Y η Y_\eta Y η is exactly the height reflection a k ↦ p − a k a_k \mapsto p-a_k a k ↦ p − a k on every bond (strand bond included): all the F F F -symbols in (21) multiply to + 1 +1 + 1 on every admissible configuration, so Y η Y_\eta Y η is a pure permutation matrix. These are the operators used throughout the numerics: the ED charges and the η \eta η -block reduction of § 3.5 use (21) directly, and the DMRG uses its matrix-product form. For the exact diagonalisation these are sparse matrices in the admissible basis; the charges of the eigenstates are read off from ⟨ ψ ∣ Y ( 2 , 1 ) ∣ ψ ⟩ \langle\psi|Y_{(2,1)}|\psi\rangle ⟨ ψ ∣ Y ( 2 , 1 ) ∣ ψ ⟩ and ⟨ ψ ∣ η ∣ ψ ⟩ \langle\psi|\eta|\psi\rangle ⟨ ψ ∣ η ∣ ψ ⟩ after diagonalising Y Y Y inside each degenerate multiplet of H H H (the lines commute with H and with each other, so this yields a simultaneous eigenbasis), and η \eta η , being a permutation of the basis, is also used to block-diagonalise H H H (§ 3.5 ). For the DMRG the same operators are written as matrix-product operators: η \eta η has bond dimension 1 (a product of local reflections), and Y ( 2 , 1 ) Y_{(2,1)} Y ( 2 , 1 ) has bond dimension ( p − 1 ) 2 (p-1)^2 ( p − 1 ) 2 , the bond carrying the pair ( b k − 1 , a k − 1 ) (b_{k-1},a_{k-1}) ( b k − 1 , a k − 1 ) needed by the next F F F -symbol, with the ring closure handled by summing over the boundary value of the bond index; they are used both to project on a charge sector (by adding a multiple of them to H H H ) and to verify the charges of the converged states.
Check on the undeformed ring. Table 1 tests (4) on the undeformed two-strand ring at p = 5 p=5 p = 5 : each of the ten lowest levels of H 1 ⊕ H ( 3 , 1 ) \mathcal{H}_1\oplus\mathcal{H}_{(3,1)} H 1 ⊕ H ( 3 , 1 ) lands in the sector predicted from its left label—including the two members of the driving doublet, ( 1 , 3 ) ⊗ ( 3 , 3 ) (1,3)\otimes(3,3) ( 1 , 3 ) ⊗ ( 3 , 3 ) in a + a+ a + and ( 3 , 3 ) ⊗ ( 1 , 3 ) (3,3)\otimes(1,3) ( 3 , 3 ) ⊗ ( 1 , 3 ) in c + c+ c + , which have different left labels and are therefore distinguished by the preserved lasso, as they must be for the deformation to couple to only one of them.
Table 1. The undeformed two-strand ring at p = 5 p=5 p = 5 (K = λ = 0 K=\lambda=0 K = λ = 0 , L = 10 L=10 L = 10 ): lowest levels, their identification in H 1 ⊕ H ( 3 , 1 ) \mathcal{H}_1\oplus\mathcal{H}_{(3,1)} H 1 ⊕ H ( 3 , 1 ) , the charges predicted by (4) from the left label, and the sector in which the level is found numerically (a + = ( φ , + ) a+=(\varphi,+) a + = ( φ , + ) , b − = ( φ − 1 , − ) b-=(\varphi^{-1},-) b − = ( φ − 1 , − ) , c + = ( − φ − 1 , + ) c+=(-\varphi^{-1},+) c + = ( − φ − 1 , + ) , d − = ( − φ , − ) d-=(-\varphi,-) d − = ( − φ , − ) in the notation of § 5.3 ). Weights of M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) : h 1 , 2 = 1 8 h_{1,2}=\frac{1}{8} h 1 , 2 = 8 1 , h 1 , 3 = 2 3 h_{1,3}=\frac{2}{3} h 1 , 3 = 3 2 , h 2 , 1 = 1 5 h_{2,1}=\frac{1}{5} h 2 , 1 = 5 1 , h 2 , 2 = 1 40 h_{2,2}=\frac{1}{40} h 2 , 2 = 40 1 , h 2 , 3 = 1 15 h_{2,3}=\frac{1}{15} h 2 , 3 = 15 1 , h 3 , 2 = 21 40 h_{3,2}=\frac{21}{40} h 3 , 2 = 40 21 , h 3 , 3 = 1 15 h_{3,3}=\frac{1}{15} h 3 , 3 = 15 1 , h 4 , 2 = 13 8 h_{4,2}=\frac{13}{8} h 4 , 2 = 8 13 .
A single crossing. On a ring crossed once by D D D the lattice Y x Y_x Y x is a single lasso on the twisted sector H D \mathcal{H}_D H D , and the same over-crossing rule applies. We checked this on the undeformed single-strand rings: at p = 5 p=5 p = 5 all twelve lowest levels of H D ( 2 , 1 ) \mathcal{H}_{D(2,1)} H D ( 2 , 1 ) (Δ = 2 15 , 3 20 , 2 5 , 11 20 , 11 15 , … \Delta=\frac{2}{15},\frac{3}{20},\frac{2}{5},\frac{11}{20},\frac{11}{15},\ldots Δ = 15 2 , 20 3 , 5 2 , 20 11 , 15 11 , … , e.g. ( 1 , 2 ) ⊗ ( 2 , 2 ) (1,2)\otimes(2,2) ( 1 , 2 ) ⊗ ( 2 , 2 ) and ( 2 , 2 ) ⊗ ( 1 , 2 ) (2,2)\otimes(1,2) ( 2 , 2 ) ⊗ ( 1 , 2 ) both at 3 20 \frac{3}{20} 20 3 ) carry exactly the charges (4) of their left label; at p = 4 p=4 p = 4 the same holds for both lines with the C 4 C_4 C 4 F F F -symbols and the junction normalization used here. A change of sign in the normalization of the ( 2 , 1 ) (2,1) ( 2 , 1 ) crossing junction flips Y ( 2 , 1 ) Y_{(2,1)} Y ( 2 , 1 ) on a ring with a single crossing and cancels on a ring with an even number of crossings; this junction normalization is distinct from the categorical distinction between C 4 C_4 C 4 and S U ( 2 ) 2 SU(2)_2 S U ( 2 ) 2 . The single-strand spectra of § 4 are thus labelled by the Verlinde values of the left label: the Y ( 2 , 1 ) = − 2 Y_{(2,1)}=-\sqrt{2} Y ( 2 , 1 ) = − 2 sector holds the fields with left label ( 1 , 2 ) (1,2) ( 1 , 2 ) , the + 2 +\sqrt{2} + 2 sector those with left label ( 1 , 1 ) (1,1) ( 1 , 1 ) and ( 1 , 3 ) (1,3) ( 1 , 3 ) , and Y ( 2 , 1 ) = 0 Y_{(2,1)}=0 Y ( 2 , 1 ) = 0 those with left label ( 2 , s ) (2,s) ( 2 , s ) .
The clean deformation operator
What the continuum asks for. The KRW flow [6 ] is triggered by a local perturbation on the defect. The field that drives it is the non-chiral defect field ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) \phi_{(1,3)}\otimes\bar{\phi}_{(3,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) of D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) , with ( h , h ˉ ) = ( h 1 , 3 , h 3 , 3 ) (h,\bar{h})=(h_{1,3},h_{3,3}) ( h , h ˉ ) = ( h 1 , 3 , h 3 , 3 ) and scaling dimension
Δ ϕ = 1 − ε , ε = 2 ( p − 1 ) p ( p + 1 ) ( Δ ϕ = 0.7 at p = 4 , 11 15 at p = 5 ) , (22) \Delta_{\phi}=1-\varepsilon,\qquad\varepsilon=\frac{2(p-1)}{p(p+1)}\quad(\Delta_{\phi}=0.7\text{ at }p=4,\ \frac{11}{15}\text{ at }p=5), \tag*{(22)} Δ ϕ = 1 − ε , ε = p ( p + 1 ) 2 ( p − 1 ) ( Δ ϕ = 0.7 at p = 4 , 15 11 at p = 5 ) , (22)
so it is relevant; it carries spin h 1 , 3 − h 3 , 3 h_{1,3}-h_{3,3} h 1 , 3 − h 3 , 3 (1 2 \frac{1}{2} 2 1 at p = 4 p=4 p = 4 , 3 5 \frac{3}{5} 5 3 at p = 5 p=5 p = 5 ), which is allowed for a field living on a line. Its mirror image ϕ ( 3 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) \phi_{(3,3)}\otimes\bar{\phi}_{(1,3)} ϕ ( 3 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) , of the same dimension and opposite spin, is also a defect field of D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) , but it sits in a different charge sector of the wrapped ( 2 , 1 ) (2,1) ( 2 , 1 ) line (§ 3.2 ): on the undeformed ring the two appear as one level at Δ = 11 15 \Delta=\frac{11}{15} Δ = 15 11 in the sector α + \alpha+ α + and one in c + c+ c + (in the notation of § 5.3 ). Among the nonidentity defect primaries in the symmetry sector of the identity, ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) \phi_{(1,3)}\otimes\bar{\phi}_{(3,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) is the only relevant one. A lattice operator that commutes with the wrapped lines, as ours does, connects the vacuum only to that sector, so the lattice flow is KRW’s single-field perturbation of D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) ; it does not preserve the exchange of the two sides of the defect, in agreement with the fact that the KRW trajectory from D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) to C C C lies off the symmetric diagonal λ l = λ r \lambda_l=\lambda_r λ l = λ r of the D ( 1 , 2 ) D_{(1,2)} D ( 1 , 2 ) flows (only the endpoints, D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) and C C C , are symmetric). On the lattice we therefore seek a local operator at the strand with a nonzero overlap with this primary. Its remaining nonidentity relevant content is excluded by the preserved symmetry; identity and irrelevant contributions are allowed and do not obstruct the intended scaling limit. The identity shifts the defect energy, while irrelevant fields produce corrections to scaling, the leading symmetry-allowed primary being ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) \phi_{(1,3)}\otimes\bar{\phi}_{(1,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) of dimension 2 h 1 , 3 2h_{1,3} 2 h 1 , 3 . Symmetry determines which fields may occur, but does not guarantee a nonzero overlap with the driving field. The overlaps below establish that overlap for our chosen operator and motivate a combination that suppresses the displayed admixtures.
Two w’s next to each other. Fix a strand x = w x=w x = w inserted after site j j j . Reading from left to right the relevant heights are (Figure 11 )
Figure 11. The heights around a strand. The strand is an extra bond c j c_j c j (red box) inserted between a j a_j a j and a j + 1 a_{j+1} a j + 1 ; it is always a domain wall (∣ c j − a j ∣ = 1 |c_j-a_j|=1 ∣ c j − a j ∣ = 1 ), so the site between a j a_j a j and c j c_j c j always carries a w w w . The two neighbouring sites may or may not be occupied, depending on the heights. e L , n L e_{\mathrm{L}},n_{\mathrm{L}} e L , n L act on the left pair (site j j j , strand), whose fusion channel is read from the end heights ( a j − 1 , c j ) (a_{j-1},c_j) ( a j − 1 , c j ) ; e R , n R e_{\mathrm{R}},n_{\mathrm{R}} e R , n R act on the right pair (strand, site j + 1 j+1 j + 1 ) with end heights ( a j , a j + 1 ) (a_j,a_{j+1}) ( a j , a j + 1 ) .
a j − 1 → site j a j → strand c j → site j + 1 a j + 1 . a_{j-1}\xrightarrow{\text{site }j}a_j\xrightarrow{\text{strand}}c_j\xrightarrow{\text{site }j+1}a_{j+1}. a j − 1 site j a j strand c j site j + 1 a j + 1 .
The strand bond a j a_j a j –c j c_j c j is always a domain wall, i.e. always carries a w (the defect anyon itself). Site j j j carries a w iff a j − 1 ≠ a j a_{j-1}\ne a_j a j − 1 = a j , site j + 1 j+1 j + 1 iff c j ≠ a j + 1 c_j\ne a_{j+1} c j = a j + 1 . Whenever a neighbouring site is occupied there are two w’s side by side—the site’s domain wall and the strand—and two w’s fuse in two ways,
w ⊗ w = 1 ⊕ adj , (23) w\otimes w=1\oplus\operatorname{adj}, \tag*{(23)} w ⊗ w = 1 ⊕ adj , (23)
where “adj” is the object of C p C_p C p with the adjoint (spin-1) fusion label, the ε \varepsilon ε -like label (η = ( 3 , 1 ) \eta=(3,1) η = ( 3 , 1 ) at p = 4 p=4 p = 4 ; τ F i b = ( 3 , 1 ) \tau_{\mathrm{Fib}}=(3,1) τ Fib = ( 3 , 1 ) at p = 5 p=5 p = 5 ). In the fusion-tree basis the channel is read off from the heights at the two ends of the pair: for the left pair (site j j j , strand) the ends are ( a j − 1 , c j ) (a_{j-1},c_j) ( a j − 1 , c j ) , for the right pair (strand, site j + 1 j+1 j + 1 ) they are ( a j , a j + 1 ) (a_j,a_{j+1}) ( a j , a j + 1 ) . A path l → m → r l\to m\to r l → m → r with ∣ m − l ∣ = ∣ r − m ∣ = 1 \lvert m-l\rvert=\lvert r-m\rvert=1 ∣ m − l ∣ = ∣ r − m ∣ = 1 has r ∈ l ⊗ w ⊗ w r\in l\otimes w\otimes w r ∈ l ⊗ w ⊗ w , and
r = l ± 2 ⇒ channel adj only , r = l ⇒ channels 1 and adj both possible . r=l\pm2\ \Rightarrow\ \text{channel adj only},\qquad r=l\ \Rightarrow\ \text{channels 1 and adj both possible}. r = l ± 2 ⇒ channel adj only , r = l ⇒ channels 1 and adj both possible .
On the r = l r=l r = l block the identity channel is a single state, the quantum-dimension-weighted superposition of the middle height. This defines the two channel projectors P 1 P^1 P 1 and $P^{\mathrm{adj}}】【。 on a triple ( l , m , r ) (l,m,r) ( l , m , r ) :
P 1 ∣ l , m , r ⟩ = δ l , r ∑ m ′ S m S m ′ d w S l ∣ l , m ′ , r ⟩ , P a d j = [ l ≠ m and m ≠ r ] − P 1 (24) P^{1}|l,m,r\rangle=\delta_{l,r}\sum_{m'}\frac{\sqrt{S_mS_{m'}}}{d_wS_l}|l,m',r\rangle,\qquad P^{\mathrm{adj}}=[l\ne m\ \text{and}\ m\ne r]-P^{1} \tag*{(24)} P 1 ∣ l , m , r ⟩ = δ l , r m ′ ∑ d w S l S m S m ′ ∣ l , m ′ , r ⟩ , P adj = [ l = m and m = r ] − P 1 (24)
where the sum runs over heights m ′ m' m ′ adjacent to both l l l and r r r , S a S_a S a and d w d_w d w are as in (9) , and [ ⋯ ] [\cdots] [ ⋯ ] is the indicator that both sites of the triple are domain walls (i.e. that there really are two w’s to fuse). The identity ∑ m a d j l S m = d w S l \sum_{m\ \mathrm{adj}\ l}S_m=d_wS_l ∑ m adj l S m = d w S l gives ( P 1 ) 2 = P 1 (P^{1})^2=P^{1} ( P 1 ) 2 = P 1 ; P 1 P a d j = 0 P^{1}P^{\mathrm{adj}}=0 P 1 P adj = 0 ; and P 1 + P a d j P^{1}+P^{\mathrm{adj}} P 1 + P adj is the projector on the whole two-w subspace. They are the height-basis images of the fusion-tree projectors ∣ w w ; 1 ⟩ ⟨ w w ; 1 ∣ |\mathrm{ww};1\rangle\langle\mathrm{ww};1| ∣ ww ; 1 ⟩ ⟨ ww ; 1∣ and ∣ w w ; a d j ⟩ ⟨ w w ; a d j ∣ |\mathrm{ww};\mathrm{adj}\rangle\langle\mathrm{ww};\mathrm{adj}| ∣ ww ; adj ⟩ ⟨ ww ; adj ∣ , drawn in Figure 10 . In the diagrammatic calculus the projector on the channel c c c is the diagram in which the two w’s are fused into a c c c line and split again, normalised by d c / θ ( w , w , c ) d_c/\theta(\mathrm{w},\mathrm{w},c) d c / θ ( w , w , c ) with θ \theta θ the theta-net of the vertex. For c = 1 c=1 c = 1 the connecting line is invisible and the diagram is the cap–cup divided by the loop value, P 1 = e / d w P^{1}=e/d_w P 1 = e / d w , whose matrix elements are the products of the two F F F -symbols F l ; m 1 l w w = S m / ( d w S l ) F^{l\ \mathrm{w}\ \mathrm{w}}_{l;m\ 1}=\sqrt{S_m/(d_wS_l)} F l ; m 1 l w w = S m / ( d w S l ) of (24) ; for c = a d j c=\mathrm{adj} c = adj the cap and the cup are joined by an adjoint line. The sum of the two projectors is the identity on two strands, i.e. the occupation n n n of (25) below, which is how P a d j P^{\mathrm{adj}} P adj is computed in practice. They commute with the wrapped lines Y ( r , 1 ) Y_{(r,1)} Y ( r , 1 ) .
Figure 10. The fusion-channel projectors of (24) as fusion diagrams (w lines in red). The projector on the channel c c c fuses the two w’s (the domain wall on the site and the strand) into a c c c line and splits them again, with the normalisation d c / θ ( w , w , c ) d_c/\theta(\mathrm{w},\mathrm{w},c) d c / θ ( w , w , c ) that makes the diagram idempotent. For c = 1 c=1 c = 1 the channel line is trivial and P 1 P^{1} P 1 is the cap–cup diagram divided by the loop value d w d_w d w , i.e. the Temperley–Lieb generator e e e up to normalisation; in the height basis the cap and the cup each contribute one F F F -symbol, F l ; m 1 l w w = S m / ( d w S l ) F^{l\ \mathrm{w}\ \mathrm{w}}_{l;m\ 1}=\sqrt{S_m/(d_wS_l)} F l ; m 1 l w w = S m / ( d w S l ) , which gives the weights S m S m ′ / ( d w S l ) \sqrt{S_mS_{m'}}/(d_wS_l) S m S m ′ / ( d w S l ) of (24) . For c = a d j c=\mathrm{adj} c = adj the cap and the cup are joined by an adjoint line (green). Their sum is the identity on two parallel strands, the occupation n n n of (25) .
The two natural building blocks. Two familiar local operators are linear combinations of these projectors: the Temperley–Lieb generator e e e [47 ] , which is d w d_w d w times the identity-channel projector, and the occupation n n n , which is 1 exactly when the site next to the strand is occupied,
e = d w P 1 , n = P 1 + P a d j . (25) e=d_wP^{1},\qquad n=P^{1}+P^{\mathrm{adj}}. \tag*{(25)} e = d w P 1 , n = P 1 + P adj . (25)
Each can be placed on either side of the strand. Written out on the four heights ( a j − 1 , a j , c j , a j + 1 ) (a_{j-1},a_j,c_j,a_{j+1}) ( a j − 1 , a j , c j , a j + 1 ) ,
n L = [ a j − 1 ≠ a j ] , n R = [ c j ≠ a j + 1 ] , n_{\mathrm{L}}=\left[a_{j-1}\ne a_j\right],\qquad n_{\mathrm{R}}=\left[c_j\ne a_{j+1}\right], n L = [ a j − 1 = a j ] , n R = [ c j = a j + 1 ] ,
e L ∣ a j − 1 , a j , c j , a j + 1 ⟩ = δ a j − 1 , c j [ a j − 1 ≠ a j ] ∑ m ′ ≠ a j − 1 S a j S m ′ S a j − 1 ∣ a j − 1 , m ′ , c j , a j + 1 ⟩ , e_{\mathrm{L}}\lvert a_{j-1},a_j,c_j,a_{j+1}\rangle=\delta_{a_{j-1},c_j}\left[a_{j-1}\ne a_j\right]\sum_{m'\ne a_{j-1}}\frac{\sqrt{S_{a_j}S_{m'}}}{S_{a_{j-1}}}\lvert a_{j-1},m',c_j,a_{j+1}\rangle, e L ∣ a j − 1 , a j , c j , a j + 1 ⟩ = δ a j − 1 , c j [ a j − 1 = a j ] m ′ = a j − 1 ∑ S a j − 1 S a j S m ′ ∣ a j − 1 , m ′ , c j , a j + 1 ⟩ ,
e R ∣ a j − 1 , a j , c j , a j + 1 ⟩ = δ a j , a j + 1 [ c j ≠ a j + 1 ] ∑ m ′ ≠ a j S c j S m ′ S a j ∣ a j − 1 , a j , m ′ , a j + 1 ⟩ (26) e_{\mathrm{R}}\lvert a_{j-1},a_j,c_j,a_{j+1}\rangle=\delta_{a_j,a_{j+1}}\left[c_j\ne a_{j+1}\right]\sum_{m'\ne a_j}\frac{\sqrt{S_{c_j}S_{m'}}}{S_{a_j}}\lvert a_{j-1},a_j,m',a_{j+1}\rangle \tag*{(26)} e R ∣ a j − 1 , a j , c j , a j + 1 ⟩ = δ a j , a j + 1 [ c j = a j + 1 ] m ′ = a j ∑ S a j S c j S m ′ ∣ a j − 1 , a j , m ′ , a j + 1 ⟩ (26)
with m ′ m' m ′ restricted to heights adjacent to the two fixed ends. Thus e L e_{\mathrm{L}} e L acts on the triple (bond, site-j j j wall, strand) and replaces the ordinary height a j a_j a j ; e R e_{\mathrm{R}} e R acts on the triple (strand, site-( j + 1 ) (j+1) ( j + 1 ) wall, bond) and replaces the strand height c j c_j c j ; n L , n R n_{\mathrm{L}},n_{\mathrm{R}} n L , n R are diagonal. The matrix of e L , e R e_{\mathrm{L}},e_{\mathrm{R}} e L , e R is the TL move w e w_e w e of (13) and Figure 5 , normalised so that e 2 = d w e e^2=d_w e e 2 = d w e : it is the same Temperley–Lieb generator that builds the bulk energy operator, now straddling the strand. None of these four operators braids anything—they act on the strand bond c j c_j c j exactly as on an ordinary bond—which is why they are local and Hermitian without any B B B -conjugation, and why they commute with the wrapped lines.
Measuring the field content. Which continuum fields does a given lattice operator O O O at the strand contain? Radial quantisation answers this directly: the low-lying eigenstates ∣ Ψ ⟩ \lvert\Psi\rangle ∣ Ψ ⟩ of the undeformed two-strand ring are the defect operators of D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) , so the matrix elements ⟨ Ψ ∣ O ∣ v a c ⟩ \langle\Psi|O|\mathrm{vac}\rangle ⟨ Ψ∣ O ∣ vac ⟩ are, up to normalisation, the coefficients of the corresponding fields in the scaling limit of O O O . Table 2 gives the absolute values of these overlaps for the natural candidates at p = 4 p=4 p = 4 , where the defect spectrum is completely known, and Table 3 repeats the exercise at p = 5 p=5 p = 5 . The tables display magnitudes; linear combinations and cancellation tests use the signed or complex matrix elements before taking absolute values. The result is clear-cut: n L + n R n_{\mathrm{L}}+n_{\mathrm{R}} n L + n R contains no driving field at all—only the identity (a constant defect tension) and the irrelevant ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) \phi_{(1,3)} \otimes\bar{\phi}_{(1,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) ; while e L + e R e_{\mathrm{L}}+e_{\mathrm{R}} e L + e R does contain the driving field but drags along the same two admixtures.
Defect-field content of the candidate lattice operators on the p = 4 p=4 p = 4 D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) ring (∣ ⟨ Ψ ∣ O ∣ v a c ⟩ ∣ |\langle\Psi|\mathcal{O}|\mathrm{vac}\rangle| ∣ ⟨ Ψ∣ O ∣ vac ⟩ ∣ , L = 10 L=10 L = 10 ; absolute overlap magnitudes). Only Φ c l \Phi_{\mathrm{cl}} Φ cl has vanishing identity and ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) \phi_{(1,3)}\otimes\bar{\phi}_{(1,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) components; these overlaps vanish to numerical precision at the tested sizes (see text).
Table 2.
Absolute overlap magnitudes, as in Table 2 , at p = 5 p=5 p = 5 (D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) two-strand ring, L = 10 L=10 L = 10 , neutral sector). The identity and the irrelevant field enter e L + e R e_{\mathrm{L}}+e_{\mathrm{R}} e L + e R and n L + n R n_{\mathrm{L}}+n_{\mathrm{R}} n L + n R in the same ratio 1.62 = φ = d w 1.62=\varphi=d_{\mathrm{w}} 1.62 = φ = d w (at p = 4 p=4 p = 4 , Table 2 , the ratio is 1.41 = 2 = d w 1.41=\sqrt{2}=d_{\mathrm{w}} 1.41 = 2 = d w ), so the combination d w ( e L + e R ) − ( n L + n R ) d_{\mathrm{w}}(e_{\mathrm{L}}+e_{\mathrm{R}})-(n_{\mathrm{L}}+n_{\mathrm{R}}) d w ( e L + e R ) − ( n L + n R ) with d w = φ d_{\mathrm{w}}=\varphi d w = φ removes both; the projector difference 2 d w ( e L + e R ) − ( n L + n R ) \frac{2}{d_{\mathrm{w}}}(e_{\mathrm{L}}+e_{\mathrm{R}})-(n_{\mathrm{L}}+n_{\mathrm{R}}) d w 2 ( e L + e R ) − ( n L + n R ) does not.
Table 3.
The clean combination. Since both contaminations sit in the same two operators, and in the same ratio (Table 2 and Table 3 ), a single combination of the two channel projectors removes both. It is the traceless combination, the projectors weighted by their inverse quantum dimensions,
P 1 d 1 − P a d j d a d j ∝ d a d j P 1 − P a d j = d w e − n , (27) \frac{P^{1}}{d_{1}}-\frac{P^{\mathrm{adj}}}{d_{\mathrm{adj}}}\propto d_{\mathrm{adj}}P^{1}-P^{\mathrm{adj}}=d_{\mathrm{w}}e-n, \tag*{(27)} d 1 P 1 − d adj P adj ∝ d adj P 1 − P adj = d w e − n , (27)
where we used (25) and d w 2 = 1 + d a d j d_{\mathrm{w}}^{2}=1+d_{\mathrm{adj}} d w 2 = 1 + d adj (from w ⊗ w = 1 ⊕ a d j \mathbf{w}\otimes\mathbf{w}=\mathbf{1}\oplus\mathrm{adj} w ⊗ w = 1 ⊕ adj ). Summing the two sides of the strand,
Φ c l = ( d a d j P 1 − P a d j ) L + ( d a d j P 1 − P a d j ) R = d w ( e L + e R ) − ( n L + n R ) (28) \Phi_{\mathrm{cl}}=\left(d_{\mathrm{adj}}P^{1}-P^{\mathrm{adj}}\right)_{\mathrm{L}}+\left(d_{\mathrm{adj}}P^{1}-P^{\mathrm{adj}}\right)_{\mathrm{R}}=d_{\mathrm{w}}(e_{\mathrm{L}}+e_{\mathrm{R}})-(n_{\mathrm{L}}+n_{\mathrm{R}}) \tag*{(28)} Φ cl = ( d adj P 1 − P adj ) L + ( d adj P 1 − P adj ) R = d w ( e L + e R ) − ( n L + n R ) (28)
At p = 4 p=4 p = 4 the adjoint is the invertible η = ( 3 , 1 ) \eta=(3,1) η = ( 3 , 1 ) , d a d j = 1 d_{\mathrm{adj}}=1 d adj = 1 and d w = 2 = 2 / d w d_{\mathrm{w}}=\sqrt{2}=2/d_{\mathrm{w}} d w = 2 = 2/ d w , so (28) is simply the projector difference P 1 − P a d j = 2 d w ( e L + e R ) − ( n L + n R ) P^{1}-P^{\mathrm{adj}}=\frac{2}{d_{\mathrm{w}}}(e_{\mathrm{L}}+e_{\mathrm{R}})-(n_{\mathrm{L}}+n_{\mathrm{R}}) P 1 − P adj = d w 2 ( e L + e R ) − ( n L + n R ) ; at p = 5 p=5 p = 5 the adjoint is the Fibonacci object τ F i b = ( 3 , 1 ) \tau_{\mathrm{Fib}}=(3,1) τ Fib = ( 3 , 1 ) with d a d j = φ d_{\mathrm{adj}}=\varphi d adj = φ and the clean field is φ P 1 − P a d j = φ ( e L + e R ) − ( n L + n R ) \varphi P^{1}-P^{\mathrm{adj}}=\varphi(e_{\mathrm{L}}+e_{\mathrm{R}})-(n_{\mathrm{L}}+n_{\mathrm{R}}) φ P 1 − P adj = φ ( e L + e R ) − ( n L + n R ) , which is not the projector difference. The bottom rows of Table 2 and Table 3 show what this buys on the undeformed D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) ring, at both values of p p p : among the three lowest primaries of the neutral sector, Φ c l \Phi_{\mathrm{cl}} Φ cl connects the vacuum to the driving field only, its identity and ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) \phi_{(1,3)}\otimes\bar{\phi}_{(1,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) components vanish, whereas at p = 5 p=5 p = 5 the naive difference P 1 − P a d j P^{1}-P^{\mathrm{adj}} P 1 − P adj retains both (Table 3 ). The vanishing of the identity component is not an accident of the table: in the topological ground state the two channels of the strand and its neighbouring wall are populated in the ratio of their quantum dimensions, ⟨ n ⟩ = d w ⟨ e ⟩ \langle n\rangle=d_{\mathrm{w}}\langle e\rangle ⟨ n ⟩ = d w ⟨ e ⟩ , so that ⟨ v a c ∣ Φ c l ∣ v a c ⟩ = 0 \langle\mathrm{vac}|\Phi_{\mathrm{cl}}|\mathrm{vac}\rangle=0 ⟨ vac ∣ Φ cl ∣ vac ⟩ = 0 holds to machine precision (10 − 15 10^{-15} 1 0 − 15 ) at every size we have checked, L ≤ 14 L\leq14 L ≤ 14 , at both p p p . This numerical observation is consistent with the interpretation of the identity component as a quantum trace, which vanishes for the channel asymmetry d a d j P 1 − P a d j d_{\mathrm{adj}}P^{1}-P^{\mathrm{adj}} d adj P 1 − P adj since d a d j d 1 − d a d j = 0 d_{\mathrm{adj}}d_{1}-d_{\mathrm{adj}}=0 d adj d 1 − d adj = 0 . We have not attempted a general proof relating this categorical trace to the interacting lattice ground-state expectation value at arbitrary size. Two qualifications are in order. First, like any lattice operator, Φ c l \Phi_{\mathrm{cl}} Φ cl also has components on the descendants of the driving field and on higher primaries of the same lasso sector (the state Φ c l ∣ v a c ⟩ \Phi_{\mathrm{cl}}|\mathrm{vac}\rangle Φ cl ∣ vac ⟩ is an eigenvector of both lassos to numerical precision at the tested sizes, § 3.3 , which fixes the sector but not a single scaling field); these are irrelevant corrections, and what the construction needs is only that the relevant content of the perturbation be the driving field alone. Second, the symmetry argument of § 2.4 says which fields a symmetric lattice operator may contain, not that every symmetric operator contains the driving field; that Φ c l \Phi_{\mathrm{cl}} Φ cl does, with the amplitudes of the tables, is a numerical fact. The meaning of (28) is transparent. On each side of the strand it asks whether the defect anyon and its neighbouring domain wall fuse to the identity or to the adjoint, weighting the answer so that the trivial (channel-blind) part cancels: it is a Kondo-type coupling [7 , 48 ] that distinguishes the two screening channels of the impurity anyon and realises the driving perturbation up to irrelevant corrections. The main structural result is that this channel asymmetry realises the KRW perturbation in the scaling limit, with the unwanted identity and leading irrelevant-primary overlaps vanishing to numerical precision at the tested sizes.
The deformed Hamiltonian. The model studied in the rest of the paper is
H = H 0 + ∑ s t r a n d s [ K Φ c l + λ ( n L + n R ) ] . (29) H=H_0+\sum_{\mathrm{strands}}\left[K\Phi_{\mathrm{cl}}+\lambda(n_{\mathrm{L}}+n_{\mathrm{R}})\right]. \tag*{(29)} H = H 0 + strands ∑ [ K Φ cl + λ ( n L + n R ) ] . (29)
with the two seams built as in § 3.2 . The coupling K K K drives the KRW flow; the second coupling λ \lambda λ multiplies the channel-blind occupation, which contains only the identity and irrelevant fields and is used as a counterterm to tune the finite-size spectrum onto the fixed point (the “improved point” of § 4 , § 5 ). Both terms commute with the wrapped lines, so every eigenstate of (29) keeps its ( Y ( 2 , 1 ) , η ) (Y_{(2,1)},\eta) ( Y ( 2 , 1 ) , η ) labels along the whole flow.
Symmetry-resolved numerics
The admissible Hilbert space grows quickly with L L L (about 2.7 × 10 6 2.7\times10^6 2.7 × 1 0 6 states at p = 4 p=4 p = 4 , L = 16 L=16 L = 16 , and 1.9 × 10 6 1.9\times10^6 1.9 × 1 0 6 at p = 5 p=5 p = 5 , L = 14 L=14 L = 14 , with two strands), so the exact symmetries are used to block-diagonalise the Hamiltonian. The wrapped invertible line Y ( p − 1 , 1 ) = η Y_{(p-1,1)}=\eta Y ( p − 1 , 1 ) = η acts on the height basis as the reflection a ↦ p − a a\mapsto p-a a ↦ p − a , i.e. as a permutation of basis states, so projecting on its η = ± 1 \eta=\pm1 η = ± 1 eigenspaces halves the problem at no cost. At p = 4 p=4 p = 4 , inside η = + 1 \eta=+1 η = + 1 the line Y ( 2 , 1 ) Y_{(2,1)} Y ( 2 , 1 ) satisfies Y ( 2 , 1 ) 2 = 1 + η = 2 Y_{(2,1)}^2=1+\eta=2 Y ( 2 , 1 ) 2 = 1 + η = 2 , hence Y ( 2 , 1 ) / 2 Y_{(2,1)}/\sqrt{2} Y ( 2 , 1 ) / 2 is a second involution that halves the block again. Each block is diagonalised separately (dense for small blocks, Lanczos otherwise) and the charges ( Y ( 2 , 1 ) , η ) (Y_{(2,1)},\eta) ( Y ( 2 , 1 ) , η ) are read back on the eigenvectors, which yields a fully sector-resolved spectrum at a fraction of the cost. The midpoint reflection combined with complex conjugation supplies an antiunitary consistency condition, rather than an additional sector label. For larger sizes (L = 24 , 32 L = 24, 32 L = 24 , 32 at p = 5 p = 5 p = 5 ) we use DMRG [49 , 50 ] on the same Hamiltonian, written as an MPO on the ring, isolating the charge sectors by adding to H H H a multiple of the wrapped lines aligned with the target charges.
Observables
Before turning to the two flows we define the quantities that will be used to characterise an interface, and how each is measured on the ring of § 3 . Throughout, L L L denotes the number of physical sites of the ring, which is also its linear size in the scaling limit.
Scaling dimensions. In radial quantisation about one strand, the eigenstates of H H H are the defect operators of the interface and their energies give the scaling dimensions [42 , 43 , 44 ]
Δ = ( E − E 0 ) L 2 π v , (30) \Delta= \frac{(E-E_0)L}{2\pi v}, \tag*{(30)} Δ = 2 π v ( E − E 0 ) L , (30)
with E 0 E_0 E 0 the ground-state energy of the ring and v v v the sound velocity of § 3.1 . For a topological line D D D the spectrum of the two-strand ring is H D × D = ⨁ k N D D k H k \mathcal{H}_{D\times D} = \bigoplus_k N^k_{DD}\mathcal{H}_k H D × D = ⨁ k N D D k H k (§ 2.3 ), whose content is fixed by the fusion rules; for a non-topological conformal interface it is not known a priori and is one of the outputs of the computation. At a topological endpoint D × D D\times D D × D contains the identity, so the ground state is the vacuum and (30) measures dimensions from it.
Charge sectors. The wrapped lines Y ( 2 , 1 ) Y_{(2,1)} Y ( 2 , 1 ) and η = Y ( p − 1 , 1 ) \eta= Y_{(p-1,1)} η = Y ( p − 1 , 1 ) commute with H H H along the whole flow (§ 3.2 ), so every level carries a pair of quantum numbers. Their eigenvalues on a state of the ( r ′ , s ′ ) (r',s') ( r ′ , s ′ ) type are the Verlinde ratios S ( r , 1 ) , ( r ′ , s ′ ) / S ( 1 , 1 ) , ( r ′ , s ′ ) S_{(r,1),(r',s')}/S_{(1,1),(r',s')} S ( r , 1 ) , ( r ′ , s ′ ) / S ( 1 , 1 ) , ( r ′ , s ′ ) , where for a defect field ϕ ( r ′ , s ) ⊗ ϕ ˉ ( r ′ ′ , s ) \phi_{(r',s)} \otimes\bar{\phi}_{(r'',s)} ϕ ( r ′ , s ) ⊗ ϕ ˉ ( r ′′ , s ) the left label ( r ′ , s ) (r',s) ( r ′ , s ) is the relevant one; the values we need are listed with the sectors in § 4 and § 5 . At p = 4 p = 4 p = 4 , Y ( 2 , 1 ) ∈ { + 2 , 0 , − 2 } Y_{(2,1)} \in\{+\sqrt{2}, 0, -\sqrt{2}\} Y ( 2 , 1 ) ∈ { + 2 , 0 , − 2 } and η = ± 1 \eta= \pm1 η = ± 1 ; at p = 5 p = 5 p = 5 , Y ( 2 , 1 ) = Y τ F i b η Y_{(2,1)} = Y_{\mathrm{\tau Fib}}\eta Y ( 2 , 1 ) = Y τ Fib η with Y τ F i b ∈ { ϕ , − ϕ − 1 } Y_{\mathrm{\tau Fib}} \in\{\phi, -\phi^{-1}\} Y τ Fib ∈ { ϕ , − ϕ − 1 } on the Fibonacci-neutral and Fibonacci-charged states respectively, and η = ± 1 \eta= \pm1 η = ± 1 .
The g g g -function. The g g g -function of an interface is its “ground-state degeneracy” g = ⟨ 0 ∣ D ∣ 0 ⟩ g = \langle0|D|0\rangle g = ⟨ 0∣ D ∣0 ⟩ [26 ] , the overlap of the interface with the vacuum on the cylinder; it equals the quantum dimension for a topological line (g η = 1 g_\eta= 1 g η = 1 , g D ( 2 , 1 ) = 2 g_{D_{(2,1)}} = \sqrt{2} g D ( 2 , 1 ) = 2 at p = 4 p = 4 p = 4 and ϕ \phi ϕ at p = 5 p = 5 p = 5 ) and decreases along any RG flow of the interface [26 , 51 , 52 ] . We extract it from the entanglement entropy, using the fact that the boundary or impurity entropy appears as an additive constant in the entanglement entropy of an interval that contains the defect [53 , 54 , 55 , 56 ] . Let S D ( ℓ ) S_D(\ell) S D ( ℓ ) be the von Neumann entropy of an interval of ℓ = 2 k + 1 \ell= 2k + 1 ℓ = 2 k + 1 bonds centred on a strand (the strand bond included) in the ground state of the ring, and S b u l k ( ℓ ) S_{\mathrm{bulk}}(\ell) S bulk ( ℓ ) the same quantity for an interval of ℓ \ell ℓ bonds in the pure chain. For an interval containing the impurity in its interior the ℓ \ell ℓ -dependence cancels in the difference and
ln g = S D ( ℓ ) − S b u l k ( ℓ ) (31) \ln g = S_D(\ell) - S_{\mathrm{bulk}}(\ell) \tag*{(31)} ln g = S D ( ℓ ) − S bulk ( ℓ ) (31)
up to corrections that vanish as ℓ , L → ∞ \ell,L \to\infty ℓ , L → ∞ . The method is calibrated on the two exactly known topological cases at p = 4 p = 4 p = 4 (§ 4 ).
The transmission coefficient. A conformal interface is characterised, beyond g g g , by how much energy it lets through. Following Quella, Runkel and Watts [57 ] (see also [58 , 59 ] ), the transmission coefficient T T T is defined by the two-point function of the stress tensor across the interface, ⟨ T ( z ) T ( w ) ⟩ \langle T(z)T(w)\rangle ⟨ T ( z ) T ( w )⟩ with z z z and w w w on opposite sides, normalised by its bulk value: T = 1 T=1 T = 1 for a topological line (which is transparent) and T = 0 T=0 T = 0 for a factorising interface (a pair of boundaries). On the lattice we use the energy current
J j = i [ h j , h j + 1 ] , (32) J_j = i [h_j,h_{j+1}], \tag*{(32)} J j = i [ h j , h j + 1 ] , (32)
the time derivative of the local energy; being the divergence of a conserved density, its scaling limit is proportional to T − T ˉ T-\bar{T} T − T ˉ with no admixture of the identity or of scalar primaries, and its bulk two-point function decays as 1 / x 4 1/x^4 1/ x 4 . Each J j J_j J j acts on a block of four consecutive sites. The measurement is made on the two-strand ring as follows. We take the ground state of the ring and evaluate the connected correlator ⟨ J i J j ⟩ c = ⟨ J i J j ⟩ − ⟨ J i ⟩ ⟨ J j ⟩ \langle J_iJ_j\rangle_c=\langle J_iJ_j\rangle-\langle J_i\rangle\langle J_j\rangle ⟨ J i J j ⟩ c = ⟨ J i J j ⟩ − ⟨ J i ⟩ ⟨ J j ⟩ for every pair of blocks with i i i in one half of the ring and j j j in the other, so that the pair straddles the strand at site L − 1 L-1 L − 1 ; blocks that overlap either strand are excluded. The same correlator is evaluated in the ground state of the pure ring of the same length, as a function of the separation d d d of the two blocks, and averaged over the position of the pair. The estimator is the ratio
T ( d ) = ⟨ J i J j ⟩ r i n g w i t h s t r a n d s ⟨ J i J j ⟩ p u r e r i n g , d = separation of the blocks i , j , (33) T(d) = \frac{\langle J_iJ_j\rangle^{\mathrm{ring\ with\ strands}}}{\langle J_iJ_j\rangle^{\mathrm{pure\ ring}}}, \qquad d=\text{separation of the blocks }i,j, \tag*{(33)} T ( d ) = ⟨ J i J j ⟩ pure ring ⟨ J i J j ⟩ ring with strands , d = separation of the blocks i , j , (33)
which divides out the non-universal normalisation of J J J , the 1 / d 4 1/d^4 1/ d 4 decay and the curvature of the ring. Several pairs share the same d d d and differ in where the two blocks sit relative to the two strands; in the continuum the ratio depends only on d d d , on the lattice the pairs scatter at the percent level. We therefore quote, for each d d d , the mean over all admissible pairs together with the standard error of that mean, σ / n p a i r s \sigma/\sqrt{n_{\mathrm{pairs}}} σ / n pairs . Two systematic effects have to be kept in mind. At small d d d one of the two blocks is necessarily close to the strand and the ratio contains the reflected image contribution; T ( d ) T(d) T ( d ) starts well above its asymptotic value and decreases with d d d until it reaches a plateau, which is why the number quoted for a given L L L is the antipodal value d = L / 2 d=L/2 d = L /2 , where both blocks are as far from both strands as the ring allows. In the DMRG calculations the ring is represented by an open matrix-product state whose two ends are joined by the long-range terms of the MPO; the truncation error is largest across that cut, and a block lying on it gives a connected correlator that is low by 1–2% at χ = 200 \chi=200 χ = 200 . Pairs with a block within two sites of the cut are therefore discarded (the cut sits at the centre of one half-ring, away from the strands). Finally the antipodal values are extrapolated linearly in 1 / L 1/L 1/ L .
Lattice β \beta β -function. The KRW flow can be read off directly on the lattice, by the finite-size scaling method introduced by Lüscher, Weisz and Wolff for asymptotically free theories [60 ] : a running coupling is defined by a finite-volume observable (there the mass gap of the O ( 3 ) O(3) O ( 3 ) σ \sigma σ -model in a box of size L L L , building on the small-volume expansion of [61 ] ), and its change from L L L to s L sL s L , the step-scaling function, is the integrated β \beta β -function of that scheme.
Developed into the Schrödinger-functional scheme [62 , 63 ] , this is how the running coupling of QCD is computed non-perturbatively [64 ] . In the statistical-mechanics literature the same idea, applied to an exponent-like finite-size observable, is the phenomenological renormalisation group of Nightingale and of Roomany and Wyld [65 , 66 ] . One needs a dimensionless running coupling, defined by a finite-volume quantity, that vanishes at the UV defect, is monotonic along the flow and takes a definite value at each fixed point. We take it to be a finite-volume Casimir energy of the single-strand ring (one D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) strand, spectrum H D ( 2 , 1 ) \mathcal{H}_{D_{(2,1)}} H D ( 2 , 1 ) of § 2.3 ), which is precisely the type of coupling used in [60 ] : the L L L -dependent part of the energy of one fixed state of the ring. At the UV point the two lowest states of H D ( 2 , 1 ) \mathcal{H}_{D_{(2,1)}} H D ( 2 , 1 ) are the mirror pair ( 1 , 2 ) ⊗ ( 2 , 2 ) (1,2) \otimes(2,2) ( 1 , 2 ) ⊗ ( 2 , 2 ) and ( 2 , 2 ) ⊗ ( 1 , 2 ) (2,2) \otimes(1,2) ( 2 , 2 ) ⊗ ( 1 , 2 ) , degenerate at
Δ k = 1 = h 1 , 2 + h 2 , 2 = ( p − 1 ) 2 + 2 4 p ( p + 1 ) ( 11 80 at p = 4 , 3 20 at p = 5 ) , (34) \Delta_{k=1}=h_{1,2}+h_{2,2}=\frac{(p-1)^2+2}{4p(p+1)}\quad\left(\frac{11}{80}\text{ at }p=4,\ \frac{3}{20}\text{ at }p=5\right), \tag*{(34)} Δ k = 1 = h 1 , 2 + h 2 , 2 = 4 p ( p + 1 ) ( p − 1 ) 2 + 2 ( 80 11 at p = 4 , 20 3 at p = 5 ) , (34)
exchanged by the swap of the two sides of the defect and distinguished by their left labels, hence by their ( Y ( 2 , 1 ) , η ) (Y_{(2,1)},\eta) ( Y ( 2 , 1 ) , η ) sectors. We track the member with left label ( 1 , 2 ) (1,2) ( 1 , 2 ) , on which Y ( 2 , 1 ) Y_{(2,1)} Y ( 2 , 1 ) takes the Verlinde value S ( 2 , 1 ) , ( 1 , 2 ) / S ( 1 , 1 ) , ( 1 , 2 ) = − 2 cos ( π / p ) S_{(2,1),(1,2)}/S_{(1,1),(1,2)}=-2\cos(\pi/p) S ( 2 , 1 ) , ( 1 , 2 ) / S ( 1 , 1 ) , ( 1 , 2 ) = − 2 cos ( π / p ) , so that it sits in the sector ( Y ( 2 , 1 ) , η ) = ( − 2 , + 1 ) (Y_{(2,1)},\eta)=(-\sqrt{2},+1) ( Y ( 2 , 1 ) , η ) = ( − 2 , + 1 ) at p = 4 p=4 p = 4 and in ( − φ , − 1 ) (-\varphi,-1) ( − φ , − 1 ) , the sector d − d- d − of § 5.3 , at $p=5, and call it the k = 1 k=1 k = 1 state (in the loop model of § 6 it is the ground state of the sector whose non-contractible loops carry the weight − 2 cos ( π / p ) -2\cos(\pi/p) − 2 cos ( π / p ) ). At integer p p p it remains the lowest state of its charge sector along the whole flow, so the conserved charges identify it at every K K K ; its mirror partner is the ground state of the ring and barely moves. Let E ( k = 1 ) ( L , K ) E^{(k=1)}(L,K) E ( k = 1 ) ( L , K ) be its energy and δ E ( L , K ) = E ( k = 1 ) ( L , K ) − E ( k = 1 ) ( L , 0 ) \delta E(L,K)=E^{(k=1)}(L,K)-E^{(k=1)}(L,0) δ E ( L , K ) = E ( k = 1 ) ( L , K ) − E ( k = 1 ) ( L , 0 ) . Writing ℓ \ell ℓ for the midpoint of two consecutive ring sizes ℓ ∓ 1 \ell\mp1 ℓ ∓ 1 (the ring is computed for even L L L , so ℓ = 7 , 9 , 11 , … \ell=7,9,11,\ldots ℓ = 7 , 9 , 11 , … ), the running coupling attached to the size ℓ \ell ℓ is the symmetric size difference
u C ( ℓ , K ) = ℓ 2 − 1 4 π v [ δ E ( ℓ − 1 , K ) − δ E ( ℓ + 1 , K ) ] , (35) u_C(\ell,K)=\frac{\ell^2-1}{4\pi v}\left[\delta E(\ell-1,K)-\delta E(\ell+1,K)\right], \tag*{(35)} u C ( ℓ , K ) = 4 π v ℓ 2 − 1 [ δ E ( ℓ − 1 , K ) − δ E ( ℓ + 1 , K ) ] , (35)
a lattice version of − ℓ 2 ∂ ℓ δ E / 2 π v -\ell^2\partial_\ell\delta E/2\pi v − ℓ 2 ∂ ℓ δ E /2 π v evaluated at ℓ \ell ℓ . Since E ( L , K ) = e b u l k L + e d e f ( K ) + 2 π v L ( x − c 12 ) + … E(L,K)=e_{\mathrm{bulk}}L+e_{\mathrm{def}}(K)+\frac{2\pi v}{L}\left(x-\frac{c}{12}\right)+\ldots E ( L , K ) = e bulk L + e def ( K ) + L 2 π v ( x − 12 c ) + … [42 , 43 , 44 ] , subtracting K = 0 K=0 K = 0 removes the bulk energy and the size difference removes the non-universal defect energy, so u C → x ( K ) − x ( 0 ) u_C\to x(K)-x(0) u C → x ( K ) − x ( 0 ) , the change of the scaling dimension of the tracked state. It vanishes at D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) by construction and takes a universal value at a fixed point: on the η \eta η line of p = 4 p=4 p = 4 the state becomes ( 1 10 , 3 20 ) \left(\frac{1}{10},\frac{3}{20}\right) ( 10 1 , 20 3 ) of the Z 2 \mathbb{Z}_2 Z 2 -twisted sector, so u C ∗ = 0.7 − 11 80 = 0.5625 u_C^*=0.7-\frac{11}{80}=0.5625 u C ∗ = 0.7 − 80 11 = 0.5625 , while at a conformal interface its value is one of the outputs. Since u C u_C u C is dimensionless, its logarithmic size derivative at fixed K K K ,
B C ( K ; ℓ 1 , ℓ 2 ) = d u C d ln ℓ ≃ u C ( ℓ 2 , K ) − u C ( ℓ 1 , K ) ln ( ℓ 2 / ℓ 1 ) , (36) B_C(K;\ell_1,\ell_2)=\frac{d u_C}{d\ln\ell}\simeq\frac{u_C(\ell_2,K)-u_C(\ell_1,K)}{\ln(\ell_2/\ell_1)}, \tag*{(36)} B C ( K ; ℓ 1 , ℓ 2 ) = d ln ℓ d u C ≃ ln ( ℓ 2 / ℓ 1 ) u C ( ℓ 2 , K ) − u C ( ℓ 1 , K ) , (36)
is a lattice β \beta β -function for the coordinate u C u_C u C : it vanishes where u C u_C u C is scale-invariant (a fixed point), and its slope at a zero is minus the RG eigenvalue of the perturbation that moves u C u_C u C away from it, so that a zero with negative slope is IR-attractive. In the UV, u C ∝ K ℓ 1 − Δ ϕ u_C\propto K\ell^{1-\Delta_\phi} u C ∝ K ℓ 1 − Δ ϕ to first order, so the curve leaves the origin with slope 1 − Δ ϕ 1-\Delta_\phi 1 − Δ ϕ . Plotting B C B_C B C against the running coordinate A C = 1 2 [ u C ( ℓ 1 ) + u C ( ℓ 2 ) ] A_C=\frac{1}{2}[u_C(\ell_1)+u_C(\ell_2)] A C = 2 1 [ u C ( ℓ 1 ) + u C ( ℓ 2 )] for successive size pairs collapses the data onto a single curve if scaling holds.
Improved points. A lattice Hamiltonian of the family (29) sits at a fixed point only up to irrelevant operators, whose amplitudes depend on ( K , λ ) (K,\lambda) ( K , λ ) and produce 1 / L 1/L 1/ L corrections to every quantity. We call an improved point a choice of ( K , λ ) (K,\lambda) ( K , λ ) at which the leading corrections are minimised, in the spirit of Symanzik improvement [67 ] , as judged by the flatness of the finite-size spectrum in 1 / L 1/L 1/ L . Different improved points are different lattice representatives of the same continuum interface.
Numerical methods. Exact diagonalisation (ED) uses the symmetry-reduced blocks of § 3.5 for L ≤ 14 L \le14 L ≤ 14 (up to 16 for a few quantities). Density-matrix renormalisation group (DMRG) calculations [49 , 50 ] treat the same ring at L = 24 , 32 , 48 L = 24,32,48 L = 24 , 32 , 48 (p = 4 p = 4 p = 4 ) and L = 24 , 32 , 48 L = 24,32,48 L = 24 , 32 , 48 (p = 5 p = 5 p = 5 ) with matrix-product states of bond dimension up to χ = 300 \chi= 300 χ = 300 ; the vacuum and the first excited state of the neutral sector are obtained by adding to H H H a multiple of the wrapped lines that penalises the other charge sectors, and excited states by orthogonalisation against the states already found. The entropies S D S_D S D , S b u l k S_{\mathrm{bulk}} S bulk and the current correlators are evaluated directly on the converged states.
Results at p = 4 p = 4 p = 4 : D ( 2 , 1 ) → η D_{(2,1)} \to\eta D ( 2 , 1 ) → η
For the tricritical Ising model the category is Ising, D ( 2 , 1 ) = σ ′ D_{(2,1)} = \sigma' D ( 2 , 1 ) = σ ′ has g = 2 g = \sqrt{2} g = 2 , and the only candidate endpoints of a flow that preserves the ( r , 1 ) (r,1) ( r , 1 ) symmetry and lowers g g g are the invertible line η \eta η (g = 1 g = 1 g = 1 ) and the identity (g = 1 g = 1 g = 1 as well). KRW do not determine the endpoint at p = 4 p = 4 p = 4 : the last model in which their TCSA shows convincing evidence of the defect C C C is M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) , while in M ( 4 , 5 ) M(4,5) M ( 4 , 5 ) the signature of the fixed point is no longer seen and, in their words, the position of the conformal fixed point may have moved a large distance or may require a higher truncation level ([6 ] , §6.2.4 and § 7 ). They do note that at p = 4 p = 4 p = 4 a factorising defect, D ( 2 , 1 ) F 1 ∣ 1 D_{(2,1)}F^{1|1} D ( 2 , 1 ) F 1∣1 , is an additional candidate with g g g below that of D ( 1 , 2 ) D_{(1,2)} D ( 1 , 2 ) ([6 ] , § 4 ). The invertible line η \eta η was not among the candidates considered by KRW. It is, however, singled out by the symmetry argument of § 2 combined with the g g g -theorem: a flow that preserves the over-crossing lasso and lowers g g g below 2 \sqrt{2} 2 can only end on a defect that commutes with the ( r , 1 ) (r,1) ( r , 1 ) lines and has g ≤ 2 g \le\sqrt{2} g ≤ 2 , which leaves η \eta η , the identity, and the factorising defects, and the lattice decides between them. The result of this section is that the endpoint is η \eta η : the KRW flow of the tricritical Ising model ends on a topological, invertible line rather than on a genuinely conformal defect. (See [11 ] for the conformal defects of the tricritical Ising model.) The bulk chain used here, the dilute A 3 A_3 A 3 model at the integrable critical point of [15 ] , belongs to the same class of Ising-fusion-category chains (labels 1 , σ , ψ 1,\sigma,\psi 1 , σ , ψ , i.e. 1 , w , η 1,w,\eta 1 , w , η in our notation, with the fusion constraint on neighbours) as the two-parameter model of Roychowdhury and Wang [24 ] . The integrable weights (13) place us directly at that point with no tuning.
Calibration of the g g g -function. We first apply (31) to the two undeformed topological rings, with strands η \eta η and σ ′ \sigma' σ ′ respectively:
η = D ( 3 , 1 ) : ln g = − 0.01 ( g = 0.99 ) , σ ′ = D ( 2 , 1 ) : ln g = + 0.335 ( g = 1.40 ) , \eta= D_{(3,1)}:\ \ln g = -0.01\ (g = 0.99), \qquad\sigma' = D_{(2,1)}:\ \ln g = +0.335\ (g = 1.40), η = D ( 3 , 1 ) : ln g = − 0.01 ( g = 0.99 ) , σ ′ = D ( 2 , 1 ) : ln g = + 0.335 ( g = 1.40 ) ,
in agreement with g η = 1 g_{\eta}=1 g η = 1 and g σ ′ = 2 = 1.414 g_{\sigma'}=\sqrt{2}=1.414 g σ ′ = 2 = 1.414 . The residual ± 0.01 \pm0.01 ± 0.01 sets the accuracy of the method at these sizes.
Endpoint of the flow. Driving D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) with the clean field Φ c l \Phi_{\mathrm{cl}} Φ cl and following ln g \ln g ln g as a function of g g g , the entropy drops from ln 2 \ln\sqrt{2} ln 2 and saturates at ln g = − 0.04 \ln g=-0.04 ln g = − 0.04 (g = 1.00 g=1.00 g = 1.00 ) at strong coupling, stable against further increase of K K K : the flow ends on an interface with g = 1 g=1 g = 1 , and the sector-resolved spectrum below identifies it as η \eta η .
Improved point. The endpoint η \eta η is invertible and topological. It corresponds to a conformal phase of (29) . We can tune away the leading finite-size effect by searching for a point with good conformal tower. The couplings are
( K , λ ) η = ( 4.82 , 0.17 ) . (37) (K,\lambda)_{\eta}=(4.82,0.17). \tag*{(37)} ( K , λ ) η = ( 4.82 , 0.17 ) . (37)
g g g -function, transmission and the leading irrelevant dimension. Figure 12 shows the three observables of § 3.6 at the improved point as functions of 1 / L 1/L 1/ L , in the same format as for C 5 C_5 C 5 below, from ED at L = 8 L=8 L = 8 –14 and DMRG at L = 24 L=24 L = 24 , 32 and 48. The DMRG uses the commuting lines to restrict the search to the neutral sector: the ground state and first excited state of H ′ = H − 1.0 Y ( 2 , 1 ) − 0.8 η H'=H-1.0\,Y_{(2,1)}-0.8\,\eta H ′ = H − 1.0 Y ( 2 , 1 ) − 0.8 η , which pushes every non-neutral sector up by at least 1.4, are the vacuum and the leading neutral excitation of H H H (their charges come out as ⟨ Y ( 2 , 1 ) ⟩ = 1.4142 \langle Y_{(2,1)}\rangle=1.4142 ⟨ Y ( 2 , 1 ) ⟩ = 1.4142 , ⟨ η ⟩ = 1.0000 \langle\eta\rangle=1.0000 ⟨ η ⟩ = 1.0000 ). At p = 4 p=4 p = 4 the local Hilbert space dimension is 3, and the Hamiltonian MPO has bond dimension 20, so χ = 160 \chi=160 χ = 160 (truncation error 2 × 10 − 9 2\times10^{-9} 2 × 1 0 − 9 ) is enough. (i) The g g g -function from (31) is flat at g = 0.99 ( 2 ) g=0.99(2) g = 0.99 ( 2 ) for L = 8 L=8 L = 8 –14 (individual values 1.002, 0.975, 0.998, 0.985) and g = 0.9965 g=0.9965 g = 0.9965 , 0.9980, 0.9993 at L = 24 L=24 L = 24 , 32, 48, i.e. at the invertible value g η = 1 g_{\eta}=1 g η = 1 within the ± 0.01 \pm0.01 ± 0.01 accuracy of the calibration,
Figure 12. p = 4 p=4 p = 4 at the improved point ( K , λ ) = ( 4.82 , 0.17 ) (K,\lambda)=(4.82,0.17) ( K , λ ) = ( 4.82 , 0.17 ) : finite-size dependence of the g g g -function (left), the transmission T T T at antipodal separation d = L / 2 d=L/2 d = L /2 (centre) and the dimension of the leading neutral irrelevant operator (right). Filled circles: ED at L = 8 , … , 14 L=8,\ldots,14 L = 8 , … , 14 ; open squares: DMRG at L = 24 L=24 L = 24 , 32, 48. Dotted grey lines: reference values (g η = 1 g_{\eta}=1 g η = 1 and g D ( 2 , 1 ) = 2 g_{D_{(2,1)}}=\sqrt{2} g D ( 2 , 1 ) = 2 ; the topological T = 1 T=1 T = 1 ; Δ ε ′ = 1.2 \Delta_{\varepsilon'}=1.2 Δ ε ′ = 1.2 ). The dashed curve in the centre panel is a power-law guide 1 + c L − 2.7 1+cL^{-2.7} 1 + c L − 2.7 through the points; the dashed line and star on the right are the linear-in-1 / L 1/L 1/ L extrapolation of the L ≥ 12 L\geq12 L ≥ 12 points.
and far from g D ( 2 , 1 ) = 2 g_{D_{(2,1)}}=\sqrt{2} g D ( 2 , 1 ) = 2 . (ii) The transmission (33) at d = L / 2 d=L/2 d = L /2 is 1.54, 1.31, 1.21 for L = 10 L=10 L = 10 , 12, 14 and 1.049, 1.021, 1.002 at L = 24 L=24 L = 24 , 32, 48, above the topological value 1 and decreasing towards it approximately as a power law, T − 1 ∝ L − 2.7 T-1\propto L^{-2.7} T − 1 ∝ L − 2.7 . A value above 1 1 1 is not possible for a conformal interface and is a short-distance artefact. We check that the estimator itself is exact for a topological line: on the undeformed η \eta η and σ ′ \sigma' σ ′ rings at L = 12 L = 12 L = 12 it gives T = 1.000 \mathcal{T} = 1.000 T = 1.000 at every separation. The approach to 1 1 1 from above is therefore consistent with a topological endpoint, in contrast with the p = 5 p = 5 p = 5 case where T \mathcal{T} T approaches a value clearly below 1 1 1 from above (§ 5.2 ). (iii) The leading irrelevant defect field in the neutral sector (Y ( 2 , 1 ) = + 2 Y_{(2,1)} = +\sqrt{2} Y ( 2 , 1 ) = + 2 , η = + 1 \eta= +1 η = + 1 ) has Δ = 1.238 , 1.245 , 1.247 , 1.248 \Delta= 1.238, 1.245, 1.247, 1.248 Δ = 1.238 , 1.245 , 1.247 , 1.248 for L = 8 L = 8 L = 8 –14 14 14 and 1.244 , 1.241 , 1.236 1.244, 1.241, 1.236 1.244 , 1.241 , 1.236 at L = 24 , 32 , 48 L = 24, 32, 48 L = 24 , 32 , 48 , extrapolating to 1.23 1.23 1.23 : it is the image of the driving field along the flow and is identified with the bulk field ϵ ′ \epsilon' ϵ ′ of dimension 1.2 1.2 1.2 , the residual ∼ 3 % \sim3\% ∼ 3% being of the same size as for the other primaries in Figure 13 .
Figure 13. p = 4 p = 4 p = 4 sector-resolved spectrum at the improved point ( K , λ ) = ( 4.82 , 0.17 ) (K,\lambda) = (4.82, 0.17) ( K , λ ) = ( 4.82 , 0.17 ) , Δ \Delta Δ versus 1 / L 1/L 1/ L , one panel per ( Y ( 2 , 1 ) , η ) (Y_{(2,1)},\eta) ( Y ( 2 , 1 ) , η ) sector, showing the three (four) lowest levels of each sector from ED at L = 8 , 10 , 12 , 14 L = 8, 10, 12, 14 L = 8 , 10 , 12 , 14 (circles) and the three lowest from DMRG at L = 24 L = 24 L = 24 (stars); lines connect the same level across sizes. Grey dashed lines: the diagonal tricritical-Ising spectrum expected on the η \eta η line, primaries and their first descendants (T T T , T ‾ \overline{T} T at 2 2 2 above the identity; the L − 1 , L ‾ − 1 L_{-1},\overline{L}_{-1} L − 1 , L − 1 pairs at 1.2 1.2 1.2 above ϵ \epsilon ϵ and at 1.075 1.075 1.075 above σ \sigma σ ; the level-2 2 2 descendants of ϵ \epsilon ϵ at 2.2 2.2 2.2 ). Open diamonds at 1 / L = 0 1/L = 0 1/ L = 0 : linear-in-1 / L 1/L 1/ L extrapolations of the leading levels through the three largest sizes. The primaries converge to within a few percent; the descendant pairs are still split at these sizes and approach their common value from both sides. Every level is consistent with the invertible-η \eta η spectrum, confirming D ( 2 , 1 ) → η D_{(2,1)} \to\eta D ( 2 , 1 ) → η .
Sector-resolved spectrum. At the IR endpoint each strand has become η \eta η , and η × η = 1 \eta\times\eta= 1 η × η = 1 , so the ring spectrum must be the untwisted, diagonal spectrum of the tricritical Ising model itself, with the Verlinde eigenvalues
Y ( 2 , 1 ) = + 2 , η = + 1 : { 1 , ϵ ′ } , Y ( 2 , 1 ) = − 2 , η = + 1 : { ϵ , ϵ ′ ′ } , Y_{(2,1)} = +\sqrt{2},\ \eta= +1:\ \{\mathbf{1},\epsilon'\}, \qquad Y_{(2,1)} = -\sqrt{2},\ \eta= +1:\ \{\epsilon,\epsilon''\}, Y ( 2 , 1 ) = + 2 , η = + 1 : { 1 , ϵ ′ } , Y ( 2 , 1 ) = − 2 , η = + 1 : { ϵ , ϵ ′′ } ,
Y ( 2 , 1 ) = 0 , η = − 1 : { σ , σ ′ } . Y_{(2,1)} = 0,\ \eta= -1:\ \{\sigma,\sigma'\}. Y ( 2 , 1 ) = 0 , η = − 1 : { σ , σ ′ } .
where ϵ , ϵ ′ , ϵ ′ ′ , σ , σ ′ \epsilon,\epsilon',\epsilon'',\sigma,\sigma' ϵ , ϵ ′ , ϵ ′′ , σ , σ ′ are the primaries of dimension Δ = 0.2 , 1.2 , 3 , 0.075 , 0.875 \Delta= 0.2, 1.2, 3, 0.075, 0.875 Δ = 0.2 , 1.2 , 3 , 0.075 , 0.875 . Figure 13 shows the three sectors at the improved point for L = 8 , … , 14 L = 8,\ldots,14 L = 8 , … , 14 (ED) and L = 24 L = 24 L = 24 (DMRG, three states per sector, obtained as the lowest states of H H H shifted by the wrapped lines aligned with the sector charges, as in § 5.3 ). Linear-in-1 / L 1/L 1/ L extrapolation of the leading levels through the three largest sizes gives
Δ → 0.08 ( σ , 0.075 ) , 0.24 ( ϵ , 0.2 ) , 0.87 ( σ ′ , 0.875 ) , 1.24 ( ϵ ′ , 1.2 ) , \Delta\to0.08\ (\sigma,\,0.075), \qquad0.24\ (\epsilon,\,0.2), \qquad0.87\ (\sigma',\,0.875), \qquad1.24\ (\epsilon',\,1.2), Δ → 0.08 ( σ , 0.075 ) , 0.24 ( ϵ , 0.2 ) , 0.87 ( σ ′ , 0.875 ) , 1.24 ( ϵ ′ , 1.2 ) ,
and the descendants behave as they should. Together with g → 1 g \to1 g → 1 this establishes D ( 2 , 1 ) → η D_{(2,1)} \to\eta D ( 2 , 1 ) → η or D ( 2 , 1 ) → 1 D_{(2,1)} \to1 D ( 2 , 1 ) → 1 at p = 4 p = 4 p = 4 . The single strand spectrum will distinguish the two cases.
Lattice β \beta β -function. We now apply the step-scaling construction of § 3.6 : the Casimir coupling u C u_C u C of (35) , built from the k = 1 k=1 k = 1 state of the single-strand ring, which at p = 4 p=4 p = 4 is the lowest state of the sector ( Y ( 2 , 1 ) , η ) = ( − 2 , + 1 ) (Y_{(2,1)},\eta)=(-\sqrt{2},+1) ( Y ( 2 , 1 ) , η ) = ( − 2 , + 1 ) (its mirror partner, the lowest state of ( 0 , − 1 ) (0,-1) ( 0 , − 1 ) , is the ground state of the ring). On the η \eta η line the tracked state becomes ( 1 10 , 3 5 ) (\frac{1}{10},\frac{3}{5}) ( 10 1 , 5 3 ) of the Z 2 \mathbb{Z}_2 Z 2 -twisted sector, so the fixed-point value is u C ∗ = 0.7 − 11 80 = 0.5625 u_C^*=0.7-\frac{11}{80}=0.5625 u C ∗ = 0.7 − 80 11 = 0.5625 . We follow the ray through the improved point ( K , λ ) η = ( 4.82 , 0.17 ) (K,\lambda)_\eta=(4.82,0.17) ( K , λ ) η = ( 4.82 , 0.17 ) , i.e. H = H 0 + K Φ c l + λ ( n L + n R ) H=H_0+K\Phi_{\rm cl}+\lambda(n_L+n_R) H = H 0 + K Φ cl + λ ( n L + n R ) with λ = 0.035 K \lambda=0.035K λ = 0.035 K (the counterterm is kept in fixed proportion to the clean coupling), for couplings 0 ≤ K ≤ 3 0\le K\le3 0 ≤ K ≤ 3 in steps of 0.1 (0.05 for K ≤ 0.2 K\le0.2 K ≤ 0.2 ), with E ( k = 1 ) ( L , K ) E^{(k=1)}(L,K) E ( k = 1 ) ( L , K ) from the symmetry-reduced ED of § 3.5 at L = 6 , 8 , … , 16 L=6,8,\ldots,16 L = 6 , 8 , … , 16 , i.e. ℓ = 7 , 9 , … , 15 \ell=7,9,\ldots,15 ℓ = 7 , 9 , … , 15 and the size pairs ( 7 , 9 ) , … , ( 13 , 15 ) (7,9),\ldots,(13,15) ( 7 , 9 ) , … , ( 13 , 15 ) .
Figure 14. Lattice β \beta β -function at p = 4 p=4 p = 4 along the ray through the improved point, H = H 0 + K Φ c l + λ ( n L + n R ) H=H_0+K\Phi_{\rm cl}+\lambda(n_L+n_R) H = H 0 + K Φ cl + λ ( n L + n R ) with λ = 0.035 K \lambda=0.035K λ = 0.035 K , from the k = 1 k=1 k = 1 state of the single-strand ring (ED, L = 6 , … , 16 L=6,\ldots,16 L = 6 , … , 16 ). Left: B C = d u C / d ln ℓ B_C=du_C/d\ln\ell B C = d u C / d ln ℓ of (36) versus the running coordinate A C = 1 2 [ u C ( ℓ 1 ) + u C ( ℓ 2 ) ] A_C=\frac{1}{2}[u_C(\ell_1)+u_C(\ell_2)] A C = 2 1 [ u C ( ℓ 1 ) + u C ( ℓ 2 )] for the size pairs ℓ = ( 7 , 9 ) , … , ( 13 , 15 ) \ell=(7,9),\ldots,(13,15) ℓ = ( 7 , 9 ) , … , ( 13 , 15 ) , for 0 ≤ K ≤ 3 0\le K\le3 0 ≤ K ≤ 3 . The black square is D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) (u C = 0 u_C=0 u C = 0 ), the dotted line the UV slope 1 − Δ ϕ = 0.3 1-\Delta_\phi=0.3 1 − Δ ϕ = 0.3 , the vertical dashed line the η \eta η value u C ∗ = 0.5625 u_C^*=0.5625 u C ∗ = 0.5625 , and the dash-dotted line the universal IR slope Δ ϵ ′ − 1 = − 0.2 \Delta_{\epsilon'}-1=-0.2 Δ ϵ ′ − 1 = − 0.2 through the zero. Right: the raw coupling u C ( ℓ , K ) u_C(\ell,K) u C ( ℓ , K ) versus K K K for ℓ = 7 , … , 15 \ell=7,\ldots,15 ℓ = 7 , … , 15 ; the dashed line is u C ∗ u_C^* u C ∗ .
Figure 14 shows the result. The collapse of the successive size pairs is excellent—from ( 7 , 9 ) (7,9) ( 7 , 9 ) to ( 13 , 15 ) (13,15) ( 13 , 15 ) the curves lie on top of one another over the whole range—and the curve has the expected structure: it leaves D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) with positive slope, rises to a maximum B C ≃ 0.065 B_C\simeq0.065 B C ≃ 0.065 near K ≃ 1.2 K\simeq1.2 K ≃ 1.2 (A C ≃ 0.3 A_C\simeq0.3 A C ≃ 0.3 ), and crosses zero with negative slope at
K ∗ = 1.900 , 1.892 , 1.900 , 1.912 , K^*=1.900,\ 1.892,\ 1.900,\ 1.912, K ∗ = 1.900 , 1.892 , 1.900 , 1.912 ,
A C ∗ = 0.516 , 0.514 , 0.515 , 0.518 for ℓ = ( 7 , 9 ) , ( 9 , 11 ) , ( 11 , 13 ) , ( 13 , 15 ) , A_C^*=0.516,\ 0.514,\ 0.515,\ 0.518 \qquad\text{for }\ell=(7,9),\ (9,11),\ (11,13),\ (13,15), A C ∗ = 0.516 , 0.514 , 0.515 , 0.518 for ℓ = ( 7 , 9 ) , ( 9 , 11 ) , ( 11 , 13 ) , ( 13 , 15 ) ,
within 8% of the η \eta η value u C ∗ = 0.5625 u_C^*=0.5625 u C ∗ = 0.5625 and drifting slowly towards it; the right panel shows the same thing in the raw coupling, which rises through u C ∗ u_C^* u C ∗ at K ≃ 2.1 K\simeq2.1 K ≃ 2.1 , overshoots to 0.62–0.69 at K = 3 K=3 K = 3 with the overshoot shrinking with ℓ \ell ℓ . The UV slope B C / A C B_C/A_C B C / A C at the smallest couplings is 0.19, 0.21, 0.22, 0.23 for the four pairs and extrapolates linearly in 1 / ℓ 1/\ell 1/ ℓ to 0.28, consistent with the predicted 1 − Δ ϕ = 0.3 1-\Delta_\phi=0.3 1 − Δ ϕ = 0.3 (the first-order coefficient of u C u_C u C carries a lattice correction that decreases with ℓ \ell ℓ ). The slope at the zero, − 0.59 -0.59 − 0.59 to − 0.63 -0.63 − 0.63 , is much steeper than the universal Δ ε ′ − 1 = − 0.2 \Delta_{\varepsilon'}-1=-0.2 Δ ε ′ − 1 = − 0.2 of the leading irrelevant operator on η \eta η : as at p = 5 p=5 p = 5 below, the corrections to scaling at these sizes are dominated by more irrelevant operators, the amplitude of ε ′ \varepsilon' ε ′ being what the improved-point tuning suppresses, and the slow L − 0.2 L^{-0.2} L − 0.2 tail would only become visible at much larger sizes. The β \beta β -function thus reproduces, from the single-strand setup alone, the flow from D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) to η \eta η with the correct UV exponent and, to within the finite-size accuracy, the correct IR value of the Casimir coupling.
Single-strand ring: η \eta η , not the identity. The two-strand spectrum above is the untwisted spectrum H D × D ‾ \mathcal{H}_{D\times\overline{D}} H D × D , and it would be the same, namely the diagonal tricritical-Ising spectrum, for either candidate endpoint with g = 1 g=1 g = 1 : the invertible line η \eta η (η × η = 1 \eta\times\eta=1 η × η = 1 ) and the trivial line (1 × 1 = 1 1\times1=1 1 × 1 = 1 ). A ring with a single strand distinguishes them. In radial quantisation its spectrum is the twisted Hilbert space H D \mathcal{H}_D H D of the endpoint line: for the identity line this is again the untwisted spectrum, with a vacuum at Δ = 0 \Delta=0 Δ = 0 ; for η \eta η it is the Z 2 \mathbb{Z}_2 Z 2 -twisted sector of the tricritical Ising model, H η = ⨁ ( r , s ) V ( r , s ) ⊗ V ‾ ( 4 − r , s ) \mathcal{H}_{\eta}=\bigoplus_{(r,s)}V_{(r,s)}\otimes\overline{V}_{(4-r,s)} H η = ⨁ ( r , s ) V ( r , s ) ⊗ V ( 4 − r , s ) , which has no vacuum and consists of the primaries ( 3 80 , 3 80 ) (\frac{3}{80},\frac{3}{80}) ( 80 3 , 80 3 ) , ( 1 10 , 3 5 ) (\frac{1}{10},\frac{3}{5}) ( 10 1 , 5 3 ) , ( 3 5 , 1 10 ) (\frac{3}{5},\frac{1}{10}) ( 5 3 , 10 1 ) , ( 7 16 , 7 16 ) (\frac{7}{16},\frac{7}{16}) ( 16 7 , 16 7 ) , ( 0 , 3 2 ) (0,\frac{3}{2}) ( 0 , 2 3 ) , ( 3 2 , 0 ) (\frac{3}{2},0) ( 2 3 , 0 ) of dimensions 0.075 0.075 0.075 , 0.7 0.7 0.7 , 0.7 0.7 0.7 , 0.875 0.875 0.875 , 1.5 1.5 1.5 , 1.5 1.5 1.5 (the second, third, fifth and sixth carrying spin ± 1 2 \pm\frac{1}{2} ± 2 1 and ± 3 2 \pm\frac{3}{2} ± 2 3 ) and their descendants. We therefore repeat the computation on a ring with one D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) strand at the same improved point, for L = 8 L=8 L = 8 –14 14 14 by ED and L = 24 L=24 L = 24 by DMRG, measuring energies from E b u l k ( L ) + e D E_{\mathrm{bulk}}(L)+e_D E bulk ( L ) + e D , where e D = 1 2 [ E 0 t w o s t r a n d s ( L ) − E b u l k ( L ) ] e_D=\frac{1}{2}[E^{\mathrm{two\ strands}}_0(L)-E_{\mathrm{bulk}}(L)] e D = 2 1 [ E 0 two strands ( L ) − E bulk ( L )] is the non-universal energy of one deformed strand extracted from the two-strand ring at the same L L L (e D = − 3.5921 e_D=-3.5921 e D = − 3.5921 at L = 14 L=14 L = 14 and − 3.5915 -3.5915 − 3.5915 at L = 24 L=24 L = 24 –48 48 48 , so it is well converged). The wrapped lines still commute with H H H and label the levels. The result (Figure 15 ) is unambiguous: there is no state at Δ = 0 \Delta=0 Δ = 0 , the ground state is ( 3 80 , 3 80 ) (\frac{3}{80},\frac{3}{80}) ( 80 3 , 80 3 ) at Δ = 0.081 → 0.08 \Delta=0.081\to0.08 Δ = 0.081 → 0.08 , and the remaining levels reproduce the twisted spectrum—the spin-1 2 \frac{1}{2} 2 1 pair at 0.69 / 0.67 → 0.68 / 0.67 0.69/0.67\to0.68/0.67 0.69/0.67 → 0.68/0.67 , split by the σ ′ \sigma' σ ′ -line into Y ( 2 , 1 ) = ± 2 Y_{(2,1)}=\pm\sqrt{2} Y ( 2 , 1 ) = ± 2 according to its left label (§ 2.3 ); ( 7 16 , 7 16 ) (\frac{7}{16},\frac{7}{16}) ( 16 7 , 16 7 ) at 0.877 → 0.87 0.877\to0.87 0.877 → 0.87 ; the spin-3 2 \frac{3}{2} 2 3 pair at 1.46 / 1.52 → 1.51 / 1.53 1.46/1.52\to1.51/1.53 1.46/1.52 → 1.51/1.53 ; and descendants at 1.69 / 1.67 1.69/1.67 1.69/1.67 and 1.03 1.03 1.03 . The single-strand ring thus confirms that the IR endpoint of the p = 4 p=4 p = 4 flow is the invertible line η \eta η and not the trivial line.
Figure 15. p = 4 p = 4 p = 4 , ring with a single D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) strand at the improved point ( K , λ ) = ( 4.82 , 0.17 ) (K,\lambda) = (4.82, 0.17) ( K , λ ) = ( 4.82 , 0.17 ) : sector-resolved spectrum versus 1 / L 1/L 1/ L , with Δ = ( E − E b u l k − e D ) L / 2 π v \Delta= (E - E_{\mathrm{bulk}} - e_D)L/2\pi v Δ = ( E − E bulk − e D ) L /2 π v measured from the bulk ground state plus the defect energy of one strand. Circles: ED at L = 8 , … , 14 L = 8,\ldots,14 L = 8 , … , 14 ; stars: DMRG at L = 24 L = 24 L = 24 (three states per sector); open diamonds at 1 / L = 0 1/L = 0 1/ L = 0 : linear extrapolations through the three largest sizes. Grey dashed lines: the η \eta η -twisted sector of the tricritical Ising model, primaries ( h , h ˉ ) (h,\bar{h}) ( h , h ˉ ) and their first descendants. There is no state at Δ = 0 \Delta= 0 Δ = 0 : the endpoint is η \eta η , not the identity line.
Results at p = 5 p=5 p = 5 : D ( 2 , 1 ) → C 5 D_{(2,1)}\to C_5 D ( 2 , 1 ) → C 5
At p = 5 p=5 p = 5 the category is F i b ⊠ Z 2 \mathrm{Fib}\boxtimes\mathbb{Z}_2 Fib ⊠ Z 2 and D ( 2 , 1 ) = τ F i b ⊗ η D_{(2,1)}=\tau_{\mathrm{Fib}}\otimes\eta D ( 2 , 1 ) = τ Fib ⊗ η has g = φ = 1.618 g=\varphi=1.618 g = φ = 1.618 . The invertible line η \eta η is still available as an endpoint with g = 1 g=1 g = 1 , but the flow does not reach it: it stops at a conformal interface C 5 C_5 C 5 with
g C 5 = 1.503 ( 2 ) ( ln g = 0.408 ) , T = 0.665 ( 10 ) , g_{C_5}=1.503(2)\quad(\ln g=0.408),\qquad\mathcal{T}=0.665(10), g C 5 = 1.503 ( 2 ) ( ln g = 0.408 ) , T = 0.665 ( 10 ) ,
to be compared with ln g = ln φ = 0.481 \ln g=\ln\varphi=0.481 ln g = ln φ = 0.481 and T = 1 \mathcal{T}=1 T = 1 for D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) . This section documents the location of C 5 C_5 C 5 , its three defining numbers with their 1 / L 1/L 1/ L dependence, its sector-resolved spectrum, and the lattice β \beta β -function of the flow.
Location of the fixed point
We use two parametrisations of the family (29) , with Φ c l = φ ( e L + e R ) − ( n L + n R ) \Phi_{\mathrm{cl}} = \varphi(e_L + e_R) - (n_L + n_R) Φ cl = φ ( e L + e R ) − ( n L + n R ) the clean field (28) at p = 5 p = 5 p = 5 . The first is the clean ray λ = 0 \lambda= 0 λ = 0 ,
H = H 0 + K Φ c l , (38) H = H_0 + K\Phi_{\mathrm{cl}}, \tag*{(38)} H = H 0 + K Φ cl , (38)
on which a single coupling K K K runs from the UV defect ( K = 0 ) (K = 0) ( K = 0 ) through C 5 C_5 C 5 . On this ray the fixed point is located by the lattice β \beta β -function of the Casimir coupling (§ 3.6 ), which for the clean ray is computed in the loop-model formulation of § 6 (identical to the RSOS chain at p = 5 p = 5 p = 5 , § 6.1 ): its zeros for the successive size pairs converge with geometrically decreasing steps (§7.1) to
K C 5 ∗ = 0.95 ( 5 ) ( λ = 0 ) (39) K_{C_5}^{*} = 0.95(5)\quad(\lambda= 0) \tag*{(39)} K C 5 ∗ = 0.95 ( 5 ) ( λ = 0 ) (39)
The second is a two-parameter tower-improved point,
( K , λ ) = ( 0.84 , 0.34 ) , (40) (K,\lambda) = (0.84,\ 0.34), \tag*{(40)} ( K , λ ) = ( 0.84 , 0.34 ) , (40)
the same clean field with an irrelevant counterterm in the channel-blind occupation, which gives the flattest finite-size spectrum and is used for the spectra of § 5.3 . The two are distinct lattice representatives of the same continuum interface, differing only in irrelevant content, and they give consistent values of g g g , T T T and Δ X \Delta_X Δ X . In the following, the g g g -function, the transmission and the leading irrelevant dimension of § 5.2 , the sector-resolved spectrum and its 1 / L 1/L 1/ L extrapolation of § 5.3 , and the Potts-interface spectrum of § 5.5 —is computed at the tower-improved point (40) , where the corrections to scaling are smallest, with ED for L ≤ 16 L \le16 L ≤ 16 and DMRG at L = 24 L = 24 L = 24 , 32 32 32 , 48 48 48 . The lattice β \beta β -function of § 5.4 requires a one-parameter family that starts at the ultraviolet defect D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) at K = 0 K = 0 K = 0 and passes through the fixed point; it is measured along the ray through the tower-improved point, λ = 0.405 K \lambda= 0.405K λ = 0.405 K , as at p = 4 p = 4 p = 4 ; the clean ray (38) , whose step-scaling curves collapse better and which fixes the location (39) of C 5 C_5 C 5 , is followed in the loop-model formulation of Sections 6 and 7. The quantities that can be compared between the two points—the conserved charges and sector assignments, the leading irrelevant dimension Δ X \Delta_X Δ X and the fixed-point value of the Casimir coupling—agree within the quoted uncertainties.
g-function, transmission and the leading irrelevant dimension
Figure 17 and Table 4 collect the three defining numbers of C 5 C_5 C 5 and their finite-size flow, combining ED (L ≤ 16 L \le16 L ≤ 16 ) with DMRG (L = 24 , 32 L = 24, 32 L = 24 , 32 at χ = 200 \chi= 200 χ = 200 and L = 48 L = 48 L = 48 at χ = 300 \chi= 300 χ = 300 ).
Table 4. g g g -function and transmission of C 5 C_5 C 5 at the tower-improved point ( K , λ ) = ( 0.84 , 0.34 ) (K,\lambda) = (0.84, 0.34) ( K , λ ) = ( 0.84 , 0.34 ) versus system size. ED for L ≤ 16 L \le16 L ≤ 16 ; DMRG at L = 24 , 32 L = 24, 32 L = 24 , 32 (χ = 200 \chi= 200 χ = 200 ) and L = 48 L = 48 L = 48 (χ = 300 \chi= 300 χ = 300 ). ln g \ln g ln g is the entropy difference (31) for the interval of 2 k + 1 2k+1 2 k + 1 bonds centred on the strand, T \mathcal{T} T the current-correlator ratio (33) at antipodal separation d = L / 2 d = L/2 d = L /2 , averaged over the admissible pairs (blocks near the MPS cut excluded), with the standard error of the mean in parentheses where more than one pair contributes. The values below the table are the linear-in-1 / L 1/L 1/ L extrapolations through L ≥ 14 L \ge14 L ≥ 14 (Figure 17 ).
Figure 17. Finite-size flow of the three C 5 C_5 C 5 observables versus 1 / L 1/L 1/ L : the g g g -function (left), the transmission T \mathcal{T} T at antipodal separation d = L / 2 d=L/2 d = L /2 (centre), and the dimension Δ X \Delta_{\mathrm{X}} Δ X of the leading irrelevant operator (right). Filled circles: ED; open squares: DMRG at L = 24 , 32 L=24,32 L = 24 , 32 (χ = 200 \chi=200 χ = 200 ) and L = 48 L=48 L = 48 (χ = 300 \chi=300 χ = 300 ; g g g and T \mathcal{T} T only); dashed line and star: linear-in-1 / L 1/L 1/ L extrapolation. Dotted grey lines mark the reference values (g = φ g=\varphi g = φ and the topological T = 1 \mathcal{T}=1 T = 1 ).
(i) The g g g -function from (31) is already flat by L = 12 L = 12 L = 12 : ln g = 0.4085 \ln g = 0.4085 ln g = 0.4085 , 0.4089 0.4089 0.4089 , 0.4084 0.4084 0.4084 , 0.4076 0.4076 0.4076 at L = 14 L = 14 L = 14 , 24 24 24 , 32 32 32 , 48 48 48 , i.e. g = 1.503 ( 2 ) g = 1.503(2) g = 1.503 ( 2 ) (the residual drift between L = 32 L = 32 L = 32 and 48 48 48 is 10 − 3 10^{-3} 1 0 − 3 in ln g \ln g ln g and is within the χ \chi χ -dependence of the entropy), well below the topological value φ \varphi φ and above the invertible value 1.
(ii) Figure 16 shows how the transmission is measured. The estimator (33) is averaged over all admissible pairs of current blocks at the same separation d d d (left panel, L = 48 L = 48 L = 48 ); the ratio depends on d d d only for d ≲ 8 d \lesssim8 d ≲ 8 , where the image contribution of the block nearest to the strand dominates, and then settles on a size-dependent plateau. The number quoted for each L L L is the antipodal value d = L / 2 d = L/2 d = L /2 : 0.801, 0.762, 0.740, 0.726 from ED (L = 10 L = 10 L = 10 –16) and 0.707, 0.6975(10), 0.685(2) from DMRG (L = 24 , 32 , 48 L = 24, 32, 48 L = 24 , 32 , 48 ). These lie on a straight line in 1 / L 1/L 1/ L (right panel) and extrapolate to T = 0.663 \mathcal{T} = 0.663 T = 0.663 through the DMRG points (0.664 through L ≥ 14 L \ge14 L ≥ 14 ; a 1 / L + 1 / L 2 1/L + 1/L^2 1/ L + 1/ L 2 fit through L ≥ 12 L \ge12 L ≥ 12 gives 0.677). We quote T = 0.665 ( 10 ) \mathcal{T} = 0.665(10) T = 0.665 ( 10 ) , the uncertainty covering the spread between fit forms. This is clearly below 1: C 5 C_5 C 5 is genuinely partially reflecting, not a topological line in disguise.
Figure 16. Measurement of the transmission coefficient of C 5 C_5 C 5 at the tower-improved point. Left: the estimator (33) , i.e. the connected current correlator across the strand divided by the bulk correlator at the same separation, versus the separation d d d of the two four-site current blocks, for L = 48 L = 48 L = 48 (DMRG, χ = 300 \chi= 300 χ = 300 ). Each point is the mean over all admissible pairs of blocks at that d d d (both blocks away from the strands and at least three sites from the cut of the matrix-product state); the standard errors of the means, ≲ 2 × 10 − 3 \lesssim2 \times10^{-3} ≲ 2 × 1 0 − 3 , are smaller than the symbols. The dashed line is the L → ∞ L \to\infty L → ∞ value. Right: the antipodal values d = L / 2 d = L/2 d = L /2 against 1 / L 1/L 1/ L with their standard errors and the linear extrapolation through the three DMRG points (dotted). The star marks the quoted value T = 0.665 ( 10 ) T = 0.665(10) T = 0.665 ( 10 ) .
(iii) The leading irrelevant defect field of C 5 C_5 C 5 , which we call X X X , is the lowest excitation in the fully neutral sector (Y τ F i b = φ Y_{\tau\mathrm{Fib}} = \varphi Y τ Fib = φ , η = + 1 \eta= +1 η = + 1 ; the sector of the identity). It is the image of the driving field along the flow, and its dimension controls the approach to the fixed point. Its gap rises past marginality with L L L and extrapolates to
Δ X = 1.12 ( 1 ) , ϵ X ≡ Δ X − 1 = 0.12 ( 1 ) . (41) \Delta_X = 1.12(1), \qquad\epsilon_X \equiv\Delta_X - 1 = 0.12(1). \tag*{(41)} Δ X = 1.12 ( 1 ) , ϵ X ≡ Δ X − 1 = 0.12 ( 1 ) . (41)
The DMRG points at L = 24 , 32 , 48 L = 24, 32, 48 L = 24 , 32 , 48 are essential here: they anchor the extrapolations far beyond the ED reach and confirm that the small-L L L trends are not finite-size artefacts.
Sector-resolved spectrum and its 1 / L 1/L 1/ L dependence
The four charge sectors of the ring are labelled by the Fibonacci charge (1 1 1 or τ F i b \tau_{\mathrm{Fib}} τ Fib , i.e. Y τ F i b = φ Y_{\tau\mathrm{Fib}} = \varphi Y τ Fib = φ or − φ − 1 -\varphi^{-1} − φ − 1 ) and the Z 2 \mathbb{Z}_2 Z 2 parity η = ± 1 \eta= \pm1 η = ± 1 ; we use the shorthand
a + = ( 1 , + ) , b − = ( τ F i b , − ) , c + = ( τ F i b , + ) , d − = ( 1 , − ) , a+ = (1,+), \qquad b- = (\tau_{\mathrm{Fib}},-), \qquad c+ = (\tau_{\mathrm{Fib}},+), \qquad d- = (1,-), a + = ( 1 , + ) , b − = ( τ Fib , − ) , c + = ( τ Fib , + ) , d − = ( 1 , − ) ,
which are the four combinations that occur (the wrapped Y ( 2 , 1 ) = Y τ F i b η Y_{(2,1)} = Y_{\tau\mathrm{Fib}}\,\eta Y ( 2 , 1 ) = Y τ Fib η takes the values φ \varphi φ , φ − 1 \varphi^{-1} φ − 1 , − φ − 1 -\varphi^{-1} − φ − 1 , − φ -\varphi − φ on them, respectively). The identity and X X X live in a + a+ a + ; the lowest state of the whole ring above the vacuum lies in b − b- b − . Table 5 lists the leading levels of each sector at the tower-improved point for L = 8 , 10 , 12 L = 8, 10, 12 L = 8 , 10 , 12 , the DMRG values at L = 24 , 32 L = 24, 32 L = 24 , 32 , and 1 / L 1/L 1/ L extrapolations, and Figure 18 shows all levels obtained from the symmetry-reduced diagonalisation together with integer-shifted guide lines, as a check that the finite-size spectrum organises into conformal towers.
Table 5. C 5 C_5 C 5 sector-resolved spectrum at the tower-improved point ( K , λ ) = ( 0.84 , 0.34 ) (K,\lambda)=(0.84,0.34) ( K , λ ) = ( 0.84 , 0.34 ) , v = 5 / 3 v=5/3 v = 5/3 . Sectors as defined in § 5.3 . DMRG column: L = 24 L=24 L = 24 for all sectors (X X X also at L = 32 L=32 L = 32 , second value); the last column is the 1 / L 1/L 1/ L extrapolation.
Figure 18. C 5 C_{5} C 5 sector-resolved spectrum versus 1 / L 1/L 1/ L at the tower-improved point ( K , λ ) = ( 0.84 , 0.34 ) (K,\lambda) = (0.84, 0.34) ( K , λ ) = ( 0.84 , 0.34 ) , one panel and one colour per charge sector. Circles: symmetry-reduced ED at L = 8 , 10 , 12 L = 8, 10, 12 L = 8 , 10 , 12 (up to 14–16 levels per sector) and the leading two levels at L = 14 L = 14 L = 14 ; stars: DMRG at L = 24 L = 24 L = 24 (ground state, X X X and X ′ X' X ′ in a + a+ a + ; the two lowest levels in each of b − b- b − , c + c+ c + , d − d- d − ) and L = 32 L = 32 L = 32 (ground state and X X X ). Lines connect the same level across all system sizes, ED and DMRG alike. Red dashed lines mark the extrapolated leading primaries of each sector (Table 5 ), grey dotted lines their integer-shifted descendants Δ + n \Delta+ n Δ + n . The levels are essentially flat in 1 / L 1/L 1/ L , and the pattern is consistent with conformal towers: in the neutral sector a + a+ a + the identity module has no state at Δ = 1 \Delta= 1 Δ = 1 and a pair of states at Δ ≃ 2 \Delta\simeq2 Δ ≃ 2 (T T T and T ˉ \bar{T} T ˉ ), while X X X and X ′ X' X ′ each come with a descendant one unit higher; in the charged sectors the leading primaries at 0.055/0.073 (b − b- b − ), 0.16/0.40 (c + c+ c + ) and 0.45/0.61 (d − d- d − ) are followed by clusters of states near Δ + 1 \Delta+ 1 Δ + 1 (e.g. the pair at Δ ≃ 1.15 \Delta\simeq1.15 Δ ≃ 1.15 above the c + c+ c + primary at 0.16, as expected for L − 1 L_{-1} L − 1 and L ˉ − 1 \bar{L}_{-1} L ˉ − 1 descendants of a primary with h , h ˉ ≠ 0 h,\bar{h} \ne0 h , h ˉ = 0 ). The finite-size drifts of the subleading c + c+ c + and d − d- d − primaries (0.35 → \to → 0.40, 0.39 → \to → 0.45) are the largest in the spectrum and are the main source of uncertainty in the identification of the higher levels.
L → ∞ L \to\infty L → ∞ : g = 1.503 ( 2 ) g = 1.503(2) g = 1.503 ( 2 ) , T = 0.665 ( 10 ) \mathcal{T} = 0.665(10) T = 0.665 ( 10 ) .
Lattice β \beta β -function
Construction. We apply the step-scaling construction of § 3.6 on the single-strand ring. At p = 5 p = 5 p = 5 the k = 1 k = 1 k = 1 state of (34) is the lowest state of the sector d − d- d − (Y ( 2 , 1 ) = − φ Y_{(2,1)} = -\varphi Y ( 2 , 1 ) = − φ , η = − 1 \eta= -1 η = − 1 ). The Casimir coupling u C u_C u C of (35) is formed from its energy at consecutive even sizes, and the β \beta β -function B C B_C B C of (36) from u C u_C u C on successive size pairs. We follow the ray through the tower-improved point, λ = 0.405 K \lambda= 0.405K λ = 0.405 K , so that K = 0.84 K = 0.84 K = 0.84 is the point (40) . The energies come from the symmetry-reduced ED of § 3.5 (η = − 1 \eta=-1 η = − 1 block) for L = 6 , … , 16 L=6,\ldots,16 L = 6 , … , 16 and from DMRG [49 , 50 ] for L = 18 , 20 , 22 L=18,20,22 L = 18 , 20 , 22 , where the d − d- d − state is made the ground state by adding to H H H the wrapped lines with the signs of its charges, H ′ = H + 0.6 Y ( 2 , 1 ) + 0.5 η H'=H+0.6Y_{(2,1)}+0.5\eta H ′ = H + 0.6 Y ( 2 , 1 ) + 0.5 η (bond dimension χ = 120 \chi=120 χ = 120 , truncation error ≲ 4 × 10 − 8 \lesssim4\times10^{-8} ≲ 4 × 1 0 − 8 ; at L = 12 L=12 L = 12 the DMRG energy agrees with ED to 4 × 10 − 7 4\times10^{-7} 4 × 1 0 − 7 ). The size pairs are ℓ = ( 7 , 9 ) , … , ( 13 , 15 ) \ell=(7,9),\ldots,(13,15) ℓ = ( 7 , 9 ) , … , ( 13 , 15 ) from ED and ( 15 , 17 ) (15,17) ( 15 , 17 ) , ( 17 , 19 ) (17,19) ( 17 , 19 ) , ( 19 , 21 ) (19,21) ( 19 , 21 ) with DMRG.
Result. Figure 19 shows B C B_{\mathrm{C}} B C against A C A_{\mathrm{C}} A C . The curve has the structure of a flow to an IR-stable fixed point: it leaves D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) with positive slope, turns over and crosses zero with negative slope; the strong-coupling regime beyond K = 1.6 K=1.6 K = 1.6 , where the curve turns back up and whose nature is not yet understood, is left for future work. Along this ray, however, the successive size pairs do not collapse: the whole curve shifts upward from pair to pair and the zero drifts,
Figure 19. Lattice β \beta β -function at p = 5 p=5 p = 5 along the ray through the tower-improved point, H = H 0 + K Φ c l + λ ( n L + n R ) H=H_0+K\Phi_{\mathrm{cl}}+\lambda(n_L+n_R) H = H 0 + K Φ cl + λ ( n L + n R ) with λ = 0.405 K \lambda=0.405K λ = 0.405 K (so that K = 0.84 K=0.84 K = 0.84 is the point (40) ), from the Casimir coupling u C u_C u C of the k = 1 k=1 k = 1 state of the single-strand ring: B C = d u C / d ln ℓ B_C=du_C/d\ln\ell B C = d u C / d ln ℓ of (36) versus the running coordinate A C A_C A C , for 0 < K ≤ 1.6 0<K\leq1.6 0 < K ≤ 1.6 (the region containing the first fixed point). Circles: ED at L = 6 , … , 16 L=6,\ldots,16 L = 6 , … , 16 (pairs ℓ = ( 7 , 9 ) , … , ( 13 , 15 ) \ell=(7,9),\ldots,(13,15) ℓ = ( 7 , 9 ) , … , ( 13 , 15 ) ); squares: DMRG at L = 18 , 20 , 22 L=18,20,22 L = 18 , 20 , 22 (pairs ( 15 , 17 ) (15,17) ( 15 , 17 ) , ( 17 , 19 ) (17,19) ( 17 , 19 ) , ( 19 , 21 ) (19,21) ( 19 , 21 ) ), at K = 0.5 , 0.75 , 0.95 , 1.1 ) K=0.5,0.75,0.95,1.1) K = 0.5 , 0.75 , 0.95 , 1.1 ) . The black square is the UV defect D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) (u C = 0 u_C=0 u C = 0 ) and the dotted line the predicted UV slope 1 − Δ ϕ = 0.267 1-\Delta_\phi=0.267 1 − Δ ϕ = 0.267 . The zero crossing with negative slope is the IR-stable interface C 5 C_5 C 5 ; it drifts to larger K K K with increasing size (see text).
K ∗ = 0.345 , 0.441 , 0.526 , 0.600 ( E D , ℓ = ( 7 , 9 ) , … , ( 13 , 15 ) ) , K^*=0.345,\ 0.441,\ 0.526,\ 0.600 \qquad(\mathrm{ED},\ \ell=(7,9),\ldots,(13,15)), K ∗ = 0.345 , 0.441 , 0.526 , 0.600 ( ED , ℓ = ( 7 , 9 ) , … , ( 13 , 15 )) ,
K ∗ = 0.663 , 0.734 , 0.793 K^{*} = 0.663, 0.734, 0.793 K ∗ = 0.663 , 0.734 , 0.793 (DMRG, ℓ = ( 15 , 17 ) , ( 17 , 19 ) , ( 19 , 21 ) \ell= (15, 17), (17, 19), (19, 21) ℓ = ( 15 , 17 ) , ( 17 , 19 ) , ( 19 , 21 ) ), with A C ∗ A_{C}^{*} A C ∗ rising from 0.12 to 0.24 and the slope at the zero between − 0.14 -0.14 − 0.14 and − 0.18 -0.18 − 0.18 . Notice the leading irrelevant operator in the neutral channel has Δ X − 1 = − 0.12 \Delta_{X} - 1 = -0.12 Δ X − 1 = − 0.12 , causing strong finite-size drift of the B C = 0 B_{C} = 0 B C = 0 point.
C5 as an interface of the three-state Potts model
The tetracritical Ising model M ( 5 , 6 ) M(5, 6) M ( 5 , 6 ) and the critical three-state Potts model are related by gauging a Z 2 \mathbb{Z}_{2} Z 2 [68 , 69 , 70 ] : the Potts model is the orbifold of M ( 5 , 6 ) M(5, 6) M ( 5 , 6 ) by the invertible line η = D ( 4 , 1 ) \eta=D_{(4,1)} η = D ( 4 , 1 ) , and conversely M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) is the orbifold of Potts by charge conjugation C C C . Under this correspondence the S 3 = Z 3 ⋊ Z 2 S_3=\mathbb{Z}_3\rtimes\mathbb{Z}_2 S 3 = Z 3 ⋊ Z 2 symmetry of Potts becomes the fusion category Rep ( S 3 ) \operatorname{Rep}(S_3) Rep ( S 3 ) on the M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) side: the sign representation is η \eta η , and the two-dimensional representation, i.e. the sum ω + ω 2 \omega+\omega^2 ω + ω 2 of the two Z 3 \mathbb{Z}_3 Z 3 lines, is the non-invertible line D ( 1 , 3 ) D_{(1,3)} D ( 1 , 3 ) of quantum dimension 2. The line D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) commutes with η \eta η and survives the gauging as the Fibonacci line τ F i b \tau_{\mathrm{Fib}} τ Fib of the Potts model (quantum dimension φ \varphi φ , associated with the h = 2 / 5 h=2/5 h = 2/5 field) [5 ] . Any η \eta η -symmetric interface of M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) therefore defines an interface of the Potts model, and since C 5 C_5 C 5 is reached from D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) by an η \eta η -even deformation, the flow D ( 2 , 1 ) → C 5 D_{(2,1)}\to C_5 D ( 2 , 1 ) → C 5 is at the same time a flow τ F i b → C 5 P o t t s \tau_{\mathrm{Fib}}\to C^{\mathrm{Potts}}_5 τ Fib → C 5 Potts of Potts interfaces. In this subsection we work out the operator content of C 5 P o t t s C^{\mathrm{Potts}}_5 C 5 Potts and identify the symmetry under which it is a stable interface.
Orbifold rule and lattice implementation
Let I I I be the gauging interface [3 , 70 ] , with M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) on one side and Potts on the other; it has quantum dimension 2 \sqrt{2} 2 and obeys I I † = 1 ⊕ η II^\dagger=1\oplus\eta I I † = 1 ⊕ η on the M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) side and I † I = 1 ⊕ C I^\dagger I=1\oplus C I † I = 1 ⊕ C on the Potts side. The Potts interface associated with C 5 C_5 C 5 is X = I † C 5 I X=I^\dagger C_5I X = I † C 5 I , and its interface operators live on a Potts ring with X X X and X ˉ \bar{X} X ˉ inserted antipodally. Fusing, X X ˉ = I † C 5 ( I I † ) C ˉ 5 I X\bar{X}=I^\dagger C_5(II^\dagger)\bar{C}_5I X X ˉ = I † C 5 ( I I † ) C ˉ 5 I , so on each of the two arcs between the M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) strands one sums over 1 and η \eta η ; of the four resulting η \eta η -line configurations only two are inequivalent, because C 5 C_5 C 5 commutes with η \eta η and an η \eta η line on one arc can be slid through C 5 C_5 C 5 onto the other. The remaining η \eta η loop of I I † I I^\dagger I I † , wrapping the ring in the time direction, acts as the projector 1 2 ( 1 + η ) \frac{1}{2}(1+\eta) 2 1 ( 1 + η ) . Hence
H X X ˉ = [ H C 5 C ˉ 5 ⊕ H C 5 η C ˉ 5 ] η -even (42) \mathcal{H}_{X\bar{X}}=\left[\mathcal{H}_{C_5\bar{C}_5}\oplus\mathcal{H}_{C_5\eta\bar{C}_5}\right]^{\eta\text{-even}} \tag*{(42)} H X X ˉ = [ H C 5 C ˉ 5 ⊕ H C 5 η C ˉ 5 ] η -even (42)
the usual orbifold prescription: the η \eta η -even states of the untwisted interface Hilbert space, plus the η \eta η -even states of the η \eta η -twisted space in which an η \eta η line is inserted between the two interfaces. The gauged theory carries the dual (“quantum”) Z 2 \mathbb{Z}_2 Z 2 of the gauging [69 , 70 ] : gauging η \eta η means summing over η \eta η lines wrapping both cycles of the torus, the line along space being the projector 1 2 ( 1 + η ) \frac{1}{2}(1+\eta) 2 1 ( 1 + η ) and the line along time defining the twisted sector, and the dual symmetry is the operator that measures the latter, the holonomy of the Z 2 \mathbb{Z}_2 Z 2 gauge field around the spatial circle. It therefore acts as + 1 +1 + 1 on the untwisted summand of (42) and as − 1 -1 − 1 on the η \eta η -twisted one; it is a symmetry because an η \eta η line in the time direction is topological and cannot be removed by the dynamics, so the two summands never mix, and gauging it in turn returns M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) . Since M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) is the orbifold of Potts by charge conjugation, this dual Z 2 \mathbb{Z}_2 Z 2 is C C C itself (the same statement in the opposite direction is that the C C C -twisted sector of Potts supplies the η \eta η -odd fields of M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) , such as ( 2 , 2 ) (2,2) ( 2 , 2 ) with h = 1 / 40 h=1/40 h = 1/40 ; the Ising model, whose Z 2 \mathbb{Z}_2 Z 2 orbifold is again Ising with the disorder field μ \mu μ of the twisted sector as the new odd spin, is the familiar prototype). The only convention is the overall sign of C C C on the twisted sector, fixed by C = + 1 C=+1 C = + 1 on the untwisted vacuum. Note that η × C 5 \eta\times C_5 η × C 5 and C 5 C_5 C 5 give the same X X X , since η \eta η is absorbed by I I I ; the two are distinct interfaces of M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) (the twisted space below contains no identity operator) but indistinguishable from the Potts side. On the lattice both spaces in (42) are two-strand rings of § 3.5 : the untwisted one is the ring used throughout § 5 , the twisted one is the same ring with a third, η = ( 4 , 1 ) \eta=(4,1) η = ( 4 , 1 ) strand placed midway between the two C 5 C_5 C 5 strands (§ 3.2 ; the η \eta η strand carries no deformation). Two charges label the resulting states. The Fibonacci charge q q q is the eigenvalue of the wrapped line Y ( 2 , 1 ) Y_{(2,1)} Y ( 2 , 1 ) (equal to Y ( 3 , 1 ) Y_{(3,1)} Y ( 3 , 1 ) on η \eta η -even states), q = 1 q=1 q = 1 for eigenvalue φ \varphi φ and q = τ F i b q=\tau_{\mathrm{Fib}} q = τ Fib for − φ − 1 -\varphi^{-1} − φ − 1 . The charge under the dual Z 2 \mathbb{Z}_2 Z 2 of the gauging—charge conjugation C C C of the Potts model—is, by the argument above, simply the ring a state comes from: + 1 +1 + 1 on the untwisted and − 1 -1 − 1 on the twisted ring. The Z 3 \mathbb{Z}_3 Z 3 charge is not diagonal in this basis, because the Z 3 \mathbb{Z}_3 Z 3 lines map the untwisted sector into the twisted one; on a Z 3 \mathbb{Z}_3 Z 3 -symmetric interface it shows up instead as exact degeneracies between C = + C=+ C = + and C = − C=- C = − states (the ω \omega ω and ω 2 \omega^2 ω 2 components of a Z 3 \mathbb{Z}_3 Z 3 doublet are exchanged by C C C , so both C C C -eigencombinations must be present at the same energy).
Check at the topological point. At K = λ = 0 K=\lambda=0 K = λ = 0 the construction must reproduce the interface operators of the Potts τ F i b \tau_{\mathrm{Fib}} τ Fib line, and it does. The merged η \eta η -even spectrum at L = 10 L=10 L = 10 is 0 0 0 ; 2 / 15 2/15 2/15 with multiplicity four (q = τ F i b q=\tau_{\mathrm{Fib}} q = τ Fib , twice C = + C=+ C = + and twice C = − C=- C = − ); 2 / 5 2/5 2/5 twice (q = 1 q=1 q = 1 and q = τ F i b q=\tau_{\mathrm{Fib}} q = τ Fib , both C = − C=- C = − ); 11 / 15 11/15 11/15 with multiplicity four (q = 1 , τ F i b q=1,\tau_{\mathrm{Fib}} q = 1 , τ Fib , C = ± C=\pm C = ± ); 0.80 0.80 0.80 twice (q = τ F i b q=\tau_{\mathrm{Fib}} q = τ Fib , C = + C=+ C = + ). The fourfold levels are the Z 3 \mathbb{Z}_3 Z 3 doublets built on the Potts spin field σ \sigma σ (h = 1 / 15 h=1/15 h = 1/15 ) and on the h = 2 / 3 h=2/3 h = 2/3 field, split evenly between C = ± C=\pm C = ± exactly as the orbifold rule demands, while the h = 2 / 5 h=2/5 h = 2/5 -type levels, which are Z 3 \mathbb{Z}_3 Z 3 neutral, appear without C C C partners.
Spectrum of C 5 P o t t s C_5^{\mathrm{Potts}} C 5 Potts . Table 6 and Figure 20 give the interface-operator spectrum at the improved point ( K , λ ) = ( 0.84 , 0.34 ) (K,\lambda)=(0.84,0.34) ( K , λ ) = ( 0.84 , 0.34 ) for L = 8 , 10 , 12 , 14 , 16 L=8,10,12,14,16 L = 8 , 10 , 12 , 14 , 16 (ED) and L = 24 L=24 L = 24 (DMRG, two leading states per sector), sector by sector. The L = 24 L=24 L = 24 points are obtained as in § 5.3 , with the η \eta η -twisted ring realised in the MPS by reflecting the heights across one site (the wrapped ( 2 , 1 ) (2,1) ( 2 , 1 ) line is closed with the same reflection); the construction reproduces the ED twisted spectra at L = 8 , 12 L=8,12 L = 8 , 12 to all digits. The lowest levels are converged to two or three digits, and the linear extrapolations of the leading levels through L = 14 , 16 , 24 L=14,16,24 L = 14 , 16 , 24 are 0.060, 0.078, 0.754 in [ τ F i b , − ] [\tau_{\mathrm{Fib}},-] [ τ Fib , − ] , 0.164, 0.404, 0.622 in [ τ F i b , + ] [\tau_{\mathrm{Fib}},+] [ τ Fib , + ] , 0.484, 0.592 in [ 1 , − ] [1,-] [ 1 , − ] and X = 1.116 X=1.116 X = 1.116 in [ 1 , + ] [1,+] [ 1 , + ] (the third τ F i b \tau_{\mathrm{Fib}} τ Fib levels from ED alone, L = 12 , 14 , 16 L=12,14,16 L = 12 , 14 , 16 ). The Z 3 \mathbb{Z}_3 Z 3 doublets of the UV point are gone: the four states at 2 / 15 2/15 2/15 have split into two C = − C=- C = − levels at 0.06, 0.08 and two C = + C=+ C = + levels at 0.16, 0.40, and no C = + C=+ C = + level has a degenerate C = − C=- C = − partner anywhere in the spectrum. The interface C 5 P o t t s C_5^{\mathrm{Potts}} C 5 Potts therefore breaks Z 3 \mathbb{Z}_3 Z 3 . The deformation φ ( 1 , 3 ) ⊗ φ ˉ ( 3 , 3 ) \varphi_{(1,3)}\otimes\bar{\varphi}_{(3,3)} φ ( 1 , 3 ) ⊗ φ ˉ ( 3 , 3 ) is charged under the lasso of D ( 1 , 3 ) D_{(1,3)} D ( 1 , 3 ) (its eigenvalue is S ( 1 , 3 ) ( 1 , 3 ) / S ( 1 , 1 ) ( 1 , 3 ) = − 1 = ω + ω 2 S_{(1,3)(1,3)}/S_{(1,1)(1,3)}=-1=\omega+\omega^2 S ( 1 , 3 ) ( 1 , 3 ) / S ( 1 , 1 ) ( 1 , 3 ) = − 1 = ω + ω 2 rather than the neutral value 2), i.e. in Potts language it is one real component of a Z 3 \mathbb{Z}_3 Z 3 doublet of interface operators on τ F i b \tau_{\mathrm{Fib}} τ Fib , and turning it on selects one of three Z 3 \mathbb{Z}_3 Z 3 -related images of the interface, in the same way that a fixed boundary condition of the Potts model selects one of the three colours. C 5 P o t t s C_5^{\mathrm{Potts}} C 5 Potts preserves the Z 2 = ⟨ C ⟩ \mathbb{Z}_2=\langle C\rangle Z 2 = ⟨ C ⟩ subgroup of S 3 S_3 S 3 (the dual of the η \eta η symmetry of the lattice model) and the Fibonacci line τ F i b \tau_{\mathrm{Fib}} τ Fib .
Table 6. Interface operators of C 5 P o t t s C_5^{\mathrm{Potts}} C 5 Potts at ( K , λ ) = ( 0.84 , 0.34 ) (K,\lambda)=(0.84,0.34) ( K , λ ) = ( 0.84 , 0.34 ) , labelled [ q , C ] [q,C] [ q , C ] and listed by increasing dimension at L = 12 L=12 L = 12 . The row marked * is DMRG at L = 24 L=24 L = 24 (χ = 120 \chi=120 χ = 120 ): the two leading states of each sector, from the untwisted ring for C = + C=+ C = + and from the η \eta η -twisted ring for C = − C=- C = − . The only Fibonacci-neutral, C C C -even operator below the marginal line is the identity; the first non-trivial one is the leading irrelevant operator X X X (1.095 at L = 24 L=24 L = 24 , extrapolating to 1.11), the same as on the M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) side (§ 5.3 ).
Figure 20. Sector-resolved interface spectrum of C 5 P o t t s C_5^{\mathrm{Potts}} C 5 Potts at ( K , λ ) = ( 0.84 , 0.34 ) (K,\lambda)=(0.84,0.34) ( K , λ ) = ( 0.84 , 0.34 ) against 1 / L 1/L 1/ L , one panel per [ q , C ] [q,C] [ q , C ] sector, in the style of Figure 18 : circles are the ED levels at L = 8 , 10 , 12 , 14 , 16 L=8,10,12,14,16 L = 8 , 10 , 12 , 14 , 16 (untwisted ring for C = + C=+ C = + , η \eta η -twisted ring for C = − C=- C = − ), stars the two leading DMRG states per sector at L = 24 L=24 L = 24 , lines connect the same level across sizes, open diamonds are linear-in-1 / L 1/L 1/ L extrapolations of the leading levels through the three largest sizes (L = 14 , 16 , 24 L=14,16,24 L = 14 , 16 , 24 ; ED only, L = 12 , 14 , 16 L=12,14,16 L = 12 , 14 , 16 , for the third τ F i b \tau_{\mathrm{Fib}} τ Fib levels), dotted grey lines mark the values at the topological UV point (τ F i b \tau_{\mathrm{Fib}} τ Fib line), and the dashed line is the marginal dimension. Levels above Δ = 1.2 \Delta=1.2 Δ = 1.2 are omitted. The only state below the marginal line in the fully neutral sector [ 1 , + ] [1,+] [ 1 , + ] is the identity; the leading irrelevant operator extrapolates to X ≃ 1.11 X\simeq1.11 X ≃ 1.11 , as on the M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) side. In [ τ F i b , − ] [\tau_{\mathrm{Fib}},-] [ τ Fib , − ] the grey dashed lines mark Δ + 1 \Delta+1 Δ + 1 for the two leading primaries, 1.06 and 1.08: at the UV point the two Δ = 2 / 15 \Delta=2/15 Δ = 2/15 states have four level-one descendants at 17 / 15 17/15 17/15 (left- and right-moving), and along the flow these become the levels at 0.99, 1.15 (shown) and 1.31, 1.38 (above the frame). On the conformal interface only a single Virasoro algebra acts, so only two of them remain L − 1 L_{-1} L − 1 descendants; the pair at 0.99/1.15 straddles the expected positions with a splitting that hardly moves between L = 8 L=8 L = 8 and 16 (0.98 → \to → 0.99 and 1.145 → \to → 1.150), as expected from mixing through the nearly marginal X X X (Δ X ≃ 1.11 \Delta_X\simeq1.11 Δ X ≃ 1.11 , corrections ∝ L − 0.11 \propto L^{-0.11} ∝ L − 0.11 ), while their mean, 1.07, agrees with the mean of the two expected descendant positions.
Stability. An interface is stable under a symmetry if none of its relevant operators is neutral under that symmetry. Reading Table 6 : the relevant operators of C 5 P o t t s C_5^{\mathrm{Potts}} C 5 Potts are the [ τ F i b , − ] [\tau_{\mathrm{Fib}},-] [ τ Fib , − ] levels at 0.06, 0.08, 0.75, 0.79, 0.99, the [ τ F i b , + ] [\tau_{\mathrm{Fib}},+] [ τ Fib , + ] levels at 0.16, 0.40, 0.62, 0.74, 0.85, and the [ 1 , − ] [1,-] [ 1 , − ] levels at 0.48, 0.59; the first [ 1 , + ] [1,+] [ 1 , + ] operator is X ∼ 1.11 X\sim1.11 X ∼ 1.11 . Hence C 5 P o t t s C_5^{\mathrm{Potts}} C 5 Potts is stable if and only if both the Fibonacci line τ F i b \tau_{\mathrm{Fib}} τ Fib and charge conjugation C C C are imposed. Imposing C C C alone leaves the relevant Fibonacci-charged, C C C -even operators at 0.16, 0.40 and 0.62; imposing τ F i b \tau_{\mathrm{Fib}} τ Fib alone leaves the Fibonacci-neutral, C C C -odd operators at 0.48 and 0.59. Since the interface breaks Z 3 \mathbb{Z}_3 Z 3 , the full S 3 S_3 S 3 cannot be used, and Z 3 \mathbb{Z}_3 Z 3 plays no role in the stability count; in a Potts lattice realisation the interface would have to be protected by the non-invertible line τ F i b \tau_{\mathrm{Fib}} τ Fib together with C C C , exactly as C 5 C_5 C 5 is by Y ( 2 , 1 ) Y_{(2,1)} Y ( 2 , 1 ) together with η \eta η on the M ( 5 , 6 ) M(5,6) M ( 5 , 6 ) side. The g g g -function and the transmission are unchanged by the orbifold, since ⟨ 0 ∣ D ∣ 0 ⟩ \langle0|D|0\rangle ⟨ 0∣ D ∣0 ⟩ and the current two-point function only involve the untwisted vacuum: g C 5 P o t t s = 1.503 ( 2 ) g_{C_5}^{\mathrm{Potts}} = 1.503(2) g C 5 Potts = 1.503 ( 2 ) and T = 0.665 ( 10 ) T = 0.665(10) T = 0.665 ( 10 ) .
The loop-model formulation at continuous p p p
The RSOS construction fixes p p p to an integer. Everything used above, however, lives in the dilute Temperley–Lieb algebra [71 ] with loop weight n = 2 cos ( π / p ) = d w n = 2\cos(\pi/p) = d_w n = 2 cos ( π / p ) = d w , and can be written in the loop (dilute O ( n ) O(n) O ( n ) ) model [35 , 72 , 73 , 74 ] , where p p p is a continuous parameter. This section does so and follows the fixed point C p C_p C p from p = 5 p = 5 p = 5 down to p = 4 p = 4 p = 4 .
Loop model with seams
Hilbert space. We use the oriented-strand (three-state vertex) representation of the dilute O ( n ) O(n) O ( n ) model [35 , 73 , 74 ] , Figure 21 (a). Each slot i i i of the ring carries s i ∈ { 0 , + , − } s_i \in\{0, +, -\} s i ∈ { 0 , + , − } : either it is empty or it is crossed by a strand of the loop gas, oriented upwards (+ + + ) or downwards (− - − ) in the time direction. A ring of L L L sites with one seam has M = L + 1 M = L + 1 M = L + 1 slots, the seam being an extra slot c c c inserted between a L − 1 a_{L-1} a L − 1 and a 0 a_0 a 0 (Figure 22 (a)) which is never empty : it is the world-line of the ( 2 , 1 ) (2,1) ( 2 , 1 ) strand. The total orientation S z = ∑ i s i S^z = \sum_i s_i S z = ∑ i s i , seam slot included, is conserved, and we work at S z = 0 S^z = 0 S z = 0 , the sector that contains every state with all strands closed into loops; at L = 12 L = 12 L = 12 the seam ring has 3 13 3^{13} 3 13 configurations of which dim = 139 152 \dim= 139\,152 dim = 139 152 survive the two constraints. Because the seam slot is always occupied there is no separate seam label to track: what remains dynamical is the orientation s c = ± s_c = \pm s c = ± of the seam strand and how it connects to the bulk strands, and that orientation changes under the dynamics, since the pair term can project the seam strand and its neighbour onto a singlet, annihilate the pair and recreate it. The seam strand carries the twist phase when it crosses the cut like any other strand.
Figure 21. The loop model in the oriented-strand representation. (a) Slot states. (b) The six face moves of (45) , drawn with time upwards; the coefficients are the RSOS weights (13) . (c) Orientation phases: a cup or cap carries q ± 1 / 2 q^{\pm1/2} q ± 1/2 according to its orientation, so the two-slot singlet is ∣ c ⟩ = q 1 / 2 ∣ + − ⟩ + q − 1 / 2 ∣ − + ⟩ \lvert c\rangle= q^{1/2}\lvert+-\rangle+ q^{-1/2}\lvert-+\rangle ∣ c ⟩ = q 1/2 ∣ + − ⟩ + q − 1/2 ∣ − + ⟩ (and ⟨ c ~ ∣ \langle\tilde{c}\rvert ⟨ c ~ ∣ likewise), and a contractible loop, summed over its two orientations, has weight ⟨ c ~ ∣ c ⟩ = q + q − 1 = n \langle\tilde{c}\vert c\rangle= q + q^{-1} = n ⟨ c ~ ∣ c ⟩ = q + q − 1 = n , which is the Temperley–Lieb relation e 2 = n e e^2 = ne e 2 = n e .
Figure 22. Seam and twist. (a) The single-seam ring: L L L site slots and one always-occupied seam slot c c c , the L − 1 L-1 L − 1 ordinary pair terms, the three-slot seam term carrying the deformation, and the cut at which the twist is applied. (b) The seam term (46) : the strand is braided in front of the content of a 0 a_0 a 0 , the ordinary h h h acts, and the strand is braided back with the inverse braid, again in front of a 0 a_0 a 0 , so that the two crossings cancel. (c) The braid on a pair of slots, (47) . (d) A twist ϕ \phi ϕ gives non-contractible loops the weight 2 cos ϕ 2\cos\phi 2 cos ϕ while contractible loops keep n n n ; the twist sector is the loop-model version of the Y ( 2 , 1 ) Y_{(2,1)} Y ( 2 , 1 ) eigenvalue, (49) .
Loop weight from orientations. Loops acquire their weight from the phases attached to their turns (Figure 21 (c)). A strand that turns by ± π \pm\pi ± π (a cup or a cap) carries q ± 1 / 2 q^{\pm1/2} q ± 1/2 according to its orientation, where q = e i π / p q = e^{i\pi/p} q = e iπ / p ; a closed contractible loop turns by ± 2 π \pm2\pi ± 2 π in total and weighs q ± 1 q^{\pm1} q ± 1 , so that the sum over its two orientations gives q + q − 1 = 2 cos ( π / p ) = n q + q^{-1} = 2\cos(\pi/p) = n q + q − 1 = 2 cos ( π / p ) = n . Concretely, with ∣ c ⟩ = q 1 / 2 ∣ + − ⟩ + q − 1 / 2 ∣ − + ⟩ \lvert c\rangle= q^{1/2}\lvert+-\rangle+ q^{-1/2}\lvert-+\rangle ∣ c ⟩ = q 1/2 ∣ + − ⟩ + q − 1/2 ∣ − + ⟩ and ⟨ c ~ ∣ = q 1 / 2 ⟨ + − ∣ + q − 1 / 2 ⟨ − + ∣ \langle\widetilde{c}\rvert= q^{1/2}\langle+-\rvert+ q^{-1/2}\langle-+\rvert ⟨ c ∣ = q 1/2 ⟨ + − ∣ + q − 1/2 ⟨ − + ∣ , the Temperley–Lieb generator and the cup and cap on two adjacent slots are
e = ∣ c ⟩ ⟨ c ~ ∣ , c u p = ∣ c ⟩ ⟨ 00 ∣ , c a p = ∣ 00 ⟩ ⟨ c ~ ∣ , e 2 = ⟨ c ~ ∣ c ⟩ e = n e , (43) e = \lvert c\rangle\langle\widetilde{c}\rvert,\qquad\mathrm{cup} = \lvert c\rangle\langle00\rvert,\qquad\mathrm{cap} = \lvert00\rangle\langle\widetilde{c}\rvert,\qquad e^2 = \langle\widetilde{c}\vert c\rangle e = n e, \tag*{(43)} e = ∣ c ⟩ ⟨ c ∣ , cup = ∣ c ⟩ ⟨ 00 ∣ , cap = ∣ 00 ⟩ ⟨ c ∣ , e 2 = ⟨ c ∣ c ⟩ e = n e , (43)
so that e e e is “cap followed by cup”. Note that ⟨ c ~ ∣ \langle\widetilde{c}\rvert ⟨ c ∣ is not the Hermitian conjugate of ∣ c ⟩ \lvert c\rangle ∣ c ⟩ (that would carry q ∓ 1 / 2 q^{\mp1/2} q ∓ 1/2 ): for complex q q q the representation is non-Hermitian, which is the price for a continuous p p p (complex-symmetric defect Hamiltonians of the same kind arise for lattice defects of the non-unitary minimal models [75 ] ), and we diagonalise with left and right eigenvectors throughout. The eigenvalues are nevertheless real on all the Kac levels , by which we mean the states that belong to the minimal model proper: those whose scaling limit lies in a representation V ( r , s ) ⊗ V ( r ′ , s ′ ) V_{(r,s)} \otimes V_{(r',s')} V ( r , s ) ⊗ V ( r ′ , s ′ ) with both labels inside the Kac table of M ( p , p + 1 ) M(p,p+1) M ( p , p + 1 ) , 1 ≤ r ≤ p − 1 1 \le r \le p - 1 1 ≤ r ≤ p − 1 , 1 ≤ s ≤ p 1 \le s \le p 1 ≤ s ≤ p , so that their dimension is h r , s + h r ′ , s h_{r,s} + h_{r',s} h r , s + h r ′ , s plus integers. These are the states that have a counterpart in the RSOS chain, with the same energies at integer p p p , varying smoothly in between. The loop model also contains non-Kac levels, with labels outside the table or outside the U q ( s l 2 ) U_q(sl_2) U q ( s l 2 ) highest-weight sector, which the RSOS restriction removes; they are discussed in § 6.2 .
Heights and orientations. The RSOS chain of § 3 is related to this vertex chain by a map of configurations. A domain wall of the height model is a point where a loop crosses the chain, and the sign of the height step across it, a i − a i − 1 = ± 1 a_i-a_{i-1}=\pm1 a i − a i − 1 = ± 1 , says on which side of the loop the height is larger. Since the inside of a closed domain wall differs from its outside by the same ± 1 \pm1 ± 1 all along the wall, this sign is constant around each loop: it is an orientation of the loop. The map s i = a i − a i − 1 ∈ { 0 , + , − } s_i=a_i-a_{i-1}\in\{0,+,-\} s i = a i − a i − 1 ∈ { 0 , + , − } therefore sends a height configuration to an oriented loop configuration, injectively (given a 0 a_0 a 0 , the heights are the running sum of the steps), its image being restricted by 1 ≤ a i ≤ p − 1 1\leq a_i\leq p-1 1 ≤ a i ≤ p − 1 and, on the ring, by ∑ i s i = 0 \sum_i s_i=0 ∑ i s i = 0 ; this is why the loop model is studied at S z = 0 S^z=0 S z = 0 , and for the seam slot it is the statement that the RSOS seam label c j = a j ± 1 c_j=a_j\pm1 c j = a j ± 1 is the orientation s c = ± s_c=\pm s c = ± of the seam strand. What the two descriptions do not share is the weight attached to an orientation. In the height representation the generators (43) act, on a triple of heights with equal ends l l l , as
e ∣ l a l ⟩ = S a S a ′ S l ∣ l a ′ l ⟩ , cup ∣ l l l ⟩ = S a S l ∣ l a l ⟩ , cap ∣ l a l ⟩ = S a S l ∣ l l l ⟩ , (44) e\lvert l\,a\,l\rangle=\sqrt{\frac{S_aS_{a'}}{S_l}}\lvert l\,a'\,l\rangle,\qquad\operatorname{cup}\lvert l\,l\,l\rangle=\sqrt{\frac{S_a}{S_l}}\lvert l\,a\,l\rangle,\qquad\operatorname{cap}\lvert l\,a\,l\rangle=\sqrt{\frac{S_a}{S_l}}\lvert l\,l\,l\rangle, \tag*{(44)} e ∣ l a l ⟩ = S l S a S a ′ ∣ l a ′ l ⟩ , cup ∣ l l l ⟩ = S l S a ∣ l a l ⟩ , cap ∣ l a l ⟩ = S l S a ∣ l l l ⟩ , (44)
with a , a ′ a,a' a , a ′ adjacent to l l l (these are the TL and cup/cap moves w e w_e w e , w ∪ ∩ w_{\cup\cap} w ∪∩ of (13) and the projector P 1 = e / d w P^1=e/d_w P 1 = e / d w of Figure 10 ). A contractible loop with outside height a a a and inside height a ± 1 a\pm1 a ± 1 thus weighs S a ± 1 / S a S_{a\pm1}/S_a S a ± 1 / S a , where the loop gas gives q ± 1 q^{\pm1} q ± 1 for the corresponding orientation, and the two agree only after the sum over the two orientations, ∑ a ′ a d j a S a ′ / S a = n = q + q − 1 \sum_{a'\,\mathrm{adj}\,a}S_{a'}/S_a=n=q+q^{-1} ∑ a ′ adj a S a ′ / S a = n = q + q − 1 (the height-basis form of ⟨ c ~ ∣ c ⟩ = n \langle\widetilde c|c\rangle=n ⟨ c ∣ c ⟩ = n ), not term by term; for a non-contractible loop they differ even then, which is where the twist introduced below enters. The height and vertex chains are therefore two representations of the same elements of the dilute Temperley–Lieb algebra rather than two gauges of one matrix. What coincides is the spectrum, sector by sector, through the standard RSOS–vertex correspondence [33 , 76 ] : the RSOS sectors are reproduced by twists of the S z = 0 S^z=0 S z = 0 sector (point (iii) below), and the vertex representation simply continues to exist when p p p is not an integer.
Hamiltonian. The pair term is the element of the dilute Temperley–Lieb algebra that defines the RSOS face operator of § 3.1 , Figure 21 (b), now in the vertex representation:
h = w ∅ 1 ∅ + w ∥ 1 ∥ + w = 1 = + w ⇄ ( hop ) + w ∪ ∩ ( cup + cap ) + w e e , H 0 = − ∑ j h j , j + 1 , (45) h=w_{\varnothing}1_{\varnothing}+w_{\parallel}1_{\parallel}+w_{=}1_{=}+w_{\rightleftarrows}(\operatorname{hop})+w_{\cup\cap}(\operatorname{cup}+\operatorname{cap})+w_e e,\qquad H_0=-\sum_j h_{j,j+1}, \tag*{(45)} h = w ∅ 1 ∅ + w ∥ 1 ∥ + w = 1 = + w ⇄ ( hop ) + w ∪∩ ( cup + cap ) + w e e , H 0 = − j ∑ h j , j + 1 , (45)
where 1 ∅ 1_{\varnothing} 1 ∅ , 1 ∥ 1_{\parallel} 1 ∥ , 1 = 1_{=} 1 = project on the pair configurations with zero, one and two strands, “hop” moves a single strand to the other slot of the pair, and cup, cap and e e e are the vertex generators (43) . The six scalar coefficients are the same functions of λ = π ( p − 1 ) / 4 p \lambda=\pi(p-1)/4p λ = π ( p − 1 ) /4 p as in (13) , now with p p p real: the first four are unchanged, and w ∪ ∩ = sin 2 λ / N 0 w_{\cup\cap}=\sin2\lambda/\sqrt{\mathcal N_0} w ∪∩ = sin 2 λ / N 0 , w e = − sin λ / N 0 w_e=-\sin\lambda/\sqrt{\mathcal N_0} w e = − sin λ / N 0 , the height-dependent square roots of (13) being the matrix elements (44) of the generators rather than part of the coefficients. In the vertex language h h h is a 9 × 9 9\times9 9 × 9 matrix acting on a pair of slots.
The seam. The pair term that straddles the seam is built exactly as in § 3.2 : the strand is braided out of the way, the ordinary h h h acts on the two physical contents, and the strand is braided back (Figure 22 (b)),
h x c y s e a m = B c y − 1 h x c B c y , x = a L − 1 , y = a 0 , (46) h^{\mathrm{seam}}_{xcy}=B^{-1}_{cy}h_{xc}B_{cy}, \qquad x=a_{L-1},\ y=a_0, \tag*{(46)} h x cy seam = B cy − 1 h x c B cy , x = a L − 1 , y = a 0 , (46)
a 27 × 27 27 \times27 27 × 27 matrix on the three slots ( x , c , y ) (x,c,y) ( x , c , y ) . The braid B B B on a pair of slots is the Kauffman skein evaluation of a crossing [77 ] , Figure 22 (c): it is 1 1 1 on two empty slots, the permutation P P P when one slot is empty (the strand simply moves past a vacancy), and
B = A 1 + A − 1 e , A = i q 1 / 2 , n = − A 2 − A − 2 = q + q − 1 , (47) B=A\mathbf{1}+A^{-1}e,\qquad A=iq^{1/2},\qquad n=-A^2-A^{-2}=q+q^{-1}, \tag*{(47)} B = A 1 + A − 1 e , A = i q 1/2 , n = − A 2 − A − 2 = q + q − 1 , (47)
when both slots carry a strand. Equation (47) is the F R F − 1 FRF^{-1} F R F − 1 braiding (15) of the RSOS seam written in the skein language, and the two blocks of B B B are related by a definite relative phase: the one-strand block is P P P with coefficient 1 1 1 , the two-strand block carries A A A . This phase is not a convention. Loops that cross the seam once pick it up, and it is fixed by requiring that the two-seam ring (a D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) strand and its conjugate) reproduce the bulk vacuum energy exactly, i.e. that the seam be topological; with the wrong phase the two-seam ring is off the vacuum by a finite amount at every L L L . The braid is unitary in the RSOS representation but not in the vertex one; h s e a m h^{\mathrm{seam}} h seam is similar to h h h and has its spectrum.
One difference of bookkeeping between the two representations should be kept in mind. In the RSOS chain the seam term changes a j a_j a j and c j c_j c j only and leaves the neighbouring heights fixed (§ 3.2 ); those heights are fusion labels, not the contents of slots. In the vertex chain there is no such spectator: (46) acts on all three slots ( x , c , y ) (x,c,y) ( x , c , y ) , and the orientation, and even the occupation, of y y y can change in one step—for instance when y y y is empty, the braid moves the seam strand onto y y y , h h h then creates a pair on ( x , c ) (x,c) ( x , c ) , and the inverse braid acts on the two occupied slots ( c , y ) (c,y) ( c , y ) . What is preserved is the total orientation S z S^z S z and the occupation of the seam slot c c c after the step.
The clean field. The deformation (28) lives on the same three slots,
Φ c l = n ( e x c + e c y ) − ( n x + n y ) , n i = 1 − ∣ 0 ⟩ ⟨ 0 ∣ i , (48) \Phi_{\mathrm{cl}}=n\left(e_{xc}+e_{cy}\right)-\left(n_x+n_y\right),\qquad n_i=1-\lvert0\rangle\langle0\rvert_i, \tag*{(48)} Φ cl = n ( e x c + e cy ) − ( n x + n y ) , n i = 1 − ∣ 0 ⟩ ⟨ 0 ∣ i , (48)
where e x c e_{xc} e x c (e c y e_{cy} e cy ) is the Temperley–Lieb generator (43) on the strand and its left (right) neighbour and n i n_i n i counts whether slot i i i is occupied; this is d w ( e L + e R ) − ( n L + n R ) d_w(e_L+e_R)-(n_L+n_R) d w ( e L + e R ) − ( n L + n R ) with d w = n d_w=n d w = n . The full Hamiltonian of the single-seam ring is then H = H 0 − h s e a m + K Φ c l H=H_0-h^{\mathrm{seam}}+K\Phi_{\mathrm{cl}} H = H 0 − h seam + K Φ cl with
H 0 = − ∑ o r d i n a r y p a i r s h H_0=-\sum_{\mathrm{ordinary\ pairs}}h H 0 = − ordinary pairs ∑ h
containing the L − 1 L-1 L − 1 ordinary pair terms (Figure 22 (a)).
Twist sectors. In the RSOS calculation the symmetry sectors were obtained by diagonalising the wrapped ( 2 , 1 ) (2,1) ( 2 , 1 ) line. The loop model has no such operator at non-integer p p p , but it does not need one: the sectors can be selected directly by a boundary condition, the twist , in three steps.
(i) Conserved orientation and a boundary-condition parameter. On the ring the total orientation S z = ∑ k s k S^z=\sum_k s_k S z = ∑ k s k is conserved, since every face move creates or destroys strands in oppositely oriented pairs. In the height language S z S^z S z is the net height change around the ring, so the periodic height configurations of the RSOS ring all have S z = 0 S^z=0 S z = 0 , and we work in that sector throughout. Because S z S^z S z is conserved we may also impose a twisted boundary condition: a cut is placed at one site of the ring, between the seam slot c c c and a 0 a_0 a 0 (Figure 22 (a)), and every matrix element in which a strand of orientation s = ± 1 s = \pm1 s = ± 1 crosses the cut is multiplied by e i s ϕ e^{is\phi} e i s ϕ . The angle ϕ \phi ϕ labels a choice of Hamiltonian, rather than an additional conserved operator; in practice only the wrap-around term ((46) ) and Φ c l \Phi_{\mathrm{cl}} Φ cl act across the cut, so the twist enters only there. Concretely, the twist is not a projection but a twisted boundary condition, and each sector is a separate Hamiltonian H ( ϕ ) H(\phi) H ( ϕ ) : a matrix element that moves an up-strand from c c c to a 0 a_0 a 0 is multiplied by e i ϕ e^{i\phi} e i ϕ and one that moves it back by e − i ϕ e^{-i\phi} e − i ϕ , a down-strand gets the opposite phases, and a cup or cap created or annihilated across the cut gets the phase of the orientation that ends up on the far side; all other terms of H H H are untouched. H ( ϕ ) H(\phi) H ( ϕ ) is complex and non-Hermitian, and for each value of ϕ \phi ϕ it is diagonalised separately in the S z = 0 S^z = 0 S z = 0 sector with biorthogonal left and right eigenvectors.
(ii) What the twist weighs. Consider the loops traced out in space-time by the transfer matrix. A contractible loop crosses the cut an even number of times with opposite orientations, so the twist phases cancel and its weight stays n = 2 cos ( π / p ) n = 2\cos(\pi/p) n = 2 cos ( π / p ) , which comes entirely from its turning phases. A loop that winds once around the ring has no net turning, so its turning phases give nothing, and without a twist it would have weight 1 + 1 = 2 1 + 1 = 2 1 + 1 = 2 ; it crosses the cut exactly once, so with the twist its two orientations give e i ϕ + e − i ϕ = 2 cos ϕ e^{i\phi} + e^{-i\phi} = 2\cos\phi e i ϕ + e − i ϕ = 2 cos ϕ (Figure 22 (d)). The twist is therefore a knob that sets the weight of non-contractible loops independently of n n n ; the choice ϕ = π / p \phi= \pi/p ϕ = π / p gives them the same weight n n n as contractible loops and is the vacuum sector. This is the standard twisted-boundary-condition construction of the loop-model sectors [76 , 78 ] .
(iii) Which ϕ \phi ϕ gives which RSOS sector. The decomposition of the RSOS torus partition function into loop (Coulomb-gas) partition functions with twisted boundary conditions is standard: summing over the heights of the plateaus between loops that wind around the ring turns the sum over configurations into a sum over the eigenvalues λ r = 2 cos ( π r / p ) \lambda_r = 2\cos(\pi r/p) λ r = 2 cos ( π r / p ) , r = 1 , … , p − 1 r = 1,\ldots,p - 1 r = 1 , … , p − 1 , of the adjacency matrix of the A p − 1 A_{p-1} A p − 1 graph, i.e. into a sum of loop partition functions with non-contractible loop weight λ r \lambda_r λ r [33 , 76 , 79 ] ; in transfer-matrix language the RSOS levels are recovered from the corresponding twisted loop representations with 2 cos ϕ = λ r 2\cos\phi= \lambda_r 2 cos ϕ = λ r after the RSOS restriction. At integer p p p , this restriction selects the appropriate Kac subquotients and removes the extra states present in the unrestricted vertex representation; the full vertex spectra must not be identified with the RSOS spectrum. The label r r r is the one of the wrapped ( 2 , 1 ) (2,1) ( 2 , 1 ) line, since that line acts on the plateau heights as the adjacency matrix and is itself a non-contractible loop of the loop model. Including the sign ( − 1 ) r + s (-1)^{r+s} ( − 1 ) r + s that the η \eta η charge contributes in ((4) ), the Kac state with holomorphic label ( r , s ) (r,s) ( r , s ) has Y ( 2 , 1 ) = ( − 1 ) r + s 2 cos ( π r / p ) Y_{(2,1)} = (-1)^{r+s} 2\cos(\pi r/p) Y ( 2 , 1 ) = ( − 1 ) r + s 2 cos ( π r / p ) , hence it lives in the twist sector
ϕ = π r p ( r + s even ) , ϕ = π − π r p ( r + s odd ) , (49) \phi= \frac{\pi r}{p}\quad(r+s\ \text{even}),\qquad\phi= \pi- \frac{\pi r}{p}\quad(r+s\ \text{odd}), \tag*{(49)} ϕ = p π r ( r + s even ) , ϕ = π − p π r ( r + s odd ) , (49)
where the second case is the sign flip r → p − r r \to p-r r → p − r produced by η \eta η : the sectors ϕ \phi ϕ and π − ϕ \pi-\phi π − ϕ are distinct in the oriented model because each strand carries one unit of phase. At integer p p p , after the RSOS restriction, these correspond to the sectors a , b , c , d a,b,c,d a , b , c , d of § 5.3 , ( Y ( 2 , 1 ) , η ) = ( φ , + ) , ( φ − 1 , − ) , ( φ − 1 , + ) , ( − φ , − ) (Y_{(2,1)},\eta) = (\varphi,+),(\varphi^{-1},-),(\varphi^{-1},+),(-\varphi,-) ( Y ( 2 , 1 ) , η ) = ( φ , + ) , ( φ − 1 , − ) , ( φ − 1 , + ) , ( − φ , − ) at p = 5 p=5 p = 5 corresponding to ϕ = π / 5 , 2 π / 5 , 3 π / 5 , 4 π / 5 \phi= \pi/5, 2\pi/5, 3\pi/5, 4\pi/5 ϕ = π /5 , 2 π /5 , 3 π /5 , 4 π /5 ; the “k = 1 k = 1 k = 1 ” sector of the next subsection, with non-contractible loop weight − 2 cos ( π / p ) -2\cos(\pi/p) − 2 cos ( π / p ) , is ϕ = π − π / p \phi= \pi- \pi/p ϕ = π − π / p . Choosing ϕ \phi ϕ thus selects the twisted vertex Hamiltonian associated with a given RSOS charge sector, subject to the RSOS restriction at integer p p p ; it is not itself a projector. The twisted Hamiltonians continue to be defined at non-integer p p p , where they organize the loop-model calculation without imposing a minimal-model quotient or requiring an RSOS wrapped-line operator.
Validation. The construction was validated at p = 5 p = 5 p = 5 against the RSOS chain: the bulk ring reproduces the vacuum energy to 10 − 13 10^{-13} 1 0 − 13 and the levels of all four sectors; the two-seam ring reproduces the identity at exactly the bulk vacuum energy (the seam is topological) and every level of the D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) spectrum with its charges; the single-seam ring reproduces the single-strand spectrum, and the running coupling u C u_C u C of (35) along the clean ray, recomputed on the RSOS chain for L ≤ 14 L \le14 L ≤ 14 , agrees with the loop model to all digits. The loop model has more states than the minimal model: its sectors are Fock-type modules of the Coulomb gas [73 , 76 , 80 ] , so besides the Kac levels they contain states that are null at integer p p p (for instance, a level at Δ = 0.20 \Delta= 0.20 Δ = 0.20 below the ( 1 , 1 ) ⊗ ( 2 , 1 ) (1, 1) \otimes(2, 1) ( 1 , 1 ) ⊗ ( 2 , 1 ) state of the a a a sector at p = 5 p = 5 p = 5 , and Δ = 1 \Delta= 1 Δ = 1 current-like states in the bulk). The Hamiltonian is non-Hermitian for all p p p , so levels are tracked with left and right eigenvectors (biorthogonal overlaps).
The running coupling at continuous p p p
For continuous p p p we use the same running coupling as at integer p p p , the Casimir coupling u C u_C u C of (35) built from the k = 1 k = 1 k = 1 state of the single-seam ring, with the velocity v ( p ) = 4 p / ( 3 ( p − 1 ) ) v(p) = 4p/(3(p - 1)) v ( p ) = 4 p / ( 3 ( p − 1 )) (so v ( 4 ) = 16 / 9 v(4) = 16/9 v ( 4 ) = 16/9 and v ( 5 ) = 5 / 3 v(5) = 5/3 v ( 5 ) = 5/3 as before) and ℓ = 7 , 9 , 11 , 13 \ell= 7, 9, 11, 13 ℓ = 7 , 9 , 11 , 13 for L = 6 , … , 14 L = 6, \ldots, 14 L = 6 , … , 14 . At non-integer p p p there are no exact charges to identify the state, so it is followed continuously in K K K in steps of 0.05 0.05 0.05 by biorthogonal overlap, never re-selected as the lowest level; at K = 0 K=0 K = 0 it is the ground state of the sector with non-contractible loop weight − 2 cos ( π / p ) -2\cos(\pi/p) − 2 cos ( π / p ) , i.e. the Kac state ( 1 , 2 ) ⊗ ( 2 , 2 ) (1,2)\otimes(2,2) ( 1 , 2 ) ⊗ ( 2 , 2 ) with Δ k = 1 \Delta_{k=1} Δ k = 1 of (34) . Two properties make it suitable at continuous p p p : the construction involves no excited level that becomes null at p = 4 p=4 p = 4 , and its energy stays real along the whole flow to C p C_p C p (it crosses the lower non-Kac level of its sector near K ≃ 0.3 K\simeq0.3 K ≃ 0.3 without mixing, also at non-integer p p p ) and only acquires an imaginary part at K ≳ 3.5 K\gtrsim3.5 K ≳ 3.5 , far beyond the fixed point. The β \beta β -function B C B_C B C of (36) is then formed from w C w_C w C on successive size pairs and plotted against A C A_C A C exactly as in § 4 and § 5 .
Results at non-integer p: the family C p C_p C p and its endpoint
The motivation for following the flow at non-integer p p p is the contrast between the two integer results: at p = 5 p=5 p = 5 the flow ends on a genuinely conformal, partially transmitting interface C 5 C_5 C 5 , while at p = 4 p=4 p = 4 it ends on the invertible topological line η \eta η , which transmits everything. Two scenarios could connect them. In the first, C 5 C_5 C 5 belongs to a one-parameter family C p C_p C p of stable fixed points that moves continuously with p p p and simply arrives at η \eta η when p p p reaches 4. In the second, the stable interface C p C_p C p disappears before p p p reaches 4: it merges with another fixed point of the same ray—the familiar annihilation of a stable and an unstable fixed point, after which the pair continues into the complex plane [81 , 82 , 83 ] —and the flow out of D ( 2 , 1 ) D(2,1) D ( 2 , 1 ) , no longer stopped by C p C_p C p , runs on to a different stable fixed point, which at p = 4 p = 4 p = 4 is η \eta η . In the second scenario C 5 C_5 C 5 and η \eta η would be unrelated interfaces, the p = 4 p = 4 p = 4 endpoint would be an accident of the annihilation happening before p = 4 p = 4 p = 4 , and one would expect the familiar signatures of a merger: a slow, walking flow near the would-be fixed point for p p p just below the merger, and a β \beta β -function whose zero is approached with vanishing slope. The loop model lets us decide between the two, because it realises the same construction at every real p p p in [ 4 , 5 ] [4, 5] [ 4 , 5 ] and gives the β \beta β -function of § 6.2 at each of them. The result, presented in this section, is the first scenario: the zero of the β \beta β -function moves smoothly from C 5 C_5 C 5 to η \eta η , the slope at the zero grows in magnitude rather than vanishing, and no second fixed point appears on the ray. The merger with η \eta η at p = 4 p = 4 p = 4 is then of a different kind, discussed in § 7.2 .
C p C_p C p for 4 ≤ p ≤ 5 4 \leq p \leq5 4 ≤ p ≤ 5
Figures 23 and 24 show B C B_C B C for p = 5 , 4.75 , 4.5 , 4.25 , 4 p = 5, 4.75, 4.5, 4.25, 4 p = 5 , 4.75 , 4.5 , 4.25 , 4 from ED at L = 6 L = 6 L = 6 –14 14 14 (ℓ = 7 , 9 , 11 , 13 \ell= 7, 9, 11, 13 ℓ = 7 , 9 , 11 , 13 ; the L = 14 L = 14 L = 14 sector has 1.17 × 10 6 1.17 \times10^6 1.17 × 1 0 6 states and was run in steps of 0.1 up to just beyond the first zero), up to slightly beyond the first zero. At every p p p the curve leaves D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) with positive slope, turns over and crosses zero with negative slope, and the successive size pairs lie close to one another:
Table 1 (PDF p. 53).
The stable fixed point C p C_p C p therefore exists for every p p p between 5 and 4: it moves smoothly to larger bare coupling (K ∗ K^* K ∗ from 0.86 to 1.87) and to a larger Casimir value (A C ∗ A_C^* A C ∗ from 0.30 to 0.52), and nothing happens to it along the way. The p = 5 p = 5 p = 5 and p = 4 p = 4 p = 4 columns are the clean-ray counterparts of the improved-ray β \beta β -functions of Sections 4 and 5. At p = 4 p = 4 p = 4 the zero K ∗ = 1.87 K^* = 1.87 K ∗ = 1.87 , A C ∗ = 0.517 A_C^* = 0.517 A C ∗ = 0.517 matches the improved-ray result 1.90, 0.516, the counterterm there being small. At p = 5 p = 5 p = 5 the clean ray is the better behaved of the two: its three pairs lie close together, the steps of K ∗ K^* K ∗ (0.064, 0.044) decrease by a factor 0.7 from pair to pair, which is the basis of the estimate K C 5 ∗ = 0.95 ( 5 ) K_{C_5}^* = 0.95(5) K C 5 ∗ = 0.95 ( 5 ) of (39) , and the slope at the zero, −0.21, is close to the improved-ray value −0.16(2); the improved ray, by contrast, drifts much more because of the K K K -proportional counterterm (§ 5.4 ). The crossing value A C ∗ = 0.27 A_C^* = 0.27 A C ∗ = 0.27 –0.30 0.30 0.30 on the clean ray is also larger than the improved-ray values, the two being on their way to the common universal value with the slow L − 0.12 L^{-0.12} L − 0.12 corrections. Finite-size effects are stronger at large p. This is expected from the leading irrelevant dimension along the ray: the slope at the zero is Δ X − 1 \Delta_X - 1 Δ X − 1 , and it grows in magnitude steadily from −0.21 at p = 5 p = 5 p = 5 to − 0.64 -0.64 − 0.64 at p = 4 p = 4 p = 4 , so the L Δ X − 1 L^{\Delta_X-1} L Δ X − 1 corrections that make the p = 5 p = 5 p = 5 crossing drift are strongly suppressed at p = 4 p = 4 p = 4 .
The endpoint at p = 4 p = 4 p = 4 : the line D ( 1 , 4 ) D_{(1,4)} D ( 1 , 4 )
The family C p C_p C p ends, at p = 4 p = 4 p = 4 , on the invertible line η \eta η . In the Kac labelling of M ( p , p + 1 ) M(p,p+1) M ( p , p + 1 ) this line is η = D ( p − 1 , 1 ) ≅ D ( 1 , p ) \eta= D_{(p-1,1)} \cong D_{(1,p)} η = D ( p − 1 , 1 ) ≅ D ( 1 , p ) , so at p = 4 p = 4 p = 4 it is D ( 1 , 4 ) D_{(1,4)} D ( 1 , 4 ) . Unlike ( p − 1 , 1 ) (p-1,1) ( p − 1 , 1 ) , the label ( 1 , 4 ) (1,4) ( 1 , 4 ) makes sense for every p p p : the generalised minimal model at c = 1 − 6 / p ( p + 1 ) c = 1 - 6/p(p+1) c = 1 − 6/ p ( p + 1 ) [84 , 85 ] has topological lines for all degenerate labels ( r , s ) (r,s) ( r , s ) , with fusion rules that factorise into an r r r -part and an s s s -part and quantum dimensions d ( r , s ) = [ r ] q [ s ] q ′ d_{(r,s)} = [r]_q[s]_{q^\prime} d ( r , s ) = [ r ] q [ s ] q ′ , q = e i π / p q = e^{i\pi/p} q = e iπ / p , q ′ = e i π / ( p + 1 ) q^\prime= e^{i\pi/(p+1)} q ′ = e iπ / ( p + 1 ) . The line D ( 1 , 4 ) D_{(1,4)} D ( 1 , 4 ) therefore exists along the whole family, commutes with all the ( r , 1 ) (r,1) ( r , 1 ) lines (it lies in the other factor of the fusion category and is mutually transparent with them), and has
d ( 1 , 4 ) = sin ( 4 π / ( p + 1 ) ) sin ( π / ( p + 1 ) ) = 1.73 , 1.57 , 1.40 , 1.20 , 1 at p = 5 , 4.75 , 4.5 , 4.25 , 4 , (50) d_{(1,4)} = \frac{\sin\left(4\pi/(p+1)\right)}{\sin\left(\pi/(p+1)\right)} = 1.73,\ 1.57,\ 1.40,\ 1.20,\ 1 \quad\text{at } p = 5,\ 4.75,\ 4.5,\ 4.25,\ 4, \tag*{(50)} d ( 1 , 4 ) = sin ( π / ( p + 1 ) ) sin ( 4 π / ( p + 1 ) ) = 1.73 , 1.57 , 1.40 , 1.20 , 1 at p = 5 , 4.75 , 4.5 , 4.25 , 4 , (50)
so it is invertible only at p = 4 p = 4 p = 4 ; at p = 5 p = 5 p = 5 the Kac identification makes it D ( 1 , 4 ) ≅ D ( 4 , 2 ) D_{(1,4)} \cong D_{(4,2)} D ( 1 , 4 ) ≅ D ( 4 , 2 ) , a non-invertible line of dimension 3 \sqrt{3} 3 . For p > 4 p > 4 p > 4 , however, D ( 1 , 4 ) D_{(1,4)} D ( 1 , 4 ) cannot be the endpoint of the flow: as shown next, it carries a pair of relevant symmetric defect operators, so it is not infrared-stable, and the measured C p C_p C p (g C 5 = 1.503 g_{C_5} = 1.503 g C 5 = 1.503 , T = 0.665 \mathcal{T} = 0.665 T = 0.665 at p = 5 p = 5 p = 5 , against the topological values 3 \sqrt{3} 3 and 1 1 1 ) is a distinct, genuinely partially reflecting interface that merges with D ( 1 , 4 ) = η D_{(1,4)} = \eta D ( 1 , 4 ) = η only at p = 4 p = 4 p = 4 .
Defect operators on D ( 1 , 4 ) D_{(1,4)} D ( 1 , 4 ) . The operators on D ( 1 , 4 ) D_{(1,4)} D ( 1 , 4 ) are the twisted sector of D ( 1 , 4 ) ⊗ D ‾ ( 1 , 4 ) D_{(1,4)} \otimes\overline{D}_{(1,4)} D ( 1 , 4 ) ⊗ D ( 1 , 4 ) , i.e. the sum over the channels c ∈ ( 1 , 4 ) × ( 1 , 4 ) = ( 1 , 1 ) + ( 1 , 3 ) + ( 1 , 5 ) + ( 1 , 7 ) c \in(1,4) \times(1,4) = (1,1) + (1,3) + (1,5) + (1,7) c ∈ ( 1 , 4 ) × ( 1 , 4 ) = ( 1 , 1 ) + ( 1 , 3 ) + ( 1 , 5 ) + ( 1 , 7 ) of the c c c -twisted spaces; those neutral under the ( r , 1 ) (r,1) ( r , 1 ) lines (chiral label ( 1 , s ) (1,s) ( 1 , s ) with s s s odd, § 4 ) are the ( 1 , s ) ⊗ ( 1 , s ′ ) (1,s) \otimes(1,s^\prime) ( 1 , s ) ⊗ ( 1 , s ′ ) with s , s ′ s,s^\prime s , s ′ odd:
Table 7.
For generic p p p the line carries a pair of relevant symmetric operators, the chiral defect fields ϕ L = ϕ ( 1 , 3 ) ⊗ 1 ˉ \phi^L = \phi_{(1,3)} \otimes\bar{1} ϕ L = ϕ ( 1 , 3 ) ⊗ 1 ˉ and ϕ R = 1 ⊗ ϕ ˉ ( 1 , 3 ) \phi^R = 1 \otimes\bar{\phi}_{(1,3)} ϕ R = 1 ⊗ ϕ ˉ ( 1 , 3 ) of weights ( h 1 , 3 , 0 ) (h_{1,3},0) ( h 1 , 3 , 0 ) and ( 0 , h 1 , 3 ) (0,h_{1,3}) ( 0 , h 1 , 3 ) , Δ = ( p − 1 ) / ( p + 1 ) \Delta= (p-1)/(p+1) Δ = ( p − 1 ) / ( p + 1 ) . They live in the c = ( 1 , 3 ) c = (1,3) c = ( 1 , 3 ) channel, and that channel is exactly what the Kac quotient removes at p = 4 p = 4 p = 4 : there ( 1 , 4 ) ≅ ( 3 , 1 ) (1,4) \cong(3,1) ( 1 , 4 ) ≅ ( 3 , 1 ) and ( 3 , 1 ) × ( 3 , 1 ) = ( 1 , 1 ) (3,1) \times(3,1) = (1,1) ( 3 , 1 ) × ( 3 , 1 ) = ( 1 , 1 ) only, so the line becomes invertible and only the diagonal operators ( 1 , s ) ⊗ ( 1 , s ) (1,s) \otimes(1,s) ( 1 , s ) ⊗ ( 1 , s ) survive, the first being ( 1 , 3 ) ⊗ ( 1 , 3 ) (1,3) \otimes(1,3) ( 1 , 3 ) ⊗ ( 1 , 3 ) with Δ = 6 / 5 \Delta= 6/5 Δ = 6/5 , the leading invariant operator on η \eta η found in § 4 . For p > 4 p > 4 p > 4 the ( 1 , 3 ) (1,3) ( 1 , 3 ) channel survives, so D ( 1 , 4 ) D_{(1,4)} D ( 1 , 4 ) (at p = 5 p = 5 p = 5 , D ( 1 , 4 ) ≅ D ( 4 , 2 ) D_{(1,4)} \cong D_{(4,2)} D ( 1 , 4 ) ≅ D ( 4 , 2 ) ) is unstable there; this is why it cannot be the endpoint of the flow for p > 4 p > 4 p > 4 , as anticipated above.
The relevant fields ϕ L , R \phi^{L,R} ϕ L , R of D ( 1 , 4 ) D_{(1,4)} D ( 1 , 4 ) are null at p = 4 p = 4 p = 4 and the line D ( 1 , 4 ) = η D_{(1,4)} = \eta D ( 1 , 4 ) = η becomes infrared stable exactly at p = 4 p = 4 p = 4 . This is the analogue of the mechanism by which states of vanishing norm drop out of the O ( n ) O(n) O ( n ) and Potts models at integer n n n or Q \mathbb{Q} Q [80 , 86 ] .
Summary and discussion
The dilute A p − 1 A_{p-1} A p − 1 RSOS chain, read as an anyon chain of C p C_p C p with the S U ( 2 ) p − 2 SU(2)_{p-2} S U ( 2 ) p − 2 fusion rules, realises the ( r , 1 ) (r,1) ( r , 1 ) Verlinde defects of M ( p , p + 1 ) M(p,p+1) M ( p , p + 1 ) exactly topologically: a defect is one more anyon threaded through the chain, coupled by the F F F - and R R R -moves of the category, and the wrapped lines commute with the Hamiltonian to machine precision. The quantum-dimension-weighted projector combination (28) between the two fusion channels of the defect anyon and its neighbouring domain wall is the clean lattice representative of the KRW field ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) \phi_{(1,3)} \otimes\bar{\phi}_{(3,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 3 , 3 ) : its overlaps with the identity and the irrelevant primary ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) \phi_{(1,3)} \otimes\bar{\phi}_{(1,3)} ϕ ( 1 , 3 ) ⊗ ϕ ˉ ( 1 , 3 ) vanish to numerical precision at the tested sizes, while descendants and higher fields can contribute irrelevant corrections.
With the clean operator the KRW flow is reproduced from first principles. At p = 4 p = 4 p = 4 , it ends on the invertible line η \eta η , an endpoint that had not been proposed before: the g g g -function drops from 2 \sqrt{2} 2 to 1.00, every level of the two-strand sector-resolved spectrum flows to the diagonal tricritical-Ising value, a ring with a single strand reproduces the Z 2 \mathbb{Z}_2 Z 2 -twisted spectrum with no vacuum, which rules out the trivial line as the endpoint. At p = 5 p = 5 p = 5 it ends on the conformal interface C 5 C_5 C 5 with g = 1.503 ( 2 ) g = 1.503(2) g = 1.503 ( 2 ) and T = 0.665 ( 10 ) \mathcal{T} = 0.665(10) T = 0.665 ( 10 ) , whose leading irrelevant field has Δ X = 1.12 ( 1 ) \Delta_X = 1.12(1) Δ X = 1.12 ( 1 ) . Gauging η \eta η maps C 5 C_5 C 5 to an interface of the three-state Potts model, τ F i b → C 5 P o t t s \tau_{\mathrm{Fib}} \to C^{\mathrm{Potts}}_5 τ Fib → C 5 Potts , which breaks Z 3 \mathbb{Z}_3 Z 3 and is stable precisely when the Fibonacci line τ F i b \tau_{\mathrm{Fib}} τ Fib and charge conjugation C C C are both imposed (§ 5.5 ). Followed at continuous p p p in the loop model (§ 6 ), the family C p C_p C p moves smoothly from C 5 C_5 C 5 to η \eta η , and the merger at p = 4 p = 4 p = 4 is not a fixed-point annihilation: C p C_p C p is the endpoint of the symmetric ϕ ( 1 , 3 ) \phi_{(1,3)} ϕ ( 1 , 3 ) deformation of the line D ( 1 , 4 ) D_{(1,4)} D ( 1 , 4 ) , whose relevant ( 1 , 3 ) (1,3) ( 1 , 3 ) -channel fields become null at p = 4 p = 4 p = 4 and decouple, so that η = D ( 1 , 4 ) \eta= D_{(1,4)} η = D ( 1 , 4 ) becomes stable exactly there (§ 7.2 ).
Figure 24. The ℓ = ( 11 , 13 ) \ell=(11,13) ℓ = ( 11 , 13 ) curves of Figure 23 (L = 10 , 12 , 14 L=10,12,14 L = 10 , 12 , 14 ) overlaid: the stable fixed point C p C_p C p moves continuously from p = 5 p=5 p = 5 to p = 4 p=4 p = 4 .
Figure 23. Loop model on the clean ray H 0 + K Φ c l H_0 + K\Phi_{\mathrm{cl}} H 0 + K Φ cl : β \beta β -function B C B_C B C of the defect Casimir coupling (35) versus A C A_C A C for p = 5 , 4.75 , 4.5 , 4.25 , 4 p = 5, 4.75, 4.5, 4.25, 4 p = 5 , 4.75 , 4.5 , 4.25 , 4 (ED, L = 6 L = 6 L = 6 –14, size pairs ℓ = ( 7 , 9 ) \ell= (7, 9) ℓ = ( 7 , 9 ) , ( 9 , 11 ) (9, 11) ( 9 , 11 ) and ( 11 , 13 ) (11, 13) ( 11 , 13 ) ), drawn up to 1.35 K ∗ 1.35 K^* 1.35 K ∗ ; the black square is D ( 2 , 1 ) D_{(2,1)} D ( 2 , 1 ) , the filled circles the zero of each pair. The last panel collects K ∗ K^* K ∗ (all pairs) and A C ∗ A_C^* A C ∗ (largest pair) versus p p p .
Comparison with conformal perturbation theory and with integrable defect flows. The transmission of C C C has been computed in conformal perturbation theory by Makabe and Watts [12 ] : for the diagonal KRW fixed point they find the reflection coefficient R C = 1 − T = 8 9 π 2 y 4 + O ( y 5 ) R_C=1-T=\frac{8}{9}\pi^2y^4+O(y^5) R C = 1 − T = 9 8 π 2 y 4 + O ( y 5 ) with y = 1 − h 1 , 3 = 2 / ( p + 1 ) y=1-h_{1,3}=2/(p+1) y = 1 − h 1 , 3 = 2/ ( p + 1 ) , which gives R = 0.284 R=0.284 R = 0.284 at p = 4 p=4 p = 4 (against 0.366 0.366 0.366 from their non-perturbative work on tricritical Ising defects [11 ] ) and R = π 2 / 72 = 0.137 R=\pi^2/72=0.137 R = π 2 /72 = 0.137 at p = 5 p=5 p = 5 , to be compared with our R = 1 − T = 0.335 ( 10 ) R=1-T=0.335(10) R = 1 − T = 0.335 ( 10 ) . As for the g g g -function, the leading term underestimates the departure from the topological value at y = 1 / 3 y=1/3 y = 1/3 , by a factor 2.4 here; the general leading-order formulas for the change of g g g , of the reflectivity and of the entanglement entropy along a defect flow driven by a spinless field are given in [13 ] . On the lattice side, the chiral KRW flows—the axes of the ( λ l , λ r ) (\lambda_l,\lambda_r) ( λ l , λ r ) plane—have been realised in the dense RSOS models by an inhomogeneous spectral parameter on one column, which keeps the model integrable and gives the running g g g -function from the thermodynamic Bethe ansatz [23 , 87 ] ; those flows end on topological lines, whose identity is now fixed in general by the functional-relation analysis of [27 ] , and the present construction is the complementary one, in which the non-chiral deformation of an exactly topological seam reaches the non-topological interface C C C .
Comparison with open string field theory. The g g g -function of the KRW defect C C C has also been computed by Kudrna, Maccaferri and Schnabl [14 ] as a level-truncated solution of open string field theory in the doubled theory M ( p , p + 1 ) ⊗ M ( p , p + 1 ) M(p,p+1)\otimes M(p,p+1) M ( p , p + 1 ) ⊗ M ( p , p + 1 ) . Their extrapolated values are
g C O S F T ≃ 1.08 – 1.10 ( p = 4 ) , g C O S F T ≃ 1.489 ( p = 5 ) , g_C^{\mathrm{OSFT}} \simeq1.08\text{--}1.10 \quad(p=4), \qquad g_C^{\mathrm{OSFT}} \simeq1.489 \quad(p=5), g C OSFT ≃ 1.08 – 1.10 ( p = 4 ) , g C OSFT ≃ 1.489 ( p = 5 ) ,
with the p = 4 p=4 p = 4 number flagged there as requiring higher levels, and the leading conformal-perturbation-theory estimate quoted in the same work, ln g C − ln g D ( 2 , 1 ) = − π 2 4 y 3 + O ( y 4 ) \ln g_C-\ln g_{D(2,1)}=-\frac{\pi^2}{4}y^3+O(y^4) ln g C − ln g D ( 2 , 1 ) = − 4 π 2 y 3 + O ( y 4 ) with y = 2 p + 1 y=\frac{2}{p+1} y = p + 1 2 , gives 1.477 at p = 5 p=5 p = 5 . These numbers are comparable with our lattice results, g = 1.503 ( 2 ) g=1.503(2) g = 1.503 ( 2 ) for C 5 C_5 C 5 and g = 1.00 g=1.00 g = 1.00 for the p = 4 p=4 p = 4 endpoint.
Acknowledgments
This work was heavily assisted by Claude Fable 5 (Anthropic): the lattice constructions, the numerical codes, the computations reported here and the drafting of the text were developed in collaboration with it, with the author directing the questions and checking the results. The work of the author is supported by the European Union (ERC, QFTinAdS, project number 101087025, PI Balt van Rees). Views and opinions expressed are, however, those of the authors only and do not necessarily reflect those of the European Union or the European Research Council Executive Agency.