Entanglement Entropy in Integrable Field Theories with Line Defects I. Topological Defect

In this paper and a companion one, we study the effect of integrable line defects on entanglement entropy in massive integrable field theories in 1+1 dimensions. The current paper focuses on topological defects that are purely transmissive. Using the form factor bootstrap method, we show that topological defects do not affect the the entanglement entropy in the UV limit and modify slightly the leading exponential correction in the IR. This conclusion holds for both unitary and non-unitary field theories. In contrast, non-topological defects affect the entanglement entropy more significantly both in UV and IR limit and will be studied in the companion paper.


Introduction
Quantum entanglement in many body systems is one of the most fundamental and fascinating phenomena in nature. It is therefore very important to quantify (define proper quantities such as entanglement entropy, Rényi entropy, mutual information ect. to quantify quantum entanglement), calculate (develop analytical and numerical methods to compute the quantities defined in the previous step for different systems) and measure (in experiments) quantum entanglement. One most commonly studied quantity for quantum entanglement is the entanglement entropy. Computing entanglement entropy in many body systems is important as well as challenging for theories and also interesting for experiments (see for example [2,3]). In the past decade, impressive progress has been made in computing entanglement entropy, especially in conformal field theories where powerful geometric methods [4,5] and holographic insights [6,7] play a crucial role. On the other hand, computation of entanglement entropy in generic massive quantum field theories remains open. Among them integrable field theories in 1+1 dimensions is special. Due to integrability, they are more tractable analytically than generic field theories. In a series of works [8][9][10][11][12], the authors studied entanglement entropy in the context of integrable field theories using the form factor bootstrap approach. The method is based on a generalized version of the traditional form factor bootstrap program [13,14] in the context of replica theory.
The form factor method is non-perturbative and is especially powerful in the IR limit where the result can be approximated to high precision by the first few terms in the spectral expansion. In this way, the authors obtained the leading corrections in the IR for the entanglement entropy, which takes a remarkably universal form that only depends on the spectrum of the theory [8].
On the other hand, higher order corrections depend on more details of the theory.
Given an integrable field theory, one can introduce additional structures such as boundaries and defects. In general, these additional structures will spoil the integrability of the theory. However, there exist special kinds of boundaries and defects that preserve integrability [15][16][17][18] which allow one to apply powerful integrability techniques. In [10], the authors investigated Ising model with integrable boundary conditions. The main goal of the current work is to study the case of integrable line defects in integrable field theories.
There are at least two reasons to consider line defects in quantum field theories. Firstly, homogeneous systems are usually mathematical idealization and it is important to estimate the effect of defects on entanglement entropy. This question has been investigated extensively for spin models and conformal field theories, see for example [19][20][21][22][23][24][25] and references therein. Secondly, defects often contain useful information about the structure of the bulk theory and thus classification and the study of the defects are of considerable theoretical interest [18,[26][27][28][29][30][31][32][33][34].
The requirement of integrability imposes very stringent constraints on the possible defects. The situation for free theories with S = ±1 and for interacting theories, namely those with non-trivial S-matrices, are quite different.
For interacting integrable field theories, the only allowed defects are either purely reflective (i.e. integrable boundaries) or purely transmissive (i.e. topological defects) [17,35]. For the latter case, the defect is topological in the sense that the stress tensor is continuous across the defect and one can move the defect freely without costing any energy as long as it does not cross any other operators. This kind of defect is simply characterized by a transmission amplitude T ij (θ). Form factors of a local operator in interacting theories can usually be worked out in practice for states with few particles and closed formulae for states with any number of particles are not known in general. Therefore, for interacting theories, we restrict ourselves to the leading corrections to the entanglement entropy [8] in the IR and show how the defect affects the universal correction of the bulk entanglement entropy for both unitary and nonunitary theories.
For free theories, integrable defects can be transmissive and reflective simultaneously. This kind of defect is characterized by a transmission amplitude T ij (θ) and a reflection amplitude R ij (θ) and is non-topological. On the one hand, there exist closed formula for the form factors of local operators (e.g. branch-point twist fields) for any number of particles [8,36,37] due to Wick's theorem. Given a closed formula for form factors of branch-point twist fields, we can expect more than just finding leading corrections in the IR. One can also consider the UV limit and make contact with the underlying CFT. In this way, one can obtain the boundary or defect entropy [38] which characterize the number of degrees of freedoms on the boundary or defect.
The effects of topological and non-topological defects on entanglement entropy are also quite different. Intuitively, topological defects are very "soft" and one does not expect them to modify the bulk entanglement entropy drastically. This is indeed the case. In the UV limit where the main contribution comes from particles with high energies, the defect becomes transparent and the UV behavior is not modified. In the IR limit, the effect of the defect enters the leading exponential correction via transmission amplitude, which also takes a universal form Bulk EE: whereT (θ, α) = T ( iπ 2 − θ, α) is the Wick rotated transmission amplitude. Here for simplicity, we give the result for theories whose spectrum consist a single particle. The result can be generalized straightforwardly to theories with more kind of particles. The leading correction takes a universal form which depends on both the spectrum of the theory and the transmission amplitude that characterize the topological defect. The transmission amplitude depends on some parameter which we denote α. This parameter can be related to the coupling constant in the Lagrangian formulation. We find that for certain range of the parameter α, the contribution from topological defect is finite while for the parameters outside this range, the effect of defect can be completely neglected also in the IR limit.
For non-unitary theories like the scaling Lee-Yang model, the calculation of entanglement entropy is much more subtle due to non-unitarity [12,39]. Nevertheless, we have the same conclusion that the UV limit is unmodified while the IR corrections are modified for some range of the parameter that characterizes the defect. The explicit expression is slightly more involved and is given in (6.31).
For non-topological defects in the free theories, on the other hand, the entanglement entropy is modified more significantly both in the UV and IR. The analysis for these defects are more technical and the results will be given in the companion paper [1].
The rest of the paper is structured as the follows. In section 2 we introduce integrable line defects in integrable field theories in 1+1 dimensions. In section 3 we review the method of form factor bootstrap in the context of entanglement entropy and generalize it to the defect case. In section 5 and 6, we consider the leading corrections to the entanglement entropy in unitary and non-unitary field theories using Sinh-Gordon theory and Scaling Lee-Yang theory as main examples, respectively. Finally we conclude in section 7.

Line defects in integrable field theories
In this section, we review line defects in integrable field theories. Our discussion mainly follows the seminal papers [16,17] and [18,40].

Line defects
We consider an infinitely long line defect in a 1+1 dimensional integrable quantum field theory (IQFT). In the Lagrangian description, the defect can be introduced by adding a term which has support only on the defect line where S B is the bulk action and L D is the defect lagrangian. We further require the interaction between bulk and defect be consistent with the existence of infinitely many conserved charges. The defect of this type preserves integrability and is referred to as an integrable defect.
In the bootstrap framework, the bulk physics is characterized by elastic scattering processes of bulk particles. In the presence of the defect, we also take into account scattering of bulk particles on the defect. Scattering processes on the defect is necessarily elastic due to integrability. A particle which scatters on the defect can either go through without changing its energy and momentum, or it can be reflected which preserves energy but flips the sign of momentum. It has been shown [17,35] that the only theories which allows both transmissive and reflective processes are the generalized free models (free boson and free fermion) with bulk scattering matrix S = ±1. For interacting integrable field theories, the defect has to be either purely reflective or purely transmissive. The first case is nothing but the integrable boundary, which has been investigated intensely starting from the seminal work of Ghoshal and Zamolodchikov [15]. The latter case is called topological defect since the stress energy tensor is continuous across the defect. Topological defects have many special properties. They can be moved freely without changing physical observables and their positions are not important as long as they do not cross other local operators. One can also fuse a topological defect with integrable boundaries to obtain new integrable boundary conditions [40]. Topological defects of integrable field theories can be obtained as chiral perturbation of topological defects of the underlying conformal field theory [34,41].

Defect algebra
It is convenient to represent the defect formally as an operator D [17]. In an IQFT, excitations can be described by the so-called Zamolodchikov-Faddeev (ZF) operators A i (θ) and A † i (θ) which satisfy a set of algebraic relations (ZF algebra) where θ 1,2 are the rapidity of the particles and the index "i, j, k, l" denote quantum numbers of the excitations. Here S kl ij (θ 1 − θ 2 ) is the bulk scattering matrix which characterize the bulk dynamics. The consistency of ZF algebra leads to Yang-Baxter equations and unitarity conditions for the scattering matrix. Similarly, the interaction between the defect and bulk excitations can be described by an algebra between the defect operator D and the ZF operators A † i (θ). Assuming the defect has no internal degree of freedom, the defect algebra reads where R ± ij , T ± ij are the reflection and transmission matrices respectively. Similar to the bulk case, the consistency of the defect algebra gives rise to constraints for R ± ij and T ± ij [17,35].

Correlation functions
Using the replica trick, the computation of entanglement entropy can be recast in terms of correlation functions of branch-point twist fields in the replica theory. In the bulk case, correlation functions of local operators can be computed by the form factor approach [13,14]. The main idea is to write correlation functions in terms of form factors via spectral expansion. The form factors are determined by solving a set of functional equations called the bootstrap axioms.
In the presence of a line defect, if the local operator does not sit on the defect, the form factor of the operator are essentially the same as the bulk case [18]. On the other hand, the presence of the defect affects the spectral expansion of correlation functions and gives rise to extra factors which are matrix elements of the defect operator.
In order to see this, it is more convenient to perform a "double Wick rotation" which exchanges temporal and spatial directions [17,42]. Originally, the defect is placed at x = 0 which is localized in space but extend in time. After the Wick rotation, we place the defect at t = t D . Then the defect can be taken into account by defining an extended operator D 1 placed at t = t D acting on the bulk states. This is similar to the Wick rotation in the boundary integrable theories which defines a boundary state |B [15]. The correlation functions can be expressed as Here we put a "D" on the l.h.s to remind that it is the correlation functions in a defect theory. The fields on the l.h.s. are the fields in Heisenberg representation whose evolution is governed by the complete Hamiltonian (bulk+boundary). On the r.h.s. the operators O i (x i , t i ) are the bulk fields whose time evolution is governed by the bulk Hamiltonian.
Similar to the boundary case, we can represent the operator D in terms of ZF operators [42,43] where summation over repeated indices is understood.
After the Wick rotation we can perform the spectral expansion on the r.h.s. and apply the form factor bootstrap approach. For our purpose, the two-point function is most relevant. We consider the case where the two operators resides on two sides of the defect with x 1 = x 2 = 0 and t 1 < t D < t 2 . Schematically, where r L = |t 1 −t D | and r R = |t 2 −t D |. F O i N are the N -particle form factor of operator O i , which can be determined by the form factor bootstrap program. Note that here we use a shorthand notation for the resolution of identity The effect of the defect comes in through the matrix element D M,N .
If we take a trivial defect D M,N = δ M,N in (2.7), we obtain the usual form factor expansion in the bulk theory. If we assume the defect to be purely reflective, then the non-vanishing matrix elements are D 2M,2N = D 2M,0 D 0,2N . Note that is nothing but the Ghoshal-Zamolodchikov boundary state. The spectral expansion of the two-point function in this case reduces to the product of two boundary expansions.
The matrix elements D 2N,0 and D 0,2N only depend on the reflection amplitude R ij . If the defect is purely transmissive, the non-vanishing elements are D N,N and hence the spectral expansion becomes where r = |t 1 − t 2 | = r L + r R and the matrix element D N,N only depends on the transmission amplitude T ij .

Entanglement entropy in the presence of a line defect
In this section, we formulate our set-up of computing the bipartite entanglement entropy in the presence of an integrable line defect. From now on, we focus on topological defects. We consider a line interval A of length r and denote its complement asĀ. There are three possibilities for the location of the line defect, as is shown in figure 3.1. In this paper, we consider the case (b) of figure 3.1. For case (a), the defect is outside the finite interval A and can be moved to infinity since it is topological, which reduces to the bulk case. The case (c) is more subtle and will be considered elsewhere. Let us denote the Hilbert space by H = H A ⊗ HĀ. For a given pure state |ψ , which in our case is the defect vacuum state |0 D , the reduced density matrix ρ A is given by and the entanglement entropy is defined as the von Neumann entropy

Replica trick
We apply the replica trick to compute the entanglement entropy [8]. One first compute Tr ρ n A for positive integer n and then analytically continue n to any value. The entanglement entropy is obtained the as the following limit In order to compute Tr ρ n A , we consider n copies of the original model and glue them cyclicly, as is shown in figure 3.2. This amounts to compute the partition function on the multi-sheeted Riemann surface with defects. The theory defined on this Riemann surface is however not local. Locality can be restored by introducing a multi-copy model formed by n identical copies of the original model with specific boundary conditions for the fields. The problem is thus reduced to the computation of the partition function for the multi-copy model on the complex plane C with proper boundary conditions for the fields. When applying the replica trick, the defect is also duplicated and present on each sheet of the Riemann surface, which we denote by D µ , µ = 1, 2, · · · , n. As a result, in the multi-copy model on the complex plane, the defect is the  fused defect from all the n copies, namely D = D 1 D 2 · · · D n . This is depicted in figure 3.3. It has been shown that the entanglement entropy in the replica theory can be computed as the two-point function of certain special local operators called branch point twist fields [8]. Note that the prescription is different for non-unitary theories which will be given in section 6. Using the branch point twist fields, we can write Here T is the branch-point twist field andT = T † is its hermitian conjugate. The constant Z n is an n-dependent non-universal constant which satisfies Z 1 = 1. The short-distance cut-off ε is chosen in such a way that ∂ n Z n | n=1 = 0. ∆ n is the conformal dimension of the branch point twist field The insertion of a topological defect does not modify the form of the two-point function in the UV limit, namely where we recall that This can be proven within the framework of form factor method, which will be present in section 4. Plugging (3.6) into (3.4) one finds that the r.h.s. is dimensionless and the entanglement entropy in the UV limit is given by Therefore the entanglement entropy is not modified by the topological defect in the UV limit. Let us denote the mass of the lightest particle of the IQFT by m and write where we have chosen the normalization such that D = 1. The entanglement entropy in the IR limit where r L , r R → ∞ can then be written as In the limit n → 1, the three parts are given by The derivatives of the first two terms are where U is a model dependent constant. Therefore the entanglement entropy takes the form in the IR limit In the following, we will compute the leading corrections for different models. The effect of topological defect will enter via the last term of (3.13).

Spectral expansion and form factors of branch-point twist fields
The main calculation of the entanglement entropy now boils down to the computation of the two-point function of the branch-point twist fields with defect operator 0|T (t 1 )DT (t 2 )|0 .
Performing the spectral expansion we can write Here the form factors F T n and FT n are the same as the bulk case and have been derived in [8]. The factors D M,N are matrix elements of the defect operator where the asymptotic states are given by and for simplicity we omit the replica indices in (3.16) and in what follows.
For the interacting theories, the defect is purely transmissive. For simplicity, we consider the theory with one type of particle such as Sinh-Gordon and assume that the defect is parity invariant and does not have internal degrees of freedom. In this case, the bulk interaction is characterized by a simple S-matrix S(θ) and the defect is characterized by a transmission amplitude T (θ). The only non-vanishing matrix elements are where µ i , ν i = 1, · · · , n are replica indices of particles with rapidities θ i , θ i respectively. Here 'permutations' means that we need to take into account all possible contractions between the incoming and outgoing particles. The contraction between particles {θ i , µ i } and {θ j , ν j } contributes a factorT (θ i )δ µ i ,ν j δ(θ i −θ j ). When the order of the particles are exchanged, we pick up the corresponding S-matrix in the replica theory [S(θ)] δµ i ,µ j since only the particles with the same replica index have non-trivial interactions. For example, the first two defect matrix elements are given by Plugging the defect matrix element (3.18) into the spectral expansion (3.15), we find that As we see the effect of purely transmissive defect is simply modifying the measure of the integral by multiplication of transmission factors and only depends on r instead of r L , r R . This confirms our claim that we can move topological defects freely within the interval. The same result can also be obtained from a slightly different point of view from the defect form factor axioms [18].

UV limit
In this section, we study the short distance behavior of the defect two-point function T (t 1 )DT (t 2 ) and show that it is the same as T (t 1 )T (t 2 ) for topological defects, namely lim r L ,r R →0 This implies that the UV limit of the entanglement entropy is not modified by topological defects.
Purely transmissive defects can be obtained from topological defects of the underlying CFT by chiral perturbation [34,41]. It has been shown in [44,45] that topological defects in CFT do not change the universal behavior of entanglement entropy. This can also be analyzed in the framework of form factor approach. On the other hand, if the defect is non-topological, then it will in general modify the UV behavior.
In order to analyze the short distance behavior of two-point functions in the form factor framework, it is useful to consider the form factor expansion of the logarithm of the two-point function [12,46,47] log Similarly, the two-point function with a purely transmissive defect has the following expansion where the function H O|µ 1 ···µ k k (θ 1 , · · · , θ n ) can be seen as the "connected parts" of the form factors and can be worked out explicitly for any k. The first few H k 's are given by [12] By studying the limit mr → 0 one can obtain the scaling dimension ∆ O of the operator O, which is given by [12,46,47] Now we analyze the same limit for the expansion with defect (4.3). We shift the integration variable θ j → θ j + θ k (j = 1, · · · , k − 1). Due to the relativistic invariance of the bulk form factor, we have = ξ cosh(θ k + τ ) 2 Note that here we assume the theory is unitary and the expansion is slightly different from that in [12] which is for the non-unitary case up to some signs. we can write the exponential factor in the second line of (4.6) as e −mr cosh(θ k +τ ) . This factor is almost 1 in the interval log(mr) < θ k + τ < − log(mr) and zero outside the interval when mr → 0, as is shown in figure 4.4. Therefore the small mr limit is determined by the large θ k behavior of the integral. Since θ k appears in every transmission amplitudeT (θ j + θ k ) and T (θ k ), we thus need to analyze the large rapidity behavior for the transmission matrixT (θ).
A purely transmissive defect is characterized by two transmission factors T ± (θ) (for parity invariant defect we have T + (θ) = T − (θ)) which describe the scattering of particles with defect from left and right, respectively. The two transmission factors are not independent and are related by defect unitarity and crossing symmetry The quantity that appear in the spectral expansion is given byT (θ) = T − ( iπ 2 − θ) = T + ( iπ 2 + θ). The first relation in (4.10) can be written asT (θ)T (θ − iπ) = 1. When θ is very large, the shift of iπ can be neglected and we can write lim θ→±∞T (θ) 2 = 1 (4.11) Therefore asymptoticallyT (θ) → ±1 by unitarity. Now let us specify to the models we study in this paper. For scaling Lee-Yang model, the transmission amplitude is given in (6.2) with the asymptotic behavior lim θ→±∞TsLY (θ) = 1. Therefore the defect is transparent in the large rapidity limit and does not modify the scaling dimension. For Sinh-Gordon model, the transmission amplitude is given in (5.2) with the asymptotic behavior lim θ→±∞TShG (θ) = ∓1. This might create some problem for the terms involve odd number of particles since the asymptotics of two limits cancel each other. However, Sinh-Gordon theory is a parity invariant theory and the operators can be classified according to the parity. For parity even (odd) operators, only terms with even (odd) number of particles contribute to the spectral expansion [48]. This property is still true in the replica theory and the branch-point twist operator is parity even, so the terms with odd number of particles are automatically zero [49]. For the terms with even number of particles, the products of transmission amplitudes are transparent in the large rapidity limit. To conclude, topological defects do not change the two-point function of branch-point twist fields in the UV and the entanglement entropy stays the same.

Unitary case : Sinh-Gordon theory
In this section, we consider the leading correction to the entanglement entropy for unitary IQFT. We take the Sinh-Gordon theory as our example. There is only one type of particle and the scattering process is characterized by the following bulk scattering matrix 3 where β is the coupling constant in the Sinh-Gordon lagrangian. The transmission amplitude is given by (see for example [40]) where κ is a parameter that characterize the defect. In the Sinh-Gordon model, the leading contribution starts with D 2,2 since the one-particle form factor is vanishing. Plugging (3.19) into the spectral expansion (3.14) and (3.15), we obtain the leading contribution 3 Note that the convention of the scattering matrix [8] is slightly different from those in [40]. We follow the convention of [8] in this paper.
where in the second line we have used the property of the two-particle form factor F T |i i+k (θ 1 , θ 2 ) = F T |j j+k (θ 1 , θ 2 ) to perform the summation over µ 2 . The explicit form of the two-particle form factors of branch-point twist fields can be found in [8]. For the leading correction, we do not need the explicit form. It is sufficient to apply the following important property of the twoparticle form factor of the branch-point twist field after we analytically continue n according to the prescription given in [8] ∂ ∂n Using (5.4), the leading correction to entanglement entropy reads If we setT (θ) → 1, we recover the leading correction to the bulk entanglement entropy. Note that our result for the leading exponential correction (5.5) is also universal in the sense that it only depends on the spectrum of the theory and the transmission amplitude of the defect. The contribution due to the presence of the defect can be extract to be For Sinh-Gordon model, we plug in the explicit transmission amplitude and obtain This quantity depends on the parameter Bκ. We plot the function f exp (θ) = e −2mr cosh θ and f D (θ, Bκ) = 1/(cosh 1 2 (θ − Bκ)) 2 with different values of Bκ in figure 5.5. We see that both functions have similar behavior, namely they are finite in some range and quickly damped away outside this range. The parameter Bκ controls the center of the "peak". If the finite range of f exp (θ) and f D (θ, Bκ) do not overlap, the product is damped and gives almost zero. On the other hand, if the finite range of the two functions overlap, then the contribution is small but finite.
In order to have some idea about the finite range, we perform some numeric analysis in what follows. Let us consider the overlap of two functions f exp (θ) and f D (θ), which is the colored region in figure 5.6. To be more explicit, we consider f exp (θ) = e −2mr cosh θ , f D (θ) = 1 cosh( 1 2 (θ − λ mr)) 2 (5.8)  where we have set Bκ = λ mr, namely we use mr as the unit for the rapidity θ. The value of λ determines the position of the peak, as is shown in figure 5.6. We denote the maximal value of f exp (θ) as I max = e −2mr and set it as the unit for the value of the two functions. The maximal overlap is obtained at λ = 0. As |λ| increases, the overlap decreases. For large enough |λ|, the curves of two functions intersect at one point 4 which we denote by θ = θ * . We denote the value of the two functions at this point by f exp (θ * ) = f D (θ * ) = η * I max . Both θ * and η * depend on λ and hence θ * = θ * (λ) and η * = η * (λ). Let us fix a small value for η * , say η 0 = 1.0 × 10 −5 and we consider the overlap to be negligible for 0 < η * < η 0 and non-negligible otherwise. The requirement for non-negligible overlap 1 > η * > η 0 fixes a finite range for |λ| < λ 0 where λ 0 > 0 and satisfies η * (λ 0 ) = η * 0 . The analytic relation η * (λ) is the solution of the following equation and is hard to write down explicitly. However, for a fixed η 0 and mr it is straightforward to find the corresponding λ 0 numerically. In the following table, we list a set of values of λ 0 for different values of mr. In 1+1 dimensions, non-unitary field theories make perfect physical sense and describe many statistical mechanics system. The computation of entanglement entropy for non-unitary IQFT's such as the scaling Lee-Yang model is more complicated than the unitary ones. First of all, the partition function on the Riemann surface is no longer given by the two-point function of branch-point twist fields. One needs to introduce a new kind of twist field which is defined by the operator product of the branch-point twist field T and the scalar field with smallest scaling dimension denoted by φ. The partition function on the Riemann surface is then given by the ratio of these two kinds of two-point functions. Secondly, non-unitarity also affects spectral expansion since many fields become non-Hermitian. The computation of entanglement entropy in non-unitary IQFT's has been investigated in [12]. In this section, we consider non-unitary IQFT's with integrable topological defect using scaling Lee-Yang model as an example.

Defect scaling Lee-Yang model
The scaling Lee-Yang model can be obtained as a relevant perturbation of the CFT minimal model M (2,5) with central charge c = −22/5. This model has only one type of neutral particle φ with mass m and the bulk scattering is characterized by the following S-matrix S(θ) = sinh θ + i sin π 3 sinh θ − i sin π 3 (6.1) In this model, one can introduce the topological defect which preserves integrability with the transmission amplitude [40] T 2) The defect model and the corresponding form factor bootstrap program has been studied extensively, see for example [18,34,40].

Replica trick and entanglement entropy
As in the unitary case, one can apply the replica trick to compute the entanglement entropy. However, the computation is more subtle due to non-unitarity. For non-unitary CFTs like Lee-Yang model, the physical vacuum state (i.e. lowest energy state) does not coincide with conformal vacuum. Due to state-operator correspondence, we can associate the physical vacuum with a local field φ with scaling dimension ∆. As expected, this field will enter the computation of EE. It is shown by careful analysis in [12,39] that the partition function on the Riemann surface for non-unitary theories is now given by where we denote the scaling dimension of the branch-point twist field T as ∆ T . The normalization φ(x 1 )φ(x 2 ) n in the denominator can be seen as the norm of the vacuum state in CFT.
The new twist field T φ is defined as We denote the scaling dimension of this new branch-point twist field as ∆ T φ which reads From (6.3) we see that in order to compute the entanglement entropy, we need to evaluate two types of two-point functions.

Form factors and defect matrix elements
In this subsection, we give the ingredient that are necessary for the the computation of the entanglement entropy. We need form factors of the twist field T φ and the fundamental field φ up to two particles. The form factor of field φ has been derived by Zamolodchikov [50].
Form factors of fundamental field φ The one-particle form factor is a constant due to relativistic invariance and is given by and the two particle form factor can be written as Here the two-point minimal form factor F min (θ, n) is given by where a(θ, n) encodes the pole for bound states a(θ, n) = cosh θ n − 1 cosh θ n − cosh 2π 3n (6.10) and f (θ, n) is given by t sinh(nt) cosh t 2 cosh t n + iθ π dt (6.11) Form factors of twist field T φ andT φ The one-particle form factor for the branch-point twist field is given by 12) The two-particle form factor is given by Form factors with generic replica indices can be expressed through the fundamental one by , · · · , θ k + 2πi(µ k − 1)) (6.14) From the inversion of copy numbers, we can related the form factors of T φ andT φ Non-unitarity of the theory modifies complex conjugation of the form factors by an extra phase factor Matrix element of defect D In the defect theory, we also need the matrix elements of the defect operator, which in our case is given by the transmission matrix (6.5). In the computation of entanglement entropy, we needT (θ) = T − ( iπ 2 − θ). The defect depends on the parameter b. Here we consider a special case b = 3 + i α (α ∈ R) where the transmission amplitudeT (θ) is real and given by 5 Our motivation for such a choice is that for generic parameter b, the correction to the entanglement entropy will be complex which has no obvious physical meaning. The one-and two-particle defect matrix element are given by the same formulas as in (3.19). The main difference is that now both D 11 and D 22 will contribute to the entanglement entropy.

Spectral expansion and entanglement entropy
In this subsection, we compute the entanglement entropy using the form factor bootstrap method. Slightly generalizing the expression in [12], we can write the entanglement entropy in terms of two-point functions as and r = |t 1 − t 2 |. The entanglement entropy in the UV limit is not modified by the topological defect and is given by [39] S A = c eff 3 log(r/ε) (6.20) where for scaling Lee-Yang model c eff = c − 24∆ = 4/5. In the IR limit, the entanglement entropy (6.18) can be written as where U is defined in [12]. We concentrate on the leading correction which is given by where A D (r, n) = dA D (r, n)/dn. Both A D (r, n) and B D (r) can be expanded in terms of form factors. We denote A D,k (r, n) and B D,k (r) the contribution of the k-particle state and A D (r, n) = A D,1 (r, n) + A D,2 (r, n) + · · · (6.23) Keeping the r.h.s of (6.22) up to two particle contributions, we obtain Now we compute the above quantities (θ)e −mr cosh θ (6.25) Note that the factor F contains n in a non-trivial way, so that A D,1 (r, 1) needs to be computed with some care and we obtain The derivative of the two particle contribution takes a simpler form due to the general result (5.4). The two-particle contribution reads (θ 1 )T (θ 2 )e −mr cosh θ 1 −mr cosh θ 2 . e −mr cosh θ The two-particle contribution reads The first three quantities (6.26),(6.28) and (6.29) take the form that is similar to (5.5) and the same analysis as the unitary case applies. For the quantity B D,2 we need to integrate over two variables. Although technically more involved, physically it is the same as the unitary case, namely the exponential factors f exp (θ 1 , θ 2 ) = e −mr cosh θ 1 −mr cosh θ 2 and the defect matrix element f D (θ 1 , θ 2 , α) =T (θ 1 )T (θ 2 ) have certain peaks in a finite region and quickly damp away outside, as is shown in figure 6.7. The position of the peak is controlled by the parameter α. If the peaks of f exp and f D overlap, the defect has finite contribution to entanglement entropy, otherwise the contribution is negligible and we recover the bulk result. Similar to the unitary case, we give numeric estimation of the peak overlapping. From our final expression (6.31), there are two types of overlaps to consider, namely the overlap of f exp,1 = e −mr cosh θ , f D,1 =T (θ) − 1 (from the first two lines) and of f exp,2 = e −2mr cosh θ , f D,2 =T (θ) 2 − 1 (from the last two lines) whereT (θ) is given by (6.17). As in the unitary case, we set α = λ mr. For the two cases, we have I max,1 = e −mr and I max,2 = e −2mr . We set η * 0 = 1.0 × 10 −5 as before for both cases, the corresponding λ 0 are given in the following tables. For the first case, we have

Conclusions
In this paper, we considered the effect of integrable purely transmissive line defect on the bipartite entanglement entropy in 1+1 dimensional integrable field theories. We show that the topological defect does not affect the universal scaling of entanglement entropy in the UV limit. In the IR limit, the line defects have finite contribution to the leading corrections of entanglement entropy within certain range of the parameters that characterize the defect and have no contribution outside this finite range. The result holds for both unitary and non-unitary theories.
As a future direction, it is very interesting to investigate the finite volume/temperature entanglement entropy for integrable field theories, first for the bulk case and then for the boundary/defect case. Although obtaining the result for generic volume/temperature might be quite challenging, a systematic low temperature/large volume expansion for the bulk case should be within reach following similar ideas proposed in [51][52][53]. The generalization to the defect case is also known [54].
In the bulk case, the leading exponential correction of the entanglement entropy is universal even for non-integrable massive field theories [55]. In this paper, we have shown that the leading exponential correction is also universal for the topological defect case. It would be interesting to see if we could generalize this result to generic non-integrable massive field theories 6 .