Heavy-ion Physics at a Fixed-Target Experiment Using the LHC Proton and Lead Beams (AFTER@LHC): Feasibility Studies for Quarkonium and Drell-Yan Production
B. Trzeciak, C. Da Silva, E.G. Ferreiro, C. Hadjidakis, D. Kikola,, J.P. Lansberg, L. Massacrier, J. Seixas, A. Uras, Z. Yang

TL;DR
This paper explores the feasibility of using the LHC's fixed-target mode with lead beams to study heavy-ion collisions, focusing on quarkonium and Drell-Yan production, and presents simulation results for potential measurements.
Contribution
It demonstrates the potential of AFTER@LHC to extend heavy-ion physics studies into new energy and rapidity domains, especially for quarkonium and Drell-Yan observables.
Findings
Feasibility of studying $Upsilon(nS)$, $J/Cpsi$, and $Cpsi(2S)$ production in new domains.
Simulations show potential for analyzing nuclear modification of partonic densities.
Supports the use of quarkonia as probes of quark-gluon plasma formation.
Abstract
We outline the case for heavy-ion-physics studies using the multi-TeV lead LHC beams in the fixed-target mode. After a brief contextual reminder, we detail the possible contributions of AFTER@LHC to heavy-ion physics with a specific emphasis on quarkonia. We then present performance simulations for a selection of observables. These show that , and production in heavy-ion collisions can be studied in new energy and rapidity domains with the LHCb and ALICE detectors. We also discuss the relevance to analyse the Drell-Yan pair production in asymmetric nucleus-nucleus collisions to study the factorisation of the nuclear modification of partonic densities and of further quarkonia to restore their status of golden probes of the quark-gluon plasma formation.
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8
Figure 9
Figure 10
Figure 11
Figure 12| (mb) | |
|---|---|
| 5.51 10-4 | |
| 5.90 10-5 | |
| (1S) | 6.87 10-7 |
| (2S) | 1.82 10-7 |
| (3S) | 7.56 10-8 |
| Drell-Yan (M 1.5 GeV/c2) | 6.66 10-6 |
| 1.07 10-1 | |
|
9.44 10-5 ()
5.86 10-5 () |
|
| minimum bias | 25.38 |
| Beam | Target gas | |||
|---|---|---|---|---|
| Pb | Xe | 30 | 165.1 | |
| (Kr, Ar, Ne) | ||||
| Xe | 2 | 3.7 | ||
| Kr, Ar, Ne | 3 , 6.5 , 150 | |||
| H | 250 | 1 |
| Yields | signal | S/B | |
|---|---|---|---|
| 1.33 | 29.0 | ||
| Xe | 1.39 | 7.8 | |
| PbXe | 4.33 | 1.8 | |
| 2.92 | 8.2 | ||
| Xe | 3.06 | 2.2 | |
| PbXe | 9.56 | 5.0 | |
| 1.37 | 10.3 | ||
| Xe | 1.44 | 2.8 | |
| PbXe | 4.49 | 6.2 |
| Yields | : 4-5 GeV/ | : 5-6 GeV/ | ||
|---|---|---|---|---|
| signal () | () | signal () | () | |
| : 2-3 | 11.5 | 2.6 | 4.8 | 6.3 |
| : 3-4 | 26.9 | 0.2 | 15.6 | 0.4 |
| : 4-5 | 15.3 | 0.3 | 8.3 | 1.0 |
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.
∎
11institutetext: B. Trzeciak 22institutetext: Institute for Subatomic Physics, Utrecht University, Utrecht, The Netherlands
C. Da Silva 33institutetext: P-25, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
E.G. Ferreiro 44institutetext: Departamento de Física de Partículas, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
C. Hadjidakis, J.P. Lansberg, L. Massacrier 55institutetext: IPNO, CNRS-IN2P3, Univ. Paris-Sud, Université Paris-Saclay, Orsay, France
D. Kikola 66institutetext: Faculty of Physics, Warsaw University of Technology, Warsaw, Poland
J. Seixas 77institutetext: LIP and IST, Lisbon, Portugal
A. Uras 88institutetext: IPNL, Lyon, France
Z. Yang 99institutetext: Tsinghua University, Beijing, China
Heavy-ion Physics at a Fixed-Target Experiment Using the LHC Proton and
Lead Beams (AFTER@LHC): Feasibility Studies for Quarkonium and Drell-Yan Production
B. Trzeciak
C. Da Silva
E.G. Ferreiro
C. Hadjidakis
D. Kikola
J.P. Lansberg
L. Massacrier
J. Seixas
A. Uras
Z. Yang
Abstract
We outline the case for heavy-ion-physics studies using the multi-TeV lead LHC beams in the fixed-target mode. After a brief contextual reminder, we detail the possible contributions of AFTER@LHC to heavy-ion physics with a specific emphasis on quarkonia. We then present performance simulations for a selection of observables. These show that , and production in heavy-ion collisions can be studied in new energy and rapidity domains with the LHCb and ALICE detectors. We also discuss the relevance to analyse the Drell-Yan pair production in asymmetric nucleus-nucleus collisions to study the factorisation of the nuclear modification of partonic densities and of further quarkonia to restore their status of golden probes of the quark-gluon plasma formation.
††journal: Few Body Systems
Contents
-
3.3 Drell-Yan pair production in heavy-ion collisions at 72 GeV
-
3.4 Acceptances of the ALICE detectors in the fixed-target mode
1 Introduction
Hadro-production experiments in the fixed-target mode have indisputably played a fundamental role in the history of quarkonium physics. This began with the co-discovery of the Aubert:1974js in 1974 at the Brookhaven National Laboratory and followed by the discovery of the Herb:1977ek and the first observation of Armstrong:1992ae at the Fermi National Laboratory. The dedicated heavy-ion program at the Super Proton Synchrotron (SPS) at CERN then uncovered many novel and unexpected features of quark and gluon dynamics, including the anomalous suppression of Abreu:2000ni in PbPb collisions, and helped discover the strong non-factorising nuclear suppression of the hadro-production at high Hoyer:1990us . Let us recall that the first observation of -pair production was carried out at the SPS as early as in 1982 Badier:1982ae . Fixed-target experiments have also played a central role in measuring Drell-Yan pair production, in particular below the bottomonium region Baldit:1994jk ; Hawker:1998ty ; Vasilev:1999fa ; Zhu:2006gx .
In this context, we find it important to remind that collisions of the proton and heavy-ion LHC beams on fixed targets open a remarkably large range of physics opportunities Brodsky:2012vg ; Koshkarev:2016ket ; Signori:2016lvd ; Koshkarev:2016acq ; Lansberg:2016gwm ; Lansberg:2016urh ; Signori:2016jwo ; Pisano:2015fzm ; Vogt:2015dva ; Fengand:2015rka ; Barschel:2015mka ; Kikola:2015lka ; Kurepinand:2015jka ; Zhou:2015wea ; Arleo:2015lja ; Lansberg:2015lva ; Brodsky:2015fna ; Massacrier:2015qba ; Anselmino:2015eoa ; Lansberg:2015kha ; Ceccopieri:2015rha ; Goncalves:2015hra ; Kanazawa:2015fia ; Lansberg:2015hla ; Massacrier:2015nsm ; Lansberg:2014myg ; Chen:2014hqa ; Rakotozafindrabe:2013au ; Lansberg:2013wpx ; Lansberg:2012sq ; Lorce:2012rn ; Rakotozafindrabe:2012ei ; Boer:2012bt ; Lansberg:2012wj ; Lansberg:2012kf ; Liu:2012vn . The AFTER@LHC project bears on them and aims at an ambitious heavy-ion, spin and hadron-structure physics programme. These opportunities result from a high luminosity – be it with an internal gas target or with an extracted beam– together with a relatively high center-of-mass system (c.m.s.) energy of 115 GeV per nucleon with a 7 TeV proton beam and 72 GeV per nucleon with a 2.76 A.TeV lead beam. Besides, the boost typical of the fixed-target mode makes forward detectors such as LHCb and ALICE nearly ideal to probe the negative region, which is essentially uncharted despite its interest.
The AFTER@LHC energy range stands half way between the c.m.s. energies of SPS and the Relativistic Heavy Ion Collider (RHIC), allowing, in particular, for detailed studies of the bottomonium production and dynamics and offering a very large physical acceptance for detectors covering between 1 to 3 units of rapidity.
In Massacrier:2015qba ; Kikola:2017hnp , we have focused on the possibility offered by the LHC nine-months-per-year proton program during which AFTER@LHC would be able to study the production of quarkonia in , and collisions with an unparalleled statistical accuracy down to . Such high-precision quarkonium-production measurements in are essential Andronic:2015wma ; Brambilla:2010cs ; Lansberg:2006dh to solve longstanding puzzles in and production and for the understanding of their behavior in the (Cold) Nuclear Matter (CNM).
AFTER@LHC can also provide extremely relevant information on the Quark-Gluon Plasma (QGP), which should be created at 72 GeV in Pb collisions. Despite the decrease in the c.m.s. energy, yields similar to that of LHC at 2.76 and 5.5 TeV and RHIC at 200 GeV are to be expected. These are two orders of magnitude larger than the ones obtained at RHIC at 62 GeV. The first results from the LHC at 2.76 and 5.02 TeV Andronic:2003zv confirmed that the pattern of quarkonium (anomalous) suppression at high energy is very intricate with subtle and dependencies. Low energy experiments, where the recombination process Andronic:2003zv are not expected to be significant, can then play a key role in understanding the underlaying physics processes. Such a role will become central if it is possible to carry out measurements of and even production and suppression in heavy-ion collisions –2 measurements thus far not performed in any other experimental configuration. The quest for the sequential suppression of quarkonia as a QGP thermometer would then become realistic again.
In addition, the measurement of the Drell-Yan-pair production in asymmetric nucleus-nucleus collisions provides a unique opportunity to test the factorisation of the initial-state nuclear effect such as the nuclear PDFs. Such a factorisation is routinely assumed at high energies but a number of non-linear effects can violate it. If these are significant, the nuclear effects from both colliding nucleus do not add up linearly and this would simply prevent us to use the information gained from proton-nucleus or electron-nucleus collisions to characterise the initial stage of nucleus-nucleus collisions. Because of a very high background at the LHC in the collider mode, such measurements are not even feasible in PbPb collisions. At the LHC, the studies of asymmetric nucleus-nucleus systems are not planned, and such a measurement is very difficult at RHIC. With the lower energy at AFTER@LHC, the studies may be within reach if sufficient different systems are used (Pb, Pb, PbXe, Xe, …) to pin down and quantify factorisation-violation effects.
The structure of the paper is as follows. In section 2, we introduce the case for heavy-ion physics for AFTER@LHC. In section 3, we discuss our simulation set-up and present performance studies for quarkonium and Drell-Yan probes. We conclude in section 4.
2 Quark-gluon-plasma studies with AFTER@LHC
2.1 General considerations
One of the main incentives to study relativistic heavy-ion collisions is the unique opportunity to probe highly-excited and dense nuclear matter under controlled laboratory conditions. If the temperature and density are high enough, Quantum ChromoDymanics (QCD), the theory of the strong interaction, predicts the existence of a new phase of matter, the so-called quark-gluon plasma (QGP), in which most of the quarks and gluons, normally confined within hadrons, are liberated. Lattice-QCD calculations indeed show that there is a rapid rise of the entropy density when the temperature reaches 160 MeV. Beyond this temperature the effective number of degrees of freedom saturates near the number of quark and gluon helicity and colour states and the entropy density approaches the Stefan-Boltzmann limit.
Relativistic heavy-ion collisions offer an ideal environment to reach the conditions of the phase transition, both in terms of the necessary temperature and volume for the system thermalisation. Even at moderately high energies such as the ones obtained at the CERN SPS, the average multiplicity in proton-proton () collisions at a c.m.s. energy GeV is on the order of 3 per unit rapidity with an average momentum close to GeV/. This leads to a volume for the system of approximately fm3 and an energy density on the order of 0.4 GeV.fm*-3*, that is roughly 3 times the density of the normal nuclear matter. Heavy nucleus-nucleus () collisions, at the same energy per nucleon, rather leads to an energy density of GeV.fm*-3* and an initial volume, before the expansion, on the order of 150 fm3.
Among the expected signatures of the QGP formation, quarkonia play an essential role. Since the seminal work of T. Matsui and H. Satz in 1986 Matsui:1986dk , the behaviour of quarkonia at high temperature has been considered a signal for deconfinement in heavy-ion collisions. The original idea behind this proposal is that Debye colour screening in a deconfined plasma would reduce the binding of quarks, and thus also affect the formation of heavy-quark bound states. This effect on quarkonia emerged in the QGP would result in the suppression of the quarkonia yields relative to those in the absence of the plasma formation.
The experimental verification of the Matsui-Satz prediction since the late 1980’s certainly is one of the most important quests in particle physics in the last 30 years. Huge progress was achieved since then at the SPS, RHIC and LHC in understanding the properties of matter at extreme conditions Braun-Munzinger:2015hba , in particular that produced in high energy proton-nucleus and nucleus-nucleus collisions. Yet, crucial aspects of the resulting system remain unclear. Besides, the mechanism of quarkonium production and suppression in Brambilla:2010cs ; Lansberg:2006dh , and Andronic:2015wma collisions is still far from understood despite significant efforts from the theoretical and experimental sides.
One of the biggest challenge is to find a good baseline to properly use heavy quarkonia to diagnose the QGP. This baseline should allow us to correct the yields for the effects characteristic of heavy-quarkonium production and evolution in hadronic matter when QGP is absent. Such effects, a priori measurable in proton-nucleus collisions, include (i) the modification of the nuclear parton densities, commonly known as shadowing, anti-shadowing, EMC effect or Fermi motion – depending on the probed kinematics; (ii) the multiple scattering of partons or of the heavy-quark pair in the nucleus before or after the hard scattering, which leads to an energy loss or the break-up of the formed quarkonium state; (iii) the interaction with other particles produced in the collision – the so-called comovers.
The interplay between these effects and the quarkonium suppression in the QGP is still under debate Andronic:2015wma . Only high precision data on different quarkonium states, in different and systems and in different kinematics regions (, and ) would clear up these debates.
Another source of complexity indeed results from a possible new charmonium-production mechanism suggested to be at work at high c.m.s. energy. Indeed, and quarks are so abundantly produced that they could (re)combine into charmonia. In other words, a deconfined nuclear state could lead to more produced in some kinematical domains.
2.2 The case for AFTER@LHC
In this context, it is important to keep in mind the central role played by fixed-target experiments in this field. In general, they provide very high precision measurements whose relevance is clear from the above discussion. The series of experiments NA38, NA50 and NA60 at the CERN SPS program have indeed shown the first tangible signs of anomalous charmonium suppression in heavy-ion collisions Abreu:1997jh . Yet, as of today, this anomalous character could not unambiguously be related to a specific phenomenon or ensemble of phenomena.
Without any surprise, AFTER@LHC, as it happened with these SPS experiments, present significant assets with respect to the present collider experiments, LHC and RHIC. Its c.m.s. energy of 72 GeV allows one to study quarkonium observables in the region of a possible phase transition where charmonium recombination is not expected to play an important role. High luminosities, inherent to the fixed-target mode, are not only necessary to increase the statistical precision but also to provide us with a better control of acceptances and efficiencies. The NA60 experiment together with experiments at HERA and FNAL provide ample proof to that.
Moreover, the NA61 experiment at SPS energies and the RHIC beam-energy-scan program are the natural complement of AFTER@LHC. Together they can provide a much awaited detailed picture of the phase transition region from SPS to RHIC energies and clarify its nature.
Taking into consideration the AFTER@LHC c.m.s. energy and luminosity, the first inclusive high-precision program for quarkonia is at reach. It is also complemented by heavy- and light-flavour studies as we detail now.
2.2.1 The AFTER@LHC quarkonium-physics case
Quarkonium excited states.
The scope of AFTER@LHC for quarkonium excited states is fully inclusive and encompasses not only the study of , and \mbox{\Upsilon}(nS) states but also the states and the corresponding associated-production channels whose study requires high luminosities, good detector performances (granularity, resolution, PID), wide acceptance and reliable efficiency controls.
It is a known fact that the and states receive significant feed-down contributions. The observed ground states and \mbox{\Upsilon}(1S) are in both cases significantly produced through feed down from higher excited states. This is also true for \mbox{\Upsilon}(2S) and \mbox{\Upsilon}(3S). Understanding quarkonium formation and suppression will most probably only be achieved once the excited-state effect are properly accounted for. This begins with careful studies of the , \mbox{\Upsilon}(2S) and \mbox{\Upsilon}(3S) yields which is clearly possible with AFTER@LHC.
In section 3, we present our projection for the study of and states in collisions at AFTER@LHC energies. These complement a previous study Massacrier:2015qba for and collisions.
Yet, a full description of the mechanisms for the formation of the quarkonium states cannot be dissociated from that of states. Direct measurements are probably possible with AFTER@LHC with a reduced combinatorial background thanks to the reduced energies compared to the LHC or RHIC. The study of topological properties of the decay distribution such as polarisation Faccioli:2012kp ; Faccioli:2008dx ; Faccioli:2014cqa discussed below can complement such direct studies without requiring to detect photons.
Quarkonium polarisation.
Quarkonium polarisation – despite being considered for a long time to be a smoking-gun signal – has so far eluded a thorough explanation despite decades of experimental and theoretical efforts. In the recent years, it was realised that both polar and azimuthal distributions should be measured altogether Faccioli:2008dx . AFTER@LHC, being a high precision experiment, can play in this respect a very important role in clarifying the evolution of the angular anisotropy parameters for the decay of heavy quarkonia from low to high energies and thus enlighten the production mechanism for these states.
Moreover, as mentioned above, the and production in nucleus-nucleus collisions necessarily impacts the interpretation of the and suppression as a probe of QGP formation. The observation of the and suppression patterns in nucleus-nucleus collisions can help to confirm or falsify the sequential quarkonium melting scenario and, therefore, discriminate between the QGP interpretation and other options.
In addition to a direct observation of the and signals in their radiative decays to and which is admittedly challenging in heavy-ion collisions, it may be possible to determine the relative yield of and states by only performing dilepton polarisation measurements. The predicted differences in the nuclear dissociation patterns of and states – owing to their different binding energies– should result in a change of the observed and polarisations from proton-proton to central nucleus-nucleus collisions. This may provide an indication for quarkonium sequential suppression in the QGP. For such observations to be conclusive, the impact of other QGP- and nucleus-related phenomena on the and polarisation should carefully be investigated alongside in collisions or on the states receiving few or no feed down, like the for instance.
2.2.2 Open charm and beauty production
The study of mesons is a natural continuation of the study of and formation and dissociation. mesons has been proposed Satz:2005hx as a baseline to understand QGP effects on quarkonium. In particular, it was proposed to measure their yield in and collisions versus rapidity to get some insights on CNM affecting heavy-quark production, which can be common to quarkonium production. Such studies can be complemented with the measurements of azimuthal anisotropies as a function of rapidity to constraint transport properties of the QGP such as the shear viscosity and the heavy-quark diffusion coefficients. The measurement of the nuclear-modification factor, for open-heavy flavour hadrons can also provide insights into the nature of the heavy-quark interaction with the surrounding matter. The data collected at RHIC at the energy of = 200 GeV Adamczyk:2014uip ; Adare:2006nq show a strong suppression at transverse momenta larger than 3 GeV/c. These results can be considered as the evidence of heavy-energy loss in the hot and dense QGP. However, the details of the interactions with the QGP are not well understood; the major difficulty is to determine the role of the gluon radiation (radiative energy loss) and of the collisional energy loss (due to collisions with other objects in the QGP). A few different models with distinctive assumptions describe reasonably well the current data Adamczyk:2014uip . New observables –preferably for open charm and beauty separately– are thus required to constrain these theoretical approaches.
AFTER@LHC will be capable of delivering high-quality data to address this issue. The measurements of nuclear modification factors for D mesons as a function of transverse momentum and rapidity with AFTER@LHC will provide extremely valuable inputs since the radiative and collisional energy loss are expected to exhibit different and dependence Kikola:2015lka . These high statistics data will also facilitate studies of heavy-flavour azimuthal correlations () which will give another handle on heavy-quark in-medium interactions. These will indeed allow for precise determinations of the transport properties of the QGP, including the energy-loss transport coefficient and the charm-quark-diffusion coefficients. Finally, the measurements of the nuclear-modification factor for -hadrons in interactions will be measurable with a good precision Kikola:2015lka ; Massacrier:2015qba . -hadron measurements via non-prompt in Pb collisions will also be at reach with a slightly reduced statistics compared to the prompt ones but with a smaller background.
2.2.3 Drell-Yan
The hadro-production of a large invariant-mass lepton pair, known as the Drell-Yan (DY) process, is a clean, precise and controllable probe of the short-distance dynamics and the partonic structure of hadrons. In particular, the DY process on nuclear targets is an ideal tool to quantify the initial-state effects taking place during such proton-nucleus reactions. Let us cite the nuclear modification of the parton densities or coherence effects. Indeed, the production of the lepton pair by an electroweak gauge boson does not suffer from final-state interactions, typically associated with the phenomena of energy loss or absorption .
Beside being a typical background for the study of quarkonia, it is also an important process per se where the factorisation of initial-state nuclear effects can be tested. With a nuclear target, the AFTER@LHC kinematic with a LHCb-like detector cover correspond to the antishadowing region for the parton in the targer () for invariant masses slightly above the charmonium family. Increasing the dilepton mass allows one to probe the EMC region, where the nuclear modification of the PDFs was historically observed for the first time. The use of Pb beam allows one to probe a region where shadowing effects are expected. Overall, the study of , Pb and Pb collisions altogether will offer an unique playground to test the factorisation of nPDF effects.
2.2.4 Soft-probe studies
The aforementioned AFTER@LHC studies on hard probes will be complemented with that of soft probes in heavy-ion collisions. In fact, the bulk of low-momentum particles produced in collisions at RHIC and at the LHC is supposed to originate from a fluid-like medium, whose transverse expansion is largely driven by the density gradients produced in the earliest stage of the collision Ollitrault:1992bk ; Kolb:2000sd . The harmonic analysis of particle production in collisions typically focuses on the Fourier decomposition of the azimuthal angle dependence over a narrow region in pseudorapidity around . The coefficients of this decomposition play an important role in understanding the collective behavior of the system. The study of both even and odd order coefficients Naselsky:2012nw provide tight constraints on the collective dynamics and should be determined for that reason.
This requires, in particular, a baseline comparison of (and ) to existing models in collisions, for two otherwise opposite assumptions - with and without a hydrodynamic phase of the evolution. Figure 1 shows predictions of charged hadrons versus pseudorapidity in semi-central PbPb collisions at GeV for a pure cascade (UrQMD) case and for the case when a hydrodynamic evolution is included (vHLLE+UrQMD) Karpenko:2015xea . AFTER@LHC will be able to continue the existing studies on at energies between the RHIC Beam Energy Scan (BES) and the top RHIC and LHC energies not only with high precision but also in a broad rapidity range. Doing so, it will provide a detailed account of the collective evolution between the two energy ranges, around the expected transition line in the QCD phase diagram, between the existing high precision SPS results and the LHC.
2.2.5 Going further
Other points that will be envisaged include:
During the ramp up in energy of the LHC, AFTER@LHC can perform an energy scan to search for a possible critical point along the phase transition line, along the same lines as the RHIC BES program. Even if this scan may be an experimental challenge owing to the possible beam instabilities during this phase, it can provide relatively high statistics allowing for direct searches of the critical point and to look for phase-transition signatures. As for the soft probes study at 72 GeV, it would be complementary to the RHIC BES II program.
A study of the chiral symmetry restoration à la NA60 can be performed in AFTER@LHC by studying the dimuon mass spectrum in for the , , systems. Different measurements can be performed, namely the change in the branching ratio between leptonic and hadronic channels and the change in the width of the mass peaks. This requires significant statistics, a high resolution and a careful control of the background processes including the charm contribution. A comparison with the STAR observation at RHIC (excess with 2 at \mbox{\sqrt{s_{NN}}}=62 GeV – without charm background subtraction) would be extremely useful. Since the NA60 results Arnaldi:2006jq were obtained with a resolution of 20 MeV at the and peaks and AFTER@LHC studies should aim at such precision to be competitive. The particle identification which could be done with a LHCb-like detector would also considerably simplify the analysis and overcome many of the drawbacks present in NA60.
Overall, AFTER@LHC can considerably enrich the existing studies devoted to the search for signatures of the chiral symmetry restoration which plays a fundamental role in explaining the origin of 95% of the meson and baryon masses. Aided by lattice studies and a considerable wealth of information collected since decades on the fundamental aspects of chiral symmetry in the strong sector Borasoy:2007yi , AFTER@LHC is in a unique position to clarify the role of the chiral-symmetry restoration in the dynamics of meson and baryon interactions and also its importance in the interplay between chiral-symmetry restoration and deconfinement transition Karsch:1994hm .
Finally, AFTER@LHC, due to its extended kinematic range can scan a much larger chemical baryon potential range in the QGP phase diagram compared to measurements at mid-rapidity which are presently available in most of the experimental results.
3 Simulations
In this section, we present feasibility studies for quarkonium and Drell-Yan production via the di-muon decay channels in heavy-ion collisions with a 2.76 TeV per nucleon LHC lead beam on a fixed target ( 72 GeV). In Ref. Massacrier:2015qba , we have presented simulations for collisions at GeV resulting from 7 TeV protons impinging on a hydrogen target with a LHCb-like detector. In the following, we follow a similar procedure which we have updated for the heavy-ion case as described below. In addition, we elaborate a little on the possibilities offered by the ALICE detectors used in the fixed-target mode.
3.1 Simulation framework
The simulations presented in this paper are performed for collisions at 72 GeV. Projections for and collisions at 72 GeV are done by applying a nuclear scaling factor to the simulations, assuming neither nuclear nor isospin effects.
Different di-muon sources are separately studied in order to keep a good control over the input distributions and the normalization of these different sources. As a signal we consider either quarkonia or Drell-Yan (which is a background in the case of quarkonium studies). The background under the signal is twofold: the correlated background – charm and beauty pair production (+ and +), and the combinatorial background, mainly muons from pion and kaon decays.
The quarkonium signals as well as the Drell-Yan and charm– and beauty–pair production ( and ) are simulated with HELAC-Onia Shao:2012iz ; Shao:2015vga which provides Les Houches Event Files Alwall:2006yp as output. These events are then processed with Pythia 8 Sjostrand:2007gs to perform the hadronisation, to account for initial/final-state radiations and to decay the resonances. On this step additional muons from the underlying Pythia event can be produced. We have checked that the contribution of combinations of these muons with muons from initial quarkonia or , pairs is negligibly small and thus is not included in the correlated background.
, , states and the continuum rates are obtained in a data-driven way as in Massacrier:2015qba . Quarkonia processed in Pythia 8 are forced to decay into the dimuon-decay channel and the simulated yields are then weighted by the corresponding cross section multiplied by the Branching Ratio (BR) Beringer:1900zz . The open–beauty simulation is performed with a Leading Order (LO) matrix element normalised to a Next-To-Leading-Order (NLO) factor found to be 1.83 Massacrier:2015qba . The Drell-Yan simulation is performed with the process at LO, with the CTEQ6L1 pdf set. In this case, the factor is fond to be 1.2. In order to decrease the time for the simulation in the di-muon invariant mass () region of interest, the simulation is done with a 1.5 GeV/ requirement.
The combinatorial background is obtained from minimum bias collisions generated with Pythia 8, using the process SoftQCD:nonDiffractive with the MRSTMCal.LHgrid LHAPDF (6.1.4) set Buckley:2014ana . The dominant source of combinatorial opposite-sign di-muon pairs are coming from or decays.
The Pythia 8 output is then processed via a fast simulation framework. This last step is performed to account for realistic detector-resolution and particle-identification performances. The detector response – momentum resolution, identification efficiency and misidentification probability with – is simulated with a detector setup similar to the LHCb detector Alves:2008zz with a pseudorapidity coverage of 2 5. The minimum transverse momentum of single muons is required to be greater than 0.7 GeV/. The considered momentum resolution is taken as 0.5 %, the LHCb reporting: 0.4 (0.6)% for a momentum of 3 (100) GeV/ Archilli:2013npa . The considered single identification efficiency of 98% is an average LHCb efficiency for muons coming from decays, for GeV/ and GeV/ Archilli:2013npa . In the case of muons from hadronic decays, the muon is not considered in the simulation if a / decays before 12 m. This corresponds to a decay before the LHCb calorimeter which can be rejected by the LHCb tracking system. If the is produced beyond 12 m or if a is misidentified with in the muon stations, the following momentum-dependent misidentification probabilities Aaij:2014jba are applied: and , for and K respectively. The di-muon efficiency is taken as: . It is assumed in these simulations that the detector performance does not decrease with the event multiplicity.
Projections for and collisions are done by applying a nuclear scaling factor, accounting for the number of binary collisions, to the cross sections obtained from the simulations. Then, the signal and background sources are normalised to the desire integrated luminosity according to the production cross section of the process (taking into account the initial phase space cuts, if any). The values of the cross sections that are obtained from the simulations at 72 GeV are reported in Table. 1, they are integrated over rapidity and . In the case of quarkonium, Drell-Yan, open charm and beauty simulations, the nuclear scaling factor for minimum-bias collisions is and for collisions it is , where and are the beam and projectile mass numbers, respectively. These scaling factors correspond to the absence of hot/cold nuclear and isospin effects. The scaling factor for the combinatorial background is and for and collisions, respectively, where is an average number of binary collisions obtained from Glauber simulations Alver:2008aq ; Loizides:2014vua ; Kikola:2015lka for minimum bias collisions for a given colliding system.
In what follows, we will study the case of Xe gas target, one of the heavier noble gases. The LHC heavy-ion runs are expected to last for about a month per year, which we take as s. In order to provide a baseline at the same energy, 2.76 TeV proton beam is delivered for about a week per year, and our estimated running time for proton-beam collisions on a hydrogen gas target is s. The expected instantaneous and yearly (as defined before) luminosities for and PbXe collisions at 72 GeV are gathered in Table 2. A month of running with the LHC Pb beam on a gas target provides an access to heavy-ion studies with a yearly integrated luminosity as high as 30 nb*-1* with an integrated luminosity of the baseline reaching 250 pb*-1*. Furthermore, studies of collisions at the same centre-of-mass energy of 72 GeV would be of particular interest. The instantaneous luminosity for collisions (with a Xe, Kr, Ar or Ne gas target) is included in Table 2. These instantaneous luminosities are expected to be achieved with a 1-m-long HERMES-like internal gas target Barschel:2015mka . Such an open-storage-cell system with a dedicated pumping system offers a possibility to reach a higher gas pressure111In Massacrier:2015qba , a SMOG-like system was considered and resulted in lower instantaneous luminosities. Given these high instantaneous luminosities, one can obtain the same yields of , , and Drell-Yan in collisions as in the case with relatively short runs. For example, in the case of the Xe system considered in this paper, a running time of 13 hours is enough, which corresponds to an integrated luminosity of 2 . For Kr, Ar and Ne cases that would be 20 hours, 42 hours and 84 hours of running, respectively. These running times correspond to integrated luminosities of 3 , 6.5 and 150 for Kr, Ar and Ne collisions, respectively 222The mentioned instantaneous luminosities correspond to an inelastic rate of 400 MHz. This may however be limited by detector capabilities. In such case, longer running times may be needed in order to achieve the quoted integrated luminosities for systems.. Table 2 also includes average number of binary collisions for minimum-bias events, for each of the considered colliding systems: , Xe and PbXe.
The projections which are presented here for Xe collisions give the highest combinatorial background out of the mentioned Xe, Kr, Ar, Ne cases with the 2.76 TeV proton beam, due to the highest number of binary collisions. The important advantage of the fixed-target mode is the possibility to change relatively easily from one target to another. This would allow us to study different systems, such as Pb, Pb and collisions at the same energy, in different years of LHC running.
3.2 Quarkonium studies at 72 GeV
For the quarkonium PbXe simulations where the full background is considered, the combinatorial one is subtracted assuming the like-sign technique. The background can be reconstructed as an arithmetic or geometric mean of all and pairs produced in an event. That is the worst-case scenario in terms of the expected statistical uncertainties from the combinatorial-background subtraction: the statistical uncertainties of the background scale as the background yield under the signal () – – and add to that of the signal when subtracted. However, this is the most straightforward method of extracting the combinatorial background in an experiment. It does not require an additional normalisation factor and thus does not suffer from additional systematic uncertainties, as the event-mixing technique for example.
As can be seen on Figure 2, high yields are expected in PbXe heavy-ion collisions with an integrated luminosity of 30 . The plot shows the di-muon invariant mass distribution in the mass range after the combinatorial background subtraction with a like-sign technique assuming 1. Thanks to the large yields and large enough signal / background ratios (), each of the state – , and – can easily be reconstructed. The dashed lines represent expected width of states assuming the LHCb-like single muon resolution: the states can clearly be separated. The other background sources, namely correlated and Drell-Yan pairs, are negligibly small in the invariant-mass range. The expected yields for , and and ratios are summarised in Table 3 for , Xe and PbXe collisions. They are integrated over the and a rapidity range of 3 5, accessible for the production. The and rapidity spectra, for each state are shown in Figure 3.
With the integrated luminosity of 30 , extremely large and yields are expected. The and , , and rapidity, , spectra are shown in Figure 3, assuming 1. and signals can be studied with a good precision over a wide transverse momentum range, up to 12 GeV/, and in the rapidity range limited only by the detector acceptance.
The left panel of Figure 4 shows the di-muon invariant-mass distribution in the and mass region and in the entire rapidity range of 2 5. As in the case, the combinatorial background is subtracted using the like-sign technique. In addition, the contributions from correlated , and Drell-Yan backgrounds are included. The correlated background dominates over the other background sources, under the and signal peaks.
The right panel of Figure 4 shows our predicted statistical precision for the nuclear-modification factors () for (red markers) and (blue markers) in PbXe (filled symbols) and Xe (open symbols) collisions at = 72 GeV. The statistical uncertainties include that of the baseline uncertainties of the and expected yields and of the like-sign combinatorial background subtraction - which is the dominant source here. can be extracted as a function of rapidity, in three ranges, with a very good precision.
Such measurements of quarkonia in wide rapidity ranges will allow one to probe the hot and dense medium resulting from these heavy-ion collisions at different temperatures. Furthermore, to gain even more information on the medium properties, such high and expected yields would be used for analyses of differential in centrality, as well as measurements of the elliptic flow, , which we leave for a future publication.
3.3 Drell-Yan pair production in heavy-ion collisions at 72 GeV
The physics program of AFTER also includes the precise measurement of the Drell-Yan pair-production process which can probe initial-state effects on quarks using different Pb collision species.
A great challenge in the Drell-Yan measurement at high energies is the large correlated background from + and +. This background is smaller at the c.m.s. energy of AFTER@LHC in Pb collisions, and can easily be further reduced with a modern vertex detector by applying a secondary vertex cuts on displaced production of and -mesons. Figure 5 separately shows the expected yields for Drell-Yan, and production in PbXe collisions as a function of the di-muon invariant mass in the range between the and (4 8 GeV/), for different rapidity ranges within 2 5. A single muon cut of 1.2 GeV/ is applied. We also stress here that going to the very backward rapidity region (2 3) – which is probably the most interesting one physics wise – renders the Drell-Yan signal even cleaner since the quark-induced processes are favoured and the background is reduced accordingly.
Due to the continuous nature of the signal di-muon invariant-mass distribution shape and to the reduced ratios, the like-sign technique does not provide enough statistical precision for the signal extraction in this case. However, one can precisely determine the high combinatorial background expected in Pb collisions with the event-mixing technique333where an arbitrarily large number of and from different events can be mixed and paired. This technique combined with a large amount of like-sign di-muon background pairs (used to normalise the shape determined by the mixed events) should provide a robust subtraction of the combinatorial background even for small values. These expected Drell-Yan yields and ratios are gathered in Table 4. In such case, the main remaining uncertainty would be of systematical origin, thus difficult to simulate. As such, the combinatorial background is not accounted for in the distributions shown in Figure 5, which should therefore be considered as idealised. On the way, we note that the event-mixing technique can also be applied to improve the extraction of the quarkonium-signal mentioned above.
We note that the also improves with increasing and when going to the most backward rapidity region. This is clear from Table 4 and Figure 6, where the signal over the combinatorial-background ratio is studied as a function of the di-muon invariant mass for the integrated rapidity region of 2 5, and for the following 3 ranges: 2 3, 3 4 and 4 5. As can be also seen on Figure 6, the S/B can be increase by increasing the cut on a single muon. In order to maximise the , the single muon minimum cut is changed from 0.7 GeV/ to 1.2 GeV/. Indeed the particles that contribute to the combinatorial background are mostly produced with lower . However, a further increase of the minimum muon starts to impact the Drell-Yan signal.
3.4 Acceptances of the ALICE detectors in the fixed-target mode
The ALICE experiment is designed to cope with the high multiplicity of heavy-ion collisions and a fixed-target mode can also be implemented for this experiment. The ALICE detector consists in a central barrel part, which measures hadrons, electrons, and photons and in a forward muon spectrometer Aamodt:2008zz ; Abelev:2014ffa . The ALICE upgrade is planned after the Long Shutdown 2 (LS2) of 2018 Abelevetal:2014cna and is designed for recording data with a rate of 200 kHz in and collisions and 50 kHz in PbPb collisions.
The central barrel covers the pseudorapidity range around mid-rapidity, while the muon spectrometer covers the forward pseudorapidity range . It should be noticed that, due to the absorber system the only available region for the installation of a fixed target in ALICE is the A-side, opposite to the muon spectrometer one.
3.4.1 Muon Spectrometer
In a fixed-target mode, the acceptance of the muon spectrometer allows for measurements in the rapidity regions of -2.3<\mbox{y_{\rm c.m.s.}}<-0.8 with a TeV proton beam and of -1.7<\mbox{y_{\rm c.m.s.}}<-0.2 with a A.TeV Pb beam, when considering the target at the nominal interaction point of the experiment.
The vertex can in principle be displaced upstream (in the A-side of ALICE) from the nominal interaction point by up to a few meters. The geometrical acceptance of the muon spectrometer as a function of the longitudinal position of the fixed target, is shown in Figure 7 separately for the tracking system of the muon spectrometer (identified as MUON, corresponding to the tracking stations installed after the hadron absorber) and the MFT (Muon Forward Tracker) which will serve as vertex tracker for the muon spectrometer after LS2 CERN-LHCC-2015-001 .
As one can notice from Figure 7, the geometrical acceptances of the spectrometer and the MFT decouple very early as the target is moved away from the nominal interaction point. For this reason, an optimised inner-tracking system would be needed, providing both vertexing capabilities and the proxy for the tracking in the spectrometer.
On the other hand, the geometrical acceptance of the spectrometer stays reasonably close to \mbox{y_{\rm c.m.s.}}\simeq 0, keeping a range of about one unit of rapidity even for longitudinal positions of the fixed target very far from the nominal interaction point.
It should be however noted that large displacements between the fixed-target position and the hadron absorber would automatically result in an increase of the combinatorial background from semi-muonic decays of pions and kaons. This aspect should be carefully accounted for when evaluating the physics performance of the muon spectrometer in a fixed-target configuration with displaced vertex.
3.4.2 Central Barrel
If the target is located at the nominal interaction point of the experiment, the central barrel covers the far backward region with \mbox{y_{\rm c.m.s.}}<-3.9 and \mbox{y_{\rm c.m.s.}}<-3.3 with a TeV proton beam and a A.TeV Pb beam, respectively. The detection and the identification of particles with the central barrel therefore seems very appealing to access far negative , down to the target fragmentation region. Several questions however need to be addressed in the case of a target far from the nominal interaction point. The tracking performance of the Time Projection Chamber, Time Of Flight and Transition Radiation Detector need to be evaluated for particles entering and traversing their volume with an angular direction different from the nominal specifications.
As in the muon spectrometer case, a new vertex detector would be needed if the current ALICE Inner Tracking System (ITS) falls outside the geometrical acceptance of interest for particle production in a fixed-target configuration. It would provide both vertexing capabilities and the proxy for the tracking in the main detectors of the central barrel.
4 Conclusions
We have outlined the broad spectrum of possibilities offered by AFTER@LHC for studies of the nuclear matter properties in and heavy-ion collisions in an energy range where a phase transition to a partonic matter in heavy-ion collisions is expected.
The heavy-ion studies with the LHC beams used in the fixed-target mode include:
- •
the correlation of the medium temperature and the quarkonium suppression as a function of their binding energies by looking for signs of a sequential suppression with as many states as possible;
- •
the precise and complete determination of the Fourier expansion of identified particle azimuthal asymmetries as a function of rapidity, down to the target fragmentation region;
- •
the determination of different mixtures of quark-gluon plasma and cold nuclear matter effects on the rapidity- and transverse-momentum-dependent particle-yield modifications, in particular heavy-flavoured particles;
- •
the verification/falsification of the factorisation of the initial-state nuclear effects on Drell-Yan-pair measurements on different systems including asymmetric nucleus-nucleus collisions;
- •
the search for net baryon singularities from a possible QCD critical point using the beam-energy-ramp phase;
- •
the search for modifications in the mass width and branching ratios between hadronic and leptonic channels of , and decays as a signature of the chiral-symmetry restoration in the same spirit of the NA60 studies.
As we have shown for a LHCb-like detector, various charmonium and bottomonium states can be measured with an unprecedented accuracy over wide and rapidity ranges both in and collisions. These studies will be accompanied by open heavy-flavour, light-flavour, and Drell-Yan measurements to deliver in-depth characteristics of the partonic matter and the phase transition. Based on our first acceptance studies, the ALICE detectors used in the fixed-target mode can provide complementary measurements in this new energy domain.
Acknowledgements
This research was supported in part by the French P2IO Excellence Laboratory and the French CNRS via the grant FCPPL-Quarkonium4AFTER, and Ministerio de Ciencia e Innovacion of Spain and Xunta de Galicia.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1(1) E 598 Collaboration, J. J. Aubert et al. , “Experimental Observation of a Heavy Particle J,” Phys. Rev. Lett. 33 (1974) 1404–1406 . · doi ↗
- 2(2) S. W. Herb et al. , “Observation of a Dimuon Resonance at 9.5-Ge V in 400-Ge V Proton-Nucleus Collisions,” Phys. Rev. Lett. 39 (1977) 252–255 . · doi ↗
- 3(3) T. A. Armstrong et al. , “Observation of the p wave singlet state of charmonium,” Phys. Rev. Lett. 69 (1992) 2337–2340 . · doi ↗
- 4(4) NA 50 Collaboration, M. C. Abreu et al. , “Evidence for deconfinement of quarks and gluons from the J / psi suppression pattern measured in Pb + Pb collisions at the CERN SPS,” Phys. Lett. B 477 (2000) 28–36 . · doi ↗
- 5(5) P. Hoyer, M. Vanttinen, and U. Sukhatme, “Violation of factorization in charm hadroproduction,” Phys. Lett. B 246 (1990) 217–220 . · doi ↗
- 6(6) NA 3 Collaboration, J. Badier et al. , “Evidence for ψ ψ 𝜓 𝜓 \psi\psi Production in π − superscript 𝜋 \pi^{-} Interactions at 150-Ge V/ c 𝑐 c and 280-Ge V/ c 𝑐 c ,” Phys. Lett. B 114 (1982) 457–460 . · doi ↗
- 7(7) NA 51 Collaboration, A. Baldit et al. , “Study of the isospin symmetry breaking the in the light quark sea of the nucleon from the Drell-Yan process,” Phys. Lett. B 332 (1994) 244–250 . · doi ↗
- 8(8) Nu Sea Collaboration, E. A. Hawker et al. , “Measurement of the light anti-quark flavor asymmetry in the nucleon sea,” Phys. Rev. Lett. 80 (1998) 3715–3718 , ar Xiv:hep-ex/9803011 [hep-ex] . · doi ↗
