Discovery reach for wino and higgsino dark matter with a disappearing track signature at a 100 TeV $pp$ collider
Masahiko Saito, Ryu Sawada, Koji Terashi, Shoji Asai

TL;DR
This paper evaluates the potential of a 100 TeV collider to detect wino and higgsino dark matter particles through disappearing track signatures, aiming to reach the thermal relic density limits.
Contribution
It presents a detailed study of the discovery prospects for wino and higgsino dark matter at a future 100 TeV collider, including detector design modifications and background analysis.
Findings
Potential to discover wino/higgsino dark matter at 100 TeV collider
Strategies to extend detection to thermal relic density limits
Proposed detector modifications to enhance sensitivity
Abstract
Within the theory of supersymmetry, the lightest neutralino is a dark matter candidate and is often assumed to be the lightest supersymmetric particle (LSP) as well. If the neutral wino or higgsino is dark matter, the upper limit of the LSP mass is determined by the observed relic density of dark matter. If the LSP is a nearly-pure neutral state of the wino or higgsino, the lightest chargino state is expected to have a significant lifetime due to a tiny mass difference between the LSP and the chargino. This article presents discovery potential of the 100 TeV future circular hadron collider (FCC) for the wino and higgsino dark matter using a disappearing-track signature. The search strategy to extend the discovery reach to the thermal limits of wino/higgsino dark matter is discussed with detailed studies on the background rate and the reference design of the FCC-hadron detector under…
| Layer # | 1 | 2 | 3 | 4 | 5 |
|---|---|---|---|---|---|
| Default layout (#1) | 25 | 60 | 100 | 150 | 270 |
| Alternative layout (#2) | 25 | 60 | 100 | 150 | 200 |
| Alternative layout (#3) | 25 | 50 | 80 | 110 | 150 |
| Layout | Default (#1) | Alternative (#2) | Alternative (#3) |
|---|---|---|---|
| wino ( = 3 TeV) | |||
| % | % | % | |
| % | % | % | |
| higgsino ( = 1 TeV) | |||
| % | % | % | |
| % | % | % | |
| Layout | Default (#1) | Alternative (#2) | Alternative (#3) |
|---|---|---|---|
| wino ( = 3 TeV) | |||
| Leading jet threshold [TeV] | |||
| threshold [TeV] | |||
| Signal yield | |||
| Background yield | |||
| Significance | |||
| higgsino ( = 1 TeV) | |||
| Leading jet threshold [TeV] | |||
| threshold [TeV] | |||
| Signal yield | |||
| Background yield | |||
| Significance | |||
| Layout | Default (#1) | Alternative (#2) | Alternative (#3) |
|---|---|---|---|
| wino ( = 3 TeV) | |||
| Leading jet threshold [TeV] | |||
| threshold [TeV] | |||
| Signal yield | |||
| Background yield | |||
| Significance | |||
| higgsino ( = 1 TeV) | |||
| Leading jet threshold [TeV] | |||
| threshold [TeV] | |||
| Signal yield | |||
| Background yield | |||
| Significance | |||
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
∎
\thankstext
e1e-mail: [email protected] \thankstexte2e-mail: [email protected]
11institutetext: International Center for Elementary Particle Physics and Department of Physics, The University of Tokyo
Discovery reach for wino and higgsino dark matter with a disappearing track signature at a 100 TeV collider
Masahiko Saito\thanksrefaddr1
Ryu Sawada\thanksrefe1,addr1
Koji Terashi\thanksrefe2,addr1
Shoji Asai\thanksrefaddr1
Abstract
Within the theory of supersymmetry, the lightest neutralino is a dark matter candidate and is often assumed to be the lightest supersymmetric particle (LSP) as well. If the neutral wino or higgsino is dark matter, the upper limit of the LSP mass is determined by the observed relic density of dark matter. If the LSP is a nearly-pure neutral state of the wino or higgsino, the lightest chargino state is expected to have a significant lifetime due to a tiny mass difference between the LSP and the chargino. This article presents discovery potential of the 100 TeV future circular hadron collider (FCC) for the wino and higgsino dark matter using a disappearing-track signature. The search strategy to extend the discovery reach to the thermal limits of wino/higgsino dark matter is discussed with detailed studies on the background rate and the reference design of the FCC-hadron detector under possible running scenarios of the FCC-hadron machine. A proposal of modifying the detector layout and several ideas to improve the sensitivity further are also discussed.
Keywords:
FCC Dark matter Supersymmetry Disappearing track
††journal: Eur. Phys. J. C
1 Introduction
Astrophysical observations of galaxies and the large-scale structure in the universe strongly indicate that dark matter predominate the matter contents in the universe. The nature of dark matter is, however, still unknown and the Standard Model (SM) of particle physics has no counterparticles acting as dark matter. Within the theory of supersymmetry (SUSY), a dark matter candidate is a neutralino, which is a mixed state of neutral supersymmetric partners of the SM U(1)SU(2) gauge bosons and the two SU(2) Higgs doublets. The lightest neutralino () is often assumed to be the lightest supersymmetric particle (LSP). If we assume that dark matter is the neutral wino (higgsino), which is a supersymmetric partner of the SM SU(2) gauge boson (Higgs boson) and that it was produced in thermal processes, the upper limit of the LSP mass is determined by the observed relic density of dark matter to be 3 (1) TeV for the wino (higgsino) LSP scenario HISANO200734 ; CIRELLI2007152 .
Wino LSP is predicted naturally by Anomaly Mediated SUSY Breaking (AMSB) models Giudice:1998xp ; Randall:1998uk and Pure Gravity Mediation (PGM) models IBE2012374 ; PhysRevD.85.095011 ; PhysRevD.87.015028 . Pure wino dark matter has been excluded by indirect searches Cohen_2013 ; Fan2013 ; PhysRevD.94.115019 when cusped dark matter profiles in the Milky Way and dwarf galaxies are assumed. However, uncertainties in the dark matter distributions in galaxies are considerably large Nesti_2013 ; Sugai2016 . When a cored profile is assumed, only a narrow range of the wino dark-matter mass around the Sommerfeld resonance at TeV is excluded Cohen_2013 . Direct dark matter search results do not exclude the wino dark matter because of the extremely small direct-detection cross section of the pure wino HISANO2010311 .
Very-pure higgsino dark matter has been excluded, in the case of very heavy gauginos, due to the absence of inelastic dark matter scattering in direct detection experiments Nagata2015 . If light gauginos are assumed, the higgsino-only dark matter is also excluded due to large mixing under the requirement of a small fine tuning Baer2018 . In such natural SUSY scenarios, the higgsinos are allowed to make up only a portion of the dark matter Baer2018 . Dark matter scenarios with mixed higgsino/gaugino states (so called “well-tempered” ARKANIHAMED2006108 ) are also getting less appealing due to direct search bounds BADZIAK2017226 ; Beneke2017 . However, if a certain amount of fine tuning is allowed, the neutral higgsinos could have a mass splitting due to non-negligible gaugino mixing, and therefore are not excluded by direct searches due to the absence of vector interactions. In this scenario, the nearly-pure higgsino-only dark matter with masses around the thermal limit, which is a subject of this paper, is not constrained Nagata2015 ; Kowalska2018 by direct searches due to the small neutralino-nucleon cross section or by indirect searches because of the negligible Sommerfeld enhancement in the annihilation cross section.
The mass difference of the neutral wino and the charged winos is dominated by the radiative contributions from SM particles when other SUSY particles are a few times heavier than the winos. Contributions from heavier SUSY particles are negligible Ibe:2012sx due to a suppression by the large masses. As shown in Ref. Ibe:2012sx , the mass difference is nearly constant (160 MeV) for the wino with masses larger than about GeV. The SUSY particle contributions to the mass difference of higgsinos are non-negligible. However, when gauginos are heavier than (100) TeV considered in this paper, the mass difference of higgsinos is almost constant ( MeV) because SM contributions dominate Nagata2015 . Due to the small mass difference, the has a long lifetime of approximately 0.2 (0.023) ns for the 3 TeV wino (1 TeV higgsino) at the thermal limit 1126-6708-2003-03-045 .
With this mass difference, the decays into a and a charged pion predominantly. The passes through the detector without any electromagnetic or strong interactions while the pion is not reconstructed as a track in the inner-tracking detector due to very low transverse momentum (). Therefore, such a long-lived chargino is observed in the detector as a short ‘disappearing’ track, which has hits in several silicon-detector layers near the interaction point and no associated hits after a certain radius, with a typical length of (1–10) cm in collider experiments PhysRevD.55.330 ; Chen:1999yf ; Asai2008185 .
In this paper, a search strategy for the long-lived charginos with a disappearing-track signature at the 100 TeV FCC-hadron machine (FCC-hh) is discussed with the reference FCC-hadron detector. The expected sensitivity is presented under the scenario of 30 ab*-1* of data collected under possible running conditions of the FCC-hh. In earlier studies on the prospects of disappearing-track searches in future colliders Low2014 ; Cirelli2015 ; Mahbubani2017 ; FUKUDA2018306 ; PhysRevD.98.035026 ; Cao2018 , effects of multiple collisions occurring simultaneously with a signal event (pileup) were assumed not to alter the results significantly by improving the analysis method or the background rate was scaled with a certain factor without detailed studies. In this paper, the background rates were estimated based on the detector layout and simulation. In addition, several ideas that could potentially improve the sensitivity further are discussed.
2 Simulation samples
Simulated samples of Monte Carlo (MC) events are used to obtain the kinematic distributions of the signal and background processes. The wino mass spectrum is calculated using Softsusy 3.7.3 Allanach:2001kg assuming the minimum Anomaly Mediated SUSY Breaking (AMSB) model with tan, positive sign of the higgsino-mass term and the universal scalar mass () of 20 TeV. This setting naturally accommodates the observed Higgs boson mass of 125 GeV without maximal mixing. The higgsino mass spectrum is calculated by Susyhit Djouadi:2006bz assuming the general minimal supersymmetric extension of the SM (MSSM) with tan. Other SUSY particle masses are set to 100 TeV to prevent the higgsino from being affected from loop effects. Only the pure wino and higgsino LSP scenarios are considered here. The production cross sections for the winos and higgsinos in 100 TeV collisions are calculated using Prospino2 Beenakker:1996ch . The kinematic distributions are simulated using MadGraph5_aMC@NLO 2.3.3 Alwall:2014 assuming the wino spectra. The same kinematic distributions are assumed for the wino and higgsino if they have the same masses. In this study, pair-production processes of charginos and neutralinos are considered: and for the wino scenario, and , and for the higgsino scenario. As the and are degenerate in the higgsino case, the difference in the cross section or decay processes between and is negligible. Therefore, only the process is generated and the cross section is doubled to account for . Minimum-bias collisions are produced using Pythia 8.230 Sjostrand:2014zea and are overlaid to evaluate contribution from non-genuine tracks caused by the random combination of hits in the inner-tracking detector, which are referred to as fake tracks hereafter. The details of the minimum-bias event overlay are described in Sect. 5.
SM background processes such as -boson production process in association with jets, top-quark pair and single top-quark production processes are simulated using MadGraph5_aMC@NLO and the generated samples are used to evaluate the selection efficiency. The kinematic distributions of jets, electrons, muons and the magnitude of missing transverse momentum () are processed using Delphes 3.4.2 deFavereau2014 with the reference FCC-hh run card to account for the detector reconstruction efficiency and energy/momentum measurement. Jets are reconstructed from particle-flow candidates using the anti- algorithm Cacciari:2008gp with a distance parameter of 0.4.
3 Inner Tracker Layout
The reconstruction of short tracks relies on the performance of the inner-most tracking system. This paper considers the reference design of the FCC-hadron detector FCC-hh_detector as a baseline. For the inner tracker, the most relevant part for the disappearing-track search is the inner-pixel layers, in particular those in the barrel region where the signal acceptance is high. In the present study three different layouts are considered for the barrel region to evaluate impact of the inner-tracker configuration.
The default layout (labelled #1) is the one in the reference design, containing four pixel layers within mm in radial distance from the beam line. In the second layout (#2) a pixel layer is added at mm to the default layout and all the other layers are unchanged. The third layout (#3) considers an additional pixel layer at mm and the inner four pixel layers are moved closer to the beam line. The - views of all the three layouts are shown in Fig. 1.111 We use a right-handed coordinate system, with the -, -, and -axis are defined by the horizontal, the vertical and the beam directions respectively. Its origin is at the nominal interaction point. The pseudorapidity is defined in terms of the polar angle by .
The outer macro-pixel layers located at mm in the barrel region or the forward tracking system are unchanged from the reference design for all the three layouts. The first five inner layers are most relevant for the short-track reconstruction as described below. Their radial positions considered in this study are summarized in Table 1.
4 Event Selection
The chargino/neutralino-pair production is characterized by the presence of significant associated with the and one or two disappearing tracks from the . The SM background is suppressed by requiring an existence of a high- disappearing track. The signal event will presumably contain high- jet(s) as well, originating from initial state radiation, and they provide the means of triggering on the event in conjunction with . The analysis therefore selects events which contain at least one high- jets, large and no isolated lepton (electron or muon) with GeV. Figure 2 shows the leading jet and distributions of signal and background events selected with the isolated lepton veto and no requirement on the jet or . The thresholds on the jet and are determined, for each configuration of the inner-tracker layouts, by maximizing the sensitivity based on the signal acceptance and background rate (summarized later in Tables 3 and 4). Although the background yield can be reduced significantly by applying high thresholds to the leading jet and , there will be no discovery power for the wino or higgsino at the thermal limit if the disappearing-track information is not used in the search.
The chargino-track reconstruction requires the particle to traverse a certain minimum number of inner-tracking layers before decay. The ATLAS search with 36 fb*-1* of Run 2 data at TeV Aaboud:2017mpt used a short-track reconstruction with only hits in the four pixel layers. As the chargino lifetime becomes short, the acceptance for the chargino-track reconstruction increases with decreasing radial position of the outermost layer necessary for the track reconstruction. On the other hand, the use of short tracks with less number of layer hits will usually increase fake tracks (discussed in Sect. 5). In order to assess the impact on the sensitivity from different requirements on the number of hits used in the track reconstruction (), the track reconstruction is performed with the requirements of and separately. When the (5), the four (five) inner-most layers with hits are used to reconstruct tracks and the fifth (sixth) layer or above is not considered in the track reconstruction. The contribution from fake tracks increases significantly at high due to higher hit occupancy. Figure 3 shows normalized distributions of fake tracks in minimum-bias collision events (details in Sect. 5.2) and chargino tracks in 3 TeV wino signal events. Given the dramatic increase of fake background with increasing , the analysis requires a candidate track to be within (further discussion is in Sect. 7).
The signal acceptance for the track requirement is estimated based on the tracker geometry and the chargino lifetime by assuming that the chargino tracks can be reconstructed at 100% of time if the charginos traverse at least four or five inner-most layers before decay. Table 2 summarizes the signal acceptances for the wino and higgsino under the three scenarios of the inner-tracker layouts.
5 Background estimation
Backgrounds for the disappearing-track signature are categorized into two components Aaboud:2017mpt . The first one is physical background, which arises from charged particles (mainly electrons or charged pions) scattered by the material in the inner tracker. The second one is unphysical background, which arises from fake tracks. The estimate of both background components is discussed separately below.
5.1 Physical background
The physical background due to particle scattering is estimated by assuming that it is solely determined by the amount of detector material passed through by particles. Under this assumption the physical background rate at the FCC-hh is obtained by scaling the background rate measured in ATLAS Aaboud:2017mpt by the ratio of the material budgets in the FCC-hh and ATLAS inner trackers ATLASMaterialBudget . A support structure of the FCC-hh tracker has not been fully defined yet, so the expected impact from particle scattering with support material is accounted for by assuming the same relative amounts of tracker layers and support material measured in ATLAS. The physical background is dominated by +jets processes, which have large and a high- track associated with the lepton from the boson decay. Therefore, the final estimate of the physical background is obtained by correcting the scaled background rate above for the cross section and kinematic selection efficiency for +jets processes at TeV.
5.2 Background with fake tracks
The fake-track background is evaluated by counting the number of reconstructed tracks with good quality using simulated minimum-bias events. The samples of minimum-bias events are generated with the average number of interactions per bunch crossing () of 200 or 500, and simulated for the FCC-hadron detector using Geant4 Agostinelli:2002hh . The two values are assumed to be representative of FCC-pileup conditions. Pileup collisions are produced with two different models of soft QCD processes: the first model considers only non-diffractive processes while the second one the mixture of diffractive and non-diffractive processes with the relative fraction determined according to the Pythia cross sections222The non-diffractive processes are set using the SoftQCD:nonDiffractive option while the diffractive ones using SoftQCD:singleDiffractive and SoftQCD:doubleDiffractive options in Pythia.. The former model predicts more fake tracks than the latter due to the absence of diffractive events, and is used as the nominal pileup model to be more conservative.
The track reconstruction is performed as follows. First, the track-seed finding is done by finding a group of inner-tracker hits within a narrow region with a size of about in and . Second, a track is reconstructed by fitting a line with a given curvature to the collection of the hits. The fit is iterated by removing hits which are far from the line. The following quality criteria are required to select chargino-track candidates: GeV, , , at least one hit in each layer and the association to the hard-scattering vertex by applying the nominal impact-parameter requirement (transverse and longitudinal impact parameters of mm and mm). The final number of fake tracks is obtained by counting the number of tracks with loose impact-parameter requirements ( mm, mm) and scaling that to the nominal impact-parameter requirement assuming a uniform distribution of and . Figure 4 shows the estimated probability of finding a fake track per event as a function of the . The fake-track finding probability increases significantly with increasing . However, the probability can be reduced by 2–3 orders of magnitude by changing the requirement from 4 to 5. The number of fake-track backgrounds is estimated by multiplying the number of SM background events satisfying the kinematic selection criteria (in Sect. 4) by these probabilities.
The time resolution of pixel detectors can be better than 50 ps by utilizing, for example, low-gain avalanche detectors Pellegrini2014 ; Lange2017 . Therefore, the time information could be used in the track fit as an additional track parameter to potentially improve the quality of the track. The reduction of fake-track background by requiring a good time quality (i.e, small in the estimation of track time obtained from the associated pixel hits) is investigated using the simulated minimum-bias samples with . In the simulation, the , , positions and the time of collisions are distributed randomly according to Gaussian probability-density functions with the standard deviations of 0.5 mm, 0.5 mm, 50 mm and 160 ps, respectively (without assuming any correlations among them). The is computed assuming the constant time resolution of 50 ps for a single layer hit. The fake-track background is found to be reduced by 96% by requiring the -probability to be larger than 0.05 for tracks reconstructed with when collisions occur in bunch crossings separated by 25 ns. The signal efficiency for this selection is approximately 95%, evaluated using a sample of single muon events.
6 Result
The discovery sensitivity with the requirement of is higher than that with due to much lower 5-layer-hit fake-track background, as indicated in Fig. 4. The fake-track background rate appears to be rather insensitive to the inner-tracker layout for a fixed requirement. Therefore, the alternative layout #2 is expected to have sensitivity in-between the default and alternative layout #3 (hence not shown in the rest of the figures) because the signal acceptance is a factor 2–4 smaller for the #2 than #3 (see Table 2).
Figure 5 shows the expected discovery significance for the wino (higgsino) LSP model with the proper lifetime of 0.2 (0.023) ns with 30 ab*-1* with the requirement of , assigning 30% systematic uncertainty in the background yields. The background uncertainty is based on the estimate in the ATLAS Run 2 analysis Aaboud:2017mpt . The time information discussed in Sect. 5.2 is not used in Fig. 5. The discovery significance is obtained using an approximate formula to calculate the significance of a signal+background hypothesis against a background-only hypothesis, taking into account the background uncertainty Linnermann2003 ; Cranmer2005 . The sensitivity is degraded significantly when the increases from 200 to 500, even if five layer hits are used in the track reconstruction. The alternative layout (#3) clearly improves the sensitivity, reaching well above the discovery for the 3 TeV wino in both pileup scenarios. For the 1 TeV higgsino the discovery is also feasible in both pileup scenarios but with less margin for the high pileup case. The discovery significance is re-evaluated using the time information and the result is shown in Fig. 6. The sensitivity is restored at the high pileup scenario and the degradation with increasing is much more relaxed. The discovery is also reached with sufficient margin for the 1 TeV higgsino with 30 ab*-1* or even less if the time information is used. The discovery sensitivities for the three tracker layouts at the pileup scenarios of and 500 (without time information) are summarized in Table 3 and 4, respectively. The reach of the FCC-hadron machine is shown in Fig. 7 as well as the HL-LHC reach and the current observed limits by the ATLAS experiment. According to the present study the FCC-hadron machine has potential for answering conclusively yes or no to the thermal production of nearly-pure wino or higgsino dark matter.
7 Conclusions and discussions
The discovery potential for the 3 TeV wino and 1 TeV higgsino scenarios has been evaluated using a disappearing-track signature in collisions at TeV with the FCC-hadron detector. The FCC-hh sensitivity is very promising to both the wino and higgsino scenarios, reaching up to the thermal limits, with inner-tracker layouts optimized for a short-track reconstruction, even in the high pileup environment of . If the time resolution of 50 ps is achievable at pixel-hit level, wino or higgsino dark matter will be confirmed or refuted up to the thermal limits with an integrated luminosity of 30 ab*-1* or smaller.
Finally, we discuss several ideas that could potentially improve the sensitivity further. The range could be expanded from (used in the present analysis) to e.g, , to improve the overall signal acceptance, then split into smaller regions to optimize the sensitivity. As shown in Fig. 3, the contribution from fake tracks, that grows very rapidly with increasing , will have to be evaluated for the determination of the regions. By tilting the sensor modules and making overlaps between them, the charginos produced at the interaction region may pass through multiple modules in a single layer. This will potentially help improving the acceptance for a short-lifetime signal without decreasing the required number of hits in the track reconstruction, therefore without increasing the fake-track background level. It would be also beneficial if the barrel pixel layers can be moved even closer to the beam line than the alternative layout #3. This has a large impact on the higgsino signal acceptance. Another interesting avenue to pursue is to use the speed of particles relative to the speed of light (), which can be measured using the time information of pixel hits, to reject background and characterize the observed disappearing-track signature. With the single muon sample, the resolution is found to be about 14% (with the pixel-hit time resolution of 50 ps) for the tracks with under the alternative tracker layouts. The estimated resolution is roughly proportional to the pixel-hit time resolution. The fake-track background will be reduced by requiring the measured to be compatible with a particle (produced from the primary collision vertex) with a certain range of momentum and mass. Because new heavy particles move more slowly in the detector than SM particles, the could be used as an additional discriminant to separate the signal from scattered SM particles.
Acknowledgements.
This work was supported by MEXT KAKENHI Grant number JP16K21730 and JSPS KAKENHI Grant Number JP18J11405.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1(1) J. Hisano, S. Matsumot, M. Nagai, O. Saito, M. Senami, Phys. Lett. B 646 , 34 (2007). DOI 10.1016/j.physletb.2007.01.012. ar Xiv:hep-ph/0610249
- 2(2) M. Cirelli, A. Strumia, M. Tamburini, Nucl. Phys. B 787 (1), 152 (2007). DOI 10.1016/j.nuclphysb.2007.07.023. ar Xiv:0706.4071 [hep-ph]
- 3(3) G.F. Giudice, M.A. Luty, H. Murayama, R. Rattazzi, JHEP 12 , 027 (1998). DOI 10.1088/1126-6708/1998/12/027. ar Xiv:hep-ph/9810442
- 4(4) L. Randall, R. Sundrum, Nucl. Phys. B 557 , 79 (1999). DOI 10.1016/S 0550-3213(99)00359-4. ar Xiv:hep-th/9810155
- 5(5) M. Ibe, T.T. Yanagida, Phys. Lett. B 709 (4), 374 (2012). DOI https://doi.org/10.1016/j.physletb.2012.02.034. ar Xiv:1112.2462 [hep-ph]
- 6(6) M. Ibe, S. Matsumoto, T.T. Yanagida, Phys. Rev. D 85 , 095011 (2012). DOI 10.1103/Phys Rev D.85.095011. ar Xiv:1202.2253 [hep-ph]
- 7(7) B. Bhattacherjee, B. Feldstein, M. Ibe, S. Matsumoto, T.T. Yanagida, Phys. Rev. D 87 , 015028 (2013). DOI 10.1103/Phys Rev D.87.015028. ar Xiv:1207.5453 [hep-ph]
- 8(8) T. Cohen, M. Lisanti, A. Pierce, T.R. Slatyer, J. Cosmol. Astropart. Phys. 2013 (10), 61 (2013). DOI 10.1088/1475-7516/2013/10/061. ar Xiv:1307.4082 [hep-ph]
