Dynamical renormalization group for mode-coupling field theories with solenoidal constraint

The recent inflow of empirical data about the collective behaviour of strongly correlated biological systems has brought field theory and the renormalization group into the biophysical arena. Experiments on bird flocks and insect swarms show that social forces act on the particles' velocity through the generator of its rotations, namely the spin, indicating that mode-coupling field theories are necessary to reproduce the correct dynamical behaviour. Unfortunately, a theory for three coupled fields - density, velocity and spin - has a prohibitive degree of intricacy. A simplifying path consists in getting rid of density fluctuations by studying incompressible systems. This requires imposing a solenoidal constraint on the primary field, an unsolved problem even for equilibrium mode-coupling theories. Here, we perform an equilibrium dynamic renormalization group analysis of a mode-coupling field theory subject to a solenoidal constraint; using the classification of Halperin and Hohenberg, we can dub this case as a solenoidal Model G. We demonstrate that the constraint produces a new vertex that mixes static and dynamical coupling constants, and that this vertex is essential to grant the closure of the renormalization group structure and the consistency of dynamics with statics. Interestingly, although the solenoidal constraint leads to a modification of the static universality class, we find that it does not change the dynamical universality class, a result that seems to represent an exception to the general rule that dynamical universality classes are narrower than static ones. Our results constitute a solid stepping stone in the admittedly large chasm towards developing an off-equilibrium mode-coupling theory of biological groups.

The recent inflow of empirical data about the collective behaviour of strongly correlated biological systems has brought field theory and the renormalization group into the biophysical arena. Experiments on bird flocks and insect swarms show that social forces act on the particles' velocity through the generator of its rotations, namely the spin, indicating that mode-coupling field theories are necessary to reproduce the correct dynamical behaviour. Unfortunately, a theory for three coupled fields -density, velocity and spin -has a prohibitive degree of intricacy. A simplifying path consists in getting rid of density fluctuations by studying incompressible systems. This requires imposing a solenoidal constraint on the primary field, an unsolved problem even for equilibrium mode-coupling theories. Here, we perform an equilibrium dynamic renormalization group analysis of a mode-coupling field theory subject to a solenoidal constraint; using the classification of Halperin and Hohenberg, we can dub this case as a solenoidal Model G. We demonstrate that the constraint produces a new vertex that mixes static and dynamical coupling constants, and that this vertex is essential to grant the closure of the renormalization group structure and the consistency of dynamics with statics. Interestingly, although the solenoidal constraint leads to a modification of the static universality class, we find that it does not change the dynamical universality class, a result that seems to represent an exception to the general rule that dynamical universality classes are narrower than static ones. Our results constitute a solid stepping stone in the admittedly large chasm towards developing an off-equilibrium mode-coupling theory of biological groups.

I. INTRODUCTION
The phenomenon whereby systems with very different microscopic details have similar critical behaviour, known as universality, finds an elegant explanation in the context of the Renormalization Group (RG) [1,2]. As a system approaches a critical point, i.e. a second order phase transition, it exhibits emergent collective phenomena on increasingly wider length-scales, which leads to the presence of a diverging correlation length. In this regime, a set of phenomenological scaling laws for near-critical systems have been proposed [3,4] in order to describe the critical behaviour through a set of critical exponents ruling the divergence of thermodynamic quantities with respect to the correlation length. In the RG framework, these scaling laws naturally arise as a consequence of the diverging correlation length, thanks to the flow in the parameter space generated by iterating an RG transformation many times. Through the concept of attractive fixed points of the RG flow, it is possible to prove that systems sharing only few general properties, as symmetries, dimensionality and range of the interactions, develop the same large-scale behaviour near the critical point, hence belonging to the same universality class. Moreover, the RG gives a constructive method to compute the critical exponents [5], making it one of the most powerful tools of the theory of critical phenomena, broadly employed to characterise both static [6] and dynamic [7] universality classes in equilibrium systems.
When we shift our attention from equilibrium to nonequilibrium phenomena, the study of emergent collective dynamics is of fundamental interest also in biological systems, in which the response to external perturbations is more efficient when individuals behave collectively [8,9]. The existence of collective behaviour in biological systems is strictly related to the presence of scale-free correlations among the system's components, as experimentally observed in bird flocks [10], sheep herds [11] and bacterial clusters [12], in strong similarity with emergent phenomena studied in condensed matter, as superfluidity, superconductivity and ferromagnetism. The experimental discovery of scale-free correlations provides some justification to the theoretical investigation of collective behaviour in biological systems using the tools developed for critical phenomena, including of course the renormalization group.
The early efforts in the study of collective phenomena in active biological systems focused on the ordered phase, namely in those groups exhibiting net group motion, the most prominent example being represented by bird flocks [13,14]. More recently it has become experimentally evident that also systems in their disordered phase, as natural swarms of insects, display collective phenomena even in absence of collective order [15]. Moreover, within natural swarms both static and dynamic scaling laws have been found to hold [15,16]. Since scaling laws are funda-mental features of near-critical systems, a renormalization group approach to swarming behaviour seems particularly appealing.
Dynamic critical systems not only have a diverging correlation length, but also a diverging relaxation time [17]. According to the dynamic scaling hypothesis, the correlation length ξ and the collective (i.e. zero wavelength) relaxation time τ are linked by the relation τ ∼ ξ z , where z is the dynamic critical exponent ruling the dynamical relaxation of the system. Experiments show that the dynamic behaviour of natural swarms is characterized by a critical exponent z ∼ 1 [16], which is very far from the value z ∼ 2 typical of the simplest dissipative ferromagnetic models, generically represented by Model A in the famous classification of Halperin and Hohenberg [7]; this is somewhat problematic, as the earliest theories for active biological systems are essentially a mix of Model A ferromagnets and Navier-Stokes equations [13]. Finding what are the correct dynamical equations that provide a value of the dynamical critical exponent z equal to the experimental value represents a very important step, not only towards understanding collective behaviour in natural swarms (which will be our focus here), but more broadly in the effort to prove that statistical field theory is a predictive tool in biophysics.
There are two possible mechanisms that can explain a drop of the value of the dynamical critical exponents in natural swarms compared to that of standard ferromagnetic systems in the Model A class. The first possibility is activity, namely the fact that individuals are self-propelled, which leads to a time-dependent network of interactions between the particles' velocities, thus driving the system out of equilibrium. Activity gives rise to density fluctuations, thus meaning that in addition to the velocity field v (x, t), also the density field ρ (x, t) must be included in the hydrodynamic description of the system [13,14]. Broadly speaking, one can say that the feedback between velocity and density fluctuations is responsible for a more efficient transport of information within the system, which may reflect into a lower value of z. Indeed, theoretical RG calculations in presence of activity have been performed in the incompressible case, showing that z is indeed lower than its equilibrium value, leading the exponent from z ∼ 2 to z ∼ 1.7 in d = 3 dimensions [18]; and yet this is still far from the experimental value z ∼ 1.
The second factor that may lower the value of z, and that has been shown to be relevant in the description of both insect swarms and bird flocks, is the presence of inertial dynamics, namely a coupling between the velocity field v (x, t) and the generator of its rotations, known as the spin s (x, t) [16,19]. In equilibrium systems, as for example planar ferromagnets, superfluids and quantum antiferromagnets (respectively Model E, F and G in the classification of [7]) this coupling (and the conservation law generated by it) radically alters the mechanism of information propagation throughout the system and it also lowers the dynamical critical exponent, which becomes z = 1.5 in d = 3 [20,21]. Even in this case, though, we are quite far from the experimental value, but then it seems likely that a calculation that keeps into account both the off-equilibrium active nature of swarms and the mode-coupling nature of their interaction, may provide a value of the dynamical critical exponent close to the experimental one, z ∼ 1. The bad news is that a fullfledged out of equilibrium field theory taking into account both these features must contain a density field ρ (x, t), a velocity field v (x, t) and a spin field s (x, t), with a large number of couplings between them and a near-prohibitive level of intricacy. Therefore, to make theoretical progress, one needs some simplifying assumptions.
The first and most natural simplification is that of incompressibility. Incompressibility in the RG study of the Navier-Stokes equations was first used in the seminal paper by Forster, Nelson and Stephen [22], while it was used first in the context of dynamical field theories for active matter in [18]. This hypothesis seems quite reasonable both in natural insect swarms and bird flocks, because experiments show that density fluctuations are indeed very limited in these systems, and they do not seem to be a crucial ingredient in their collective behaviour [15,23]. Incompressibility is obtained by enforcing a solenoidal constraint on the velocity field, that is by adding to the dynamics the condition, ∇ · v = 0. This allows to neglect the (constant) density field, whose fluctuations are suppressed thanks to the solenoidal constraint. Incompressible swarms in the fully dissipative case (i.e. without velocity-spin coupling) have been studied in [18], but we need to move to the mode-coupling case, of velocity and spin, and here some new problems arise. In the case of a theory with one single field (the velocity), the solenoidal constraint is simply enforced by projecting the force acting on the velocity onto the direction orthogonal to the wave vector k in Fourier space [18,22]; in this way, the constraint is conserved along the time trajectories. But in systems with mode-coupling dynamics incompressibility cannot be achieved in this simple way, as the forces now act on the spin, rather than on the velocity, and the spin is not subject to the solenoidal constraint! Hence we are in a paradoxical situation: the assumption that was supposed to radically simplify the calculation, reducing from three to two the number of fields, creates itself a new challenges.
Instead of giving up the crucial simplification of incompressibility, we choose to try and solve the problem of how to impose a solenoidal constraint in a mode-coupling theory; and because this new problem has nothing to do with the complications due to off-equilibrium dynamics, we tackle the problem in the case of a zero-activity equilibrium system, namely a system living on a fixed-network. Indeed, surprising as this may seem, the problem of how to impose the solenoidal constraint on a mode-coupling theory has never been investigated so far. Once we will have solved this problem, we will be in a far better shape to push the calculation out of equilibrium, still retaining the incompressibility condition. Notice that in absence of activity it does not make much sense to insist in calling the primary field "velocity". Hence, we will use as our primary field ψ, the field describing the average direction of motion in the vanishing speed limit. The virtue of this approach is that any future calculation performed in the active case must have the present calculation as a limit in the zero-speed case.
To summarize, in the present work we will derive and study the coarse-grained dynamic equations for a critical field ψ which is coupled to the generator of its rotations, s, in the presence of the solenoidal constrain, ∇ · ψ = 0. We shall refer to this theory as Solenoidal Model G (SMG), because in absence of the solenoidal constraint this field theory is known in literature as Model G [7]: this classic model for equilibrium quantum antiferromagnets will therefore be our stepping stone towards a theory for off-equilibrium natural swarms. We will show that the solenoidal constraint leads to the emergence of an additional non-linear term in the equation of motion of the spin, as a consequence of the suppression of the ψ mode longitudinal to the wave-vector k. We will also prove that the equations we obtain are eigenstates of the renormalization group, by calculating perturbative corrections up to one loop. Moreover, the dynamic behaviour will be shown to be consistent with the static behaviour of a ferromagnet with dipolar interactions [24,25]. Finally, we will find as dynamic critical exponent z = d 2 , leading to the same dynamic universality class as the unconstrained theory [7,20]. Therefore, while incompressibility affects the static universality class, it does not modify the dynamic universality class, suggesting that the critical behaviour of homogeneous systems is not affected by imposing incompressibility.
In section II we will review the biological origin of the mode-coupling dynamics and the reason for which incompressible constraint is imposed. In section III the coarse-grained theory will be derived, starting form the static Hamiltonian and the Poisson-bracket relations between the various quantities. In section IV we will derive one-loop corrections to the couplings of the theory and verify that the equations of motion are self-consistent. Self-consistency is proven by showing that no new RGrelevant interaction is generated and that the dynamic renormalization reproduces the correct static behaviour. Finally, in section V we will derive the dynamical critical exponent.

A. Dynamic scaling
The spatio-temporal statistical behavior of a collective system can be described through the connected spacetime correlation function of the velocity, C(r, t), which expresses how much the velocity fluctuations of some individual in the system influence and are influenced by the velocity fluctuations of another individual at distance r and after a time t [26]. For t = 0 this is just the static (i.e. equal times) correlation function, whose decay range defines the correlation length, ξ. Conversely, for r = 0 and t = 0, this is the singe-particle auto-correlation function. In general, when both r and t are non-zero, the correlation function is quite a loaded concept, entailing the fundamental relationship between relaxation in time and relaxation in space. In the case of critical systems, such relationship acquires a particularly illuminating form, which goes under the name of dynamic scaling [27,28]. Dynamic scaling states that the dynamic correlation function C, when expressed as a function of wave-vector and frequency, takes the following simplified scaling form, where ξ is the correlation length and where the static correlation function C 0 has in turns the scaling form, and where the characteristic frequency at scale k is given by, In the relations above, Ω, F 0 and F are well-behaved scaling functions, whose explicit form is inessential to capture the gists of dynamic scaling [27,28]; η is the critical exponent for the static correlation function (normally called anomalous dimension [29]). The fundamental meaning of dynamic scaling is that in critical systems space and time do not scale independently from each other, but they are linked by the dynamic critical exponent z. The space-time correlation function has a very simple form, as its whole spatio-temporal dependence goes through the product kξ. When k = 0, the collective relaxation time τ of the system is linked to the correlation length ξ through the relation, This phenomenon is known as critical slowing down and it represents a consequence of the fact that the time τ needed to decorrelate a spatially correlated region grows with the region's size, namely with the correlation length ξ, making the latter the only relevant scale also at a dynamic level.
Notably, natural swarms of insects have been found to obey dynamic scaling [16], with a dynamical critical exponent z quite close to 1 (z = 1.2, obtained by looking at the collapse of the correlation functions [16], is the safest determination). On one hand, because the validity of scaling laws is one of the hallmarks of critical systems, the experimental evidence of static and dynamic scaling in natural swarms [15,16] suggests that the swarming phase can be theoretically described as a near-ordering phase of classic ferromagnetic theories. On the other hand, the value z ∼ 1 is definitely anomalous, as the standard exponents in Model A class of statistical systems [7] is z ∼ 2, with small corrections at two loops. This state of affairs suggest that a critical dynamical universality class, different from Model A, must be found for natural swarms. This new universality class will need to take into account the two main features that affect the collective behaviour of these systems, namely activity and mode-coupling dynamics. Let us illustrate the first factor through the most venerated hydrodynamic theory describing active matter.

B. The hydrodynamic theory of Toner and Tu
In most collective biological systems, from bacterial clusters up to insect swarms, bird flocks and animals herds, individuals are self-propelled thanks to their ability to convert energy into systematic movement [30]. The presence of metabolic processes, through which an energy supply from the environment is provided, guarantees a constant speed. Activity means that the interaction network over which the velocities of the individuals are interacting with each other evolves in time, representing the principal non-equilibrium feature of these systems.
The hydrodynamic theory developed by Toner and Tu [13,14], based on the discrete model proposed by Vicsek and collaborators [31], paved the way to a theoretical analysis of collective behaviours in the presence of activity. Together with the velocity field v, also the density number field ρ is an hydrodynamic variable of the system, since activity allows local density fluctuations. The velocity field obeys to a dynamical behaviour that is a crossover between that of a Landau-Ginzburg O(n) model, reflecting the presence of an effective alignment among the animals velocities, and a Navier-Stokes dynamics, reflecting the fact that velocity and density are coupled as in a standard fluid; the density evolves according to a continuity equation, since the total number of individuals is assumed to be fixed. The essential terms of the Toner and Tu (TT) theory are given by [32], In Eq. (5), the r.h.s. encloses Model A dynamics [7], with the gaussian random white noise θ and the hamiltonian force − δH δv , to which a pressure force is added contrasting an infinite compressibility. Here the effective hamiltonian H takes the Landau-Ginzburg structure, namely Depending on the value of r 0 , two phases can be identified: an ordered phase for r 0 < r c and a disordered phase for r 0 > r c , where r c is the critical value of the parameter r. The main difference with equilibrium dynamics is represented by the γ v term on the l.h.s. of Eq. (5), which is the advection term 1 typical of the Navier-Stokes structure. Since its development, the TT theory has been able to explain many fundamental aspects of collective behaviours in active systems, from the existence of an ordered phase also in 2 dimensions, to the presence of linear sound modes [13,33]. The success of this remarkable theoretical effort is due to the minimal description given by the TT equations of motion, which tie together activity and effective alignment. However, experimental evidence has shown how this theory is not able to fully describe some aspects of biological systems', in particular information propagation in flocks [19], and dynamical correlation functions in swarms [16]. These discrepancies arise as consequence of the structure of Eq. (5), which is essentially that of an overdamped Langevin equation for the velocity, Moreover, RG analysis of the near-ordering dynamics of the TT theory in the incompressible case predictss a dynamical critical exponent z 1.7 for d = 3 [18], significantly different from the value z ∼ 1 observed in natural swarms.

C. Restoring inertia: the mode-coupling theory
The first hints about the necessity of restoring an underdamped dynamics came from experiments on bird flocks. In order to explain the information propagation dispersion law within turning flocks, the generator of the rotations of the velocity vector has to be conserved during the turn [19]. In analogy with quantum mechanics, this generator of the rotations in the internal space of the velocity, which must not be confused with the generator of spatial rotations (the standard angular momentum), takes the name of spin.
A second piece of empirical evidence that inertial dynamics was missing from the original Vicsek model and Toner-Tu theory came from the calculation of the dynamical correlation functions in real swarms [16]. The overdamped dynamics for the velocity typical of the Vicsek model and TT theory generates a classic exponential decay of the dynamical correlation functions; however, in natural swarms this is not the case, and the correlation function displays a non-exponential form typical of underdamped inertial systems. This feature too can be included at the theoretical level by incorporating in the dynamics the generator of the rotations of the velocity.
The discrete model describing this kind of spin-velocity mode-coupling behaviour is known as the Inertial Spin Model (ISM) [34], and it is characterized by an underdamped structure for the equations of motion; where the alignment force exerted on the focal particle by its neighbours does not act directly on the velocity, but it is mediated by the generator of rotations, thus restoring an inertial behaviour.
Adding a conservation law leads to a modification of the hydrodynamic behaviour of a system, which now must include the dynamics of the spin density and the effects of a conserved total spin on the other variables. In the present case, the conservation of the spin leads a mode-coupling between the order parameter and the spin density field. This structure is not new to physical systems: this is what happens in the case of planar magnets, superfluid helium and isotropic anti-ferromagnets (Models E, F and G in [7] respectively). The equations of motion for these models, in terms of the order parameter ψ and the spin s, take the following form for a 3 dimensional order parameter (Model G) where θ and ζ are gaussian random white noises, while the Hamiltonian has the usual Landau-Ginzburg form for ψ, plus a Gaussian non-interacting kinetic term in the spin, (9) and (10) the two terms proportional to g, namely ∂ t ψ ∼ δ s H and ∂ t s ∼ δ ψ H, represent the Hamiltonian conservative dynamics enforced by the Poissonbraket relation where αβγ is the Levi-Civita antisymmetric symbol, meaning that s is rotating ψ. Moreover, the equation of motion of s can be written in the form of a continuity equation, reflecting the fact that due to the rotational symmetry, the total spin S(t) = d d x s(x, t) is conserved 2 .

D. Incompressible flow and fixed network assumptions
In order to join the two ingredients described in the last two sections, namely activity and inertia (i.e. modecoupling dynamics) into one single dynamical field theory, 2 In the biophysical context, the spin is not strictly conserved.
Thus, the presence of a dissipative term ∂ts ∼ −ηδsH is allowed in Eq. (10). In this system a crossover between a conservative and a dissipative dynamics has been observed [20,21], in which the former regulates the behaviour of finite-size systems as swarms.
As we will see later on, the solenoidal constraint will violate the spin conservation, but does not generate a dissipative behaviour, thus allowing us to work at η = 0 one needs to write equations for three coupled fieldsdensity, velocity and spin -with several non-linear couplings among them, giving rise to a diagrammatic RG proliferation that it is impossible to keep under control without any simplifying assumption. One such simplification is that of incompressibility. The huge advantage of working under the hypothesis of incompressibility is that the density field drops out the theoretical description [18,22], thus reducing the number of fields that have to be studied from three to two, hence dramatically simplifying the theoretical investigation. One may worry that incompressibility introduces some non-local interactions in the system, then potentially changing the critical exponents with respect to the compressible case. However, numerical simulations of the standard compressible theory have found the dynamic scaling exponent to be in perfect agreement with theoretical RG results for incompressible theory [35], thus indicating that the effects of activity on the dynamical critical exponent are independent from whether incompressibility is enforced or not. This is very useful indeed, as it means that incompressibility can be used as a simplifying tool, without changing critical dynamics.
A second remark about incompressibility is in order. In the standard compressible theory, the presence of density fluctuations in the system, as a consequence of activity, makes the phase transition of Vicsek-like models a firstorder (i.e. discontinuous) transition [36]. In incompressible systems, however, the absence of density fluctuations makes the phase transitions always second-order (continuous) [18], thus scaling behaviour is observed at all sizes and all the complications about phase separation are avoided. Of course, this is only reasonable because in real natural swarms of insects (as in bird flocks), density fluctuations are negligible and play no determinant role in ruling the collective dynamics of the system [15,37]. Hence, incompressibility is not only a simplifying hypothesis, but also a sound theoretical description of actual empirical data.
Thanks to the continuity equation (6), the condition of incompressibility reduces to imposing a solenoidal constraint on the velocity field, As we shall see, the introduction of a solenoidal constraint on the primary field within a mode-coupling dynamics is far from trivial. The only encouraging thing is that the complications arising are completely unrelated to activity, so that it seems sound to first find a conceptually consistent way to impose the solenoidal constraint on a theory at equilibrium, and then to use that results to make progress off-equilibrium in the future. Therefore, in this work activity will be neglected by assuming the adjacency network to be fixed in time. Since in the microscopic description of Vicsek-like systems each individual has a fixed speed v 0 , the fixed-network approximation can be formally seen as a limit where v 0 vanishes, which is equivalent to freeze the position of each particle. While in this limit the velocity field acquires a singular behaviour, since no particle is actually moving, the coarse-grained direction of motion ψ, defined by the relation still has a smooth behaviour when v 0 → 0, and therefore represents the ideal candidate to be the order parameter in the fixed network approximation. Since the solenoidal constraint given by Eq. (13) holds also for ψ at any finite value of v 0 , it must hold also in the fixed network approximation.

III. THE SOLENOIDAL MODE-COUPLING THEORY
The first problem we have is that the presence of the solenoidal constraint, forces the order parameter to have the same dimensionality as the space in which the theory is defined. Thus, since in the following sections we will need to perform an RG expansion near d = 4, we must work with the generalization of Model G to arbitrary dimensions, known as the Sasvari-Schwabl-Szepfalusy (SSS) model [38,39]. 3 In this model, the spin s is an anti-symmetric d × d matrix, with one independent component for each possible plane around which a rotation can be performed. The Poisson-bracket relation between s and ψ becomes, where repeated greek-letter indices are intended to be summed, if not otherwise indicated, and the tensor I is the identity tensor in the space of d × d anti-symmetric matrices, 4 given by The static properties of the SSS model are fully determined by its effective Hamiltonian, which reflects the fact that ψ is a critical field with O (d) symmetry group and local interactions, having the static critical properties of the Landau-Ginzburg universality class, while s is a non-critical massive field. The effective Hamiltonian of the system takes the following natural generalization The equations of motion of the critical SSS model can be derived starting from the Poisson-bracket relation (16) and the effective hamiltonian (18) by following a procedure coming from the works of Mori et al. [40] and Zwanzig [41] (see [42] and [43]). These equations are given by where θ and ζ are two gaussian white noises, with variance The λ 0 ∇ 2 terms in the stochastic parts of the dynamics of s ensure that the spin is conserved. 3 Given that we will be calculating the properties of the Solenoidal Sasvari-Schwabl-Szepfalusy model, strictly speaking we should use the nomenclature "SSSS model"; however, chiefly for aesthetic reasons, we prefer not to do that, and we will rather stick to Solenoidal Model G (SMG) even in generic dimension d. 4 The factor 1 2 in the definition of I arises as a consequence of the fact that, when s is represented as an anti-symmetric matrix, each independent component appears twice.

A. Static critical behaviour of the solenoidal theory
We will now review, for the benefit of the reader, the classic results of the universality class of dipolar ferromagnets [24,25], since the static behaviour of the ψ field in SMG belongs to this class. From now on, if not explicitly indicated, we shall work with the Fourier transformed fields, defined by the relation where we introduced the notation In Eq. (23) the cutoff Λ is the value of the wave-vector above which fluctuations have no physical meaning. Since a field theory, in the context of biophysics, is obtained by coarse-graining a discrete model, fluctuations cannot occure on distances shorter than the microscopic length a, therefore meaning that the cutoff can be taken to be Λ ∼ a −1 . In Fourier space the solenoidal constraint reads An equivalent way in which this can be formulated is by requiring that where we have introduced the transverse projection operator The static behaviour of a field characterized by the Hamiltonian of the SSS model, given in Eq. (18), to which a solenoidal constraint is applied, is described by the static universality class of isotropic dipolar ferromagnets [24,44]. The universality class of dipolar ferromagnets describes the more generic class of critical theories described by an Hamiltonian with a structure as that of Eq. (18), in which the ψ (k) mode is taken to be non-critical. The renormalization group analysis of isotropic dipolar ferromagnets [24] shows that the non-criticality of the ψ (k) mode leads to its full suppression in the longwavelength limit, meaning that the RG stable fixed point describes a solenoidal-constrained theory. To order ε = 4 − d, the RG recursive relation for the ferromagnetic coupling constant is [24], In dimension d < 4 the stable fixed point u * ruling the critical behaviour of dipolar ferromagnets is given by [24] u * = 2 17 ε (for the sake of simplicity we have set to 1 the volume of the d-dimensional unit sphere). The fact that at the stable fixed point the ψ (k) mode is suppressed [24,44] means that the solenoidal theory is robust with respect to weak violations of the constraint, i.e. small non-critical fluctuations of the ψ (k) mode. Even though the change of static universality class is very interesting at the theoretical level, it must be said that the new critical exponents are so close to those of the Landau-Ginzburg universality class that the difference is often experimentally difficult to observe [44]. A two-loop estimation of the critical exponents 2ν, ruling the behaviour of the mass r while the transition is approached, and η, modifying the spatial dependence of the correlation function, are reported in Table I. Here it is clear the little difference between the scaling behaviour of dipolar ferromagnets class and Landau-Ginzburg class.

B. Dynamics of the solenoidal theory
In this crucial section, the dynamical equations of the mode-coupling theory subject to a solenoidal constraints will be derived. We will use the Poisson relations, the effective Hamiltonian of the SSS model, and the classic Mori-Zwanzig formalism [40,41].
Because the field ψ is subject to the constraint, the most natural representation in which one would like to derive the equations of motion is, whereê = k/|k| is the unitary vector identifying the direction of k, while theê ⊥,i are orthogonal unitary vectors spanning the space perpendicular to k, in such a way that e ·ê ⊥,i = 0 andê ⊥,i ·ê ⊥,j = δ ij ; indeed, within this set of coordinates, the constraint simply reads, leaving only the (d − 1) independent modes ψ ⊥,i (k) to take care of. The advantage of this notation is that we can formulate the constrained field theory in terms of the (d − 1) independent transverse modes and not in terms of a constrained set of d cartesian coordinates, ψ α , to which it is not clear how to apply the Mori-Zwanzig procedure. However, the explicit form of both the effective hamiltonian H and of the Poisson-brackets are given in terms of the cartesian coordinates, ψ α , while their form in terms of the ψ ⊥,i would be extremely cumbersome. What we shall do, then, will be to first obtain the equation of motions for the ψ ⊥,i and then to go back to the standard ψ α fields by using the chain rule. In doing that, something new will pop out in the spin equation.
The dynamic behaviour of the constrained variables ψ α is given, in terms of the independent variables ψ ⊥ , by Following the Mori-Zwanzig procedure, the equations of motion for the (d − 1) independent fields ψ ⊥,i and s take the following form where the t dependence is always understood even when not made explicit for reasons of space. In Eqs. (33) and (34) H is the Hamiltonian given by Eq. (18) in which ψ is set to 0;θ, ζ are white gaussian noises with variance respectively given by 2Γ 0 δ αβ and 4k 2 Λ αβγν . The tensor Λ αβγν = Λ αβγν (λ 0 ; k) is function of the diffusive coefficient λ 0 and potentially also of the wave-vector k; in the field theory without constraint it takes the form Λ αβγν = λ 0 I αβγν ; however, since the solenoidal constraint generates an anisotropy in Fourier space for the order parameter, we expect that anisotropic effects can affect also the spin dynamics; therefore, we generalize Λ αβγν by allowing λ 0 to take different values for the longitudinal and transverse components, namely by taking, where P is the generalization of the projection operator P defined in Eq. (26) acting on the space of 2-indices antisymmetric tensors, and takes the following form In order to find the equations of motion of ψ α and s αβ , we will proceed by making explicit the terms in Eqs. (33) and (34) exploiting the chain rule between the fields ψ ⊥ and ψ α . The chain rule applied to the variations of the Hamiltonian with respect to ψ ⊥,i reads, while the Poisson-bracket relation between s and ψ ⊥,i can be written as where in the last equality of both Eq. (37) and (38) we used the relations, Thanks to Eq. (37) and (38) we can write the mode-coupling term of the dynamics of the spin in the following way where we used the following relation, which comes directly from the definition of the unitary vectorsê ⊥,i Thanks to the relations we found in Eq. (37), (38) and (40) we can write Eq. (33) and (34) in the following form By substituting Eq. (42) in Eq. (32) and by using the explicit Poisson-bracket relation given in Eq. (16), the equations of motion for the standard Cartesian components ψ α and for s αβ can finally be written, Here the hamiltonian H can be taken to be the same of SSS model, since all the terms involving derivatives of H with respect to ψ are already projected. Moreover, the two gaussian random white noises ζ and θ have variance

Effects of the constraint
The solenoidal constraint has a double effect on the equations of motion. The first, and maybe the most trivial, is that the equation of motion for ψ is projected orthogonally to k, as it happens in incompressible field theories with no mode-coupling interaction [18,22]. The second effect is less obvious, and it is represented by the presence of a projection operator P νρ (q) in the mode coupling interaction of the spin dynamics. The existence of this projector is a consequence of the fact that, in the presence of a solenoidal constraint, the conservative Hamiltonian force is not simply, but rather, The linear part of the force is not affected by this new projector, but the non-linear terms are. This can be seen by writing explicitly the new force, where in the first linear term we used the fact that P νρ (q) ψ ρ (q) = ψ ν (q). The linear part of the force contributes to the dynamics of s with the same term as Model G [7] The factor q 2 − (k − q) 2 arises as a consequence of mode-coupling and it vanishes as k → 0, thus conserving the total spin S(t) = d d x s (x, t) = s (k = 0, t). But now, thanks to the presence of the projector, also the non-linear term contributes to the dynamics of the spin! More precisely, it does so through a novel dynamical interaction term given by, This is a completely new term, which mixes the static ferromagnetic interaction (the coupling constant u 0 ) with the dynamic mode-coupling interaction (the coupling constant g 0 ); such vertex is absent in the non-constrained theory, since when P νρ is substituted by δ νρ , as in the non-solenoidal case, the r.h.s. of Eq. (52) vanishes. In honour of the historic NYC diner at Houston and Ludlow, we call this new interaction, the Katz vertex. As we will demonstrate later on, the Katz vertex is crucial in order to keep closed and self-consistent the RG calculation and to recover the correct static critical exponents. The Katz interaction does not vanish when k → 0, meaning that the equation of motion of s cannot be written as a continuity equation anymore and thus that the total spin S(t) is no longer conserved. However, in the following section we will show thats no spin dissipation is generated by the Katz vertex, suggesting that the violation of the spin conservation is equivalent to a generalized precession of the total spin vector.

C. Field theoretical description
In order to set up the renormalization of the theory we will follow the procedure proposed by Martin, Siggia, Rose [45], Janssen [46] and De Dominicis [47] to write stochastic differential equations as a field theory formulated using path integrals. Thanks to this procedure, the behaviour of a field φ governed by a stochastic differential equation with a deterministic evolution operator F and a gaussian noise θ can be described through a field-theoretical action that correctly reproduces the statistics, i.e. the correlation and response functions, of Eq. (53), where 2L αβ is the variance of the gaussian noise, whilê φ is an auxiliary field. From now on, if not explicitly mentioned, we will work in Fourier space for both the spacial and temporal dependence, s wherek = (k, ω) and k = k . The presence of the additional fieldφ is the cost which has to be paid in order to exploit standard path integral formulation, which will allow us to use the standard rules of static renormalization and write the perturbative series in terms of Feynman diagrams. This auxiliary field takes the name of response field, since the propagator, namely the response function, can be written as [48]: The Martin-Siggia-Rose (MSR) action for the stochastic differential equations (44) and (45) while S I , which takes contributes from the non-linear dynamic terms and represents the interaction part of the action, is given by where we have introduced the new parameters, which are a sort of dynamical generalization of mass and ferromagnetic coupling constant. In Eq. (60) the following new tensors have been introduced Q αβγν (k) = P αβ (k) δ γν + P αγ (k) δ βν + P αν (k) δ βγ (62)

Free theory
The starting point to build the perturbative expansion of the equations of motion is the free, or Gaussian, dynamic theory, obtained by setting to 0 all the dynamic non-linear couplings, namely g 0 and u 0 . From the gaussian part of the action, given by Eqs. (58) and (59), we can easily derive the expressions for the bare propagators and correlation functions for the effective field theory, which are given by: The subscripted 0 on thermal averages indicate that they are computed within the non-interacting theory, namely by setting u 0 = g 0 = 0. The tensors G and C are given by In Eq. (66), (68), (67) and (69) we have, In the diagrammatic framework, bare propagators and correlation functions are represented in the following way where the arrows in the propagators aways point in the direction of the response field.

Non-linear terms: the vertices
The four terms that compose S I represent the nonlinear interactions in the equations of motion. Two of them involve one fieldψ, since they derive from the equation of motion of ψ, while the other two involve one fieldŝ, since they derive from the equation of motion of s. In the diagrammatic framework these interactions are represented by vertices, in which different lines merge together, each representing one of the fields involved in the interaction. Here we are representing with a solid line the fields ψ andψ, with wavy lines the fields s and s. Moreover, an arrow is used to recognize which legs represent a response field.
The first vertex -namely interaction -involvingψ, represents the mode coupling non-linearity, (75) This interaction represents a purely dynamic interaction, since it is proportional only to the dynamic coupling g 0 .
The second vertex involvingψ derives from the ferromagnetic quartic interaction of the static hamiltonian, ensuring that the field ψ relaxes towards the static equilib-rium distribution, and therefore is proportional to static coupling u 0 . It is represented by the term (76) The other two vertices involve one fieldŝ, and both derive from the mode-coupling interaction in the equation for s. The first comes from the linear part of the "force" defined in Eq. (49), representing a purely dynamic interaction proportional to g 0 , and it takes the usual form as in the non-constrained theorŷ Here the factor (p 2 − q 2 ), coming from the cross product structure of the mode-coupling interaction, guarantees that this interaction vanishes when k = 0. This is a consequence of the symmetry of the non-constrained theory, which conserves the total instantaneous spin S(t) = s(k = 0, t).
The last interaction term is the Katz vertex, given by the novel non-linear contribution peculiar of the solenoidal theory, discussed at the end of Sec. III B. This interaction mixes static and dynamic terms, since it represents the effects of the static quartic interaction on the dynamics of s, mediated by the mode-coupling dynamic interaction. Therefore, the Katz vertex is proportional to the product of the static coupling u 0 and the dynamic coupling g 0 and it takes the following form, The Katz vertex arises as a consequence of the solenoidal constraint, since the tensor K αβγνστ vanishes in the nonconstrained theory. At variance with the vertex (77), Katz causes a violation of the spin conservation due to the fact that the order parameter has lost the O(d) symmetry as a consequence of the solenoidal constraint, meaning that it does not vanish when k = 0. However, the fact that the total spin is not conserved does not mean that it is dissipated. In fact, we will show in the following sections that no perturbative corrections dissipating the spin arise after the shell integration as a consequence of the presence of the Katz vertex, reinforcing the hypothesis according to which s is an hydrodynamic slow-variable of the system. Note that, were this not the case, the RG flow would lead to the stable fixed point of solenoidal Model A [20].

IV. RENORMALIZATION GROUP CALCULATION
The key idea behind the renormalization group is that, under the assumption of an infinitely large correlation length, scaling laws and critical exponents can be obtained by looking at how the parameterss of a theory change by changing the length-scale at which the system is observed [1]. The RG itself consists in a transformation through which a set of equations describing the dependence of the couplings from the length-scale, namely the β-functions, can be derived.
We will use Wilson's momentum shell approach [2], which can be performed by following a two-steps procedure: i) integrating out the short wavelength details, hence decreasing the cutoff in momentum space; ii) rescal-ing space and time, so to formally reinstate the same original cutoff. The first step is carried out by marginalizing the probability distribution of the fields over the modes in the momentum shell Λ/b < k < Λ, where Λ is the cutoff. Here b is a parameter larger than 1, but close to it, while Λ is the cutoff of the theory. The effect of this integration is to modify the parameters of the theory, and shift the cutoff from Λ to Λ/b. To compare the coefficients of the new theory with the bare ones, space, and consequently time, must be rescaled in order to restore the original cutoff. The iteration of this procedure defines a flow in the parameter space, from which information on the critical behaviour can be obtained.
In the Gaussian theory, namely when u 0 = g 0 = 0, the shell integration is harmless since modes at different wavelength are independent [5]; hence, only the rescaling step remains, giving to all parameters, fields, and coupling constants their naive (or engineering) scaling dimension. However, in the interacting theory, where nonlinear interactions are present, during the shell integration the coupling between long and short wavelength modes generates perturbative corrections to the parameters of the theory, which have the effect to correct the scaling dimensions. In order to explicitly compute these corrections, one performs a perturbative expansion of the shell integrals in terms of the parameter ε = d c − d, where d c is the upper critical dimension, namely the dimension above which mean-field theory is exact. This expansion method, known as ε-expansion, is nowadays a well established procedure in the context of perturbative RG techniques [5,17,49,50].

A. Renormalization group equations
In this section we will explicitly show how the renormalization group changes the parameters of the model. As stressed out before the RG procedure unfolds in two steps. In the first step we integrate out the small length scale (large momenta) fluctuations, namely fluctuations with momenta Λ/b < k < Λ. The effect of this integration is twofold: i) it changes the cut-off of the theory, since now only fluctuation k < Λ/b are allowed; ii) it changes the value of the parameters of the model, which acquire corrections due to the coupling between low and high momenta fluctuations, We omitted the tensorial structure of the action for easier reading. We remark the fact that all the corrections are proportional to the volume of the momentum shell, which is always proportional to ln b. The standard way to carry out this task, and to compute the corrections to the bare parameters of the model, is using perturbation theory; to be more precise we will compute the corrections δP using a Feynman diagram expansion. The second step consist in re-scaling momenta, frequencies and fields 5 , The effect of this second step is to restore the original cut-off of the theory. Moreover after this step all the parameters of the model acquire a naive scaling factor, corresponding to their naive/engineering scaling dimension: where z is the dynamic critical exponents, which determines how the order parameter relaxes close to the critical point. After these two steps we completed the RG step, and we ended up with a theory defined by a new set of parameters: How the parameters change iterating this procedure defines the renormalization group flow, and the fixed point of this flow rules the critical dynamics of the system.
Even though we apparently have six equations for six parameters, it is possible to reduce the complexity of the problem by introducing an appropriate set of four effective coupling constants and effective parameters, namely, 86) that, through equations (85), are regulated in a closed manner by the following four equations, 6 We point out that if we set λ 0 = λ ⊥ 0 ≡ λ the set of effective parameters (86) coincides with the one of the standard, unconstrained Model G [7]. Once we iterate the RG transformation, we obtain the final recursive RG flow equations, The fixed point value of these parameters will determine the value of the critical exponents of the model. We notice that the naive scaling dimension of the effective couplings f 0 and u 0 is = 4 − d, which suggest that their fixed point will be of order . In the next sections we shall compute all the contributions arising from the shell integral at first order in (one loop), using the Feynman diagrams technique.
Of course, the equations above are rather useless without a determination of the various perturbative corrections δΓ, δg, . . . , which we now calculate.

B. Self-energies
We start computing the perturbative corrections to the Gaussian parameters of the action, namely m 0 , Γ 0 , λ ⊥ 0 , λ 0 . Diagrammatically these perturbative corrections are given by the two-field vertex function, also referred as self-energies: The self-energies Σ and Π are given by the following Feynman diagrams expansion: which corrects the action as follows: Therefore Σ αβ corrects Γ 0 and r 0 , while Π αβγν corrects λ ⊥ 0 and λ 0 . These diagrams can be written as integrals, using the standard Feynman diagrammatic rules, and read 6 In fact, u 0 is not a new effective coupling constant! We are just going back to the original static ferromagnetic coupling constant; we are sorry for this back-and-forth, but it was quite inevitable.
Π αβγν (k) = g 2 0 I αβσµ I ρτ γν wherep + =p+k 2 , whilep − =p−k 2 . The integration over the frequency is performed explicitly, as ω has no cutoff, while the integration in the wave vector is performed by using a thin-shell approximation, which is valid for b ∼ 1; in this way we get: Here we evaluate Σ αβ and Π αβγν a the relevant order in the momenta, namely up to order k 2 . From (96), (97) it is possible to read the perturbative correction to the parameters Γ 0 , m 0 , λ ⊥ 0 , λ 0 :

Absence of spin dissipation
It is important at this point to emphasize a key result: the self-energy Π αβγν (k, ω) vanishes when k → 0. This fact implies that the renormalization group is not generating a dissipative term for the spin, namely it is not generating a linear term in the spin equation of motion which is finite at k = 0. This result is strictly related to the particular structure of the mode coupling vertex of the spin equation of motion (77), and in particular to the fact that this vertex vanishes at zero external momentum k. We believe that this result, which we proved here only at one loop level, could be valid at all orders in perturbation theory. Indeed, the most general diagram that can generate a dissipation is given by, where the blob represents the sum of all 1-particle irriducible diagrams compatible with the given external legs, namely the renormalized mode coupling spin vertex, in which all the possible diagrammatic corrections (at all orders) are taken into account. If the renormalized spin mode coupling vertex vanishes at zero external momentum k, then this diagram is zero too at k = 0, implying that no dissipation is generated. Therefore, as far as the structure of the spin mode coupling vertex is preserved under RG, no spin dissipation can be generated. Even though in this work we explicitly showed that the structure of this vertex is preserved under RG only at one loop level, there are some hints suggesting that this result could remain valid at all order in perturbation theory. It can be shown, by computing the Fokker-Planck equation [7] from the Langevin equations (44) and (45), that the probability density of the system approaches the equilibrium Gibbs-Boltzmann stationary state, P ∼ exp (−H). This non-perturbative result relies on the specific structure of the mode coupling vertices, which means that if the structure of the mode coupling interactions (77) and (75) were different, the system would not have a stationary equilibrium distribution. As this result is not perturbative, it is reasonable to think that the RG does not violate it; indeed, it would be very strange if under coarse-graining an equilibrium model flowed to an out-ofequilibrium one. Therefore, we believe that the structure of the mode coupling vertices is preserved by the RG at all orders in perturbation theory. Finally, since an essential requirement to generate a spin dissipation is that the RG changes the structure of the mode coupling vertices, we believe that no dissipation is generated at all orders in perturbation theory. Note also that this result ensures that in the presence of a small bare dissipation, the same crossover observed in [20,21] between a conservative and a dissipative dynamics is recovered.
Finally, we remark once again that the absence of spin dissipation does not imply that the total spin is instantaneously conserved: due to the presence of the novel Katz vertex (52), which is not zero at k = 0, conservation is broken, even though it is still true that the mean value of the total spin vector is conserved, suggesting that a generalized spin precession occurs: similarly to a standard angular momentum without dissipation but in presence of an external force, which performs a periodic precession with constant time average, so the total spin in presence of the solenoidal constraint has a non-dissipative time dynamics that conserves its mean value.

C. Mode coupling vertex corrections
The corrections to the mode coupling constant g 0 are given by the following vertex functions, These two vertex functions correct the action as follows: We are interested in how these vertex function change the value of the mode coupling constant, hence we can compute Vψ ψs at the zeroth order in the momenta, and Vŝ ψψ at the second order in the momenta (because the g 0ŝ ψψ term in the action is of second order in the wave vector).
A first key consistency check: perturbative renormalization vs symmetry generator The coupling g 0 is the parameter conjugated to the generator of the rotational symmetry (the spin), hence it plays a central role in the definition of the Poisson structure and for this reason it cannot take perturbative contributions from the RG calculation. Therefore, we expect both Vψ ψs and Vŝ ψψ to be zero, and fortunately this is indeed the case in our calculation. From the technical point of view, however, the fact that these vertex functions are zero is extremely nontrivial and it is worth showing, as it is a vital consistency check of the calculation and in particular of the necessity of the new Katz vertex. The vertex function Vψ ψs is given by the following Feynman diagrams: To compute the corrections to g 0 we must compute Vψ ψs at the zeroth order it the momenta. As it happens in the non-constrained case [7], these four diagrams, at zero external momentak,q, sum up to zero.
Vψ ψs αβγν (0, 0) = 0 The vertex function Vŝ ψψ αβγν (k,q), on the other hand, is given by the following nonzero Feynman diagrams: The first two diagrams (the triangles) cancel each other, exactly as in the non-solenoidal case. The other two dia-grams are, on the other hand, specific to the solenoidal case. The first diagram, which vanishes in the nonconstrained theory [21], is nonzero when the solenoidal constraint is present, due the suppression of the longitudinal ψ mode, (108) where the ellipses stand for higher order in the momentum expansion, representing corrections to RG-irrelevant interactions. This vertex corrections not only would give to g 0 a perturbative correction due to the shell integration, but it would generate a novel interaction term too, since it has not the same tensorial structure as the original interaction (77). Therefore, if no other diagram canceling it were present, the RG would not have a closed structure, meaning that the equations of motion would not be eigenstates of the RG transformation, since the shell integration generates new interaction terms that were not present in the bare theory. This unpleasant scenario, in which new relevant interactions arise during the RG flow, is avoided by the key presence of a new Feynman diagram formed by a bubble connection of two lines of the Katz vertex, namely, where, as before, the ellipses stand for higher order in the momentum expansion. The presence of this diagram is fundamental, since it exactly cancels the contributions of diagram (108), therefore curing the anomalies that the latter carries and making the solenoidal RG calculation self-consistent. We therefore consider the following diagrammatic equation a key result of our calculation: Moreover, although the Katz vertex contributes to the dynamic behaviour of s also at k = 0, namely violating the conservation of the total spin, the last diagrams vanishes at k = 0. The fact that the simplest one-loop diagram that can be constructed starting from the Katz vertex does not give any perturbative contribution at vanishing momenta suggests that no spin dissipation should arise, not even at higher orders in the ε-expansion. This is related to the fact that the Katz vertex violates the conservation of the spin leading to a precession of it, but keeping fixed its average value.

D. Ferromagnetic vertex corrections
In contrast to the mode coupling vertex, the ferromagnetic coupling does have perturbative corrections due to the shell integration. The coupling J 0 is corrected by the four-field vertex function Vψ ψψψ , that corrects the action as follows: At one loop, only the following Feynman diagrams contribute to the vertex function Vψ ψψψ in a non-trivial way: The first term is the classic fish diagram of the standard ferromagnetic theory; the second diagram is generated by joining a mode-coupling vertex with the Katz vertex; the last three diagrams are of purely mode-coupling origin.
Computing Vψ ψψψ at the zeroth order in the external momenta, gives, so that the perturbative corrections to the coupling J 0 are,

A second key consistency check: statics vs dynamics
The dynamical RG calculation (at equilibrium) must of course contain in itself the static RG calculation; more specifically, if a coupling constant is present also in the static case, its dynamical renormalization must be exactly the same as its static renormalization. This is clearly the case for the ferromagnetic coupling, which is perfectly welldefined also within a purely static framework. Therefore, in this section we show that this consistency between statics and dynamics is achieved by our calculation.
First of all we recall that the actual ferromagnetic coupling, namely the coupling constant that appears in the static Hamiltonian, is u 0 = J 0 /Γ 0 (see (61)). Hence, the static ferromagnetic coupling u 0 gets perturbative corrections both from J 0 and Γ 0 , and from equations (115) and (98) we have, Hence, the static coupling u 0 does not receive any perturbative corrections from the dynamic coupling g 0 , which is healthy. But the crucial check is whether the recursive relation we get for u from the dynamic RG is the same as the static one, equation (27). Fortunately, it is. From equation (88) and (117) we obtain, which is exactly the same as the static RG recursive equation (27). We stress that this key consistency is recovered in an extremely nontrivial way; in particular, the Katz vertex plays a crucial role. The cancellation of the dynamical coupling g 0 in the perturbative correction of the static coupling is achieved through the following diagrammatic identity: where the l.h.s of the equation is computed at zero external momenta. The diagram D 1 is the product of the interplay between the mode-coupling vertex and the Katz vertex, which is therefore crucial in recovering the correct static behaviour.
There is a second, and subtler, consistency check re-lated to the renormalization of the ferromagnetic coupling constant. The coupling u 0 not only appears in front of the ferromagnetic vertex, but -due to the static-dynamic coupling induced by the solenoidal constraint -it also appears in front of the Katz vertex, Vŝ ψψψψ , which is indeed proportional to g 0 u 0 (see equation (78)); we know that g 0 does not acquire perturbative corrections, hence any diagrammatic correction to the Katz vertex must be billed to u 0 ; but u 0 has been already corrected by its natural ferromagnetic vertex Vψ ψψψ , in the static-compliant way that we have just seen, equation (117). Hence, it seems we have two potentially independent corrections to u 0 , one coming from the bona fide ferromagnetic vertex, and a second one from the Katz vertex! If these diagrammatic corrections were different from each other, we would have a serious problem, as there would be a bifurcation of the ferromagnetic interaction, with highly dubious physical interpretation, not to mention the impossible recovery of the equilibrium static results. Once again, fortunately, the calculation does not disappoint us, even though in a very nontrivial way. The Katz vertex function Vŝ ψψψψ has (at one loop) the following non-vanishing diagrammatic contribution, Computing this Feynman diagram we get, and therefore the correction to u 0 that we obtain from the Katz vertex is exactly the same as from the ferromagnetic vertex, namely δu = −17/2u 0 , which saves the day.

V. THE CRITICAL DYNAMICS OF SOLENOIDAL MODEL G
We now have all that we need to finally calculate the dynamical critical exponent z at one loop in a modecoupling theory subject to a solenoidal constraint, that is in solenoidal Model G.
directions of the fields, both in real and Fourier space, meaning that transverse and parallel modes must be equal and therefore λ = λ ⊥ = λ , namely x = 1. However, at the new stable fixed point of the solenoidal Model G described by Eq. (129) x * = 1/2, meaning that the anisotropy due to the suppression of the longitudinal ψmode leads to a different dynamic behaviour of the two s-modes s ⊥ and s , in such a way that λ * = 2λ ⊥ * . This result directly follows from the fact that the perturbative corrections δλ ⊥ and δλ , given in Eqs. (100) and (101) respectively, obey the relation It is not yet clear to us whether this factor 2 between λ ⊥ and λ can be guessed through a direct analysis of the equations of motion, or if it valid only in the long wavelength and long time dynamic behaviour. In any case, since the diagram contributing to δλ ⊥ and δλ does not involve neither propagators nor correlators of the field s, we may observe that this result is a pure consequence of suppression of the ψ modes. In fact, even if we had naively left λ 0 = λ ⊥ 0 = λ 0 at a bare level, the RG transformation would have led to two different perturbative corrections δλ ⊥ and δλ , meaning that the infrared behaviour of this theory has two different diffusive coefficients for the s ⊥ and s modes.

C. The dynamical critical exponent z
To find the dynamic critical exponent, following [7] and [17], we require that the kinetic coefficient of the primary field, Γ, is non-singular at the RG fixed point, thus ensuring that the effective RG theory has a non-singular characteristic time scale. This amounts to imposing the condition, and thus we need to write explicitly the recursive RG equation for the kinetic coefficient; that equation can be found in (85), complemented by its perturbative corrections, equation (98), thus giving, By using the obvious expansion, b x (1 + y ln b) = b x+y , we finally obtain from (131) and (132) that the general expression for the dynamic critical exponent is, and once we use the one-loop values of the parameters at the stable fixed point (129), we obtain Despite the difference of all effective parameters and coupling constants at the stable fixed point, this is exactly the same dynamic critical exponent as the standard unconstrained Model G [7]. This result is somewhat surprising. The solenoidal constraint does change the static universality class: the static critical exponents are different from the Landau-Ginzburg class, and define the novel dipolar ferromagnet class (see Table I and [24]). According to the common wisdom in the theory of critical phenomena, we would expect a change also in the dynamic universality class, as universality is normally broader at the static level than at the dynamic level. For example, Ising-like ferromagnets, have the same static critical exponents, while the dynamical critical exponents varies depending on whether the order parameter is conserved (Model B) or not (Model A). Here, we see something different: the dynamic universality class does not change due to the solenoidal constraint. Even though we have derived this conclusion perturbatively, this is probably a non-perturbative result due to the great power of the symmetry; the lack of diagrammatic renormalization of the coupling constant conjugated to the generator of the rotations, g 0 , leads to the recursive relation, If we now ask that the central charge of the symmetry does not change at all, we get z = d/2, at the non-perturbative level.
The fact that the solenoidal constraint does not alter the dynamical critical exponent is an encouraging result, since it suggests that using the incompressibility condition to simplify the dynamical equations of active matter with mode coupling interactions is a reasonable approximation, as it does not change dramatically the dynamical behaviour of the system.

VI. CONCLUSIONS
We have studied the effects of a solenoidal constraint on the critical dynamics of a field ψ with O(d)-symmetry in the presence of mode-coupling interaction with the generator of the rotational symmetry s, which we called spin; more succinctly, we have studied solenoidal Model G. The presence of the constraint leads to the suppression of the ψ-mode parallel to the wave-vector k, namely ψ (k) = k · ψ (k) /|k|, violating the O(d)-symmetry and modifying the static behaviour. The equations of motion of the constrained theory have been derived starting from the symmetries and Poisson-bracket relations between the hydrodynamic variables, namely the order parameter ψ and the spin s. We performed a one-loop renormalization group calculation to investigate the long wave-length and long time behaviour in the critical region. The closed structure of the RG transformation and the consistency of the RG flow with the static behaviour of dipolar ferromagnets provide a self-consistent proof that no RG-relevant interaction has been omitted and that the equations of motion we derived are correct.
Two main dynamic effects arise as a consequence of the solenoidal constraint. The first, and most predictable one, is the projection on the plane orthogonal to k of the equation of motion for the order parameter ψ. On the contrary, no similar projection of the equation for s can be performed in order to obtain the spin dynamics; instead, the suppression of ψ (k) leads to a novel nonlinear interaction -the Katz vertex -combining the effect of the static ferromagnetic coupling of the field ψ and the mode-coupling dynamic interaction. The presence of this new mixed interaction is the second, less intuitive, effect of the constraint, which adapts the spin dynamics to the presence of the constraint by making the static quartic coupling contributing to the torque-like interaction ∂ t s ∼ g 0 ψ × δ ψ H. Moreover, the Katz vertex contributes to the time-derivative of the spin also at zero wave vector, therefore violating the conservation of the total spin.
The lack of conservation of the spin is not something strange: the order parameter is not rotational-invariant, due to the solenoidal constraint, and therefore the generator of its rotations is not a conserved quantity. It is however crucial to understand whether the spin is dissipated or not, since the presence of a dissipation generated by the RG would make the spin stop being an hydrodynamic variable, therefore suppressing any inertial behaviour in the critical region. The torque-like nature of the Katz vertex, which is the only dynamic term violating the spin conservation, indicates that this violation gives rise to a generalized precession of the total spin, rather than a dissipation. At one loop, the perturbative expansion confirms this interpretation, since the self-energy of the spin Π does not contain any perturbative corrections at k = 0. Moreover, we show that the presence of any dissipative term in the linear dynamics of the spin can arise only if the dynamic mode-coupling vertex in the equation of s did not vanish when k = 0, which seems not to be the case for this theory.
Our RG calculation passed several nontrivial consistency checks. First of all, the fact that the equations of motion appear to be eigenstates of the RG, in the sense that the shell integration step does not generate new interaction terms, ensures that we did not miss any relevant coupling in the problem description. Secondly, the RG recursive relations we found for the dynamic theory reproduce the behaviour of dipolar ferromagnets, ensuring that the static behaviour is correctly reproduced by the dynamics. We must remark that both these results directly follow from the presence of the new Katz vertex, in absence of which the theory would not describe correctly the dynamics of the system; therefore meaning that the this new non-trivial interaction plays a crucial role in making the dynamic behaviour of the spin compatible with the constraint.
From the study of the RG recursive relations, the dynamic behaviour has been shown to be characterized by a critical exponent z = d 2 , which is the same as the nonconstrained theory. This result was somewhat unexpected. In general, static properties are more robust compared to dynamical properties; models with different dynamical critical behaviours are often characterized by the same static behaviour, such as Model A and Model G of [7]. On the contrary, in our model the solenoidal constraint changes the static universality class, leaving unchanged the critical dynamics of the system. This suggests that the dynamic critical behaviour of homogeneous systems does not change when a solenoidal constraint, i.e. incompressibility, is enforced; this indicates that we can try and understand their dynamic critical behaviour by studying their incompressible version. Homogeneous systems are governed by equations of motion in which incompressibility is, in general, not required, but where density fluctuations, and therefore density-velocity couplings, are negligible. Incompressibility, though, completely suppresses density fluctuations and therefore represents a stronger requirement than homogeneity. Moreover, requiring incompressibility, hence imposing a solenoidal constraint on the velocity field, generates long-range interactions that could change the properties of a system; while this is indeed the case for the static behaviour of our theory, the long-range interactions are not sufficient to modify its dynamic universality class. Theoretical evidences of this fact have already been discussed in homogeneous active systems [35], therefore suggesting that the solenoidal constraint does not significantly affect the critical dynamic behaviour in the presence of neither activity nor mode-coupling terms. This result is very encouraging, as it allows to study the homogeneous phase of the offequilibrium Inertial Spin Model under a incompressible hypothesis, where the absence of the density field leads to a great simplification of the calculation.
Our result is an important stepping stones towards developing an RG theory for natural swarms. The Katz interaction vertex in the spin dynamics that we derived in this work, which emerges as an effect of the solenoidal constraint, must certainly characterize also the equations of motion of an incompressible out-of-equilibrium field theory, in which terms coupling the order parameter to its generator of rotations are present. It would have been extremely difficult to derive the Katz vertex had one tackled the problem directly in the presence of activity. Despite this step forward, though, the complexity of the present calculation suggests that significant more theoretical efforts will be needed to carry out the full-fledged out-of-equilibrium mode-coupling RG study of natural swarms.
We warmly thank M. Testa, C. Castellani, and U. Tauber for fruitful discussions. L. Di Carlo also thanks I. Basile for interesting exchanges about the work. This work was supported by ERC Advanced Grant RG.BIO (contract n.785932) to AC, and ERANET-CRIB Grant to AC and TSG. TSG was also supported by grants from CONICET, ANPCyT and UNLP (Argentina).