{'fraction_non_alphanumeric': 0.04318159664684667, 'fraction_numerical': 0.011989472658153815, 'mean_word_length': 4.373254189944134, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 23, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Non Local Global Symmetries of a Free Scalar Field in a Bounded Domain

BMS symmetries have been attracting a great deal of interest in recent years. Originally discovered as being the symmetries of asymptotically flat spacetime geometries at null infinity in General Relativity, BMS symmetries have also been shown to exist for free field theories over Minkowski spacetime. In wanting to better understand their status and the underlying reasons for their existence, this work proposes a general rationale towards identifying all possible global symmetries of a free field theory over Minkowski spacetime, by allowing the corresponding conserved generators not to be necessarily spatially local in phase space since fields and their conjugate momenta are intrinsically spatially non local physical entities. As a preliminary towards a separate study of the role of asymptotic states for BMS symmetries in an unbounded Minkowski spacetime, the present discussion focuses first onto a 2+1 dimensional free scalar field theory in a bounded spatial domain with the topology of a disk and an arbitrary radial Robin boundary condition. The complete set of global symmetries of that system, most of which are dynamical symmetries but include as well those generated by the local total energy and angular-momentum of the field, is thereby identified.


Introduction
Consider a free field, whether classical or quantum, over a background spacetime of specific geometry, say that of a Minkowski spacetime. Besides the symmetries of that geometry -thus those of the Poincaré group in the instance of a Minkowski spacetime -what are all the possible global symmetries of such a field dynamics?
Of course these symmetries include the global symmetries of the background spacetime geometry. However there are far more symmetries possible still, certainly in the case of a free field. As is well known -a result which came as a surprise when it was first discovered -in Einstein's General Relativity an asymptotically flat spacetime possesses asymptotic global symmetries at null infinity that include the Poincaré group but are in fact much larger with an infinite (countable) set of conserved quantities, namely the generators of the so-called BMS symmetries of supertranslations and super-rotations [1,2,3,4,5,6]. More recently it has been established that free fields in a Minkowski spacetime share the same BMS symmetries as well [7,8,9]. In particular in Refs. [8,9] it was shown that these BMS symmetries are related to conserved charges and generators which are bilinear in the phase space fields but are, as a matter of fact, spatially non local.
Because of their wave-like dynamics Lorentz covariant fields are intrinsically spatially non local yet with a causally consistent time evolution. On the other hand, by definition symmetries of a dynamical system are transformations of its degrees of freedom that leave its equations of motion form-invariant, namely that map solutions into solutions to the same equations of motion. However dynamical systems may possess symmetries larger than those that leave their Lagrangian action invariant (up to total time derivatives), namely symmetry transformations which, when acting on the associated phase space (rather than on the system's configuration and conjugate momentum spaces separately), leave the system's Hamiltonian action invariant (up to total time derivatives). Such symmetry transformations mixing configuration space degrees of freedom with their conjugate momenta correspond generally to so-called dynamical symmetries 1 . Obviously Noether's (first) theorem (for global symmetries) applies likewise to the phase space parametrisation of a dynamics in terms of its Hamiltonian action, allowing for the identification of the Noether conserved charges generating global symmetries over phase space [10] through Poisson brackets. Consequently the general setting appropriate to studying the global symmetries of a system is that of its canonical Hamiltonian formulation.
In the case of a free field, say a bosonic one, its action being quadratic and its equation of motion being linear in the field, hence also its Hamiltonian equations of motion, the general transformations of the phase space fields leaving these equations invariant, hence also the Hamiltonian action, are a priori linear in the phase space fields. The corresponding (Noether charge) generators thus need to be bilinear in these fields. However because of the intrinsic spatially non local character of field dynamics, these generators need not necessarily be local in space (as are for instance the total energy, momentum and angular-momentum of the field when conserved), and may indeed be spatially bilocal, namely be given by a double spatial integral of a bilinear form in the phase space fields. This is thus the rationale to be implemented in the present paper when wanting to identify all the possible global symmetries of a free field dynamics. Specifically the analysis to be presented applies to the simplest case possible, namely that of a free and real bosonic scalar field in a 2+1 dimensional Minkowski space-time (generalisation to Minkowski spacetimes of arbitrary dimension is of course feasible). Furthermore in order to better understand the possible role of asymptotic behaviour at infinity for the existence of BMS symmetries (in the spirit of their original discovery), in the present work a spatially bounded domain will first 2 be considered, namely that of a disk of finite radius R > 0, D(R), in order to at least preserve -besides the invariance under global time translations (thus the conservation of the field's total energy) -the invariance under spatial rotations around a fixed center as well (and thus the conservation of the field's total angular-momentum). In other words the considered 2+1 dimensional spacetime has the topology of a plain cylinder R × D(R). However the spatial boundary condition to be imposed on the scalar field will be as general as possible in a manner consistent with rotational invariance, namely in terms of a radial Robin boundary condition of uniform Robin parameter along the disk boundary, allowing for a smooth transition between a Dirichlet and a Neumann boundary condition if required. This choice of boundary condition is to be implemented in the field dynamics through a specific boundary term in the Lagrangian action of the free field.
With the purpose as well of specifying notations and establishing results for later use, Section 2 discusses the free field dynamics based on an action principle, its Hamiltonian formulation and eventually its canonical quantisation. Section 3 then implements the rationale outlined above, to identify the complete ensemble of global symmetries of the considered free field dynamics, specifically for the flat disk spatial topology and geometry. The comprehensive identification of the complete global symmetry group is achieved in Section 4. While some final comments are presented in the Conclusion.

Free Field Dynamics and its Quantisation
Given the background 2+1 dimensional Minkowski spacetime of topology R × D(R), in natural units (c = 1 = ) let us thus consider a real scalar field φ(t, r, θ), (r, θ) being polar coordinates in the disk (0 ≤ r ≤ R and 0 ≤ θ < 2π). The dynamics of that field derives from the following choice of Lagrangian action which is purely quadratic in the field, where µ ≥ 0 is a real mass parameter while λ ∈ R is the real, constant and dimensionless Robin parameter determining the radial Robin boundary condition for the equation of motion of the free scalar field. In the above expression the first line is the bulk contribution to the action while the second line is a boundary contribution. Note that besides its invariances under global time translations and spatial rotations around the center of the disk D(R), this dynamics is also invariant under parity and time reversal transformations.

Classical Solutions
Obviously given the above action the variational principle, which requires the action to be stationary only up to a total time derivative, implies the usual Klein-Gordon wave equation for the scalar field in the disk, but subjected as well to the radial Robin boundary condition of parameter λ, Note that the value λ = 0 reproduces a Neumann boundary condition, while the limits λ → ±∞ imply a Dirichlet boundary condition for φ(t, r, θ). The Klein-Gordon equation being linear in the field, through a separation of variables it suffices to find a basis of solutions to construct its general solution by superposition. The lack of spatial translational invariance in the present instance however, implies that the usual travelling plane waves are not available for such a basis. A basis of solutions needs rather to be identified in terms of stationary waves. Given the invariances of the dynamics under global translations in time and global rotations in the disk with a 2π periodicity in the θ dependence, through a separation of variables such stationary waves of angular-momentum ℓ ∈ Z and positive angular frequency ω ℓ > 0 may be considered in the following (complex) form, Direct substitution into the Klein-Gordon equation then requires, where . This is the differential equation defining the real Bessel and Neumann functions of order ℓ, J ℓ (x) and N ℓ (x). Since the field φ(t, r, θ) must remain regular throughout the disk, only the Bessel functions of the first kind, J ℓ (x), may be retained as solution for the radial dependence of the field, namely f ℓ (r) = J ℓ (k ℓ r) up to a constant normalisation factor. Let us then also recall the parity property However the spectrum of wave number values k ℓ is restricted by imposing the Robin boundary condition, which leads to the following equation in k ℓ (J ′ ℓ (x) denotes the derivative of J ℓ (x) relative to its argument x), The properties of the positive roots of this equation in x R have been discussed in detail in Refs. [11,12]. For generic values of the Robin parameter λ and a given value of ℓ, the positive roots form an infinite discrete set of non degenerate values, x ℓ,n , increasing with n = 1, 2, · · · , the smallest of which is strictly positive, x ℓ,1 > 0. Correspondingly the spectrum of wave numbers solving the Robin boundary condition is thus given as k ℓ,n = x ℓ,n /R, with ω ℓ,n = k 2 ℓ,n + µ 2 . As a function of R and for a given value of ℓ, it is only at most for a finite set of isolated values for λ (if not a single such value) that the smallest root, x ℓ,1 , may be doubly degenerate. In this work we shall assume that such a degenerate circumstance is not encountered nor that accidental degeneracies in the roots x ℓ,n would occur for different values of |ℓ| (in which case otherwise the set of global symmetries to be identified for the free field theory would become slightly enlarged as compared to our final identification of the global symmetries of this system).
However, note that because of parity invariance, or equivalently time reversal invariance of the field, these roots are doubly degenerate for opposite non vanishing values of ℓ, as also follows form the fact that J −ℓ (x) = (−1) ℓ J ℓ (x).
In conclusion the full spectrum of solutions to the Laplacian eigenvalue problem on the disk, ∆ϕ(r, θ) = −k 2 ϕ(r, θ) with k ∈ R and k ≥ 0, subjected to the rotationally invariant radial Robin boundary condition of parameter λ, namely, is given by the following basis of functions, with x ℓ,n = k ℓ,n R, where N ℓ,n are (possibly complex) constant normalisation factors still to be determined. Requiring these functions to be orthonormalised in the disk, namely implies for the normalisation factors 3 , Finally because of the well known representation our choice of phase convention for the normalisation factors N ℓ,n is such that Besides the above orthonormality properties the functions ϕ ℓ,n (r, θ) also possess the following completeness relation as a genuine basis of functions in polar coordinates on the disk obeying the radial Robin boundary condition of parameter λ, Note that because of the properties under the change of sign in ℓ pointed out above, and the phase convention for N ℓ,n specified above, we also have, in a way similar to an analoguous property obeyed by the usual spherical harmonics on the 2-sphere. 3 Using the differential equation obeyed by Jn(kx) (with n ∈ N) it is possible to establish that . Orthogonality of the functions ϕ ℓ,n (r, θ) for different values of n (and a same value for ℓ) then follows by relying on the Robin boundary condition identity defining the roots x ℓ,n .
In conclusion the complete and general solution to the Klein-Gordon equation subjected to the considered radial Robin boundary condition is thus of the form, where the choice of normalisation is made for later convenience, while a † (ℓ, n) ≡ a * (ℓ, n) and a(ℓ, n) are complex valued integration constants corresponding to the excitation amplitudes and phases of each of the possible excitation modes (ℓ, n) of standing waves in the disk, which for the quantised field will correspond to the (normalised) Fock algebra creation and annihilation operators for the different quanta possible for this free field theory.
Within the Hamiltonian formulation of the system (see the next Subsection), in polar coordinates in the disk the momentum conjugate to the field, π(t, r, θ), is given as π(t, r, θ) = r∂ t φ(t, r, θ), namely, The relevant inverse relations are then, for the corresponding mode amplitudes,

Hamiltonian Formulation
Given the polar coordinate parametrisation of the Lagrangian action in (1) the momentum conjugate to the field configuration φ(t, r, θ), defined by π(t, r, θ) = δL/δ(∂ t φ(t, r, θ)), is obtained as, for which the Poisson brackets are canonical, with the only non vanishing (equal time) bracket being Correspondingly the canonical Hamiltonian of the system is given as, inclusive of the boundary term, in terms of which the Hamiltonian first-order action is expressed as It should be noted that in view of its definition the conjugate momentum field obeys a rotationally invariant radial Robin boundary condition of parameter (λ + 1), namely, These are thus the boundary conditions to be considered when solving the associated Hamiltonian equations of motion, of which the solutions will be expressed in terms of the basis of functions ϕ ℓ,n (r, θ) and rϕ ℓ,n (r, θ) (see (15) and (16)). Note that given the equation obeyed by the roots x ℓ,n = k ℓ,n R, we have indeed as it should, Through a careful treatment of the boundary term contribution in the evaluation of Poisson brackets with the above Hamiltonian H, it may readily be checked that the first-order and linear Hamiltonian equations of motion read as of which, given the appropriate Robin boundary conditions, the solutions are provided by the expressions (15) and (16). Note that as it should, these same equations of motion also derive directly from the variational principle applied to the above Hamiltonian action S H , inclusive of spatial boundary contributions that need to vanish as well.
Using the inverse relations in (17) and (18) as well as the orthonormality properties of the basis functions ϕ ℓ,n (r, θ), one may also readily establish that the above canonical Poisson brackets for the phase space fields φ(t, r, θ) and π(t, r, θ) translate into the following (only non vanishing) Poisson brackets for the mode amplitudes, Given the expression (21) and using the Robin boundary condition for φ, through partial integrations it is straightforward to establish the following purely bulk representation for the Hamiltonian which is spatially local in phase space and quadratic in the phase space fields, Then using the fact that, ∆ − µ 2 ϕ ℓ,n (r, θ) = 1 r ∂ r r∂ r + 1 r 2 ∂ 2 θ − µ 2 ϕ ℓ,n (r, θ) = −ω 2 ℓ,n ϕ ℓ,n (r, θ), ω 2 ℓ,n = k 2 ℓ,n + µ 2 , (28) as well as the orthonormality properties of the basis functions ϕ ℓ,n (r, θ), a straightforward evaluation establishes the following representation for the conserved total energy of the field, in which none of the factors a(ℓ, n) and a † (ℓ, n) have been permuted with one another.
Besides the global symmetry under constant translations in time for which the conserved Noether charge is the total energy or the Hamiltonian H, the system also possess a global symmetry under constant rotations of the disk around its center at r = 0, for which the conserved Noether charge is the total angular-momentum of the field, L. A simple analysis identifies its spatially local phase space expression which is quadratic in the phase space fields, A straighforward substitution in terms of the mode amplitudes then finds,

The Quantised Field
In terms of the results established above within the Hamiltonian formulation of the system, its canonical quantisation is straightforward. Its Hilbert space of quantum states is spanned by the tensor product of the Fock spaces generated from a Fock vacuum by the Fock algebras of annihilation and creation operators, a(ℓ, n) and a † (ℓ, n), obeying the Fock algebras (with = 1), a(ℓ 1 , n 1 ), a † (ℓ 2 , n 2 ) = δ ℓ 1 ,ℓ 2 δ n 1 ,n 2 I, Time dependence of the system is generated through the total quantum Hamiltonian, which we choose to be given by the normal ordered form of the classical quantity, namely, In particular in the Heisenberg picture the time dependent quantum field operatorsφ(t, r, θ) and π(t, r, θ) remain then given by the same expressions as in (15) and (16), now in terms of the Fock operators.
Likewise the total angular-momentum operator of the system is represented by the normal ordered expression,L Hence the action of a † (ℓ, n) onto any quantum state of the system adds one more quantum of energy ω ℓ,n and angular-momentum ℓ in the standing wave mode of the field of quantum numbers (ℓ, n) related to the root x ℓ,n of the equation expressing the radial Robin boundary condition for the angular-momentum mode ℓ. Because of the parity and time reversal invariance property that x −ℓ,n = x ℓ,n , note well that the energy spectrum of the system is doubly degenerate under the exchange ℓ ↔ −ℓ for ℓ = 0. In other words one may as well write,

Bilocal Global Dynamical Symmetries and their Charges
Given the rationale outlined in the Introduction, let us now consider the possibility that quantities which are bilocal and bilinear in the phase space fields, hence bilinear in the creation and annihilation operators, could generate (as the would-be associated conserved Noether charges) global symmetries of the free field dynamics. In a compact notation such that "1" or "2" stand for the pairs of coordinates (r 1 , θ 1 ) and (r 2 , θ 2 ), respectively, such candidate conserved charges are of the general form where the triplet of time independent functions f (1; 2) ≡ f (r 1 , θ 1 ; r 2 , θ 2 ), g(1; 2) ≡ g(r 1 , θ 1 ; r 2 , θ 2 ) and h(1; 2) ≡ h(r 1 , θ 1 ; r 2 , θ 2 ) are real functions to be identified. In particular note that the functions f (1; 2) and g(1; 2) are symmetric under the exchange 1 ↔ 2, while h(1; 2) does not possess any specific symmetry under that exchange, Through Poisson brackets with the fields φ(t, r, θ) and π(t, r, θ) such charges Q(t) generate the following linear transformations mixing the phase space fields, In a form more convenient for later use, these transformations may also be expressed as, δφ(t, r, θ) = R 0 dr 1 r 1 2π 0 dθ 1 g(r, θ; r 1 , θ 1 ) 1 r 1 π(t, r 1 , θ 1 ) + 1 r 1 h(r 1 , θ 1 ; r, θ) φ(t, r 1 , θ 1 ) , 1 rr 1 f (r, θ; r 1 , θ 1 ) φ(t, r 1 , θ 1 ) + 1 r h(r, θ; r 1 , θ 1 ) 1 r 1 π(t, r 1 , θ 1 ) .(42) A first series of restrictions on the functions f , g and h is implied by the necessity that the transformed fields still obey their relevant Robin boundary conditions of parameter λ or (λ + 1) as the case may. Consequently one must require, Note that since φ(t, r, θ) and π(t, r, θ)/r both obey the radial Robin boundary condition with a same parameter λ, as a matter of fact so do the functions f (1; 2)/(r 1 r 2 ), g(1; 2) and h(1; 2)/r 1 , which is equivalent to the above boundary conditions for the triplet of functions (f, g, h). Different strategies may be considered in order to determine under which conditions the charges Q(t) generate global symmetries of the field dynamics. By definition a symmetry leaves the equations of motion form-invariant; through a direct substitution into the Hamiltonian equations of motion of the above transformations δφ and δπ one may find specific restrictions that the functions f , g and h would have to meet. However there would then still remain the necessity to check under which conditions the composition (through Poisson brackets or commutators) of transformations generated by such charges Q(t) for different such triples of functions (f, g, h) closes on the same class of transformations.
A rather more efficient approach considers directly the Hamiltonian action of the system, S H . If the above transformations define a global symmetry leaving the equations of motion forminvariant, necessarily they leave the action invariant possibly up to a total time derivative. If a total time derivative is indeed thereby induced it could be that the associated Noether charge, even though conserved under its complete equation of motion 4 , possesses an explicit time dependence as a function over phase space [10], which in turn may lead to a classical central extension of the algebra spanned by the Noether generators of such global symmetries of the system 5 .
Let us thus consider the Hamiltonian action (22) with the Hamiltonian given as in (27) (without a boundary term), and determine its linearised variation under the transformations δφ and δπ in (39) and (40). Through a series of careful integrations by parts in r and in θ, and by exploiting the Robin boundary conditions that the triplet of functions (f, g, h) must satisfy, the linearised variation of the Hamiltonian action reduces to the expression, In order that this variation reduces solely to a total time derivative contribution, in addition to those properties already listed above the triplet of functions (f, g, h) must satisfy as well the following further properties, Incidentally recall that f (1; 2) needs to be symmetric under the exchange 1 ↔ 2. Hence we have established so far that Note that the functions g(1; 2), h(1; 2)/r 1 , and f (1; 2)/(r 1 r 2 ), all satisfy radial Robin boundary conditions of parameter λ, in both their dependencies in r 1 and r 2 . Given that the functions ϕ ℓ,n (r, θ) provide a basis of functions defined on the disk D(R) and obeying this very Robin boundary condition, it proves relevant to expand in that basis the symmetry transformation functions just listed in order to solve the different conditions they must meet, in the following form, g(r 1 , θ 1 ; r 2 , θ 2 ) = ℓ 1 ,n 1 ℓ 2 ,n 2 α(ℓ 1 , n 1 ; ℓ 2 , n 2 ) ϕ ℓ 1 ,n 1 (r 1 , θ 1 ) ϕ ℓ 2 ,n 2 (r 2 , θ 2 ), and, 1 where the complex coefficients α(ℓ 1 , n 1 ; ℓ 2 , n 2 ) and β(ℓ 1 , n 1 ; ℓ 2 , n 2 ) are to be restricted in order to solve for all the necessary properties that the functions g(1; 2) and h(1; 2)/r 1 must possess, while each of the double sums runs over ℓ ∈ Z and n = 1, 2, 3, · · · .
As a consequence all the coefficients β(ℓ 1 , n 1 ; ℓ 2 , n 2 ) vanish identically unless ℓ 1 = ℓ = −ℓ 2 and n 1 = n 2 , in which case one finds, with ℓ ∈ Z and n = 1, 2, 3, · · · , β(ℓ, n; −ℓ, n) = iβ(ℓ, n), with β(ℓ, n) ∈ R and β(−ℓ, n) = −β(ℓ, n), where β(ℓ, n) is a collection of arbirary real constants (note that β(0, n) = 0). Thus finally, Having identified the triplet of functions (f, g, h) a direct substitution into (47) determines the explicit form for the general generator of all global (and dynamical) symmetries of the free field φ(t, r, θ) in terms of its creation and annihilation operators. A patient but straightforward calculation then establishes that, once brought into normal ordered form in the case of the quantum operator, Note that this operator is such that Q † (t) = Q(t), as it should. Furthermore, since one obviously has [Q(t),Ĥ] = 0 (for which the degeneracy property ω −ℓ,n = ω ℓ,n is crucial), even though our approach accounted for the possibility that the generator of the general global symmetry could possess an explicit time dependence as a function defined over phase space, it turns out that in fact the conserved quantity Q(t) does not possess any explicit time dependence.

The Complete Global Symmetry Group
The algebra of all global symmetries of the 2+1 dimensional free scalar field theory on the disk D(R) is thus spanned by the operators N ℓ,n and Q ℓ,n , respectively the number operator and the angular-momentum flipping operator for each of the (ℓ, n) modes of the stationary waves of the field and its conjugate momentum inside the disk. It should be clear these transformations do indeed transform a solution to the Klein-Gordon equation into another solution to the same Klein-Gordon equation, and this without changing the total energy of that field configuration since ω −ℓ,n = ω ℓ,n . Furthermore in particular, note that the conserved total energy and angularmomentum of the system are part of that large dynamical symmetry, which is much larger than the finite dimensional global symmetry group of the underlying spacetime geometry with the topology of R × D(R), namely time translations and disk rotations of which, when acting on the field and its conjugate momentum the generators are, respectively 6 , H = ℓ,n ω ℓ,n N ℓ,n ,L = ℓ,n ℓ N ℓ,n . (61) The commutator algebra generated by all these conserved Noether charges of course closes, and is given by, In particular it also follows that, Ĥ , Q ℓ,n = 0, L , Q ℓ,n = 2ℓ Q ℓ,n , as it should of course, while furthermore, [N ℓ 1 ,n 1 , a(ℓ 2 , n 2 )] = −δ ℓ 1 ,ℓ 2 δ n 1 ,n 2 a(ℓ 1 , n 1 ) , N ℓ 1 ,n 1 , a † (ℓ 2 , n 2 ) = +δ ℓ 1 ,ℓ 2 δ n 1 ,n 2 a † (ℓ 1 , n 1 ), [Q ℓ 1 ,n 1 , a(ℓ 2 , n 2 )] = −δ ℓ 1 ,ℓ 2 δ n 1 ,n 2 a(−ℓ 1 , n 1 ) , Q ℓ 1 ,n 1 , a † (ℓ 2 , n 2 ) = +δ ℓ 1 ,−ℓ 2 δ n 1 ,n 2 a † (ℓ 1 , n 1 ).
In order to identify the symmetry group associated to the above algebra generated by N ℓ,n and Q ℓ,n let us take as a clue the expressions for the total Hamiltonian and angular-momentum operators in (35) and (36). Clearly quantum sectors with different values for n = 1, 2, 3, · · · are all decoupled form one another, and so are the sectors with different values of |ℓ| for ℓ ∈ Z. However given a fixed value for n, the sectors with opposite values of ℓ = 0 and −ℓ are coupled to one another through the action of the global symmetries generated by N ℓ,n and Q ℓ,n .
Let us thus first consider the sector with ℓ = 0 and a given value for n = n 0 . Since such a sector does not contribute to the total angular-momentumL, is left invariant by the sole symmetry generator N 0,n 0 = Q 0,n 0 = a † (0, n 0 )a(0, n 0 ) and contributes to the total energy asĤ(0, n 0 ) = ω 0,n 0 N 0,n 0 , any such sector is equivalent to that of an ordinary one-dimensional harmonic oscillator of angular frequency ω 0,n 0 , of which the global symmetry is the U(1) 0,n 0 phase symmetry generated by the number operator N 0,n 0 . Hence the complete global symmetry of all ℓ = 0 sectors of the free scalar field in the disk D(R) is ∞ n 0 =1 U (1) 0,n 0 . Consider now a specific non vanishing and positive value of ℓ = 1, 2, 3, · · · as well as a given value for n = n ℓ . Let us then adapt the notations as follows, in order to emphasize the analogy to be highlighted, and consider the following combinations of the symmetry generators for the sectors (ℓ, n ℓ ) and (−ℓ, n ℓ ), Obviously 2ω ℓ,n T 0 is the total contribution of the sectors (ℓ, n) and (−ℓ, n) to the total energy of the system, while 2ℓT 3 is that to its total angular-momentum. Once again T 0 generates a global U(1) ℓ,n ℓ phase symmetry for these two sectors of the system. However one has furthermore or equivalently, in which one recognises the SU(2) Lie algebra. Indeed the Fock algebras (a ± , a † ± ) are precisely those of the two dimensional spherically symmetric harmonic oscillator in the Euclidean plane (in the helicity basis), which is well known to possess a SU(2) dynamical symmetry with the above SU(2) algebra. Thus the sectors (ℓ, n ℓ ) and (−ℓ, n ℓ ) with ℓ = 0 of the free scalar field in the disk D(R) possess a global and dynamical SU(2) ℓ,n ℓ symmetry, which is not spatially local in phase space.
Consequently one comes to the final conclusion that the complete global symmetry group of the free scalar field in the disk D(R) -which for most of it is a dynamical and spatially non local symmetry group -is identified as the following infinite countable symmetry group The finite dimensional symmetry group of the underlying spacetime geometry, namely the direct product of global time translations and space rotations, is a specific subgroup of the above, with the following two abelian generators which are also spatially local in phase space, To make the above points somewhat more explicit, let us also consider now the finite symmetry transformations generated by each of the generators N ℓ,n , Q ℓ,n and Q † ℓ,n = Q −ℓ,n to check that indeed they map solutions into other solutions to the equations of motion, by redefining and mixing the mode amplitudes a(ℓ, n) and a † (ℓ, n) and thus the latter's contributions to the field and its conjugate momentum, φ(t, r, θ) and π(t, r, θ), in the Heisenberg picture.
Beginning with the Hermitian number operators N ℓ,n = N † ℓ,n , let us consider a specific but otherwise arbitrary choice of values ℓ 0 ∈ Z and n 0 ∈ N + and the number operator N ℓ 0 ,n 0 = a † (ℓ 0 , n 0 )a(ℓ 0 , n 0 ) corresponding to the mode (ℓ 0 , n 0 ) of the field. Since the only sector (ℓ, n) and its Fock algebra which does not commute with N ℓ 0 ,n 0 is that of the mode (ℓ 0 , n 0 ), all finite and global symmetries generated by that number operator through the action of the unitary operator U (α) ≡ e iαN ℓ 0 ,n 0 , α ∈ R,

Conclusions
In view of the double degeneracy in the energy spectrum under the exchange ℓ ↔ −ℓ for ℓ = 0 as made explicit in (35) and (36), and given the hindsight gained through the present analysis, it would appear rather obvious that indeed the system possesses as finite global (and dynamical) symmetries the group identified above, namely However that there do not exist any further global symmetries is less obvious, a conclusion which requires an approach such as the one having been used herein and based on the rationale outlined in the Introduction. This is not to say that there do not exist other conserved quantities for the system. Rather it means that there do not exist conserved quantities other than N ℓ,n and Q ℓ,n that would also be generators of additional continuous global symmetries of the system.