Off-shell cubic hypermultiplet couplings to N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 2 higher spin gauge superfields

We construct manifestly 4D, N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 2 supersymmetric and gauge invariant off-shell cubic couplings of matter hypermultiplets to the higher integer spin gauge N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 2 multiplets introduced in arXiv:2109.07639. The hypermultiplet is described by an analytic harmonic 4D, N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 2 superfield q+ with the physical component spins s=120\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathbf{s}=\left(\frac{1}{2},0\right) $$\end{document} and an infinite number of auxiliary fields. The cubic coupling constructed has the schematic structure q+ℋ̂s++q+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {q}^{+}{\hat{\mathrm{\mathscr{H}}}}_{(s)}^{++}{q}^{+} $$\end{document}, where ℋ̂s++\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\hat{\mathrm{\mathscr{H}}}}_{(s)}^{++} $$\end{document} is a differential analytic operator of the highest degree (s− 1) accommodating the massless gauge N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 2 multiplet with the highest spin s. For odd s the gauge group generators and couplings are proportional to U(1)PG generator of the internal SU(2)PG symmetry of the hypermultiplet and so do not exist if SU(2)PG is unbroken. If this U(1)PG is identified with the central charge of 4D, N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 2 supersymmetry, a mass for the hypermultiplet is generated and the odd s couplings vanish in the proper massless limit. For even s the higher-spin gauge transformations and cubic superfield couplings can be defined for both massive and massless (central-charge neutral) hypermultiplets without including U(1)PG generator. All these features directly extend to the case of n hypermultiplets with the maximal internal symmetry USp(2n) × SU(2).


Introduction
The theory of higher spin fields attracts vast attention due to its intimate relationships with both quantum field theory and (super)string theory. In particular, it was conjectured [1] that the superstring theory describes the "massive" phase of the massless higher spin theory corresponding to the spontaneous breakdown of the infinite-rank gauge symmetries of the latter, with masses of the higher-spin excitations being generated through some generalized Higgs-type mechanism. Though the precise role of the higher-spin fields in fundamental interactions is still unclear, it is expected that the relevant particles could dominate the string interactions at high energies (see, e.g., [2]). Like in other instances, the study of supersymmetric versions of various higher-spin models is of pivotal importance, as supersymmetry can shed more light on the hidden structures and properties of the theory. Since the higher-spin theory, like the superstring theory, predicts infinite towers of fermions and bosons related by supersymmetry, there arises the natural problem of describing these JHEP05(2022)104 sets in the appropriate field-theory language, where fields are combined into supermultiplets and further into superfields. Thus, the issue of describing supersymmetric theories of higher spins in the superfield approach acquires high actuality.
The superfield formulations have a great privilege to provide an explicit off-shell realization of supersymmetry. They are most preferable for constructing invariant actions since the transformations of given supersymmetry do not depend on the choice of a specific model. In four dimensions, there were developed the exhausting superfield formulations of N = 1 theories in general and chiral 4D, N = 1 superspaces (see, e.g., [3,4]), as well as of N = 2 theories in the harmonic superspace and its analytic subspace with half the original Grassmann coordinates [5][6][7].
In a recent paper [8] we have constructed the completely off-shell manifestly N = 2 supersymmetric superfield extension of arbitrary 4D integer-spin free massless theory. The construction was based on the harmonic superspace method [5][6][7], which so confirmed its universal role as the most adequate and convenient approach to 4D, N = 2 supersymmetric field theories. The theory constructed yields just N = 2 off-shell supersymmetrization of the free massless higher-spin field theory pioneered by Fronsdal [9] and Fang and Fronsdal [10].
The natural next step is to set up interacting extensions of the free theory. The construction of interactions in the theory of higher spins is a very important, but highly difficult task. The only fully nonlinear constructions known at the moment are limited to the Vasiliev's theory in (A)dS d [11][12][13] (for a review see [14][15][16][17][18]) and the 3D Chern-Simons higher spin gravity formulations [19][20][21][22] (for a review see [23]). As for the interacting higher spin theories in the Lagrangian approach, one of the possible obstructions here is provided by no-go theorems imposing significant restrictions on the corresponding interactions (for a review and references see, e.g., [24]). The statement of the no-go theorems relevant to our consideration is that the constraints on the low-energy scattering in flat space-time seemingly forbid massless particles with spins s > 2 to participate in interactions. However, in many cases no-go theorems may be evaded by relaxing various original assertions.
Let us turn back to 4D, N = 2 supersymmetric theory. The most general lower spin 4D, N = 2 matter multiplet is hypermultiplet with the physical component spins s = ( 1 2 , 0). The off-shell hypermultiplet requires an infinite number of auxiliary fields and admits an unconstrained formulation only in N = 2 harmonic superspace [5]. So if one intends to study off-shell N = 2 cubic vertices, it is unavoidable to use the harmonic superspace approach.
In this paper we construct, for the first time, the off-shell manifestly N = 2 supersymmetric cubic couplings 1 2 − 1 2 − s of an arbitrary integer higher spin s gauge N = 2 multiplet to the hypermultiplet matter in 4D, N = 2 harmonic superspace. 2 The couplings are linear in the analytic gauge superfields, bilinear in the analytic hypermultiplet q + superfields and are written as integrals over the analytic harmonic superspace. An interesting new peculiarity is the necessity of breaking the Pauli-Gürsey SU(2) PG symmetry of the free hypermultiplet action down to U(1) for odd spins s and the explicit presence of this U(1) generator in the gauge transformations and the relevant superfield cubic actions. If this U(1) generator is identified with the central charge of rigid 4D, N = 2 supersymmetry, it induces a mass for the hypermultiplet and only even spin couplings survive in the massless limit. As regards the higher-spin gauge superfields, we deal with the same linearized actions for them as in ref. [8].
The cubic couplings of the hypermultiplet to N = 2 higher spin gauge superfield constructed in the present paper are based on the gauge principle: the interactions are resulted from localization of global symmetry. We begin with the free hypermultiplet theory in harmonic superspace and find for this theory an infinite set of higher-derivative global symmetry transformations. The localization of these transformations in the harmonic analytic superspace requires introducing the corresponding higher spin gauge harmonic superfields with the gauge group already described in [8]. As a result, we arrive at the manifestly N = 2 supersymmetric cubic interactions of the hypermultiplet with the higher spin gauge superfields, in the first order in the latter. 3 The paper is organized as follows. Section 2 is devoted to a brief description of a completely off-shell invariant action for the free massless N = 2 gauge theory with an arbitrary maximal integer spin s of the supermultiplet. In section 2.2 we introduce the analytic N = 2 invariant differential superfield operators associated with the higher-derivative JHEP05(2022)104 hypermultiplet gauge transformations. Just these operators are the basic ingredients of our construction. In section 3 we present the detailed derivation of the manifest form of 1 2 − 1 2 − s couplings of N = 2 higher spins to the hypermultiplet for s = 1, 2, 3, 4, as well as the corresponding global and local transformations of the hypermultiplet. A generalization to arbitrary even and odd N = 2 spins is the subject of section 4. In section 5 we summarize the results, briefly describe a generalization to an arbitrary number of hypermultiplets (with the internal symmetry USp(2n)) and discuss possible directions of further development of the material presented.

JHEP05(2022)104
An important ingredient of the HSS formalism is the harmonic derivatives D ++ and D −− which have the following form in the analytic basis Here we used the standard notations for the partial derivatives with respect to harmonic variables: The crucial difference between the derivatives D ++ and D −− is that D ++ preserves analyticity, while D −− does not. All superfields except for the hypermultiplet are assumed to be x 5 -independent, while the action of ∂ 5 on the hypermultiplet is identified with the generator of some U(1) isometry of the free hypermultiplet action (see section 3).
For what follows it will be useful to present the "passive" transformations of various partial derivatives under N = 2 supersymmetry: (2.11)

Analytic prepotentials and invariant action
In this section we summarize the main results of [8]. We present all the necessary information about description of N = 2 higher spin supermultiplets in harmonic superspace. Detailed examples of spin 2 and spin 3 description were given in [8], here we will at once deal with the general structure of arbitrary spin s theory. The N = 2 gauge multiplet of highest integer spin s (s ≥ 2) 5 is accommodated by two real analytic bosonic superfields and two conjugated analytic spinor superfields These transformations can be used to choose the appropriate WZ gauge and then show that the physical multiplet involves spins s, s − 1 2 , s − 1 2 , s − 1 : h ++α(s−1)α(s−2)+ = h ++α(s−1)α(s−2)+ . (2.14) Here the fields (2.15) are physical and describe the bosonic spin s, a doublet of the fermionic s − 1 2 spin fields and the bosonic spin s − 1. All other fields are auxiliary.
For discussion of the residual gauge transformations and further features of the component reduction we refer the reader to [8], where all details in the s = 2 and s = 3 cases were presented. It was shown there that the residual gauge transformations completely coincide with the corresponding transformations for the physical bosonic [9] and fermionic fields [10].
To construct the invariant linearized actions, one must introduce the negatively charged potentials (2.16)

JHEP05(2022)104
They are defined as solutions of the harmonic flatness conditions: These extra non-analytic potentials will not be present in the Lagrangian of cubic interactions with hypermultiplet and are needed only for constructing the invariant actions of the gauge superfields. The bosonic analytic prepotentials (2.12a), (2.12b) and the negatively charged potentials (2.16) have non-standard transformation laws under N = 2 global supersymmetry: 7 Spinor potentials have the standard N = 2 superfield passive transformation rules: 8 As the building blocks of N = 2 supersymmetric actions, one defines non-analytic N = 2 singlet superfields, It is straightforward to be convinced that indeed δ G ±±··· = 0. The newly introduced superfields satisfy the harmonic flatness conditions and transform under the gauge group as

JHEP05(2022)104
where Λ's are composed out of the analytic gauge parameters λ's and the coordinates θ −α ' The N = 2 gauge invariant action, up to a normalization factor, has the universal form for any s: where κ s is a higher spin analog of Newton's constant. The N = 2 supersymmetry of (2.22) is manifest, while the gauge invariance can be checked by bringing the gauge variation of S (s) to the form and proceeding further by taking into account that the ++-potentials in (2.19) are linear in the negatively charged coordinates with the analytic coefficients. The straightforward calculation immediately yields that δ λ S (s) = 0. The dimension of the coupling constant κ s is fixed by the dimension of the analytic prepotentials. In mass units, the "engineering" dimension of the analytic gauge potentials is can . In what follows, it will be more convenient to deal with the gauge superfields of engineering dimensions.
Summarizing, here we have repeated the main results of [8]. We presented the transformation laws of the spin s superfields under both the gauge group (2.13) and N = 2 global supersymmetry (2.18a), (2.18b) and then wrote down the invariant action (2.22) which is valid for any s.

Gauge transformations via differential operators
In this section, we define analytic, N = 2 supersymmetry-preserving differential operators which encode all the analytic prepotentials of the spin s higher spin supermultiplet. We will start with the case of N = 2 Einstein supergravity (corresponding to s = 2) and then generalize this construction to an arbitrary higher spin N = 2 supermultiplet. These operators will play the key role in the construction of cubic couplings of hypermultiplet to higher spin gauge multiplets.

Differential operators for gauge transformations in N = 2 Einstein supergravity
One of the underlying principles of Einstein N = 2 supergravity in harmonic superspace [5] is the covariantization of flat harmonic derivative (2.7) Analytic superfields h ++αα , h ++μ+ and h ++5 are unconstrained analytical prepotentials of N = 2 Einstein supergravity. One can easily check that the operatorĤ ++ (2) is invariant under N = 2 global supersymmetry (2.11) accompanied by (2.18a) and (2.18b) specialized to the case s = 2. Note that the operatorĤ ++ (2) and its arbitrary spin s generalizationŝ H ++ (s) are dimensionless and the engineering dimensions of the component gauge superfields are uniquely fixed by the dimensions of the related (higher-order) derivatives.
Fundamental gauge group of Einstein N = 2 supergravity (its "minimal" version [77,78]) is the following group of superdiffeomorphisms of the analytic harmonic superspace: (2.26) It will be useful for future application to realize these transformations on the gauge superfields as the "active" ones using the differential operator: Hereafter we use the notations M = (αα,μ+, 5). By construction, the full covariant harmonic derivative is invariant under the action of gauge group of supergravity: From (2.28) one can easily deduce the active transformation laws for the analytic prepotentials (2.13) The linearized form of the transformation law (2.29) reads: So we conclude that the analyticity preserving N = 2 supersymmetry invariant differential operatorĤ ++ (2) has the very simple linearized transformation law realized as the commutator of the flat harmonic derivative with the differential operator involving all the superfield gauge parameters which are contracted with the corresponding partial derivatives. In the case of spins s ≥ 3 there are no realizations of the gauge groups on the coordinates, so just the active form of gauge transformations proves to be relevant.

Differential operators for gauge transformations in N = 2 spin s higher spin theory
Now we will extend s = 2 construction of N = 2 supersymmetry invariant analytic differential operatorĤ ++ (2) (2.24) to general higher N = 2 spins s. It will be the necessary ingredient for construction of interactions of hypermultiplet with higher spins.
First of all, we will define the analytic N = 2 invariant differential operator that is a natural extension of the differential operator (2.24) for s = 2: (2.32) Here s ≥ 2 and it is assumed that the corresponding spinor indices in h ++α(s−2)α(s−2)M are always symmetrized with those in M , e.g., The corresponding operator involving the analytic gauge parameters is constructed as a natural generalization of (2.27): with the same convention concerning the spinor indices hidden in the index M , e.g., . In these notations, the gauge transformations (2.13) are realized as δ λĤ Note, that here, in contrast to the spin 2 case (2.31), we added s − 2 vector derivatives ∂ s−2 α(s−2)α(s−2) in (2.34). These extra derivatives are necessary both for reproducing the correct transformation properties of the analytic prepotentials from (2.34) and for ensuring the compatibility with the rigid N = 2 supersymmetry. Indeed, it is easy to check that only the operatorĤ ++ (s) is N = 2 supersymmetric, δ Ĥ ++ (s) = 0. For s ≥ 3, neither the operator (2.32) nor any its contraction with ∂ αα , such that some amount of the free spinor indices is left, are invariant under N = 2 supersymmetry (2.11), (2.18).
One can also construct a useful N = 2 superfield: where we used the notation where The superfield Γ ++ (s) is also invariant under N = 2 supersymmetry (2.11) and (2.18): We will use the operator (2.35) and the superfield (2.37) for construction of gaugeinvariant couplings with the hypermultiplet. Note that there are only one N = 2 supersymmetric invariant operatorĤ ++ (s) and only one invariant superfield Γ ++ (s) which are linear in the analytic prepotentials. This will significantly restrict the possible cubic couplings with the hypermultiplet.

Hypermultiplet in harmonic superspace
The N = 2 hypermultiplet free action has the form [5][6][7] We can write the action in the two equivalent forms: in terms of the pseudo-real analytic superfield q +a (ζ) with the Pauli-Gürsey SU(2) PG doublet indices (a = 1, 2) or in terms of the complex superfields q + andq + . These two representations are related by: Hereafter, we will use the first form of the action since the manifest SU(2) PG symmetry crucially simplifies the calculations. A peculiarity of our theory is the presence of the derivative ∂ 5 and of the corresponding gauge superfields. The option ∂ 5 q +a = 0 is most general (though the choice ∂ 5 q +a = 0 is also admissible). In order to avoid any integration over x 5 , we impose the standard Scherk-Schwarz condition that ∂ 5 q +a coincides (up to a phase factor) with the action of some U(1) PG ⊂ SU(2) PG . Without loss of generality we assume 9 We use the standard definition of the analytic superspace integration measure [5]:

JHEP05(2022)104
The parameter m is a mass of the hypermultiplet, so we deal with the massive hypermultiplet in the general case. It is easy to check that The action (3.1) is invariant under the rigid N = 2 supersymmetry, where we need to take into account also the transformation of x 5 (2.6) for m = 0. Because of (3.5), it does not affect the N = 2 supersymmetry transformation of the Lagrangian, though modifies N = 2 transformation of q +a . The internal symmetry of the single hypermultiplet action is SU(2) PG × SU(2) aut in the massless case (∂ 5 q +a = 0) and U(1) PG × SU(2) aut in the massive case (∂ 5 q +a = 0). The equation of motion for the free massive hypermultiplet is: On shell, the analytic harmonic superfield q +a is reduced to: Eq. (3.7) also implies the massive equations of motion for the physical fields. We will not need their explicit form. The massless hypermultiplet corresponds to setting m = 0 in (3.3) and (3.8).
In the sequel we will derive the coupling of the hypermultiplet to the gauge supermultiplets of integer spins described in the HSS approach by analytic superfields (2.12a) and (2.12b).

Guiding principles
Before turning to our basic subject, we will formulate a few generic á priori restrictions on the structure of possible interactions: • Analyticity. We require that coupling must be analytic because hypermultiplet is described by an analytic superfield and the well known couplings of the hypermultiplet to the gauge and supergravity N = 2 theories preserve the analyticity. Moreover, the harmonic analyticity is one of the major and crucial features of the harmonic superspace approach and it proved its power while having constructed N = 2 Fronsdal theory in the harmonic superspace [8].
• N = 2 supersymmetry. The manifest rigid N = 2 supersymmetry is a necessary ingredient of our construction and it is one of the main general reasons for employing harmonic superspace. Since the hypermultiplet has the N = 2 transformation law (3.6) the gauge higher-spin superfields must also appear in a way preserving N = 2 supersymmetry. This requirement significantly limits the possible form of cubic interaction. Using the results of section 2.2.2 the gauge superfields can appear only JHEP05(2022)104 within two possible terms, the first one involving the supersymmetry-invariant differential operator (2.35), and the second one involving the supersymmetry-invariant superfield (2.37): Here A ab and B ab = B ba are some matrices, which in general can break the Pauli-Gürsey SU(2) PG symmetry. The possible structure of these matrices is fully specified by the gauge symmetry (which, in turn, is determined by the appropriate global symmetry of the q + action through the gauging procedure).
• Gauge invariance. The most crucial property is gauge invariance in the leading order in the gauge superfield.
As usual, the gauge transformations should be defined through gauging of the global symmetry transformations and so should yield global symmetry upon restricting gauge parameters to constant values.
We know the linearized gauge transformation laws of the basic quantitiesĤ ++ (s) (2.34) and Γ ++ (s) (2.39). So in the leading order in the gauge parameters, the gauge variation of the full action should be composed of the two possible terms: These variations must be canceled by the appropriate variations of the free hypermultiplet action induced by the higher-spin gauge transformations of the hypermultiplet.
Thus, the general strategy for building cubic interactions of the hypermultiplet with the given spin s supermultiplet should be the following. As the starting point, one finds out the appropriate global symmetry of the free hypermultiplet action realized on the superfield q +a . Secondly, one considers the most general gauging of this global symmetry and singles out the necessary combinations of the gauge transformations of the hypermultiplet ensuring cancelation of the gauge superfield variations (3.10a) and (3.10b). Finally, from the last two steps, one determines the matrices A ab and B ab , and, as a result, derives the sought cubic couplings.
It is worth emphasizing that the general gauging procedure described above makes it possible to gauge all global symmetries realized on the hypermultiplet. Thus it is capable to yield all admissible local first-order Noether couplings of the hypermultiplet to higher-spin gauge fields (see the relevant discussion also in section 5).
Here we adapt it to the generic form applicable to higher spins.

JHEP05(2022)104
Rigid symmetry. The free hypermultiplet action (3.1) is invariant under the U(1) rigid transformation which can be realized as a "shift" with respect to x 5 with the generator (3.3). Now we will gauge this global symmetry.
Gauging. The spin 1 gauge transformations read: where λ 5 (ζ) is an arbitrary analytic gauge parameter. The formalism of section 2.2.2 is adapted to this degenerate case as The transformation of the hypermultiplet action (3.1) under (3.11) reads: If D ++ λ 5 = 0, the gauge transformation (3.11) becomes a rigid internal symmetry of the action, in the full agreement with gauging procedure. In order to make the action (3.1) invariant under gauge transformations (3.11), one needs just the compensating gauge superfield h ++5 (ζ) which is introduced through the substitution (3.12) Using the gauge transformations δh ++5 = D ++ λ 5 , one can choose the Wess-Zumino gauge for h ++5 as: which yields just the off-shell field content of N = 2 Maxwell multiplet. The engineering dimension of h ++5 is −1 and the passing to the gauge superfield h ++ with the zero canonical dimension is accomplished as h ++ = mh ++5 . The gauge invariant and N = 2 supersymmetric action of the spin 1 gauge superfield has the standard form, Note that the action (3.15) is exactly invariant under the gauge transformations (3.11) and (3.13) in consequence of the relation In the case of higher spin couplings (s > 2), no such an exact gauge invariance is present in the leading order in gauge superfields in view of absence of the analogous relations. So it is the peculiarity of the spin 1 coupling only.

JHEP05(2022)104
One more notable feature of the simplest case considered is the following. The relations (3.11) and (3.15) imply that if the hypermultiplet is x 5 independent, the mass parameter m vanishes in (3.15) and a non-trivial cubic coupling of such a massless hypermultiplet to spin 1 gauge field is seemingly impossible. To avoid this obstacle and get a possibility to construct the cubic coupling for both massive and massless theory, one can resort to the following reasoning. The main point is that in general it is not obligatory to identify the central charge ∂ 5 with the generator J of U(1) PG ⊂ SU(2) PG to be gauged. The basic condition for such a generator is whence, without loss of generality, This corresponds just to replacing ∂ 5 to J in the gauge transformation law (3.11): Hence the hypermultiplet can still be chosen ∂ 5 -neutral (∂ 5 q +a = 0) and thus massless. The interaction with Maxwell multiplet will have the same form as in (3.15), where we have to replace m by 1 and h ++5 by the dimensionless h ++ . Note that in the massive case, with ∂ 5 q +a ∼ i(τ 3 ) a b q +b = 0, it is necessary to identify ∂ 5 with the gauged U(1) PG generator J (up to an unessential constant of mass dimension, ∂ 5 ∼ mJ), in order to preserve the gauge invariance of the total q +a action. In the massless case J still commutes with all N = 2 supersymmetry generators (as SU(2) PG generators do) but has no any relation with the central charge.
Let us point out once more that the N = 2 Maxwell theory analytic potential and its coupling to the charged hypermultiplet can be constructed without any reference to the fictitious coordinate x 5 and its interplay with the generator J. The way we have followed here is just a simple illustration of the general procedure outlined in the previous subsection. As we will see soon (in the subsections 3.5 and 4.1), the above procedure with the extra coordinate x 5 and the matrix generator J perfectly well works for all odd spins s ≥ 3 .
We finish this subsection with the two comments. First, in the gauge-covariant derivative (3.12) one can remove the mass-generating background in D ++ just by redefining 10 This means that, in the gauged s = 1 theory for the single hypermultiplet, there is no actual difference between massive and massless cases: the mass parameter m appears only as a coupling constant in the minimal interaction (3.15) and can be removed (or set equal to 1) by the proper rescaling of h ++5 (the same phenomenon can be easily traced also in the component formulation).
Secondly, perhaps a more convincing explanation why it is useful to introduce the matrix generator J ⊂ su(2) PG besides ∂ 5 is as follows. As we saw in the previous section, JHEP05(2022)104 the differential operators relevant to the spin s case are always of the order (s − 1), so it would be more natural to describe the spin 1 by the differential operators of the zeroth order, than of the first order as in (3.13). Such a description, in accord with the general scheme of section 2, amounts to The flat harmonic derivative D ++ is still defined by the general expression (2.7), with the central charge operator ∂ 5 . The relevant gauge transformation of q +a is now postulated as in (3.21), with the SU(2) PG symmetry matrix generator J properly realized on q +a . The covariantization of (3.22), when it acts on q +a , is then accomplished as The standard condition [D ++ (1) , J]q +a = 0 still admits the solution ∂ 5 q +a = 0 , which yields a massless hypermultiplet. For ∂ 5 q +a = 0 , the same condition necessitates the relation ∂ 5 ∼ J, yielding the massive hypermultiplet. As we argued above, for the spin 1 coupled to the hypermultiplet, the massive and the massless Lagrangians are related through a redefinition of the gauge superfield h ++ . However, no such an equivalence between massive and massless Lagrangians is valid for higher odd spins s ≥ 3. It is just the description outlined here that directly extends to the s ≥ 3 case.
Here we repeat the relevant construction in the form most convenient for further generalizations.
Rigid symmetry. The free hypermultiplet possesses the following global supetranslational symmetry It involves five constant bosonic parameters λ αα , λ 5 , four constant spinor parameters λα i , such that λ ±α = λα i u ± i , and it can be treated as a copy of the rigid N = 2 supersymmetry transformations in their active form. However, we will gauge just it, leaving N = 2 supersymmetry still rigid, so that the latter forms a semi-direct product with the gauge extension of (3.25). Recall that in the previous subsection we have already introduced λ 5 transformations in order to describe spin 1 supermultiplet after their gauging. The symmetry (3.25) is an extension of this ∂ 5 symmetry, such that its gauging generates the multiplet of minimal N = 2 Einstein supergravity.

JHEP05(2022)104
OperatorΛ rig commutes with the harmonic derivative (we assume, as earlier, that all operators act on analytic superfields): The Lagrangian (3.1) is invariant, up to total derivative, under these transformations: (3.28) Gauging. In this case there are two possibilities for gauge transformations of the hypermultiplet: The relevant first-and zeroth-order differential operators are the particular s = 2 case of the general operators (2.33) and (2.40). Here λ αα (ζ), λ +α (ζ) and λ 5 (ζ) are arbitrary analytic gauge parameters. The first type of transformations, δ 1 q +a , corresponds to the direct gauging of the above supertranslations, while the second type δ 2 q +a can be treated as a special gauging of the constant-parameter rescaling of the hypermultiplet (which is not invariance on its own). Under (3.29) the hypermultiplet action (3.1), up to a total derivative, transforms as: The transformation (3.30) leads to the following variation of the Lagrangian: where we made use of the evident property q +a (D ++ Ω (2) )q + a = 0 . (3.33) Thus the total gauge variation of the free Lagrangian reads: This variation vanishes for the constant parameters, so the relevant transformations provide the evident symmetry of the action.
To couple spin 2 theory to the hypermultiplet we use the differential operator (2.24), Respectively, the gauging is accomplished as: The variation ofĤ ++ (2) in the second term cancels (3.34). So the action (3.37) provides the gauge invariant coupling of the spin 2 to the hypermultiplet in the leading (first) order in gauge superfields. Now we could try to restore the next orders. The δ 1 and δ 2 variations of the hypermultiplet q +a in the second term of (3.37), up to a total derivative, generate the following new term:δ 1 So, for the full variation of the gauge superfield part we obtain: Thus, in contrast to the spin 1 case, where we had the exact gauge invariance without any modification of the transformation law of the relevant operatorĤ ++ (1) , in the spin 2 case one is led to modify the gauge transformation ofĤ ++ (2) to achieve the complete gauge invariance.
The deformation (3.41) of the spin 2 transformation law means that the linearized action of N = 2 supergravity (2.23) must also be modified by including nonlinear terms in the action, so as to achieve the full gauge invariance for both the pure supergravity gauge superfield action and the hypermultiplet couplings. It is an essential difference from the spin 1 case where the pure gauge part of the spin 1 total action is given by the action (3.17) and does not require any correction terms. A complete nonlinear harmonic superfield action for N = 2 Einstein supergravity has been constructed in [79]. We emphasize that the discussion of any nonlinear aspects of this kind are beyond the scope of the present paper.
To summarize, we have restored, from the transformation law for hypermultiplet and gauge fields, the well known complete coupling of the hypermultiplet to the minimal N = 2 supergravity multiplet:

JHEP05(2022)104
It is worth noting that it is easier to check the gauge invariance of the analytic superspace action with the Lagrangian (3.42), using the passive form of the gauge transformations Our last remark is that it is consistent to choose ∂ 5 q +a = 0 in the spin s = 2 case. Under this restriction, the prepotential h ++5 drops out from the hypermultiplet action which so corresponds to the massless hypermultiplet. However, the h ++5 term should still be present in the gauge action (2.22) with s = 2 for ensuring the gauge invariance of the latter.

Spin 3 coupling
In this case we will also start with the appropriate rigid transformations. Our previous consideration suggests that for any s the q +a variation of the highest order in derivatives involves the same differential operator as in (2.35) Here we used the notation (2.33) for the differential operatorΛ αα in the s = 3 case. The bosonic parameters λ (αβ)(αβ) , λ αα , as well as the coefficients of the harmonics in the fermionic parameters λ (αβ)α± = λ (αβ)αi u ± i ,λ (αβ)α± =λ (αβ)αi u ± i , are assumed to be coordinate-independent. So we face 9 + 4 = 13 bosonic parameters and 2 · 3 · 2 = 12 fermionic parameters. Assuming that the operators act on the analytic superfields, we check that: Using this property and integration by parts one can derive the relevant variation of the free Lagrangian: 11 (3.47) We observe that the variation of the Lagrangian is not reduced to a total derivative, which means that the transformations (3.44) do not constitute a global symmetry of the JHEP05(2022)104 free hypermultiplet. Thus we need to look for some alternative definition of the rigid transformations of the hypermultiplet in the spin s = 3 case.
Option 2. This is a generalization of the spin 1 rigid transformations (3.3). We take the rigid transformations in the form: As in the spin s = 1 case, we have introduced here a new generator Jq +a which á priori is not obliged to be collinear to ∂ 5 appearing in D ++ . To substantiate the choice (3.48), we note that the relevant modification of the important condition (3.46), when applied to the hypermultiplet superfield q +a , implies, apart from (3.46), also the condition Assuming that J is identified with some U(1) PG ⊂ SU(2) PG , like ∂ 5 , eq. (3.50) leaves us with only two possibilities: the second option being related to the fact that no mutually commuting two generators can be found in SU (2). Without loss of generality, Jq +a in (3.51a) can be chosen the same as in (3.51b). The variation of the free Lagrangian under (3.48) is vanishing, Gauging. The gauge transformation operators in this case arê and all gauge parameters are unconstrained analytic superfields. There are few admissible transformations laws, which generalize global symmetry (3.48). It is convenient to choose the following ones as independent:

JHEP05(2022)104
We vary L +4 free by the transformations (3.55)-(3.58) with arbitrary coefficients (up to a common rescaling), using the integration by parts and the relations like In particular, up to a total derivative, Then, properly fixing the numerical coefficients, we single out two appropriate combinations of the variations (3.55)-(3.58), (3.60) such that they have the necessary form (3.10a), (3.10b) for s = 3: whereĤ as well as the transformation law for the superfield Γ ++ (3) , (3.67)

JHEP05(2022)104
It will be also useful for the future generalizations to rewrite the transformation laws δ λ q +a and δ ξ q +a as:  12 In the present paper we limit ourselves to the invariances only in the leading order in gauge superfields.
Last but not least. As we saw, the consistent minimal coupling of the hypermultiplet to the gauge N = 2 spin s = 3 superfields is possible only provided Jq +a = i(τ 3 ) a b q +b = 0, which implies that this coupling necessarily breaks SU(2) PG down to U(1) PG , the generator of which is further identified with J. In accord with (3.51), the hypermultiplet can still stay massless [option (a), with ∂ 5 q +a = 0], or massive [option (b)], such that the operator J is proportional to the central charge ∂ 5 which is not vanishing in the second case. These two possibilities are essentially different because setting ∂ 5 q +a = 0 in the second case not only makes the hypermultiplet massless but also eliminates all its couplings to N = 2 spin s = 3 superfields. As distinct from the simplest s = 1 case, in the s = 3 case it is impossible to relate the massless and massive hypermultiplet Lagrangians by any redefinition of the gauge N = 2 potentials. The same features are characteristic of all N = 2 odd spins (see below). On the contrary, the hypermultiplet couplings to the higher spin N = 2 spin multiplets with even s can be defined for both massive (∂ 5 q +a = 0) and massless (∂ 5 q +a = 0) cases on equal footing, without insertions of the SU(2) PG generators in the transformation laws and cubic superfield coupling. Note that the full symmetry of the odd spin s case is always U(1) PG ×SU(2) aut , as distinct from the maximal symmetry SU(2) PG ×SU(2) aut of the even spin case (which can be attained for ∂ 5 q +a = 0).

Spin 4 coupling
Rigid symmetry. Similarly to other cases, the rigid symmetry to be gauged is given by the maximal-degree differential operator: The variation of the action L +4 free under the transformation (3.70) is easily checked to vanish (modulo total derivatives), So we have picked up the appropriate symmetry of the free hypermultiplet action. Now we will gauge it.
Gauging. As before, we promote the constant parameters to arbitrary analytic superfields. One can construct six independent transformation laws with the parameters λ M and derivatives thereof. It will be convenient to choose as a basis the following gauge JHEP05(2022)104 variations: Here we used the general definition (2.33) forΛ αβαβ and (2.40) for Ω αβαβ . The derivatives act freely on all objects to the right, in accord with our previous conventions. The crucial difference from the spin s = 3 case is the absence of the derivative ∂ 5 (or the generator J) in these transformation laws. As a result, they are non-trivial for both the ∂ 5 q +a = 0 and the ∂ 5 q +a = 0 cases. In fact, they are a generalization of the spin s = 2 gauge transformations. The gauge variations of L +4 free are performed straightforwardly, integrating by parts at the intermediate steps. We present only the final answers These variations still admit, through integration by parts, some other forms sometimes more convenient for calculations, e.g.,

JHEP05(2022)104
(the variation (3.82) is reproduced by taking off all derivatives from D ++ q + a in (3.86), except for the first term).
Invariant Lagrangian. Using the formulas for variations and summing them with undetermined coefficients, one can single out their unique combination which can be canceled by the gauge transformation of N = 2 invariant analytic differential operatorĤ ++(αβ)(αβ) , constructed out of the analytic gauge prepotentials: (3.88) The factor 1 2 was introduced for ensuring that δ λ q +a is reduced to (3.70) in the rigid limit. Thus the Lagrangian describing the coupling of spin 4 gauge supermultiplet to the hypermultiplet, in the leading order in the former, reads: The Lagrangian L +4 gauge is gauge invariant (in the leading order) and completely N = 2 supersymmetric. In contrast to the spin s = 3 case, the option of adding terms with Γ ++ (4) is absent in the s = 4 theory.
For further generalizations, it is useful to rewrite δ λ q +a as: (3.90) or, in a more explicit form, We stress once more that the Lagrangian (3.89), like its spin 3 counterpart (3.67), is invariant under gauge transformation only in the leading order in the gauge prepotentials. As distinct from the spin 2 case, no clear way is seen how to modify the linearized transformation laws (2.13) so as to restore the full spin s = 4 gauge invariance of the hypermultiplet coupling (if exists).

General odd spins
In this section, we will assume that s is odd integer. As before, we start with rigid transformations: and Here all the parameters are coordinate-independent. We have s 2 + (s − 1) 2 bosonic parameters and 2s(s − 1) fermionic ones. There is valid the property: The variation of the free hypermultiplet Lagrangian (modulo total derivatives) is zero due to the symmetry of (τ 3 ) ab : To check this, one needs to integrate s − 1 times by parts in the second term.
Gauging. Like in the spin 3 case, one could start with the most general combination of gauge transformation and then follow the strategy of section 3.2. But due to the strong constraints (3.9) on the form of interaction, it will be sufficient to guess, from the very beginning, the gauge transformations of hypermultiplet which lead to the gauge-invariant coupling with higher spins. The interactions obtained in this way will be unique. Fortunately, the sought gauge transformations prove to be a direct generalization of those we met while considered the spin 3 and spin 4 examples.
In the generic case of the odd spin supermultiplet, the proper generalization of the spin 3 gauge transformations (3.68) is as follows The variation of the hypermultiplet Lagrangian can be calculated in full analogy with (3.81): α(s−2)α(s−2) q +a q +b . (4.7)

JHEP05(2022)104
Now we are led to find a transformation that would cancel the first term in (4.7). It is a direct generalization of the second term in (3.68): The relevant variation of the free hypermultiplet action, in full analogy with the spin 3 case, is (4.9) Using (4.9), one can cancel the first and last terms in (4.7). Collecting all terms, we obtain for the full variation: (4.10) Using the transformation law (2.34), one can cancel these terms by passing to the gauge superfield-modified Lagrangian: The relevant action is a direct generalization of the corresponding action for s = 3 defined by the Lagrangian (3.67). The Lagrangian (4.11) is also invariant only up to the leading order in gauge prepotentials.
Let us now turn to the ξ-transformations (3.69). A direct generalization of the latter is The relevant variation of the free Lagrangian has the form: and it can be easily canceled by adding Γ ++ (s) (2.37), so the most general Lagrangian reads: (4.14) Here ξ is an arbitrary real parameter. To summarize, for an arbitrary odd N = 2 spin s the gauge transformations of the hypermultiplet are given by δ λ,2 q +a and δ (s) ξ , with the variations being defined in (4.6), (4.8) and (4.12). The action which is gauge invariant up to the first order in gauge superfields is given by eq. (4.14). The hypermultiplet can be massless (∂ 5 q +a = 0 , Jq +a = 0) or massive (∂ 5 q +a = 0 , Jq +a = 0 , ∂ 5 q +a = mJq +a ).

General even spins
In this section, we will generalize the spin 4 transformation laws (3.90) to an arbitrary integer even spin s.
Let us start with the definition of rigid symmetry of the free hypermultplet for even higher spins: δ (s) The group-parameter structure of coincides with that already given in (4.3), the crucial difference being the absence of the generator J (and, respectively, of the matrix (τ 3 ) a b ) in the transformation law (4.15). It is a direct generalization of the rigid symmetry (3.70) pertinent to the spin s = 4. Here we face s 2 +(s−1) 2 constant bosonic parameters and 2s(s−1) constant fermionic parameters. One can check that this transformation, in the complete analogy with (3.73), indeed provides a symmetry of the free hypermultiplet action. Now we will gauge this symmetry, using the even higher-spin s gauge supermultiplet.
Gauging. Like in the previous cases, one can define many local generalizations of the rigid transformation (4.15), with the proper analytic superfield parameters and derivatives thereof. However, as a consequence of the strong restrictions (3.9) on the possible structure of the interaction, it is enough to explicitly guess some kind of implementation of the gauge group, such that it meets the restrictions just mentioned. The interaction constructed in this way will be most general. The option leading to the desired result is a direct generalization of the spin 4 gauge transformations (3.90): To explicitly check that the transformation (4.17) is the needed one, we divide it into two terms. The first part of the variation, transforms the free hypermultiplet Lagrangian (up to a total derivative) as The second part of the transformation (4.17),

JHEP05(2022)104
gives rise to the following transformation of the free hypermultiplet Lagrangian: (4.21) So the full variation of the Lagrangian under (4.17) allows one to determine the possible interaction terms: As a consequence, the gauge invariant Lagrangian has the form: Like in all previous cases, it is the covariantization of the free hypermultiplet action under (4.17) through extending the harmonic derivative D ++ by the differential operator H ++ (s) (2.35). The highest derivative term in this operator have the degree (s − 1). Once again, the action corresponding to (4.23) is invariant only to the first-order in gauge superfields. For the time being, we do not know how to achieve the complete invariance of such an action.
It is worth noting that the group of rigid symmetries of the free hypermultiplet is much wider than (5.1). For any spin there are also transformations of the form:

JHEP05(2022)104
To be convinced of this, consider a variation of the free hypermultiplet action of the hypermultiplet with respect to (5.3): After integration by parts, one can rewrite this variation as: (5.6) Then the requirement of invariance of the action amounts to the condition (5.4).
For even spins s ≥ 2, this does not provide new possibilities as compared to (5.1), while for odd spins s ≥ 3 new rigid symmetries come out. These extra rigid symmetries generalize to the case of the hypermultiplet the symmetries suggested in ref. [37]. Perhaps, it would be of interest to explore their gauging. Leaving aside the detailed treatment of this interesting problem, note that such a gauging would imply the introduction of new gauge superfield prepotentials with a nontrivial index structure. This should also lead to a change in the structure of the action of higher spins and to the appearance of new indices for the component fields. So gauging of such symmetries could yield some "non-abelian higher spin theory". To the best of our knowledge, even for the bosonic higher spins, the interactions of this kind had never been seriously explored. The only relevant remark is contained in ref. [37], where it is claimed that the gauging of such symmetries "run into difficulties".
A similar (though different) kind of "non-abelian" higher spin invariances of the free hypermultiplet action (3.1), also existing for odd spins s = 1, 3, . . ., and generalizing those of ref. [37], is provided by the transformations δ (s) where the parameters are c-numbers and for s = 1 yield just rigid SU(2) PG transformations. So gauging of such symmetries should result in some higher-spin generalizations of Yang-Mills theory associated with the group SU(2) PG . So much for these new opportunities which we hope to study elsewhere.

2.
As the next crucial step after defining the global symmetries (5.1), we gauged them and explicitly presented the gauge transformations of 4D, N = 2 analytic hypermultiplet superfield. They constitute a family of transformations which differ by the highest degree s of the involved derivatives. These transformations are consistent with the analyticity and include the same differential operators as the gauge transformations of the higher-superspin analytic prepotentials (2.34). These transformations have the universal form:

JHEP05(2022)104
When the parameters in (5.8) are independent of the coordinates, these transformations reproduce the global symmetry group (5.1) of the free hypermultiplet action (3.1).

3.
As the final result, in section 3 we presented N = 2 gauge invariant cubic couplings 1 2 − 1 2 − s of the hypermultiplet to higher-spin s N = 2 gauge supermultiplets, which we described in section 2. These couplings can be brought in the universal form for all spins: Thus, the interactions for even and odd spins differ by the presence, in the odd spin case, of an additional interaction with the dimensionless coupling constant ξ. At the superfield off-shell level, this second interaction has the form clearly distinct from the basic interaction. The most important difference between these two interactions is that the ξ-terms contain no derivatives on the hypermultiplet superfield. The interplay of these two terms at the on-shell component level will be discussed elsewhere. Note that after passing to the gauge superfields of the canonical dimension, the constants κ s defined in (2.22) appear in (5.9) in front of the interaction terms, and so the interactions with hypermultiplets vanish in the limit κ s → 0 .

4.
The presence of the U(1) PG generator J of the internal SU(2) PG symmetry in the transformations (5.8) and couplings (5.9) results in an essential difference in the treatment of the odd and even spin s cases. For odd s the relevant higher-spin transformations of the hypermultiplet and its couplings to the gauge fields exist only provided SU(2) PG is broken, with ∂ 5 q +a = 0 , Jq +a = 0 in the case of massless hypermultiplet and ∂ 5 q +a = 0 , Jq +a = 0 , ∂ 5 q +a ∼ mJ in the case of massive hypermultiplet. The unbroken SU(2) PG symmetry implies ∂ 5 q +a = Jq +a = 0, so the gauge transformations of the hypermultiplet and the relevant gauge invariant couplings in this case can be defined only for the even spins s, the hypermultiplet being massless.
It is the appropriate place here to discuss how our consideration can be extended to the case of few hypermultiplets. The free Lagrangian of n hypermultiplets can be written in the manifestly USp(2n) invariant form as L +4 free,n = 1 2 q +A D ++ q + A , q + A = Ω AB q + B , A = 1, 2, . . . , 2n , (5.11)

JHEP05(2022)104
where Ω AB = −Ω BA is USp(2n) invariant constant 2n × 2n symplectic metric. Up to a numerical factor, this Lagrangian can be rewritten in an equivalent complex form as L +4 free,n ∼q +a D ++ q + a − D ++q+a q + a , a = 1, 2, . . . , n , q + A = (q + a , −q +a ) . (5.12) In this form it is manifestly invariant under the group U(n) = SU(n) × U(1) ⊂ USp(2n), with respect to which q + a andq +a transform in the fundamental and cofundamental representations, while the transformations completing U(n) to USp(2n) are realized as δq +a = C ab q + b (and c.c.), with C ab = C ba . In this general case it is natural to identify J with the common phase U(1) ⊂ U(n) , Jq + a = iq + a , Jq +a = −iq +a . (5.13) It breaks USp(2n) down to U(n) ⊂ USp(2n). 13 So in this case the residual group acting on the bosonic fields f i a is U(n) × SU(2) aut . All the higher-spin gauge transformations derived earlier for a single hypermultiplet and the relevant cubic couplings to the gauge superfields can be directly transferred to this general case. Note that for generic n > 1 there are much more possibilities to choose a pair of mutually commuting internal symmetry generators to be identified with J and, by the Scherk-Schwarz mechanism, with ∂ 5 .

5.
The results on the cubic N = 2 supersymmetric couplings of hypermultiplet to higher spins obtained in the present paper could be further extended along several directions (besides those already mentioned in the item 1): • The natural next step is the construction and investigation of 4D, N = 2 higherspin supercurrents of the hypermultiplet; • Of primary interest is also the study of the component structure of the interactions constructed. The bosonic physical fields f Ai of n hypermultiplets (4n independent fields) are transformed according to the bi-fundamental representation of USp(n) × SU(2) aut in the massless case or U(n) × SU(2) aut in the massive one, and it is tempting to examine how the bosonic subsector of the N = 2 higher spin gauge group acts on these fields; • An interesting task is to explore the relationship of the N = 2 couplings presented here with the known N = 1 cubic interactions. To accomplish this, it is necessary to reduce our N = 2 harmonic superspace to N = 1 superfields; • The approach under consideration can be applied to interactions of N = 2 higher spin theory with the hypermultiplet defined on AdS and more general superconformally-flat superbackgrounds. The hypermultiplet action in 5D AdS harmonic superspace was constructed in [81,82] and can be directly dimensionally reduced to 4D. It is still unknown how to extend the harmonic superspace JHEP05(2022)104 construction of the higher spin theories given in [8] to 4D AdS framework. The linearized action for N = 2 supergravity on the AdS background was built in [83]. However, this formulation was never generalized to higher spins.