On the Zakharov-Mikhailov action: $4$d Chern-Simons origin and covariant Poisson algebra of the Lax connection

We derive the $2$d Zakharov-Mikhailov action from $4$d Chern-Simons theory. This $2$d action is known to produce as equations of motion the flatness condition of a large class of Lax connections of Zakharov-Shabat type, which includes an ultralocal variant of the principal chiral model as a special case. At the $2$d level, we determine for the first time the covariant Poisson bracket $r$-matrix structure of the Zakharov-Shabat Lax connection, which is of rational type. The flatness condition is then derived as a covariant Hamilton equation. We obtain a remarkable formula for the covariant Hamiltonian in term of the Lax connection which is the covariant analogue of the well-known formula"$H=Tr L^2$".


Introduction
Classical integrable field theories in two dimensions are characterised by the existence of a Lax connection which is on-shell flat and depends meromorphically on an auxiliary complex curve, typically the Riemann sphere. Determining whether a given field theory is integrable or not is usually a very difficult problem since there is no systematic way of constructing a suitable Lax connection, if one exists.
Over the past couple of years, however, two closely related general frameworks have emerged for explaining the algebraic and geometric origins of Lax connections in 2d integrable field theories. From a Hamiltonian perspective, an origin was proposed in [V1], and further developed in [DLMV2], based on Gaudin models associated with untwisted affine Kac-Moody algebras and the representation theory of such algebras. From a Lagrangian perspective, an origin was proposed by Costello and Yamazaki in [CY], following earlier work on integrable spin chains in [Cos1,Cos2,W,CWY1,CWY2], based on the introduction of surface defects in 4d Chern-Simons theory. A much older connection between Lagrangians for (hierarchies of) integrable field theories in 2d and field theories of Wess-Zumino-Witten type was pioneered in [N].
In the Hamiltonian formulation of integrable field theories, there is an important dichotomy between so-called 'ultralocal' and 'non-ultralocal' theories. This distinction is based on the classical r-matrix formalism [S1, S2], specifically on whether or not the r-matrix of the given integrable field theory is skew-symmetric [Se,M1,M2].
The affine Gaudin model perspective on integrable field theories was specifically developed in [V1] to address the problem of quantisation of non-ultralocal theories. Note that a related approach was used in [Z] to treat ultralocal field theories as Gaudin models associated with affine Kac-Moody algebras. On the other hand, it was demonstrated on examples in [CY] that both ultralocal and non-ultralocal field theories can be described from the perspective of 4d Chern-Simons theory. In the non-ultralocal case, further examples were shown in [DLMV3] to fit within this framework, and more recently a very general family of new non-ultralocal integrable field theories was also constructed using this approach in [LV] following [BSV]. By performing a Hamiltonian analysis of 4d Chern-Simons theory, it was shown in [V2] that in the case of non-ultralocal field theories this frameworks is, in fact, intimately related to that of affine Gaudin models. By contrast, a Hamiltonian analysis of the class of ultralocal theories from the perspective of 4d Chern-Simons theory has so far not been considered. The main purpose of this paper is to initiate such a study.
In fact, very recently, an independent line of research emerged in [CS1] where the classical r-matrix structure was derived in the context of covariant Hamiltonian field theory. In that setting, a covariant Poisson bracket replaces the standard Poisson bracket and the r-matrix determines the Poisson algebra satisfied by the whole Lax connection (a 1-form) and not just by its spatial component, called the Lax matrix.
Such results have been established successfully for ultralocal theories with rational r-matrix (nonlinear Schrödinger and modified Korteweg-de Vries) and trigonometric r-matrix (sine-Gordon). However, the generalisation to non-ultralocal theories has resisted all attempts so far. In particular, the famous example of the principal chiral model, which is intrinsically non-ultralocal, does not seem to be easily amenable to this covariant formalism. Nevertheless, a certain reduction of the principal chiral model dynamics can be reproduced by an ultralocal integrable field theory [FR], for which an action was obtained in [AHP]. This model can be seen as a special case of a large class of models with Lax pairs of Zakharov-Shabat type which derive from an action first introduced by Zakharov and Mikhailov [ZM]. Our observation is that this general class of models admits an ultralocal r-matrix structure of rational type and is therefore suited for a covariant Hamiltonian treatment.
The main goal of the present work is to begin exploring the covariant Hamiltonian structure of certain ultralocal integrable field theories which can be obtained from the 4d Chern-Simons perspective, using the Zakharov-Mikhailov class of models as our guiding example. The covariant approach to integrable field theories initiated in [CS1] is in contrast with the long tradition of analysing the standard Hamiltonian formulation of integrable field theories and may offer new insights when it comes to their (covariant) quantisation. An interesting by-product of our approach is the interpretation of the flatness condition of the Lax connection as a covariant Hamilton equation associated with a covariant Hamiltonian which we derive from the Zakharov-Mikhailov action.
In Section 2, we show that the Zakharov-Mikhailov action of [ZM] can be derived from 4d Chern-Simons theory. Since, in our case, the meromorphic 1-form ω appearing in the 4d Chern-Simons action is ω = dz, it has a double pole at infinity so we follow a similar approach to [BSV] by first regularising the action of 4d Chern-Simons theory. We then couple minimally the 4d gauge field A to a collection of Lie group valued fields {φ m } N 1 m=1 and {ψ n } N 2 n=1 localised along surface defects. In Section 3, we derive the covariant Poisson algebra satisfied by the Lax connection of the Zakharov-Mikhailov class of models. We also present the covariant Hamiltonian of the theory and derive a remarkable formula connecting it to the Lax connection. This formula represents the covariant analogue of the well-known formula relating a traditional Hamiltonian H with the Lax matrix L which we write schematically as "H = Tr L 2 ". We show that the flatness condition of the Lax connection takes the form of a covariant Hamilton equation, thus giving it a new interpretation in this context. The results of this section rely heavily on the variational bicomplex formalism as presented in [D,Chap. 19] and on ideas developed for instance in [K]. For a detailed account geared specifically towards the implementation of these ideas in 2d integrable field theories, we refer the reader to [CS1].
We shall work with the general linear group GL N with Lie algebra gl N of N × N matrices, following [ZM], but we expect our results to hold more generally for any semisimple Lie algebra. We denote the trace by Tr : gl N → R.
Let X := Σ × CP 1 . We shall use the notation for gl N -valued p-and q-forms on X, where p, q = 0, . . . , 4, u A , v B ∈ gl N and A, B are multi-indices with |A| = p and |B| = q so that {dx A } and {dx B } denote bases of 1-forms for the space Ω p (X) and Ω q (X), respectively.
2.1. Regularised 4d Chern-Simons action. Since the 2d integrable field theory we want to describe is ultralocal, we consider the meromorphic 1-form ω = dz. The Lagrangian of the corresponding 4d Chern-Simons theory is given by where CS(A) := Tr(A ∧ dA + 2 3 A ∧ A ∧ A) denotes the Chern-Simons 3-form and A is a gl N -valued 1-form on X which we can decompose as Note that there is no need to include a dz-component since this would drop out from the Lagrangian (2.1). The components A ξ , A η and Az are taken to be smooth functions away from the set of marked points {a m } N 1 m=1 and {b n } N 2 n=1 , but it will be important for later to allow A ξ and A η to be singular at those points. Specifically, we will assume that these components can be written locally as A ξ = (z −a m ) −1 B m,ξ near a m for m = 1, . . . , N 1 and as A η = (z − b n ) −1 B n,η near b n for n = 1, . . . , N 2 , where B m,ξ and B n,η are smooth functions on X. One easily checks that, despite the presence of these singularities, the Lagrangian (2.1) remains locally integrable near Σ × {a m } for m = 1, . . . , N 1 and near Σ × {b n } for n = 1, . . . , N 2 . However, since the 1-form dz has a double pole at z = ∞, the 4-form dz ∧ CS(A) is not locally integrable near Σ × {∞}. For this reason we need to regularise the action which we do following [BSV]. First, note that Here we have used the fact that Tr(A ∧ A ∧ A ∧ A) = 0 for any 1-form A ∈ Ω 1 (X, gl N ) by the cyclicity of the trace.
We can now rewrite the Lagrangian (2.1) of 4d Chern-Simons theory as where the first term is exact but has a double pole at infinity while the second term only has a simple pole and is therefore locally integrable near Σ×{∞}. We therefore define the regularised action of 4d Chern-Simons theory by dropping the exact term above and keeping only the second term, namely we set Note that we can continue to assume that A has no dz-component, namely it can be expressed as in where in the second line we have used the fact that dF for any g ∈ C ∞ (X, G). Indeed, the transformation of the curvature F (A) under a gauge transformation (2.5) is given by conjugation F ( g A) = gF (A)g −1 , hence the result follows by the invariance of the trace.

2.2.
Adding surface defects. We would like to modify the action (2.3) by adding to it terms which couple the 4d bulk field A to a collection of 2d fields localised on the surface defects Σ × {a m } and Σ × {b n } for m = 1, . . . , N 1 and n = 1, . . . , N 2 . We shall make use of the embedding ι . Specifically, to each marked point a m for m = 1, . . . , N 1 we associate a Lie group valued field φ m ∈ C ∞ (Σ, GL N ) which we think of as living on the surface defect Σ × {a m }. Likewise, to each of the marked points b n for n = 1, . . . , N 2 we associate a Lie group valued field ψ n ∈ C ∞ (Σ, GL N ), living on the surface defect Σ × {b n }. Let us also fix constant non-dynamical elements U Following the discussion of order defects in [CY], we now couple the 4d gauge field A to the collection of 2d fields φ m and ψ n on the surface defects by replacing the regularised 4d Chern-Simons action (2.3) with Here d Σ denotes the de Rham differential on Σ.
To maintain the gauge invariance of the action under (2.5) after introducing the surface defects, we need to let the 2d fields transform as It is straightforward to check that the extended action (2.6) is then gauge invariant 2.3. Bulk equations of motion. Consider the variation A → A + ǫa of the action (2.6), for some arbitrary a = a η dη + a ξ dξ + azdz ∈ Ω 1 c (X, gl N ) of compact support. This reads n for all n = 1, . . . , N 2 . As we will see below, the 2d action obtained from our 4d action with defects effectively gives equations of motion for U m and V n which are valued in the Lie algebra gl N . Without any particular model in mind, it is a matter of taste at this stage whether one wants to interpret the fields of the theory as being those Lie algebra elements or the group elements φ m and ψ n . In the former interpretation, the phase space is thus the (co)adjoint orbit through U n . In the first term on the right hand side above we have dropped a boundary term which vanishes since a ∈ Ω 1 c (X, gl N ) is of compact support. The curvature F (A) is given in components by (recall that (2.4) ensures that we can take A with no dz-component) Note that the last term does not contribute to the equation of motion. The Az equation of motion is then given by which will become the zero curvature equation for the Lax connection. On the other hand, the A η and A ξ equations of motion respectively read where the δ-functions, satisfying the property for any x ∈ C and any smooth function f on X, come from the fact that the surface defect terms are localised at z = a m or z = b n .
2.4. Lax connection. Given the resemblance of (2.9) with the zero curvature equation satisfied by the Lax connection, we would like to turn A into the Lax connection itself. There are two obvious issues with this.
The first main issue is that A has an additional dz-component compared to the Lax connection L = L η dη+L ξ dξ. We can eliminate this problem by focusing on field configurations with no dz-component. This will break some of the gauge invariance since we must now impose that (2.5) does not re-create any dz-component in the gauge field. In other words, we impose that (2.12) Az = 0,∂gg −1 = 0.
An obvious way to ensure the latter condition is to take g ∈ C ∞ (Σ, GL N ), i.e. g no longer depends on CP 1 . These residual gauge transformations will correspond to gauge transformations in the 2d theory. The next difference between A = A η dη + A ξ dξ and a Lax connection is that the former depends smoothly on CP 1 , with singularities at the marked points a m and b n of the form described in §2.1, while the latter is meromorphic on CP 1 . This issue is resolved by focusing again on a subset of gauge fields which satisfy the equations of motion (2.10). Having fixed Az = 0, these now reduce to Using the identity ∂zz −1 = −2πiδ(z) we deduce that a solution of the above is These expressions coincide with those for the dξ and dη-components U and V of the Lax connection from [ZM,(2) & (6)].
Note that if we have U 0 = d ξ hh −1 and V 0 = d η hh −1 for some h ∈ C ∞ (Σ, GL N ), cf. [ZM,(5)], then we can set them both to zero in (2.13) using a gauge transformation with g = h −1 . This would have the effect of fixing the residual gauge symmetry down to the global transformations and the Lax connection (2.13) would then have no constant term, i.e. it would take the form We will, however, keep the residual gauge symmetry for the remainder of this section, which will become the gauge symmetry in the 2d action.
2.5. Defect equations of motion. We may also consider the variation of the action (2.6) with respect to the 2d defect fields φ m and ψ n . Consider the variation φ m → e ǫαm φ m for arbitrary α m ∈ C ∞ (Σ, gl N ) with m = 1, . . . , N 1 and ψ n → e ǫβn ψ n for arbitrary β n ∈ C ∞ (Σ, gl N ) with n = 1, . . . , N 2 in the action (2.6). This gives Taking into account the solution (2.13) this then leads to the equations of motion for m = 1, . . . , N 1 and n = 1, . . . , N 2 . These coincide with [ZM,(4)] (noting that there is a sign mistake in [ZM,(4)]). Of course, the equations (2.14) are nothing but the residues of the zero curvature equation (2.9) at a m and b n respectively, taking into account the solution (2.13).
2.6. The Zakharov-Mikhailov action. We now substitute the solution (2.13) for A, which we write as L = L η dη + L ξ dξ since it correspond to the Lax connection, into the action (2.6). The 4d Chern-Simons action term becomes where in the second equality we used the fact that L only has components along dη and dξ which implies, in particular, that Tr(dL ∧ L ∧ L) = 0. Using the explicit form (2.13) of L we find Substituting this into (2.15) and performing the integral over dz ∧ dz by using the property (2.11) of the δ-function we find On the other hand, substituting the solution (2.13) for A into the two surface defect contributions to the action, namely (2.7), we obtain Combining together (2.16) and (2.17), and recalling that This coincides with the Zakharov-Mikhailov action [ZM,(10)] up to an overall sign. 2.7. Example. The simplest non-trivial example of the Zakharov-Mikhailov action is obtained by taking N 1 = N 2 = 1. In this case we only have two fields φ 1 and ψ 1 which we denote simply as φ and ψ. Moving to a gauge where V 0 = U 0 = 0, as described in §2.4, and choosing U (0) = −Λ and V (0) = Λ for some fixed constant matrix Λ, the action (2.18) takes the simple form where we have introduced the coupling 2ν := a 1 − b 1 . This action coincides with that of the so-called linear chiral model constructed in [AHP,(3.20)]. The latter can be seen as a generalisation to an arbitrary Lie algebra (here written only for gl N ) of the model proposed by Faddeev and Reshetikhin in [FR] as an ultralocal reduction of the SU (2) principal chiral model. More precisely, the Faddeev-Reshetikhin model is defined by replacing the non-ultralocal Poisson bracket of the SU (2) principal chiral model by an ultralocal one. However, the latter is degenerate and therefore the Faddeev-Reshetikhin model can only reproduce a reduction of the original principal chiral model dynamics, in which the Casimirs of the ultralocal Poisson bracket have been set to constants. In the next section we will derive the covariant Poisson algebra of the Lax connection (2.13) which in the present two-point case generalises the ultralocal algebra for the Lax matrix of the linear chiral model found in [AHP,(3.5)].

Covariant Poisson bracket and r-matrix for the 2d theory
In this section, we will rely heavily on the calculus in the variational bicomplex as presented in [D]. Informally, we introduce two differentials: d is the "horizontal" differential, and acts as the usual exterior differential, while δ is the "vertical" differential that acts only with respect to the fields. We consider (p, q)-differential forms that have a vertical degree p and a horizontal degree q. For instance, L = L η dη + L ξ dξ is a (0, 1)-form, or a horizontal 1-form, and Ω (1) below (3.2) is a (1, 1)-form. For details on how this is used in deriving the r-matrix structure of the covariant Poisson bracket of the Lax connection of a 2d integrable field theory, or more generally an integrable hierarchy, we refer the reader to [CS1,CS3].
We proceed in four steps: to begin with, we derive the multisymplectic form of the theory by considering the variation of its Lagrangian volume form, as established in [D, (19.5.2)] and then used in [CS1]. We can then define the covariant Poisson bracket of certain horizontal forms, called Hamiltonian, using the multisymplectic form. We then show that the Lax form associated with the Zakharov-Mikhailov theory is Hamiltonian and compute its covariant Poisson bracket structureà la Sklyanin [S1, S2], thus exhibiting its r-matrix structure. Finally, we construct the covariant Hamiltonian for the 2d theory, which is the covariant analogue of the usual Hamiltonian obtained by performing the Legendre transformation with respect to both independent variables, and we interpret the zero-curvature equations as covariant Hamilton equations.
3.1. The multisymplectic form. Our starting point is the Lagrangian volume form associated with (2.18), where from now on we shall drop the inessential overall minus sign compared to [ZM,(10)]. However, throughout this section we shall work in the gauge where U 0 = V 0 = 0, so we start from n ψ −1 n . In particular, we have δU m = [δφ m φ −1 m , U m ] and δV n = [δψ n ψ −1 n , V n ]. We also note the identities To show these we need, in particular, to use the fact that δd = −dδ along with the cyclicity of the trace. Combining the above we then find As expected, the first term reproduces the Euler-Lagrange equations in the form (2.14), recalling that we are working in the gauge where U 0 = V 0 = 0. On the other hand, the last term on the right hand side of (3.1) allows us to identify the form which in turn yields the multisymplectic form Ω := δ Ω (1) of the model as (3.3) The multisymplectic form Ω = ω (ξ) ∧ dξ + ω (η) ∧ dη provides the covariant symplectic structure of a field theory. Its coefficients ω (ξ) and ω (η) contain the pull-back to the group of the Kostant-Kirillov forms for the orbits through U where denotes the interior product of a vector field with a form. Let F and G be two Hamiltonian forms . We define their covariant Poisson bracket as follows where q is the horizontal degree of F . Notice that the vector field X F in (3.4) will generally not be unique since Ω may have a non-trivial kernel. Nevertheless, the covariant Poisson bracket (3.5) is seen to be independent of the choice of vector fields X F and X G for both of the Hamiltonian forms F and G. We remark that the covariant Poisson bracket is non-trivial and well defined when the Hamiltonian forms F and G are either both horizontal 1-forms, or one is a horizontal 1-form and the other one is a 0-form (i.e. a function). Our objective is to compute the covariant Poisson bracketà la Sklyanin for the Lax connection L = L η dη + L ξ dξ corresponding to the solution (2.13) for A, in the gauge where U 0 = V 0 = 0. Specifically, let E ij be the canonical basis for gl N and write the Lax connection in this basis as where from now on we shall show the explicit dependence on the spectral parameter. To compute the covariant Poisson brackets between any two components of the Lax connection, we first need to show that these are Hamiltonian 1-forms.
For this we shall need the following useful identities. If M is any GL N -valued field with components M ij , i, j = 1, . . . , N and C is any non-dynamical matrix (meaning δC = 0), then we have In particular, we can use these with M = φ n , C = U Then a direct calculation shows that satisfies δL ij (z) = X ij (z) Ω. Therefore all the components L ij (z) for i, j = 1, . . . , N of the Lax connection are Hamiltonian 1-forms, as required. We shall write the covariant Poisson bracket of the Lax connection L using the standard tensorial notation L 1 := L ⊗ 1 and L 2 := 1 ⊗ L as 3.3. The r-matrix structure. We now turn to the computation of the components on the right hand side of (3.8). We have Introducing the permutation operator P 12 : for any gl N -valued field M with components M ij for i, j = 1, . . . , N , and noting that for any distinct z, w, a ∈ C we have the identity we may rewrite the covariant Poisson bracket (3.8) as In other words, we have the announced result that the Lax connection satisfies the following Poisson algebra with respect to the covariant Poisson bracket {| , |}, where r 12 (z) := P 12 /z is the rational r-matrix. The fact that we have been working with the Lie algebra gl N was convenient for writing the GL N -valued fields φ n and ψ n in components. However, the above derivation can be adapted to hold more generally for any semisimple Lie algebra, working in a basis of the latter.
3.4. The covariant Hamiltonian. Following [D,Lemma 19.5.9], the covariant Hamiltonian related to L ZM is found to be equal to This can be reexpressed directly in terms of the Lax connection as (3.10) This is a rather remarkable formula extending to the present covariant context the familiar formula "H = Tr L 2 " for extracting a Hamiltonian from a Lax matrix in many finite-dimensional integrable systems, such as the Gaudin model or Calogero-Moser system. In fact, the formula (3.10) is very reminiscent of the expression for the Hamiltonian in non-ultralocal integrable field theories described by Gaudin models associated with affine Kac-Moody algebras [V1,DLMV2].
3.5. Flatness of the Lax connection as a covariant Hamilton equation. It was shown for the first time in [CS1] for certain 2d integrable field theories (nonlinear Schrödinger, sine-Gordon, modified Korteweg-de Vries) that the zero curvature equation (2.9) is a covariant Hamilton equation for the Lax connection L associated with the density of the covariant Hamiltonian. By this we mean that, if we define the "covariant flow" of L by in analogy to what one would do in the traditional Hamiltonian formalism, then since we have we can conclude that dL(z) + L(z) ∧ L(z) = 0. The main steps in the derivation of the crucial equality (3.11) are as follows. First, we have by definition Second, we find where we have used the identity (3.6b) in the second equality and (3.9) in the last equality. Substituting the above into (3.12) we obtain (3.11).

Conclusion and outlook
In this paper we derived the Zakharov-Mikhailov 2d action from the 4d Chern-Simons action in the presence of certain surface defects. At the 2d level, the covariant Poisson algebra of the Lax connection was shown to possess a classical r-matrix structure of rational type, thereby recasting the pioneering results of Sklyanin [S1, S2] into a covariant Hamiltonian context. So far, this had only been shown for the sine-Gordon model [CS1] and the entire AKNS hierarchy [CS2,CS3]. There are a number of tantalising questions and possible further directions following this work.
Some of the models (e.g. deformed Gross-Neveu models) considered in the series of papers [By1,By2,By3,ABW] seem to be cousins of the models of Zakharov-Mikhailov type studied here. It would be natural to expect that the covariant Poisson algebra of the Lax connection also holds for these models. Whether this could be achieved by relating them to the present Zakharov-Mikhailov construction is an interesting problem. The extension of the covariant Poisson algebra structure to an entire hierarchy, as obtained in [CS3], is based on the notions of Hamiltonian multiform and multi-time Poisson bracket introduced in [CS2]. In turn, these are based on the idea of Lagrangian multiforms [LN] which provide a generalised variational principle that is able to capture the integrability properties of classical field theories. The Zakharov-Mikhailov action was analysed from this point of view and embedded into a Lagrangian multiform in [SNC]. It is an intriguing problem to understand how such a multiform could effectively arise from a higher dimensional theory, in parallel to the present situation where a single 2d Lagrangian is derived from a 4d one.
The Poisson algebra of the Lax matrix of a non-ultralocal 2d integrable field theory was derived in [V2] by performing a Hamiltonian analysis of the 4d Chern-Simons action for a general 1-form ω. It would be interesting to similarly derive the covariant Poisson algebra of the Lax connection in the present ultralocal setting for which ω = dz. This would involve performing a covariant Hamiltonian analysis of the 4d action (2.6) in order to rederive the covariant Poisson bracket obtained in §3 directly from the 4d Chern-Simons theory.
We showed that the gauge transformations in the Zakharov-Mikhailov action arose as special types of gauge transformations in the 4d Chern-Simons theory for which the gauge transformation parameter g ∈ C ∞ (Σ, G) is independent of the spectral parameter. This is the crudest way of ensuring (2.12) but we believe that a more appropriate condition would be to require that g is (sectionally) holomorphic in order to make a connection with the theory of dressing transformations [ZS]. In other words, it would be interesting to understand if dressing transformations in the 2d integrable field theory can also be understood as arising from gauge transformations in 4d Chern-Simons theory by allowing g ∈ C ∞ (X, G) to depend also on CP 1 as long as the pole structure of the Lax connection (2.13) remains unchanged under such gauge transformations.
In the ultralocal setting considered in the present paper, the regularised 4d Chern-Simons action is easily seen to be gauge invariant. Therefore any defect terms added to the action, as in (2.6), should be gauge invariant themselves. By contrast, in the non-ultralocal setting one needs to impose boundary conditions on the bulk field A at the disorder defects, which are located at the poles of ω [CY]. This is necessary in order to ensure that the action is gauge invariant [BSV]. Alternatively, the gauge invariance can be ensured by introducing new fields living on the surface defects, called the edge modes, and coupling these to the bulk field A [BSV]. It would be very interesting to explore the possibility of combining these two approaches by adding further gauge invariant defect terms to the 4d Chern-Simons action with edge modes. This would have the interesting effect of coupling, in the sense of [DLMV1,DLMV2], ultralocal integrable field theories to a non-ultralocal one.