Entanglement Content of Quantum Particle Excitations II. Disconnected Regions and Logarithmic Negativity

In this paper we study the increment of the entanglement entropy and of the (replica) logarithmic negativity in a zero-density excited state of a free massive bosonic theory, compared to the ground state. This extends the work of two previous publications by the same authors. We consider the case of two disconnected regions and find that the change in the entanglement entropy depends only on the combined size of the regions and is independent of their connectivity. We subsequently generalize this result to any number of disconnected regions. For the replica negativity we find that its increment is a polynomial with integer coefficients depending only on the sizes of the two regions. The logarithmic negativity turns out to have a more complicated functional structure than its replica version, typically involving roots of polynomials on the sizes of the regions. We obtain our results by two methods already employed in previous work: from a qubit picture and by computing four-point functions of branch point twist fields in finite volume. We test our results against numerical simulations on a harmonic chain and find excellent agreement.

In two recent works [32,33] we have investi- gated the entanglement entropy of gapped systems in zero-density excited states, where excitations have a (quasi-)particle interpretation. We found a very intuitive mathematical structure for the entropy differences between the ground state and such excited states. This structure allowed for simple closed formulae in particular limits. At the basis of our results laid the observation that the additional entanglement entropy due to such excitations is equal to the entanglement entropy of very simple states, involving a small number of qubits with amplitudes which represent the probabilities of finding particles in various system's regions. This idea of localized excitations and their one-to-one correspondence with increased entanglement has turned out to be surprisingly general. So far it has been shown to work for systems without interactions and in some instances of interacting models and particular excited states. We also found numerical evidence of the validity of our results for the harmonic chain and lattice (2 dimensions) and beyond the universal quantum field theory regime. In the present work we extend this intuition, with precise results for the free boson theory, to the entanglement entropy of disconnected regions and the (replica) logarithmic negativity. In one dimension, the configuration that we will be considering throughout this manuscript is illustrated by Fig. 1. Here L is the full length of the system and 1 , 3 are the lengths of subsystems A and B respectively. The length of subsystem C is then L´ 1´ 3 , which includes the mid-region of length 2 .
The von Neumann and Rényi entanglement entropies are measures of the amount of quantum entanglement, in a pure quantum state, between the degrees of freedom associated to two sets of independent observables whose union is complete on the Hilbert space. For a subsystem consisting of two disconnected regions as in Fig. 1, the Hilbert space may be seen as having the structure Consider a bipartition where the two sets of observables correspond to the local observables in the two finite-size complementary regions, A Ť B and C. Let the system be in a state |Ψy L . Then the von Neumann entropy associated to region A Ť B is where ρ A Ť B " Tr C p|Ψy LL xΨ|q is the reduced density matrix associated to subsystem A Ť B. One may obtain the von Neumann entropy S Ψ 1 p 1 , 3 , Lq as a limiting case of the sequence of nth Rényi entropies defined as thanks to the property lim nÑ1 S Ψ n p 1 , 3 , Lq " S Ψ 1 p 1 , 3 , Lq . (1.4) The logarithmic negativity of the same state |Ψy L [40,41,42,43,44,45] is defined as: (1.5) The geometry is the same as in Fig. 1 but the density matrix ρ T B A Ť B differs from ρ A Ť B , thanks to the operation T B , which represents partial transposition on subsystem B. Crucially, after such transposition, the matrix ρ T B A Ť B is no longer guaranteed to be positive-definite: specifically, it may have positive and negative eigenvalues. The operation Tr A Ť B |¨| involved in (1.5) represents the "trace norm" of ρ T B A Ť B , namely the sum of the absolute values of its eigenvalues. Similarly as for the von Neumann entropy, it is possible to define the logarithmic negativity as a limit of a more general set of quantities depending on an additional parameter n as: (1. 6) such that E Ψ p 1 , 3 , Lq :" lim nÑ1 E Ψ n p 1 , 3 , Lq, from n even . (1.7) The meaning of the last equation is that one sets n " 2m for m P N, and then analytically continues in m towards m " 1{2 (under appropriate conditions in the complex m plane making this analytic continuation unique). This definition was first proposed in [36,37]. From now on we will call these functions replica logarithmic negativities. Note that the (replica) logarithmic negativity is a measure of the entanglement between the non-complementary regions A and B whereas the von Neumann and Rényi entropies defined earlier are measures of the entanglement between the two complementary regions A Ť B and C. Therefore they are distinct measures capturing different information about the state.
From the definitions above it is clear that in order to compute these measures of entanglement we need to have access to the reduced density matrix ρ A Ť B and its partially transposed version, or at least to their eigenvalues. This may be achieved numerically for some systems, but if the system admits a quantum field theory description, then is also accessible analytically by relating the nth powers ρ n A Ť B and pρ T B A Ť B q n to specific correlators of local fields in quantum field theory [8,9,12]. In [12] these fields were termed branch point twist fields and characterized as symmetry fields in an n-copy version of the model under study. Branch point twist fields have since been systematically employed in the computation of measures of entanglement in massive QFTs and also in CFT. The entanglement measures considered in this paper, in the setup of Fig. 1, were studied systematically in the ground state of CFT in [34,35,36,37,46] and also partially in massive integrable QFT [47,48].
In the present work we want to extend our study to zero-density excited states. More precisely, we want to study the increment of the entanglement entropies of two disconnected regions and of the replica logarithmic negativities in a zero-density excited state, compared to the ground state. We define the new variables, , for i " 1, 2, 3, and r :" 1´r 1´r3 . (1.8) In this paper we will consider the scaling limit, defined as L, 1 , 2 , 3 Ñ 8 , with r, r 1 , r 2 , r 3 finite, and 0 ď r, r 1 , r 2 , r 3 ď 1 . (1.9) In terms of branch point twist field correlators, the increment of the Rényi entropies of two disconnected regions is given by: ∆S Ψ n prq :" S Ψ n pr 1 , r 3 q´S 0 n pr 1 , r 3 q (1.10) " lim LÑ8 1 1´n log « L xΨ|T p0qT pr 1 LqT ppr 1`r2 qLqT ppr 1`r2`r3 qLq|Ψy L L x0|T p0qT pr 1 LqT ppr 1`r2 qLqT ppr 1`r2`r3 qLq|0y L ff , where S 0 n pr 1 , r 3 q are the entropies in the ground state. As the notation suggests, we will see later that ∆S Ψ n prq is a function of r only. The increments of the replica logarithmic negativities are given by: ∆E Ψ n pr 1 , r 3 q :" E Ψ n pr 1 , r 3 q´E 0 n pr 1 , r 3 q (1.11) " lim LÑ8 log « L xΨ|T p0qT pr 1 LqT ppr 1`r2 qLqT ppr 1`r2`r3 qLq|Ψy L L x0|T p0qT pr 1 LqT ppr 1`r2 qLqT ppr 1`r2`r3 qLq|0y L ff , and are functions of r 1 and r 3 only. As standard, T is the branch point twist field,T is its hermitian conjugate and |0y L is the ground state in the compactified space of Fig. 1. As is wellknown [12] the branch point twist fields are symmetry fields associated with cyclic permutation symmetry in an n-copy version of the theory under consideration. From this definition it follows that T "T if n " 2, and we will see later that, indeed, the results for the replica logarithmic negativity equal those for the Rényi entropy (up to a sign) for n " 2. When employing the branch point twist field technique, we will therefore be computing ratios of four-point functions in the infinite volume limit. Exact explicit formulae for the four-point functions above are generally inaccessible even in the ground state and/or CFT. Remarkably though, the increments computed here, in this particular scaling limit, turn out to be the logarithms of order kˆn polynomials (where k is the number of excitations) in the parameters r 1 , r 3 (or just in r for the entanglement entropy) with positive, integer coefficients. This paper is organized as follows: In section 2 we summarize the main results of the paper. In section 3 we show how these results may be obtained from a qubit interpretation. That is, we show that the formulae for the increment of entropies (negativities) in excited states of QFT equal the entropies (negativities) of a very special class of states involving a finite number of qubits. In section 4 we obtain the same results from a branch point twist field computation in the massive free boson theory. In section 5 we present some numerical results for the harmonic chain. We conclude in section 6. In Appendix A we summarize some formulae for the Rényi entropies and replica negativities for particular values of n and k. In Appendix B we derive the logarithmic negativities of states of two and three identical excitations from the exact diagonalization of the reduced and partially transposed density matrix. In Appendix C we extend our qubit results for the Rényi entropies of two disconnected regions to the case of any number of disconnected regions. In Appendix D we summarize some useful results for the residua of certain types of first, second and third order poles which play a role in our branch point twist field computation.

Summary of the Main Results
Consider a state consisting of a single particle excitation. Let us denote the entropy increments of such a state by ∆S 1 n prq with r as defined earlier (1.8). We find that This is identical to the result found in [32,33] for a connected region of length r and indeed, further results for the von Neumann and other types of entropy, for excited states consisting of a higher number of excitations, identical or not, also take the same form as found in [32,33]. We will not report all of these here again. We report however the formula for a k-particle excited state consisting of identical excitations, as it will feature later on. This is given by This means that the entanglement entropy increase due to the presence of a discrete number of excitations depends only upon the size of the regions A and B relative to that of the whole system, and not on their connectivity. Indeed, in Appendix C we show that the same conclusion can be drawn for a subsystem consisting of any number of disconnected regions.
Results are more interesting for the (replica) logarithmic negativity. First, we find that it is a function of r 1 , r 3 , where each parameter now enters independently. For a state consisting of a single particle excitation the increment of the replica logarithmic negativity is given by the simple expression: where r.s denotes the integer part. The increment of the logarithmic negativity can be obtained by analytically continuing this expression from n even to n " 1. Evaluation of the sum above 3) for n " 2, 4, 6 and 8 (the higher the value of n, the more negative the value of the function). Right: The function (2.5). As we can see, there is a change of curvature and sign when taking the limit n Ñ 1. A similar interesting analytic behaviour near n " 1 was found in [37] for the compactified free boson. Recall that 0 ď r 1`r3 ď 1, which restricts the domain of definition of the functions shown. A cross-section of these surfaces is provided in Fig.4 (left).
For a state consisting of k distinct excitations (particles with distinct momenta), the result is simply k times the above, just as for the Rényi entropies. The case of identical excitations (identical momenta) is more interesting. Consider an excited state of k identical excitations. The increment of the replica logarithmic negativity is given by: where the coefficients A p,q are defined as follows: and σ n 0 pqq represents the set of integer partitions of q into n non-negative parts. Note that the coefficients are zero whenever any of the arguments of the factorials in the denominator becomes negative and this selects out the partitions that contribute to each coefficient for given values of p and q. As should be, formula (2.3) is the k " 1 case of (2.6). Indeed the coefficients inside the sum (2.3) are nothing but the number of partitions of p into n parts, p of which are 1 and n´p of which are 0, with the constraint that there are no consecutive 1s.
Since the coefficients A p,q are rather non-trivial, it is not easy to perform the sums in (2.6) explicitly and the analytic continuation leading to the logarithmic negativity is rather involved. We have not been able to obtain closed formulae for the logarithmic negativity for all k but have obtained expressions for k " 1, 2, 3. For k " 2 it is given by and ∆ 0 " r 4`1 0r 2 r 1 r 3`7 r 2 1 r 2 3 , ∆ 1 "´2r 6´3 0r 4 r 1 r 3´6 9r 2 r 2 1 r 2 3`2 0r 3 1 r 3 3 . (2.10) The case k " 3 has been studied in Appendix B and the logarithmic negativity in this case is given by (B.34). The explicit polynomials (2.6) for a few values of n and k are listed in the tables of Appendix A. The function (2.6) for n " 2, 4, 6 and 8 and k " 2. Right: The function (2.8). Again, there is a change of curvature and sign when taking the limit n Ñ 1. Recall that 0 ď r 1`r3 ď 1, which restricts the domain of definition of the functions shown. A cross section of these surfaces is provided in Fig.4 (right).
For mixed excited states consisting of k i particles of identical momentum p i with i " 1, . . . m and p i ‰ p j for i ‰ j we have that the Rényi entropies and the (replica) logarithmic negativities can be expressed in terms of the building blocks given above, namely: Note that the functions (2.5), (2.8) give nice simple illustrations of the non-trivial nature of the analytic continuation from n even to n " 1. Although for integer n the functions (2.3) and (2.6) are simple polynomials in r 1 , r 3 , the logarithmic negativity involves roots which appear to be of increasing order as k is increased. As explained in Appendix B, we have obtained the analytic continuation of (2.6), for k " 2 giving (2.8) and for k " 3 giving (B.34), by explicitly computing the eigenvalues of the partially transposed reduced density matrix obtained from a qubit interpretation.
The marked difference between functions (2.5), (2.8) and (2.3), (2.6) can also be seen in Figs. 2 and 3 where they are shown to exhibit opposite sign and curvature. The gradual change in curvature as n is reduced is particularly clear from the cross-section views given in Fig. 4.
The expression (2.6) with (2.7) suggests a combinatorial interpretation for our results and indeed such interpretation can be made concrete. A one-to-one relationship between (2.6) and the generating function of the number of simply-connected graphs with 2nk nodes, nk edges and certain connectivity conditions will be proven in [38].
Let us for now present two approaches whereby these results can be derived: from an interpretation in terms of qubits, and from a direct calculation within massive QFT.

Qubit Interpretation
It was shown in [32,33] that exact expressions for the increment of entanglement entropy in a zero-density excited state and a bipartition involving only simply-connected domains could be derived from a so-called "qubit interpretation". This interpretation was first suggested by the result (2.1) and the realization that it is identical to the Rényi entropy of a two-qubit state of the form |Ψ qb y " ? r|01y`?1´r|10y: the entanglement entropy increase depends only on whether the particle is localized in region A Ť B with probability r, or on region C, with probability 1´r. The strength of this picture lays on the fact that it can be extended naturally to states consisting of multiple excitations. Their increase of entanglement entropy is also the entropy of multi-qubit states, with coefficients which also have a simple probabilistic interpretation. For a single connected region, this interpretation has been numerically shown to hold even in the two-dimensional harmonic lattice [33].
This interpretation worked so well for the entanglement entropy of simply-connected domains and appears so natural that one may expect it to still hold for disconnected regions and for the logarithmic negativity. Clearly, it suggests that results should be independent of the topology of the regions. This is in agreement with our first main formula (2.1) (and its multiple-particle generalisation) and holds as well in the case of a region composed of any number of disconnected parts (see Appendix C). As we see below, this interpretation also allows us to obtain the result (2.6) for the the replica logarithmic negativities.

Single Particle Excitation
Let us consider the simplest case, that of an excited state of a single excitation. We can simply write the following three qubit state: where the first qubit represents the presence (1) or not (0) of a particle in region A, the second qubit represents the same for region C, and the final qubit likewise for region B. Tracing over the mid-qubit we have that and In matrix form we have ?
where the first row and first column refer to the states involved in (3.2)-(3.3). The eigenvalues of ρ A Ť B are: λ 1 " 0 , λ 2 " 0 , λ 3 " 1´r , λ 4 " r , (3.5) and Note that in the latter case, the last eigenvalue is clearly negative. This means that, for a one-particle excitation we have: which is what we expected, and E 1 n pr 1 , r 3 q " logˆr n 1`r n 3`p r`?r 2`4 r 1 r 3 q n 2 n`p r´?r 2`4 r 1 r 3 q n 2 n˙. (3.8) Interestingly, for n integer, even or odd, there is no square-root dependence of the polynomial above (the square-roots always cancel). Indeed, it can be equivalently written as given in (2.3). However, the logarithmic negativity itself involves the square root in (3.8). From the eigenvalues above this gives which is the result (2.5) reported earlier. As we can see, (3.7) and (3.8) are in general rather different functions. However, it is easy to show that the polynomials inside the logarithm coincide for n " 2.

Multiple Particle Excitations
The approach developed above, where we compute the eigenvalues of the reduced (partially transposed) density matrix, has been applied also to the cases of two and three identical excitations, leading to the results (2.8) and (B.34). The explicit computation is reported in Appendix B. For higher number of excitations, this exact diagonalization approach becomes quickly intractable. Here, we instead present an approach based on evaluating the replica trace more generally. It directly leads to the results (2.1) and (2.6). However, in the latter case, we do not know in general how to extract the analytic continuation to n " 1 from even values, necessary for the logarithmic negativity.

Distinct Excitations
Having considered the particular case above, it is worth noting that the qubit states associated with multiple distinct excitations are by construction factorized states. For instance for two distinct excitations we have the state: qb y " r 1 |1100000y`r|001100y`r 3 |000011y`?r 1 r 3 p|100001y`|010010yq ? r 1 r p|100100y`|011000yq`?rr 3 p|000110y`|001000yq where the first two qubits live in region A, the next two in region C and the final two in region B and |Ψ p1q qb y is the state (3.1). Thus the trace is also factorized and can be written as From this factorization it is easy to see that both the Rényi entropies and the logarithmic negativity will be twice their value for a single excitation. This generalizes to k-particle states consisting of distinct excitations.

Identical Excitations
Consider instead a k-particle state consisting of k identical excitations. Its associated qubit state can be written as: where σ 3 0 pkq represents the set of integer partitions of k into three non-negative parts, and we recall once more the definitions (1.8).
The state may be interpreted as follows: each state |k 1 k 2 k 3 y is a state of k 1 excitations in region A, k 2 excitations in region C and k 3 excitations in region B. The square of the corresponding coefficient c 2 k 1 ,k 2 ,k 3 is the associated probability that this configuration occurs if we were to place randomly and independently, with uniform distribution, k particles on the interval r0, 1s covered by three non-intersecting subintervals of lengths r 1 , r and r 3 .
Provided that all vectors |k 1 k 2 k 3 y are normalized to one, then the vector |Ψ pkq qb y is also a unit vector, xΨ From this expression it is then possible to explicitly construct the reduced density matrix and its partially transposed version as: Here the sums run over all non-negative integers, and whenever the constraint k 1`k2`k3 " k is violated, the corresponding coefficient c k 1 ,k 2 ,k 3 is zero by definition. From this constraint we know in fact that 0 ď k 2 ď k so we could have restricted the summation range k 2 P N 0 much more. However we will write it like this for now, for simplicity, and discuss the summation ranges more precisely at the end of our calculation. With these results we can now evaluate the matrix elements of the n-th powers of the reduced density matrices above. These are given by: This means that the factor where we defined q :" ř n j"1 k j 3 . This finally allows us to rewrite (3.22) as and σ n 0 pqq represents the set of integer partitions of q into n non-negative parts. We have relabelled k j 3 :" k j , and the range of the sums in p and q is determined by the condition of Z p,q ‰ 0. From the definition (3.24), 0 ď p ď k. Regarding the values of q, we know that q can not be negative (by definition) so q ě 0. Its maximum value is obtained if k´p´k j " 0 for all j. This corresponds to q " npk´pq. In fact the sum over q can be rewritten as with r 1`r3 " 1´r so that S k n pr 1 , r 3 q ": S k n prq " as reported in (2.2). Therefore, the entanglement entropy depends only on the parameter r and is given by exactly the same expression as found in [32,33]. In other words, in the qubit picture, the entanglement entropy depends only on the overall size of regions and not on whether or not they are connected. This implies that the same result should also hold for more than two disconnected regions. Indeed this can be shown by similar methods. The proof is presented in Appendix C.

Results for Replica Logarithmic Negativity
A similar analysis can be carried out for the replica logarithmic negativity. Starting with (3.22) the second delta function gives the condition: and this eliminates the sums over k j 1 with j " 1, . . . , n. We then have sums over k j 2 and k j 3 left but we can also eliminate one of these by implementing the second set of delta-functions together with the conditions above. This gives the constraints, We may regard this equation as a first order difference equation for the sequence k j 2 . The solution to such an equation is the sum of the solution to its homogenous version (a constant) and a particular solution of the full equation which can be worked out by inspection to be´k j 3´k j`1

. The general solution is then
where γ is an arbitrary constant. With this we also have that k j 1 " k´γ`k j 3 . We can now evaluate the product where q :" ř n j"1 k j 3 and p " k´γ. Relabelling k j`1

3
:" k j we then find as given in (2.6) and (2.7). The range of sums in p and q is fixed by selecting out those contributions for which A p,q ‰ 0. This requires that the arguments of the factorials in the denominator remain non-negative, which in turn restricts the type of partitions that can contribute to the sum over k 1 , . . . , k n . Consider the sum in p. The range of this sum can be determined easily from the relation (3.32). This implies that 0 ď γ ď 2k. Together with the definition of p this gives´k ď p ď k. This guarantees that all arguments of the factorials in the denominator remain non-negative.
The range of values of q can also be determined as follows. The lower limit is easy to establish as whenever p ă 0 the partitions contributing to A p,q must have k j ě´p. Thus, the smallest value of q giving a non-vanishing contribution corresponds to taking all k j "´p for all j which gives q "´np. On the other hand, if p ě 0 then the smallest value q can take is zero corresponding to all k j " 0. This fixes the lower bound to maxp0,´npq. Let us now consider the upper bound. Given a certain p, the largest value q can take corresponds to having k´p´k j´kj`1 " 0 for all j, or k j`kj`1 " k´p. Writing we have obviously that the right hand side gives npk´pq whereas the left hand side gives 2 ř n j"1 k j " 2q. Therefore, for generic parity of n and k´p, we obtain q " r npk´pq

Computation from Branch Point Twist Fields
In this section we evaluate the leading volume contribution to the increment of Rényi entropies of two disconnected regions and the increment of the replica logarithmic negativity using branch point twist field methods in the free massive boson theory. This provides a formal derivation, from QFT methods, of the results presented in the previous sections, which were argued for based on the qubit picture. We follow the same strategy as for the computation presented in [33], generalized to the study of the four-point functions entering the definitions (1.10) and (1.11). All the techniques that we use in this section (branch point twist fields, finite volume form factors, doubling trick) have been exhaustively reviewed in [33] and we refer the reader to this paper for further details. Here we will just enumerate the main ideas and techniques needed.

Techniques and Main Formulae
As usual in the context of branch point twist fields we will be working with a replica free boson theory, consisting of n copies of the original model. We then first implement the doubling trick [49] on the replica free boson theory to construct a replica free massive complex boson theory. The doubling introduces a new U p1q internal symmetry on each replica. This makes it possible to diagonalize the action of the branch point twist fields on asymptotic states. Associated to these new U p1q symmetries there are n twist fields T p with p " 1, . . . , n in terms of which the branch point twist field can be expressed as: where p labels the U p1q sectors. The form factors of the U p1q twist fields will constitute the building blocks of all our computations so we will recall the main formulae here. They can be computed by the standard techniques [50,51,52] and have been known for a long time [53,54]. The formulae that we will need here are: where θ 12 :" θ 1´θ2 . This is the normalized (by the vacuum expectation value) matrix element of a sector-p U p1q twist field between the vacuum in that sector and a two-particle state. The state consists of a complex free boson (with creation operator a : p pθ 1 q and rapidity θ 1 ) and its chargeconjugate boson (with creation operator b : p pθ 2 q and rapidity θ 2 ). The dependence on the sector is determined by the form factor equations, in particular the crossing relation f n p pθ`2πiq " e 2πip n f n p pθq. Since the theory is free, higher particle form factors can be obtained by simply employing Wick's theorem. For the complex free boson they have the structure F p,n pθ 1 , . . . , θ m ; β 1 , . . . , β m q :" p x0|T p p0qa : p pθ 1 q¨¨¨a : p pθ m qb : where σ are all elements of the permutation group S m of m symbols and each two-particle form factor can be associated with a particle-antiparticle contraction. Note that any other form factors are zero (e.g. if the numbers of particles and antiparticles are different). The second important ingredient in our calculation is understanding the structure of the Hilbert space of the theory in the replica doubled theory. That is, the description of the asymptotic states in each sector of the theory in a basis upon which the U p1q twist fields act diagonally. The solution to this problem was discussed at length in [33] and can be schematically presented as follows. First, that is, the creation operator of a real boson on copy j is mapped after doubling into a linear combination of creation operators on copy j associated with a complex free boson (+) and its charge conjugate (-). As expected, the branch point twist field action on these operators is such that they are mapped to operators on copy j`1 of the theory. However, there exists a new basis ta : p pθq, b : p pθqu with p " 1, . . . , n on which the branch point twist field acts diagonally, with each U p1q field (4.1) acting on a particular sector. For example, in can be shown where the operators a : p pθq, b : p pθq satisfy also a complex boson algebra. We will from now on perform computations on this basis, where we can employ the matrix elements of U p1q twist fields seen above. Although this simplifies the computation of form factors enormously, the price to pay is a more complicated structure for the excited states which must now be also expressed in this new basis (see e.g. (4.16) for an example).
The third and final issue to consider is the finite volume extension of the above. In particular, rapidities will be quantized in finite volume and the quantization conditions will depend upon the sector p the corresponding creation/anhilation operator is acting on, and upon an index α that parametrizes the periodicity conditions for the fields Φ j px`Lq " Φ j`α pxq, where Φ j pxq is the complex bosonic field with associated creation operators in (4.4). A set of (generalized) Bethe-Yang equations [55,56] can be written as where "˘is the U p1q charge of the particle and J P Z. The sum over a complete set of states in sector p with quantization condition α can be written as The summation runs over the integer sets tJ˘u " J1 , . . . , Jm`, J1 , . . . , Jḿ´( and the rapidities satisfy the quantization conditions (4.6) as: We can then define the complete sum over all sectors The branch point twist fields intertwine between the different quantization sectors. Denoting the corresponding Hilbert space by H α we can write The excited states |Ψy L we consider in this calculation are elements of the trivial quantization sector H 0 " H n . This, combined with the properties of the branch point twist fields means that the four point functions of interest may be spanned as ) where x 1,2,3 denote the positions of the branch point twist fields which is related to the original lengths in Fig. 1 as (4.13) Employing (4.1) the four-point functions then factorize into products of contributions from every sector, involving the finite volume matrix elements of U p1q twist fields. These matrix elements are related to infinite volume form factors (4.3) in a simple well-known fashion [57,58] which involves dividing each infinite volume form factor by the square root of the product of the associated particle energies (in general, for interacting theories, this would be the density of states) of the excitations involved, and employing the crossing property of form factors. This can be summarized as F p,n pθ 1 1 , . . . , θ 1 s 1 , β 1`i π, . . . , β q`i π; β 1 1 , . . . , β 1 up to exponentially decaying corrections Ope´µ L q and Epθq " m cosh θ. An important feature of the two-particle form factors (4.2) is their behaviour near kinematic poles

Four Point Functions in Single-Particle Excited States
Let us focus on the calculation for a single-particle excited state denoted by |1y L . If the excitation has rapidity θ, the state in the U p1q basis has the form where C n ptN˘uq are the (known) coefficients containing all the phase factors of the basis transformations (see subsection 4.1.3 in [33] for an explicit example). The rapidity is the solution of the quantization condition (4.6) h p, 0 pθq " mL sinh pθq " 2πI 0 . (4.17) The summation runs over the integer sets tN˘u " N1 , N1 , . . . , Nǹ , Nń ( subject to the condition n ÿ p"1 ÿ "˘N p " n . (4.18) Following previous considerations, the four point functions are with the different sector contributions F p´Np ,Ñp¯" p;L x0| ra p pθqs Np rb p pθqs Nṕ T p p0q I p 1T p px 1 q I p 0 (4.21) In sector n the U p1q twist fields coincide with the identity operator, hence the contributions from this sector are just the normalization of the state F n´Nn ,Ñn¯"F n´Nn ,Ñn¯" Nǹ !Nń !δ Nǹ ,Ñǹ δ Nń ,Ñń . (4.23) The strategy of the calculation from this point on is to evaluate the matrix elements of the twist fields using the finite volume form factor formula (4.14), turn the summation for quantum numbers into contour integrals, and by contour manipulation extract the leading volume contribution. First, we focus on the increment of the Rényi entropy of two disconnected regions. We will later see that the logarithmic negativity can be computed in a very similar way.

Entanglement Entropy of two Disconnected Regions
Let us compute the ratio of correlators (1.10) for a single-particle excited state. We need to evaluate (4.19) with (4.21). We will start by focussing on the evaluation of (4.21) through the introduction of the complete sums I p 0 and I p 1 . We first introduce some notations. We denote the rapidities of the complete sets of states, starting from left to right, by β i , α i ,β i and their numbers by m , k ,m , respectively. As usual "˘. The rapidities satisfy the following Bethe-Yang quantization conditions (4.6) We introduce also the following notation for the various sets of rapidities tθu N " tθ, . . . , θu loooomoooon N , tβ u m " tβ 1 , β 2 , . . . , β m u , tβ u m " tβ 1`i π, β 2`i π, . . . , β m`i πu , (4.27) and similarly for the α i andβ i rapidities. Using the finite volume form factor formula (4.14) and these notations, (4.21) takes the form where P pθq " m sinh θ are the particle momenta. We have also used the symmetry property of the form factors (since we are dealing with free bosons, the form factors are symmetric in all rapidities). We stress that at this stage all the rapidities in the formula are solutions of the appropriate Bethe-Yang equations. This will not be the case anymore when we rewrite our formulae in terms of contour integrals. We will start by transforming the tβ u set by writing where we now consider the function h p, 1 pβ i q as the function of β i defined in (4.24) so as to ensure that the integrand has a pole exactly when the Bethe-Yang equation is satisfied. The function f on the left hand side of the equation is to be taken from the J i summand in (4.28), and the contour is the sum of small contours around the solutions of the Bethe-Yang quantization condition (4.6) h p, 1 pβ i q " 2πJ i , (note that these all lie along the real line) We deform the contour into a contour encircling the real axis with positive orientation C Ô and subtract the residua of the form factors' kinematic poles. The form factors (4.3) have kinematic poles, when either β i " θ in the first form factor or β i " α j for some j, in the second form factor. However, if the quantum number I j " I 0 , the quantum number of the excited state, then α j " θ and a second order pole can result from the combination of the kinematic singularities of the two form factors. This leads to a different residue, hence we should proceed with care. In order to correctly account for the different contributions, we separate the quantum number sums into a group that coincide with I 0 and the rest where we used the symmetry property of the form factors and the invariance of the integral under relabelling of variables tα u k . With these considerations, the contour after the deformation is where α j ‰ θ, and the multi-contour for all β i rapidities is (4.33) We can expand the integration multi-contour by substituting (4.32) into (4.33) and employing the same symmetries of the integrand mentioned above. As a result there are several terms of the multi-contour that give the same residue and we can write where G pm Ô , m θ , m α q is a combinatorial factor arising from the power expansion, m θ is the number of contours about the rapidities θ, and m α is the number of rapidities in the set tγ u m α which is the subset of tα u k , itself being involved in a particular term of the expansion (4.31). The nonnegative integer summation indices are constrained by m " m Ô`m θ`m α , (4.35) Note that although, in general, the order in which the integrals over the various contours are performed matters, in the expansion (4.34) we can obviate this by employing the fact that all such orderings are equivalent under relabeling of rapidities and that in this case all such relabelings are equivalent due to the symmetries of the free boson form factors. Let us focus on the residua of the contour C θ (that is the poles with β i " θ). As mentioned before, if M ‰ 0, the product of the two first form factors produces both second and first oder poles. We denote by m θ,1{2 the number of first order poles arising from the first/second form factor, and by m θ,d the number of second order poles. They satisfy the condition m θ " m θ,1`m θ,2`m θ,d . (4.36) Evaluating the residua of these poles produces several types of factors. In the phase factors, the rapidity is trivially replaced by θ or γ i (recall that the rapidities tγ u m α are a subset of tα u k as defined in (4.34)). There are constant factors coming from the combination of the form factor residua and the denominators, simplified by using the Bethe-Yang equations (4.24) and (4.25). These factors contain derivatives in the case of higher order poles. Moreover, in case of first order poles, the contracted rapidity θ or γ j appears in the non-singular form factor. In the case of C θ , this also means taking the regular part of that form factor, defined as the part that does not contain any two-particle form factors depending on β i´θ . The singular terms thus omitted are already accounted for in the double order pole residua. The form factors minus such omitted terms are what we later denote by an additional subscript "reg" (see Eq. (4.43)). There are combinatorial factors coming from counting all possible contractions of β i with θ or γ j rapidities inside a given form factor, that we discuss in detail later. The residua of the singular two-particle form factor(s) are given by R θ,1 " 1 , R θ,2 " e 2πi p n , R α " e 2πi p n , R θ,d " LE pθq g n p pr 1 q , (4.37) where R θ,1 , R θ,2 and R α correspond to first order poles and R θ,d corresponds to second order poles. r 1 " x 1 L " 1 L and g n p prq " 1´r`re 2πi p n , (4.38) the function that already played an important role in the computation of the entanglement entropies of one interval [33]. The same consideration is valid for theβ rapidities.m θ,1{2 denotes the number of first order poles resulting from the third/fourth form factor having a singularity, andm θ,d denotes the number of second order poles. The residua of the singular two-particle form factor(s) are given byR θ,1 " 1 ,R θ,2 " e 2πi p n , where r 3 " x 3´x2 L " 3 L . We stress that, as earlier, only the second order poles contribute a factor proportional to the volume LE pθq.
Successive computation of the residua associated with rapidities β i about any of the rapidities tγ u m α andβ i about rapidities tγ um α results in a sum of contributions where a "reshufling" of the rapidities tα u k has taken place. These will now appear in the first or second and in the third or fourth form factors, in various combinations. To reflect this, we introduce the following partition of the rapidities where tα ij u k ij denotes the subset of rapidities in tα u k that feature as arguments of the ith and jth form factors after the residue evaluation. k ij are the cardinalities of these sets satisfying k " k 13`k 14`k 23`k 24 , m α " k 13`k 14 ,m α " k 14`k 24 . With this partitioning, a generic sum over k turns into where also tI u k " tI 13 u k 13 Ť tI 14 u k 14 Ť tI 23 u k 23 Ť tI 24 u k 24 . The new form of (4.10) after all these manipulations is given by Aptα u k , tβ u m Ô , tβ um Ô q Bptα 13 u k 13 , tα 14 u k 14 , tα 23 u k 23 , tα 24 u k 24 q Dpθ, m θ , M , N p ,Ñ p q Rpm θ,1 , m θ,2 , m θ,d , m α ,m θ,1 ,m θ,2 ,m θ,d ,m α q with Aptα u k , tβ u m Ô , tβ um Ô q :"  Dpθ, m θ , M , N p ,Ñ p q :" e ix 1 m θ P pθq e ipx 2´x1 qM P pθq e ipx 3´x2 qm θ P pθq e´i x 3Ñ p P pθq rLE pθqs and, finally where the "reg" subscript of the form factors indicates they do not contain any contractions within the sets of θ rapidities (e.g. as explained earlier).G in (4.47) contains the combinatorial factors from the summation over contours and from the combinatorics of the residua. The prime over the summation of m θ,i ,m θ,i for i " 1, 2, d indicates certain constraints due to the limited total number of θ rapidities that can take part in the residue calculation, namely The next step of the calculation is to transform the quantum number sums for tI u into contour integrals as well where the contour is a combination of small contours about the Bethe-Yang solutions (4.25). As before, we deform the contour into one encircling the real axisC Ô , however, we need to subtract the residua of poles at α i " θ, since these poles are not included in the original contour. That is It is important to note thatC Ô must be chosen such as to run closer to the real axis than the contour C Ô , for the tβ u m Ô and tβ um Ô rapidities, to avoid capturing undesired residua.
The denominator e ih p, 0 pα i q´1 of (4.49) is singular at α i " θ and the form factors (which would be part of the function f pα i , . . .q in the numerator) can also have kinematic singularities at this point. Since any given α ij rapidity appears in two form factors it follows then that we can have first, second, and third order poles. Using the residue formulae in Appendix D, it is clear that only third order poles can generate contributions that are proportional to mL. This is the reason why we will only need to consider third order poles in order to obtain the largevolume leading contribution to the entanglement entropies. To find the leading contribution, let us denote by s ij the number of third order poles evaluated by rapidities belonging to the set tα ij u k ij . These numbers are constrained by the number of α ij and θ rapidities. Clearly s ij ď k ij . In addition, we have the following less trivial constraints The maximum power of the volume is then ∆ " 0. This corresponds to the situation when all dependency on the θ rapidities has been cancelled by the evaluation of residua. Hence there are no other poles to consider other than third order ones (as these will produce subleading contributions). Using the results of Appendix D, the residua from third order poles are (up to factors depending on rapidities other than θ) and r 2 " x 2´x1 L " 2 L . The origin of the function (4.58) is explained in Appendix D. After relabelling k ij´s ij Ñ k ij and using (4.55), the combinatorial factor arising from the evaluation of the residua, that is from counting the number of ways of picking the rapidities for the residua and the different contractions, has the simple form with Ω n p pN, r 1 , r 2 , r 3 q "g pr 1`r2 q g pr 1`r2`r3 q m minpm,mq ÿ s"maxpM´N`m`m,0q rg pr 1`r2`r3 qg pr 2 qs s rg pr 1`r2 qg pr 2`r3 qs´s s! pm´sq! pm´sq! pN´M´m´m`sq! . (4.61) Here the number of sums has been greatly reduced thanks to the constraints (4.55). The four remaining sums correspond to the four independent variables in (4.56) which can been relabelled for convenience (s 13 Þ Ñ s, M Þ Ñ M, m θ,1 Þ Ñ m,m θ,1 Þ Ñm and N p Þ Ñ N ). It is easy to show that the sums in m,m and s are nothing but the multinomial expansion of 1 pN´M q! r1´A´B`ABCs N´M with A "g pr 2`r3 q g pr 1`r2`r3 q , B "g pr 1`r2 q g pr 1`r2`r3 q , C "g pr 1`r2`r3 qg pr 2 q g pr 1`r2 qg pr 2`r3 q , (4.62) The last three sums are then simply and Ω n p pN, r 1 , r 2 , r 3 q " " g n p pr 1 q g n p pr 3 q´"1´e 2πi p n ı 2 r 1 r 3 Note that, from (4.63) we see that the dependence on the parameter r 2 drops out. Thus the results do not depend on the distance between regions A and B as long as 0 ď r 2 ď 1 in the scaling limit. This leads to the following final result for the four-point function Recall that this is exactly the same result as for the single region entanglement entropy [33] with region length r 1`r3 . Therefore, as observed earlier for the qubit states, the increment of the entanglement entropy due to the presence of one excitation is independent of the connectivity of the region under consideration. This can be shown more generally both with qubits (see Appendix C) and using form factors.

Replica Logarithmic Negativity
To calculate the leading order large-volume contribution to (4.22) we have to go through exactly same steps as presented in the previous section. For the logarithmic negativity the last two form factors are "exchanged" hence the quantization condition in the complete set of states inserted between the last two fields is now different. This only affects the values of some residua. The changed residua areR θ,2 " e´2 πi p n , The final result for (4.22) is lim LÑ8Fp`Np , Npp ;L x0| T p p0qT p px 1 qT p px 2 q T p px 3 q |0y p;L "˜ź "˘N p !Ω n p`N p , r 1 , r 2 , r 3˘¸`O`L´1˘, (4.67) whereΩ n p pN, r 1 , r 2 , r 3 q " " g n p pr 1 q g ń p pr 3 q´"1´e 2πi p n ı " 1´e´2 πi p n ı r 1 r 3 with the functionĝ n p pr 1 , r 3 q " 1´r 1´r3`r1 e 2πi p n`r 3 e´2 πi p n . (4.69) As for the Rényi entropies we find that the functionΩ n p pN, r 1 , r 2 , r 3 q does not ultimately depend on the parameter r 2 . The final result for the logarithmic negativity for the one-particle state becomes

Multi-Particle States
Following the notations introduced in [33], a general multi-particle state in the U p1q basis has the form |k 1 , . . . , k m y L "¨m where k q is the multiplicity of rapidity θ q in the multi-particle state. Each rapidity θ q with q " 1, . . . , m is a solution of the Bethe-Yang equation (4.25) with quantum number I 0,q . All such quantum numbers are distinct for different values of q. D kq n ptN q,˘u q are the numerical coefficients involved in expressing a k q -particle state with coinciding rapidities in the U p1q basis (they are related to the C n coefficients in (4.65)-(4.70) and shown in Eq. (4.43) of [33]). The integer sums obey the selection rules ÿ "˘n ÿ p"1 N q, p " nk q , @q . The four-point functions factorize into the product of U p1q sector contributions as before in (4.19) and (4.20), but now they depend on all the rapidities F p´t Np u, tÑp u¯" p;L x0|˜m ź q"1 ra p pθ q qs N q,p rb p pθ q qs N q,ṕ¸T ‰Ñ q,ṕ¸| 0y p;L , (4.73) andF p´t Np u, tÑp u¯" p;L x0|˜m ź q"1 ra p pθ q qs N q,p rb p pθ q qs N q,ṕ¸T are the increments for a k-particle state of identical particles and the integer sums still obey (4.72). Recall that r 1`r3 " 1´r and (4.75) is invariant under r Þ Ñ 1´r. We cross-checked these final results against those obtained from the qubit picture for several multi-particle states and found perfect agreement.

Numerical Results
In this section we present some numerical results for the Rényi entropies and replica logarithmic negativities of a harmonic chain. Recall that this is a discrete theory whose continuum limit is the massive free boson. Therefore we expect that for appropriate parameter choices we should find agreement with our predictions of the previous sections. The dispersion relation is where ∆x is the lattice spacing, m is the mass and p is the momentum of the excitation. The QFT regime corresponds to ∆x¨p ! 1 where we recover the relativistic dispersion relation. The numerical procedure was explained in detail in Appendix A of [33] and is based on exact diagonalization. In all the plots we have chosen: ∆x " 0.02, m " 1, and L " 20 . Numerical values of the 2nd Rényi entropy for r 1 " r 3 " R 2 as a function of r " 1´R and a single excitation. These are given for four different choices of the momentum of the excitation: p " 0 (circles), p " π (triangles) and p " 2π (squares) and 4π (rombi). Right: Numerical values of the 2nd (circles), 3rd (triangles) and 4th (squares) Rényi entropies for p " 4π (single excitation) compared to analytic formula (continuous curves).

Rényi Entropies
We have evaluated numerically the Rényi entropies for a subsystem consisting of two disconnected regions of the same length r 1 " r 3 . The subsystem is A Ť B with A " r0, RL{2s and B " rL{2, L{2`RL{2s, and the values of R are taken from 0 to 1 by increments of 0.1. The total subsystem's length is therefore 1` 3 " RL, ranging from 0 to L. Numerical values of the 2nd Rényi entropy of a state of two identical excitations (circles) and two distinct excitation (triangles). The momenta are p 1 " p 2 " 4π for identical excitations (k " 2) and p 1 " 3π, p 2 " 4π for two distinct excitations. In all cases agreement with the analytic results (continuous curves) is excellent.
In Fig. 5 we present results for the Rényi entropies of one excitation. The figure on the left explores the dependence on the momentum of the excitation; this dependence is shown to be very small, and the numerical results to agree very well with the analytical prediction. At p " 0 the disagreement is slightly larger because in this case only the correlation length m´1 " 1 controls the accuracy. The relative disagreement is however higher for the small value R " 0.1 (omitting the trivial value R " 0) than it is at R " 0.5, where correlation length effects should be smallest. The entanglement is higher than predicted at the small value R " 0.1, where the regions A and B are small, presumably because one observes additional entanglement due to the correlations within distances comparable to the correlation length; while it is smaller at R " 0.5, presumably because one looses entanglement due to correlations between the regions A and B which are now adjacent. We observe that the numerical results for p " π (which is equal to 10¨2π{L, giving a De Broglie's wavelength of 2), 2π (De Broglie's wavelength of 1) and 4π (De Broglie's wavelength of 0.5) are virtually indistinguishable from each other and from the analytic curve (2.1) here presented as the solid curve. The figure on the right in Fig. 5 shows very good agreement with the analytic predictions for three values of n and relatively high momentum (De Broglie's wavelength of 0.5). Similarly good agreement is found for states consisting of two identical (or not) excitations, whose Rényi entropies are presented in Fig. 6.

Replica Logarithmic Negativity
We have evaluated numerically the replica logarithmic negativity between regions A and B with various ratios of lengths. We have chosen A " r0, f RLs and B " rf L, f L`p1´f qRLs, where the ratio of lengths |A|{|B| " f {p1´f q is controlled by the parameter f taking values 1{2, 2{5 and 1{3, giving r 1 " r 3 " R{2 (f " 1{2), r 1 " 2R{5, r 3 " 3R{5 (f " 2{5), and r 1 " R{3, r 3 " 2R{3 (f " 1{3) respectively. The overall scale is again R, ranging from 0 to 1 by increments of 0.1.  Figure 7: Results for k " 1 and p " 4π with n " 2 (left) and n " 3 (right). For each case we present three curves corresponding to r 1 " r 3 " R 2 (circles), r 1 " 2R 5 and r 3 " 3R 5 (triangles), and r 1 " R 3 , r 3 " 2R 3 (squares). For n " 2 all three sets of data are identical, since in this case the negativity only depends on the sum r 1`r3 " R and equals minus the 2nd Rényi entropy. The symbols are numerical data and the curves are obtained from the formula (2.  Figure 8: Left: Results for k " 1 and p " 4π with n " 2 (dots), n " 3 (triangles), and n " 4 (squares). The dashed curve is the logarithmic negativity, showing once again the gradual change in curvature as n Ñ 1. The symbols are numerical data and the solid curves are obtained from the formula (2.3). The dashed curve is obtained from (2.5). Right: The replica logarithmic negativity for n " 3 and a state of two distinct excitations of momenta p 1 " 3π and p 2 " 4π. The circles give data for r 1 " r 3 " R 2 and the triangles correspond to r 1 " 2R 5 and r 3 " 3R 5 . Agreement with the analytic curves is again excellent. These correspond to twice the replica logarithmic negativity of a single excitation.
The regions' lengths are therefore 1 " f RL and 3 " p1´f qRL, whose sum is RL ranging from 0 to L. We have used different values of n and different types of excited state (with distinct and equal excitations). Here we concentrate on higher momenta, which are still well within the QFT regime but also well within the quasiparticle regime, indeed finding excellent agreement with our analytic predictions. In Figs. 7 and 8 (left) we present results for k " 1 and n " 2, 3 and 4 and different ratios r 1 {r 3 . Fig. 8 (right) and Fig. 9 explore some configurations for states with k " 2, 3 and states of equal/distinct excitations. In particular we find confirmation of the statement (2.12) that the replica logarithmic negativity of a state of k distinct excitations is k times the replica logarithmic negativity of the state of one excitation. : Left: n " 3 for a state of two equal particles k " 2 with p " 4π (triangles) and for a state of two distinct particles with momenta p 1 " 3π and p 2 " 4π (circles). The solid lines are the analytic formulae for ∆E 2

Conclusion
In two previous publications [32,33] we developed a methodology for the evaluation of the Rényi entropies (and related quantities) in zero-density excited states of one-dimensional free massive quantum field theory. Three methods were used: a finite-volume form factor approach based on the use of branch point twist fields and U p1q twist fields, a qubit picture showing that the same excess entanglement can be obtained from simple qubit states, and, finally, a numerical approach based on exact diagonalization for a harmonic chain. These methods have been employed again in the current work to extend results to the Rényi entropies of any number of disconnected regions and to the replica logarithmic negativity, a measure of entanglement for non-complementary regions. Our results demonstrate that all three approaches can still be successfully employed to evaluate the increment experienced by these more mathematically complex quantities in zero-density excited states. In the current paper we have focussed on the case of the free boson. It would be interesting to extend our results to the massive Majorana fermion.
We find that the increment of the Rényi, von Neumann, and single copy entropies of multiple disconnected regions of lengths i takes exactly the same form as for a single connected region of length , up to the replacement Þ Ñ ř i i . More precisely, in the scaling limit (1.9) the increment of the Rényi and related entropies is independent of the connectivity of the regions. It depends only upon their overall length.
For the replica logarithmic negativities we find that their increment is the logarithm of a symmetric polynomial of the variables r 1 and r 3 . By employing a qubit picture we have obtained closed formulae for all such polynomials (2.6). The numerical coefficients involved are given by (2.7) and point towards a combinatorial interpretation of the result. Indeed, we show in [38] that such polynomials are partitions functions for counting of particular families of graphs. This is also true for the Rényi entropies which may be recovered from the replica logarithmic negativity by taking either r 1 or r 3 to zero.
From the qubit interpretation of our results we are also able to generalize our conclusions to any number of disconnected regions and find that both the increments of Rényi entropies and logarithmic negativities do not depend on the connectivity of the regions, just on the overall size of each subsystem.
This work, together with [32,33] and the follow up paper [38], provides a complete understanding of the most popular measures of entanglement in zero-density excited states of free one-dimensional massive quantum field theory. Some of the results have also been shown to hold for interacting models [33], and all results hold as well in free higher dimensional theories [32,38], the most general proof being given in [38]. The main challenge now is to consider interacting models more generally and to show under which conditions (if any) these results still are valid in the presence of interactions. We expect to study this problem in a future work.

A Logarithmic Negativity for Identical Particles
The tables below give examples of the sort of polynomials that are obtained for the replica logarithmic negativities E k n pr 1 , r 3 q for a state of k identical particles. Recall that, as previously, r :" 1´r 1´r3 and r 1 , r 3 are the scaled versions of regions A and B sizes in Fig. 1. As indicated, the polynomials for n " 2 can be rewritten entirely in terms of r only. They are the same as those entering the expressions for the second Rényi entropies ∆S k 2 prq obtained in [32,33]. This is because for n " 2 we have T "T and so the two ratios of four-point functions in (1.10) and (1.11) are identical.

B Eigenvalues of the Reduced (Partially Transposed) Density
Matrix for k " 2 and k " 3.
It is instructive to repeat the analysis of subsection 3.1 for the case of two and three identical excitations.
(B.6) Thus, the Rényi entropies are given by above. In this contribution the coefficients A˘1 ,q once again are related to counting partitions into 0s and 1s with no consecutive 1s. The sum pλ t 7 q n`p λ t 8 q n`p λ t 9 q n is found (by inspection) to have the following structure: which corresponds to the p " 0 contribution in (2.6). The coefficients are positive and integer and they can be systematically computed from the explicit formulae for the eigenvalues λ t 7 , λ t 8 and λ t 9 , and/or by employing some useful properties of roots of third order polynomials as described for instance in [59]. The first few coefficients are given by and, finally, A 0,n " " 2 n for n odd 2`2 n for n even . (B.17) Interestingly, all the coefficients above are positive integers for n integer and they can be shown to agree with the formula (2.7), even though showing this analytically is relatively involved as there is a sum over partitions to perform. The replica logarithmic negativities are then given by E 2 n pr 1 , r 3 q " logpσ 1 pnq`σ 2 pnqq, (B.18) As usual, we are interested in the analytic continuation of this result to n " 1. In this case, as we know all roots explicitly we can simply compute the logarithmic negativity using its formal definition as the sum of the absolute values of the eigenvalues of the partially transposed, reduced density matrix. We have that lim nÑ1 σ 1 pnq " For the remaining eigenvalues it is possible to use properties of cubic roots to show that in our case, lim nÑ1 σ 2 pnq " namely, it can be shown that λ t 7 ă 0 whereas λ t 8 and λ t 9 are positive. ∆ and ∆ 0 were defined in (2.9), (2.10). The logarithmic negativity is then as given by (2.6).

(B.33)
In summary, where now H 1 " H 1 b¨¨¨b H α . Again, we can define an orthonormal basis t|ky P H : k " pk 0 , . . . , k α q P N α`1 0 u such that k m is the number of excitations in region R m and ř α m"1 k m " k for a state of k identical excitations. Let where r m is the scaled length of region R m . We define the qubit state: |Ψ pkq α y " ÿ k"tk 0 ,...,kαuP σ α`1 where σ α`1 0 pkq represents the set of integer partitions of k into α`1 non-negative parts. It is easy to extend (1.2) to the case of α disconnected regions. Indeed, if we now introduce S k n pr 1 , . . . , r α q :" exp tp1´nqS k n pr 1 , . . . , r α qu where S k n pr 1 , . . . , r α q is the Rényi entropy for the state (C.2) it is easy to see that this can be written as: (C.7) Again the multinomial coefficients constrain the sums. The presence of q! in the denominator means that q ě 0. We know also that k i m must be non-negative for all i " 1, . . . , n and m " 1, . . . , α´1. Furthermore the only non-zero terms in the sums are given by k´q´ř α´1 m"1 k i m ě 0 and thus q ď k. In summary, for the same reasons as in the two region case 0 ď q ď k.
We can re-write (C.7) as: S k n pr 1 , . . . , r α q " notice it takes the same form as (2.2) (with r replaced by r 0 ). Therefore, the Rényi (and related) entropies of this particular class of qubit states depend only on the relative size of the two parts in the bipartition and not on whether or not they are connected.

C.2 Replica Logarithmic Negativities of Disconnected Regions
A very similar computation can be performed for the replica logarithmic negativities. The starting point is the assumption that regions A and B are now disconnected, namely so that A consists of a number α and B of a number β of disconnected regions. Let regions A i and B i have scaled lengths given by r i 1 and r i 3 , respectively. Then our results for the (replica) logarithmic negativities will still hold up to the identifications: and r " 1´r 1´r3 . (C.14) In the qubit picture this can be shown in a very similar way as for the Rényi entropies in the previous section, so we do not present the computation here.

D Form Factors and Residue Formulae
Let upβq and vpβq be generic functions of β such that vpβq has a single zero at β " θ. The corresponding residua of the following first, second, and third order poles are given by Res β"θ u pβq pβ´θq v pβq " u 1 pθq v 1 pθq´u pθq v 2 pθq 2 rv 1 pθqs 2 , (D.2) Res β"θ u pθq These generic formulae have been employed in order to obtain (4.57). In particular, the functioñ gprq defined in (4.58) results from the evaluation of residua of third order poles of the type and the denominator factors are cancelled by the form factor residues in the calculation.