Field Theories with Conformal Carrollian Symmetry

Conformal Carrollian groups are known to be isomorphic to Bondi-Metzner-Sachs (BMS) groups that arise as the asymptotic symmetries at the null boundary of Minkowski spacetime. The Carrollian algebra is obtained from the Poincare algebra by taking the speed of light to zero, and the conformal version similarly follows. In this paper, we construct explicit examples of Conformal Carrollian field theories as limits of relativistic conformal theories, which include Carrollian versions of scalars, fermions, electromagnetism, Yang-Mills theory and general gauge theories coupled to matter fields. Due to the isomorphism with BMS symmetries, these field theories form prototypical examples of holographic duals to gravitational theories in asymptotically flat spacetimes. The intricacies of the limiting procedure lead to a plethora of different Carrollian sectors in the gauge theories we consider. Concentrating on the equations of motion of these theories, we show that even in dimensions d=4, there is an infinite enhancement of the underlying symmetry structure. Our analysis is general enough to suggest that this infinite enhancement is a generic feature of the ultra-relativistic limit that we consider.


Effective theories and singular limits
Physics is all about length scales, and most theories that we use to describe Nature are effective theories, including the spectacularly successful Standard Model of particle physics. Effective theories are often formulated by going to a certain sub-sector of a higher theory. Hence, often there is a lot to be learnt by focussing on a sub-sector of a particular theory. One of the very obvious examples that come to mind in this regard is the non-relativistic sector of well defined relativistic theories, which become important for most of the physics around us. Galilean physics is obtained from relativistic physics by taking the speed of light to infinity. Group theoretically, the symmetries contract from the relativistic Poincare group to the Galilean group as one takes c → ∞. Various counter-intuitive things happen in this singular limit. The light-cones open up, the spacetime metric degenerates and the Riemannian structure of spacetime manifold is lost. New mathematical structures emerge, viz. the spacetime manifold becomes what is called a Newton-Cartan manifold (for a brief introduction see e.g. [1]), and often the underlying symmetries get enhanced. Some of these ideas have come into light only recently and this is a very fertile area of current research [2][3][4][5][6][7][8][9][10][11], some of which have been motivated by understanding the notion of holography for non-relativistic systems [12][13][14][15][16].
Out of mathematical curiosity, one could ask if similar new structures may emerge when one takes the speed of light to go to zero instead of infinity. This peculiar c → 0 limit has been dubbed the Carrollian limit in literature [17]. This defines a contraction of the Poincare group different from the Galilean one, and the group so generated is called the Carrollian group. In this paper, we would be interested in constructing field theoretic examples that exhibit Carrollian symmetry. A selection of previous related work in this direction includes [18][19][20][21][22].
Purely mathematical curiosity is perhaps not very good justification for research in theoretical physics, although many very important breakthroughs, like the formulation of Yang-Mills theories that later have led to the understanding of the weak and strong nuclear forces, have come via this path. Our investigations in this paper are driven by some very strong motivations, which we elaborate on below.

Holographic duals to flatspace
It is more than 20 years since the advent of the famed AdS/CFT correspondence [23], which has given a firm footing to the idea of holography [24,25]. This has opened up an avenue of studying quantum gravity by looking at quantum field theories and vice versa. Our physical world is, however, clearly not AdS. For many applications, especially astrophysical ones, the universe can be well approximated by an asymptotically flat spacetime. It is thus of great importance to extend the notion of holography from its original setting in asymptotically AdS spacetimes to flat backgrounds.
A natural way to construct a holographic quantum field theory for a general gravitational theory is to consider the symmetry structure at the boundary of the spacetime in which the gravitational theory lives. These symmetries at the boundary are formally given by the Asymptotic Symmetry Group (ASG) and its associated algebra, the Asymptotic Symmetry Algebra. In mathematical terms, given a set of boundary conditions, Asymptotic Symmetry Group = Group of all allowed diffeomorphisms Group of trivial diffeomorphisms . (1.1) One can then propose that the dual field theory lives on the asymptotic boundary of the spacetime and inherits the symmetry of the ASG. This analysis of the ASG of AdS 3 famously led Brown and Henneaux [26] to two copies of the Virasoro algebra which are of course the symmetries of a 2d conformal field theory. This could be looked upon as the principle precursor to the discovery of AdS/CFT. Recently, similar asymptotic symmetry analysis for theories of higher spin in AdS 3 [27] has led to the uncovering of higher spin dualities [28]. Earlier, similar ideas were used to propose a holographic dual to de Sitter spacetimes [29]. We wish to extend this analogy to the case of asymptotically flat spacetimes. For Einstein gravity in 3 and 4 dimensional Minkowski spacetime, the asymptotic symmetries at the null boundary of spacetime are given, not by the Poincare group, but by the infinite dimensional Bondi-Metzner-Sachs (BMS) groups. The associated algebras are The structure at null infinity in 3d is IR × S 1 . L n form the diffeomorphisms of the circle at I + , while the M n 's are angle-dependent translations called supertranslations. The central charge 3/G to the [L, M ] commutator is the value for Einstein gravity [30]. One could in general have a non-zero central extension to the [L, L] commutator. For Einstein gravity, this is zero. But one could add a topologically massive term to generate a non-zero central extension for theories of gravity with higher derivative interactions [31].
In 4d, I + becomes IR × S 2 . The structure on S 2 enhances to two copies of the Witt algebra, very much like 2d CFTs [32] and the algebra becomes One cannot find physical boundary conditions (e.g. without excluding gravitational radiation and the Schwarzschild black holes in 4d) so as to limit these infinite dimensional groups to the finite Poincare group in dimensions 3 and 4 1 . 1 It is important to mention here that another infinite dimensional extension of the BMS4 exists where the infinite dimensional group is taken to be the semi-direct product of the supertranslations with Dif f (S 2 ), i.e. the group of all smooth diffeomorphisms on the conformal sphere S 2 [33]. We shall however not be interested in this or its possible generalisations (or the lack of it) to higher dimensions. This is however not true in d > 4. Here with boundary conditions strict enough one can reduce the ASG to ISO(d − 1, 1) [34][35][36]. Interestingly though, one can find looser boundary conditions so that there are infinite enhancements in d > 4 [37]. This begs the question as to which set of boundary conditions is the more physically relevant.
Of late, there has been a resurgence in the study of physics in flat-spacetimes, initiated by Strominger and his collaborators. An incomplete list of interesting work in this direction includes [38][39][40][41][42][43][44][45][46][47][48][49][50]. A very intriguing picture has emerged linking asymptotic symmetries with soft theorems in QFT in asymptotically flat spacetimes and to memory effects. The interested reader is pointed to [51] for a review of the current developments in this field. The asymptotic BMS symmetries arise at null infinity and it is thus very natural that these should play an important role in scattering theory in flat space. Weinberg's soft graviton theorem [52] relates the S-matrix element of a theory of quantum gravity to another which differs from the original one by the addition of a soft graviton. Recently it has been shown that this soft-graviton theorem actually arises out of a supertranslation Ward identity [40]. Weinberg's theorem does not depend on the dimension of spacetime where the quantum theory of gravity lives. So, if we were to assume that the relation between asymptotic symmetries and soft theorems exist in all dimensions, then these infinitely extended supertranslation symmetries should also exist in all dimensions. It is thus perhaps of greater physical significance to consider looser boundary conditions in d > 4 which generate an infinite dimensional asymptotic symmetry algebra.
If we are to draw inspiration from AdS/CFT, holography in asymptotically flat spacetimes should involve these infinite dimensional symmetry algebras and the putative dual theories should be non-gravitational quantum field theories living on the null boundary and invariant under the infinite extended BMS algebras. A natural avenue to explore flat holography is to investigate the singular limit where the bulk theory goes from AdS to flat space, i.e. taking the radius of AdS to infinity. It can be shown that this leads to an ultra-relativistic contraction of the boundary CFT [53]. So these conformal versions of Carrollian theories that we discussed in the beginning are putative duals of flat space. It has also been shown that the conformal Carrollian symmetries are isomorphic to BMS symmetries [54,55].
A relativistic CFT in d dimensions has a finite symmetry algebra, viz. so(d, 2). The ultra-relativistic contraction of this also obviously leads to a finite symmetry algebra iso(d, 1). This, clearly, is also the process of taking radius of AdS d+1 to infinity to get to flat space, which has (d + 1) dimensional Poincare symmetry. Even before speaking of infinite extensions, it is clear that these d-dimensional field theories are exotic theories. For example, for asymptotically flat spacetimes in 4 dimensions, we are seeking a field theory in d = 3 that is invariant under iso (3,1). These are clearly not usual relativistic field theories and as we have stated just above, one of the natural ways to generate them is to consider Carrollian limits of relativistic CFTs.
The case of d = 2 is special. Here, as is very well known, the relativistic conformal algebra enhances to two copies of the infinite dimensional Virasoro algebra, which under the Carrollian limit, generate the ultra-relativistic algebra (1.2): L n = L n −L −n , M n = L n +L −n , → 0. (1.4) Here is the speed of light in the 2d theory. In holographic terms, the identification is as follows: = G/ , with being the AdS radius and G is the Newton's constant. This is thus also the flatspace limit. The ultra-relativistic contraction and its link to the large radius limit of AdS is explored in detail in [53]. The symmetries in both the parent AdS theory and its contracted flat space version remain infinite dimensional and there is a lot to be learnt by exploiting this singular limit. This has been used to understand features of flat holography in 3 bulk and 2 boundary dimensions. An incomplete list of references include [31,[53][54][55][56][57][58][59][60][61][62][63][64][65][66][67][68]. The interested reader is pointed to [19,69] for a summary of work in this direction. Some of the more recent interesting developments include [70][71][72][73].
In boundary dimensions d > 3, the story is very different. The limit from AdS yields just finite dimensional symmetries. But, in flat space, we need more,-we need to generate infinite dimensional supertranslation symmetries in order to relate to Weinberg's theorem. Starting out with finite symmetries on the boundary via the process of contraction that we have just outlined above, we thus need a mechanism to generation these infinite supertranslations. Following earlier work [19], we describe an algebraic process which achieves this. We will focus specifically in d = 4 and find that generic field theories generated by this process of contraction naturally have these infinite dimensional symmetries that emerge in this singular limit.

Other potential applications
Taking the speed of light to zero focuses us on a subsector of a parent theory that is very highly energetic. One could think of doing this on the worldsheet of the relativistic string. This singular limit, interestingly, is also the tensionless limit of string theory [74] where a greater symmetry structure has been conjectured to exist and this has been in focus for its relation to recent higher spin holographic theories [75]. The residual symmetry on the tensionless string worldsheet after fixing the equivalent of the conformal gauge turns out to be (1.2) [76][77][78]. The construction of tensionless strings has been recently revisited in connection to these symmetry structures [79][80][81]. It is very likely that when one is interested in taking the tensionless limit of membranes, higher BMS or equivalently conformal Carrollian field theories would have a decisive role to play [78].

Outline of this paper
In this paper, we systematically construct Carrollian field theories which are further invariant under the Conformal Carrollian group or equivalently the BMS group (in one higher dimension). We start from well-known relativistic conformal field theories and show that in the ultra-relativistic limit, the contraction of these field theories become invariant under the infinitely-extended conformal Carrollian group. In Sec. 2, we begin with a review of the algebraic aspects of conformal Carrollian symmetry. We discuss the process of contraction that gets one from the relativistic conformal algebra to first the finite conformal Carrollian Algebra (CCA) and then proceed to give the CCA an infinite dimensional lift. We then provide a detailed survey of the representation theory of this algebra.
In Sec. 3, we begin our construction of field theories invariant under the CCA. We start with scalar fields and outline the process of finding the symmetries of the limiting field theory by focusing on equations of motion (EOM). We are thus discussing on-shell symmetries throughout this work. We then move on to fermions and then to the Carrollian version of Yukawa theory. We focus on dimensions d = 4. In each case, the study of EOM reveals an infinite symmetry enhancement of the underlying theory to the (infinitely) extended conformal Carrollian group.
In Sec. 4, we revisit the construction of Carrollian electrodynamics and the emergence of the Electric and Magnetic limits here [18,19]. We then add matter to these theories. We showcase the example of Carrollian scalar electrodynamics with the most general scalings of the (massless) matter fields and show how consistency of the theory reduces the allowed parameter space of scalings to a particular patch in each of the Electric and Magnetic limits. All these allowed values of scalings give rise to EOM which are invariant under the infinite CCA. We show this for a theory with an arbitrarily chosen set of parameters. We then add fermions to the U (1) theory and perform a similar analysis and find infinite symmetries of the EOM.
In Sec. 5, we move on to non-Abelian theories with Carrollian invariance. First we review pure Carrollian Yang-Mills theory following earlier work [19]. We then go on to add fermionic matter to these theories and investigate the various limits arising from the parent theory. The (gory!) details of the SU (2) theory are worked out in Appendix A. In the main body of the paper, we present the general analysis of SU(N ) Carrollian Yang-Mills theory coupled to fermionic matter. We end up with numerous possible sectors in the UR limit which we systematically analyse. A number of these sectors can be discarded for the absence of kinetic terms in some equations of motion. We again obtain that in all of these sectors, the symmetries of the EOM get enhanced infinitely and the CCA turns out to be the relevant symmetry algebra. This enhancement of symmetries in the ultra-relativistic limit of conformal field theories thus seems to be generic and we propose that this would be true at the classical level for all dimensions and all conformal field theories.
We conclude with a summary of our results and a number of possible directions of future work. Appendix B contains details of the cousin of the Carrollian Yang-Mills theory, viz. the Galilean SU(N ) Yang-Mills theory and associated additional fermionic matter. We present the similarities and differences between the two theories in Appendix C.

Conformal Carrollian Symmetry
We begin with some algebraic preliminaries. We will first consider how to obtain conformal Carrollian algebras from relativistic conformal algebras and then go on to a detailed review of the representation theory which we would use throughout the rest of the paper.

Contraction and algebra
The ultra-relativistic limit of a relativistic conformal field theory is obtained by performing an Inönü-Wigner contraction on the generators of conformal group. Let us consider a d dimensional CFT. The limit that one needs to take in terms of the spacetime coordinates is: . This means we are sending the speed of light, c → 0. Under this limit, the relativistic conformal algebra contracts and the generators of the contracted algebra, which is the finite conformal Carrollian algebra (CCA), are written as: The non vanishing brackets of the finite CCA generators are given by: where J ij , H, P i , B i are the rotation generators in (d − 1) spatial dimensions, time translation, spatial translations, and Carrollian boosts respectively. These form the Carrollian group, which is the ultra-relativistic analogue of the Galilean group. The conformal generators are dilatation D, and the Carrollian versions of temporal and spatial conformal generators K, K i . The set of the generators {J ij , P i , D, K i } close to form a so(d − 1) subalgebra of the finite conformal Carroll algebra, which unsurprising is iso(1, d − 1). This so(d − 1) subalgebra can be looked upon as the conformal generators of the S d−1 at the null boundary of flat spacetime, which as stated earlier, has a structure IR × S d−1 .
In [19], it was shown that it is possible to give the CCA infinite dimensional extensions. These are different for different dimensions. The infinite extensions in the d = 2 case yields (1.2) and for d = 3 gives (1.3). Here we are interested in the d = 4 construction. We have three spatial directions and one temporal one. Of course, for d ≥ 4, there is no infinite enhancement of the (relativistic) conformal symmetries on the sphere at infinity like in the 2d and 3d examples. The proposed infinite extension for the generators is thus only in the supertranslation part: The finite part of the algebra described above fit into this scheme as follows: The infinite part of CCA algebra are given by the lie brackets [19,22]: The other commutators (e.g. for the other components of P i or K i ) follow similarly. The supertranslations form an abelian sub-algebra among themselves: Generalisation to higher dimensions: The above construction has a natural generalisation to dimensions d > 4. Again, we would look to give an infinite lift to only the supertranslation generators: The conformal structure on the S d remains finite dimensional and generated by the so(d) algebra mentioned above. The algebra of generators follows immediately and is a direct generalisation of the d = 4 case: The infinite supertranslations can also be written as: and the above algebra can be written in a compact form for d ≥ 4 [22]:

Representation Theory
We now discuss the building blocks of infinite conformal Carrollian algebra, the representation theory based on highest weights. Unless otherwise stated, our construction would be for d = 4. As with the algebra described above, it is expected that the representation theoretic aspects would also have a natural generalisation to higher dimensions. This subsection is primarily intended to be a detailed review of earlier work [19] along with an extension in case of different integer and half-integer spins. We aim to study the conformal Carrollian limit of different free and interacting field theories in d = 4. Each of these field theories can be classified based on their spins. Hence, a natural choice of constructing the representation theory is to label the states by the eigenvalues of dilatation and rotation operators. Thus, the states will be labelled with definite scaling dimension ∆ and spin j. This particular choice of representation is known as the scale-spin representation. Below we describe in detail the scale-spin representation of spin 0, spin 1 2 and spin 1 fields.
Scale-spin representation: As just stated above, we would label the states by their scaling dimension ∆ and spin j.
where, J 2 is the sum of square of the three rotation generators in three spatial dimensions. We propose that similar to 2d CFT, in these conformal Carrollian theories, there exists a state-operator correspondence: This is not strictly necessary, but helps us in translating from operator language to states. It follows that In (2.15) the generators have been promoted to operators and the brackets imply commutators henceforth. Next, we introduce the notion of primaries in a way similar to CFTs. We demand that the spectrum of these field theories (i.e. the value of ∆) is bounded from below. Following [19], the conformal Carrollian primary is defined as an operator annihilated by all lowering operators. This means for the finite algebra: where Σ ij is the spin operator. The spatial and time translation of a generic field at any spacetime point is generated by the Hamiltonian H and momentum operators P i as , it is clear that the primaries are not eigenstates of Carrollian boosts B i . Hence, following previous work [19] we use the Jacobi identity to examine the action of boosts on the primaries: The most general transformation is Here, ϕ, {φ, χ}, {A t , A k } are primaries of different spins and r, {f, f }, {a, b} are some constants yet to be determined. For the purposes of this paper, we will truncate the above expansions at spin 1. In principle, we could have mixing with higher spins as well.
Now we review the action of the finite and infinite CCA generators on the primaries at any general spacetime point. A CCA primary Φ(t, x) is related to its own version at the origin by Hence, for any finite CCA generator O We employ Baker-Campbell-Hausdorff formula along with the algebra in (2.3) to recast the above expression for finite and infinite generators in a simplified form as The equations (2.24) are the action of finite and infinite symmetry generators on a CCA primary. Below we explicitly write down the transformations for scalars, fermions and spin-1 vector bosons.
The representation of CCA for scalars is given by {∆, j, r}. These would be fixed by input from dynamics, e.g. when we will demand that the Carrollian scalar theory emerges as a limit of the relativistic scalar field theory, as we do in the next section, we would be able to fix the value of r to be zero.
Again, as stated before, merely invariance under the underlying symmetries is not enough to fix these values. We would need input from dynamics.
Vectors: The representation is given by the quadruple {∆, j, a, b} Having detailed the features of the scale-spin representation, we are now in a position to deal with conformal Carrollian field theories. Our construction of these theories in the next few sections will crucially depend on the algebraic details outlined in this section. We will start with relativistic field theories and through a process of contraction, generate Carrollian field theories. We will then use the representation theory developed in this section to examine the symmetries of the equations of motion of these Carrollian theories.
We should emphasise that all Carrollian theories need not to be constructed in this way and there should be an independent way of arriving at these peculiar field theories without recourse to any limiting procedure. The advantage of constructing theories in this particular way from the ultra-relativistic limit of relativistic conformal field theories is that the values of the constants undetermined by symmetries that we discussed in this section would be fixed by the limit.
-11 -3 Carrollian field theories: spin 0 and spin 1 2 In this section we will begin our analysis of constructing conformal Carrollian field theories for spin-0 and spin-1 2 fields. We will be looking at the conformal Carrollian versions of free scalars, free fermions and the interacting Yukawa theory.

Scalar Fields
The simplest relativistic field theory is of course the theory of a massless spin-0 real scalar field. The relativistic Lagrangian of massless scalars and the corresponding equation of motion are invariant under conformal transformations in all dimensions. While dealing with Carrollian field theories, we would also be initiating our discussion with the spin-0 field. We wish to see whether the finite conformal Carrollian symmetry remains in the ultra-relativistic limit of real scalar. Afterwards, we would extend our analysis in search of the infinite enhancement of conformal Carrollian symmetry.
The equation of motion of a relativistic massless real scalar field is We take appropriate spacetime scaling to obtain the ultra-relativistic limit of massless scalar field.
Plugging eq(3.2) in (3.1), the Carrollian scalar equation of motion can be expressed as: We will now find the symmetries of Carrollian scalar fields at the level of its equation of motion.
The philosophy that we would be employing to understand the symmetries of the equations of motion of a particular theory is the following. Suppose we have a generic equation: where D is some differential operator acting on a field Φ(t, x) of arbitrary spin, and J is a source term. We wish to examine the symmetries of this equation under a symmetry group G, whose associated algebra has generators Q α . We will say that the equation is invariant under this symmetry when: where ε α is the infinitesimal parameter of the particular transformation.
Fixing the values of the representation theory from the limit We now need to fix the values of the constants {∆, j, r} which remained unfixed in our symmetry analysis by input from the dynamics. We will follow the limit from the relativistic massless scalar theory in order to do this. The relativistic field Φ had a scaling dimension The relativistic generator does not change in form in the ultra-relativistic limit. So its action on the field is also unchanged. We thus will work with the input We are dealing with a scalar field. So, of course, we have j = 0.
For fixing the boost labels, we start with the relativistic boost generator The action of the relativistic boost operator on a generic relativistic field Φ is where Σ i0 describes analogue of the spin for the spatial rotation. As, we are dealing with a scalar field, Σ i0 = 0. When we take the ultra-relativistic limit, as stated in eq(3.2), we get Hence, the action of Carrollian boost generator on a Carrollian scalar field is obtained as, From the representation theory, the action of Carrollian boost on Carrollian scalar fields using equation (2.25) can be written as Comparing the above equation with (3.11), we fix the value of boost label 'r' for the Carrollian scalar field ϕ to be zero.

Symmetries of the EOM
Under finite conformal Carrollian generators, the scalar field equation transforms as: Hence, real scalar field respects the finite conformal Carrollian symmetry. Furthermore, it also has the infinite conformal Carrollian symmetry in d = 4. 14) It can also be checked that the Carrollian scalar obeys conformal Carrollian invariance in all dimensions.

Fermionic fields
In this section, we will be looking at the ultra-relativistic limit of spin-1 2 field or free massless Dirac fermions in 4 dimensions. The Dirac fermion Ψ can be written as two component Dirac spinors φ and χ: We will be working in the Pauli-Dirac representation in which the gamma matrices are given by The σ i 's are the Pauli matrices. The relativistic free Dirac equation is In terms of the spinors (φ, χ) this becomes

Carrollian scalings for Fermions
For the ultra-relativistic massless fermionic theory, we scale the spinor fields as Taking the scaling (3.18) on the relativistic equations of motion, the ultra-relativistic Dirac equation can be recast as, We could have considered a general scaling for the spinors: In this scaling the ultra-relativistic Dirac equation takes the form We wish to see an interaction term apart from the usual kinetic term in the EOM. This can be achieved if s + 1 − r = 0, or, r + 1 − s = 0. (3.22) We have used the second condition, without the loss of generality This condition justifies the scaling in (3.18), if we put r = 0. Considering, any higher value of r will also result in the same Carrollian Dirac equations (3.19). It is worth mentioning here that choosing the first condition in (3.22), will give the Carrollian Dirac equations where the role of φ and χ will be exchanged among each other in (3.19).

Fixing coefficients
Now, we will fix the values of the undetermined constants {∆, f, f } by taking the ultrarelativistic limit on the parent relativistic dilatation and boost transformation. We start with (3.9) and write down the action of Σ i0 on the relativistic fermion Ψ.
where we have lowered the indices in the final result. Also, we mention here the action of relativistic spin operator Σ ij on Ψ, Here, σ i 's are the Pauli Matrices and [σ i , σ j ] = 2i ijk σ k . For a generic spacetime point, we obtain the transformation of the relativistic spinors under relativistic boost as in (3.9) , the transformation of the spinors under Carrollian boost becomes: Taking input from the dynamics (3.28), we are able to fix the values of the boost labels in the representation theory to {f = 0, f = − 1 2 }. Similarly, following the same argument as in (3.6) and (3.7), we fix the scaling dimension ∆ for the Carrollian spinors to be For 4 dimensions, ∆ takes the value of 3 2 .

Symmetries of EOM
Next, we are interested in checking the invariance under conformal Carrollian algebra. Under scale transformation, Similarly, under the spatial part of the Carrollian special conformal transformation (SCT) Here, we rewrite the representation of the massless Carrollian fermion theory for conve- It is thus clear that the equations are invariant under spatial SCT, for these particular values. It can be easily checked that the equations are invariant under the temporal SCT. Under M m 1 ,m 2 ,m 3 , the equations of motion are trivially invariant.
So, we have seen that the EOM for Carrollian fermions are also invariant under the infinite CCA. It should be emphasised here that we have focussed on d = 4. The analysis is expected to hold for all higher dimensions as well.

Yukawa Theory
In the previous sections, we described in detail, the theory of Carrollian version of the free spin-0 and spin-1 2 fields. We are now equipped to analyse ultra-relativistic interacting theories. In this section, we will construct the Carrollian version of Yukawa theory.
We will consider relativistic massless Yukawa theory and then consider its ultrarelativistic limit. The parent theory consists of a scalar field ϕ and a fermion Ψ. Massless Yukawa theory is classically conformally invariant in d = 4. The action for relativistic massless Yukawa theory is The equations of motion are We will now decompose the Dirac spinor into (φ, χ) and use Pauli-Dirac representation of gamma matrices. The relativistic EOM can be written down in terms of (φ, χ) as: Now we want to see the ultra-relativistic limit of massless Yukawa theory and check for its symmetries.

Carrollian scalings for Yukawa theory
Along with the usual coordinates scaling, the scalar and fermion are scaled as We consider this general scaling in order to make sure that we are not missing out on interesting interaction terms by restricting our scalings. Plugging the scaling (3.37) in the equations of motion (3.36) we obtain the scaled equations: For consistency of the theory, when the fermions are turned off, the EOM (3.38) need to reduce to free Carrollian scalar EOM and similarly, when the scalar is turned off the EOM should reduce to the free fermionic EOM. These put restrictions on the allowed values of (α, β) in our scalings above. The constraints on α, β are summarised in Table[1]. Combining all the constraints, we finally arrive at the following relations: Hence the possible values of α, β are Interestingly, the parameter space of different values of α, β splits into two distinct cases: This yields an interesting interacting theory. We will be primarily interested in this and will call this Carrollian Yukawa theory.
• Case 2: (Reproduces free theory) (α = −1, β = α + 1): This just gives back the free equations of the scalar and fermions. Since there are no interaction terms, we will not be interested in these class of scalings.

Symmetries of EOM
As just mentioned, we will choose α = −1, β = 0. Now we investigate the symmetries for the equations (3.40). For this, we use the values of the constants of the individual free scalar and fermionic theories 2 . They are given as The transformations of the EOM under scale transformation are given by Similarly, under special conformal transformation (K l ), we have So we have constructed our first example of an interacting conformal Carrollian theory in this paper and have checked that its EOM have infinitely enhanced symmetries for d = 4.
2 Since we are concerned with the classical invariance of the EOM in this work, interaction terms would not affect the values of these constants and hence using the constants obtained from the free theories is a justified choice. 3 Invariance under K is trivial.
-18 -Having investigated Carrollian scalars and fermions, the natural next step is to consider spin-1 theories. We begin with a brief review of Carrollian electrodynamics [18,19] and then we go on to add massless scalars and fermions to the U (1) theory and investigate its different sectors and the associated symmetry structure of the equations of motion in each sector, like in the previous section.

Carrollian Electrodynamics: A quick look back
This section is intended to be a quick review of the earlier work [19] before stepping onto the next part of this paper. The Carrollian limits of Abelian (U (1) Electrodynamics) and Non-abelian (SU (N ) Yang-Mills) gauge theories were discussed in detail in [19]. We will summarise the details of Electrodynamics in this current section. The relativistic Electrodynamics and Yang-Mills theory are classically conformally invariant in d = 4. The starting point of studying gauge fields in ultra-relativistic limit concerns breaking the Lorentz invariance by treating the spatial and temporal part of a four vector differently, in addition to the scaling of the underlying spacetime. The corresponding scaling results in two distinct sectors: The existence of two different Carrollian sectors of electrodynamics can be traced to the fact that there is no non-degenerate metric for the entire Carrollian spacetime. For this precise reason, the covariant and contravariant vectors form two different inequivalent characterisation of the spin-1 theory. On taking the UR limit on the spacetime, contravariant vectors transform as We associate this with the magnetic sector (4.1b): On the other hand, covariant vectors transform as We can write the transformation of the covariant vectors (after an overall multiplication by a factor of which does not change anything) as This is the electric sector. The dynamics of the two sectors are governed by different equations of motion. These EOM for the electric and magnetic sector are respectively given as: As we have done previously, we can check for the symmetries of the EOM in these sectors, based on the representation theory discussed in Sec [2.2]. We will also need input from dynamics, which as before, can be obtained by looking carefully at the limit from the relativistic parent theory. Both sets of EOM turn out to have finite and infinite conformal Carrollian symmetry. For the detailed analysis of the symmetries of the EOM, the reader is pointed to [19]. We will use this knowledge of the pure gauge theory to go beyond and study Carrollian gauge fields coupled to matter in the next section.

Adding scalars
Scalar electrodynamics is the theory of U (1) gauge field coupled to a complex charged scalar field. The Lagrangian density of the massless scalar electrodynamics theory is given by where F µν is the field strength defined as with e being the coupling parameter and A µ being the U (1) gauge field. ϕ is the massless complex scalar field charged under the U (1) gauge group.
The gauge transformation are given as, Here, α is a parameter of gauge transformation. The equations of motion of the massless complex scalar field ϕ and the gauge field A µ is given as, Massless scalar electrodynamics preserves relativistic conformal symmetry in d = 4. The complex scalar field ϕ can be expressed in terms of two real scalar fields ϕ 1 and ϕ 2 .
The equations of motion (4.9) in terms of the components are:

Carrollian scalar electrodynamics
As we have just seen, in Carrollian limit, we get two different sectors of electrodynamics depending on the scaling of gauge field. Similarly, the scalar electrodynamics will also exhibit these two different sectors, viz. the Electric and Magnetic sectors.
Electric sector: In this sector, we scale the gauge field in Carrollian scalar electrodynamics in the fashion similar to electric sector of electrodynamics. Along with usual spacetime scaling, the scalings of the gauge field and the scalar field are chosen as (we are considering the most generic situation by scaling the two real scalar fields individually): where p 1 , p 2 are two arbitrary constants which will be determined later. The scalar and gauge field equations in the electric sector are then given by: This of course are just the equations with the scalings put in. The actual EOM will be determined when the appropriate scalings have been performed and the terms with higher powers of as compared to the leading term are set to zero.
The above equations must reduce to free scalar and electrodynamics equations when the gauge and the scalar field are respectively turned off. This leads to constraints on p 1 , p 2 . The constraints that we obtain for electric sector is: (4.14) The allowed values of p 1 , p 2 are shown in Fig[1]. The marked points represent different sectors within electric limit of Carrollian scalar electrodynamics.

Different interacting sub-sectors
Restricting ourselves to integral and half integral values of (p 1 , p 2 ), we find that there are five distinct interacting sub-sectors within the Electric sector of Carrollian scalar electrodynamics. These different sectors for different (p 1 , p 2 ) are: • Case 1 : p 1 = 1, p 2 = 0 : (4.15b) Key: Points in black: free Carrollian theories, Coloured points: non-trivial interacting sectors.
• Case 2: • Case 4: p 1 = 0, p 2 = 1 : We have demarcated the distinct sectors with coloured dots in Fig[1]. The points shown in black reduce to free Carrollian scalar and electric sector of Carrollian electrodynamics. The values of p 1 , p 2 which result in same sectors are connected with each other in the same figure. It is of interest to note here that the non-trivial interacting sub-sectors within the Electric sector of Carrollian scalar electrodynamics all lie along the corners of the allowed region. This is a feature that will be true for all cases we consider after this, where we look at the parameter space of these allowed field theories.

Symmetries of EOM
Now we want to check the symmetries the above equations posses. We will examine an arbitrary sector, say (p 1 = 0, p 2 = 1) to discuss the symmetry. It will require the following inputs from the representation theory, (4.20) In this sector the equations of motion are: Let us now check for the invariance under scale transformation: Hence the scale invariance and special conformal invariance is preserved in d = 4. There is invariance under M m 1 ,m 2 ,m 3 for all the equations.
It can be shown in an identical fashion that in all the other sectors mentioned in our previous analysis, the EOM would also have this infinite dimensional conformal Carrollian invariance. Magnetic sector: Carrollian scalar electrodynamics also has a magnetic sector, named after the particular sector in Carrollian electrodynamics with magnetic field dominating over the electric field. In the magnetic sector, the scaling are given by The equations of motion in this limit are: Studying the equations we find the constraints on p 1 , p 2 as, The allowed values of p 1 , p 2 are easily realised in Fig[2]. The marked points represent different sectors within magnetic limit of Carrollian scalar electrodynamics.

Different interacting sub-sectors
We again find 5 distinct sub-sectors within the Magnetic sector. These different possible sub-sectors are: • Case 1: p 1 = 1, p 2 = −1: • Case 3: • Case 4: p 1 = −1, p 2 = 1: • Case 5: p 1 = − 1 2 , p 2 = 3 2 : The black points in Fig[2] reduce to free Carrollian scalar and magnetic limit of Carrollian electrodynamics. The non trivial sectors having interaction terms are marked in different colours in Fig[2]. The values of p 1 , p 2 which result in same sectors are connected with coloured lines in the same figure.

Symmetries of EOM
Before discussing the symmetries, we require the values of the constants {∆, ∆ , r, a, b} from the representation theory.
We are choosing a sector randomly, p 1 = 1 2 , p 2 = − 1 2 . The equations in this sector are: We will now check for invariance of the above equations. Under scale transformations, we have Similarly, under K i , the transformation is given as The right hand side of all of the above expression vanishes as {∆ = 1, ∆ = 1}. The equations of motion also have infinite M m 1 ,m 2 ,m 3 symmetry.
We can perform the same symmetry analysis for all the possible subsectors and we find that the emergent infinite dimensional symmetry exists in all of these sectors.

Adding Fermions
The relativistic theory of U (1) gauge field A µ coupled to fermionic fields Ψ exhibit conformal symmetry in d = 4. The Lagrangian density of this theory is given as where D µ is the covariant derivative and F µν is the electromagnetic field strength: The relativistic equations of motion are: We will decompose the Dirac fermion in two component spinors as in Sec [3.2]. The relativistic equations of motion get recasted as, Next, we are interested in studying the Carrollian limit of the same theory.

Carrollian electrodynamics with fermions
We have seen that the gauge fields can be scaled in two ways, electrically and magnetically.
Here, we will also scale the spinors accordingly. For easier understanding, we will be using the conventional names of the different limits in a way similar to Carrollian electrodynamics.
Electric sector: In this sector, the scaling of gauge fields and fermions are given as: Here α and β are two arbitrary constants which will be determined later. In this sector, the equations of motion become: The above equations must reduce to the free Carrollian equations of the electric sector of electrodynamics and fermions separately. These constraints put some restrictions on the arbitrary constants α, β. This can be seen in Table[2].
Combining all the constraints we finally arrive at the following relations:

Different sub-sectors in the theory
We will now write down different sectors for the Electric sector of Carrollian electrodynamics coupled to fermions: • Case 1: α = 0, β = 1: • Case 2: α > 0, β = α + 1: We see that in Case 2, the equations (4.46) reduces to Carrollian equations of electrodynamics and free Dirac fermions.

Symmetries of EOM
The next task is to find the invariance of these equations under conformal Carrollian symmetries. For that, we will only consider Case 1, since these equations contain interaction terms. For this case, the values of the constants are mentioned here: The invariance under scale transformation is given by Similarly, the invariance under M m 1 ,m 2 ,m 3 can be worked out in a similar way.
We see that the Electric sector comes out to be invariant under finite and infinite conformal Carrollian symmetries in d = 4. Magnetic sector: The scaling of the fields in this sector are given by

Free equation Constraints
In this limit, the equations of motion becomes Comparing the above equations with the free equations of magnetic sector of Carrollian electrodynamics and free fermions we obtain restrictions on the arbitrary constants α, β mentioned in Table 3. Analysing all the constraints we get the following relations: Hence the possible values of α, β are:

Symmetries of EOM
The values of the constants under Carrollian symmetries are given as (4.56) We will consider Case 1 as the representative sector to examine the symmetries. The same analysis can be repeated for Case 2, with similar results. The transformations of the equations of motion under scale transformations are given by Putting in the values of {∆, ∆ }, we find that the equations are invariant. Similarly, under K l , we have Again, putting in the values of the constants, the EOM turn out to be invariant. Under the infinite dimensional M m 1 ,m 2 ,m 3 , the invariance of the EOM are straight-forward.
So, we can conclude that the Magnetic sector of Carrollian U (1) gauge field coupled to fermion respects the conformal Carrollian symmetry in d = 4.

Lessons from this section
For the reader who does not wish to be buried in the details of the rather tedious analysis we have done so far in this section, we will provide a brief recapitulation of the main results. In this section, we have carried out a detailed analysis of Carrollian Electrodynamics theory with scalars and fermions. In case of Carrollian Electrodynamics, there were only two sectors, viz. the electric and magnetic sectors. However, the inclusion of matter into the theory led to many diverse subsectors inside each of the two above mentioned sectors. Massless scalars led to the splitting of each sector into 5 different interacting sub-sectors. The inclusion of fermions made the splitting a bit more restrictive and there were two subsectors in each of the Electric and Magnetic sectors.
Next, we checked for the symmetries of the theories under conformal Carrollian generators. We found that the field theories possess finite as well as infinite conformal Carrollian symmetries in d = 4. The infinite extension of symmetries makes Carrollian Electrodynamics with matter prominently distinct from its relativistic counterpart.

Carrollian Yang-Mills theory and massless fermions
So far we have constructed Carrollian field theories with Abelian gauge symmetry. We have found that these theories respect finite and infinite conformal Carrollian symmetry in d = 4. At this point, it is natural to ask if we can extend our analysis to field theories with non-Abelian gauge symmetry as well. The ultra-relativistic limit of SU (N ) Yang-Mills has been described in detail in [19]. In this section, we would construct Carrollian SU (N ) Yang-Mills theory coupled to massless fermionic matter in d = 4.

Carrollian Yang-Mills Theory: A brief review
This section is going to be a quick review of pure Carrollian Yang-Mills (YM) theory presented in [19]. The field content of Carrollian YM theory is the gauge potential Here T a (a = 1, . . . , D) are the Lie algebra generators and the dimension of the gauge group is D. In [19], it was shown that there are a large number of new sectors arising in the Carrollian limit, as opposed to the two sectors in Carrollian Electrodynamics. The reason lies in the projection of individual gauge fields into either individual electric or magnetic sectors.

Different Sectors
Following [19], we divide the theory into D + 1 sectors, each characterised by a vector: Ξ a (p) denotes its a'th component. We choose a particular sector (say, p 0 -th) denoted by the vector Ξ (p 0 ) . In p 0 -th sector Ξ (p 0 ) , the gauge fields scale as If Ξ a (p 0 ) = 1, the gauge fields transform electrically and we denote them by the index α, β: On the other hand, for Ξ a (p 0 ) = 0, the gauge fields transform magnetically and we denote them by the index I, J: It is seen from (5.2) that I, J run from 1, 2 . . . , D − p 0 and α, β are in the range D − p 0 + 1, . . . , D. The Carrollian Yang-Mills equations are: The Carrollian scalar equation is: The vector one is: Scalar equation: Vector equation: It is worth noting here that the equation (5.8) does not have a kinetic term. Hence, these sectors would not be considered when we are looking at sensible Carrollian theories. This argument turns out to be a little too hasty. Notice also that the interaction piece vanishes when there is only one Magnetic leg (I = J = 1). Hence, in [19], the next higher order term in the series of was considered and this led to the following equation: So, this is a mixed leg that we would consider as a valid sector.
To summarise, the sectors of pure Carrollian Yang-Mills theory that we interested in are: • Purely Magnetic sector: The equations of relevance here are {(5.6), (5.7)}; • Purely Electric sector: The purely electric sector survives. This is the one where we consider {(5.9) (with I = 0), (5.10)}. There are no magnetic legs, and hence we have to put A I = 0 in these equations. This leads to: • Mixed sector (with J = 1): Only one mixed sector survives. This corresponds to p 0 = D − 1 and equations {(5.9) (with I = 1), (5.10) (with I, J, K = 1)}.
It was also shown in [19] that these the sectors in Carrollian Yang-Mills theory respect finite and infinite conformal Carrollian symmetries in d = 4.

Yang-Mills with Fermions
The Lagrangian density of the relativistic SU (N ) Yang-Mills with fermions is given by where D µ ≡ ∂ µ − igT a A a µ is the non abelian gauge covariant derivative and F a µν = ∂ µ A a ν − ∂ ν A a µ + gf abc A b µ A c ν is the non abelian field strength tensor. The label a is the colour index, m is an internal symmetry index and f abc are the structure constants of the underlying gauge group with generators following the algebra [T a , T b ] = if abc T c . The equations of motion for this theory are The relativistic equations of motion under the decomposition of Dirac fermion (3.15) are The relativistic Yang-Mills coupled to fermions have conformal symmetry in d = 4.

Carrollian SU (N ) Yang-Mills with Fermions
The simplest example of a non-Abelian group is, of course, SU (2). We have discussed the details of the Carrollian SU (2) Yang-Mills with fermionic matter extensively in Appendix A, to which the interested reader is pointed. Here we proceed to present the general analysis for Carrollian SU (N ) Yang-Mills with matter. The discussion can also be generalised immediately to any other compact gauge group. Let us now move to the scaling of fields in SU (N ) case. The scaling of gauge fields is very similar to the scaling in [19] and is briefly explained in the previous section. Next step will be to take scaling on the fermions (φ m , χ m ) where m = 1, . . . , N . The scaling is taken as, We will look at the equations of motion of different sectors of SU (N ) Carrollian Yang-Mills coupled to fermions: The equations of motion in this case are given as For simplicity, we will denote the LHS of first equation by A and second by B.
• D − p 0 + 1 ≤ a ≤ D: The equations of motion in this case are given as We will denote the LHS of the third and fourth equation by C and E. When there is only one magnetic leg I = J = 1, the first term in (5.19b) goes to zero by the antisymmetric property of the structure constant f α IJ , (f α 11 = 0). For that reason we consider the next order term in (that is 0 ), and the equation becomes The Dirac equations are same in both the cases and they are given as Hence, in this section we are left with the equations A, B, C, E, F, G, H depending on different sectors of Carrollian SU (N ) Yang-Mills coupled to fermions. Please note that we are yet to specify a m 's and b m 's, which we will be doing in the next section.

Detailed discussion of different sectors
At this juncture, it is important to remind ourselves of the discussion in the pure Carrollian YM case. There we had discarded sectors which had no kinetic terms. This obviously would also continue to hold when we put in matter fields. Clearly, a sector where the equation E (5.19b) holds, is not a "nice" sector. In the analysis that follows, we will continue to investigate all sectors, keeping in mind, that if we were to focus on "nice" sectors of Carrollian YM coupled to matter, all the sectors which have a non-trivial realisation of (5.19b) would be discarded. In case of SU (2) Carrollian Yang-Mills with fermions there are three gauge fields and 4 different sectors: EEE, EEM, EMM, MMM. Similarly for SU (N ) case, we can have D + 1 distinct sectors. We would organise D + 1 sectors depending on the numbers of magnetic legs.
Number of Magnetic legs J ≥ 1: We start with the case in which one and more gauge field transform magnetically. For J = 1, the equation E does not hold. Hence instead of E, we get Magnetic leg J ≥ 1

From equation Constraints
Equation A of Gauge field A I (i) a m + a n + 1 ≥ 0 (xiv) a n − b m + 1 ≥ 0 Table 4. Constraints on a m , b m in J ≥ 1 sector.
The rest of the equations A, B, C, G, H remains same with J = 1. Now we want to see the constraints on a m and b m . In Table[ 4], we find the constraints by comparing with the free equations. Next, we analyze the range of a m and b m . From equation (xii) and (xiv) we can see, (xiv) : a n − b m + 1 ≥ 0 =⇒ a n − a m ≥ 0.  Table 5. Values of a m , b m in J ≥ 1 sector.
We are yet to specify the range in which a m belong. They are given by For J = 1 (vii) : a m + a n ≥ 0 =⇒ a m ≥ 0 (5.26) Combining equations (5.23)-(5.26), we get the Table [5]. Writing the equations of motion in each sector separately, we get For the case J > 1, the equations of motion are written as, • Case 2: (a m > − 1 2 , b m = a m + 1): Both of the above mentioned sectors contain no kinetic term in (5.27c),(5.28c). But (5.27b),(5.27c),(5.28b) and (5.28c) drop off completely when α = 0, or in other words when there is no electric leg. It corresponds to the purely magnetic limit of SU (N ) Yang-Mills. Hence, for the purely magnetic limit we get two individual sectors which differ by a fermionic interaction term gφ † m T I mn φ n within themselves (see (5.27a) and (5.28a)).
For the case J = 1, the equations of motion are written as, • Case 1: (a m = 0, b m = 1) : • Case 2: (a m > 0, b m = a m + 1): From the above analysis we see that there are (2+2) different sectors in the purely magnetic (α = 0) and only one magnetic leg (J = 1) limit respectively.
Number of Magnetic legs J = 0: This sector defines the case in which all the gauge fields transform electrically. The equations of motion are given as Similar to the previous case, we will be writing down Table [6] to have a clear understanding of the constraints on a m and b m . The relevant equations that we need to look at are (ii), (v) and (vii) in Table [6].
(ii) : a m + a n ≥ 0, (v) : a n − a m + 1 ≥ 0, Unlike J ≥ 1 sector, we get no more constraints to specify a m . Hence there will be many sectors describing J = 0 case. At this point, we want to count how many sectors actually arise for the purely electric case of Carrollian SU (N ) Yang-Mills coupled to fermions considering the constraints (5.32).
Purely Electric sector (J = 0)   Table[ 7]. Any other combination of these a m 's, which obeys the constraints (5.32), gives back no new sector.

Gauge Transformation for Carrollian SU (N ) Yang-Mills with Fermions
Let us consider the transformations of gauge fields and fermions under SU (N ) gauge group: where Θ is the gauge parameter. We will now see gauge transformations in ultra-relativistic limit of SU (N ) Yang-Mills theory. We scale the gauge transformation parameter for the electric and magnetic legs as We would be plugging (5.34) into (5.33) along with the usual ultra-relativistic scaling of spacetime, gauge fields and fermions. Our aim is to find the value of p, q which keep (5.33) finite in the limit → 0. They are given by Thus we arrive at the following gauge transformation: Before finding the invariance, we write the values of the constants in representation theory. We should keep in mind that, here, all the gauge fields scale electrically.
The invariance under scale transformation are given as Under special conformal transformation, Hence, all the above mentioned equations in J = 0 sector, have infinite conformal Carrollian symmetry in d = 4. Next, we will examine the sections where J = 1 and J > 1. We are choosing {a m = 0, b m = 1} for all m. The relevant equations in these two representative sectors are, The values of constants in representation theory are given by, The invariance under scale transformation is Under special conformal transformation, we have The above analysis shows that both the J ≥ 1 sectors have finite as well as infinite conformal Carrollian symmetries in d = 4.

Brief summary of the current section
Let us summarise the results that we have obtained in this section. We started out with a brief review on SU (N ) Carrollian Yang-Mills theory. We explained the scaling of gauge fields (A a t , A a i ) in detail and wrote down the equations of motion. These equations were invariant under conformal Carrollian generators including the infinite supertranslations in d = 4. For further details about the invariance of these equations, the reader is referred to [19]. One of the upshots of the analysis was that the number of "nice" sectors with EOM having kinetic pieces, boiled down to just 3, viz. the purely electric and magnetic sectors, and one mixed sector with one magnetic leg.
Then, we moved on to SU (N ) Carrollian Yang-Mills coupled with fermions. The scaling of gauge fields gave rise to (D + 1) sectors, of which only 3 were "nice", as just mentioned. For the sake of completeness, we did not discard the other sectors outright. We arranged the (D + 1) sectors depending on the number of magnetic legs, denoted by J. For J = 1 and J > 1, the equations (A, B, C, G, H) (see eq.(5.18a) -(5.21b)) were held perfectly. For J > 1, equation E (5.19b) was true. This was the problematic equation in the pure YM case, which contained only interaction pieces. So in principle, all these sectors with J > 1, except the purely magnetic sector, should be discarded on grounds of having a non-trivial realisation of this equation. For J = 1, we observed that the equation E did not hold and we had to consider the next order in , captured by equation F (5.20). Each of the J = 1 and J > 1 sectors had two different subsectors within them (see eq. (5.27) -(5.30)).
We also considered J = 0 case, the purely electrical sector. The absence of the constraint a n = a m for all m, n's into the scaling of fermions (that appeared in the previous cases, see eq. (5.24)) resulted in a more diverse class of subsectors (see Table[ 7]).
We looked at the gauge transformations and the invariance of equations of motion in each subsector. We took a particular representation for each of the subsectors in J = 0, ≥ 1 cases. We investigated the symmetries of the Carrollian YM theory coupled to fermions at the level of equations of motion. Again, like in the previous sections, we found the emergence of finite as well as infinite conformal Carrollian symmetries d = 4 in all of these subsectors. This emergence of infinite enhancement of conformal Carrollian symmetries in d = 4, thus, seems to be a generic feature of all field theories that arise as a ultra-relativistic limit of classically conformally invariant theories in d = 4.

SU (3)
SU (4) SU (5) + fermion + fermion + fermion + fermion In this paper, we have systematically constructed conformal Carrollian field theories in d = 4 by following the ultra-relativistic limit of relativistic CFTs in d = 4. We have started with free scalars, fermions and reviewed previous constructions of Carrollian gauge theories. Then we have constructed arbitrary interacting theories by adding matter to general Carrollian gauge theories. The intricacies of the limit resulted in a multitude of different sectors of these gauge field theories. We showed that in all these sectors the equations of motion possess an infinite dimensional symmetry. We again stress here that this enhancement of symmetries in conformal field theories in the ultra-relativistic limit seems to be a generic feature of this limit. We also saw that algebraically the process of infinite extension for d > 4 was very similar to d = 4. Taking this analogy further, we propose that if one starts with any relativistic conformal field theory in any dimension and constructs its Carrollian version by taking c → 0 limit, the equations of motion of the ultra-relativistic theory so constructed would be invariant under the infinite conformal Carrollian algebra.
Our analysis in this paper, as stated earlier, has been a direct follow-up to earlier work [19]. The underlying mathematical formulation of the investigation of symmetries from the point of view of equations of motion actually is very similar in nature to the investigations of non-relativistic systems and the emergence of infinite Galilean conformal symmetry in theories derived as a limit of relativistic conformal field theories [4][5][6]. Given the similarity of these analyses, it is important to compare and contrast the two different singular limits of gauge theories. In [6], the non-relativistic counterpart of the field theories studied in the current paper were constructed, namely Galilean scalar, fermions, Yukawa theory and interacting theories like electrodynamics and SU (2) Yang-Mills coupled to matter. This analysis has been extended for the case of SU (N ) Yang-Mills coupled to matter in appendix B, which also contains a brief introduction to aspects of the non-relativistic symmetries for the uninitiated reader. The important question of the connection these two opposite limits (Galilean and Carrollian) in terms of their similarities and differences is discussed in detail in Appendix C.
The infinite enhancement of symmetries in the ultra-relativistic limit, that has been the central theme of our present paper, seems intimately related to the geometry of the spacetime that these field theories live on. If we start with relativistic CFTs that live in flat spacetimes, the Carrollian limit leads to a degeneration of the metric structure of the Minkowskian manifold: The relativistic conformal structure is more restrictive as it involves space and time on the same footing. In the ultra-relativistic case, the (pseudo) Riemannian structure (IR 1,d ) turns into a fibre bundle structure with a IR d base and temporal one-dimensional IR t fibres. Now there is more "wiggle-room" and the conformal structure gets split into that of the base and the fibre. This degeneration of the spacetime metric into a purely spatial part (for the base) and a purely temporal part (for the fibre) seems to be at the heart of these infinite dimensional enhancements in the ultra-relativistic limit. This also should be true in the non-relativistic limit, where similar infinite enhancements are noticed. We wish to investigate these geometrical aspects involving Carrollian and its dual Newton-Cartan structures more carefully in the near future.
As we stated in the introduction, the ultimate aim of this programme is to make connections with holography in asymptotically flat spacetimes. As a specific case, we would like to understand what happens to the dual field theory when we take a flat limit on the best known example of the AdS/CFT duality, viz. the duality between type IIB superstring theory on AdS 5 ×S 5 and N = 4 Supersymmetric Yang-Mills (SYM) theory in d = 4. In the bulk, the radii of the AdS 5 and S 5 are equal and hence taking a limit on one induces a limit on the other, thereby giving us type IIB superstring theory on a 10d flat spacetime. On the boundary, the limit would lead to a Carrollian version of N = 4 SYM. In future work, we would attempt an understanding of SYM theories with various supersymmetry building up to this example. The analysis done in this paper would provide the basic ingredients for this investigation. It seems very unclear at present what these various sectors arising in the Carrollian limit would correspond to in the holographic setting. The infinite dimensional emergent symmetry, on the other hand, has tantalising hints of a new integrable sector in the parent gauge theories and relatedly holography.
One of the stumbling blocks of our analysis in this paper, and our programme as a whole, has been the lack of an action principal. For the case of Carrollian electrodynamics, this has been partially addressed and one has an action for the Electric sector [22], and the same symmetries have been understood by a rigorous canonical analysis 4 . In ongoing work [82], it has been possible to show that one can have a similar action principal for scalars coupled to Carrollian electrodynamics for a particular sector. Generalisation of this to the other possible sectors and then to the case of Yang-Mills theories and SYM is the obvious goal of this line of work. With an action, quantities of interest like correlation functions can be derived by functional methods and this would also pave the way to a possible quantisation of the Carrollian theories and investigations into anomalies to understand whether the infinite dimensional symmetries survive in the quantum regime. This work is also partially supported by the following grants: DST-INSPIRE faculty award, SERB Early Career Research Award (ECR/2017/000873), DST-Max Planck mobility award, DST-BMWF India-Austria bilateral grant, Royal Society International Exchange grant, SERB extra-mural grant (EMR/2016/008037).

Appendices A Carrollian SU (2) Yang-Mills with Fermions
In the main body of the paper, we have dealt with a general SU (N ) Carrollian Yang-Mills theory coupled to fermionic matter. It may be of interested for the reader to understand the elementary details and hence in this appendix, we present the complete details of the SU (2) Carrollian theory with fermions.
In SU (2) Yang-Mills with fermions, we have three colour indices for the gauge fields and a doublet of fermions in the fundamental representation. The relativistic equations of motion are given by (5.16). In the Carrollian limit, we scale the fermions as, In [19], the detailed analysis of Carrollian SU (2) Yang-Mills is carried out. There are four different sectors in case of SU (2) Yang Mills: EEE (purely Electric), MMM (purely Magnetic), EEM and EMM (2 mixed sectors). Following the same trail, in case of SU (2) Yang-Mills coupled to fermions, we scale the gauge fields in the similar way. Let us discuss the different sectors in details below.

A.1 EEE sector
In this sector, all gauge fields are scaled electrically. The scaling is given as where m = 1, 2. The fermions are scaled as in (A.1). The equations of motion in this sector are given as The above equations must reduce to free equations of a particular field in absence of rest of the fields. Also, if we replace T a by an identity matrix and set the structure constant to be zero, we should get back the equations of U (1) Carrollian electrodynamics coupled to fermions. For example, suppose we are looking at the limit EMM of the SU (2) gauge fields

EEE sector
Free equation coupled to fermions. In this sector A 1 scales electrically and A 2 , A 3 scales magnetically. So we demand that if we replace T a by an identity matrix, the A 1 equation must reduce to (4.45) or (4.46) while the rest of the gauge fields must follow (4.54) or (4.55). The same goes for the fermions as well. If we turn off the gauge fields in the equations obtained in EMM case, we must get back the free fermion equations. If we keep the gauge field, replace T a by an identity matrix, we should get the fermion equation from either (4.45),(4.46) or (4.54), (4.55)(depending on how the particular gauge field coupled to the fermion scales). The above conditions impose constraints on the constants {α m , β m }. Please see Table [8].
There is no equality in some of the constraints in the Table [8] to reduce it to Carrollian electrodynamics with fermions case. Constraint (ii),(vi) and (vii) gives The rest of the constraints give no new information. In Mathematica, we can do a Regionplot [3] to plot (A.4) and find the possible values of α 1 , α 2 . Each intersecting point in Regionplot [3] corresponds to a sector in the EEE limit. Some of the possible values are: We will now discuss different sectors within EEE limit. All the intersecting points in Plot [3] do not give unique sectors. We are enlisting below only the different sectors. . Allowed values of (α 1 , α 2 ) in EEE sector.

Gauge Transformation for EEE sector
Let us consider relativistic gauge transformations of gauge and fermions fields as We will now see gauge transformations in this sector. We scale the gauge transformation parameter for EEE sector by Θ a → p Θ a . (A.12) We would be plugging (A.12) into (A.11) along with the usual ultra-relativistic scaling of spacetime, gauge fields and fermions. Our aim is to find the value of p which keeps (A.11) finite in the limit → 0. It is given by Thus we arrive at the following gauge transformation for EEE sector of Carrollian SU (2) Yang Mills coupled to fermions: To see the invariance under these transformation, we have to take one set of values of the coefficients (α m , β m ) along with → 0 limit and plug them into (A.14).

Checking symmetries of EEE sector
Next, we will find the invariance of the equations of motion. For that we will choose a representative sector, that is α 1 = 0, α 2 = 0, β 1 = 1, β 2 = 1. The equations are Before checking invariance under Carrollian symmetry, we remind the readers, the values of constants in representation thoery: EEM sector Table 9. Constraints on α n , β n in EEM sector.

From equation Constraints
• For spinors φ m , χ m : Following the same reasoning mentioned above we can find the constraints on (α m , β m ) in Table [9].
From the allowed values of (α 1 , α 2 ) we can construct only two different sectors.

Gauge Transformation for EEM sector
For the EEM sector, we scale the gauge transformation parameter as, We want to find p, q in the same way mentioned before. Now we scale Θ a 's differently for A 1,2 and A 3 . The values comes out to be The gauge transformation in this sector is given by The invariance of each subsector under these transformations can be seen when we take one set of values of the coefficients (α m , β m ) along with → 0 limit and plug them into (A.31).

Gauge Transformation for EMM sector
For the EMM sector, we scale the gauge transformation parameter as, We scale Θ a 's differently for A 1 and A 2,3 . The values of p, q comes out to be The gauge transformations in this sector is given by The invariance under these transformations can be seen when we take one set of values of the coefficients (α m , β m ) along with → 0 limit and plug them into (A.45).

Gauge Transformation for MMM sector
For the MMM sector we scale the gauge transformation parameter as, We want to find the value of q that keeps the gauge transformation (A.11) of fermions and gauge fields finite along with the ultra-relativistic spacetime contraction. The value of q comes out to be q = 2, Θ a → 2 Θ a . (A.55) We also scale the gauge fields magnetically and the fermions in the usual way. Thus we arrive at the gauge transformation of MMM limit of Carrollian SU (2) Yang Mills coupled to fermions: When we take one set of values of the coefficients (α m , β m ) along with → 0 limit and plug them into (A.56), the equations remain invariant in that subsector.
Checking symmetries of MMM sector The representative limit for this sector is chosen to be α 1 = 0, α 2 = 0, β 1 = 1, β 2 = 1. The equations are The values of the constants are The invariance under scale transformation is Under SCT, we have The equations are also invariant under M m 1 ,m 2 ,m 3 .

B Galilean SU (N ) Yang-Mills Theory with matter
In this section we extend the analysis of [6], to our non-relativistic counterpart of SU (N ) Yang-Mills coupled to matter. It is a further generalisation of Galilean SU (2) Yang-Mills with matter described in [6].

Galilean Conformal Algebra
We start with a brief review on the Galilean Conformal Algebra (GCA). GCA is obtained from conformal algebra by performing an Inonu-Wigner contraction. In the process of going to the Galilean framework, the underlying Lorentz symmetry breaks down. The spacetime scale differently, as The generators of GCA can be found by imposing the limit (B.1) on relativistic conformal algebra. The generators in this limit emerge as where L (−1,0,1) = H, D, K (H, D and K are the Galilean Hamiltonian, dilatation and temporal SCT) and M (−1,0,1) i = P i , B i , K i (P i , B i and K i represent spatial momentum, Galilean boost and spatial SCT). The full GCA is given as The interesting point to be noted here is that the above algebra closes, even if the index n runs over all integers. Hence, it is concluded that GCA has an infinite extension in all spacetime dimensions, which is very unlike relativistic conformal algebras for d > 2.

Galilean Yang-Mills
In [6], a detailed analysis of Galilean conformal field theories was performed including free as well as interacting theories. The last example presented in this paper was Galilean SU (2) Yang-Mills theory coupled to fermions. The various sectors of O(1500) emerged as sublimits of EEE, EEM, EMM and MMM. However, using some consistency requirements the possible sectors were brought down to 19 (including both interacting and free sectors). This significant reduction in the number of sectors hinted towards something deeper coming into the scenario, which is unfortunately yet to be explained. Another interesting observation was that, the equations of motion in all the sub-limits of SU (2) YM with fermions turned out to be invariant under GCA. The purpose of this appendix is to generalise the existing calculations of SU (2) case [6] to non-relativistic limit of SU (N ) Yang-Mills theory coupled to fermions. We would look into the structure of equations of motion and try to find the constraints on parameters a m , b m appearing in the theory. Ultimately, we would check the invariance of a representative sector under GCA.

C Comparison between Galilean and Carrollian CFTs
The reader familiar with earlier work [4][5][6] will find many similarities in the methods and indeed the final results of Carrollian and Galilean conformal field theories. Here we attempt to make a comparison between these field theories, which can be thought of as two distinct sub-sectors of a parent relativistic theory.
The underlying algebras, the CCA and the GCA, are obviously different as they appear from different inequivalent contractions of the relativistic conformal algebra. For the finite algebras, arising solely out of contractions, the number of generators are the same. As we have seen, the infinite extensions in both cases have different flavours and the algebraic structure of the infinite dimensional algebras are quite different. The GCA has the same structure for all dimensions, which the CCA differs significantly depending on dimensions (and whether or not one can give the super-rotation part an infinite lift).
Moving on to the explicit examples of the field theories, we find quite a few important distinctions. We will begin the discussion with our first example, the scalar field theory. In [6], we found that in Galilean limit (B.1) of ∂ µ ∂ µ ϕ = 0, we get ∂ i ∂ i ϕ = 0. (C.1) The above expression (C.1) contains spatial derivative only. If we take Carrollian limit (3.2) on scalar equation, the result comes out to be ∂ t ∂ t ϕ = 0, (C.2) which only involve time derivatives. This is a generic feature and obviously has to do with the underlying contraction where in the non-relativistic case, space is scaled down and as a result, spatial derivatives scaled up as opposed to the Carrollian case where time and time-derivatives are scaled down and up respectively. For our second example, we will look at massless Dirac field theory. In Galilean case, the equations are given as where φ, χ are two component spinors in Pauli-Dirac representation. In Carrollian case, the equations are given by (3.19). In both cases, the scaling of the spinors are taken to be φ → φ, χ → χ. Again, as expected, in the NR case, we get an equation which contain only spatial derivatives whereas in the UR case we have the time derivative acting on the spinor.
Let us now move on to Electrodynamics. The relativistic EOM are: To take the non-relativistic limit on (C.4), we used the non-relativistic scaling of the coordinates and scaled the gauge field components as Electric sector: A t → A t , A i → A i (C.5a) Magnetic sector: A t → A t , A i → A i (C.5b) The reason that we get these two sectors is due to that fact that the contravariant and covariant vectors behave like two independent quantities (because the metric becomes degenerate). This can also be seen clearly as follows where we have usedx i → x i ,t → t in the intermediate steps. We see that (C.6a) can be associated with electric sector, whereas (C.6b) denotes the magnetic sector. In case of ultrarelativistic limit, the gauge fields are scaled in a similar fashion (C.5). But contravariant vectors gives rise to magnetic sector (see (4.2)) and the covariant ones give the electric sector. This is because of opposite scaling of spacetime coordinates (x i → x i ,t → t) in contrast to Galilean case. The equations of motion in Galilean case and Carrollian case are also different. Nonrelativistic limit taken on the relativistic equation yields Electric sector: Magnetic sector: For Carrollian limit, we have Electric sector: Magnetic sector: We now want to focus on the differences that occur when we check for the symmetries of the EOM in both cases. We have worked in the scale-spin representation. The scalespin representations of GCA and CCA are determined by the set {∆, j, a, b}. The constants a, b are boost labels on the fields, exclusive to each spin. For Galilean Electrodynamics (GED), the values of the constants are given by (a e , b e ) = (1, 0) for electric limit and (a m , b m ) = (0, 1) for magnetic limit 5 . For Carrollian Electrodynamics case, the value of the constants are given by (a e , b e ) = (0, 1) for electric limit whereas for magnetic limit, the values are (a m , b m ) = (1, 0). We see that the boost labels are interchanged between the sectors when we compare both limits. For example, the values of constants in electric sector of GED are same as that in magnetic sector of Carrollian ED.
Moving on to gauge fields coupled to matter, we also see some differences. In GED coupled to fermions, the magnetic limit is devoid of any interactions between the fields In Carrollian ED coupled to fermions, only the sectors with (α > 0) of electric limit is devoid of the interaction terms and can be seen explicitly from the equations Another point of difference is that while taking generalised scaling of fermions (φ, χ) as φ → α φ, χ → β χ, (C.11) the minimum value of α for electric limit turns out to be α = −1 and for magnetic limit α = − 1 2 in Galilean case, whereas in Carrollian case, the values are given as α = 0 (for electric limit) and α = − 1 2 (for magnetic limit). These minimum values of α show that any value below this, fails to reproduce the free equations (kinetic terms).
Lastly, we consider the example of SU (N ) YM theory. For the case, D−p 0 +1 ≤ a ≤ D in SU (N ) Carrollian YM (section 5.2), we have two separate scalar equations because of the magnetic index J. For J = 1, the equation is denoted by F (5.20) and for J > 1, we obtain E (5.19b). Correspondingly, we get two different constraints on a m for SU (N ) Carrollian YM coupled to fermions. They are given as a m ≥ 0 for J = 1 and a m ≥ − 1 2 for J > 1. (C.12) The sectors where the E hold, are not "nice" sectors, as these equations have only interaction pieces and no kinetic terms. This means that we need to discard all those sectors where E has a non-trivial realisation. In fact, the only surviving sector with J > 1 is the purely magnetic one as E drops off there. In SU (N ) Gallilean YM, we do not have any such divisions because we do not get any equation similar to E. We only have single scalar equation for D − p 0 + 1 ≤ a ≤ D case and it is given as ∂ i ∂ i A α t = 0. On the one hand, this is a simplification in terms of the underlying equations. On the other hand, the absence of the equivalent of equation E means that one has a huge number of possible sectors to deal with in the Galilean YM theories.
For the convenience of the reader, we have listed the points we discussed above in an extensive table in the next page.

4.
We will write some of the commutators to show the difference between two algebras. They are written down as Similarly, for CCA we have

5.
We can have representations labelled by the scale and boost in all dimensions.
Here the scale-boost representation is only possible in d = 2 case. For d > 2 case, there are no non-trivial representation that are labelled by weights under scale and boost.
In case of scalars, the equation of motion is given as ∂ i ∂ i ϕ = 0. Here, we have no time derivatives.
For scalars, the equation is given as ∂ t ∂ t ϕ = 0. We see the absence of spatial derivatives.

7.
In fermion case, one of the equations just contains the spatial derivatives and is given by iσ i ∂ i φ = 0 Whereas, in this case, one of the equations contains only the time derivative and is given as i∂ t φ = 0.

8.
In GED case, the contravariant vectors are associated with electric sector, whereas the covariant vectors give rise to the magnetic sector.
In Carrollian ED case, the contravariant vectors are associated with the magnetic sector, whereas the covariant vectors denotes the electric sector.

9.
The constants that defines the scale-spin representation are given by a, b. For GED, the values of the constants are given by a e = 1, b e = 0 for electric limit and a m = 0, b m = 1 for magnetic limit.
In Carrollian ED case, the values of the constants are given by a e = 0, b e = 1 for electric limit whereas for magnetic limit, the values are a m = 1, b m = 0.
10. In GED coupled with fermions, the magnetic limit is devoid of any interactions between the fields.
In Carrollian ED coupled with fermions, the electric limit (only the sectors with α > 0) is devoid of the interaction terms.
11. In SU (2) GYM, the self interactions among gauge fields are present in EEE, EMM sector.
In SU (2) Carrollian YM, the self interactions among gauge fields are present in EEM, MMM sector.
12. In SU (N ) GYM, we do not have any such divisions because we do not get any equation similar to E. We only have a single scalar equation for the D − p 0 + 1 ≤ a ≤ D case.
For the case D − p 0 + 1 ≤ a ≤ D in SU (N ) Carrollian YM, we have two separate scalar equation because of the magnetic index J. For J = 1, the equation is denoted by F and for J > 1, we have E. Due to such division, we get two different constraints on a m for SU (N ) Carrollian YM coupled with fermions. For J = 1, we have a m ≥ 0 and a m ≥ − 1 2 for J > 1.
13. Number of sectors in GYM is large as there are no sectors we discard for nonexistence of a kinetic term The number of sectors in Carrollian YM is more limited due to the existence of sectors where there are EOM like E. This class of EOM only contain interaction terms. So all sectors with non-trivial realisations of this class of EOM are discarded.