Phenomenology at the LHC of composite particles from strongly interacting Standard Model fermions via four-fermion operators of NJL type

A new physics scenario shows that four-fermion operators of Nambu-Jona-Lasinio (NJL) type have a strong-coupling UV fixed point, where composite fermions F (bosons Π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Pi $$\end{document}) form as bound states of three (two) SM elementary fermions and they couple to their constituents via effective contact interactions at the composite scale Λ≈O\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda \approx {{{\mathcal {O}}}} $$\end{document}(TeV). We present a phenomenological study to investigate such composite particles at the LHC by computing the production cross sections and decay widths of composite fermions in the context of the relevant experiments at the LHC with pp collisions at s=13\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}=13$$\end{document} TeV and s=14\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}=14$$\end{document} TeV. Systematically examining all the different composite particles F and the signatures with which they can manifest, we found a vast spectrum of composite particles F that has not yet been explored at the LHC. Recasting the recent CMS results of the resonant channel pp→e+F→e+e-qq¯′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$pp\rightarrow e^+F \rightarrow e^+e^- q\bar{q}'$$\end{document} (as F=E\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F=E$$\end{document}), we find that the composite fermion mass mF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_F$$\end{document} below 4.25 TeV is excluded for Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}/mF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_F$$\end{document} = 1. We further highlight the region of parameter space where this specific composite particle F can appear using 3 ab-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {ab}^{-1}$$\end{document}, expected by the High-Luminosity LHC, computing 3 and 5 σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} contour plots of its statistical significance.


Introduction
The parity-violating gauge symmetries and spontaneous/ explicit breaking of these symmetries for the hierarchy pattern of fermion masses have been at the center of a conceptual elaboration that has played a major role in donating to mankind the beauty of the Standard Model (SM) and possible scenarios beyond SM for fundamental particle physics. A simple description is provided on the one hand by the composite Higgs-boson model or the Nambu-Jona-Lasinio (NJL) model [1] with effective four-fermion operators, and on the other by the phenomenological model of the elementary Higgs boson [2][3][4][5][6][7]. These two models are effectively equivalent for the SM at low energies. After a great experimental effort for many years, using pp collision data at √ s = 7, 8 TeV at the Large Hadron Collider (LHC), the ATLAS [8] and CMS [9] collaborations have shown the first observations of a 125 GeV scalar particle in the search for the SM Higgs boson [10,11]. This far-reaching result begins to shed light on this most elusive and fascinating arena.
Recently, in the Run 2 of the upgraded LHC, studies on √ s = 13 TeV pp collision data are performed by ATLAS and CMS to search for new (beyond the SM) resonant and/or non-resonant phenomena [12][13][14][15][16]. These studies are continuously pushing up exclusion bounds on the parameter spaces of many possible scenarios beyond SM [17][18][19][20]. Among these models, are of particular interest composite-fermion scenarios that have offered a possible solution to the hierarchy pattern of fermion masses [21][22][23][24]. In this context [25][26][27][28][29], the assumption is that SM quarks "q" and leptons " " are assumed to be bound states of some not yet observed fundamental constituents generically referred to as preons and to have an internal substructure and heavy excited states F of masses m * F that should manifest themselves at the high energy compositeness scale . Exchanging preons and/or binding quanta of unknown interactions between them results in effective contact interactions of SM fermions and heavy excited states. While different heavy excited states have been considered in literature [30][31][32], below, we take as a reference the case of a heavy composite Majorana neutrino, N , for which the interaction Lagrangian would be (g * / ) 2q L γ μ q LN γ μ L . Its theoretical studies and numerical analysis have been carefully elaborated in [33,34]. Moreover, an experimental analysis of √ s = 13 TeV pp collisions at the LHC of the process pp → N → qq of the dilepton (dielectrons or dimuons) plus diquark final states has been carried out by the CMS collaboration [35] excluding the existence of N for masses up to 4.60 (4.70) TeV at 95% confidence level, assuming m N = .
In this article, we present phenomenological studies of new composite states according to a scenario recently proposed in Refs. [43,44] that relies on the four-fermion operators (interactions) of the NJL type and has escaped the spotlight of the LHC searches so far. The four-fermion interactions beyond SM considered in this new model are motivated by the theoretical inconsistency [36][37][38] between the SM bilinear Lagrangian of chiral gauged fermions and the natural UV regularization of unknown dynamics or quantum gravity, that implies quadrilinear four-fermion interactions (operators) of the NJL type, or Einstein-Cartan type [39][40][41][42], at high energies. On the basis of SM gauge symmetries, four-fermion operators of SM left-and right-handed fermions (ψ L , ψ R ) in the charge sector "Q i " and flavor family " f " can be written as From the point of view of an effective theory, these effective operators are attributed to the new physics at the high energy cutoff cut . The effective coupling G (1) has two fixed points: (i) the weak-coupling infrared (IR) fixed point and (ii) the strongcoupling ultraviolet (UV) fixed point. In the scaling domain of IR fixed point of the weak four-fermion coupling G at the electroweak scale v ≈ 239.5 GeV, effective operators (1) give rise to SM physics with tightly composite Higgs particle via the NJL mechanism, and also offers possible solution to the hierarchy pattern of fermion masses [43,[45][46][47][48][49][50][51][52][53][54][55][56]. The heaviest top quark mass is generated by the spontaneous breaking of SM gauge symmetries in the top sector (in Eq. 1) [59] with a tt bound state as a candidate for the SM Higgs particle, and three Goldstone bosons becoming the longitudinal modes of massive gauge bosons W ± and Z 0 . The reason why only the top sector undergoes the condensation is due to the energetically favorable ground state of NJL interac-tions (Eq. 1), as shown in Refs. [46,60]. Other SM fermion masses are generated by the explicit breaking of SM gauge symmetries due to the CKM flavor mixing of three SM generations [43]. Most importantly, the measured Higgs mass makes it possible to uniquely determine the solutions of the renomalization group equations for the form factor and quartic interaction of the composite Higgs particle [69][70][71]. The extrapolation of these solutions to TeV regime implies that the composite Higgs is a tighly bound state and a strongcoupling dynamics occurs. In the scaling domain of UV fixed point of the strong four-fermion coupling G at the composite scale ∼ O (TeV), composite fermions (bosons) form as bound states of three (two) SM elementary fermions and they couple to their constituents via effective contact interactions [44,[61][62][63][64][65][66][67][68].
We focus on the composite particles arising from fourfermion operators of NJL type, with massive (m F ) com-  [44,[69][70][71]. The effective coupling between the composite fermion (boson) and its constituents is given by the following contact interaction, which describes composite particle F f ( f ) production and decay: where (g * / ) 2 is a phenomenological parameter, and one can choose g 2 * = 4π so that 4π/ 2 is a geometric crosssection in the order of magnitude of inelastic processes forming composite fermions (Fig. 1). Whereas, (F / ) 2 is the Yukawa coupling between composite boson (Fig. 2) and two fermionic constituents, and (g * /F ) 2 relates to the form factor of composite boson. The composite fermion is in fact a bound state of an SM fermion and composite boson, namely The composite scales and F can only be experimentally determined like the electroweak scale v. The composite-fermion (-boson) mass m F , m ∝ and the proportionality is of the order of unity.
Analogously to composite-fermion scenarios mentioned above where new particles originate from preons (for more details see Refs. [21][22][23][24][25][26][27][28][29][30][31][32]) the present scenario in the domain of UV fixed point has two model-independent properties that are experimentally relevant: (a) the existence of composite fermions; (b) the existence of contact interactions, in addition to SM gauge interactions, which represents an effective approach for describing the effects of the unknown internal dynamics of compositeness. However, the present scenario is not only conceptually, but also consequently and quantitatively rather different from the previous composite- Fig. 1 A lepton , two quarks q (u-type) andq (d-type) form a composite fermion F via the contact interaction (dark blob) P L ,R (g 2 * / 2 ), where P L ,R = (1 ∓ γ 5 )/2. The thin solid line represents an SM elementary fermion, and the thick double line represents a composite fermion F. By a crossing symmetry applied to the lepton line → † (dashed line) the same diagram describes a 2 → 2 production process qq → † F Fig. 2 We show the Feynman diagrammatic representation for the contact interaction between the composite boson and its constituent quarks, where the thin solid line represents an SM elementary fermion, double wave line represents a composite boson, and the blob represents an interacting vertex (F / ) 2 P L ,R fermion scenarios. In fact, the composite fermions are formed as bound states of SM fermions, not preons, by strong fourfermion interactions of SM fermions at high energies and they further have different contact interacting processes. Therefore, it deserves more detailed phenomenological studies to reveal new features of the present scenario that are relevant to LHC experiments. This is the aim of this article and we find that the model foresee a large number of new composite particles that could appear in signatures not yet investigated and hence of great interest for the ongoing LHC physics program related to searches of physics beyond SM.
The model parameters in Eqs. (2) and (3) are unique for all SM fermions f and composite fermions F and bosons together with their interacting channels and we aim to study them for detailing the complete phenomenology of F, for all the corresponding flavors " f ". In Sect. 2 composite fermions' constituents and effective contact interactions among them are discussed considering the model in Eq. (1) with contact interactions of Eqs. (2) and (3). The production cross sections and decay widths of these composite fermions are calculated in Sect. 3, while in Sect. 4 the search for F, for all its flavors, is outlined deriving the final states and their topology that are relevant for its discovery at the LHC. It turns out that there is a wide range of new physical states that deserve dedicated searches at the LHC in order to investigate the entire phase space in which F can manifest. In Sect. 5, we take advantage of the aforementioned heavy composite Majorana neutrino N experimental studies [33,34] in the channel pp → N e − → e − e − qq to determine some constraints on the model parameters. We further compute 5σ contour plots of the statistical significance and highlight the region of parameter space where F can appear in the same channel using 3 ab −1 , as an example of the sensitivity to this model for a particular flavor of F. Finally, we summarize the work with some closing remarks in Sect. 6.

Four-fermion operators and contact interactions
In this section we describe the four-fermion operators and contact interactions that are relevant for the study of the phenomenology of the composite fermions at pp or ep collisions, including the LHC, which will be detailed in Sect. 3.

Composite fermions F
We consider, among four-fermion operators (1), the following SM gauge-symmetric and fermion-number conserving four-fermion operators, being the SM doublet i L = (ν e L , e L ) and singlet e R with an additional right-handed neutrino ν e R for leptons; ψ Lia = (u La , d La ) and u a R , d a R for quarks, where the color a, b and SU L (2)-isospin i indexes are summed over. Equation (4) or (5) is for the first family only, as a representative of the three fermion families. The SM left-and right-handed fermions are mass eigenstates, their masses are negligible in TeV-energy regime and small mixing among three families encoded in G is also neglected [43].
In Eqs. (4) and (5), each four-fermion operator has the two possibilities to form composite fermions, listed in Tables 1  and 2. Up to a form factor, E (N ) or D (U ) indicate a composite fermion made of an electron (a neutrino) or a down quark d (an up quark u) plus a color-singlet quark pair , and its superscript for electric charge. In Eq. (4), there are four independent composite fields F:   Table 3 Composite fermions F R , their constituents and SM quantum numbers Ra ) † γ 0 . They carry SM quantum numbers t i 3L , Y , and Q i = Y + t i 3L , which are the sum of SM quantum numbers (t i 3L , Y, Q i ) of their constituents, i.e., the elementary leptons and quarks in the same SM family [44], listed in Tables 3 and 4, so that the contact interactions in Eq. (2) are SM gauge symmetric.
The contact interactions for the production and decay of a composite fermions F are: In the case of Eq. (4) and Table 1, and their Hermitian conjugates, In the case of Eq. (5) and Table 2, and their Hermitian conjugates, These are relevant contact interactions for phenomenological studies of possible inelastic channels of composite-fermion production and decay in pp or ep collisions.

Composite bosons 0,±
From the four-fermion interaction in Eqs. (4) or (5), it is possible to form composite bosons and their Hermitian conjugates. Such normalized composite boson field has the same dimension [energy] of elemen-tary boson field. The composite boson carries the quantum numbers that are the sum over SM quantum numbers of its two constituents, see Table 5. These are pseudo composite bosons 0,± , analogous to charged and neutral pions π 0,± in the low-energy QCD. As shown in Fig. 2, the effective coupling between composite boson and its two constituents can be written as an effective contact interaction, where g Y = (F / ) 2 . By the appropriate normalizing the composite boson with the form factor (g * /F ) 2 in Eqs. (23)(24)(25), the effective contact interaction in Eqs. (26)(27)(28) can be expressed as a dimensionless Yukawa coupling g Y , whose value, corresponding to F value, can be different for composite bosons in Eqs. (23)(24)(25), but we do not consider such difference here.

Contact interaction of composite fermion and boson
In the view of the composite fermion being a bound state of a composite boson and an SM fermion, using compositeboson fields in Eqs. (23)(24)(25), we rewrite V † in Eqs. (11)(12)(13)(14) as follows: and their Hermitian conjugates V in Eqs. (7)(8)(9)(10), as shown in Fig. 3. In the same way, we rewrite V † in Eqs. (19-22) as follows: R can decay into composite bosons ± and 0 , which decay then to SM fermions, following the contact interactions in Eqs. (26)(27)(28) at the leading order of tree level. However, we shall consider other decay channels at the next-to-leading order, such as neutral composite boson decay to two SM gauge bosons 0 u,d →G +G .

Contact interaction of 0 composite boson and gauge bosons
Analogously to π 0 → γ γ , the massive 0 u,d composite boson can also decay into two gauge bosons [44] : via the contact interaction where g and g represent the couplings of gauge bosons A μ and A ν to the SM quarks u and d with different SU L (2)isospin i = u, d. Actually, this effective contact interaction (41) is an axial anomaly vertex, as a result of a triangle quark loop and standard renormalization procedure in SM. It should be mentioned that two gauge bosons can be two gluons that possibly fuse to Higgs and other particles in the final states.

Phenomenology of the composite fermions in pp collisions
In this section we study the phenomenology of the composite fermions in pp collisions. We first outline their production and decay mode and then calculate their cross section and decay width. This study leads us to the discussion on the search for F that will be discussed in the next section.

Production and decay of F
As already specified in Sect. 2, the composite fermion F can be where E, N , D, U stand for the charge sector Q = −1, 0, −1/3, 2/3 respectively, and the corresponding composite fermions of the higher generation of families in the SM. The kinematics of the processes is derived in the center of mass frame of pp collisions and virtual processes of F are not considered. If the energy √ s in the parton center of mass frame is larger than composite fermion masses, the following resonant process can occur: where the Standard Model fermion f is produced in association with the corresponding composite fermion F. We note that f in Eq. (43) can be e, ν, u, d, and the corresponding Standard Model fermions of the higher generation of families. The kinematics of final states is simple in the center of mass frame of pp collisions. If we neglect the quark-family mixing, the previous process can manifest at parton level as: The composite fermion F can decay through two different channels:f plus two quarks, via the interactions in Eqs. (11)(12)(13)(14) and Eqs. (19-22); orf plus a composite boson , via the interactions in Eqs. (29)(30)(31)(32) and Eqs. (33)(34)(35)(36), wherē f indicates a fermion that is the antiparticle of f . Then the composite fermion F decays as:  52) and (53) and the filled circles (black) represent the results from our implementation of the model in CalcHEP. We find good agreement. On the right panel, we plot the decay width of composite fermion F as a function of its mass m F for the case m F / = 0.8. Again, we observe a good agreement between the expectation from a CalcHEP simulation and the analytical result based on Eq. (54) The full decay chain is: It appears clear, considering all the possible flavors of f and F, that a large range of final states is possible. The cross section of the process pp → f F, the decay branching ratio of F, and the final states and their topologies are discussed below.
We consider the production cross section for the composite fermions F in pp collisions expected at the LHC collider according to Feynman's parton model. The QCD factorization theorem allows to obtain any hadronic cross section (e.g. in pp collisions) in terms of a convolution of the hard partonic cross sectionsσ , evaluated at the parton center of mass energy √ŝ = √ τ s, with the universal parton distribution functions f a (x,Q) which depend on the parton longitudinal momentum fractions x, and on the factorization scaleQ: The factorization and renormalization scale Q is generally fixed at the value of the mass that is being produced. The parametrization of the parton distribution function is NNPDF3.0 [72] and the factorization scale has been chosen asQ = m F . The right panel of Fig. 4 shows the agreement between analytical calculations based on Eqs. (52) and (53) for the case of the composite fermion F, and the results of simulations with CalcHEP where the model with four-fermion interactions has been implemented. We remark the quite good agreement that validates our model implementation in CalcHEP.

Decay widths
Analytical calculations, in the similar way as the first term in Eq. (5) of Ref. [33], yield to the width of the 3-body process F →f qq Note that at TeV energy scales, composite fermions are massive (m F ) Dirac fermions, whereas all SM elementary fermions are treated as massless Dirac fermions of four spinor components, consisting of right-and left-handed Weyl fermions of two spinor components. Alternatively, the decay width F has also been evaluated via CalcHEP, and numerical results are completely in agreement with analytical one in Eq. (54), see the left panel of Fig. 4. The decay width of the composite fermion F in the process F →f can easily be computed from the effective contact lagrangian in Eqs. (29)(30)(31)(32)(33)(34)(35)(36) and the total width is The decay width of the boson to two quarks is simply calculated by using the effective contact Lagrangian in Eq. (26) and (27), For the case that equals to + or − composite boson, → qq of Eq. (57) is the only decay channel, see Eqs. (23) and (26). The 0 u,d composite bosons, instead, can also decay to two gauge bosonsGG (37)(38)(39)(40), according to the contact interaction (41), the corresponding decay widths are [44] : where θ W is the Weinberg angle, and the number of colors N c = 3. Total decay rate tot ( 0 u,d →GG ) is the sum over all contributions from Eqs. (58)(59)(60)(61). The total 0 u,d -decay rate reads where ( 0 d → qq ) is given by Eq. (57). Based on these results, we calculate the exact branching ratios of different channels for different parameters of the model.

Branching ratios
The branching ratios of the 0 u,d decay to two quarks qq , and the 0 u,d decay to two gauge bosonsGG , Whereas, the branching ratios of the composite fermion F decay to f , The branching ratios of the direct decay F →f qq , and indirect decay F →f →f qq , The sum of these two branching ratios gives the total branching ratio B(F →f qq ) of F decaying tof qq , For the case ± , B( ± → qq ) = 1. For the case 0 u,d , is given by Eq. (64), and the branching ratios of decay F →f 0 u,d →fGG is given by As a result, the cross sections of these processes are: and All channels of composite fermion F production and decay give the same results at this level of approximation by using contact interactions only. For the processes with the e + e − qq final state in pp collisions, the total cross section is approximately given by and the total width is      The calculation of these quantities will be given in the next sections.

Search for F at the LHC
After having discussed the production and decay of F and its cross-section, width, and decay branching ratio, we now examine these results in terms of parameters of the model and derive the possible final states and their topologies, highlighting their impact to the current program of beyond SM searches at the LHC.

Branching ratios and topology of F with respect to model parameters
In order to present the branching ratios of different possible channels in terms of parameters of the model, we are bound to discuss physically sensible parameters to explore. This model has four parameters that can be rearranged to three dimensionless parameters for a given value: The ratio m F / < 1 (m / < 1) of the composite fermion (boson) mass and the basic composite scale gives us an insight into the dynamics of composite fermion (boson) formation. In addition, as the parameters m and F represent the same dynamics of composite boson formation we use the F /m ratio. Finally, to take into account the feature that a composite fermion F is composed by a composite boson and an elementary SM fermion, we adopt the ratio m /m F < 1 as a parameter. As a result, for given √ s and values, we have three parameters m F / , F /m , and m /m F to represent the results of the possible branching ratios. Figure 5 shows three sets of plots for the branching ratios of F →f qq and F →f (left column), → qq and → GG (center column), and the large-radius jet: full decay chains F →f qq , F →f →f qq , and F →f →f GG (right column). These branching ratios are plotted with respect to the m /m F ratio and for the values of F / m = 0.8, F / m = 5, F / m = 15. Note that the branching ratios of F → f qq and F → f (left column) do not depend on m F / , which can be seen from Eqs. (54)(55)(56). For the ± decay, the channel (57) is unique, so the branching ratio is one, independent of m F / . Whereas the 0 u,d decays also to G G, see Eqs. (64) and (65), and the branching ratio depends on m F / . However, in the regime of m /m F we consider, the branching ratio of 0 u,d → G G is very small and negligible, compared with the branching ratio of 0 u,d → q q. Therefore, our results of branching ratio of F decay presented in Fig. 5 are independent of the parameter m F / . As a result, regarding the branching ratio and topologies of F decay, the model effectively depends only on two parameters F /m and m /m F . However, the cross section of F production depends on the parameter m F / , as indicated in Eq. 52.
In Fig. 5, it is shown that the branching ratios F → f qq direct (67) and F → f → f qq indirect (68) tend to swap each other for different values of F /m . Increasing the coupling F (and the ratio F /m ) of composite boson to its two constituents (qq), see contact interactions (26,27,28), the branching ratios F → f → f qq indirect (68) becomes dominant over the branching ratios F → f qq direct (67). We thus consider two reference cases of F /m : We further outline that the value 0.2 for m /m F that separates the boosted from the resolved topology when F decays through a is indicative and may vary based on the mass of F. We have verified, using CalcHEP, that for a mass of F of 1 TeV, the two quarks decaying from the originating from F are indeed within a R(q,q ) below 0.8, which determines the cone size of a large-radius jets suitable for boosted jet at the LHC experiments. For a mass of F of 7 TeV instead, the value 0.2 is lowered to 0.15 to guarantee R(q,q ) < 0.8, as for higher masses of F its width increases.
So far we have presented the studies relying on the first SM generation, i.e. considering the electron flavor of F (F = E). However, these discussions and calculations are straightforwardly generalized to the second and third SM generations, and in the text we keep the general notation F instead of E. At the leading order of only contact interactions being taken into account, the formulae of cross-sections, decay rates and branching ratios of the second and third generations are the same as those of the first generation, however composite fermions F and their productions, decay channels and rates are different, depending on the values of their parameters m F , m and F , which vary from one SM generation to another. The reasons are that the effective couplings ∼ g 2 * / 2 and ∼ g Y of contact interactions can be different, due to the modifications from the flavor mixing and gauge interactions. However, we neglect these modifications in this article, and expect the variations of masses m F , m , decay constant F to be small. In this sense, there are two basic parameters F /m and m /m F for each SM generation to Table 7 List of all the signatures foreseen by the model for the different flavors of f , F, their topology and final states. In the table " j" indicates any reconstructed jet that originates from a u, d, s quark and "J " a large-radius (cone-size of 0.8) jet that is reconstructed from 2 quarks produced with low angular separation. The notation w/ indicates the decay of F through a . The LHC searches for these signatures are reported along with some comments on the features that are typical of this model and have not been exploited in the referenced LHC searches. In the table, "n/a" indicates the case in which an LHC search interesting the corresponding final state has not been found  [87,88] N identification Resolved w/ → qq ν(ν(qq )) [87,88] N , identification Boosted νν J [92] 2-prong, no V boson tag, boosted deacy Here, we consider one single channel (ν N ) for all the ν neutrinos of the Standard Model and one single channel ( j J ) for the u, d, s quarks. We notice that the flavors cC is taken separately, because of the improving performances in c-tagging algorithms at the LHC [93,94] and dedicated searches for new physics with c quarks in the final state [95,96]. Moreover, we distinguish the three topologies (resolved with and without and boosted) explained in the previous paragraph, so that we have in total 24 different signatures that have to be considered in order to pursue a comprehensive search of F. In Table 7 we outline these signatures based on the flavors f and F, the possible topology, the corresponding final states, and the LHC search that, to the best of our knowledge, could be more sensitive to searching for F. In the last column, we further report on features of F and its decay that have not been exploited directly in the cited LHC searches and could be used to improve possible future searches dedicated to F. We especially remark that the F quark flavors appear to be completely unexplored yet and we urge on the importance of carrying out specific analyses at LHC to investigate it. This is certainly noteworthy and can have a relevant impact on the beyond the SM physics program of the LHC. We notice that, despite possible and with a peculiar signature, the channel f F → f (f ) → f (f (GG )) is negligible since the decay 0 → (GG ) is only relevant for F / m = 0.1 when, however, the decay F → f ( ) has a branching ratio close to zero. Because of this we do not include this case in Table 7 and we point out that this case would become of interest in the case F is found in one of the possible signatures mentioned above, to study the nature of the new particle.
We put emphasis on the case of the composite fermions F = N 0 ,N 0 , N + , N − for the final state ννqq , where νν stands for the pairs of the SM left-handed neutrino ν e L and/or sterile right-handed neutrino ν e R , as the latter is a candidate of dark-matter particles.
Finally, we acknowledge that the final state ff qq is relevant for a wide range of the parameter space of the model and thus will consider it in the next section and, in particular, the case of f = e and F = E to derive limits on the model parameters based on existing LHC results.

Bounds on the model parameters
In this section we provide a discussion of the bounds of the model parameters taking as a reference the case F = E and the final state eeqq that has been shown to be sensitive to a wide portion of the parameter space of the model in the previous section. For this purpose, we recast the 95% confidence level (C.L.) experimental upper limit on σ ( pp → eeqq ) using a recent analysis [35]  We also performed a study about the potential of a dedicated analysis in the High Luminosity LHC (HL-LHC) conditions (center of mass energy of 14 TeV and luminosity of 3 ab −1 ). We used CalcHEP to generate the processes and DELPHES [73] to simulate the detector effects. In order to separate the signal from the background, we selected events with p t e 1 ≥ 180 GeV, p t e 2 ≥ 80 GeV, p t j 1 ≥ 210 GeV, m ee ≥ 300 GeV ( p t is the transverse momentum, e 1 the leading electron, e 2 the subleading electron, j 1 the leading jet and m ee the invariant mass of the two electrons). Then we evaluated the reconstruction and selection efficiencies for signal ( s ) and background ( b ) as the ratio of the selected and the total generated events. From these efficiencies, the signal and background cross sections (σ s , σ b ) and the integrated luminosity (L), it is possible to evaluate the expected number of events for the signal (N s ) and the SM background (N b ) and finally the statistical significance (S): The S = 5 contour curve is shown by the upper (solid) line in Fig. 7. It can be used to get indications about the potential for discovery or exclusion with the experiments at the HL-LHC, showing that there is a wide region of the model phase space where the existence of the composite fermions can be investigated; for the case m F / = 1 we can reach masses up to ≈ 6.2 TeV. We notice that, despite having considered the case of F = E in this section, it could be inferred that the cross section for F, which should approximately be the same for all its flavors, is sufficient for the F to appear at the LHC with the statistics already collected at the LHC experiments, and that is expected by the HL-LHC. Based on this result, we recommend that the physics program of the LHC consider the new particles foreseen by this model and their signatures in its investigations.

Summary and remarks
In the weak coupling regime the effective four-fermion operators of NJL-type possess an IR-fixed point, rendering the elegant Higgs mechanism of the SM of particle physics at low energies. In the strong coupling regime, on the other end, these operators could possess an UV-fixed point, giving rise to composite fermions F (bosons ) composed by SM fermions as bound states of three (two) SM elementary leptons or quarks, and to their relevant contact interactions with them at high energies O(TeV). We study, for the first time for this model, the spectrum of composite particles and contact interactions in quark-lepton and quark-quark sectors in relationship to their phenomenology at the LHC in order to unveil their discovery potential. The cross sections and decay rates of composite particles are calculated based on the LHC physics from pp collision at high energy TeV scale. We find out that a comprehensive investigation of the model presented here can be effectively achieved, for given √ s and values, by considering only two parameters: F / m and m /m F . Based on these results, we exhaustively examine all the possible F states and the signatures with which they can manifest at the LHC, according to different F / m and m /m F . Interestingly, we find that there is a broad variety of new composite particles that could manifest in signatures that have escaped the realm of the searches at the LHC. We summarize these cases in Table  7. They can offer an unprecedented discovery potential of physics beyond the SM and we urge on the importance for the LHC experiments to include such searches in their ongoing physics program.
In order to set bounds on the model parameters, we derive constraints for the particular case where F has electron-like flavor. We analyzed the particular processes giving e + e − qq final state by using the recast of the experimental upper limit by the CMS collaboration on the cross-section σ ( pp → eeqq ). We determine that a composite fermion F of mass m F below 4.25 TeV can be excluded for = m F . At the same time, we compute 3σ and 5σ contour plots of the sta-tistical significance and highlight the phase space in which F can manifest using 3ab −1 , foreseen at the high luminosity LHC (HL-LHC). This result shows that, even for final state traditionally considered at the LHC experiments, there is a vast range of model parameters to which a dedicated search can be sensitive to the F composite fermions. We encourage such efforts in future investigations in light of peculiar features of F not yet exploited at the LHC searches and highlighted in Table 7. We are preparing the next article presenting the investigation of phenomenology at the LHC of composite bosons.
Both composite bosons ( ) and fermions (F) are mass eigenstates and have definite SM quantum numbers, so that the CKM flavor mixing and the Feynman diagrammatic representations of SM perturbative gauge interactions can be easily implemented, see Eqs. (4.8)-(4.11) in Ref. [44]. At the leading order of contact interactions considered in this article, the contributions from all CKM flavor mixing and gauge interactions to composite particle masses are neglected. It should be mentioned that these effective couplings between composite particles and SM gauge bosons have two main effects. First, they give the possible decay channels of composite particles to final states involving SM gauge bosons, like two gauge bosons as in Sect. 2.4. The branching ratios of these channels are negligible in the composite fermion production and decay studied in this article. However, they could be important in the composite boson formation and decay, which is under examination in a separate effort. Second, these effective couplings between composite particles and SM gauge bosons should affect the well-measured SM quantities in the IR regime, like the electroweak oblique parameters or the decay width of the Z boson, as well as the Higgs physics. These issues are indeed important and necessary, and will be addressed in a future study. So far only some general and qualitative discussions have been presented [43,44], showing that the composite Higgs boson is a tight bound state of tt, as if it was an elementary particle, and possible corrections to the SM quantities are small.
It is an interesting question to see how these phenomenologies can possibly account for some recent results obtained in both space and underground laboratories. The cosmic rays pp collisions might produce composite fermions F = E that decay into electrons and positrons. This may explain an excess of cosmic ray electrons and positrons around TeV scale [74][75][76][77]. In addition, recent AMS-02 results [78] show that at TeV scale the energy-dependent proton flux changes its power-law index. This implies that there would be "excess" TeV protons whose origin could be also explained by the resonance of composite fermions F = N due to the interactions of dark-matter and normal-matter particles. These composite fermions should appear as resonances by high-energy sterile neutrinos inelastic collisions with nucleons (xenon) at the largest cross-section, then resonances decay and produce some other detectable SM particles in underground laboratories [79]. Similarly, in the ICECUBE experiment [80], we expect events where the neutrinos change their directions (lower their energies) by their inelastic collisions to form the resonances of composite fermions N at a high energy scale (≈ TeV). It is worthwhile to mention that in the IR domain of this model there are effective coupling vertexes of the SM gauge boson W and the right-handed currents: g RνR γ μ e R W + μ or g Rū R γ μ d R W + μ , where g R ∝ (v/ ) 2 [81,82]. The parity symmetry is restored at the scale [43,44] and none of additional intermediate gauge bosons W R or W is present. The recent phenomenological studies of this effective coupling in the quark sector can be found in Ref. [83]. In the lepton sector, these effective contact interactions relate to the darkmatter physics of right-handed neutrinos ν R . Similarly to the analogy between the Higgs mechanism and BCS superconductivity, the composite-particle counterparts in condensed matter physics have been recently discussed [84].  Fig. 8 The branching ratios of all the possible decays of F for values of F /m between 1 and 15, increasing in step of 1, to complement the cases reported in Fig. 5 (right column)