Asymptotic Completeness in a Class of Massive Wedge-Local Quantum Field Theories in any Dimension

A recently developed n\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n$$\end{document}-particle scattering theory for wedge-local quantum field theories is applied to a class of models described and constructed by Grosse, Lechner, Buchholz, and Summers. In the BLS-deformation setting we establish explicit expressions for n\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n$$\end{document}-particle wave operators and the S\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S$$\end{document}-matrix of ordered asymptotic states, and we show that ordered asymptotic completeness is stable under the general BLS-deformation construction. In particular, the (ordered) Grosse–Lechner S\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S$$\end{document}-matrices are non-trivial also beyond two-particle scattering and factorize into 2-particle scattering processes, which is an unusual feature in space-time dimension d>1+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d > 1+1$$\end{document}. Most notably, the Grosse–Lechner models provide the first examples of relativistic (wedge-local) QFT in space-time dimension d>1+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d > 1+1$$\end{document} which are interacting and asymptotically complete.


Introduction
Asymptotic completeness (AC) of a quantum field theory (QFT) is the requirement that incoming and outgoing multi-particle scattering states Ψ ± should be dense in the Hilbert space H of a given QFT model. This conceptual question was already raised early on during the development of the quantum theory of fields [Gre61;Ru62]. But even after many decades of research, mathematical results on AC of interacting QFT models have remained rather scarce (e.g. [CD82;DG99;IM06;Le08]) and there are various obstacles of physical and technical nature (cf. [DyG13]). In fact, the first complete proof of asymptotic completeness of an interacting QFT has been established rather recently: It was obtained by Lechner [Le06;Le08] for certain two-dimensional models with factorizing S-matrix, such as the Sinh-Gordon model.
In the present work we provide, to the best of our knowledge, the first full proof of AC in a class of interacting wedge-local QFT models on Minkowski space-times of dimension d > 1 + 1. The models in question, found by Grosse, Lechner, Buchholz and Summers in [GL07; BS08; BLS11], may not be local, but have the weaker property of wedge locality. 1 The constructive procedure developed in [BLS11], which we will call a BLS-deformation, preserves the general structure of a relativistic wedge-local theory and in addition introduces interaction. For example, the BLS-deformation of a free field theory, which we will call a GL-model, has a non-trivial twoparticle scattering matrix [GL07;BS08]. However, collision processes involving n ≥ 3 particles were not investigated in these works, since n-particle scattering theory was not available in the wedge-local setting back then. Such a scattering theory has been developed meanwhile in [Du18], so that the question of asymptotic completeness can now be posed and positively answered for these models. More than that, by adapting and generalizing methods from [BLS11; DT11] we will establish explicit expressions for the general effect of BLS-deformations on asymptotic states and scattering data. In this manner we also prove stability of asymptotic completeness under BLS-deformations for general wedge-local models with massive particles. To construct n-particle scattering states according to [Du18], we define Haag-Ruelle operators B kτ (f k ), τ ∈ R, 1 ≤ k ≤ n, which create the desired one-particle states from the vacuum Ω ∈ H . Then outgoing and incoming velocity-ordered scattering states are (1) The momentum-space configuration of such scattering states is specified by a family of regular Klein-Gordon solutions f k . For wedge-ordered velocity configurations (as defined below in (2)), the wedge-local Haag-Ruelle theorem establishes convergence in (1) and Fock structure of the resulting asymptotic states [Du18]. This result holds in general wedge-local theories with massive particle spectrum, given that all B kτ (f k ) are constructed from a common localization wedge W.
For the concrete case of the standard wedge W = W R = {(t, x) ∈ R d : |t| < x 1 }, the velocity ordering of outgoing states (τ → +∞) reads Here V f k := {(1, v m (p)) ∈ R d : p ∈ suppf k } denote the supports of the Klein-Gordon wave packets with respect to velocity v m (p) := p/ω m (p). Further, the precursor ordering relation in (2) is defined by V f k ≺ WR V f k−1 :⇐⇒ V f k−1 − V f k ⊂ W R and corresponds to the geometrical requirement that all relative velocities v k−1 −v k ∈ V f k−1 −V f k yield a positive directed separation of wave packets f k for t → +∞ with respect to the space-like opening direction of W R , i.e.
(v k−1 − v k ) 1 > 0 for all 2 ≤ k ≤ n. In the present work we prove, in particular, for GL-models that the velocity-ordered scattering states constructed according to (1) and (2) for each wedge W span dense sets H ± W in the full Hilbert space H of the interacting model. In this case, we will say that a wQFT model has the property of (ordered) asymptotic completeness.
We now briefly explain the basic ideas of the present analysis of deformed scattering data and asymptotic completeness in non-technical terms. In the deformation construction from [BS08; BLS11] the new deformed model is generated by observables constructed as formal spectral integrals In the terminology of [BS08;BLS11], the (smooth) wedge-local operators A of the initial model are 'warped' with respect to the spectral measure dE = dE (H,P ) of the energy-momentum operators P = (H, P ). Presently, Q ∈ R d 2 denotes a fixed parameter matrix satisfying certain geometrical properties [GL07]. Later it will be chosen as a mapping depending on the localization wedge of the operator A, as prescribed by Grosse and Lechner, to yield wedge-locality and Poincaré covariance (if applicable) of the deformed model. Proceeding now on the formal level of (4), let us apply the wedge-local n-particle scattering theory from [Du18] and compare the scattering states where B kQτ (f k ), 1 ≤ k ≤ n, obtained from (4) are a corresponding family of creation operators associated to the deformed model. Inserting the definition of the warped convolutions appearing in (6) we see directly that On a heuristic level we can thus write Ψ + Q as a superposition of scattering states constructed in the undeformed model Ψ + 0;p1,...,pn := lim τ →∞ α Qp1 (B 1τ (f 1 ))α Qp2 (B 2τ (f 2 )) . . . α Qpn (B nτ (f n ))Ω which can be related to scattering states of the usual form (5), with operators modified by space-time translations depending on p 1 , . . . , p n ∈ R d . As the action of space-time translations on the one-particle states can be made explicit, namely α y (B τ (f ))Ω = B τ (f y )Ω, where f y (x) = f (x − y), we may expect on heuristic grounds that, also for a general deformed wedge-local model, the corresponding subspaces of scattering states of the deformed and undeformed model will coincide, that is, H + Q,≺W = H + 0,≺W ⊂ H .
Here H denotes the full Hilbert space of the interacting model, and similarly H − Q,≺W = H − 0,≺W . This shows that asymptotic completeness of ordered scattering states is stable under BLSdeformations. Asymptotic completeness of the Grosse-Lechner models follows from AC of the free massive scalar field theory with respect to ordered scattering states. Indeed, using that for local QFTs the space of wedge-ordered scattering states H + 0,≺W coincides with the full space of scattering states H + 0 , and exploiting AC of the free field, H + 0 = H 0 , we obtain for the GL-model Hence these wedge-local QFT models are asymptotically complete.
It is the main contribution of the present paper to make the above heuristic arguments leading from (7) to (9) mathematically rigorous. Our analysis is based on oscillatory integral techniques from [BLS11] and the recent wedge-local n-particle scattering theory [Du18]. We also profit from the previous analysis of a scattering in BLS-deformed massless two-dimensional models by Dybalski and Tanimoto [DT11]. This case is simplified by the absence of dispersion, which reduces the problem to two-body scattering of left-and right-movers. Hence wedge-local scattering theory with n ≤ 2 particles suffices. This has in fact been developed and studied much earlier [BBS01; Le03; GL07; BS08] via direct application of standard Haag-Ruelle methods. In particular, the swapping strategy introduced in [Du18] to establish scattering theory for arbitrary particle numbers is not needed for this specific class of 1+1-dimensional models.
To conclude this introduction we should remark that physical intuition familiar from local QFT makes (10) very plausible on one hand. On the other, such intuition may fail in the much larger class of general wedge-local models. Interesting examples can already be found among recently constructed free product models [LTU17,Sec. 5]. There it was shown that ordered twoparticle scattering states are not sufficient for two-particle asymptotic completeness. In contrast to a previous work on clustering properties [Sol14], where scattering theory is also mentioned briefly, our results do not assume locality of the underlying wQFT model and in this regard the results of the present work are much more widely applicable within the general wedge-local framework.
2 Preliminaries on Wedge-local QFT (wQFT) and Scattering Theory 2.1 Operator-algebraic Framework for Wedge-local QFT We will work in an operator-algebraic setting of wedge-local quantum field theory on Minkowski space-time R d , and our results are valid in arbitrary spatial dimension s := d − 1. The family of wedge regions is defined as the orbit under the action of the Poincaré group P = R d ⋊ L.
A wedge-local quantum field theory model is specified by mathematical objects (A, α, H , Ω), where H is the Hilbert space of pure states containing the vacuum as a distinguished unit vector Ω ∈ H . The wedge-localization of observables is described by a family of von Neumann algebras A(W) ⊂ B(H ) associated to wedge regions W. Poincaré symmetry acts on the wedgelocal algebras A(W) by a given group of isomorphisms α λ and we denote by λ = (x, Λ) ∈ P ↑ + = R d ⋊ L ↑ + the elements of the proper orthochronous Poincaré group. In this paper we are working mostly with space-time translations of some given operator A ∈ A(W) by x ∈ R d , also denoted by α x (A).
Guided by physical intuition one asks that these objects satisfy wedge-local variants of the Haag-Kastler postulates, which are concerned with the algebraic and representation-theoretic properties of A. Firstly, for any choice of wedge regions W, W 1 , W 2 one has Here the Minkowski causal complement W ′ = (ΛW R + x) ′ = −ΛW R + x of W is also a wedge region and A(W) ′ denotes the commutant of A(W) relative to B(H ).
On the representation-theoretic side it is further assumed that translations are unitarily implemented on the vacuum Hilbert space H by a strongly continuous s + 1-parameter group, α x (A) = U (x)AU (x) * . The representing unitaries are generated by the energy-momentum operators via U (x) = U (t, x) = e itH−ix·P , whose joint spectral resolution in terms of projectionoperator-valued measures will be denoted ∆ −→ E (H,P ) (∆) and abbreviated as E(∆) for Borel sets ∆ ⊂ R d . Focusing on the analysis of scattering, we impose the following standard assumptions concerned with the vacuum representation and its one-particle spectrum, This summarizes the axiomatic operator-algebraic formalism for wedge-local QFTs providing the basis of our present investigations. Our choice of framework serves the purpose of accommodating scattering-theoretic reasoning and capturing the requirements of the wedge-local Haag-Ruelle theory [Du18]. In the literature also the closely related concept of a causal Borchers triple (R, U, Ω), corresponding to R := A(W R ) in the present framework, is often studied. Such a framework is equivalent to the above setting including Poincaré covariance (HK3 ♯ ). For further historical background and other aspects of wedge-local Quantum Field Theory, we refer to [BS08; Le15; BLS11; Du18].
Regarding the scattering-theoretic analysis, Poincaré covariance (HK3 ♯ ) is not essential. The important properties for us are the existence of isolated mass shells (HK6 ♯ ) and the wellestablished wedge duality condition (HK2 ♯ ), which strengthen (HK6) and (HK2), respectively. As such, these will be standing assumptions of this work. Whereas (HK6) merely demands positivity of the Hamiltonian in any Lorentz frame, the Haag-Ruelle mass gap condition (HK6 ♯ ) physically amounts to non-triviality of the one-particle subspace H 1 := E(H m )H . The upper mass gap M > m has the same technical purpose as in traditional Haag-Ruelle scattering theory: it enables the efficient separation of one-particle states from the remaining energy-momentum spectrum.

Warped Convolutions and the Grosse-Lechner model
Our scattering-theoretic analysis is concerned with a general class of massive wedge-local QFT models, which are constructed via a certain deformation method introduced by Buchholz, Lechner, and Summers [BLS11]. The Grosse-Lechner models were introduced and studied from a non-commutative geometry perspective in [GL07]. Here we adopt the operator-algebraic approach from [BS08; BLS11] and study, in particular, the GL-models as constructed by applying BLS-deformations to the free scalar field.
The starting point of the BLS-construction is a general wedge-local model (A 0 , α, H , Ω). With κ ≥ 0 and, in dimension d = 3 + 1 additionally η ∈ R, as parameters of the deformation, a family W −→ Q W of warping matrices is now defined according to [BS08; GL07] by either for the warping matrix of the reference wedge W R = {x ∈ R d : x 0 < x 1 } and its translates. For general wedges given by W = ΛW R + x the warping matrix is defined by Poincaré-covariance 2 as Q W := ΛQ WR Λ −1 . These warping matrices are antisymmetric with respect to the scalar product defined by the Lorentzian metric g with signature (+, −, . . . , −).
Let now Q ∈ R d×d . The warped convolution of an operator A is denoted by and usually intuitively but formally written as a spectral integral with operator-valued integrand. Buchholz, Lechner and Summers [BLS11] have given a rigorous definition by interpreting (12) as an integral of oscillatory-type. To this end, one starts with operators A on which space-time translations act smoothly and defines A Q in terms of its action on vectors from the dense domain of finite energy states. We denote by C ∞ the algebra of all regular operators A ∈ B(H ), defined by the requirement that R d ∋ x −→ α x (A) = U (x)AU (x) * is arbitrarily often differentiable with respect to the operator norm topology. By standard mollification arguments, the regular subalgebra A 0r (W) := C ∞ ∩ A 0 (W), for any wedge W, is weakly dense in the full wedge algebra A 0 (W) of the initial model.
Definition 2 (warped convolution of a regular operator [BLS11]). The warped convolution A Q of a bounded regular operator A ∈ C ∞ with respect to a warping matrix Q is defined by its action on finite-energy vectors Ψ ∈ D as limit of strong integrals of oscillatory type, with η ∈ S (R d × R d ), such that η(0, 0) = 1, serving as a regularizing function.
Here the technical assumptions of regularity of A and the consideration of only finite energy states are needed for the existence of the limit ǫ → 0. In addition to existence of this limit, it is shown in [BLS11] that A Q extends to a bounded operator, and is independent of the choice of f . Let us summarize some helpful properties of warped convolutions from the literature. For standard results about Bochner integrals we refer to [Zaa67, Ch. 6 §31].

Lemma 3. [BLS11]
The warped convolution A Q of any regular operator A ∈ C ∞ , given by (14), is well-defined. In particular it does not depend on the specific choice of the regularizing function f . The deformed operator extends to a bounded operator A Q ∈ B(H ), which satisfies (v) Warping commutes with space-time translations, (vi) If the underlying model satisfies Poincaré-covariance (HK3 ♯ ) one has further It is important to emphasize that the definition (14) of A Q is only mathematically rigorous when working with regular operators A ∈ A 0r . To obtain the full warped wedge-local QFT model in the sense of Section 2.1, we additionally pass to the weak closure of the set of warped regular operators.
Definition 4 (warped wedge-local model [BS08;BLS11]). Let Q W be given by (11) for W = W R and let Q W := ΛQ WR Λ −1 for general wedges W = ΛW R + a. The warped wedge-local algebras are defined as where the bicommutant is taken with respect to B(H ). We further note that [BLS11] study wedge-local quantum field theories in the framework of a (causal) Borchers triple (A(W R ), α, Ω). There one specifies only an observable algebra for the single wedge region W R . The full family of wedge algebras is then obtained using Poincaré symmetry (HK3 ♯ ). This symmetry of the wQFT model is of course an additional requirement, which is natural from the perspective of physics and a useful constraint for constructive efforts. In particular, it is a non-trivial feature that BLS-deformations can be used to construct Poincaré covariant models. However, the BLS-deformation method also applies to the present framework, where Poincaré covariance is optional. The proofs of Lemma 3 and Theorem 5 can be extracted without significant modifications from the results of [BLS11].

Scattering States in Wedge-local QFT
We briefly review the n-particle Haag-Ruelle scattering theory for wedge-local QFTs [Du18]. To construct one-particle states via the wedge-local Haag-Ruelle method, one chooses operators A ∈ A(W) such that E(H m )AΩ = 0. The existence of such operators follows from (HK5). Let χ ∈ S (R d ) be a Haag-Ruelle auxiliary function supported within a sufficiently small neighborhood of the mass shell H m , disjointly from the remaining energy-momentum spectrum. Then the one-particle problem is solved by Namely, in terms of the Fourier transform defined in the relativistic unitary convention bŷ we have BΩ = (2π) d/2χ (H, P )AΩ by spectral calculus and translation-invariance of the vacuum. Hence BΩ is in the one-particle space H 1 := E(H m )H , as a consequence of the support of the Haag-Ruelle auxiliary functionχ intersecting the energy-momentum spectrum supp E (H,P ) of the theory only on subsets of the mass shell H m (by construction and using also the mass gap assumption (HK6 ♯ ), see e.g. [A,Ch. 5] or [Dyb17, Sec. 2.1]).
Proceeding towards the n-particle problem we define similarly for a family of A j ∈ A(W), 1 ≤ j ≤ n, the operators B j := A j (χ). Further we consider (positive-energy) Klein-Gordon solutions as scattering-theoretic comparison dynamics, with relativistic dispersion ω m (k) := √ k 2 + m 2 for mass m > 0, and recall that we abbreviate s := d − 1. For usual technical reasons, the wave packets are assumed to be regular, that is,f j ∈ C ∞ c (R s ). By a stationary phase analysis, these regular Klein-Gordon solutions vanish rapidly in all space-and time-like directions away from the classical propagation cones We note that these cones describe the scattering geometry of the single-particle wave packets.
The construction of n-particle states is now accomplished by means of wedge-frame adapted Haag-Ruelle creation-operator approximants Here, the boost Λ ∈ L ↑ + specifies an auxiliary Lorentz frame used for the construction. The simplest case is A j ∈ A(W) for W = ±W R in two dimensions, and in higher dimensions also for spatial rotations of W R . Then one can simply take Λ = 1, corresponding to equal-time hyperplanes. For general wedges W it is technically preferable to choose this boost from Regarding the causal geometry of wedges, we recall here that in dimension d ≥ 2 + 1 a general wedge region W can be written as W = λW R := ΛW R +x in terms of the standard Rindler wedge W R and some Poincaré transformation λ = (x, Λ) with Lorentz transformation part Λ ∈ L ↑ + and translation x ∈ R d . Dimension d = 1 + 1 is special in that, after restricting to the proper orthochronous Poincaré group, the set of all wedges splits into the disjoint orbits of the right wedge W R and left wedge W L : The choice of Λ enters in the space-time localization of the operators B Λ j,τ (f j ), and keeping track of the latter is important for the Haag-Ruelle method. Namely, the localization depends on τ and the wedge W of localization of A j , translated to Here x ∈ R s } is the standard space-like hyperplane at τ = 1 and V Λ fj is the velocity support of f j with respect to the Lorentz frame specified by Λ. For further details we refer to [Du18].
Geometrically, scattering situations are concerned with phenomena at very large distances, and it is convenient to introduce the centering of W denoted by W c := ΛW R . As centered wedges W c are convex cones with W c + W c ⊂ W c , the precursor relation defined for non-empty regions O 1 , O 2 ⊂ R d , is transitive, and Poincaré covariant as a partial ordering in the sense that O 1 ≺ W O 2 ⇐⇒ λO 1 ≺ ΛW λO 2 , λ = (x, Λ) ∈ P ↑ + . To construct n-particle scattering states we consider ordered configurations of wave packet velocities for the outgoing limit τ → +∞, and the reversed ordering for the incoming scattering limit τ → −∞, respectively. In these cases the main result of [Du18] establishes convergence of the scattering-state approximants and Fock structure of the limits, where all underlying wedge-local operators A j ∈ A(W) are localizable in a common wedge W. This approach is distinct from conventional Haag-Ruelle theory, where convergence and Fock structure proofs use that the commutators vanish rapidly in norm for j = k when τ → ±∞. Such stronger estimates are obtained from locality of the QFT and disjoint velocity supports of the Klein-Gordon solutions. Applying this reasoning in a wedge-local context leads to a more restrictive setup. Firstly we should take Klein-Gordon solutions f, f ⊥ must be chosen such that the velocity supports are not merely disjoint, but satisfy the stronger geometrical ordering property Then one obtains by wedge-locality and the decay properties described in (19) and (22) a rapid decay for large |τ |, Here it should be especially noted that the geometrical configuration depends on whether the outgoing (τ > 0) or incoming (τ < 0) regime is considered. Pairs of opposite wedge configurations were used in the previous constructions of two-particle scattering states in the wedge-local context [GL07;BS08], and also using wedge-local operators in local QFT context [BBS01]. In the multi-particle generalization (26) opposite wedges appear only indirectly. Namely, we work with one-particle states that can be generated from the vacuum within two opposite wedges W, In this case we call Ψ swappable (with respect to the wedge W). Swappable vectors with W ⊥ := W ′ are dense in the full Hilbert space H as a consequence of wedge duality (HK2 ♯ ), see [Du18, App. B]. Swappable one-particle states are obtained by projecting swappable vectors onto the one-particle space H 1 . In this way we obtain from (16) and (20) an oppositely localized pair of Haag-Ruelle creation-operator approximants for each one-particle state (1 ≤ k ≤ n) where both expressions involving Haag-Ruelle operators are τ -independent by construction, so that τ → ±∞ limits can be dropped. Although scattering states can be constructed without explicit use of the swapped operators B ⊥Λ kτ (f k ), they are the main tool for establishing the wedge-local Haag-Ruelle theorem: , Ω) be a wedge-local quantum field theory satisfying wedge duality (HK2 ♯ ) and the mass gap condition (HK6 ♯ ). Let Λ ∈ L ↑ + and Ψ j 1 = E(H m )A j Ω with A j ∈ A(W) be swappable one particle states, and define B k := A k (χ) with an auxiliary function χ as in (16).
(i) For regular positive-energy Klein-Gordon solutions f j satisfying the scattering state approximants Ψ Λ n (τ ) : converge rapidly in norm for τ → ∞. More precisely, for any N ∈ N there exists a C N > 0 such that )Ω, constructed both with respect to the same wedge W, satisfy Here the right hand side is again τ -independent by construction.
Analogous statements hold for incoming scattering states, assuming opposite ordering.
For the later discussion of the scattering matrix, let us note that all asymptotic data in wedge-local models, including wave operators, must be defined depending on a localization wedge W from which the scattering states have been prepared. This is an unusual feature of wedge-local models. But we note that, of course, the W-dependence can be trivial. This happens in 1 + 1 dimensions, where our results apply as well. In this case there are only two centered wedges W R and W L , whose associated observables and scattering states can be related by swapping symmetry. In general this certainly does not imply that the models will be local. Presently the status of existence or non-existence of local observables even in 1+1-dimensional GL-models appears to be still open. In the present formalism for higher dimensions, the existence of a local QFT model underlying the wedge-local model under consideration implies a certain trivial W-dependence of scattering states, which is discussed in Section 5.
The possible more general wedge-dependence of scattering states is an interesting feature of wedge local quantum field theories in higher dimensions. In particular, we will use the formalism introduced in [Du18] to describe this wedge dependence in a more transparent manner for twoparticle and n-particle scattering reactions. One of the main aims of the present paper is to illustrate this wedge dependence in BLS-deformed wQFT models and in particular for the special case of the Grosse-Lechner models (see Sections 3.3 and 5).

N -Particle Scattering in BLS-Deformed wQFTs
Using the results and notation described in the previous sections, we can now state our results in precise form. Let (A 0 , α, H , Ω) be a given wedge-local quantum field theory satisfying wedge duality (HK2 ♯ ) and the mass gap condition (HK6 ♯ ), in addition to (HK1)-(HK6). Let (A Q , α, H , Ω) be the model constructed by BLS-deformation with some fixed warping parameter Q.
Our main aim will be to prove asymptotic completeness of the deformed model. On the technical side, this will be achieved by establishing a direct relation between a scattering state of the deformed model, and the scattering states of the initial model, Here we write , as before, and A Q k ∈ A Q (W) for 1 ≤ k ≤ n with some fixed centered wedge W = ΛW R . In the following we will for simplicity suppress dependencies on the fixed reference frame boost Λ and on the wedge in our notation, writing for , when these dependencies can be clearly seen from the context.

Wave Operators in BLS-deformed wQFTs
Our introductory heuristic considerations (see equations (7) and (8)) suggested that the deformation preserves outgoing and incoming particle numbers. This makes our approach feasible. As we can in principle deform general models, such as P (φ) 2 , the scope of our result is not restricted to particularly simple or integrable deformed models. However, we can infer similarly that the deformed scattering state (35) cannot be directly related to a scattering state in the undeformed model with the same simple product form (36). This is addressed here by using the formulation of the scattering data in terms of wave operators W ± W , as defined in [Du18, Sec. 5]. To recall the formal construction of the wave operators we start by introducing the full (unsymmetrized) Fock space by over the one particle space Here it is convenient to define the velocity ordering ≻ W also as partial order on H 1 , by writing where the velocity support V Λ Ψ1 of a one-particle vector Ψ 1 is defined analogously to (19) and (22), replacing the support off j by the spectral support of the momentum operator in the state Ψ 1 . Natural domains for the wave operators in the wedge-local setting are the velocity-ordered Fock spaces Γ ≻ (H 1 ) and Γ ≺ (H 1 ), which are defined as the closures of the spans of finite linear combinations, for outgoing-and incoming scattering states, respectively. For technical purposes we analogously define such algebraic spans Γ ≻ 0 (H ′ 1 ) and Γ ≺ 0 (H ′ 1 ) for dense subsets H ′ 1 ⊂ H 1 .
Our main technical results concern the wave operators W ± Q,W of the deformed model, and we will now discuss their construction as given in [Du18]. Here we will also include some additional details which are of particular importance for us here. First we recall that the wave operators are constructed by means of Theorem 6. In particular this means that, on the technical side, we are in fact working with a smaller subset of one-particle states, namely all which can be written in the form We call such Ψ 1 swappable (with respect to W) one-particle states of bounded energy, and we let H W 1c be the (non-closed) linear space spanned by them. Alternatively, using spectral calculus, one sees that the above generating set 3 of Ψ 1 ∈ H W 1c can also be characterized by the existence Proposition 7. In a wQFT satisfying wedge duality (HK2 ♯ ), the density H W 1c = H 1 holds for any wedge W.
Remark 8. Concerning this technical detail, there is one subtlety with respect to the standard construction of bosonic and fermionic Fock spaces, which should be pointed out. Namely, we expect the resulting Fock space to be independent of possibly different choices of H ′ 1 , as the choice of this dense subset of the one-particle space is merely technically motivated. Yet, for ordered Fock spaces the choices of H ′ 1 are not completely arbitrary due to a possible interplay with ordering conditions. In the present context one can express this requirement on H ′ 1 as follows: we say that H ′ 1 ⊆ H 1 is momentum resolving iff for any Ψ ∈ H 1 there exists a sequence (Ψ n ) n∈N ⊂ H ′ 1 such that Ψ n → Ψ and supp E (H,P ) Ψ n ⊆ supp E (H,P ) Ψ. We note that from Proposition 7 and (42) it follows by choosingf suitably that the linear spaces H W 1c are momentum resolving.
By construction, the energy-momentum operators of the initial and deformed model coincide.
In particular, the one-particle spaces of the two models are identical. Hence we can directly compare the scattering data of the two models on the present abstract level in terms of the respective wave operators, as they are defined on the same ordered Fock spaces for both the deformed and initial model. For a direct comparison of the two wave operators, we will make use of a strenthened technical result on the density of one-particle states.
Corollary 9. Any state Ψ 1 ∈ H W 1c can be generated by a pairÃ ∈ A r (W) andÃ ⊥ ∈ A r (W ⊥ ) of regular operators.
Proof. Starting from a pair A ∈ A(W), A ⊥ ∈ A(W ′ ) and regular wave packetf from Proposition 7 we write where χ ∈ C ∞ c (R d ) is a compactly supported function withχ(ω m (p), p) non-vanishing for all p ∈ suppf . Thenf ′ (p) :=χ(ω m (p), p) −1f (p) defines a new regular Klein-Gordon wave packet andχ(H, P )AΩ = (2π) −d/2 A(χ)Ω =: A ′ Ω (see eq. (17)) shows that the vector part is obtained from the vacuum as image of a regular operator A ′ ∈ A r (W) for some wedgeW ⊃ W + supp χ. An analogous calculation yields A ′⊥ ∈ A r (W ⊥ ) with W ⊥ ⊃ W ′ + supp χ. We can then arrangẽ A ∈ A r (W) by translating the above constructed operator back to W, and also replacing A ⊥ and W ⊥ with the respective translates by the same vector. These translations can be compensated by absorbing a corresponding phase into the wave packet.
Definition 10 (wave operators of initial and deformed model). Let (A, α, H , Ω) be a wedge-local quantum field theory. The incoming and outgoing wave operators W ± 0,W associated to centered wedge regions W = ΛW R , Λ ∈ L, in the initial model are the maps defined by via linear and continuous extension from product states in Γ . Similarly the wave operators of the deformed model (A Q , α, H , Ω) are For simplicity we will refer to W ± Q,W as the deformed wave operators. We note that, strictly speaking, this terminology is justified only in retrospective after the results of the present paper are established. Namely, our main result shows the wave operators W ± Q,W , as defined using the general construction in the deformed model, can be regarded as a "deformation" of the wave operators W ± 0,W from the underlying "undeformed" model.
We note that P i can be regarded as i-th unordered second quantization of the energy-momentum operator P . Essential self-adjointness of P i on the domain of vectors of finite particle number follows from standard arguments (see e.g. [RS1, Sec. VIII.10]).

Proof of the Wave Operator Identity
We start with some preparations. First we will recall how to re-express the warped convolutions as convergent H -valued integrals. This is done by using standard oscillatory integral methods, similarly to the discussions from [BLS11; DT13].
Lemma 12 (oscillatory integral method). Let R 2d ∋ (x, y) −→ Ψ(x, y) ∈ H be a map with uniformly bounded derivatives for all multi-indices β ∈ N 2d 0 up to order |β| ≤ 4d. 4 Then we have Here D reg (∂ x , ∂ y ) := s j=0 D j (∂ xj , ∂ yj ) 2 is a product of auxiliary, mutually commuting partial differential operators defined by In particular, the limit (49) exists and it is independent of the choice of the regularizing function Proof. The operators D j are constructed such that their formal adjoints satisfy Inserting this identity into the oscillatory integral with finite ǫ > 0 twice for every 1 ≤ j ≤ d, using integration by parts, and writing η ǫ (x, y) := η(ǫx, ǫy), we get In this form the ǫ → 0 limit can now be carried out: it follows by explicit calculation and induction that where the norm is defined by From (53) we now obtain an integrable majorant for the rewritten integral (52), by estimating H ) and, for 0 < ǫ < 1, η ǫ C 4d (R 2d ) ≤ η C 4d (R 2d ) using the definition of η ǫ and the chain rule. Hence, by dominated convergence, it is sufficient to verify pointwise convergence of the integrand. To this end we write In each term of the product rule expansion of the remainder R ǫ there is at least one derivative with respect to x or y acting on η ǫ . Thus this remainder is proportional to ǫ and vanishes for ǫ → 0. The first term converges pointwise to D reg (∂ x , ∂ y )Ψ(x, y) and from this the claim (49) follows.
For our scattering theoretic purposes it is useful that this method can be readily extended to obtain similar integral representations by iterating Lemma 12. In this manner we can further strengthen the decay of the integrand for large x, y ∈ R d .
Lemma 13. Let M ∈ N and let R 2d ∋ (x, y) −→ Ψ(x, y) ∈ H be a map with uniformly bounded derivatives for all multi-indices β ∈ N 2d 0 up to order |β| ≤ 4M d. Then we have The integrand satisfies the integrable norm bounds Next, we use the convergent integral representation to check that warped convolutions can be exchanged with the smearing operations used to define the creation-operator approximants in the deformed model.
, and let f be a regular Klein-Gordon solution. Then and Proof. We can apply Lemma 12 to warped convolutions, as the corresponding integrand Ψ(x, y) = α Qx (A)U (y)Ψ is arbitrarily often differentiable for A ∈ C ∞ and Ψ ∈ D, with ∂ α x ∂ β y Ψ(x, y) ≤ C αβ for suitable constant depending on the multi-indices α, β ∈ N d 0 . By translation covariance of warped convolutions from Lemma 3 (v) we have Using the decay of χ ∈ S (R d ) and estimate (53), the above integrand has integrable norm with respect to the product Lebesgue measure d 3d (x, y, z). Hence the order of integrations can be exchanged by Fubini's theorem. Taking the x-and y-dependent translation operators outside the inner strong integral, we obtain (58). The proof of (59) is analogous.
Here, P 1 = P ⊗ 1 is the energy-momentum operator acting on the first argument only, and similarly P 2 = 1 ⊗ P acts on the second argument.
Before we continue let us remark that this deformed tensor product preserves the ordered subspaces in the sense that if Ψ, Ψ ′ ∈ Γ ≻ (H 1 ), with Ψ ≻ Ψ ′ , then also Ψ ⊗ Q Ψ ′ ∈ Γ ≻ (H 1 ). We cautiously note that this deformed tensor product is clearly linear in its arguments and associative, but not commutative (unless Q = 0). This is somewhat reminiscent of the general structure of the Zamolodchikov-Faddeev relations in integrable models, in contrast to canonical commutation relations (see e.g. [Le03]). Further, ⊗ Q is in general not mixed-associative in combination with ordinary tensor products. That is, Ψ 1 ⊗ Q (Ψ 2 ⊗ Ψ 3 ) = (Ψ 1 ⊗ Q Ψ 2 ) ⊗ Ψ 3 . We also note that the definition is consistent with the fact that on Fock spaces Ψ ⊗ Ω = Ψ = Ω ⊗ Ψ are identified for any Ψ.
We can now establish the main lemma for proving the wave operator identity from Theorem 11.
With incoming ordering V 1 ≺ . . . ≺ V n we have analogously for the incoming limit τ → −∞ Proof. We consider only the case τ → ∞. We first note that Ψ kτ := B kQτ (f k ) . . . B nQτ (f n )Ω ∈ D for 1 ≤ k ≤ n, due to the compact energy-momentum transfer of the B jQτ (f j ), 1 ≤ j ≤ n (see [Du18] Lemma 7). Thus Lemmas 14 and 12 apply and we obtain where the constant M ∈ N will be chosen below. We rewrite the vector part of the integrand using that U (x 1 )α Qy1 (B 1τ (f 1 )) = α Qy1+x1 (B 1τ (f 1 ))U (x 1 ) as with x = (x 1 , . . . , x n ), y = (y 1 , . . . , y n ) ∈ R nd , and we abbreviate Q := Q W . Now we split the integration in (64) into an integral over the region and its complement, where the constant ρ > 0 is for now fixed but arbitrary, and will be specified later. The vector resulting from the integral over R ↑ ρτ will be denoted Ψ ↑ 1τ . By means of applying (57) iteratively for all pairs (x k , y k ), 1 ≤ k ≤ n, we verify that with an appropriately strong power M ∈ N of the regularizing operator D reg , the remainder integral over the complement R ↓ ρτ := R 2nd \ R ↑ ρτ becomes small for large τ . That is, In the last step we estimated The derivative norm of the vector part was then bounded by expanding it into individual differentiated terms ∂ α x ∂ β y Ψ x,y,τ , which in turn can be expanded by the product rule into terms with differential operators acting on the translated creation-operator approximants B kτ (f k ). Due to the assumed norm differentiability of A k , the differentiated B kτ (f k ) can be rewritten as creation operator approximants constructed using the differentiated A k . These will be denoted here byB kτ (f k ) andÃ k , respectively. Using this rewriting, the terms from the expanded differentiated vector parts from (67) can each be bounded using the standard estimate B kτ (f k ) ≤ C χ,f k Ã k (1 + |τ | s/2 ) to obtain α z1 (B 1τ (f 1 )) . . . α zn (B nτ (f n ))Ω ≤ n k=1 B kτ (f k ) ≤ C(1 + |τ | ns/2 ). Taking all these terms together we obtain the last step from (68), where the new constant C depends on M , on all wave packages and on the norms of derivatives of the A k operators up to order 4M nd. To obtain the decay estimate for sufficiently large τ we proceed to estimate Here constant factors such as the convergent integrals can be absorbed into the constant C.
To bound the supremum we note that by definition of R ↑ ρτ we have for any point (x, y) in the complement R 2nd \ R ↑ ρτ at least one vector with |x 0 k * | + |x k * | > ρτ or |y 0 k * | + |y k * | > ρτ . This in turn implies that at least one coordinate |x k * ,j * | or |y k * ,j * |, respectively, is larger than ρτ / √ 4s. As all other factors in the supremum are bounded from above by one, we obtain Choosing M ∈ N large enough, this contribution becomes arbitrarily small for large τ > 0.
In particular, for establishing (62) it is sufficient to consider the limit of Ψ ↑ τ . To address this convergence, letρ > 0 denote the minimum over all such constants from Lemma 17 for the families of (B jτ (f j )) 1≤j≤n , whereB jτ (f j ) = ∂ βj xj α xj (B jτ (f j )) xj=0 stand for all possible combinations of derivatives of these operators up to order |β j | ≤ 4M nd, β j ∈ N d 0 . It is used here that ∂ βj xj α xj (B jτ (f j ))Ω = ∂ βj xj α xj (B ⊥ jτ (f j ))Ω for all multi-indices β j ∈ N d 0 , which shows by setting x j = 0 that all suchB jτ (f j )Ω are swappable, as needed to apply Lemma 17 to these differentiated families. Further the derivatives with respect to y k can also be written as multiples of such differentiated operators by using the chain rule. Thus we set ρ :=ρ · (1 + Q ) −1 /(n + 1) < ∞. With this choice we have for any 1 ≤ k ≤ n and (x, y) ∈ R ↑ ρτ that the vectors z k := Qy k + x 1 + . . . + x k are contained in the double cone of radius kρ |τ | + Q ρτ ≤ρ |τ | for 1 ≤ k ≤ n. Considering the asymptotically dominant part we can obtain a τ -uniform integrable bound by applying (69) and noting that instead of the coarse bound (71) we can now estimate Here we make use of Lemma 17, noting that this uniform bound applies to each of the terms resulting from the expansion by our choice of ρ. Due to the restriction to R ↑ ρτ we can estimate the differentiable vector norm by means of Lemma 17, after expanding it into individual terms with fixed derivatives acting on the creation-operatorsB kτ (f k ). As before these can be written as Haag-Ruelle-type operators in terms of the differentiated A k , so that Lemma 17 applies. To summarize we note that the proof strategy here is in fact analogous to the above method used for the outside region R ↓ ρτ . However in R ↑ ρτ the use of the clustering bound of Lemma 17 is geometrically permitted and yields the much stronger τ -uniform estimate on the vector part.
By dominated convergence we obtain where we already used that the characteristic function 1 R ↑ ρτ (x, y) → 1 pointwise for τ → ∞. Here the Haag-Ruelle limit can be exchanged with the regularizing differential operators by explicit computation: we expand everything into differentiated translated Haag-Ruelle operators, for which the right hand side of (75) converges to the scattering state generated by the differentiated and translated operators, which are again Haag-Ruelle type creation-operator approximants. Collecting the Haag-Ruelle limits again after performing them, we have by linearity of the wave In the last equality we have written the one-particle state translations again in groups using the second quantized translations on the unordered Fock space. In this form we can now apply the warped-convolution-type integral representation of the deformed tensor product (Proposition 20). For this purpose we note that the wave operator is bounded on the range of the integrand and can therefore be taken outside the strong integral. Thus the introduction of the regularizing differential operators can be undone by iterative application of Lemma 13 and we obtain lim τ →∞ Regrouping the convergent integrals, using Fubini and continuity of the tensor product, and applying Proposition 20 iteratively, we obtain deformed tensor products, as claimed in (62). We note that after replacing all tensor products by Q-deformed tensor products the right-associative grouping from (77) becomes inessential and can be dropped. The proof of the statement for the incoming limit τ → −∞ is analogous.
We used the following auxiliary result concerning the norm of scattering-state approximants involving certain restricted translations of each operator.
Before proving this lemma let us recall the following useful technical result from [Du18], which concerns the approximation of Haag-Ruelle creation-operator approximants by wedgelocal operators.
Lemma 18 ([Du18] Lemma 9). Let A ∈ A(W). For any τ ∈ R and δ > 0 the corresponding where the constants C δ N depend on f , A and χ.
For the proof of Lemma 17 we use a corresponding version of the commutator estimate, which will be formulated as a separate lemma to be proven first. For this estimate certain translations of oppositely localized pairs are admitted, similarly to the corresponding translations appearing in Lemma 17.
, for a pair of opposite wedges W, W ⊥ , χ ∈ S (R d ), and let f, f ⊥ be regular Klein-Gordon solutions ordered by V f ⊥ ≺ W V f . Then there exists a constant ρ > 0 and for any N ∈ N a constant C N > 0 such that for any τ > 0 and x, y ∈ C ρτ , This rapid decay extends to the case that either or both operators are replaced by their adjoint. Further these estimates also hold for τ < 0 given the opposite ordering In the proof presented here we will focus on the arguments needed to generalize the corresponding commutator estimate of Corollary 10 from [Du18] to the present statement. More details can be found in Appendix A of [Du18].
Proof. We only consider the outgoing case τ > 0 and note that it is sufficient to prove (80) for all τ ≥ τ 0 with some fixed τ 0 > 0. Let us begin by assuming for simplicity that y = 0. For δ > 0 we obtain families of wedge-local approximants B and analogously B . To obtain (80) by means of the wedge-locality of these two approximating operators we have to choose δ > 0 and subsequently ρ > 0 sufficiently small so that the localization regions are space-like separated. These causality considerations can be simplified for large enough τ by absorbing any finite translations into the growing double cones and rewriting as a region involving a single growing double cone: Choosing for simplicity ρ = δ and assuming τ > τ 0 := (|x W | c + |x W ⊥ | c )/δ we obtain By the ordering assumption we have that V f − V f ⊥ is a compact subset of the open set W c . In particular there exists an Thus we can choose ρ = δ := ǫ/6 and it then follows that the two sets M τ 1 and M τ 2 are space-like separated for τ > τ 0 .
We now obtain from locality that for all τ ≥ τ 0 and x ∈ C ρτ we have [B τ )] = 0, which implies the uniform commutator estimate by expanding )·(1 + |τ |) s/2 ≤ C ′ N τ −N by estimating the commutator via the two operator norms, using Lemma 18 for its first and the standard polynomially growing norm estimate for its second argument, and analogously for the other non-vanishing commutator.
Concerning the extensions to adjoints and τ < 0 we note that the commutator estimates involving adjoints follow using the same approximation argument, noting that (18) also directly yields wedge-local approximants of adjoint operators by the C * -property of the operator norm. For τ < 0 the geometric situation in (81) is inverted, but the same arguments work if the opposite ordering V f ⊥ ≻ W V f holds.
Proof of Lemma 17. For the outgoing case, let ρ = ρ n > 0 be the minimum over the constants ρ k,j from Lemma 19 for the pairs B kτ (f k ) and B ⊥ jτ (f j ), 1 ≤ k < j ≤ n. For n = 1, estimate (78) follows directly from the fact that Ψ x1 1 := α x1 (B τ (f ))Ω = U (x 1 )B τ (f )Ω does not depend on the translation vector x 1 due to translation invariance of the norm and the vacuum, neither on τ by construction of the Haag-Ruelle operators.
For the case n ≥ 2, let Ψ x n n (τ ) := α x1 (B 1τ (f 1 )) . . . α xn (B nτ (f n ))Ω with x n := (x 1 , . . . , x n ) ∈ R nd . To simplify the notation we will drop the obvious wave packet dependences and write B x k kτ := α x k (B kτ (f k )), so that Ψ . To give an inductive argument we can write using the swapping property The first term can be bounded using the clustering property of the Haag-Ruelle operators ( [Du18] Prop. 8 (vi)), which yields Here the estimate is obtained uniformly for all x 1 , . . . , x n ∈ C ρnτ by using that x 1 , . . . , x n−1 ∈ C ρnτ ⊆ C ρn−1τ due to the induction hypothesis, and bounding the first factor as for n = 1. Note that here we use ρ n ≤ ρ n−1 , which follows directly by the definition where ρ k,j denote the constants from Lemma 19 for the pairs B kτ (f k ) and B ⊥ jτ (f j ), 1 ≤ k < j ≤ n.
Finally, the commutator term from (84) is also seen to be bounded in τ (in fact, rapidly decreasing), uniformly for x 1 , . . . , x n ∈ C ρnτ , by means of Lemma 19: we expand the big commutator in (84) into a sum of terms of vacuum expectation values of operators of the form Each of those terms is now estimated using that B x k * kτ , B ⊥xn nτ and B x k kτ , B ⊥xn nτ can for 1 ≤ k < n and x k , x n ∈ C ρnτ by definition of ρ n be bounded by C M (1 + τ ) −M . Choosing M ∈ N sufficiently large we absorb the growth of the simple norm estimates B xj * jτ = B xj jτ ≤ C j (1 + |τ | s/2 ) used for the remaining operators with 1 ≤ j ≤ n in the respective terms of the expansion. Together the desired uniform bound on (84) is obtained.
Proof. First we can infer from Lemma 12 that the limit in (88) exists and is independent of η within the specified restrictions. We rewrite the expression under the limit on the right-hand side via the spectral calculus of the energy-momentum operators as Here the Fubini theorem for exchanging the order of integrations applies, where integrability with respect to the product measure follows from the bounded energy-momentum of Ψ 1 , Ψ 2 and the rapid decay of η. Here P denote the energy-momentum operators on the full Fock space and P 1 (Ψ 1 ⊗ Ψ 2 ) := (P Ψ 1 ) ⊗ Ψ 2 acts only on the second component. Let us denote the inner scalar integral from (89) by I ǫ (p, q). By the uniqueness result of Lemma 12 we can proceed to concretely choose η(x, y) := e −|x| 2 e −|y| 2 e , where |·| e denote the Euclidean norm. Then an elementary calculation (see Proposition 21) shows that I ǫ (p, q) −→ e −ip·Qq pointwise for ǫ → 0, and in addition I ǫ (p, q) is bounded uniformly in p, q for small ǫ > 0. By dominated convergence we obtain that (89) yields for ǫ → 0 where we used that P = P 1 + P 2 and P 1 · QP 1 = 0.
Proposition 21. For any p, q ∈ R d and any warping matrix Q ∈ R d×d we have Proof. We let ε := ǫ 2 and use anti-symmetry of Q with respect to the Minkowski scalar product to write Let us express the Minkowski products in terms of ordinary scalar products involving the Minkowski metric g. Then we can calculate component-wise by Fubini, with Here we have substituted x = x ′ µ := g µν x ν . By elementary calculation one obtains from which the claim follows.
Proof of Theorem 11. We consider the outgoing case for a fixed wedge W. By definition the linear combinations of ordered product states are dense in Γ ≻W (H 1 ). Further by Proposition 9 such states can be approximated with arbitrarily small error by vectors generated by swappable Haag-Ruelle approximants Here the further restriction from swappable states to swappable states generated by regular operators is used to assure that the warped convolutions are well defined and we have just seen that the linear combinations of the vectors (97) are also dense in Γ ≻W (H 1 ) even with this additional smoothness requirement. By Lemma 3 (i) it is clear that they yield the same one-particle vectors B kQτ (f k )Ω = B ⊥ kQτ (f k )Ω = Ψ k 1 for all 1 ≤ k ≤ n. For simplicity of notation we will be dropping the tilde for the remainder of the proof.
By definition of the wave operator we now obtain where the last equality follows from Lemma 16. By induction we obtain that Hence the claimed identity holds on a dense subspace and thus by continuity of the wave operators and S Q on the full domain Γ ≻W (H 1 ) of W + Q,W . The argument for W − Q,W is analogous.

Scattering Data and Wedge-Transition Matrix Elements
The two-particle S-matrix of higher-dimensional deformed or GL-type models was worked out in [GL07; BS08] for any fixed wedge W. Even if the initial and thus also the deformed model were Poincaré covariant, these authors observed that the two-particle S-matrix is in fact not fully Lorentz covariant in higher dimensions d > 1 + 1. In [Du18] it was proposed to make this observation more precise by specifying the dependence of wave operators and all related asymptotic data on the localization wedge W of the Haag-Ruelle operators explicitly. Due to translation covariance of the S-matrix this reduces to a dependence on a wedge modulo translations, or equivalently, a dependence on a centered wedge.
Definition 22. Let W f , W i be centered wedges. Following [Du18] we define S-matrices of the initial and deformed wQFT model, respectively, as maps between the incoming ordered Fock space Γ ≺W i (H 1 ) and the outgoing space Γ ≻W f (H 1 ) by Similarly we define wedge-transition maps between two final or two initial states, respectively, as for any two centered wedges W 1 , W 2 .
As a direct consequence of Theorem 11 we obtain similar expressions for the S-matrices and wedge-transition matrix elements in deformed wedge-local models.
Corollary 23. Let W f , W i , W 1 , W 2 be arbitrary centered wedges. The S-matrices and wave operators of a BLS-deformed wQFT model can be expressed in terms of the corresponding objects of the undeformed model: Proof. We obtain the expression for the S-matrix (i) from the short computation The calculations for the wedge-transition formulas (ii) are analogous.
Remark 24. Expressions (i) and (ii) here simplify further by noting that it follows immediately from the definition of S QW and its restrictions S where in the last equality the definition of Q W was used. This further implies that these deformation maps are unitary, by writing and analogously S ≻W QW (S ≻W QW ) * = 1. Let us note that we can further identify the scattering data for opposite wedges by making use of the swapping property. For this consideration it is immaterial whether we are in the deformed or undeformed model, so we will drop the corresponding indices from the wave operators. We start with an outgoing scattering state given by By an analogous argument as in (84) and (87) we write where the commutator terms are rapidly decreasing faster than any polynomial in τ > 0. Iterating this swapping argument we obtain In the last equality we already used that by swapping B ⊥ kτ (f k )Ω = B kτ (f k )Ω and that by the definition of the precursor relation and from the fact that W ′ = −W for any centered wedge W. Let us therefore define Z : Γ u (H 1 ) → Γ u (H 1 ) by its action on n-particle states, for any n ∈ N and Ψ 1 1 , . . . , Ψ n 1 ∈ H 1 and we note that Z is a unitary and self-adjoint involution. Summarizing these considerations we obtain: Proposition 25. For any centered wedge W we have as subspaces of Γ u (H 1 ), (ii) ZΓ ≻W (H 1 ) = Γ ≻ W ′ (H 1 ), and analogously ZΓ ≺W (H 1 ) = Γ ≻W (H 1 ).
Proposition 26. For any centered wedge W the wave operators associated to complementary wedges can be identified by This identification appears to take us somewhat away from the localization properties appearing in the Haag-Ruelle construction. Yet it shows that the scattering matrices for the two simplest choices W f = W i and W f = W ′ i contain the same scattering theoretic data. These choices correspond to ones made in [GL07;BS08] for the analysis of two-particle scattering. In these works asymptotic two-particle states are constructed following more closely the methods of Haag-Ruelle theory from local QFT. Hence a mixed localization is used, where The independence of the operator order is reminiscent of Haag-Ruelle scattering theory for bosonic local QFT. In the general wedge-local scattering theory it is a special feature appearing for the case of two particles. It was already remarked in [Du18] that this definition of twoparticle scattering states captures the same information as our analysis restricted to the level of two-particle states. This can be seen by swapping the corresponding sides, Thereby we have not only illustrated the symmetry of Proposition 26 at the two-particle level.
We also see that the compatibility of our analysis with earlier calculations from [GL07; BS08] is in fact a corollary of an intermediate step of the proof of Proposition 26.

Asymptotic Completeness of BLS-Deformed wQFT
Given a wedge-local model (A, α, H , Ω), we can now proceed to our second main objective and study the completeness of asymptotic states in the common Hilbert space H of the initial and BLS-deformed wQFT model. We will use a notion of asymptotic completeness for wQFT which directly generalizes the corresponding standard definition from local QFT.
Definition 27. We say that a wave operator W ± W of a wQFT model for a centered localization wedge W is asymptotically complete iff the subspace of velocity-ordered scattering states respectively, is dense in H . We say that a wQFT model satisfies the property of ordered asymptotic completeness (more precisely, asymptotic completeness with respect to velocityordered scattering states) iff both W + W and W − W are asymptotically complete for any wedge W. The results of Section 3 give an explicit representation of the deformed wave operators W ± Q,W in terms of the undeformed wave operators W ± 0,W . This directly yields a general result regarding the stability of asymptotic completeness of wedge-local theories under BLS-deformations.
Taking the closures yields the equivalence of ordered asymptotic completeness of deformed and initial model for outgoing states. The argument for the incoming case is analogous.
In the following final section we will discuss the application of our results to GL-type models, which are constructed by applying BLS-deformations to a free field. For the scattering theoretic analysis of these models we have to work with only wedge-ordered states, as stated explicitly in (114) and dictated by the scope of the wedge-local scattering theory. In GL-type models we will see this restriction to wedge-ordered states is inessential for the particle interpretation and still yields a dense set of scattering states. This may of course be expected on grounds of the bosonic statistics of the underlying free theory. Let us also note that there are wedge-local models, which do not satisfy ordered asymptotic completeness at the two-particle level. Examples of such models have been obtained by applying a von Neumann operator-algebraic free product construction to the free field [LTU17].

Application to Grosse-Lechner Models
In the work of Grosse and Lechner [GL07] an interesting class of wedge-local models in any spacetime dimension d ≥ 1 + 1 is constructed in a wedge-local variant of the Wightman framework. A closely related class of models is obtained in the operator-algebraic framework by applying the BLS-deformation construction to the standard scalar free field [BS08; BLS11]. The scalar field has a canonical wedge-local description given by the von Neumann algebras act on the bosonic Fock space H = Γ b (H 1 ) over the scalar one-particle space H 1 = L 2 (R s ), and we write Ω F ∈ H for the Fock vacuum. The net is covariant with respect to the standard second quantized scalar representation of the proper orthochronous Poincaré group. For any non-vanishing warping matrix Q ∈ R d×d of the form (11) we will call the BLS-deformed theory (A Q , α, H , Ω) a Grosse-Lechner model. Let us now apply our general results from Sections 3 and 4 to evaluate the scattering theoretic content of these models and establish their ordered asymptotic completeness. It is a well-known fact that in the case of the free field also the conventional local Haag-Ruelle scattering theory applies (see e.g. [A;Dyb17]). It yields wave operators W ± : Γ b (H 1 ) → H , defined on the bosonic Fock space. For the free model we have H = Γ b (H 1 ) and the wave operators are trivial in the sense that W + = W − = 1. 5 This also directly implies asymptotic completeness of the free theory in the conventional sense with respect to the standard Haag-Ruelle scattering theory. That is, W ± Γ b (H 1 ) is dense in H and in these cases we will say similarly that W ± , respectively, are asymptotically complete.
For the free field, the present velocity-ordered wedge-local formalism for scattering theory applies as well. The wedge-local wave operators W ± 0,W can be efficiently determined from their local counterparts in cases for which the latter exist. For this purpose we define the embeddings which map wedge-ordered n-particle vectors to the corresponding bosonic symmetrized tensor product. In terms of the norm on Γ b (H 1 ) we get Here it is used that, as a consequence of ordered velocity supports, the one-particle states Ψ k 1 ∈ H 1 (1 ≤ k ≤ n) are pairwise orthogonal. Thus I ≻W /≺W are well defined by linear continuous extension from (117) and yield isometries on Γ ≻W /≺W (H 1 ).
Theorem 29. The wedge-local wave operators of a local quantum field theory with isolated mass shell are well-defined and can be expressed using local Haag-Ruelle wave operators W ± : Further, the S-matrix and wedge-transition maps in the case of a local quantum field theory are given by where S f i = (W + ) * W − denotes the usual scattering matrix from local Haag-Ruelle theory.
Proof. It is sufficient to establish (119), from which the other statements follow. Here we note that for the wedge-transition matrix identities (121) also the Fock structure (W + ) * W + = 1 = (W − ) * W − of scattering states in the local Haag-Ruelle theory is used. We consider the case of the outgoing wave operator. We may restrict to one-particle vectors of the form Ψ k = B kτ (f k )Ω defined in terms of local operators A k ∈ A(O) and regular Klein-Gordon solutions f k for 1 ≤ k ≤ n. In a local QFT model such vectors yield a dense subset of the one-particle space by standard arguments. We further assume that the one-particle states (and similarly the wave packets) are velocity ordered Ψ 1 ≻ W Ψ 2 ≻ W . . . ≻ W Ψ n . Lastly, these states trivially satisfy the swapping property with A ⊥ k := A k for some overlapping wedge-regions W ⊃ O, W ⊥ ⊃ O. Thus, we have on one hand by definition of the wedge-local wave operators that On the other hand, we obtain from standard local Haag-Ruelle theory that lim τ →∞ B 1τ (f 1 ) . . . B nτ (f n )Ω = W + a * (Ψ 1 ) . . . a * (Ψ n )Ω F = W + I ≺W Ψ 1 ⊗ Ψ 2 ⊗ . . . ⊗ Ψ n . (123) Equating we obtain (119) on a total subset and thus by continuity on Γ ≻W (H 1 ). The argument for the incoming case is analogous.
Using that the spectrum of the momentum operators is Lebesgue absolutely continuous on the one-particle space, I ≻W /≺W are in fact surjective. Such continuity clearly holds in the free example (116). In a general context it is known to follow from locality and the spectrum condition [BF82, Prop. 2.2], or, from Poincaré covariance [Mai68].
Proposition 30. For a local QFT model the following statements are equivalent: (i) W + is asymptotically complete, (ii) W + W is asymptotically complete for one wedge W, (iii) W + W are asymptotically complete for all wedges W, and analogously for the completeness of incoming wave operators.
Proof. The equivalence of (i) and (ii) for any wedge W follows from I ≻W /≺W being a surjective isometry and (119). As W was arbitrary, this implies (iii).
Corollary 31. The wedge-local wave operators of the free scalar field are given by In particular, the free scalar field satisfies the property of ordered asymptotic completeness.
Theorem 32. The Grosse-Lechner models are asymptotically complete with respect to velocityordered scattering states for any warping matrix Q (as defined in (11)).
Proof. We consider a GL-model for any fixed warping matrix Q. Theorem 29 shows that the wedge-local wave operators W ± 0,W of the scalar free field are asymptotically complete for any wedge W. By Theorem 28 we obtain ordered asymptotic completeness of the wave operator W ± Q,W for any W and, thereby, ordered asymptotic completeness of the considered GLmodel.
To conclude these investigations, let us note the explicit n-particle scattering matrix for the case of the Grosse-Lechner models. We will focus on the two cases of equal or opposite initial and final wedges, which correspond to the earlier investigations from [GL07;BS08]. The case W f = W ′ i can perhaps be regarded as slightly more natural due to the coincidence of Γ ≺W (H 1 ) = Γ ≻ W ′ (H 1 ) as subspaces of Γ u (H 1 ), which is a consequence of the equivalence of the ordering relations ≻ W ′ and ≺ W .
Proposition 33. Let W be a fixed initial wedge. Then for the case of W f := W ′ we obtain the Grosse-Lechner S-matrix as unitary operator on Γ ≺W (H 1 ) given by Further we have for any ordered n-particle state Ψ n = Ψ 1 1 ⊗ . . . ⊗ Ψ n 1 ∈ Γ ≺W (H 1 ) that In particular, the Grosse-Lechner S-matrix is non-trivial and factorizing.
This follows directly from Corollary 23 and Theorem 29, together with the triviality of the scattering matrix of the free field. Further we also used that ≻ W ′ and ≺ W are equivalent, giving (S We note that the same result is obtained when using Proposition 26 and (125).

Conclusions
In this paper we showed stability of ordered asymptotic completeness under BLS-deformations in wedge-local QFT. We concluded that the Grosse-Lechner model is interacting and asymptotically complete in any spacetime dimension. We also showed that this model has a factorizing S-matrix, which is an unusual feature in higher dimensions. Although we focussed on models obtained by BLS deformations, our approach should also apply to wedge-local factorizing models in two dimensions. Even in the cases in which strict locality is still open, such as the non-linear sigma models [AL17], our strategy may give asymptotic completeness and factorization of the S-matrix. We leave detailed analysis of these problems to future investigations.
Another natural direction is a generalization of our results to theories of massless particles. The BLS-deformations remain valid for such theories, but so far neither interaction nor asymptotic completeness have been studied for d > 1 + 1. Such an investigation would require a massless variant of wedge-local scattering theory from [Du18]. It may be difficult to develop such a theory at the same level of generality as its local counterpart [Bu77], since the decay of correlations in massless wedge-local models may be very slow. Also the energy bounds [Bu90], which simplify more recent constructions of massless scattering states [DH15; AD17], are not available for wedge-local theories. But under some natural assumptions on the decay of correlations massless wedge-local scattering theory appears to be within reach. It should apply, in particular, to the massless Grosse-Lechner model.
We mention as an aside, that such a scattering theory could also help to understand recent computations of infraparticle scattering amplitudes in four-dimensional string-local models [GRT21]. An apparent breakdown of unitarity of scattering amplitudes was found in this reference after adapting a formula from a two-dimensional context [DM21]. This problem may have its roots in our limited understanding of collisions of massless Wigner particles in string-and wedge-local theories.