A First-Quantized Model For Unparticles and Gauge Theories Around Conformal Window

We first quantize the action proposed by Casalbuoni and Gomis in [Phys. Rev. D \textbf{90}, 026001 (2014)], an action that describes two massless relativistic scalar particles interacting via a conformally invariant potential. The spectrum is a continuum of massive states that may be interpreted as unparticles. We then obtain in a similar way the mass operator for a deformed action in which two terms are introduced that break the conformal symmetry: a mass term and an extra position-dependent coupling constant. A simple Ansatz for the latter leads to a mass operator with linear confinement in terms of an effective string tension $\sigma\,$. The quantized model is confining when $\sigma\neq0$ and its mass spectrum shows Regge trajectories. We propose a tensionless limit in which highly excited confined states reduce to (gapped) unparticles. Moreover, the low-lying confined bound states become massless in the latter limit as a sign of conformal symmetry restoration and the ratio between their masses and $\sqrt\sigma$ stays constant. The originality of our approach is that it applies to both confining and conformal phases via an effective interacting model.


I. INTRODUCTION
It is known that some asymptotically free gauge theories with N f light fermion flavours and SU (N ) gauge group have a conformal window, i.e., there exists an energy range in which the beta function vanishes. This conformal window occurs at various values of N f and N , depending on the fermion representation, see [1] for an extensive list of examples obtained by inspection of the two-loop beta function. For example, SU (3) gauge theory with 12 light quarks flavours in the fundamental representation is conformal, as confirmed by a five-loop and nonperturbative calculations [2,3]. SU (2) gauge theory with 2 quark flavours in the adjoint representation can also be mentioned [4,5], or 2 quark flavours in the two-index symmetric representation [6]. Note that conformal windows are expected to be present for other gauge groups than SU (N ), like Sp(2N ) and SO(N ) [7]. Many evidences showing the existence of a conformal window for specific gauge theories have been found by resorting to lattice QCD methods, see e.g. the reviews [8][9][10][11][12]. General algorithms actually exist for SU (N ) theories with quarks in arbitrary representations [13].
To our knowledge, no effective model -we mean a simple enough action so as to allow for analytical calculations -has been proposed to mimic confining gauge theories when they approach the conformal window starting from a confining phase. Our starting point has been proposed in Ref. [14], where an action describing two scalar relativistic particles with conformal invariant interactions has been presented. The latter action reads with α a dimensionless parameter, x µ i = x µ i (τ ) the parametric equations for the two particles in the (D + 1)-dimensional Minkowski spacetime with metric η = diag(− + · · · +) in inertial coordinates, τ an evolution parameter and r µ := x µ 2 − x µ 1 . For any vector with components v µ , v 2 := η µν v µ v ν =: v · v is the squared norm in the Lorentzian sense. In Section II we first review and quantize the action S C . In particular, while reproducing some steps of the canonical analysis of [14], we also detail some issues that were not presented therein for the sake of conciseness but that we need for the quantization of the model. We link the spectrum obtained after quantization to unparticles, originally introduced as a nontrivial scale-invariant sector in low-energy effective field theories [15], see also [16,17].
The model is generalized in two ways that break conformal invariance. In Sec. III we introduce a mass term for the interacting particles. The quantization leads to unparticle spectrum with a mass gap. Then a confining interaction is introduced in Sec. IV via a change of the form α 2 → α 2 U 2 (r 2 ) . The quantization of our model is then performed and the spectrum is analytically computed. The transition from the confining phase to the conformal phase is finally studied and unparticles are shown to emerge from the confined spectrum. The idea that unparticles may be "hidden" in a gauge theory's conformal window has been developed in [18]. We finally argue that the present model provides a toy model to illustrate the latter proposal.

Let us review the conformally-invariant action presented by Casalbuoni and Gomis in
Refs. [14,19]. The aim of their presentation was to show the existence of higher-spintype conserved currents in their model while ours is rather to prepare the ground for the quantization of our model.
In the case of two massless particles, the authors of [14] start from the action with e i (i = 1, 2) the einbeins, and add an interaction term that couples the two massless particles in a conformally-invariant way: where we recall that r 2 := r µ r ν η µν =: r · r and where we have chosen the potential to be repulsive. In the above action, the coupling constant α is dimensionless. Since the action is manifestly Poincaré invariant it suffices to prove that it is invariant under dilations and special conformal transformations. As for the dilations, it is easy to see that the transformations x µ i → λx µ i , e i → λ 2 e i preserve the action. The invariance under special conformal transformations is established by defining them as a succession of a Poincaré translation followed by an inversion and another Poincaré translation, where the inversion is the transformation The interested reader may find additional information in the following lecture notes. The conformal invariance of the action is then readily checked.
Note that one can eliminate the two auxiliary variables e 1 (τ ) and e 2 (τ ) by virtue of their own field equations, which results [14] in the incarnation (1) of the action. We mention for completeness Ref. [20] in which it is shown that action (3) may be modified in a simple way so that it models two particles interacting conformally in Snyder space.
The action (3) does not depend on the derivativesė i , hence there are two primary constraints: where π i is the conjugate variable to e i and where the symbol "≈" denotes a weak equality, i.e., an equality which holds on the constraint surface. One then derives the canonical Hamiltonian associated with the Lagrangian action (3): as well as the total Hamiltonian √ e 1 e 2 r 2 + λ 1 π 1 + λ 2 π 2 .
At this stage it is worth saying that the same analysis can be performed even if α is a function α(r): Although breaking conformal symmetry, this case is important for our purpose and will be discussed in Sec. IV.
Coming back to the Hamiltonian (11), sinceλ 2 is an arbitrary function, one may take e 1 π 1 + e 2 π 2 as our new primary first-class constraint. Its Poisson brackets with π 1 and π 2 is weakly zero while it is strongly zero with φ 1 and φ 2 . Furthermore, H T is first-class since it does not explicitly depend on time. From the fact that the bracket of two first-class functions is first-class and the computation one derives that γ 2 is another first-class constraint. To summarise, there are two first-class (FC) constraints (one is primary and the other is secondary) and two second-class (SC) constraints (one is primary and the other is secondary): FC : γ 1 := e 1 π 1 + e 2 π 2 (primary) , γ 2 := e 1 φ 1 + e 2 φ 2 − Cπ 1 (secondary) , (13) SC : χ 1 := π 1 (primary) , One has The number f of degrees of freedom is given by the number of phase-space variables minus twice the number of first-class constraints minus the number of second-class constraints. The number of degrees of freedom for the conformally-invariant interacting system (3) is therefore given by f = 4(D + 1) + 4 − 4 − 2 = 4D + 2 , which differs from the counting obtained for the non-interacting system that produces f = 4D . If one considers the free limit α → 0, the limiting value should therefore not be considered. This is also clear from the form (1) of the action. In fact, one can show that the Dirac conjecture is not satisfied by the constrained system at hand. It turns out that first-class constraints are all gauge symmetry generators if four conditions established in Chapt. 3 of [21] are respected. One of these conditions is that the secondary second-class constraints should not appear in the Poisson bracket of the first-class constraints with the primary second-class constraint. With γ a the first-class constraints and χ α the second-class constraints, in general one has that and the condition mentioned above is that the matrix elements C β =β 1 aα 1 must be equal to zero, where α 1 and β 1 refer to the primary constraints. In our case, this condition means that the quantities C φ 1 σ i π 1 should be null quantities. However, showing that C φ 1 σ 2 π 1 = 0 . Therefore, at least one condition imposed in [21] is not satisfied, which implies that the Dirac conjecture is not true in our case. The first-class constraints are not necessarily all generators of gauge transformations and one must use the chain algorithm of [49] to determine these generators; see e.g. Appendix A for some details about it.
The field equations are obtained by taking the Poisson bracket of the dynamical variables with the total Hamiltonian (11). They explicitly reaḋ The second class constraints can be dropped by adopting the Dirac bracket defined in terms of the inverse of the matrix where (χ α ) α=1,2 = (π 1 , φ 1 ) denotes the two second class constraints and D := . One can then consider the reduced phase space where one strongly sets π 1 and φ 1 to zero. The first-class constraint γ 2 ≈ 0 therefore leads to the constraint φ 2 ≈ 0 , since the einbein e 2 is required to be non vanishing. By using that φ 1 is strongly set to zero, the equation φ 2 ≈ 0 yields a constraint where the einbeins e 1 and e 2 have been eliminated [14]: This equation is one of the main results of [14], that terminates their canonical analysis. We pursue their analysis in order to derive a more tractable constraint in view of the quantization. We express (24) in the relative (−) and centre-of-mass (+) variables that satisfy the canonical Poisson bracket relations In the new coordinate system and once the second class constraints (χ α ) α=1,2 = (π 1 , φ 1 ) have been strongly set to zero, the resulting first-class constraint (24) reads In these coordinates, the canonical equations of motion (18)-(22) read: Obviously, the total momentum p µ + is preserved by the dynamics, which simply reflects the invariance of the system under constant spacetime translations. Since Eq. (27) is not totally tractable yet, we will completely fix the gauge by finding two gauge-fixing conditions C a = 0 such that the bracket matrix with entries M ab := {C a , γ b } is non-degenerate, effectively resulting in a second-class system. We propose the following gauge-fixing conditions and readily check that the bracket matrix is invertible, as we wished. In the gauge C 1 = 0 and in the Lorentz frame where the (preserved) total momentum is p µ + = (M, 0, 0, 0) , we have p µ − = (0, p − ) and p 2 1 = p 2 2 : there is a balanced energy distribution among the two particles. Although this case will not be investigated here, it has been shown in [22] that massless bound states can also be considered by using similar methods.
The first two field equations of (31) indicate that one can further impose the condition which leads toq i + = 0 andq 0 − = 0 . Indeed, let us assume e 1 = e 2 and therefore e 1 = 1 = e 2 on account of C 2 = 0 . The two conditions e 1 = 1 = e 2 can obviously be reached by virtue of the primary constraints (5). We call this gauge the unit einbein gauge. Eq. (29) shows that q − .p + is a constant that one considers to be zero in order to have q 0 − = 0 in the Lorentz frame adopted. From the first two equations (31), the constraint e 1 = e 2 is consistent provided is zero when e 1 − e 2 = 0 . It is straightforward to check that it is indeed the case from the fact that q − .p + is zero, as we have just shown. We therefore set q − = (0, q − ) . The three-vector q + is set to zero, in accordance with translation invariance of the system. From the equations of motion and in the gauge chosen, it is clear that the only dynamical variables are q − and p − .
In the gauges and Lorentz frame we have chosen, Eq. (27) becomes Instead of extracting the square root of the above equation, with the ambiguity in which branch to pick, we recall that the remaining first-class constraint γ 2 ≈ 0 now reads φ 2 ≈ 0 that leads, in the gauges we have chosen: Finally one is led to the following dispersion relation

II.2. Quantization
In a Schrödinger quantization scheme, Eq. (37) defines the eigenequation The spherical symmetry of this operator allows to work with hyperspherical coordinates with Y ,ma (Ω D ) the spherical harmonics in D dimensions, ∈ N, m a ∈ Z and a = 1, . . . , (D − 1). Explicit forms can be found for example in [23]. More precisely, the squared mass operator is a Schrödinger Hamiltonian with repulsive inverse-squared potential: Hence the mass spectrum is continuous and the eigenstates are scattering states: where µ ∈ R + and where J λ+ D 2 is a Bessel function of the first kind. The generalized angular-momentum index λ is defined by where λ > 0 guarantees a solution that is regular at the origin. Because α 2 > 0 , the radial function actually never vanishes at the origin. The interested reader may find in [24] a detailed discussion of the inverse-squared potential in quantum mechanics.
Our model in the conformal phase contains a continuum of bosonic states with arbitrary mass M ≥ 0. We identify this continuous spectrum with unparticles, originally introduced as a nontrivial scale-invariant sector in low-energy effective field theories [16]. As discussed in [25][26][27][28], such an unparticle sector arises from a continuum of scalar fields with arbitrary mass; an effective Lagrangian of the form L ∼ ∂ µ φ −δ ∂ µ φ with δ > 0 is then found for a scalar unparticle [25,26]. An originality of the present work is to provide a realization of unparticles as binary states of two interacting massless particles through action (3).

III. MASSIVE PARTICLES AND GAPPED UNPARTICLES
We first recall the action for a massive relativistic particle in Minkowski spacetime: where the variable e is required to be nonvanishing. The above action has a smooth massless limit m → 0 .
The whole constraint analysis that we have reviewed from [14] and presented above can be repeated in the case one adds a mass term to the original action. By considering the action S (m) that clearly breaks conformal invariance through the presence of the mass terms, it is straightforward to see that both the number and nature of the constraints remain unchanged compared to the massless case reviewed in great details above. There remains two first-class constaints {γ In the massless case one had the constraint (35). In presence of the mass terms in (44), following the same procedure, one is led to the dispersion relation and, after quantization, to the spectrum (41) with M 2 = µ 2 +4m 2 . The continuous spectrum is bounded from below by M = 2m .
It has previously been noticed that coupling an unparticle to an electroweak sector leads to an unparticle spectrum with mass gap via a kind of Brout-Englert-Higgs mechanism [29].
Note also that unparticle models with a mass gap are nothing but hidden-valley models [30]. We can finally mention Ref. [31], in which gapped unparticles emerge as continuous Kaluza-Klein modes of a five-dimensional model with a brane. The action (44) may be seen as a first proposal to generate gapped unparticles from two massive interacting particles. It is worth mentioning that such an action might be used in the modelling of near-threshold neutral charm meson molecules, since it has been argued in [32] that such states (like the X(3872) which is very close to the D * 0D 0 threshold) may be seen as unparticles.

IV. FROM THE CONFINING TO THE CONFORMAL PHASE
As starting point we propose the action Both the mass term and the extra interaction term with the function U (r 2 ) break conformal invariance. One obviously recovers action (3) in the limit m → 0 and by setting U = 1 .
The replacement of α 2 by the position dependent coupling α 2 U 2 (r 2 ) breaks conformal invariance but do not lead to drastic changes in the canonical constraint analysis. Although there appear terms proportional to U (r 2 ) from the Poisson bracket of the total Hamiltonian with the secondary constraints, one is still able to ensure that the latter constraints are preserved during dynamical evolution by identically solving the equation φ 3 = 0 for the function λ 1 . The latter function takes a more complicated form than the one given in the massless case reviewed in Section II. Here we obtain where the function C is given in (10). The additional terms proportional to the derivative of the function U (r 2 ) make no difference for the rest of the canonical analysis. Proceeding as in the previous sections, it is straightforward to find the following dispersion relation: We now use the Ansatz where the term in q 4 − can be seen as the first nontrivial term in the power expansion of any function U 2 (q 2 − ): a term in q 2 − would only redefine m. As we will show in the following, (49) mimics a linear confinement, typical of (3 + 1)−dimensional Yang-Mills theories in their confining phase. We therefore assume that the action (46) with potential (49) is a relevant effective model to describe binary states in gauge theories around their conformal window, just as a Nambu-Goto Lagrangian is a relevant effective model for light mesons (ρ, f 0 ,. . . ) in the confined phase, see e.g. [33,34].
As in Section II, we may set q 0 − = 0 , so that q 2 − = q 2 − . With the Ansatz (49), the mass operator (48) is now a D−dimensional harmonic oscillator with arbitrary angular momentum: Its spectrum reads [23] with L α n the generalized Laguerre polynomials and λ given by Eq. (42). The spectrum contains massive states showing Regge trajectories, i.e., M 2 ∼ or n at large or n. Such a behaviour is observed experimentally in light meson spectroscopy, see e.g. [35] and references therein. For this reason the massive states we observe in the confined phase will be referred to as "mesons".
An more general ansatz of the form U 2 (q 2 − ) = 1 − 4 α 2 δ 2 q 2 − + 4 α 2 σ 2 q 4 − leads to the mass spectrum The mass scale δ could be used to fine-tune the ground-state mass (n = = 0). For example, the value δ 2 = 2σ(1 + (D − 2) 2 + α 2 ) + 4m 2 leads to a massless ground state. However, such a value causes U 2 (q 2 − ) to be negative for some values of q 2 − since δ 2 > ασ. In our model, the existence of a massless ground state demands to drop the positivity of the potential.
Let us now focus on the low-lying confined spectrum (n and finite) when approaching the conformal window, that is in the limit σ → 0 and m → 0. According to (50), the mass of low-lying states will go to zero as already suggested in [36][37][38]: Light mesons are expected to become massless as a signal of conformal symmetry restoration. Their masses scale as √ σ, which is a behaviour observed in lattice QCD in the case of SU (2) QCD with one adjoint Dirac quark flavour: The ratios M/ √ σ are found to be constant for the lightest (pseudo)scalar states (mesons and glueball) as the fermion mass goes to zero in order to restore conformal invariance [39]. Moreover, it is observed that the mass ratios of two meson masses are constant near the conformal window as suggested by Eq. (50). This feature has been observed in SU (2) gauge theory with two adjoint quarks [40].
It has to be noticed that the methodoloy used in most of the lattice QCD studies of theories with conformal window is to choose N , N f and the quark representation such that the conformal window is a priori reached. The quark mass m starts from a nonzero value and the conformal window is reached by taking the limit m → 0; the masses of bound states scale as m and [42] Γ(n + 1) The spectrum (50) approaches to our unparticle sector (41) as n → +∞ up to the rescaling R → √ n R. The harmonic oscillator functions are indeed normalized to unity, which leads to the vanishing of R as σ → 0 since scattering states can only be normalized to δ(µ − µ ) in principle.
Our limit (53) is actually consistent with the results of [43] showing that, in field theory, an unparticle sector can be generated by a tower of massive states with mass M 2 n = ∆ n when the mass spacing parameter ∆ goes to zero.

V. CONCLUDING COMMENTS
The action (46) appears to be an interesting toy model to describe the transition from confining to conformal phases of a Banks-Zaks-type gauge theory [15] in terms of binary states. It predicts that the low-lying bound states in the confining phase become massless when approaching the conformal window with a universal behaviour, the ratios of masses and √ σ being a constant. Highly excited radial states give rise to a unparticle sector in the conformal phase. The unparticle spectrum has a mass gap or not depending on whether the conformal symmetry is broken or not by a mass term. A schematic drawing of this picture is given is Fig. 1.
Notice that the radial wave equation we derived in the present paper also appears from a bottom-up AdS/CFT perspective, see e.g. [44]. Interestingly, in the paper [45] where they consider a U (1)-invariant gauge theory in AdS 4 propagating a vector field and a complex scalar field, a nontrivial profile for the scalar field near the conformal infinity of AdS 4 introduces two dimensionful constants m and v that break conformal symmetry of the dual theory. Various limits where m → 0 and v → 0 relate a discrete spectrum with a massless Goldstone boson to a continuous Banks-Zaks-type spectrum for the scalar conformal operator in the dual CFT 3 . The approach of [45] enables a simultaneous discussion of spontaneous and dynamical symmetry breaking in a CFT at strong coupling. Some bridges may presumably be drawn between our model and the AdS/QCD framework, although it is out of the scope of the present paper. The interested reader may find a detailed discussions of mesons in AdS/QCD in the review [46].
One may finally wonder to what kind of gauge theory our toy model best applies. Similarities with SU (2) QCD with one adjoint quark [39] have been commented in the text. However, our model is a priori better suited to model gauge theories with scalar matter than with fermionic matter, the variables x i (τ ) being then identified with two "particles" of scalar matter. The existence of a conformal window is not limited to gauge theories with fermionic matter, as discussed in Appendix B. For example, a lattice study of N f = 5 scalar SU(2) Yang-Mills theory is affordable with current computers and could be performed to check the present model. We hope that such results will become available in the future.
Although we did not explicitly mention it, the present work and quoted references focus on zero-temperature gauge theories. The interplay between conformal window crossing and confinement/deconfinement transition may lead to new interesting phenomena, as discussed in [47]. For example, a correspondence has been found between conformal QCD with N f = 7 and QCD with N f = 2 at T ≈ 2 T c (T c is the deconfinement temperature). To what extent the action (46) can be generalized at finite-temperature and bring relevant information is a problem that we leave for future works.

ACKNOWLEDGMENTS
We thank Zhenya Skvortsov for a discussion about unparticles.
Appendix A: The chain algorithm for gauge generators We start from the constraint system (13)- (14) and seek for the generators of gauge transformations. According to the chain algorithm [49] for the determination of the gauge generators, there is just one of them because we have only one primary first-class constraint.
Also, the consistency algorithm stops after the secondary constraints: the desired generator consists in the sum of up to two terms. The constraint γ 1 itself cannot be a generator since {γ 1 , H T } is secondary. The generator is built with the first class constraints: G =ε(τ )γ 1 + ε(τ )(aγ 1 + bγ 2 ) , where a and b are arbitrary functions at this stage. Our goal is to fix these functions in order to give G the properties of a generator of gauge transformation. Since the rules of the chain algorithm have to be respected, one imposes: The bracket {γ 1 , H T } = γ 2 implies that b must be equal to −1 . Then, one knows that {γ 2 , H T } =λ 2 γ 2 , thus one concludes that the other condition implies that a must be equal It can be then useful to notice the following relation: ⇒ė 1 π 1 +ė 2 π 2 = C π 1 +λ 2 γ 1 .
Applying the methodology of [1] but neglecting chiral symmetry issues, one may search for a conformal window in theories such that β 0 < 0 and β 1 > 0, i.e., with a number of flavours such that 34C 2 (adj) 2 T (R) (C 2 (adj) + 6C 2 (R)) the coupling constant g * at which it is observed being equal to g * 2 = − 16π 2 β 0 β 1 . The equation (B1) admits nontrivial solutions. A simple example is the case of matter in adjoint representation, for which a conformal phase appears if 5 < N f < 22.