arxiv: v1 [hep-ph] 14 Sep 2018

Transcription

1 From ± roduction aymmetry at the LHC to romt ν τ at IceCube Victor P. Goncalve Intituto de Fíica e Matemática, Univeridade Federal de Pelota (UFPel), Caixa Potal 354, CEP 96-9, Pelota, R, Brazil Rafa l Maciu la and Antoni zczurek Intitute of Nuclear Phyic, Polih Academy of cience, Radzikowkiego 5, PL-3-34 Kraków, Poland arxiv:89.544v [he-h] 4 e 8 he decrition of the heavy meon roduction at large energie and forward raiditie at the LHC i fundamental to derive realitic rediction of the romt atmoheric neutrino flux at the IceCube Obervatory. In articular, the romt tau neutrino flux i determined by the decay of meon roduced in comic ray - air interaction at high energie and large value of the Feynman- x F variable. Recent datafrom the LHCb Collaboration indicate a roduction aymmetry for and meon, which cannot be exlained in term of the tandard modelling of the hadronization roce. In thi aer we demontrate that thi aymmetry can be decribed auming an aymmetric trange ea ((x) (x)) in the roton wave function and taking into account of the charm and trange fragmentation into meon. Moreover, we how that the trange quark fragmentation contribution i dominant at large - x F (.3). he romt ν τ flux i calculated and the enhancement aociated to the trange quark fragmentation contribution, diregarded in reviou calculation, i etimated. PAC number: Keyword: Introduction. he exerimental reult obtained in recent year by the LHC, the Pierre Auger and IceCube Neutrino Obervatorie have challenged and imroved our undertanding of Particle Phyic. he dicovery of the Higg boon [] comleted the tandard Model (M), which i now a comlete and elf-conitent theory. On the other hand, the detection of atrohyical neutrino by the IceCube Neutrino Obervatory et the begining of the neutrino atronomy []. In addition, the data from the Pierre Auger Obervatory rovide an unique oortunity to tet Particle Phyic at energie well beyond current accelerator [3]. uch reult motivated, in articular, the develoment of new and/or more recie aroache to decribe the erturbative and nonerturbative regime of the Quantum Chromodynamic (QC). One examle i the recent imrovement in the decrition of the heavy meon roduction in hadronic colliion at the LHC, directly influenced by the need to contrain the magnitude of the romt neutrino flux, which i crucial for a recie determination of the comic neutrino flux at the IceCube [4 7]. In what follow we will exlore thi direct connection between the LHC and IceCube reult and rovide more recie rediction for the roduction at the LHC and the romt tau neutrino flux at the IceCube. he atmoheric neutrino are roduced in comic-ray interaction with nuclei in Earth atmohere[8]. At low neutrino energie (E ν 5 GeV), thee neutrino arie alo at Univerity of Rzezów, PL Rzezów, Poland Electronic addre: barro@ufel.edu.br Electronic addre: rafal.maciula@ifj.edu.l Electronic addre: antoni.zczurek@ifj.edu.l from the decay of light meon (ion and kaon), and the aociated flux i denoted a the conventional atmoheric neutrino flux [9]. On the other hand, in the energy range 5 GeV < E ν < 7 GeV, it i exected that the romt atmoheric neutrino flux aociated with the decay of hadron containing heavy flavour become imortant []. In the articular cae of the tau neutrino ν τ flux, it i dominated at low energie by the conventional atmoheric flux, via ν µ ν τ ocillation. On the otherhand, fore ν > 4 GeV, thi contribution become negligible and the romt ν τ flux i determined by the decay of meon, which have a letonic decay channel τν τ with a branching ratio of a few ercent, with the ubequent τ decay that alo contribute to the flux []. A recie determination of the romt ν τ flux i fundamental to identify the tau neutrino of comic origin, which i conidered another imortant ignature of the comic ray origin of the highet neutrino flux. A demontrated in Ref. [6], the romt neutrino flux i determined by the heavy meon roduction at high energie and very forward raiditie. herefore, the decrition of the roduction in the kinematical region robed by the LHCb Collaboration i a requiite to obtain a recie rediction of the romt ν τ flux. uring the lat year, the LHCb Collaboration releaed a large et of data aociated with the and B meon roduction. he data for the tranvere momentum and raidity ditribution are, in general, quite well decribed by theoretical aroache. On the other hand, the decrition of the exerimental data for the charge roduction aymmetrie [, 3] till i a challenge for the great majority of the theoretical aroache. In general, thee roduction aymmetry are interreted a ariing during the nonerturbative roce of hadronization a imle-

2 mented e.g. in the PYHIA Monte Carlo, which i baed on the Lund tring fragmentation model. However, thi aroach fail to decribe the recent for the ± roduction aymmetry, which have found evidence of a nonzero aymmetry. In articular, in contrat with PYHIA that redict a oitive value, the exerimental data indicate that the aymmetry i negative. Very recently, two of u have rooed in Ref. [4] an alternative aroach to decribe the aymmetry reent in the and roduction []. he baic idea i that ubleading contribution for the arton fragmentation are nonnegligible at the LHC energie and that the aymmetry come from the inherent aymmetry of the u and d valence ditribution in the incident roton. In thi aer we extend the aroach for the ± roduction and demontrate that the can be decribed if we aume that the trange quark ea in the roton i aymmetric, with. uch aymmetry i redicted e.g. by the Meon Cloud Model (ee. e.g. Ref. [5 9]) and i not excluded by the recent data and by QC global analyi. In reality, the trange ea in the roton i only oorly known, with it behavior being determined in a greatextentbythe neutrino-induced I dataoncharm roduction obtained by the CCFR/NueV and NOMA exeriment [, ]. Recent analyi of the LHC data for the W ± roduction imroved our undertanding of the trange ditribution, eecially at mall - x, but the exitence or not of an aymmetric trange ea i till an oen quetion []. A a conequence, auming that our aroach for the ubleading arton fragmentation i correct, the reult for the ± aymmetrycan be conidered a a firt ignature that in the roton. Finally, the reult reented in Ref. [4] and in what follow indicate that the ubleading contribution to the arton fragmentation are nonnegligible and have a large imact on the Feynman - x F ditribution of the heavy meon roduced in hadronic colliion. A dicued e.g. in Ref. [5, 6], thi ditribution determine the romt neutrino flux. herefore, it i exected that the romt ν τ flux will be enhanced by thee ubleading contribution. One of the goal of thi aer i to etimate thi enhancement and rovide realitic rediction for the tau neutrino flux that are baed on a formalim that i in agreement with the recent. roduction at the LHC. At high energie the charm quark are roduced dominantly by gluon - gluon interaction via the gg c c ubroce and are believed to hadronize to - meon mainly through the c fragmentation roce. herefore, at leading order, we exect an identical amount of and meon, which imlie that the charge aymmetry defined by ) = σ ) σ ) σ )σ ) () will be zero at thi aroximation. Conequently, the aymmetrie are exected to be generated by ubleading artonic ubrocee, initial tate aymmetrie and/or a ditinct decrition of the hadronization roce. Here we extend the aroach rooed in Ref. [4], which exlain the / aymmetry in term of the unfavored fragmentation function, which are reonible for light quark/antiquark fragmentation to meon, for roduction. he fact that u d in the incident roton, naturallyleadtothe / aymmetrywhen the ubleading contribution for the fragmentation i taken into account. In contrat, in the cae of roduction, the incluion of the ubleading fragmentation, aociated to the and tranformation, imlie ) = if =. herefore, our interretation of the i that the ± aymmetry arie due to an aymmetry in the trange quark ea. A exlained before, uch a behavior i redicted by ome theoretical model and i not excluded by the recent QC global analyi of the exerimental data. In articular, the CEQ Collaboration ha erformed a dedicated tudy of the trange arton ditribution of the roton in Ref. [3] and obtained that the exerimental data i quite well decribed alo auming that. In the reent aer we hall include the dominant g g c c ubroce a well a the / g / g and g / g / term with the trange arton from Ref. [3]. Moreover, we will include the ubleading / ± fragmentation, which i decribed in term of the robability of tranition of a trange quark into a meon (P ). Of coure, the leading fragmentation i aociated with the c and c tranitionwith anaociatedtranition branching of 6 % [4]. A in Ref. [4] the gluon - gluon channel will be calculated within k -factorization aroach, with the cro ection being exreed in term of the unintegrated gluon ditribution function (UGF) and off-hell matrix element for the gg c c ubroce. In the reent aer we ue the Kimber-Martin-Rykin (KMR) recrition for UGF [5] which wa hown to allow good decrition of charm roduction at the LHC [7]. In contrat, the ( )g ( )g and g( ) g( ) rocee are calculated within a leading-order collinear aroach, with the regulation of mall tranvere momenta region being done a in Ref. [4]. Moreover, the c fragmentation i decribed uing the Peteron model (with ε =.5), while the and fragmentation function are arametrized a in Ref. [4] uing the revered-peteron and triangular fragmentation function. he only free arameter in our aroach i P. In rincile, it i exected to be larger than the value obtained in Ref. [4] for the u/d tranition, due to the larger ma of the trange quark. However, a thi quantity i aociated with a nonerturbative roce, it i not oible to calculate it value from firt rincile. In what follow, we will contrain P uing the for the charge roduction aymmetry ).

3 3.5. < y < 3. = 7 ev.5. < y < 3. = 8 ev.5 Peteron(-z) (olid) riangle (dahed) Peteron(-z) (olid) riangle (dahed) P = trangene aymmetry bet-fit CEQ P = trangene aymmetry bet-fit CEQ < y < 3.5 = 7 ev.5 3. < y < 3.5 = 8 ev.5 Peteron(-z) (olid) riangle (dahed) Peteron(-z) (olid) riangle (dahed) P = trangene aymmetry bet-fit CEQ P = trangene aymmetry bet-fit CEQ < y < 4.5 = 7 ev < y < 4.5 = 8 ev.5 Peteron(-z) (olid) riangle (dahed) Peteron(-z) (olid) riangle (dahed) P = trangene aymmetry bet-fit CEQ P = trangene aymmetry bet-fit CEQ6 FIG. : aymmetry obtained with our aroach for = 7 ev (left anel) and = 8 ev (right anel) together with the LHCb exerimental data [3]. In Fig. we reent our reult for the aymmetry for = 7 ev (left anel) and = 8 ev (right anel) uing the revered-peteron and triangular fragmentation function for the tranition. Rather reaonable agreement with the i obtained auming P = 7. In other word, a mall value for the unfavoured fragmentation function i ufficient to decribe the. he data tatitic i till too low to erform a more detailed fit and/or to dicriminate between the two model for the ubleading fragmentation. However, the reult indicate that the aymmetry in trangene in the roton wave function, a decribed in the CEQ arametrization, i able to generate the correct ign for ), in contrat to PYHIA [3], a well a the enhancement of the aymmetry at larger raiditie. In Fig. (left anel) we reent the reulting rediction for the tranvere momentum ditribution of for the different range of the meon

4 4 (µb/gev) dσ/d -m k -fact. KMR CEQ6L ± BR(c ) = 6% - ) X (GeV) = 7 ev. < y <.5; m =.5 < y < 3.; m = 3. < y < 3.5; m = 3.5 < y < 4.; m = 3 4. < y < 4.5; m = 4 (µb) dσ/dx F BR(c ) =.6 =.7 P -3 - k -fact.: c c (olid line) -4 - aymmetric: (dah-dotted) -5 aymmetric: (dahed) ± meon x F = 7 ev Peteron FF CEQ6 FIG. : Left: ranvere momentum ditribution of for different range of raidity. he [6] are hown for comarion. Right: Feynman - x F ditribution for the ± roduction. raidity (y ) robed by the LHCb Collaboration [6]. A quite well agreement with the i achieved without free arameter. We have verified that the contribution of the ubleading fragmentation for the - ectra i mall ( 5%) in the kinematical range robed by the LHCb Collaboration. In contrat, the behavior of the raidity and Feynman - x F ditribution are ignificantly modified at large value of y and x F. In Fig. (right anel) we demontrate that the aymmetry in the trange ea imly different behavior for the x F - ditribution of the and meon at intermediate x F. More imortant, while at mall - x F the conventional contribution dominate, at large - x F the ituation i revered. One ha that the contribution aociated to the g g channel become dominant for x F.5. In articular,forx F.3, thechannelinitiatedbytrange quark, uually diregardedin the analyi of the roduction, are dominant. uch value of x F correond to raiditie larger than thoe robed by the LHCb detector. However, a demontrated in Ref. [6], thi i exactly the kinematical region that determine the behavior of the romt neutrino flux. Conequently, the reence of the ubleading contribution for the roduction i exected to have direct imact on the rediction of the romt tau neutrino flux. Promt ν τ flux at the IceCube. One of the current goal of the IceCube Obervatory i the meaurement of tau-neutrino [8], which are conidered an indeendent robe of the comic neutrino. uch an exectation i trongly motivated by the fact that for comic neutrino the decay of charged ion generated in atrohyical ource imlie a ratio ν e : ν µ : ν τ = : : at the Earth, while for atmoheric neutrino thi ratio i exected to be tyically ν e : ν µ : ν τ = : :. [5, ]. A a conequence, the background aociated to atmoheric tau neutrino i uually redicted to be trongly reduced in comarion to the other flavour, with the meaurement of a tau neutrino being conidered a direct robe of comic neutrino. However, reviou analyi have diregarded the ubleading contribution dicued in thi aer. It i the aim of the reent tudy to make a realitic etimate of the romt ν τ ν τ flux when uing the current information from the LHC. In order to etimate the romt ν τ ν τ flux, we will cloely follow the rocedure decribed in detail in Ref. [6, 7]. We will calculate the romt tau neutrino flux uing the emi-analytical Z-moment aroach [], where a et of couled cacade equation for the nucleon, heavy meon and leton (and their antiarticle) fluxe i olved, with the equation being exreed in term of the nucleon-to-hadron (Z NH ), nucleon-tonucleon(z NN ), hadron-to-hadron(z HH ) andhadron-toneutrino (Z Hν ) Z-moment. hee moment are inut in the calculation of the romt tau neutrino flux aociated with the roduction of a meon and it decay into a ν τ in the low-and high-energyregime. We will focu on vertical fluxe and will aume that the comic ray flux φ N can be decribed by the H3a ectrum rooed in Ref. [9], with the incident flux being rereented by roton. A in Ref. [] we will include in our calculation the contribution of neutrino roduced in the direct ν τ decay a well a thoe generated in the chain decay τ ν τ. he contribution for the romt ν τ flux aociated to the decay of meon heavier than t i negligible [5] and will not be included in our analyi. In Fig.3 (left) we how the flux of the romt ν τ ν τ onent, flux caled by Eν 3. In addition to the conventional com- aociatedtoheavyquarkroductionbyagluongluon fuion and rereented by the olid black line, we

5 5-3 E 3 Φ ν [GeV cm - - r - ] -4-5 IceCube Conv. Conv. ubleading (ymmetric) Conv.ubleading (Aymmetric) Ratio IceCube ymmetric Aymmetric H3a E ν (GeV) E ν (GeV) FIG. 3: Left: he romt atmoheric tau neutrino (ν τ ν τ) flux a a function of the neutrino energy. Right: Ratio between the full calculation (conventional ubleading contribution) and the conventional one. how the reult which include in addition the ubleading fragmentation auming a ymmetric (dahed red line) or an aymmetric (dot - dahed blue line) trange ea in the roton wave function. he ubleading mechanim lead to a ignificant enhancement of the highenergy romt τ-neutrino flux, which can be quantified calculating the ratio between the full rediction (conventional ubleading) and the conventional one. he energy deendence of thi ratio i reented in Fig.3(right), which demontrate that the incluion of the ubleading contribution imlie an enhancement of a factor 3.5 in the romt ν τ flux in the kinematical range robed by the IceCube Obervatory. In articular, for a aymmetric trange ea, needed to decribe the, we redict that the enhancement will be of a factor of about 3.5 at E ν = 5 GeV. he enhancement trongly increae for E ν > 8 GeV due to the fater decreaing of the gluon ditribution in comarion to the trange one at large value of Bjorken x variable. Concluion. In the reent aer we rooe, for the firt time, the decrition of the roduction aymmetry for and meon in term of an aymmetry in the trange ea of the roton aociated to the incluion of the ubleading fragmentation mechanim or. We have ued aymmetric ditribution derived in a global analyi of different exerimental data and that can be exlained within meon cloud icture of the nucleon. We have demontrate that a mall value for the trange fragmentation function into meon i ufficient to decribe the. uch a new contribution, diregarded in reviou tudie, become dominant at large value of x F, that i the kinematical range that determine the romt atmoheric ν τ flux at the IceCube Obervatory. We have etimated the imact of thi contribution and demontrated that it imlie an enhancement by a factor larger than 3 in the kinematical range robed by the IceCube Obervatory. Our reult indicate that a future exerimental analyi of romt tau neutrino at the IceCube can be ueful to robe the underlying mechanim of roduction at high energie and forward raiditie at the LHC. We have found recently that the roduction of τ- neutrino wa dicued very recently in the context of intrinic charm in the nucleon [3] and the beam dum fixed target exeriment HiP at CERN[3]. No reference to the LHCb aymmetry wa done there. Acknowledgment. hi tudy wa artially uorted by the Polih National cience Center grant EC- 4/5/B//58, by the Center for Innovation and ranfer of Natural cience and Engineering Knowledge in Rzezów and by the Brazilian funding agencie CNPq, FAPERG and INC-FNA (roce number /4-5). [] G. Aad et al. (ALACollaboration), Phy. Lett. B76, ();. Chatrchyan et al.(cm collaboration), Phy. Lett. B 76, 3 (). [] M. G. Aarten et al. (IceCube Collaboration), cience 34, 4856 (3). [3] A. Aab et al. (Pierre Auger Collaboration), Phy. Rev. Lett. 7, no. 9, 9 (6). [4] M. V. Garzelli,. Moch and G. igl, JHEP 5, 5 (5); A. Bhattacharya, R. Enberg, M. H. Reno, I. arcevic and A. tato, JHEP 56, (5);

6 6 R. Gauld, J. Rojo, L. Rottoli and J. albert, JHEP 5, 9(5); R.Gauld, J. Rojo, L.Rottoli,.arkar and J. albert, JHEP 6, 3 (6); F. Halzen and L. Wille, Phy. Rev. 94, 44 (6); R. Laha and. J. Brodky, Phy. Rev. 96, 3 (7); M. V. Garzelli et al. [PROA Collaboration], JHEP 75, 4 (7); M. Benzke, M. V. Garzelli, B. Kniehl, G.Kramer,.Moch andg.igl, JHEP 7, (7). [5] A. Bhattacharya, R. Enberg, Y.. Jeong, C.. Kim, M. H. Reno, I. arcevic and A. tato, JHEP 6, 67 (6). [6] V. Goncalve, R. Maciu la, R. Paechnik and A. zczurek, Phy. Rev. 96 (7) 946. [7] A. V. Giannini, V. P. Goncalve and F.. Navarra, Phy. Rev. 98, no., 4 (8). [8]. Gaier and F. Halzen, Ann. Rev. Nucl. Part. ci. 64, (4) [9] M. Honda,. Kajita, K. Kaahara,. Midorikawa and. anuki, Phy. Rev. 75, 436 (7). [] P. Gondolo, G. Ingelman and M. hunman, Atroart. Phy. 5, 39 (996). [] L. Paquali and M. H. Reno, Phy. Rev. 59, 933 (999) [] R. Aaij et al. (LHCb Collaboration), Phy. Lett. B 78, 9 (3). [3] R. Aaij et al. (LHCb Collaboration), J. High Energy Phy. 8, 8 (8). [4] R. Maciula and A. zczurek, Phy. Rev. 97, no. 7, 74 (8). [5] H. Holtmann, A. zczurek and J. eth, Nucl. Phy. A569, 63 (996). [6] G. Q. Feng, F. G. Cao, X. H. Guo and A. I. ignal, Eur. Phy. J. C 7, 5 (). [7] E.R. Cazaroto, V.P. Goncalve, E.. Navarra and M. Nielen, Phy. Lett. B74, 8 (3). [8] X. G. Wang, C. R. Ji, W. Melnitchouk, Y. alamu, A. W. homa and P. Wang, Phy. Lett. B 76, 5 (6); Phy. Rev. 94, no. 9, 9435 (6). [9] M. Gohariour, Nucl. Phy. A 973, 6 (8). [] A. O. Bazarko et al. (CCFR Collaboration), Z. Phy. C 65, 89 (995). [] G.P. Zeller et al. (NueV Collaboration), Phy. Rev. 65, 3(R) (), [] A. M. Cooer-arkar and K. Wichmann, Phy. Rev. 98, no., 47 (8) [3] H. L. Lai, P. M. Nadolky, J. Pumlin,. tum, W. K. ung and C.-P. Yuan, J. High Energy Phy. 4, 89 (7). [4] M. Liovyi, A. Verbytkyi and O. Zenaiev, Eur. Phy. J. C 76, no. 7, 397 (6). [5] G. Watt, A.. Martin and M. G. Rykin, Eur. Phy. J. C 3, 73 (3). [6] R. Aaij et al. (LHCb Collaboration), Nucl. Phy. B 87, (3). [7] R. Maciula and A. zczurek, Phy. Rev. 87, no. 9, 94 (3). [8] M. G. Aarten et al. (IceCube Collaboration), Phy. Rev. 93, no., (6); M. Uner [IceCube Collaboration], Po ICRC 7, 974 (8). [9]. K. Gaier, Atroart. Phy. 35, 8 (). [3] W. Bai and M. Hall Reno, arxiv: [3] M. Anelli at al. (HiP Collaboration), arxiv: