Testing the origin of the “f1(1420)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_1(1420)$$\end{document}” with the K¯p→Λ(Σ)KK¯π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{K}} p \rightarrow \Lambda (\Sigma ) K{\bar{K}} \pi $$\end{document} reaction

We study the K¯p→YKK¯π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{K}} p \rightarrow Y K{\bar{K}} \pi $$\end{document} reactions with K¯=K¯0,K-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{K}} = {\bar{K}}^0, K^-$$\end{document} and Y=Σ0,Σ+,Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Y=\Sigma ^0, \Sigma ^+, \Lambda $$\end{document}, in the region of KK¯π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K{\bar{K}} \pi $$\end{document} invariant masses of 1200-1550\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1200-1550$$\end{document} MeV. The strong coupling of the f1(1285)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_1(1285)$$\end{document} resonance to K∗K¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^* {\bar{K}}$$\end{document} makes the mechanism based on K∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^*$$\end{document} exchange very efficient to produce this resonance observed in the KK¯π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K{\bar{K}} \pi $$\end{document} invariant mass distribution. In addition, in all these reactions one observes an associated peak at 1420 MeV which comes from the K∗K¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^* {\bar{K}}$$\end{document} decay mode of the f1(1285)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_1(1285)$$\end{document} when the K∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^*$$\end{document} is placed on shell at higher invariant masses. We call the attention to the possibility that the peaks observed in other reactions where the “f1(1420)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_1(1420)$$\end{document}” is observed have a similar origin.


Introduction
Kaon beams are called to play a relevant role in studies of hadron dynamics. So far kaon beams in the GeV range are available at the Japan Proton Accelerator Research Complex (J-PARC) [1,2]. Kaon beams of low energy from φ decay are available in the Phi-Factory at Frascati [3,4]. The COM-PASS collaboration at CERN can also select events coming from kaon induced reactions [5][6][7]. In addition the planned K L factory at Jefferson Lab has passed the first steps for approval [8,9]. The main aim of this latter factory is the study of strange hadron spectroscopy, but its potential for non-strange hadron studies is also relevant. In this sense, having this facility in mind, a proposal was made in Ref. [10] to produce the f 0 (980) and a 0 (980) resonances inK induced reactions on protons, which should shed light on the much debated nature of these resonances. a e-mail: liangwh@gxnu.edu.cn (corresponding author) b e-mail: oset@ific.uv.es In the present work we make a proposal that shows the usefulness of such facilities to provide relevant information on the nature of the " f 1 (1420)" resonance, which is catalogued in the PDG [11] as a standard resonance. This resonance has been observed in many very high energy reactions, as proton or pion induced, e + e − , γ γ and J/ψ decays, and exceptionally in one experiment with a K − induced reaction at 32.5 GeV/c in Serpukhov [12,13]. One of the peculiar properties of this resonance is that it is only seen in the KK π decay mode, supposed to be K * K + c.c., with a small fraction seen in the a 0 (980)π channel and a very small fraction to φγ . However, the right to be catalogued as a standard resonance was challenged in Ref. [14], where it was proved that in such reactions a peak appeared naturally at 1420 MeV as the K * K + c.c. decay mode of the f 1 (1285) resonance, which showed two peaks in the KK π decay channel, one at the nominal f 1 (1285) mass and another one around 1420 MeV, where the K * K channel with K * on shell becomes open. The small fraction of a 0 (980)π decay was also shown in Ref. [14] to correspond to a triangle singularity decay mode of the f 1 (1285).
It is interesting to mention how the π 0 a 0 (980) mode is generated in Ref. [14]. A triangle mechanism is considered where the f 1 (1285) decays to K * K , K * decays to π 0 K and then theK K fuse to give the a 0 (980) resonance. This mechanism was shown to have a triangle singularity peaking at M inv (π 0 a 0 ) 1412 MeV, and it explained the π 0 a 0 spectrum observed in the experiment of Ref. [15]. In the same experiment [15], a γ φ decay mode of the f 1 (1420) is also reported, albeit with small statistics since the branching ratio is around 0.3%, about an order of magnitude smaller than the branching ratio for π 0 a 0 (980) (4%). It is not our purpose to study this process here, but it is pertinent to mention that, similarly to the π 0 a 0 (980) decay mode, one can get a γ φ decay mode with a similar triangle mecha-nism, f 1 (1285) → K * K , followed by K * → K γ and then KK fusing to give the φ. The mechanism develops a triangle singularity, at M inv (γ φ) 1384 MeV (see Eq. (18) of Ref. [16]), which is converted into a broad bump because of the width of the K * . Provided one can get improved statistics in the future, such an idea would be worth pursuing, but for the moment we note the fact that such a decay mode is unavoidable with the dynamics of the f 1 (1285) coupling strongly to K * K .
In the present work we show that the peak at 1420 MeV can also be produced in the present and future lower energy Kaon Facilities with reactions that one has under control theoretically, such that their experimental search and comparison with the theoretical predictions can shed much valuable light on the origin and nature of the " f 1 (1420)" peak. For this purpose we propose several reactions, all tied among themselves, for which we make evaluations of the cross section, with the assumption that the f 1 (1420) is only a manifestation of the K * K + c.c. decay mode of the f 1 (1285). We make absolute predictions of the cross section for the three reactions, and the strength of this cross section around the f 1 (1285) and f 1 (1420) peaks results as a consequence of the coupling of the f 1 (1285) to K * K +c.c. and dynamical details of the production mechanisms in theK -induced reaction. We should note that a strong coupling of the f 1 (1285) to K * K + c.c. components emerges in pictures where this resonance is dynamically generated from the pseudoscalar-vector meson interaction [17][18][19][20]. However, independent of this information, the strong coupling of the f 1 (1285) to K * K + c.c. has experimental backing from the explicit experimentalK K π decay mode of the f 1 (1285) (see a theoretical description of this decay mode in Ref. [21]).

Formalism
The f 1 (1285) axial vector meson has the quantum numbers I G (J PC ) = 0 + (1 ++ ). Using chiral Lagrangians for pseudoscalar-vector interaction [22] and a chiral unitary approach in coupled channels, the axial vector mesons emerge as a consequence of the interaction [17][18][19][20]. In particular the f 1 (1285) is the cleanest example, appearing in a single channel K * K −c.c., with the mass at the right position using a cut off to regularize the loops that provides masses for the other resonances in agreement with experiment as an average. The masses of the other resonances are not provided as accurately as the one for the f 1 (1285) in that approach. On the other hand, the important partial decay widths of this resonance have been studied and are described successfully in this approach, the π 0 a 0 (980) and π 0 f 0 (980) decays in Ref. [23] and the KK π decay in Ref. [21].
The f 1 (1285) state is given, with the isospin convention (K + , K 0 ), (K 0 , K 0 ) and C(K * + ) = −K * − , etc., by The coupling of the state to this K * K − c.c. combination in s-wave obtained from the residues of the t f 1 , f 1 amplitude at the pole with z the complex energy, and z R the complex pole position, is given by We shall take the second value in our calculations. The t f 1 ,K * + K − amplitude is then given by with the polarization vector, similarly with the other components.
In view of this, the K − p → Y f 1 → Y K πK decay proceeds via the mechanism shown in Fig. 1.
We should note that this mechanism is the same one used in Ref. [24] to study the K − p → f 1 (1285) reaction, but in that work the process stops in the f 1 production of Fig. 1 without looking into the specific decay channel of K πK , hence, the f 1 (1420) excitation was not addressed in that work.
We shall specify the K 0 π + K − decay channel of f 1 , and, hence, only the (b) and (g) diagrams of Fig. 1 contribute.
Apart from the f 1 coupling to K * K , we need the coupling of K * to K π which is given by the standard Lagrangian with P, V the standard SU(3) pseudoscalar and vector meson matrices On the other hand, we need the K * B B vertex. This vertex is given in terms of the Lagrangian [25,26] with the SU(3) baryon matrix B given by It is worth noting that while Eq. (8) is not provided by chiral Lagrangians, the combination of this Lagrangian together with Eq. (5) and vector exchange between mesons and baryons gives rise exactly to the lowest order mesonbaryon chiral Lagrangians [27] when q 2 → 0 in the propagator of the exchanged vector [28,29].
Since we will work with relatively lowK energies, we will take the spatial components of the K * , i , neglecting the 0 component, which has been shown to be a very good approximation even for relatively large K * momenta [30]. The relevant couplings are then given by with , for pK * − → ; 1, for pK * 0 → + .
A different B BV vertex is used in Ref. [24] following Refs. [31,32], given by where M, M , p, p are the masses and momenta of the incoming and outgoing baryons, with α = 1.15 and κ = 2.77 according to Refs. [33,34] and a slightly different value of g. The magnetic σ μν term of Eq. (12) is usually neglected in chiral dynamics calculations since the momentum transfer is small, but in this case it is not negligible and we shall take it into account. The coupling of Eq. (12) seems much bigger than that of Eqs. (10), (11) but they are accompanied with form factors which are normalized differently. In the formalism of Eq. (12) the form factor is normalized to unity when q 2 = ( p − p) 2 = m 2 K * , while in Eqs. (10), (11) it is normalized to unity for q 2 = 0. We shall take a form factor of the type 2 2 − q 2 for Eqs. , (10), (11); and we shall see the difference between the results with the two approaches, which we will accept as uncertainties of our results.
With the ingredients discussed above, the amplitude for the diagram of Fig. 1b is given bỹ We should also sum coherently the contribution of the diagram of Fig. 1g. Since K * + andK * 0 are different particles, in the limit of small K * width these diagrams do not interfere. However, with the finite K * width there is an interference and we take it into account. As a consequence, and evaluating explicitly the ū|γ i |u , ū|σ iν |u matrix elements we obtain at the end, for angles close to the forward direction, where an angular average in [( p K − p π ) · p ] 2 and [( p K − p π )· p ] 2 has been done, and linear terms in cosines of angles are neglected in the interference term. The magnitude p π in the interference term of Eq. (15) is the pion momentum in the f 1 (1285) rest frame, given by The cross section for the reaction is given by where the angle θ is defined between K − and f 1 in the K − p rest frame. The limits of the integration for M inv (K 0 π + ) with fixed M inv (K − π + ) are given in the PDG [11], and pK is the initialK momentum.
Since the peak of the M inv (K 0 π + K − ) distribution around 1420 MeV is sensitive to the tail of the f 1 (1285) resonance, we do a more refined picture for the f 1 (1285) width than taking the nominal width as a constant. We take the important π 0 a 0 (980) decay channel explicitly, with a 38% branching ratio, and also the K * K decay width, which is relevant around 1420 MeV. Hence we take with , where S(s v ) is the K * spectral function with m K * and K * the nominal K * mass and width [11], and f 1 , K * K (s v ) the K * K decay width of the f 1 with K * mass √ s v given bỹ with .
In addition, we also take K * in Eq. (15) as energy dependent, with where with π K standing for K 0 π + or K − π + depending on the term in Eq. (15). We can tune the parameter to get an important datum which is the integrated cross section for f 1 (1285) production in the K − p → f 1 (1285) at √ s = 3010 MeV, σ = 11 ± 3 μb [35]. The process can be described by the diagram of Fig. 2.
Proceeding as before, we find now where p f 1 , pK are the f 1 andK momenta, and now (

Results
The first step is to use Eq. (23) to obtain the cross section for K − p → f 1 (1285). By using Eq. (12) for the B BV coupling and a value of = 1250 MeV, we find σ = 11.07 μb. This result is compatible with those of Ref. [24] where a bigger is demanded because extra powers of 2 −m 2 K * 2 −q 2 are used in the form factor. If we use Eq. (23) together with Eqs. (10), (11) and the form factor of Eq. (13) normalized to 1 at q 2 = 0, then we need a value = 1900 MeV and we get σ = 11.24 μb. We shall use these two versions of the B BV vertex and evaluate the cross section for K 0 π + K − production of Eq. (16), accepting the differences as uncertainties. The results will be given with the coupling of Eq. (12).
In Fig. 3 we show dσ dcosθ dM inv (K 0 π + K − ) for cosθ = 1 as a function of M inv (K 0 π + K − ) for theK 0 p → + K 0 π + K − at twoK 0 energies corresponding to √ s = 2850 MeV and √ s = 3010 MeV. We observe a clear peak at the f 1 (1285) nominal mass, corresponding to the K πK decay of the f 1 [11], but interestingly we find also a peak around M inv (K 0 π + K − ) of 1420 MeV. This peak comes from the f 1 (1285) → K * K when the K * becomes on shell. We should see this peak as a combination of two factors: The increasing of M inv (K 0 π + K − ) allows the intermediate K * to get on shell, increasing the cross section, but the tail of the f 1 (1285) tends to reduce the cross section with increasing M inv . The result of it is the peak seen at 1420 MeV which is hence the manifestation of the K * K decay mode of the f 1 (1285) and does not come from any new resonance. We can make an estimate of the width of the 1420 MeV peak by taking a background below the curve, as an experimental analysis would do. A rough estimate can be done drawing a straight line between the minimum at 1350 MeV and the last point of the distribution. We find a width There is no need to show results for the K − p → 0 K 0 π + K − reaction, since neglecting the 0 , + mass difference the cross section is just 1 2 of the former one (see Eq. (11)). At this point, it is interesting to see what happens if we use the B BV coupling of Eq. (11). The results are very similar but about a factor of two smaller. We accept this as uncertainty in our results. Yet, the important thing is that the shape of the cross section is practically identical, with about the same ratio of the strength at the 1420 MeV and 1285 MeV peaks.
In Fig. 4 we show the cross section for the K − p → K 0 π + K − reaction for cosθ = 1 and √ s = 2850 MeV, 3010 MeV.
The results are very similar to those shown before except that they are about a factor of 6 smaller in size, according to Eq. (11). Once again, the peak at 1420 MeV appears as before and the ratio of the strengths at the peaks is also very similar to that found before. We have also evaluated the results with the B BV input of Eqs. (10), (11) and the form factor normalized to 1 at q 2 = 0 and, as before, we find the same shape for the cross section and a size about 1 2 the former one. The size of the cross sections obtained is relatively large, and the fact that there are results in Ref. [35] guarantees its measurability in present facilities.
It is worth making some observation about possible contamination of background from other sources. In order to minimize this potential contribution, one should bear in mind that the characteristics of the peak at 1420 MeV that we obtain is that it corresponds to K * K , concretely K * + K − andK * 0 K 0 since, as discussed above, the peak comes from  Fig. 1b, g. The other important feature is that the K * K are produced in relative s-wave. These two conditions can serve as a fitter of other possible background contributions that could distort the picture that we have obtained.
When arriving to this point it is worth making some observations. In all cases we observe that the peak of the f 1 (1285) is always bigger than for the 1420 MeV peak. This is not a universal feature independent of the reaction. Indeed, the main ingredients are the width of the f 1 (1285) and the phase space for K πK production mediated by a K * . However, in dσ/[d cosθ dM inv ( f 1 )], for fixed cosθ , as we have done, thē K * propagator depends on M inv ( f 1 ) and becomes weaker as M inv ( f 1 ) increases. We have then taken a different magnitude dσ/[dt dM inv ( f 1 )] (dt = 2 p K − ,ini p d cosθ ) and have fixed t at the value of the maximum M inv in Fig. 3 for √ s = 3010 MeV. The results are shown in Fig. 5.
What we see in Fig. 5 is that while the ratio of strengths in the peaks in dσ/[d cosθ dM inv ( f 1 )] is about 3.7 in Fig. 3, this ratio is 2.6 in Fig. 5. We simply want to draw the lesson that the relative strength of the peaks depends on the reaction and the quantity chosen to see the peaks. In this sense, the investigation of different reactions, and the present work shows one such example, is most useful to settle the issue of the meaning of the 1420 MeV peak seen in different reactions.
In this context it is worth making a survey of the different works from where the f 1 (1420) resonance has been claimed. The different works look at the mass distribution of the KK π channel to see the f 1 (1420) peak. Many of the old experiments have poor statistics and some of them do not see a signal for the f 1 (1285). This is not the case in more modern experiments. In general the experiments based in π p scattering at very high energy tend to see similar strength for the f 1 (1285) and f 1 (1420) peaks [36][37][38]. However, there are exceptions. In Ref. [39] the signal for f 1 (1285) is bigger than for f 1 (1420), and with some cuts the f 1 (1420) sig-nal becomes clearer and a bit magnified. In Ref. [40,41] the strength of the f 1 (1285) is significantly bigger than the one of the f 1 (1420), and in Ref. [42] the f 1 (1285) signal is also magnified with some cuts. In pp induced reactions the situation is similar, many of them show a similar strength for the two peaks [43][44][45][46][47]. But, again, there are exceptions. In Ref. [48] the signal for the f 1 (1420) is a bit larger than for the f 1 (1285). In Ref. [49] the signal for the f 1 (1285) is bigger than for f 1 (1420) for a range of the chosen longitudinal momentum, while for another range the strengths are similar. There are also data in pp annihilation. In Ref. [50], producing KK πππ, a similar strength for the two peaks is seen for some KK π spectra, but in the K 0 S K 0 S π 0 channel the f 1 (1285) signal is clearly seen but the f 1 (1420) disappears if a cut of π + π − π 0 is made in the ω region. In γ γ collisions in Ref. [51] the strength of the f 1 (1285) is about one third from the one of the f 1 (1420), but one must note that the f 1 (1285) falls at the beginning of the acceptance region. Another example is the Z decay studied in Ref. [52]. The bare data do not show any signal of either resonance, but when a cut is made in the KK invariant mass the two peaks appear with about the same strength, a bit larger for the f 1 (1420).
It is not clear what to conclude about the shape and strength of the two peaks, except that they depend much on the reaction and cuts made in certain observables, which means that this issue is very much tied to the dynamics of each process and each of them would deserve a separate study.
In our reaction we found the strength of the f 1 (1285) larger than the one of the f 1 (1420). We have tried to see if this is due to uncertainties of the model and for this we have investigated two sources of possible modifications. The first one is to allow the two kaons in Fig. 1 rescatter, giving rise to a triangle diagram that according to Ref. [16] develops a triangle singularity at 1420 MeV (see Eq. (18) of Ref. [16]). Technically this is implemented using the same formalism as in Ref. [14], where rescattering of the KK to give πη was considered, substituting the KK → πη amplitude by the KK → KK one. We find very small corrections, of the order of 10 ∼ 20%. This might be surprising in the presence of the triangle singularity, but without entering into further details, we should blame for it the memory that the process has of the Schmid theorem [53] (see Ref. [54] for a detailed discussion of the theorem and implications). Another aspect that we investigated was possible modifications of the f 1 (1285) amplitude in the region of 1420 MeV. For this we used directly the amplitude generated in Ref. [18] from the chiral unitary approach in the 1420 MeV region and, again, we found small changes of the same order of magnitude. In summary, the relative strength that we get from the f 1 (1285) peak to the peak at 1420 MeV does not change qualitatively.
While, without studying in detail the different reactions where the f 1 (1420) resonance has been claimed, we cannot dismiss it from the findings of our work, yet, the unavoidable peak that we get at 1420 MeV related to the K * K decay width of the f 1 (1285) should be taken into account in any work aiming at describing these experiments. On the other hand, we propose some reactions where, with moderate uncertainties, we can get the cross sections and predict the strength for these two peaks. Performance of the experiment and comparison with the predictions will be an important step to learn more about the f 1 (1420) peak and its possible origin as we have described in this work, or a consequence of a different dynamics, including its possible nature as an ordinary mesonic resonance.

Conclusions
We have evaluated the cross section for theK p → Y K 0 π + K − reactions, withK =K 0 , K − and Y = , 0 , + in the region of invariant masses around M inv (K 0 π + K − ) ∈ [1200, 1550] MeV. The fact that the f 1 (1285) couples very strongly to the K * K components has as a consequence that the mechanism with K * exchange leads to the production of that resonance in these reactions, which has been already observed in the K − p → f 1 (1285) reaction. The novelty presented here is that, together with the f 1 (1285) excitation, the reaction shows a peak at 1420 MeV in all these reactions, which would be very interesting to observe experimentally. We state that such a mechanism should be taken into account in any study trying to understand the peak at 1420 MeV in the different reactions where the f 1 (1420) resonance has been claimed, to see if by itself it can explain the experimental features observed, in which case the arguments to advocate a new resonance would disappear. On the other hand, it is interesting to propose new reactions, prior to experiment, making predictions along the lines discussed here. The reactions proposed here go in this direction and their experimental implementation and comparison with the predictions done here certainly will help unveiling the nature of the f 1 (1420).
This project has received funding from the European Union Horizon 2020 research and innovation programme under grant agreement No 824093 for the STRONG-2020 project.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: There are no external data associated with the manuscript.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .