Wedge-local fields in integrable models with bound states

Recently, large families of two-dimensional quantum field theories with factorizing S-matrices have been constructed by the operator-algebraic methods, by first showing the existence of observables localized in wedge-shaped regions. However, these constructions have been limited to the class of S-matrices whose components are analytic in rapidity in the physical strip. In this work, we construct candidates for observables in wedges for scalar factorizing S-matrices with poles in the physical strip and show that they weakly commute on a certain domain. We discuss some technical issues concerning further developments, especially the self-adjointness of the candidate operators here and strong commutativity between them.


Introduction
In recent years, we have seen many interesting developments in constructing models of quantum field theory (QFT) in the operator-algebraic approach [25]. Among them, one of the most important contributions was the construction of 1 + 1 dimensional scalar quantum field theories with factorizing S-matrices [27,29]. In that work, Lechner took a large family of analytic functions that satisfy certain conditions and constructed quantum field theories that have these functions as the two-particle S-matrix. This was remarkable because, although physicists conjectured several properties and computed many interesting quantities, a mathematically consistent construction of the models in axiomatic approaches had not been obtained before, and also because these models include many S-matrices that had not been considered by physicists, as the arguments exploit only a few properties of the S-matrices and do not depend on specific expressions or computations. Now, it was clear at the first step of the construction [27] that this method does not apply directly if the function S(ζ ) of the input has simple poles in the physical strip, 0 < Im ζ < π. On the other hand, some integrable QFTs are believed to have S-matrices with simple poles and they should correspond to bound states of elementary particles. The purpose of this paper is to extend the step of [27] to the cases with poles in the physical strip.
In the models constructed in [27,29], if n particles are incoming, then after the scattering process n particles are outgoing. Furthermore, the scattering process of n particles can be systematically constructed from that of 2 particles. In this case, the S-matrix is said to be factorizing. This property is expected for so-called integrable models, which have infinitely many conserved currents. In some cases, for example the sine-Gordon model, Fröhlich and Seiler constructed Euclidean Green's functions [23] and proved that the S-matrix is nontrivial. But, to the authors' knowledge, it is unknown whether it is factorizing. More general integrable models have more complicated Lagrangians, and the understanding of these models is limited to the perturbation theory.
An alternative approach to integrable models is the form factor program [3,41]. In this program, instead of quantizing the fields, one conjectures the S-matrix from the symmetry of the Lagrangian. Then the matrix components of local fields (form factors) are computed and then the n-point function should be reconstructed from these matrix components. Although many interesting quantities have been computed, a convergence proof of the form factor expansion of the n-point functions is still lacking in almost all cases [3].
The recent operator-algebraic approach, initiated by Schroer [40], considers so-called wedge-local fields, rather than point-like quantum fields. Wedge-local fields are observables localized in an infinitely extended, wedge-shaped region. Indeed, this infinite extension allows one to take operators that are very simple in the momentum space [27,33]. Thereafter, the algebra of local, finitely extended observables should be defined as the intersection of two algebras corresponding to left-and right-wedges [29]. It is noteworthy that the existence proof of local observables avoids explicit computations of quantum fields, but it is reduced to a certain phase space property called modular nuclearity [13,14]. The resulting net of local observables reproduces the S-matrix of the input; therefore the inverse problem for a class of S-matrices has been solved. Now, let us recall that the S-matrices treated in [27,29] must have analytic components in the physical strip. In the operator-algebraic approach, there have been many other constructions of algebras corresponding to wedges [1,7,8,15,21,24,30,32,43,44] and some related constructions [6,31,35] which take an analytic function as an input, but it was always assumed that the function has no pole in a certain region. On the other hand, a pole in the physical strip of the S-matrix components is considered to correspond to a bound state. Indeed, some integrable models, e.g.the sine-Gordon model, are believed to have S-matrix with poles in the physical strip. These S-matrices have not yet been treated in the operator-algebraic approach. For an S-matrix S without poles, Lechner constructed a pair of operator-valued distributions φ, φ and proved that they commute when they are smeared by test functions supported in the left-and right-wedges, respectively. This computation involves a shift of the integral contour of a function containing S. If one takes the same construction for S with poles in the physical strip, the shift of the integral contour yields the residues of S at these poles. Hence, the fields φ, φ themselves cannot be wedge-local.
The same problem has appeared also in the form factor program [3]. Solutions of the form factor equations should represent the matrix components of local operators, and their commutators should vanish when they are spacelike separated. Indeed, such a formal proof has been given first by Smirnov for S-matrices without poles [41]. A crucial part of the proof is again done by shifting integral contours, which is invalid when the S-matrices have poles. For S-matrices with poles, they added further properties to form factors which make a correspondence between the poles and "bound states". Quella showed that these new properties allow one to cancel the residues which come from first-order poles [36]. Higher poles are discussed by Babujian et al. [5] and local commutativity appears to be formally maintained, at least for some specific models (the so-called Z (N )-Ising models).
In this paper, for a certain S-matrix with poles in the physical strip, we introduce a new field φ = φ + χ , by adding an operator χ , which we call the bound-state operator, to the field φ of Lechner. It is the commutator of χ with its reflected operator χ that cancels the contribution of the residues coming from the commutator between φ and φ , mentioned above. The operator χ has a formal integral expression in terms of Zamolodchikov-Faddeev operators z, z † with complex arguments.
To prove wedge-commutativity in operator-algebraic approach, one considers the reflected field φ using the action of the CPT operator and shows that its commutator with φ vanishes for test functions supported in right-and left-wedges, respectively, in the weak sense, namely, the matrix elements of the commutator between suitable vectors vanish. Strong commutativity remains to be proven. Yet, by assuming the existence of nice self-adjoint extensions, we give an argument for the Reeh-Schlieder property. We also argue that our wedge-local fields are non-temperate polarization-free generators [10]. Besides, we classify the scalar S-matrices with poles in the physical strip that satisfy the requirements of our analysis. The general form of the S-matrix that we obtain is essentially given by a certain subclass of S-matrices known from [28] multiplied with a universal model independent factor which has poles in the physical strip. This especially includes the S-matrices of the Bullough-Dodd model [22]. Further, we clarify how the (weak) wedge localization of the field φ is related to the properties of form factors [4]. We conjecture a generalization of the characterization theorem of local observables in [12] to models with bound states. In [12] local observables are expanded into a series in terms of Zamolodchikov-Faddeev operators, where the expansion coefficients are related to the form factors of the observable. We sketch the outline of a proof that, if its expansion coefficients of an operator fulfill a set of conditions, slightly modified from [12], then it formally commutes with our fields φ( f ); hence they are local to each other.
The paper is organized as follows: In Sect. 2 we recall the results of Lechner [27,29] for S-matrices analytic in the physical strip, and we introduce our general notation. In Sect. 3 we summarize the properties of scalar S-matrices with poles in the physical strip in the models under investigation, and we construct the wedge-local fields φ, φ . A proof of weak wedge-commutativity is also given in this section. In Sect. 4 we show that our wedge-local fields are compatible with the form factor program in the sense explained above. Further, we explain how the form of the bound-state operator can be deduced by formal arguments. In Sect. 5 we present our conclusions and open problems. Appendix A is dedicated to the classification of scalar S-matrices fulfilling the properties introduced in Sect. 3. Appendix C comments on the problem of finding suitable self-adjoint extensions of the field φ.

Background: Haag-Kastler nets and wedge-local field in two dimensions.
Here we review the motivation to study wedge-local fields, the main objects which we construct in this paper.
In the operator-algebraic approach to quantum field theory (QFT), a model of QFT is realized as a net of algebras of observables.  [25].
If one has a Wightman field φ with certain regularity conditions, then one can construct the corresponding net by defining A(O) := {e iφ( f ) : supp f ⊂ O} , where M means the set of bounded operators commuting with any element of M. The double commutant M is the smallest von Neumann algebra which includes M. Actually it is required that φ( f ) and φ(g) have commuting spectral projections for f, g with spacelike separated support, and this follows from the regularity condition. In this way, a Haag-Kastler net is considered as the operator-algebraic formulation of quantum field theory.
One of the difficulties in constructing Haag-Kastler nets lies in the infiniteness of the family {A(O)} which needs to comply with the axioms. Instead, Borchers observed that for d = 2, actually the whole net can be recovered from the single von Neumann algebra A(W R ) associated with the (right-)wedge-shaped regions W R := {a ∈ R 2 : a 1 > |a 0 |} and the spacetime symmetry U (under a condition called Haag-duality). Furthermore, by the Tomita-Takesaki theory of von Neumann algebras [42], it is enough to know the restriction of U to the translation subgroup R 2 [9].  (D a,b )) . Then one can show that this "net" A, a collection of algebras, satisfies isotony and locality. Furthermore, the representation T extends to a representation U of the Poincaré group which makes A covariant and is still invariant. In this way one obtains a "net" (A, U, ), where the only missing property is that is Hence, in the operator-algebraic approach, the construction of Haag-Kastler nets can be split into two steps: (1) to construct Borchers triples, (2) to prove the cyclicity of . In the following, we exhibit an attempt to (1). For this purpose, we construct wedgelocal fields. Wedge-local fields (φ, φ ) is a pair of operator-valued distributions such that [φ( f ), φ (g)] = 0 if supp f ⊂ W L and supp g ⊂ W R . Then it is natural to expect that M = {e iφ ( f ) : supp g ⊂ W R } , together with appropriate U and , gives a Borchers triple (actually, this last step is our open problem).

Zamolodchikov-Faddeev algebra and wedge-local fields.
We consider quantum field theories in 1 + 1 dimensional Minkowski space (and with the convention x · y = x 0 y 0 − x 1 y 1 ) with factorizing S-matrices, which are characterized by the particle spectrum and the two-particle scattering function as the main input in the theory. In the following, we introduce the mathematical framework and the notation we will use to describe these models, following [27,29]. In particular, this subsection is meant to be an overview on previous work on the topic [27] and to be an introduction to Sect. 3, where we will consider models with scalar two-particle S-matrices which have poles in the physical strip R + i(0, π).
2.2.1. Single-particle space, S-symmetric Fock space, space-time symmetries. We parametrize the momentum of a single particle with mass m > 0 by the rapidity θ : The two-particle scattering function is generally a complex-valued meromorphic function S : R + i(0, π) → C with a certain number of symmetry properties. These symmetry properties are the well-known properties of unitarity, hermitian analyticity and crossing symmetry (see for example [26]) typically fulfilled by any twoparticle scattering function in a local integrable quantum field theory. We summarize these properties in Sect. 3.1.
With these functions we can define the Hilbert space H of the theory. In the case of models with only one species of massive scalar particle, the single particle Hilbert space is H 1 = L 2 (R, dθ). We define the n-particle Hilbert space H n as the subspace of S-symmetric functions in L 2 (R n ) = H ⊗n 1 and the Hilbert space of the theory as H := ⊕ ∞ n=0 H n with H 0 = C . We introduce the orthogonal projection P n := 1 n! σ ∈G n D n (σ ) thus we can write H n = P n H ⊗n 1 , and we denote with D the dense subspace of H of vectors with finite particle number.
Furthermore, the Hilbert space H is endowed with a unitary representation of the proper orthochronous Poincaré group, denoted by P Note that this action of the (proper) Poincaré group preserves the S-symmetrized spaces H n . We will adopt the following convention for the one-particle wave function associated with g ∈ S (R 2 ) [27]: If g is supported in W R , then g + (θ ) has an entire analytic continuation, which is bounded in R + i(−π, 0). Furthermore, let g j (x) := g(−x). Then, if g is real, (g j ) + (θ ) = g − (θ ) and g + (ζ ) = g − (ζ ). Finally, the proper Poincaré group acts on R 2 and also on the space of test functions naturally, which we denote by g (a,λ) (the space-time reflection acts by g → g j ). One can easily check that it is compatible with the action on the one-particle space: (g (x,λ) ) + (θ ) = U 1 (a, λ)g + (θ ).

Borchers triples for analytic S-matrices
For the class of two-particle scattering functions S(θ ) which are analytic in the physical strip θ ∈ R + i(0, π), local observables associated with wedge-regions, say with the standard left wedge W L , can be constructed by following an argument due to Schroer [39] and Lechner [27]. Specifically, they define a quantum field φ as We note that this reduces to the free field if S(θ ) = 1. In general, the field φ shares many properties with the free field as shown in [28, Proposition 4.2.2]. In particular, it is defined on the subspace D of H of vectors with finite particle number and it is essentially self-adjoint on D for real-valued f (we denote its closure by the same symbol φ( f )). It has the Reeh-Schlieder property, it solves the free Klein-Gordon equation and it transforms covariantly under the representation U (x, λ) of the proper orthochronous Poincaré group. The only exception is the property of locality. The field φ(x) is not localized at the space-time point x in the usual sense, but rather in an infinitely extended wedge with tip at x, W L + x. To make this more precise, we introduce the "reflected" Zamolodchikov-Faddeev operators, It has been shown in [27,Proposition 2] that the two fields φ, φ are relatively wedgelocal, in the sense that the commutator [e iφ( f ) , e iφ (g) ] is zero for any real-valued test functions f, g with supp f ⊂ W L and supp g ⊂ W R , Hence, we can interpret φ, φ as observables measurable in the wedges W L , W R , respectively. This result can be obtained by computing the commutators of z # with z # as shown in [28, Lemma 4.2.5] and by shifting a certain integral contour which critically uses the analyticity of the two-particle scattering function S(θ ) in the physical strip θ ∈ R + i(0, π).
It should be remarked [27,Proposition 2] that also the properties of the test functions f, g play an important role in the proof of wedge-locality. More specifically, the proof uses the fact that if f ∈ S (W L ) (similar arguments apply to g as well) then its Fourier transform f + fulfills certain analyticity, boundedness and symmetry properties in the strip R + i(0, π).
Starting from the fields φ, φ one can then define the corresponding von Neumann algebras. The right-wedge algebra is given by and other wedge algebras are defined by using the action of translations and reflection, The most important consequence of commutativity between φ and φ is that is separating for M. The strictly local observables can then be recovered from the intersection of a right and left wedge algebras, as explained in Sect. 2.1. Furthermore, Lechner proved [29] that the algebras A(O) are nontrivial, namely, are different from C1. In order to show this, he proved modular nuclearity condition, a spectral property of the modular operators associated with the theory. He showed that under a certain regularity condition on the scattering function S for all regions O, at least with a minimal size [2], the vacuum state is cyclic and separating for the algebra A(O).
For the class of two-particle scattering functions S(θ ) which are not analytic in the physical strip θ ∈ R + i(0, π), the fields φ( f ), φ ( f ) fail to be wedge-local and the situation becomes more complicated, as we will see in Sect. 3.2.

Poles in the S-matrix and bound states.
We now suppose that the two-particle scattering function S(ζ ) has poles in the physical strip ζ ∈ R + i(0, π). Physically, these poles are related to the notion of "bound state", which here is interpreted as the "fusion" of two bosons, and, although the S-matrix is factorizing, their components are analytically related to each other. Let us recall the physicists' arguments (see, for example, [19]). This section only tries to physically motivate the conditions on the Smatrices in Sect. 3.1.
Two elementary particles can be fused into a bound state if the total momentum of the two particles, say "1" and "2", lies on the mass shell of a third particle, say "b" [36], namely where m b is the mass of the third particle and Therefore, one can find ζ b such that the momenta of the particles are related by where ζ 1 , ζ 2 and ζ b are the (possibly complex) rapidities of the two fusing bosons and of the bound particle, respectively.
To determine the rapidities of the particles involved and the position of the pole in the rapidity complex plane, we essentially need to solve Eq. (1). In preparation for Sect. 3, we will do that in a simpler case by considering a system with only one species of particle. In that case, the fusion process becomes quite simple: Two bosons of the same species fuse to form another boson of the same species, meaning also that the masses of the particles are equal.
To solve Eq. (1), we make the ansatz that the difference of the rapidities of the fusing bosons is purely imaginary, that is, Hence, we can parametrize the rapidities of the two fusing particles and of the bound particle as ζ 1 = θ + iλ 1 , ζ 2 = θ + iλ 2 and ζ b = θ , with θ real. Using this parametrization for the momenta of the particles, one can show that there is a unique solution (up to addition of 2πi) to Eq. (1) given by λ 1 = π 3 and λ 2 = − π 3 . Demanding that the difference λ 1 − λ 2 is in the physical strip, we obtain λ = 2π 3 . In the physical literature, one associates bound states with the poles of the S-matrix [19]. Specifically, we assume that the two-particle scattering function S(ζ ) has a simple pole at the point iλ (the so called s-channel pole), and we denote its residue by R := Due to the property of crossing symmetry (see (S3) in Sect. 3.1) of the scattering function, S(ζ ) has another pole at πi 3 (the so called t-channel pole) with residue R := res ζ = πi 3 S(ζ ). Moreover, one can show that as a consequence of hermitian analyticity (S2) and unitarity (S1) of the scattering function, the residue R is purely imaginary and one has R = −R again by crossing symmetry of S. We assume that Im R > 0. Correspondingly, unitarity (S1) also implies that S(ζ ) has zeros at the points ζ = − 2πi 3 , − πi 3 . Finally, we will assume that except for the two simple poles at 2iπ 3 , πi 3 there are no other poles of S in the physical strip. An example of a 1 + 1-dimensional integrable model with scalar S-matrix fulfilling these properties is the Bullough-Dodd model [22].
A particle in these models is considered as the bound state of two of the same species, and the two-particle S-matrix S satisfies a nontrivial equation, the so-called bootstrap equation (S4).

Scalar Two-Particle S-Matrices
Here we present our wedge-local fields associated with scalar S-matrices with poles in the physical strip. We specify the properties of these S-matrices and give examples. Then an important operator χ( f ), which "binds f and the state", is introduced. We use this operator in order to construct our wedge-local fields.

Properties of scalar S-matrix with poles.
Let us consider the simplest case where the one-particle space has multiplicity one. Since we are interested in S-matrices with poles in the physical strip, the bound state of two "elementary" particle must be the same particle, therefore must have the same mass. As discussed in Sect. 2.3, one observes that the only possibility for the pole in s-channel is 2πi 3 . By crossing symmetry, the S-matrix must have another pole at πi 3 . Furthermore, the residue of S at 2πi 3 must be a positive multiple of the imaginary unit i. This is interpreted as a consequence of hermiticity of the Hamiltonian [18]. 1 Summarizing, we assume that our two-particle S-matrix S is defined on R + i(0, π) except for poles indicated below, has therefore L ∞ -boundary values at R and R + πi, and has the following properties (c.f. [33]). Actually, we will see in Appendix A that (S6) follows from (S1)-(S5). Note also that we consider bosonic particles, yet the observables are represented on the S-symmetric Fock space (see below).
We introduce η = i √ 2π |R|, in accordance with the literature [36], recalling that R is the residue of S at ζ = 2πi 3 . Examples Although specific expressions of S are not needed in our main construction, we present here a family of examples, in order to show that the above set of axioms is not empty.
The simplest examples of S which are believed to be associated to the Bullough-Dodd model [22] are the following: satisfies all properties of the S-matrix but positivity of residue. Indeed, we have The remaining factor also satisfies (S1)-(S4) except for positivity of residue, namely, it has no pole in the physical strip and satisfies has a positive residue at 2πi 3 , and therefore, satisfies (S5). (S6) is now straightforward. From the above computations, it is also clear that a product where 0 < B k < 2, B k = 1 and n is odd, satisfies all the properties of S-matrix. When n > 1, no Lagrangian is known for such two-particle S-matrix. We will completely classify all S-matrices which comply with (S1)-(S6) in Appendix A.

The bound-state operator.
For a test function f supported in W L , let us introduce an unbounded operator χ( f ) on the S-symmetric Fock space H. This operator will preserve the particle number and will be interpreted as the operator which "makes a bound state". Actually, in a model with scalar S-matrix, a bound state particle of two "elementary particles" is again the same elementary particle as there is only one species of particle [22].
Let us denote its component on H n by χ n ( f ). Firstly, χ 0 ( f ) annihilates the vacuum .
For ξ ∈ H 1 = L 2 (R, dθ), we say that ξ(θ) has an L 2 -bounded analytic continuation on a strip parallel to R (e.g. is an analytic function on that strip (with boundary value ξ(θ) at Im ζ = 0) such that for each fixed − < α < 0 (respectively 0 < α < ) the function θ → ξ(θ + iα) is an L 2 -function in θ , with uniform L 2 -bound in α. The action χ 1 ( f ) on H 1 is given as follows: where η and R are given in Sect. 3 where Dom(χ ( f )) is the algebraic direct sum of Dom(χ n ( f )), hence is a subspace of D and, of course, the domain of the product AB of possibly unbounded operators A, B is given by {ξ ∈ Dom(B) : Bξ ∈ Dom(A)}. We will discuss the question of self-adjoint extensions of these operators in Appendix C.
The operator , it holds that D n (σ )P n = P n and the n-particle component χ n ( f ) can alternatively be expressed as follows: Note that, if n is S-symmetric and in the domain of χ 1 ( f )⊗1⊗· · ·⊗1, n (θ 1 , . . . , θ n ) has a meromorphic continuation in θ k and it holds that Therefore, each term in the last expression of χ n ( f ) above can be further written as follows for k ≥ 2 and for n = P n n : where we reordered the variables in the third equality, and we used the bootstrap equation in the last equality. Note that, although S-factors have poles, these poles are cancelled by the zeros of n by the definition of the domain of χ n ( f ) in Eq. (2) and the whole expression remains L 2 (this is the meaning that n ∈ Dom(χ n ( f )). In this expression, the one-particle component χ 1 ( f ) acts on each variable of n up to a correction of S factors. As the one-particle action (2) realizes the idea that the state of one elementary particle ξ is fused with f + into the same species of particle, as in Sect. 2.3, we might call it the "bound state operator". Similarly, for a test function g supported in the right wedge W R , we introduce the reflected bound state operator χ (g): For a real g, this operator is indeed related to χ by the CPT operator J : To see this, let us consider the one-particle components. Since J 1 acts as the complex conjugation, it takes an analytic function in the lower strip to an analytic function in the upper strip, therefore the domains of χ (g) and of J χ(g j )J coincide. Recall that for a real g, As J n commutes with P n , we have χ n (g) = J n χ n (g j )J n . Since the whole operators χ(g) and χ (g) are defined as the direct sum, the desired equality follows.
This operator χ (g) has an alternative expression as χ( f ) does: where ρ k = τ n−k+1 τ n−k+2 · · · τ n−1 are the cyclic permutations and Let us check basic properties of χ( f ). As χ (g) is defined similar to χ( f ), they share a number of properties. The following propositions have the obvious counterparts for χ (g). Proof. First let us look at χ 1 ( f ). By its definition, it is densely defined. To see that it is symmetric, we take two vectors ξ, η ∈ Dom(χ 1 ( f )) which have smooth Fourier transforms with compact support. Such vectors form a core for χ 1 ( f ). Indeed, we can , where x f is the multiplication operator by f + θ + πi 3 and 1 ξ(θ) = ξ(θ − 2πi) with the obvious domain. As it 1 implements the real shift ( it 1 ξ )(θ ) = ξ(θ + 2π t), the Fourier transform diagonalizes 1 , and the vectors above form a core. Then, ξ, η are the Fourier transforms of compactly supported smooth functions, therefore they are rapidly decreasing. Recall that f + (θ + iα) is rapidly decreasing as well, if 0 ≤ α ≤ π . Furthermore, η(ζ ) is analytic in ζ . Now we compute: where we used the Cauchy theorem in the second equality and f As this equality holds on a core of χ 1 ( f ), so does it on the whole domain and we obtain the symmetry of χ 1 ( f ). Let us check that χ n ( f ) is densely defined. As P n is bounded, it is enough to see that (χ 1 ( f ) ⊗ 1 ⊗ · · · ⊗ 1)P n is densely defined. The range of P n is the set of S-symmetric functions, while the domain of χ 1 ( f ) ⊗ 1 ⊗ · · · ⊗ 1 is the functions which has an L 2 -bounded analytic continuation to πi 3 in the first variable. Let us take an arbitrary set {ξ 1 , . . . ξ n } of n vectors in the domain of χ 1 ( f ). Then, is of course S-symmetric, but has poles which come from the S-factors. Such vectors are a dense subspace of P n H n . Note that all these poles come from the poles of S at θ j − θ k = πli 3 , l = 1, 2 and 1 ≤ k < j ≤ n. In order to compensate these poles, we can multiply it (n − 2)! times by where α < 0 or π < α. This is a bounded invertible symmetric function in the real variables θ θ θ . Therefore, the multiplication operator M C n by C n preserves H n and its image of a dense subspace is again dense. The functions in the image have an analytic continuations in the first (actually any) variable to πi 3 and are S-symmetric, therefore in the domain of (χ 1 ( f ) ⊗ 1 ⊗ · · · ⊗ 1)P n . Now the symmetry of χ n ( f ) follows from the symmetry of χ 1 ( f ) and a general fact (ABC) * ⊃ (C * B * A * ). The whole operator χ( f ) is by definition the direct sum of these operators, therefore densely defined and symmetric as well.
In Sect. 2.2, we saw that there is an action U of the Poincaré group on H. Let us check that χ( f ) is covariant with respect to U . L and (a, λ) ∈ P ↑ + such that a ∈ W L . Then it holds that Ad U (a, λ)(χ ( f )) ⊂ χ ( f (a,λ) ).
Next we take a pure translation (a, 0), where a ∈ W L . As U 1 (a, λ) * multiplies by e −ia· p(θ) , which has a bounded analytic continuation in R+i(−π, 0), the domain of χ 1 ( f ) is preserved. Again, recall the one-particle action (U 1 (a, 0) f (a,0) )), it holds that As we have f (a,0) )ξ )(θ ), and this is the covariance with respect to translations.

Formal expression
We argue that our operator χ( f ) defined above is formally equivalent to This formal expression preserves the particle number, because it consists of a creation operator z † and an annihilation operator z. In this argument, we do not pay much attention to domains. We will not use this formal expression later in proofs, yet it is interesting to observe that the operator χ( f ) has such a simple expression. We know that z( f ) = Pa( f )P and z † ( f ) = Pa † ( f )P, where P = n P n is the projection from the unsymmetrized Fock space to the S-symmetric Fock space. Let us take an S-symmetric vector n ∈ H n . The unsymmetrized annihilation operator a(ζ ) substitutes the first variable by a fixed number ζ . As n is an S-symmetric function with n variables, then after the action of a, it is still an S-symmetric function with n − 1 variables. This means that we can remove several factors of P in the above formal expression and obtain Let us look at the integrand formally. For a fixed θ , the action of the annihilation operator gives Thereafter, the action of the creation operator gives the following: , θ 2 , . . . , θ n and after the integration and the multiplication by −iη, we get If we look at the action on the first variable, this is exactly nχ 1 ( f ). Recall that χ n ( f ) = n P n (χ 1 ( f ) ⊗ · · · ⊗ 1)P n . As n was arbitrary, we get finally which is the desired expression. Similarly, we have

The wedge-local fields.
We define the (left) wedge-local field by The domain of φ( f ) includes that of χ( f ), therefore the domain of φ( f ) coincides with the latter. The reflected field is given by Proposition 3.3. Let f be a real test function supported in W L . Then the operator φ( f ) defined above has the following properties. f (a,λ) ). Proof. (1) We saw in Proposition 3.1 that χ( f ) is symmetric, and we know from [27, is symmetric as well.
Of course, similar properties hold for φ (g). Compared with [27, Proposition 1], our fields are of a subtler nature. The operator φ( f ) does not preserve its domain (see Appendix C), neither is essentially self-adjoint.

Weak commutativity
Here we prove our main result. As the domains of φ( f ) and φ (g) are subtle, we cannot always form the products φ( f ) φ (g) and φ (g) φ( f ), and the commutator [ φ( f ), φ (g)]. Instead, we consider the weak form of commutation. Consider the following linear space of vectors: The intersection Dom( φ( f )) ∩ Dom( φ (g)) has an interesting property. The n-particle component n of a vector has an analytic continuation in θ 1 in the negative imaginary part, but by S-symmetry, it can be analytically continued in any variable (with possible poles at the poles of S). At the same time, it has analytic continuation in all the variables in the positive imaginary part, as it is in Dom( φ (g)). Therefore, n is not only the L 2 -boundary value of an analytic function, but it has a continuous value at θ j ∈ R. Now, by (S6), S(0) = −1 and n is S-symmetric, therefore n (θ 1 , . . . , θ n ) has a zero at θ j − θ k = 0. In the expression below, j S θ − θ j + πi 3 n (θ, θ 1 , . . . , θ n−1 ), the S-factors have simple poles at θ − θ j = 0, and they are cancelled by the zeros of n (θ, θ 1 , . . . , θ n ). Therefore, it still has pointwise meaning. We conjecture that the additional condition that it has an L 2 -bounded analytic continuation in θ should be automatic, and the linear space above should be simply the intersection Dom(χ ( f )) ∩ Dom(χ (g)). Yet, as the domains are a subtle question, we content ourselves in the present work by showing the weak commutativity on this space. One can easily check that it is dense in H by considering the vectors constructed in Proposition 3.1. For vectors , in the above space, we show weak commutativity, which is, From weak commutativity and some estimate of operators, strong commutativity can follow (e.g. [20]). Unfortunately, we have not yet been able to establish such estimates and in this work we restrict ourselves to these weak relations. See also Appendix C.
Theorem 3.4. Let f and g be real test functions supported in W L and W R , respectively. Then, for each , in the linear space above, it holds that Proof. Recall that the vectors , have finitely many non-zero components and our operators acts componentwise: In this proof, it should be kept in mind that and are already S-symmetric.
Let us compute the commutator [ φ( f ), φ (g)], which by the above expands into several terms that we will compute individually.
We can actually compute this commutator in the operator (not the weak) form, as we do not encounter the problem of domains. Recall that, if g is real, J 1 g − = g + .
Let us look at the expression we derived early in Sect. 3.2, in which χ n ( f ) is written as a sum of n terms. Its first n − 1 terms act in the first n − 1 variables in the same way as n − 1 terms in χ n−1 do. On the other hand, z (g + ) acts only on the n-th variable θ n . Therefore, the only remaining term in this commutator is the one which comes from the n-th term in the expression, which is For the later convenience, we further rewrite this using the S-symmetry of n : where in the second equality we used the bootstrap equation, in the third equality the assumed property of explained before this Theorem 3.4 and Lemma B.2, and the last equality follows because g is real.
This can be computed in a similar way as before. By paying attention that z reduces the number of variables and shifts the indices of the remaining ones, the result is: where we used the assumed domain property of , Lemma B.2, the crossing symmetry of S, and that f is real. This coincides with the result of the commutator [χ( f ), z (J 1 g − )] up to a sign, therefore they cancel each other.

The commutators [z † ( f + ), χ (g)] and [χ( f ), z † (g + )]
One can show that these two commutators cancel each other by taking adjoints and repeat the computation as in the commutators before. More precisely, before we have computed weak commutators such as K , Z − Z * , K , where K corresponds to either χ( f ) or χ (g) and Z is the annihilation operators z( Here we must compute numbers such as K , Z * − Z , K , which is complex conjugate to the previous one.

The commutator [φ( f ), φ (g)]
This part has been essentially done in [27, P.13]. The difference is that S has two poles in the physical strip, therefore we obtain residues when we shift the integration contour. This commutator preserves the particle number, therefore, it suffices to compute its action at fixed n. Let us take , with only n particle components. By considering that the poles of S are simple, the result is (note that we switched f and g from [27]) where we used that the residue of S at πi 3 is −R, which follows from the crossing symmetry. In order to justify the second equality, we should note that if any pair of θ k 's does not coincide, then the integrand has a simple pole and we have the equality. The complement of such θ θ θ 's has Lebesgue measure zero, therefore we obtain the equality. Actually, although this expression looks unbounded as a function of θ k because of the poles of S, it is actually bounded since these poles at θ k − θ j = 0 cancel each other. (This is not surprising, since φ( f ) and φ (g) are bounded on H n .)

The commutator [χ( f ), χ (g)]
Now we come to the most important part of the computations. As both χ( f ) and χ (g) preserves each n-particle space, we assume again that , have only n-particle components.
Let us recall the expressions of χ( f ), χ (g) we derived in Sect. 3.2: Therefore, the term χ (g) , χ ( f ) of the commutator is Actually, ρ k here can be replaced by any permutation σ such that σ (1) = k. Let us choose the transposition τ 1,k : (1, k) → (k, 1). τ 1,1 coincides with the unit element of S n . Similarly, we can use the transposition τ n− j+1,n instead of ρ j . Then it is clear that τ 1,k and τ n− j+1,n commute unless k = n or j = n or k = n − j + 1. For pairs of k and j such that the two transpositions commute, the scalar product reduces to the following: where the last equality follows by Lemma B.1 (by writing χ 1 ( f ) = √ 2π |R|x f 1 6 1 , etc.). If k = n, j = 1, n (respectively k = 1, n and j = n), one has τ n− j+1,n τ 1,n = τ 1,n− j+1 τ n− j+1,n (resp. τ 1,n τ 1,k = τ 1,k τ n,k ), and we can reduce these contributions to the same value (3). The case k = j = n is also easy and gives the same contribution. We will see that they get cancelled by the term The remaining terms are those with k = n − j + 1. Here, χ 1 ( f ) and χ 1 (g) act on the same variable: where we used the assumption that the product of , and the S factor remains L 2 , Lemma B.2 and the condition (S5) which implies that R = i|R|. Note that this is the second term in the commutator [χ( f ), χ (g)], therefore it gets another minus sign. The expression so obtained is equal to the first contribution from the commutator [φ( f ), φ (g)] with reversed sign and coupling with , therefore, they cancel each other.
Let us examine the term The computation is similar to the previous case and we obtain the following. As before, for pairs of k and j such that k = n − j + 1, one has the contributions: and this cancels the contribution in Eq. (3).
The remaining terms in χ( f ) , χ (g) are: . . . , θ k , . . . , θ n ) The last expression is equal to the second contribution from the weak form of the commutator , [φ( f ), φ (g)] up to the sign, therefore, they cancel each other. Altogether, we have seen that all the terms in [ φ( f ), φ (g)] cancel.

Reeh-Schlieder property
As the fields φ( f ) and φ (g) do not preserve their domains, especially one cannot iterate them on the vacuum more than once (see Appendix C), the Reeh-Schlieder property does not hold for polynomials of these fields. Instead, if we assume the existence of nice self-adjoint extensions, we can argue as below. But as we do not have such extensions, we refrain from stating it as a theorem. Let us suppose that for each g there is a self-adjoint extension of φ (g), which we denote by the same symbol, such that φ is covariant with respect to U . Suppose also that, for each f , φ( f ) has a nice self-adjoint extension, such that φ (g) and φ( f ) strongly commute. We consider the von Neumann algebra and we have to show that M = H. Actually, as M is an algebra of bounded operators containing the identity operator 1, we can freely use the fact that M = MM .
Take first the one-particle space. We have φ (g) = g + ∈ H 1 and this is in the above closure because i φ (g) = d dt e it φ (g) and is in the domain of φ( f ) [37, Theorem VIII.7]. By the one-particle Reeh-Schlieder property (e.g. [34,Theorem 3 Note that H 1 ∩Dom( φ (g)), which includes H 1 ∩Dom(χ (g)), has a dense subspace. For any such vector ξ ∈ H 1 ∩Dom(χ (g)), it holds i φ (g)ξ = d dt e it φ (g) ξ ∈ (C ⊕H 1 ⊕ H 2 ) and their projection to H 2 is dense in H 2 (as it comes from the action of z † ). We can subtract the (C ⊕H 1 )-component since M is an algebra, and obtain that H 2 -component belongs to MM = M . Therefore, it follows that The rest is shown by induction: assume that (C ⊕ · · · ⊕ H n ) ⊂ M . By differentiation, we obtain i φ (g) = d dt e it φ (g) ∈ (C ⊕ · · · ⊕ H n+1 ) and we can extract the H n+1 -component. Such H n+1 -components are obtained by z † , thus form a dense subspace in H n+1 . Namely, we showed that (C ⊕ · · · ⊕ H n+1 ) ⊂ M , which completes the induction.

(Non-)temperateness of the fields
Let us assume that we have a Haag-Kastler net. Borchers et al. [10] called an operator G which is affiliated to the wedge-algebra and generates a one-particle state from the vacuum a polarization-free generator. A polarization-free generator is said to be temperate if there is a translation-invariant dense domain of G such that, for any vector in that domain, a → GU (a, 0) is strongly continuous and polynomially bounded. The existence of temperate polarization-free generators restrict drastically the possibility of interaction.
We argue that the closure of our field φ( f ) (which we denote by the same symbol) has no such polynomially growing vector. Indeed, let be a vector in the domain of the closure of φ( f ). Even if is only in the domain of the closure of φ( f ) = φ( f ) + χ( f ), the operator φ( f ) is continuous on each space of fixed particle number H n , thus each n-particle component n of is in the domain of χ n ( f ), and φ( f ) can be computed componentwise. If we look at the one-particle component of χ( f )U (a, 0) , we get If a tends to the negative spacelike direction, this grows exponentially. As we remarked, the contribution from φ( f ) is bounded, hence this shows that φ( f )U (a, 0) is not polynomially bounded.

The Form Factor Program and Polarization-Free Generators
In this section, we discuss the connection of our approach to the form factor program [3,41]. We will also present how to derive our new term χ( f ) under certain assumptions. Both in the form factor program and in the operator-algebraic approach, we have to construct local observables. In the form factor program, the matrix elements of local operators are considered and they are subject to several conditions, the form factor axioms, and one has to find solutions (form factors) of these conditions. On the other hand, in the operator-algebraic approach, it has been proved [10] that, in any model, there are polarization-free generators, namely operators which are localized in a wedge-region and generate one-particle states from the vacuum. Therefore, if a Haag-Kastler net exists for an integrable QFT, there must be polarization-free generators. As they are localized in a wedge, they must commute with operators localized in the causal complement of the wedge, especially with local operators which come from form factors. Here we perform such formal computations.

Compatibility with the form factor program. Operator expansion for analytic S-matrix
In [12], we investigated the structure of strictly local observables in integrable quantum field theories in 1+1 dimensions in the case where the theory has only one species of massive scalar particle and the two-particle scattering function is analytic in the physical strip.
Specifically, we extracted more information on the properties of these local observables A by expanding them into a series of normal-ordered strings of z † , z, It was shown in [11] that this expansion holds for every operator or quadratic form A in a certain regularity class and it is independent of the localization region of A. It is similar to the well-known form factor expansion, although it is not identical to it; for a pointlike field A, the definition of the coefficients f [A] m,n (θ θ θ, η η η) formally agree with the form factors for certain regions of the arguments. More specifically, following [3, Section 3] [36, p.11] we have that (θ 1 , . . . , θ n ) for θ n > · · · > θ 1 and where F A denotes the form factor of A. For the other coefficients f [A] m,n (θ θ θ, η η η), there is an explicit expression in term of vacuum expectation values of A [11], and it was shown in [12] that if A is a general observable localized in a bounded region, then the f [A] m,n are boundary values of a common meromorphic function, i.e., one has Choosing the bounded region to be a standard double cone of radius r > 0 around the origin, these F k (ζ ζ ζ ) (ζ ζ ζ ∈ C k ) fulfill the following properties, which we write informally here (the exact statement of these properties can be found in [12]).
(F1) They are meromorphic on C k . In particular, except for possible first-order poles at ζ n − ζ m = iπ ("kinematic poles"), they are analytic on the tube region Im ζ 1 < · · · < Im ζ k < Im ζ 1 + 2π . (F2) They are S-symmetric, that is, for all complex arguments ζ 1 , . . . , ζ k , the following relation holds: (F3) They are S-periodic, i.e., (a similar property holds also in the other variables.) (F4) The value of their residue at ζ k − ζ 1 = iπ is given by (and a similar formula holds for the residues at the other kinematic poles.) The remaining properties concern the bounds that these functions fulfill, namely (F5) θ θ θ → F k (θ θ θ + i0 0 0) is square integrable.
(F6) At real infinity the growth behaviour of F k is essentially given by These relations between form factors and local observables can be made precise and suitable variations of (F1)-(F6) holds for F k if and only if A is localized in a double cone (see [12,Theorem 5.4]).
The functions F k have a rich pole structure: Apart from the so-called kinematic poles, they have further singularities on hyperplanes in C k , which are due to the poles of the scattering function "outside the physical strip".

Cases with poles in the physical strip
In the present paper, the scattering function has additional poles in the physical strip. We will therefore expect that properties (F1)-(F6) need to be modified to take into account these extra poles of the S-matrix. We call these new properties (P1)-(P6). In particular, for (P1) we demand that the functions F k are meromorphic with the pole structure as in (F1) and additional poles at ζ n − ζ m = 2πi 3 with n > m. The presence of these additional poles, which are a consequence of the poles of the S-matrix inside the physical strip, is a well-known feature in the form factor program [4].
Further, we expect that properties (P2) and (P3) can stand unmodified from (F2) and (F3), while (P5), which is required for the well-definedness of the expansion (4), will possibly need small modifications from (F5); these depend on details of the domain of definition of (4) in a theory with bound states, which we do not enter here.
We will also need to properly adapt (F6) to account for the extra poles of F k (ζ ζ ζ ) at ζ n − ζ m = 2πi 3 with n > m, yielding a new condition (P6). However the essential new ingredient to be added to [12,Definition 5.3] is described as follows: (P7) The F k have first order poles at ζ n − ζ m = 2πi 3 , where 1 ≤ m < n ≤ k and one has res the residues at the other poles can again be inferred from (F2).
A similar formula can be found in the form factor program [4] up to a constant factor, where it characterizes the form factors of point-like localized fields in models with bound states.

Sample computations
We would like to show that if the form factors We compute the commutator , [A, χ (g)] = , Aχ (g) − , χ (g)A using the actions of (4) and of χ (g) derived in Sect. 3.2 on the vectors , . We obtain , Aχ (g) noting that in the expansion (4) the coefficients which contribute to the matrix elements , Aχ (g) and , χ (g)A are F 0 and F 2 . We observe that the -prescription in the argument of F 2 in (5a) was introduced in [12], even if F 2 (ζ ζ ζ ) is actually analytic at ζ 2 − ζ 1 = iπ . Further, the vector (A ) 1 (ζ 1 ) = dζ 2 F 2 (ζ 1 , ζ 2 ) (ζ 2 ) in (6a) is no longer analytic at ζ 1 = πi 3 due to the pole of the integrand at ζ 2 − ζ 1 = 2πi 3 , so that χ as given in Sect. 3.2 cannot be applied to (A ) 1 (ζ 1 ). However, by the symmetry of χ (g), we can apply it to the vector and then shifting the integral contour, we obtain the boundary value indicated in (6a).
The terms (5b) and (6b) depending on F 0 cancel each other. We shift the integral contour in (5a) in the variable ξ by − iπ 3 + i0 (taking care to not cross the pole hyperplane of F 2 (ζ ζ ζ ) at Im(ζ 2 − ζ 1 ) = 2π 3 ). For this the integrand is required to decay at infinity which should follow from the growth properties of g ± and from (P6).
Hence, we find We consider now the commutator , , where we inserted the expression of φ (g).
Using the action of (4) and of z , z † on the one-particle vectors , , and noting that to these matrix elements the coefficients in (4) which contribute are F 1 and F 3 , one computes where we used that (J g − )(ξ ) = g − (ξ ); for the direction of the boundary values in F 3 and F 1 see [12,Definition 4.3].
By shifting the integral in (8c) in the variable ξ from R to R + iπ , we find (8a) up to residues (again bounds at infinity enter). Specifically, the analytically continued function 3 . Using (P7), the residue at the first pole is given by This residue gives a contribution to the difference of the two integrals (8c) and (8a), namely (8a) + (8c) where the orientation of the residue is fixed using [12,Lemma 3.4].
In computing the second residue + res we can assume θ = ρ since it suffices to compute a meromorphic function on any open set. Note that the two S-factors are analytic near the pole hyperplane ρ+iπ −i −ξ = 2πi 3 and only F 3 has a pole there. In the second equality we used (P7) and the Bootstrap equation (S4). The sign of the residue is fixed by using [12,Lemma 3.4]. This residue gives a contribution to the difference of the integrals in (8c) and (8a), namely Summarizing our results, we have where we used that g + (θ ) = g − (θ ± iπ). In (7a) the coefficient F 2 (ζ ζ ζ ) approaches the pole in the direction ζ 2 − ζ 1 = 2πi 3 + i , instead in (10b) it approaches the pole in the direction ζ 2 − ζ 1 = 2πi 3 − i . So the difference of these boundary values is again a residue that we can compute using (P7) and [12,Lemma 3.4], This residue gives a contribution to the difference of (7a) and (10b), namely Likewise, the difference of the boundary values in the integrals (7b) and (10a) gives the residue We now consider (10c) and (10d). When shifting the integral contour in (10d) in the variable ξ from R to R + iπ , F 1 is analytic, but the factor S(ζ ) has poles at ζ = 2πi 3 , πi 3 , the residues of which we denote by R, R , respectively. Using [12,Equation (3.14)], the difference of the boundary values in the integrals (10c) and (10d) gives the following residue Comparing this with (12) and (11), we find −(13) = (12) + (11), where we used that R = −R, η = i √ 2π |R|, R = i|R| and g + (ρ) = g − (ρ ± iπ). Hence, the matrix element , [A, φ (g)] vanishes.

A guesswork for the bound-state operator.
Our operator χ( f ) appears to be new and it is worth presenting how one can "find" it by formal arguments. Again in this subsection we ignore all subtleties with domains and assume that we can compute everything termwise. Borchers et al. proved in [10] that, for any Haag-Kastler net with mass eigenvalues, one can construct polarization-free generators, namely (unbounded) operators which create one-particle states from the vacuum and are affiliated to the von Neumann algebra associated to wedges. In particular, if we could construct a Haag-Kastler net for a given scattering function S, then there would have to be polarization-free generators. On the other hand, the expansion (4) in Sect. 4.1 should exist, possibly in a modified form. Let us assume that B is a polarization free generator. As B generates only a one-particle state from the vacuum, by a (termwise) straightforward computation, one realizes that f [B] m,0 must vanish except m = 0, 1. We may assume that B is symmetric, then also f [B] 0,n survive only if n = 0, 1.
Let h be a real test function such that [27] has only (m, n) = (0, 1), (1, 0) components. As we know that this cannot be wedge-local, let us consider the simplest variation of it. The (0, 0) component is a scalar and has no effect on the commutation relations. Assume B has a (1, 1) component, namely that B = z † (h + ) + X + z(J 1 h − ), that X preserves the one-particle space and X = X * = 0.
Suppose that we have the Haag-Kastler net. If B is affiliated to A(W L − (0, r )), then it must commute with operators A which arise from form factors as in Sect. 4.1. Especially, for a one-particle vector , the commutator , [A, B] must vanish. We may assume that both A and B are symmetric, and thus h + = J 1 h − . The vector is annihilated by z(J 1 h − ) and X , therefore, the only remaining terms are where in the 3rd equality we used S-periodicity of F 2 , and the relations between form factors of A and J A * J in [16,Proposition 8.9] (which is also used in the 6th equality). In the 4th equality, we used the residue formula (P7) and then in the 5th equality we shifted the integral contour by − πi 3 , where we assumed that does not have poles (again, this is a guesswork). Now, if we assume the Reeh-Schlieder property for the net A, there must be sufficiently many local observables A, in particular, there are many such A's that H 1 is spanned by the one-particle components of A 's. As we require that the above commutator should vanish, it must hold that X = χ(h) . Namely, we correctly guessed the one-particle action of χ(h). We know that this choice works, as we saw in Sect. 3.3.
Note that χ(h) cannot be written exactly in the form of (4), due to its subtle domain property. Recall [16] that the expansion (4) is available only for certain regular operators. We were fortunate that this guesswork resulted in a correct answer, but in more complicated models it fails [17]. On the other hand, if the Haag-Kaslter net for the model existed at all, there would have to be polarization-free generators [10]. This failure might be a consequence of the complicated domains, or one would have to add higher component in the expansion (4). After all, we are not excluding wedge-local fields which are more complicated, yet have better domain properties.

Conclusions and Outlook
In the present work, our aim was to extend Lechner's construction of two-dimensional models of quantum field theory with a factorizing scattering matrix to models with bound states, which are associated with poles of the S-matrix in the physical strip. We considered scalar S-matrices with only one pair of poles in the physical strip, corresponding to two bosons of the same species which fuse into another boson of the same species.
The properties of the new class of S-matrices include the presence of these extra poles in the physical strip (and the values of the corresponding residues of the S-matrix at these points) and the Bootstrap equation, in addition to the properties that can be found in [28]. Quantum integrable models with underlying S-matrix fulfilling this new set of properties include the Bullough-Dodd model.
For this class of S-matrices we constructed wedge-local fields φ( f ) by adding the bound-state operator χ to the wedge-local field of Lechner. However, with the simple domain taken here, the field φ( f ) is expected not to be self-adjoint and we have not shown strong commutativity. Rather, we have shown that the field φ( f ) weakly commute with its reflected field φ (g) on a common domain.
Furthermore, we gave a partial answer to the question of compatibility of our construction with the form factor program. In particular, we conjectured that a suitable variation of the characterization theorem for local observables in [12] holds in our class of models: we argued explicitly that this is true at least in one-particle states.

Open problems
One of the major open problems is to show the strong commutativity of the fields φ( f ) and φ (g) in order to obtain a Borchers triple. This would mean not only proving the existence of self-adjoint extensions of the two fields, but also selecting extensions that strongly commute. This is a highly non-trivial task for which the results available in the literature can offer only partial answers. Some results on the construction of extensions have already been obtained by one of the authors [45].
Following Buchholz and Lechner [14], the natural next step would then be to show the Bisognano-Wichmann property, i.e., the geometric action of the modular group. Thereafter, one would like to establish modular nuclearity condition [29] in order to prove that the wedge algebras so generated are split, and therefore the non-triviality of the double-cone algebras.
Another interesting problem would be an extension of our construction, valid for the moment only for scalar S-matrices, to a larger class of integrable models with bound states, such as the Z (N )-Ising model [5] and the sine-Gordon model [46] for a certain range of the coupling constant. Indeed, in the Z (3) model one can find a multi-component field which fulfills weak wedge-commutativity at least on the level of formal computation [17]. In the Z (N ) (with N > 3) and sine-Gordon models, this will hold only for certain components of the fields. The construction involves a multi-component generalization of the bound-state operator χ , which will now modify the multi-component field of Lechner-Schützenhofer [33]. The matrix-valued analogue of the properties of S introduced in Sect. 3 will need to include the Yang-Baxter equation, which will play an important role in the proof of wedge-commutativity.
Let us look at the bootstrap equation (S4), where the boundary value is represented . The right-hand side is meromorphic, therefore, S must have a meromorphic continuation in a neighborhood of R as well. Yet, by (S1), the boundary value has modulus 1, hence it is actually continuous in a neighborhood of R. In general, an inner function can have essential singularities at the boundary. It follows from the observations above that the only possible essential singularities in the picture of R + i(0, π) is θ = ±∞, as in [35,Appendix A]. By the uniqueness of the factorization, the properties (S2)-(S3) must be satisfied by the constant, the singular and the Blaschke factors separately. Therefore, c = ±1. Consider next the singular part S ∞ . It admits an integral representation as [38,Theorem 17.15] and it has essential singularities at points in the support of the measure. By considering the symmetry (S2) and the possible essential singularities as above, it has the form S ∞ (θ ) = e ia(e θ −e −θ ) , a ≥ 0. By a straightforward computation, this satisfies (S4) and S ∞ ( πi 3 ) = S ∞ ( 2πi 3 ) > 0. Let us then turn to the Blaschke factor S Blaschke . It can be written as the following infinite product: where α n ∈ R + i(0, π) and c n = − |β n | β n 1−β n 1−β n , is a constant, in accordance with [38,Theorem 15.21], where β n = e αn −i e αn +i (and c n = 1 if α n = πi 2 (and hence β n = 0) by convention). These α n 's are exactly the zeros of S Blaschke . By (S2), α n is purely imaginary or it must appear in a pair with −α n (including multiplicity). Next, by (S3), Im α n = π 2 or it must appear in a pair with iπ − α n . In each case, it is straightforward to see that the constant factors c n cancel each other and S Blaschke (ζ ) = n e ζ −e αn e ζ −e αn . As S and S ∞ satisfy the bootstrap equation (S4), so does c · S Blaschke . Therefore, if there is any α n such that Im α n < π 3 , then 0 < Im α n + π 3 < π and α n + πi 3 must be another zero. If π 3 < Im α m , then by combining (S3) and (S4), we may assume that there is an α n such that Im α n < π 3 . We summarize these observations. Let α n with Re α n = 0 appear in the product. Then, if Im α n = π 3 (if Im α m = 2π 3 , by (S3) we may assume that there is α n such that Im α n = π 3 ), then the factor (e ζ − e α n )(e ζ − e −α n )(e ζ − e πi−α n )(e ζ − e πi+α n ) (e ζ − e α n )(e ζ − e −α n )(e ζ − e πi−α n )(e ζ − e πi+α n ) must appear. If Im α n = π 3 , then we may assume that Im α n < π 3 and (e ζ − e α n )(e ζ − e −α n )(e ζ − e πi−α n )(e ζ − e πi+α n ) (e ζ − e α n )(e ζ − e −α n )(e ζ − e πi−α n )(e ζ − e πi+α n ) × (e ζ − e α n + πi 3 )(e ζ − e −α n + πi 3 )(e ζ − e πi−(α n + πi 3 ) )(e ζ − e πi+α n + πi 3 ) (e ζ − e α n − πi 3 )(e ζ − e −(α n + πi 3 ) )(e ζ − e πi−α n + πi 3 )(e ζ − e πi+α n + πi 3 ) must appear. Both of these factors satisfy (S1)-(S3), and also (S4) by straightforward computations. Furthermore, these factors take positive numbers if ζ is purely imaginary: each of these products satisfies (S2), hence it takes non-zero real value on the imaginary axis, therefore, it has a fixed sign, but at ζ = 0 it can be directly checked that it takes 1, hence must take positive values on the imaginary axis. Finally, if S Blaschke has a zero at α n such that Re α n = 0, then the factors f B k (ζ ) of Sect. 3.1 must appear. These factors satisfies (S1)-(S4), and f B k 3 ) < 0. Now, as the Blaschke factor is determined by zeros, S Blaschke satisfies (S4) by itself, which implies that c also must satisfy (S4), namely c = 1. In order to assure (S5), S Blaschke can have an odd number of factors as f B k (ζ ) (infinite is impossible, because S(ζ ) would have an essential singularity at ζ = 0, which we have excluded), so that the residues of the factor f 2 3 gets multiplied by a negative number in the whole two-particle S-matrix S(ζ ) = f 2 3 (ζ ) · S ∞ (ζ )S Blaschke (ζ ) (recall that f 2 3 has only simple poles). Now it is obvious that S Blaschke (0) = 1.

B. Lemmas on Shifting Integral Contour
Here we state some easy facts explicitly for the sake of clarity.
Lemma B.1. Let A 1 and A 2 be positive self-adjoint operators which strongly commute. If , ∈ Dom(A 1 ) ∩ Dom(A 2 ), then for any 0 ≤ ≤ 1, we have that , ∈ Dom(A 1− 1 A 2 ) and it holds that This might seem obvious, and indeed the proof is not too complicated, yet one should note that A 1 A 2 etc. are ill-defined.
Proof. Consider the joint spectral decomposition with respect to A 1 and A 2 . Any vector is a (possibly vector-valued) L 2 -function on this spectral space of two variables a 1 , a 2 , such that A 1 and A 2 act by multiplication with e a 1 , e a 2 , respectively. Now, the hypothesis that ∈ Dom(A 1 ) means that (a 1 , a 2 ) and e a 1 (a 1 , a 2 ) are both L 2 . This implies that e 2(1− )a 1 | (a 1 , a 2 )| 2(1− ) is L 1 1− . Similarly, if ∈ Dom(A 2 ), then e a 2 (a 1 , a 2 ) is L 2 and hence, e 2 a 2 |ξ(a 1 , a 2 )| 2 is L 1 . The first claim of the present lemma is that ∈ Dom(A 1− 1 A 2 ), which is equivalent to e 2(1− )a 1 +2 a 2 | (a 1 , a 2 )| 2 being L 1 . This is a consequence of Hölder's inequality and of the estimates above.
The desired equality is obtained simply by writing it as the spectral integral: a 1 , a 2 )e a 2 (a 1 , a 2 ) = da 1 da 2 e a 1 +(1− )a 2 (a 1 , a 2 )e (1− )a 1 + a 2 (a 1 , a 2 ) Proof. and can be considered as vectors in L 2 (R, dθ) ⊗ K and in the domain of the operator A ⊗ 1, where (A it ξ )(θ ) = ξ(θ + t ). One can either apply Lemma B.1, or essentially repeat the proof of Proposition 3.1 by scaling the variable and putting a constant instead of f + .

C. Comments on the Domains
We have seen that, if f is a real test function supported in W L , then χ( f ) is a symmetric operator. Then the question arises whether χ( f ) is (essentially) self-adjoint and, if not, what the appropriate extension is.

There are many self-adjoint extensions
We claim that the answer to the first question is in general no. Indeed, if one considers a similar operator which is a product of a multiplication operator by an analytic function and an imaginary shift, then its properties depend very much on the analytic function. More precisely, the deficiency indices of the product operator are tightly related to the zeros of the analytic function and its decay rate as |Re ζ | → ∞ [45]. In our case, as f + is strongly decreasing, we expect that χ 1 ( f ) = x f 1 12 y h · y * h 1 12 , which are self-adjoint due to the theorem of von extensions, which does not appear in the cases with analytic S-matrices, is rooted in the essential properties of bound states, therefore worth a serious investigation.
We plan to systematically study this issue in [45].