Research Article  Open Access
Pseudorapidity Distribution of Charged Particles and Square Speed of Sound Parameter in or Collisions over an Energy Range from 0.053 to 7 TeV
Abstract
Pseudorapidity distributions of charged particles produced in protonproton () or protonantiproton () collisions over an energy range from 0.053 to 7 TeV are studied by using the fourcomponent Landau hydrodynamic model. The results calculated by the model are in agreement with the experimental data of the UA5, PHOBOS, UA1, P238, CDF, ALICE, and CMS Collaborations which present orderly from low to high energies. According to the distribution widths of different components, the values and some features of square speed of sound parameter for “participant” and “spectator” quark components are obtained. It is shown that the speed of sound for “participant” quark components agrees approximately with that for “spectator” quark components in the error ranges. The present work is useful for studying nucleusnucleus collisions in the related energy range.
1. Introduction
In the laboratory conditions, the only way for creating and studying the interacting systems with hadronic or partonic degrees of freedom at extremely high energy and density is to investigate heavyion (nucleusnucleus) collisions at very high centerofmass energies. The relativistic heavy ion collider (RHIC) in the United States and the large hadron collider (LHC) in Switzerland have been built, respectively [1, 2]. Particularly, the LHC was originally designed to accelerate two protons up to total energy of 14 TeV, and it will accelerate heavy ions to collide at centerofmass energy () of 5.5 TeV. Such high energy collisions offer a new environment which leads to a new significant extension of the kinematic range in longitudinal rapidity and transverse momentum [3]. It also provides a new chance for us to understand systematically the particle statistical behavior, production process, interaction mechanism, and correlation phenomena.
The multiplicity and (pseudo)rapidity distributions of finalstate particles can often be used to test different theoretical models and wonderful ideas. The rapidity distribution and the pseudorapidity distribution of charged particles are very important quantities for the study of mainly particle production mechanism in high energy nucleonnucleon and nucleusnucleus collisions. The studies of and in protonproton () or protonantiproton () collisions not only can provide reference and baseline for heavyion collisions, but also can reflect the main particles contribution to heavyion collisions. In  or  collisions at present colliders, the centerofenergies reach a region of ultrahigh energies. Meanwhile, the highest energy nuclear experiments can be performed. Most of highdensity experimental data can help us realize the creation of a new state of matter, namely, the quark gluon plasma (QGP) [4, 5], which is a thermalized system consisting of strong coupling quarks and gluons in a very small region. It is believed that the QGP was the main matter in our early universe shortly after the Big Bang.
Many models have been introduced in the field of high energy collisions, for example, the hydrodynamics model [6], the thermodynamics model [7], the hadron resonance gas model [8], the fireball model [9, 10], and so forth. Generally speaking, a given model treats the different collisions by different considerations, and different models approach the same collisions in different methods. In our previous work [11], based on the participantspectator model [12] and Landau hydrodynamic model [13–17], the pseudorapidity distributions of charged particles produced in nucleusnucleus collisions at the RHIC and LHC have been analyzed, and the values of square speed of sound () have been obtained. We notice that the related data for  or  collisions are absent in our previous work [11].
To make up the limitation of our previous work [11], in this work, we will use the same method to study the pseudorapidity distributions of charged particles produced in  or  collisions over an energy range from 0.053 to 7 TeV and to extract the values of speed of sound parameter for the “participant” and “spectator” quark components which correspond, respectively, to the participant and spectator nucleons in nucleusnucleus collisions. The quoted experimental data are taken from the UA5 ( collisions at 0.053, 0.2, 0.546, and 0.9 TeV) [18], PHOBOS ( collisions at 0.2 and 0.41 TeV) [19], UA1 ( collisions at 0.54 TeV) [20], P238 ( collisions at 0.63 TeV) [21], CDF ( collisions at 0.63 and 1.8 TeV) [22], ALICE ( collisions at 0.9 and 2.36 TeV) [23], and CMS ( collisions at 0.9, 2.36, and 7 TeV) [24, 25] Collaborations.
2. The Model
The picture of participantspectator model [12] is the base of many theoretical models such as the nuclear fireball model [26], the twocomponent model [27], the nuclear overlap model [28], and the multisource thermal model [29, 30]. In the framework of participantspectator model [12], the projectile and target nuclei penetrating through each other in collisions. Then, a central participant region, a target spectator region, and a projectile spectator region are expected to form in the collisions. Further, we can divide the participant region into two parts: a target participant region and a projectile participant region. In the multisource thermal model [29, 30], these regions can extend to different cylinders in rapidity space. In addition, the contributions of leading nucleons can fall under the spectator regions. In high energy nucleusnucleus collisions, both the participant and spectator regions consist of nucleons.
In high energy hadronhadron collisions, the situation is different. If we use the picture of participantspectator model [12], both the participant and spectator regions will be quarks. Similarly, we have four source components: the target spectator quark component, the target participant quark component, the projectile participant quark component, and the projectile spectator quark component, which are presented orderly from low to high in the rapidity space. According to [31, 32], if we use the Landau hydrodynamic model [13–17] to describe the hadronhadron collisions, only the central pseudorapidity region can be fitted. This means that the Landau hydrodynamic model [13–17] gives a narrow pseudorapidity distribution, and the contributions of leading nucleons are not included in the model. To describe the pseudorapidity distribution in a wide region, we can apply the model for each quark component in the hadronhadron collisions.
The pseudorapidity distribution of charged particles produced in each quark component is given by a Gaussian function [17, 33] with where is the normalization constant, is in the units of GeV, denotes the speed of sound, denotes the proton mass in the units of GeV/ with , and is the pseudorapidity distribution width.
Strictly speaking, there are some differences between the rapidity distribution and the pseudorapidity distribution. In fact, the differences between the two distributions at high energies can be neglected due to the small values [34]. If we need to distinguish the two distributions, we should calculate them severally or make a complete transformation [35]. In the present work, we do not distinguish the two distributions for the purpose of simpleness. Equation (1) can be rewritten as the function of due to (2). The finalstate pseudorapidity distribution is a weighted sum of four equations (3) with different rapidity shifts and weights.
For  or  collisions, the rapidity shift () of target spectator is opposite to the rapidity shift () of projectile spectator, and the rapidity shift () of target participant is opposite to the rapidity shift () of projectile participant. The weight () of target participant is equal to that of projectile participant, and the weight of target spectator is equal to that of projectile spectator.
3. Comparisons with Experimental Data
Figures 1(a)–1(d) show the pseudorapidity distributions of charged particles produced in inelastic  collisions at = 0.053, 0.2, 0.546, and 0.9 TeV, respectively. The symbols in the range of represent the experimental data of the UA5 Collaboration [18], and those in the range of are symmetrical reflection at the midpseudorapidity. The curves are our calculated results with the fourcomponent Landau hydrodynamic model. The values of parameters obtained by fitting the experimental data are given in Table 1 with the values of ( per degree of freedom), where and denote the square speeds of sound for the participants and spectators, respectively. One can see that the calculated results are in agreement with the experimental data at the four energies. Both values of for the participants and spectators are approximately equal to each other in the error ranges, and they do not depend obviously on in the considered energy range.

(a)
(b)
(c)
(d)
The pseudorapidity distributions of charged particles produced in inelastic  collisions at and 0.41 TeV measured by the PHOBOS Collaboration [19],  collisions at TeV measured by the UA1 Collaboration [20],  collisions at TeV measured by the P238 Collaboration with ranging from 1.5 to 5.5 [21], and  collisions at and 1.8 TeV measured by the CDF Collaboration with [22] are presented in Figures 2(a)–2(f), respectively. The symbols and curves represent the same meanings as those in Figure 1. By fitting the experimental data, the obtained values of and are given in Table 1. Once again the fourcomponent Landau hydrodynamic model describes the experimental data well. Although the energy has a large change, the values of for the participants and spectators are approximately the same in the error ranges, and they have no obvious change with changing the energy.
(a)
(b)
(c)
(d)
(e)
(f)
Figure 3 gives the pseudorapidity distributions of charged particles produced in nonsinglediffractive  collisions at = 0.9, 2.36, and 7 TeV. The symbols in the left panel represent the experimental data obtained by the CMS Collaboration [24, 25], and those in the right panel stand for the experimental data of the ALICE Collaboration [23]. The curves are our calculated results. The obtained values of and are given in Table 1. We notice that the model describes the experimental data at the LHC energies. Both values of for the participants and spectators seem to be the same in the error ranges, and they do not depend obviously on the energy.
(a)
(b)
(c)
(d)
(e)
To see clearly the dependences of and on , the values of and in Table 1 are displayed in Figure 4. Different symbols represent the results extracted from different experiments as marked in the figure. The open and closed symbols represent the values of and , respectively. One can see that both the and fall mainly into a range consisting of two lines: and . These results are in agreement with the hadron resonance gas model which results in, respectively, and for the existing regions of hadron resonances including and excluding pions at the temperature of ~85 MeV and for the two existing regions at the temperature of ~190 MeV [36, 37]. The present results are also in agreement with the lattice quantum chromodynamics theory which gives and 0.31 at the temperatures of 125 and 400 MeV, respectively [38].
Both values of for the participants and spectators are approximately the same in the error ranges and they do not depend obviously on over an energy range from 0.053 to 7 TeV. From Figure 4, one can also see that seems to be less than at TeV energies. Because the fitted errors are large for both , it is hard to determine which one is less than the other.
The correlations between and as well as and are presented in Figures 5(a) and 5(b), respectively, where is in the units of GeV. The symbols represent the values of rapidity shifts and which are obtained in Figures 1–3 and listed in Table 1. The lines are our fitted results which are described by with and with , respectively. One can see the linear relation existing between the rapidity shifts and . With the increase of the logarithmic centerofmass energy, both rapidity shifts increase slightly.
(a)
(b)
The correlation between and is shown in Figure 6, where is in the units of GeV. The symbols represent the values of contribution ratio of the projectile (or target) participant quark component which are obtained in Figures 1–3 and listed in Table 1. The line is our fitted result described by with . One can see the linear relation existing between the contribution ratio and . Particularly, the contribution ratio of participants decreases and that of spectators increases with the increase of the logarithmic centerofmass energy.
4. Conclusions
To conclude, the pseudorapidity distributions of charged particles produced in  collisions at , 0.41, 0.9, 2.36, and 7 TeV and in  collisions at , 0.2, 0.54, 0.546, 0.63, 0.9, and 1.8 TeV have been studied in the present work. The calculated results of the fourcomponent Landau hydrodynamical model satisfactorily describe the experimental data of the UA5, PHOBOS, UA1, P238, CDF, ALICE, and CMS Collaborations. In the model, we have considered the application of the participantspectator model [12] in hadronhadron collisions. Then, the interacting system of hadronhadron collisions is divided into four components: the target and projectile “spectator” quark components and the target and projectile “participant” quark components. For each component, the Landau hydrodynamical model [13–17] is applied, which results in the pseudorapidity distribution being a Gaussian type.
Based on the above successful descriptions of the pseudorapidity distributions of charged particles produced in  or  collisions over an energy range from 0.053 to 7 TeV, the square speeds of sound, , for both the participants and spectators are obtained and found to be the same in the fitted error ranges. Both are approximately independent of in the considered energy range and fall mainly into the range of 0.12–0.30 which is in agreement with hadron resonance gas model which gives at the temperature of ~85 MeV and at ~190 MeV [36, 37]. The values of obtained in the present work are also in agreement with the lattice quantum chromodynamics theory which gives in the temperature range of 125–400 MeV [38].
In the process of extracting the speed of sound, the rapidity shifts and contribution ratio of the participant and spectator quark components are naturally obtained. There are linear relations existing between and , and , as well as and . Both rapidity shifts increase slightly with the increase of the logarithmic centerofmass energy. The contribution ratio of participants decreases and that of the spectators increases with the increase of the logarithmic centerofmass energy.
5. Discussions
The basis for Landau’s equation (essentially a Gaussian pseudorapidity distribution) is the applicability of Landau hydrodynamics, namely, the creation of an initial state of dense Lorentz contracted matter that expands later according to the equations of 2 + 1 dimensional relativistic hydrodynamics. It seems that such an initial state is not created in hadronhadron collision at intermediate or not too high energy, and the evolution of the interacting system would not be covered by hydrodynamics. However, hadronhadron collision at high or ultrahigh energy can be regarded as a small dense system comparing with nucleusnucleus collisions. This small dense system can produce a multihadron finalstate state and can be approximately described by Landau’s equation.
The square speed of sound extracted from the present work is in agreement with that from nucleusnucleus collisions at RHIC and LHC energies in our previous work [11]. The similar values of speed of sound in both hadronhadron and nucleusnucleus collisions provides the insight into some common laws and universal mechanisms in multihadron productions [31, 32]. We believe that there are more similar quantities in multihadron productions in the “elementary” and “complex” collisions.
We would like to point out that the present work is the first one which extracts the square speed of sound in hadronhadron collisions from the pseudorapidity distributions of charged particles over a very wide energy range. What we use is the fourcomponent Landau hydrodynamic model which is based on the participantspectator model [12] and the Landau hydrodynamic model [13–17]. Our work can provide references and baselines for nucleusnucleus collisions at high energies. We hope to extract the square speed of sound in protonproton and leadlead collisions at the highest energy at the LHC in the near future.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was partly finished at the State University of New York at Stony Brook, USA. One of the authors (FuHu Liu) thanks Professor Dr. Roy A. Lacey and the members of the Nuclear Chemistry Group of Stony Brook University for their hospitality. The authors acknowledge the supports of the National Natural Science Foundation of China (under Grant nos. 10975095, 11247250, and 11005071), the China National Fundamental Fund of Personnel Training (under Grant no. J1103210), the Open Research Subject of the Chinese Academy of Sciences LargeScale Scientific Facility (under Grant no. 2060205), the Shanxi Scholarship Council of China, and the Overseas Training Project for Teachers at Shanxi University.
References
 S. Chatrchyan, V. Khachatryan, A. M. Sirunyan, A. Tumasyan, W. Adam, T. Bergauer et al., “Indications of suppression of excited $\gamma $ states in PbPb collisions at $\sqrt{{s}_{NN}}=2.76$ TeV,” Physical Review Letters, vol. 107, no. 23, Article ID 052302, 15 pages, 2011. View at: Google Scholar
 A. A. Alves Jr., L. M. A. Filho, A. F. Barbosa, I. Bediaga, G. Cernicchiaro, G. Guerrer et al., “The LHCb detector at the LHC,” Journal of Instrumentation, vol. 3, no. 8, Article ID S08005, 217 pages, 2008. View at: Google Scholar
 R. Aaij, C. A. Beteta, A. Adametz et al. et al., “Measurement of prompt hadron production ratios in pp collisions at $\sqrt{s}=0.9$ and 7 TeV,” European Physical Journal C, vol. 72, Article ID 2168, 19 pages, 2012. View at: Publisher Site  Google Scholar
 K. J. Eskola, K. Kajantie, P. V. Ruuskanen, and K. Tuominen, “Scaling of transverse energies and multiplicities with atomic number and energy in ultrarelativistic nuclear collisions,” Nuclear Physics B, vol. 570, no. 12, pp. 379–389, 2000. View at: Google Scholar
 J.P. Blaizot, “Theoretical overview: towards understanding the quarkgluon plasma,” Journal of Physics G: Nuclear and Particle Physics, vol. 34, no. 8, pp. S243–S251, 2007. View at: Publisher Site  Google Scholar
 B. M. Waldhauser, D. H. Rischke, J. A. Maruhn, H. Stöcker, and W. Greiner, “Influence of the nuclear bulk properties and the MIT bag constant on the phase transition to the quark gluon plasma,” Zeitschrift für Physik C: Particles and Fields, vol. 43, no. 3, pp. 411–417, 1989. View at: Publisher Site  Google Scholar
 G. F. Bertsch and S. Das Gupta, “A guide to microscopic models for intermediate energy heavy ion collisions,” Physics Reports, vol. 160, no. 4, pp. 189–233, 1988. View at: Google Scholar
 F. Karsch, K. Redlich, and A. Tawfik, “Hadron resonance mass spectrum and lattice QCD thermodynamics,” The European Physical Journal C, vol. 29, no. 4, pp. 549–556, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 W. Y. Zhang, “Jets induced in emulsions and cloud chambers by cosmic ray particles of energy (${10}^{11}{10}^{14}$ eV),” Acta Physica Sinica, vol. 17, no. 8, pp. 271–295, 1961. View at: Google Scholar
 G. D. Westfall, J. Gosset, P. J. Johansen et al., “Nuclear fireball model for proton inclusive spectra from relativistic heavyion collisions,” Physical Review Letters, vol. 37, no. 18, pp. 1202–1205, 1976. View at: Publisher Site  Google Scholar
 L.N. Gao, Y.H. Chen, H.R. Wei, and F.H. Liu, “Speed of sound parameter from RHIC and LHC heavyion data,” Advances in High Energy Physics, vol. 2013, Article ID 450247, 8 pages, 2013. View at: Publisher Site  Google Scholar
 R. J. Glauber, “Highenergy collision theory,” in Lectures of Theoretical Physics, W. E. Brittin and L. G. Dunham, Eds., vol. 1, pp. 315–414, Interscience, New York, NY, USA, 1959. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 L. D. Landau, “Multiple production of particles under collision of rapid particles,” Izvestiya Akademii Nauk: Series Fizicheskikh, vol. 17, p. 51, 1953. View at: Google Scholar
 L. D. Landau, “Multiple production of particles under collision of rapid particles,” in Collected Papers of L. D. Landau, D. TerHaarp, Ed., pp. 569–585, Pergamon Press, Oxford, UK, 1965. View at: Publisher Site  Google Scholar
 “On multiple production of particles during collisions of fast particles,” in QuarkGluon Plasma: Theoretical Foundations, J. Kapusta, B. Mueler, and J. Rafelski, Eds., p. 283, Elsevier, Amsterdam, The Netherlands, 2003. View at: Google Scholar
 S. Z. Belen'kji and L. D. Landau, “Hydrodynamic theory of multiple production of particles,” Il Nuovo Cimento, vol. 3, no. 1, supplement, pp. 15–31, 1956. View at: Publisher Site  Google Scholar
 E. V. Shuryak, “Multiparticle production in high energy particle collisions,” Yadernaya Fizika, vol. 16, no. 2, pp. 395–405, 1972. View at: Google Scholar
 G. J. Alner, K. Alpgård, P. Anderer et al., “UA5: a general study of protonantiproton physics at $\sqrt{S}=546$ GeV,” Physics Reports, vol. 154, no. 56, pp. 247–383, 1987. View at: Google Scholar
 B. Alver, B. B. Back, M. D. Baker, M. Ballintijn, D. S. Barton, R. R. Betts et al., “Chargedparticle multiplicity and pseudorapidity distributions measured with the PHOBOS detector in Au+Au, Cu+Cu, d+Au, and p+p collisions at ultrarelativistic energies,” Physical Review C, vol. 83, no. 2, Article ID 024913, 24 pages, 2011. View at: Google Scholar
 G. Arnison, A. Astbury, B. Aubert, C. Bacci, R. Bernabei, A. Bézaguet et al., “Charged particle multiplicity distributions in protonantiproton collisions at 540 GeV centre of mass energy,” Physics Letters B, vol. 123, no. 12, pp. 108–124, 1983. View at: Publisher Site  Google Scholar
 R. Harr, C. Liapis, P. Karchin et al., “Pseudorapidity distribution of charged particles in pp collisions at $\sqrt{s}=630$ GeV,” Physics Letters B: Nuclear, Elementary Particle and HighEnergy Physics, vol. 401, no. 12, pp. 176–180, 1997. View at: Google Scholar
 F. Abe, D. Amidei, G. Apollinari, G. Ascoli, M. Atac, P. Auchincloss et al., “Pseudorapidity distributions of charged particles produced in $\stackrel{}{p}p$ interactions at $\sqrt{s}=630$ and 1800 GeV,” Physical Review D, vol. 41, no. 7, pp. 2330–2333, 1990. View at: Publisher Site  Google Scholar
 K. Aamodt, N. Abel, U. Abeysekara, A. A. Quintana, A. Abramyan, D. Adamová et al., “Chargedparticle multiplicity measurement in proton–proton collisions at $\sqrt{s}=0.9$ and 2.36 TeV with ALICE at LHC,” The European Physical Journal C, vol. 68, no. 12, pp. 89–108, 2010. View at: Publisher Site  Google Scholar
 V. Khachatryan, A. M. Sirunyan, A. Tumasyan, W. Adam, T. Bergauer, M. Dragicevic et al., “Transverse momentum and pseudorapidity distributions of charged hadrons in pp collisions at $\sqrt{s}=0.9$ and 2.36 TeV,” Journal of High Energy Physics, no. 2, Article ID 41, 35 pages, 2010. View at: Publisher Site  Google Scholar
 V. Khachatryan, A. M. Sirunyan, A. Tumasyan, W. Adam, T. Bergauer, M. Dragicevic et al., “Transversemomentum and pseudorapidity distributions of charged hadrons in pp collisions at $\sqrt{s}=7$ TeV,” Physical Review Letters, vol. 105, no. 2, Article ID 022002, 14 pages, 2010. View at: Google Scholar
 J. Gosset, H. H. Gutbrod, W. G. Meyer et al., “Central collisions of relativistic heavy ions,” Physical Review C, vol. 16, no. 2, pp. 629–657, 1977. View at: Publisher Site  Google Scholar
 L. Grandchamp and R. Rapp, “Charmonium suppression and regeneration from SPS to RHIC,” Nuclear Physics A, vol. 709, no. 1–4, pp. 415–439, 2002. View at: Publisher Site  Google Scholar
 S. Eremin and S. Voloshin, “Nucleon participants or quark participants?” Physical Review D, vol. 67, no. 6, Article ID 064905, 3 pages, 2003. View at: Google Scholar
 F.H. Liu, N. N. A. Allah, and B. K. Singh, “Dependence of black fragment azimuthal and projected angular distributions on polar angle in siliconemulsion collisions at 4.5A GeV/c,” Physical Review C, vol. 69, no. 5, Article ID 057601, 4 pages, 2004. View at: Publisher Site  Google Scholar
 F.H. Liu and J.S. Li, “Isotopic production cross section of fragments in ${}^{\text{56}}\text{F}$e+p and ${}^{\text{136}}\text{X}$e(${}^{\text{124}}\text{X}$e)+Pb reactions over an energy range from 300A to 1500A MeV,” Physical Review C, vol. 78, no. 4, Article ID 044602, 13 pages, 2008. View at: Publisher Site  Google Scholar
 E. K. G. Sarkisyan and A. S. Sakharov, “Relating multihadron production in hadronic and nuclear collisions,” The European Physical Journal C, vol. 70, no. 3, pp. 533–541, 2010. View at: Publisher Site  Google Scholar
 E. K. G. Sarkisyan and A. S. Sakharov, “Multihadron production features in different reactions,” in Proceedings of the 35th International Symposium on Multiparticle Dynamics and the Workshop on Particle Correlations and Femtoscopy, vol. 828 of AIP Conference Proceedings, pp. 35–41, Kroměříž, Czech Republic, August 2005. View at: Publisher Site  Google Scholar
 C. Alt, T. Anticic, B. Baatar et al. et al., “Energy dependence of particle ratio fluctuations in central Pb+Pb collisions from $\sqrt{{s}_{NN}}=6.3$ to 17. 3 GeV,” Physical Review C, vol. 79, no. 4, Article ID 044910, 11 pages, 2009. View at: Google Scholar
 G. Wolschin, “Pseudorapidity distributions of produced charged hadrons in pp collisions at RHIC and LHC energies,” Europhysics Letters, vol. 95, no. 6, Article ID 61001, 6 pages, 2011. View at: Publisher Site  Google Scholar
 C. Y. Wong, Introduction to HighEnergy HeavyIon Collisions, World Scientific, Singapore, 1994.
 A. Tawfik and H. Magdy, “Hadronic equation of state and speed of sound in thermal and dense medium,” http://arxiv.org/abs/1206.0901. View at: Google Scholar
 A. Tawfik and T. Harko, “Quarkhadron phase transitions in the viscous early universe,” Physical Review D, vol. 85, no. 8, Article ID 084032, 19 pages, 2012. View at: Google Scholar
 Sz. Borsányi, G. Endrődi, Z. Fodor et al., “QCD equation of state at nonzero chemical potential: continuum results with physical quark masses at order ${\mu}_{2}$,” Journal of High Energy Physics, vol. 2012, no. 8, Article ID 53, 14 pages, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 YaQin Gao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.