PANDA Phase One

The Facility for Antiproton and Ion Research (FAIR) in Darmstadt, Germany, provides unique possibilities for a new generation of hadron-, nuclear- and atomic physics experiments. The future antiProton ANnihilations at DArmstadt (PANDA or $\overline{\rm P}$ANDA) experiment at FAIR will offer a broad physics programme, covering different aspects of the strong interaction. Understanding the latter in the non-perturbative regime remains one of the greatest challenges in contemporary physics. The antiproton-nucleon interaction studied with PANDA provides crucial tests in this area. Furthermore, the high-intensity, low-energy domain of PANDA allows for searches for physics beyond the Standard Model, e.g. through high precision symmetry tests. This paper takes into account a staged approach for the detector setup and for the delivered luminosity from the accelerator. The available detector setup at the time of the delivery of the first antiproton beams in the HESR storage ring is referred to as the \textit{Phase One} setup. The physics programme that is achievable during Phase One is outlined in this paper.


Introduction
The Standard Model (SM) of particle physics has to date successfully described elementary particles and their interactions. However, many challenging questions are yet to be resolved. Some of these are being studied at the high energy frontier at e.g. the LHC at CERN. A different approach is the high precision/high intensity frontier provided by exclusive measurements of hadronic reactions at intermediate energies. This will be exploited in the upcoming PANDA experiment at FAIR, where antiprotonproton and antiproton-nucleus interactions serve as diagnostic tools. The PANDA physics programme consists of four main physics domains: a) Nucleon structure b) Strangeness physics c) Charm and exotics and d) Hadrons in nuclei, as illustrated in Fig. 1.
The theory describing the strongly interacting quarks and gluons is Quantum Chromodynamics (QCD) [1]. At high energies, or short distances, the strong coupling α s is sufficiently weak to enable a perturbative treatment i.e. pQCD. Quarks act as free particles due to asymptotic freedom, an inherent property of QCD [2], and the predictions from pQCD have been rigorously and successfully tested in experiments [3]. At low and intermediate energies, α s increases and pQCD breaks down. The strongly interacting quarks and gluons are confined into hadrons within a radius of ≈1 fm. A quantitative description of the strong interaction at the scale where quarks and gluons form hadrons and up to the onset of pQCD, belongs to the most challenging questions in contemporary physics. This manifests itself in the nucleon whose inherently nonperturbative properties such as the spin [4,5] and, partly  controversial, the mass [6] remain objects of intense discussions and research. Understanding the former requires detailed knowledge about the distribution and motion of the quarks and gluons inside the hadrons. These can be quantified by e.g. electromagnetic structure observables such as form factors and parton distributions.
The mass of a purely light-quark system such as the nucleon, is to a very large extent generated dynamically by the strong interaction via the QCD intrinsically generated scale Λ QCD , rather than the Higgs mechanism. Nature is close to the chiral-limit case of massless quarks. Explaining the mass of nucleons and other hadrons requires a detailed theoretical understanding of the low-energy aspects of QCD, which goes hand in hand with the experimental determination of the hadronic excitation spectrum. In particular, it is illuminating to study hadrons whose building blocks have different masses, from the massless gluons on one hand, to heavy quarks, e.g. charm, on the other. In the latter case the charm quark(s) can be viewed as a near static color source(s) surrounded by the strongly interacting light quarks, a scenario that leads to additional, new forms of matter.
Glueballs, suggested by theorists since more than 40 years [7], constitute one extreme since they consist solely of massless gluons. Hence, 100% of the glueball mass is dynamically generated by the strong interaction. However, unambiguous evidence for their existence has not yet been found. Also the mass of hybrids [8], states which consist of massive quarks and massless gluons, is expected to be largely generated dynamically.
The other extreme are "pure" quark systems containing heavier quarks, e.g. strange or charm. The experimentally well-established hyperons are baryons just like the nucleons, but contain one or several heavier quarks. Strange systems provide a bridge between nucleons, composed of essentially relativistic and non-perturbative quarks, on one side, and the fairly non-relativistic systems containing heavy charm or beauty quarks on the other. With a strong coupling at the charm scale of α s ≈ 0.3, the validity of perturbative QCD becomes questionable. In general, we note that the limit of applicability of perturbative QCD is by itself a complex subject. However, it is a reasonable approximation to describe states and processes in terms of quark and gluon degrees of freedom. Meson-like systems with hidden charm (cc) show interesting features; in particular the experimentally discovered states above the open-charm production threshold that do not fit into the naive conventional quark-antiquark picture and, thereby, likely must have a more complicated structure [9][10][11]. These states are often referred to as the XYZ states, whereby Z refers to isospin one cases, Y to those states that are formed in e + e − annihilations with spinparity of J P C =1 −− , and X to the other remaining resonances. Even for the conventional nucleon, the existence of nonperturbative intrinsic charm, a first-principle property of QCD, has been proposed [12][13][14], but experimentally not firmly confirmed.
At the next level of complexity, where nucleons form nuclei, a long-standing question is how the nuclear force emerges from QCD. The short-distance structure of nuclei, studied in hadronic interactions with atomic nuclei, can shed light on this issue. At high energies, the strong interaction is predicted to be reduced due to colour transparency [16]. At low energies, hadrons are implanted in the nuclear environment and form bound systems with finite life-time. Those could be hypernuclei where one (or several) nucleon(s) in a nucleus is replaced by a hyperon. Studies of hypernuclei could shed light on the longstanding hyperon puzzle of neutron stars. Here, hyperonnucleon and hyperon-hyperon interactions give rise to hyperon pairing which can suppress the cooling of neutron stars [17].
Finally, the validity and limitations of the SM itself remain an open question at the most fundamental level. One example is the matter-antimatter asymmetry, or baryon asymmetry, of the Universe, that cannot be explained within the SM. Unless fine-tuned in the Big Bang, the baryon asymmetry should be of dynamical origin, referred to as baryogenesis [19]. This would however require e.g. CP violating processes to an extent that so far have not been observed experimentally.
To summarise, despite the many successes of the SM, many unresolved puzzles remain. Various efforts from both, theoretical and experimental frontiers are in progress or planned in the near future to address these puzzles [18]. In this paper, we highlight PANDA, a future facility that will exploit the annihilation of antiprotons with protons and nuclei to shed light on the mysteries behind the fundamental forces in nature. PANDA has unique features associated with the usage of antiprotons and a versatile detector that provide a complementary discovery potential with respect to other facilities with the capability to carry out precision studies in the field of particle, hadron, and nuclear physics. In this paper, we outline these features and give an overview of the PANDA physics objectives with emphasis on the programme foreseen for the first phase of operation of PANDA, in the following referred to as Phase One. The structure of the paper is as follows. First, we elaborate on the advantages of antiprotons as a probe. Next, we give a detailed presentation of the PANDA experiment in general and the Phase One conditions in particular. We go through each one of the PANDA physics sections and discuss their underlying purpose and aims, the present experimental status and the potential for PANDA Phase One. Finally, we conclude each part by providing a discussion on its impact and long-term perspectives in which we also briefly outline additional follow-up aspects for the subsequent phases of PANDA.

Opportunities with antiprotons
The intense and precise antiproton beam foreseen in PANDA has many advantages: -The cross sections of hadronic interactions are generally large. -Individual meson-like states can be produced in formation without severe limitations in spin and parity combinations. -Baryons with various flavour, spin and parity can be produced in two-body reactions. -The annihilation process could proceed via gluons and, in that case, providing a gluon-rich environment.
In the following, we elaborate on these points in more detail.
The cross sections associated with antiproton-proton annihilations are generally several orders of magnitude larger than those of experiments using electromagnetic probes. This enables excellent statistical precision already at the moderate luminosities available in Phase One (∼10 31 cm −2 s −1 ). In particular, hadrons composed of strange quarks and hadrons with explicit gluon degreesof-freedom are abundantly produced as demonstrated at a multitude of previous experiments at LEAR, CERN [20].
Hadronic reactions can be divided into two classes: formation and production. In formation, the initial systems fuse into one single state. The line shape of such a state can be determined from the initial system, using a technique called resonance energy scan. The beam momentum is changed in small steps thereby varying the centre-of-mass energy in the mass region of the state of interest and the production rate is measured. Each resulting data point is a convolution of the beam profile and the resonance cross section according to Fig. 2. The true energy-dependent cross section (green dashed line) is determined by the effectively measured cross section (solid blue line) based on the measured yields (markers) and the beam momentum spread (red dotted line). The smaller the momentum spread of the beam, the more precise the measurement of the resonant line shape will be. In formation, the possible quantum numbers of the formed state depend on the probes. In e + e − annihilations, processes in which the formed state has the same quantum numbers as the photon, i.e. J P C = 1 −− , are strongly favoured. States with any other quantum number are strongly suppressed and these therefore have to be produced together with a system of recoiling particles, i.e. in production, or from decays of the 1 −− state. The disadvantage of production with recoils is that the state of interest needs to be identified by the decay products. As a consequence, the mass resolution is limited by the detector resolution, which is typically several orders of magnitude worse than the beam momentum spread. In antiproton-proton annihilations, any state withqq-like, or non-exotic, quantum numbers can be created in formation. With a cooled antiproton beam, like the one foreseen for PANDA, the centre-of-mass energy resolution is excellent reaching a precision that is expected to be about 85 keV at a centreof-mass energy of 3.8 GeV at Phase One. Experiments of this kind are therefore uniquely suited for precision studies of masses, widths and line-shapes of meson-like states with non-exotic quantum numbers that are different from 1 −− . A prominent example of this is the hidden-charm χ c1 (3872) 1 with J P C = 1 ++ , that we will discuss further in Section 6.2.2. Furthermore, PANDA is unique in its capability to probe resonances with high spin. These are difficult to produce using electromagnetic probes, as well as in decays of e.g. B mesons.
Baryons and antibaryons can be produced in twobody reactionspp →B 1 B 2 . The final state baryons can carry strangeness or charm provided theB 1 B 2 system is flavour neutral. In particular for multi-strange hyperons, this is an advantage compared to meson or photon probes, where strangeness conservation requires that the hyperon is produced with the corresponding number of associated kaons. As a result, the final state comprises at least three pseudo-stable particles, which complicates the partial-wave analysis necessary in hyperon spectroscopy. Two-body reactions on the other hand, in particular close to the kinematic threshold, typically involve few partial waves. Furthermore, spin observables and decay parameters can be accessed in a straight-forward way in twobody reactions. This enables production dynamics studies as well as charge conjugation parity (CP) symmetry tests in the strange sector. The particle-antiparticle symmetric final state minimises systematic uncertainties. In principle, the aforementioned advantages apply also for baryon-antibaryon production in e + e − colliders. However, the typically much smaller cross sections result in low production rates. The resulting data samples are therefore smaller and in order to obtain sufficiently many events, methods such as missing kinematics or single-tag analysis are common. This however limits the possibility to reduce the background and achieve good resolution. Inpp experiments, one can obtain large data samples also in exclusive analysis, which increases the discovery potential.

Resonance Cross Section
Thepp → X process includes quark-antiquark annihilations, which result in gluons. Therefore, antiprotonproton annihilation provides a gluon-rich environment, where states with a gluonic component are likely to be produced if they exist. Gluon-rich environments exist also in radiative decays of charmonia and in central hadronhadron collisions. However, in radiative decays, reconstruction of the properties of the resonant state of interest relies solely on detector information since the pro-cess is not a formation process. As a result, the resolution is limited by the detector. The same is true for central hadron-hadron collisions, where the final state consists of the scattered hadrons and the produced resonance. The spin and parity of the resulting multi-particle final state is complicated to reconstruct without assumptions about the underlying production mechanism. This in turn leads to model-dependent ambiguities. The processpp → X, where X refers to a single resonance, is less complicated in this regard.
The momentum range (1.5−15 GeV/c), momentum precision (∆p/p≈10 −4 −10 −5 ) and intensity (∼10 31 cm −2 s −1 ) of the antiproton beam in PANDA is tailored for strong interaction studies. PANDA will give access to the mass regime up to 5.5 GeV/c 2 whereby recently new and interesting forms of hadronic matter have been observed (XY Z states), it can study the hadron-antihadron formation close to their production threshold, such as for the reactionspp → ΛΛ, ΣΣ, ΞΞ, ΩΩ, DD, D sDs , Λ cΛc , Σ cΣc , Ξ cΞc , Ω cΩc , and it has the resolution to measure the line-shape of states very accurately.

The PANDA experiment at FAIR
The PANDA experiment is one of the four pillars of FAIR [21]. PANDA will be a fixed-target experiment where the antiproton beam will impinge on a cluster jet (pp orpA) or pellet target (pp) or target foils (pA). The High Energy Storage Ring (HESR) [22] can provide antiprotons with momenta from 1.5 GeV/c up to 15 GeV/c. The physics goals of PANDA outlined in this paper require a detector system with nearly full solid-angle coverage, high-resolution tracking, calorimetry and particle identification over a broad momentum range as well as vertex reconstruction.
The success of the physics programme will depend not only on the detector performance but also on the quality and intensity of the antiproton beam. Antiprotons are produced from reactions of 30 GeV/c protons on a nickel or copper target. The source of these protons will be a dedicated high-power proton Linac followed by the existing SIS18 synchrotron and the new SIS100 synchrotron. Produced antiprotons are focused by a pulsed magnetic horn and selected in a magnetic channel at a momentum of around 3.7 GeV/c. After phase-space cooling in the Collector Ring (CR), packets of about 10 8 antiprotons are transferred to the HESR for accumulation and subsequent acceleration or deceleration necessary for measurements in PANDA. In this mode of operation, the HESR is able to accumulate up to 10 10 antiprotons from 100 injections within a time span of 1000 s. In a later stage of FAIR, the accumulation will take place in a dedicated ring, i.e. the Recuperated Experimental Storage Ring (RESR), allowing for up to 10 11 antiprotons to be injected and stored in the HESR. An important feature of the HESR is the versatile stochastic cooling system operating during accumulation and target operation. It is designed to deliver a relative beam-momentum spread (∆p/p) of better than 5 · 10 −5 . Furthermore, it includes a barrier bucket cavity that compensates for the mean energy loss in the thick target and that fine-tunes the absolute beam energy. This enables precise energy scans around hadronic resonances and kinematic thresholds. The centre-of-mass resolution will be about 50 keV, which to date is unreachable by other accelerators using different probes.

Staging of the experiment
The PANDA experiment will follow a staged approach in the construction of the detector and in the usage of the antiproton beam. It comprises four phases, briefly outlined below.
The first phase, Phase Zero, started in 2018 and it refers to physics activities where PANDA detectors and analysis methods are used at existing and running facilities. One example is the usage of PANDA tracking stations in the upgraded HADES at GSI [23], another is the deployment of parts of the PANDA calorimeter for experiments with A1 at MAMI [24].
The installation of the first major detector components of PANDA, including the two spectrometer magnets, will follow Phase Zero. This installation phase will be completed with a commissioning of the detectors using a proton beam at the HESR. The start of Phase One will be marked with the usage of antiprotons together with the commissioned detectors. The corresponding physics programme is outlined in this paper. During Phase One, the HESR will be capable of accumulating at most 10 10 antiprotons in 1000 s. The luminosity is expected to rise gradually from about 10 30 cm −2 s −1 to the maximum of 2 × 10 31 cm −2 s −1 (at 15 GeV/c) during Phase One. The cluster jet target available from the beginning allows the use of hydrogen, deuterium or noble gases as target materials. The available PANDA detector of Phase One will be referred to as the start setup and includes most of the major components as shown in Fig. 3. A description of the various available detector components will be given in section 3.2. The total integrated luminosity for Phase One is expected to be about 0.5 fb −1 .
The detector will be completed according to the final design in Phase Two. The main components beyond the start setup are the detector for charged particle identification in the forward region (e.g. forward RICH detector) and the completion of the Gas Electron Multiplier (GEM) tracker and forward trackers. Moreover, a pellet target system will become available. The corresponding setup will be referred to as the full setup. In Phase Three, the RESR will be available at FAIR which provides an increase in luminosity at HESR by a factor of approximately 20.

The PANDA start setup
To achieve the full physics potential of PANDA, the complete set of detector systems is needed. In Phase One, not all of these will be available and the focus is therefore on reactions with large expected cross sections and good signal-to-background ratios as well as relatively small multiplicities of final-state particles.
In this section, we primarily describe the hardware systems to be installed as part of the start setup. The PANDA detector consists of two main parts: -The Target Spectrometer (TS) for the detection of particles at large scattering angles (> 10 • in the horizontal direction and > 5 • in the vertical direction). The momentum measurement of charged particles is based on a superconducting solenoid magnet with a field strength of 2 T. -The Forward Spectrometer (FS) for particles emitted in the forward direction (±5 • vertically and ±10 • horizontally). The momentum measurement is based on a dipole magnet with a bending power of up to 2 Tm.
The magnet system is described in Ref. [25]. Both spectrometers are integrated with devices to perform tasks such as high resolution tracking, particle identification (PID), calorimetry and muon detection. The internal target operation of PANDA will employ a cluster jet target that can be operated with hydrogen as well as heavier gases. With hydrogen, an average luminosity of 10 31 cm −2 s −1 can be reached in the experiment [26].  Figure 3: Schematic overview of the start setup of PANDA. The various tracking detectors are indicated in red, the components for particle identification in blue, the electromagnetic calorimeters in green, and the target system in black.

The Target Spectrometer
The beam-target interaction point will be enclosed by the Micro Vertex Detector (MVD) that will measure the interaction vertex position. It will consist of hybrid silicon pixels and silicon strip sensors. The vertex resolution is designed to be about 35 µm in the transverse direction and 100 µm in the longitudinal direction. Moreover, the MVD significantly contributes to the reconstruction of the transverse momentum of charged tracks [27]. The Straw Tube Tracker (STT) will surround the MVD with the primary purpose of measuring the momenta of particles from the curvature of their trajectories in the solenoid field. The low-mass (1.2% of a radiation length) STT detector will consist of gas-filled straw-tubes arranged in cylindrical layers parallel to the beam direction. From these straws, a resolution better than 150 µm in the transverse x and y coordinates can be achieved. Some straw tube layers will be skewed with respect to the beam direction which enables an estimation of the z coordinate along to the beam. The z resolution will be approximately 3 mm. The STT will also contribute to the charged particle identification by measuring the energy loss dE/dx. Details of the STT can be found in Ref. [28]. The PANDA Barrel DIRC [29], surrounding the STT, will cover the polar angle region between 22 • and 140 • . The DIRC will be surrounded by a barrel-shaped Time of Flight (TOF) detector consisting of scintillating tiles read out by silicon photomultipliers. The expected time resolution [30], better than 100 ps, will allow precise particle identification of low momentum particles. In addition it is used to assign a precise time information to reconstructed tracks to combine several tracks to one event of the same antiproton-target interaction. The high time resolution and the long mean time between events during the phase-one operation of PANDA will ensure that the mixing of events can be suppressed to a negligible amount. The electromagnetic calorimeter (EMC), that will measure the energies of charged and neutral particles, will consist of three main parts: The barrel, the forward endcap and the backward endcap. The expected high count rates and the geometrically compact design of the target spectrometer require a fast scintillator material with a short radiation length and small Molière radius. Lead-tungstate (PbWO 4 ) fulfils the demands for photons, electrons and hadrons in the energy range of PANDA. The signals from the lead-tungstate crystals are read out by large-area avalanche photodiodes, except in the central part of the forward endcap where vacuum photo-tetrodes are needed for the expected higher rates. The EMC also plays an important role in the particle identification. In particular for electron/positron identification, it can sup-press background from charged pions with a factor of about 1000 for momenta above 0.5 GeV/c. A detailed description of the detector system can be found in Ref. [31]. The laminated yoke of the solenoid magnet, outside the barrel EMC, is interleaved with sensitive layers to act as a range system for the detection and identification of muons. Rectangular aluminium Mini Drift Tubes (MDT) are foreseen as sensors between the absorber layers. Details of this system are described in Ref. [32]. Downstream of the target, within the TS, a system of GEM foils will be located.
The GEM planes will offer tracking of particles emitted with polar angles below 22 • , a region that the STT in the target spectrometer will not cover. In the start setup, two GEM stations will be installed. Part of the particles that pass the GEM tracking detector will be further registered by the Forward Spectrometer (FS) rather than the TS.

The Forward Spectrometer
The FS detector systems are conceptually similar to those of the TS, but will have a planar geometry instead of a cylindrical one. The detector planes will be arranged perpendicular to the beam pipe and thereby measure the deflection of particle trajectories in the field of the dipole magnet. Downstream of the GEMs, two pairs of straw tube tracking stations are foreseen for the start setup [33]. One will be placed in front of the dipole magnet and the other inside its field. Particle identification will be provided by the Forward TOF wall consisting of scintillating slabs. The signals from the latter will be read out by photomultiplier tubes offering a time resolution better than 100 ps [34]. Forward-going photons and electrons will be detected and identified by a Shashlyk-type calorimeter with high resolution and efficiency. The detection is based on lead-scintillator sandwiches read out with wave-length shifting fibres passing through the block and coupled to photomultiplier tubes. The system is described in detail in Ref. [35]. At the end of the FS, a muon range system is placed using sensors interleaved with absorber layers similar to the TS.

Luminosity determination
The luminosity at PANDA will be determined by using elastic antiproton-proton scattering as the reference channel. Since the Coulomb part of the elastic scattering can be calculated precisely and dominates at small momentum transfers, the polar angle of 3-8 mrad is chosen for the measurement. The track of each scattered antiproton and therefore the angular distribution of the tracks will be measured by the luminosity detector made of four layers of thin monolithic silicon pixel sensors (HV-MAPS) [36]. An absolute precision of 5% for the time integrated luminosity is expected and a relative precision of 1% during the energy scans.

Data acquisition
The PANDA data acquisition concept is being developed to match the complexity of a next-generation hadron physics experiment. It will make use of high-level software algorithms for the on-line selection of events within the continuous data stream. This so-called software-based trigger system replaces the more traditional hardwaredriven trigger systems that have been a common standard in the past. In order to handle the expected Phase One event rate of 2 MHz, every subdetector system is a selftriggering entity. Signals are detected autonomously by the sub-systems and are pre-processed in order to transmit only the physically relevant information. The online event selection occurs in computing nodes, which first perform event-building followed by filtering of physical signatures of interest for the corresponding beam-target settings. This concept provides a high degree of flexibility in the choice of trigger algorithms and hence a more sophisticated event selection based on complex trigger conditions, compared to the standard approach of hardwarebased triggers.

The simulation and analysis framework
The feasibility studies presented in this paper have been carried out using a common simulation and analysis framework named PandaROOT [37]. This framework provides a complete simulation chain starting from the Monte Carlo event generation, followed by particle propagation and detector response, signal digitisation, reconstruction and calibration, and finally the physics analysis. PandaROOT is derived from the FairROOT framework [38] which is based on ROOT [39]. FairROOT offers a large set of base classes which enables a straight-forward customisation for each individual detector setup. It offers an input-output manager, a run manager, database handling, an event display and the Virtual Monte Carlo (VMC) interface which allows to select different simulation engines. In addition, it uses the task system of ROOT to combine and exchange different algorithms into a simulation chain.
The first part in the simulation chain is the event generation. Here, the initial interaction of the antiproton beam with the target material is simulated using a Monte Carlo approach. Different generators exist for different purposes. Dedicated reactions and their subsequent decays are generated by the standard signal generator Evt-Gen [40]. For the generic background, the Dual Parton Model (DPM) [41] and the Fritiof (FTF) model [42] can be chosen. Both include all possible final states and are tuned to an exhaustive compilation of experimental data. For detector-and software performance studies, the BoxGenerator creates single types of particles within user-defined momentum and angular ranges.
The generated particles are propagated through a detailed detector model, simulating the reactions with the detector material and possible decays in-flight. For this purpose, Geant3 and Geant4 are available to the user.
The level of detail in the virtual detector description varies between the different subdetectors but all active components, as well as most of the passive material, are included. Separate descriptions are prepared for the start setup and the full setup. From this stage, the energy deposit, the position and the time of a given interaction in a sensitive detector element is delivered as output, all with infinite resolution. Real data will however consist of electronic signals with finite spatial-and time resolution. Therefore, the digitisation converts the information from the particle propagation stage into signals that mimic those of a real experiment. This includes noise and effects from discriminators and electronics. For some detector systems, the final electronics is not yet defined. In those cases, the digitisation procedure is based on realistic assumptions.
In the reconstruction, the signals from the digitisation stage are combined into tracks. The procedure is divided into two steps: a local and a global part. In the local part, detector signals in a given tracking subdetector are combined into tracklets. Furthermore, the signal information is converted back to physical quantities such as position, energy deposit and time. In the global reconstruction, the tracklets from different tracking detectors are combined into tracks. Different algorithms are applied in the barrel part and the forward part. The track finding is followed by track fitting using a Kalman filter, where effects from different particle species and materials are taken into account. PANDA simulations thereby achieve a momentum resolution of about 1%.
At the particle identification stage, the information from the dedicated PID detectors and the EMCs are associated with a charged track based on the distance between the predicted flight path and the hit position in the detector. Hits in the EMC without a corresponding charged track are regarded as neutral particles. The probabilities for various particle types of the different subdetectors are then combined into an overall probability of a given particle species.
The selection of events for partial or complete reaction channels, referred to as Physics Analysis, is performed based on the combined tracking, PID and calorimetry data using the Rho package, an integrated part of Panda-ROOT. With Rho, various constrained fits such as vertex fits, mass fits and tree fits are available.

Nucleon structure
Hadron structure observables provide a way to test QCD and phenomenological approaches to the strong interaction in the confinement domain. Electromagnetic probes are particularly convenient and have been used extensively over the past 60 years. The structure is parameterised in terms of observables like form factors or structure functions.
Electromagnetic form factors (EMFFs) quantify the hadron structure as a function of the four-momentum transfer squared q 2 . At low energies, they probe distances of about the size of a hadron. EMFFs are defined on the whole q 2 complex plane and for q 2 < 0, they are referred to as space-like and for q 2 > 0 as time-like. Space-like EMFFs are real functions of q 2 and can be studied in elastic electron-hadron scattering. Assuming one-photon exchange (OPE) being the dominant process, protons and other spin-1/2 particles are described by two EMFFs: the electric G E (q 2 ) and the magnetic G M (q 2 ) form factor. In the so-called Breit frame, these are the Fourier transforms of the charge and magnetisation density, respectively. Time-like EMFFs are complex and can be studied using different processes in different q 2 regions. In the following, we consider baryons, denoted B, B 1 and B 2 . For unstable baryons, the low- p with m l the mass of the lepton l) can be probed by the reaction pp → + − π 0 . For all types of baryons, the high-q 2 region At high energies, corresponding to distances much smaller than the size of a hadron, individual building blocks are resolved rather than the hadron as a whole. Here, the factorisation theorem applies, stating that the interaction can be factorised into a hard, reaction-specific but perturbative and hence calculable part and a soft, reaction-universal and measurable part. In the spacelike region, probed by deep inelastic lepton-hadron scattering, the structure is described by parton distribution functions (PDFs) [43], generalised parton distributions (GPDs) [44][45][46][47][48][49][50], transverse momentum dependent parton distribution functions (TMDs) [51], and transition distribution amplitudes (TDAs) [52,53]. These nonperturbative objects are complementary tools to explore the structure of the nucleon at the partonic level. They extend the information given by the EMFFs and provide more detailed descriptions of the spatial and momentum distributions of the constituent partons and the spin structure. In the time-like region, they can be accessed experimentally in hard hadron-antihadron annihilations. Detailed studies to access πN TDAs at PANDA in the reactionspp → γπ 0 → e + e − π 0 andpp → J/Ψ π 0 → e + e − π 0 have been presented in Refs. [54,55]. For these measurements, as well as for the TMD studies, the designed high luminosity of PANDA is needed to accumulate reasonable statistics. The counterparts of the GPDs in the annihilation processes are the generalized distribution amplitudes (GDAs). They can be measured in the hard exclusive processespp → γγ [56] andpp → γM [57,58], where M could be a pseudo-scalar or vector meson (e.g. π 0 , η, ρ 0 , φ). Cross section measurements as a function of s and t that allow us to test the theoretical models are expected to start in Phase Two of PANDA. In the following, we focus on EMFFs, in line with the emphasis of Phase One.

State of the art
Elastic electron-proton scattering has been studied since the 1950s [59]. During the first decades, unpolarised electron-nucleon scattering was analysed using the Rosenbluth separation method [60]. Modern facilities, offering high-intensity lepton beams and high-resolution detectors, gave rise to a renewed interest in the field [61,62]. In particular, the polarisation transfer method [63] applied by the JLab-GEp collaboration (see [62] and references therein) revealed the surprising result that the ratio G E /G M decreases almost linearly with Q 2 = −q 2 . This result is in contrast to previous measurements based on a Rosenbluth method using unpolarised elastic ep scattering which do not reveal such a dependency. The correction of the unpolarised elastic ep cross section by the two-photon exchange (TPE) contribution has been suggested to solve this discrepancy [64]. The TPE correction does not impact the polarization transfer extraction of G E /G M in a significant way. The large amount of high-quality data inspired extensive activity also on the theory side, from which we have learned about the importance of vector dominance at low q 2 [65,66].
Until recently, measurements in the time-like region have not achieved precisions comparable to the corresponding space-like data, partly because most e + e − colliders have been optimised in different q 2 regions [67,68]. Inpp annihilation experiments, the clean identification of e + e − pairs has been a challenge. Among the few experiments that so far have provided a separation between G E and G M of the proton, the results at overlapping energies disagree. The ratio R = |G E |/|G M |, accessible from the final state angular distribution, has been measured below q 2 = 9 (GeV/c) 2 by PS170 at LEAR [69], BABAR [70] and more recently by BESIII [71][72][73] and CMD-3 [74]. The PS170 and BABAR data differ up to 3σ, while the BESIII and CMD-3 measurements have large total uncertainties. In the limit |q 2 | → ∞, the space-like and the time-like form factors should approach the same value as a consequence of the Phragmén-Lindelöf theorem [75]. Experimentally, the onset of this scale has not been established (see Ref. [68] for a recent review). In measurements just below |q 2 | = 20.25 (GeV/c) 2 , the time-like magnetic form factor is about two times larger than the corresponding space-like one. A recent analysis of BaBar data above |q 2 | = 20.25 (GeV/c) 2 , indicates a decreasing difference, but the uncertainties are large [70].
In 2019, the BESIII collaboration measured the Born cross section of the process e + e − →pp and the proton EMFFs at 22 centre-of-mass energy points from q 2 = 4 (GeV/c) 2 to q 2 = 9.5 (GeV/c) 2 with an improved accuracy [76], comparable to data in the space-like region. Uncertainties on the form factor ratio |G E |/|G M | better than 10% have been achieved at different q 2 values below 5 (GeV/c) 2 . The BESIII data on the proton effective form factor confirm the structures seen by the BABAR Collaboration. These structures are currently the subject of several theoretical studies [77][78][79].
The PANDA experiment aims to improve the current situation of the time-like EMFFs by providing data in a large kinematic region between 5.08 (GeV/c) 2 and ∼ 30 (GeV/c) 2 . Precisions in this region of at least a factor 3 better than the current data, as well as measurements in the unphysical region below (2M p ) 2 , whereby M p is the mass of the proton, are called for to constrain the theoretical models and to resolve the aforementioned issues.

Potential of Phase One
The PANDA experiment in Phase One offers the opportunity to measure the proton form factor in the process pp → + − , ( = e, µ) over a wide energy range, including the high |q 2 | region [80,81]. Thepp → µ + µ − reactions can be studied for the first time. The interest forpp annihilation into heavy leptons (µ and τ ) has been discussed in several theory studies [82][83][84]. In contrast to thē pp → e + e − process, thepp → µ + µ − reaction has the advantage that corrections due to final state radiation are expected to be smaller. Measuring both channels should therefore allow the formalism for radiative corrections to be tested.
Furthermore, the unphysical region of the proton EMFFs can be accessed through the measurement of thē pp → + − π 0 process [85][86][87]. These measurements by PANDA are unique and will provide the possibility to test models for this process that contain EMFFs [88].

EMFFs inpp → e + e −
A previous simulation study of the processpp → e + e − within the PandaROOT framework demonstrates the excellent prospect of nucleon structure studies with the PANDA design luminosity [80]. The simulations were performed applying an integrated luminosity of 2 fb −1 for each energy-scan point and the full PANDA setup corresponding to the conditions that are foreseen for Phase Three of PANDA. A new, dedicated simulation study with the Phase One conditions has recently been performed at q 2 = 5.08 and 8.21 (GeV/c) 2 (p lab = 1.5 and 3.3 GeV/c, respectively). The difficulty of the measurement is related to the hadronic background, mostly annihilation with the subsequent production of two charged pions. This reaction has a cross section about five to six orders of magnitude larger than that of the production of a lepton pair. In the energy scale of the PANDA experiment, the mass of the electron is sufficiently close to the pion mass for this to be an issue. Therefore, the signal and the main background reactions have very similar kinematics. The signal events are generated according to the differential cross section parameterised in terms of proton EMFFs from Ref. [89] with setting R = |G E |/|G M | = 1. The same event selection criteria as in Ref. [80] were applied. The output of the PID and tracking subdetectors as EMC, STT, MVD, and barrel DIRC have been used to separate the signal from the background. These resulted in signal efficiencies of 40% at p lab = 1.5 GeV/c and 44% at p lab = 3.3 GeV/c. The suppression factor of the main background processpp → π + π − was found to be of the order ∼ 10 8 . The proton form factors |G E |, |G M |, and their ratio R = |G E |/|G M | are extracted from the electron angular distribution, after reconstruction and efficiency correction. The proton effective form factor |G e eff | is extracted from the determined cross section of the signal (σ) integrated over the electron polar angle. The resulting precision for different q 2 is summarised in Table 1 and shown in Fig. 5, together with existing experimental data. Systematic uncertainties arise due background contamination and uncertainties in the luminosity measurement. These effects can be quantified by MC simulations. From these we conclude that the proton EMFFs can be measured with an overall good precision and accuracy. At low q 2 , the signal event yield is relatively large. However, at higher q 2 , the cross section of the process reduces significantly which leads to a smaller event yield and thus larger statistical uncertainties for a given integrated luminosity. Previous studies show that the efficiency at larger q 2 is sufficient for precise cross section measurements [80].

EMFFs inpp
An independent Monte Carlo simulation study of thē pp → µ + µ − reaction has been carried out at q 2 = 5.08 (GeV/c) 2 . The di-muon channel provides a clean environment, where radiative corrections from final state photon emissions are reduced thanks to the larger mass of the muon. However in case of muons, the suppression of the hadronic backgroundpp → π + π − is more challenging. Muon identification is mainly based on the information from the Muon System, since other subdetectors show less separation power which complicates the background separation considerably. Monte Carlo samples of 10 8 events were generated for the background processpp → π + π − . They were used for the determination of the background suppression factor and for the calculation of the pion contamination, which will remain in the signal events after the application of all selection criteria. The separation of the signal from the background has been optimised through the use of multivariate classification methods (Boosted Decision Trees). The event selection is described in Ref. [81]. A background rejection factor of 1.2 × 10 −5 was achieved, resulting in a signal-to-background ratio of 1:8. The total signal efficiency is 31.5%. Due to the insufficient background rejection, the pion contamination needs to be subtracted from the signal and the corresponding angular distributions by Monte Carlo modelling and subsequent subtraction. This has been taken into account in our feasibility studies. The angular distributions from the pion contamination are reconstructed with both the expected magnitude and shape. The sensitivity of the EMFFs to the shape was investigated and from that, the systematic uncertainty was estimated. The ratio R, and consequently |G E | and |G M |, were extracted from the angular distribution of the muons after background subtraction and efficiency correction. The results are summarised in Fig. 5 and Table 1. The uncertainty of the signal cross section is dominated by the luminosity uncertainty. The simultaneous but independent measurements of the effective EMFFs G e eff and G µ eff from the e + e − final state and the µ + µ − final state, respectively, enable a test of radia- Table 1: Results from simulation studies ofpp → e + e − andpp → µ + µ − . tive corrections that are applied in the e + e − channel. The expected uncertainty in the ratio G e eff /G µ eff is estimated to be 3.2% already during Phase One. It should be noted that although the uncertainties from radiative corrections are not yet taken into account, these are expected to contribute with only a small fraction to the total uncertainty.

EMFFs inpp
Some information about the unphysical region can be obtained from thepp → e + e − π 0 process, when studied in different intervals of the pion angular distribution. In the time-like region, the EMFFs are complex, hence they have a relative phase. This phase is generally inaccessible for protons in an experiment with an unpolarised beam or target. However, the cross section ofpp → e + e − π 0 channel can provide some information, as outlined in Refs. [85,87].
The validity of the theoretical models used to describe the cross section of the processpp → e + e − π 0 needs to be tested experimentally. Since PANDA has almost 4π coverage, the measurement of the final state angular distributions in the processespp → e + e − π 0 andpp → γπ 0 will provide a sensitive check of these models. The EMFFs extracted at threshold viapp → e + e − π 0 andpp → e + e − or e + e − →pp can be compared and used as an additional test. We note that the process γp → pe + e − may give access to the EMFFs of the proton in the unphysical region as well. The feasibility to extract EMFFs from such challenging measurements has so-far not been demonstrated as discussed in theoretical studies in Refs. [98,99].
For an ideal detector (100% acceptance and efficiency) and an integrated luminosity of 0.1 fb −1 , the expected count rate for this reaction for q 2 < 2 (GeV/c) 2 has been found to be up to 10 5 events in different intervals of the pion angular distribution [54,88]. This number is about a factor two larger than the corresponding value forpp → e + e − at q 2 = 5.08 (GeV/c) 2 . The large expected count rate ofpp → e + e − π 0 and the clean separation between this channel and background [88], indicate good prospects for EMFF measurements in the unphysical region already in PANDA Phase One. Full simulation studies to investigate the possibility to extract the proton EMFFs in this region at PANDA are currently being carried out.

Impact and long-term perspective
The simulation studies presented in the previous sections show that PANDA will improve the precision of the proton EMFF measurements for q 2 > 5.08 (GeV/c) 2 . The measurement of the effective form factor can be extended to higher |q 2 | values in Phase Three of the experiment when the luminosity reaches its design value. This enables systematic comparisons of space-like and time-like EMFFs at large |q 2 | and hence, the onset of the convergence scale of the space-like and time-like form factors may be deduced. Furthermore, the foreseen PANDA studies of thepp → µ + µ − are unique. Since the effects from final state radiation are negligible for muons, this channel provides an important cross check of thepp → e + e − results. Finally, in PANDA, the unphysical region of the proton EMFF will be accessed for the first time through thepp → e + e − π 0 process.
In general, the relative phase between the electric and the magnetic form factors is inaccessible in unpolarised cross section measurements. To measure the phase, either a polarised antiproton beam and/or a polarised proton target is required. The feasibility of implementing a transversely polarised proton target in the PANDA detector is under investigation. If feasible, the PANDA experiment will offer a first direct measurement of the relative phase between G E and G M .
In the next stages of PANDA, we will be able to extend our studies to the hard exclusive processes and perform complementary measurements of the nucleon structure objects.

Physics with strangeness
The key question in hyperon physics is "What happens if you replace one (or several) light quark(s) in the nucleon with one (or several) heavier one(s)?". Strangeness serves as a diagnostic tool for various phenomena in subatomic physics: 1. Hyperons provide a new angle to the structure and excitations of the nucleon, since the strange quark is sufficiently light to relate the knowledge about hyperons to nucleons and vice versa. 2. Hyperon decays, where the spin is experimentally accessible, provide an ideal testing ground for CP violation and thereby searches for physics beyond the SM at the precision frontier. Furthermore, it can give clues about Baryogenesis [19]. 3. In hypernuclei, strangeness provides an additional degree of freedom which plays a key role in understanding e.g. neutron stars [100]. 4. Enhancement of strangeness in relativistic heavy-ion collisions was one of the first proposed signals of Quark-Gluon Plasma [101].  Figure 5: Expected total precisions, indicated by the red and blue error bars, on the determination of (a) the proton form factor ratio, (b) the proton effective form factor, (c) the proton electric form factor, and (d) the proton magnetic form factor, from the present simulations for PANDA Phase One as a function of q 2 . Also shown are data from PS170 [69], BaBar [70,90], BESIII [71][72][73]76], CMD-3 [74], E835 [91], Fenice [92], E760 [93], DM1 [94], DM2 [95], CLEO [96], and ADONE73 [97]. The results indicated by PANDA MC study are based on an integrated luminosity of 0.1 fb −1 for each point in q 2 . Number 1 will be explored with PANDA Phase One within the subtopics hyperon production and hyperon spectroscopy. Number 2, i.e. hyperon decays will be studied extensively in Phases Two and Three. However, a good understanding of the production mechanism has proven crucial to decay measurements [102] and the planned hyperon production studies within Phase One are therefore an important milestone in the search for CP violation in baryon decays. Number 3 will be investigated during Phases Two and Three within our programme for hadrons in nuclei. Number 4 is currently studied at ALICE [103] and is not within the scope of PANDA. However, preci-sion studies of strangeness production in elementarypp reactions contribute to a more general understanding of strangeness production, which can be useful also in more complex reactions at higher energies. The same is true for the planned studies of hyperon-antihyperon pair production inpN reactions. These will provide information on absorption and rescattering of hyperons as well as antihyperons under well-defined conditions in cold nuclei. In this chapter, we discuss the subtopics hyperon production and hyperon spectroscopy in the context of what can be achieved at Phase One. In section 7.1 we also discuss antistrange hadrons in nuclei. Production plane of the pp → ΛΛ reaction. The y-axis of the Λ decay frame is perpendicular to the production plane. The z-axis is in the direction of the outgoing Λ with respect to origin in the centre-of-mass frame.

Hyperon production
The scale probed in a hadronic reaction is influenced by the energy and, therefore, by the mass of the produced quarks. The strange quark mass is m s ≈ 100 MeV/c 2 which corresponds to the scale where quarks and gluons form hadrons. Therefore, the relevant degrees of freedom are unclear -quarks and gluons, or hadrons? It is challenging to solve QCD in this energy regime. Guidance by experimental data is needed to improve the theory such that quantitative predictions can be accurately tested. As an intermediate step phenomenological models are developed which are constrained by experimental data. Exclusive hyperon-antihyperon production provides the cleanest environment to constrain these models. Phenomenological models based on quark-gluon degrees of freedom [104], meson exchange [105] and a combination of the two [106] have been developed for single-strange hyperons. The quarkgluon approach and the meson exchange approach have also been extended to the multi-strange sector [107][108][109].
Here, the interaction requires either annihilation of two quark-antiquark pairs, or in the meson picture, exchange of two kaons. This means that the interactions occur at shorter distances which make double-strange production more suitable for establishing the relevant degrees of freedom. The clearest difference between the quark-gluon picture and the kaon exchange picture is typically found in the predictions of spin observables e.g. polarisation and spin correlations. Understanding the mechanism of hyperon production is also important in order to correctly interpret experimental data on other aspects of hyperons. One example is that of the recent theoretical and experimental studies of the hyperon structure in e + e − → ΛΛ. In Ref. [110], the time-like form factors G E and G M were predicted, including their relative phase ∆Φ = Φ(G E ) − Φ(G M ) that manifests itself in a polarised final state. Different potential models were applied, usingpp →ΛΛ data from PS185 [111] as input. In the model predictions for the channel e + e − → ΛΛ, the total cross section and the form factor ratio R = |G E /G M | differ very little for different potentials. However, the relative phase ∆Φ and hence the Λ polarisa-tion showed large sensitivity. New data from BESIII [112] provide an independent test of the ΛΛ potentials. Another example is hyperons and antihyperons in atomic nuclei, since understanding the elementarypp →Ȳ Y reactions is crucial in order to correctly interpret data frompA collisions.
Spin observables are straight-forward to measure for ground-state hyperons thanks to their weak, self-analysing decays. This means that the decay products are preferentially emitted along the spin direction of the parent hadron. Consider a spin 1 2 hyperon Y decaying into a spin as illustrated in Fig. 6a, where α Y [3] is the asymmetry parameter of the hyperon decay related to the interference between the parity conserving and the parity violating decay amplitudes. The polarisation P y is related to the production dynamics, hence it depends on the centre-of-mass (CMS) energy / beam momentum and on the hyperon scattering angle. In strong production processes, such as pp →Ȳ Y , with unpolarised beam and target, the polarisation can be non-zero normal to the production plane, spanned by the incoming antiproton beam and the outgoing antihyperon as shown in Fig. 6b. Spin correlations between the produced hyperon and antihyperon are also accessible [113] and from these, the singlet fraction can be calculated, i.e. the fraction of the hyperon-antihyperon pairs that are produced in a spin singlet state. Additional information can be obtained from hyperons that decay into other hyperons, e.g. the Ξ. In the sequential decay Ξ − → Λπ − , Λ → pπ − , the additional asymmetry parameters β and γ of the Ξ − hyperon are accessible via the joint angular distribution of the Λ hyperons and the protons [114,115]. For spin 3 2 hyperons, e.g. the Ω − , the spin structure is more complicated. Only considering the polarisation parameters of individual spin 3 2 hyperons, we find that spin 3 2 hyperons produced in strong processes like pp → Ω + Ω − have seven non-zero polarisation param-eters. Three of these can be extracted from the Λ angular distribution in the Ω − → ΛK − decay [116]. The remaining four parameters can be obtained by studying the joint angular distribution I(θ Λ , φ Λ , θ p , φ p ) of the Λ hyperons from the Ω − decay and the protons from the subsequent Λ decay [115].

State of the art
The PS185 collaboration has provided a large set of highquality data on single-strange hyperons [111,117] produced in antiproton-proton annihilation. One interesting finding is that theΛΛ pair is produced almost exclusively in a spin triplet state. This can be explained by the Λ quark structure: the light u and d quarks are in a relative spin-0 state, which means that the spin of the Λ is carried by the s quark. Various theoretical investigations reproduce this finding [104][105][106], but no model has yet been formulated to describe the complete spin structure of the reaction. The extension of models into the double-strange sector [107,108] and even the triple-strange Ω [109], have not been tested due to the lack of data. For Ξ − and Ξ 0 frompp annihilations, only a few bubble-chamber events exist [118], whereas no data at all are available related to triple-strange hyperon production since no studies have been carried out so far. As a result, further progress of this field is still pending. New data on the spin structure of pp → Y Y for ground-state multi-strange and singlecharmed hyperons would therefore be immensely important for the development of a coherent picture of the role of spin in strangeness production.

Potential of Phase One
Previous studies of mainly single-, but also a few double strange hyperon-antihyperon pairs produced in antiproton-proton annihilations show remarkably large cross sections within the PANDA energy range [117]. This means that large hyperon data samples can be collected within a reasonable time even with the reduced luminosity of the Phase One setup. Simulation studies of exclusive hyperon production, using a simplified Monte Carlo framework, were performed and presented in detail in Refs. [116,119,120,122]. New, dedicated simulation studies of hyperon production have been performed for this review: GeV/c and p beam = 6.0 GeV/c.
GeV/c and p beam = 7.0 GeV/c.
The beam momenta for the single-strange hyperons were chosen in order to coincide with those of other benchmark studies. For the double-strange Ξ − , the chosen beam momenta coincide with the hyperon spectroscopy campaign (4.6 GeV/c, see Section 5.2) and the χ c1 (3872) line-shape campaign (7 GeV/c, see Section 6.2.2). In these new simulation studies, a realistic PandaROOT implementation of the Phase One conditions was used, though with some simplifications due to current limitation in the simulation software: i) ideal pattern recognition, with some additional criteria on the number of hits per track in order to mimic a realistic implementation of the track reconstruction ii) ideal PID matching, to reduce the run-time. It was however shown in Ref. [116] that the event selection can be performed without PID thanks to the distinct topology of hyperon events: since the hyperons have relatively long life-time (10 −10 s) they travel a measurable distance before decaying. This provides a challenge in the tracking but also makes the background reduction very efficient. Around 10 6 events were generated forΛΛ andΞ + Ξ − [121,122], whereas 10 4 events forΣ 0 Λ [123]. The larger event samples in theΛΛ andΞ + Ξ − cases enable studies of spin observables. In the case ofΣ 0 Λ, only a general feasibility study of cross section and angular distribution measurements has been carried out so far. TheΛΛ andΣ 0 Λ final states were modelled using parameterisations based on data from Refs. [111,124], where it was found that single-strange antihyperons are very strongly forward-going in the CMS of the reaction. TheΞ + Ξ − final state has never been studied and was therefore generated both with an isotropic angular distribution and with a forward-peaking distribution. The results were found to differ only marginally.
The particles were propagated through the Panda-ROOT detector implementation and the signals were digitised, reconstructed and analysed. The signal events were selected by requiring all stable (p,p and γ) or pseudostable (π + and π − ) particles to be found: -ΛΛ: p, π − ,p and π + .
To reduce the number of background photon signals, additional energy cuts were applied to identify the photon from theΣ 0 decay [123]. The Λ andΛ, that appear in all channels, were identified by combining the reconstructed pions and protons/antiprotons and applying vertex fits and mass window criteria on the combinations. Furthermore, the decay vertex of the Λ/Λ was required to be displaced with a certain distance from the interaction point. To identifyΣ 0 or Ξ − /Ξ + , the Λ/Λ candidates were combined with the photons or remaining pions. In the case ofΛΛ andΣ 0 Λ final states, four-momentum conservation was applied using kinematic fits constraining the momenta and total energy of the final-state particles with the initial antiproton-proton momenta and energy to further reduce the background. Since the Ξ − decays sequentially, a more elaborate method including a decay tree fitter was applied [121,122].
The resulting signal efficiencies are given in Table 2, that also includes the results from the Ξ * study described in Section 5.2.2. The expected rates of reconstructed events are calculated based on the Phase One luminosity of 10 31 cm −2 s −1 and cross sections from Refs. [111,124] (ΛΛ andΣ 0 Λ) and Ref. [109] (Ξ + Ξ − ). The signal-to- background ratios (S/B) were obtained by simulating 10 7 events at each energy, generated with the Dual Parton Model [41].
In this work, we have also investigated the feasibility of reconstructing spin observables such as the polarisation and spin correlations using the methods outlined in Ref. [116]. For the analysis, the pp → ΛΛ,Λ →pπ + , Λ → pπ − sample was used, containing 157000 signal events surviving the selection criteria. A sample of this size can be collected within a few hours with the Phase One luminosity. The simulated events were weighted according to an input polarisation function P y = sin 2θ Λ and the spin correlation distributions C ij = sin θ Λ (i, j = x, y, z). Symmetry implies P Y = -PȲ which means that the extracted polarisation from Λ andΛ can be combined for better statistical precision.
The reconstruction efficiency was accounted for using two different, independent methods: i) regular, multidimensional acceptance correction as in Ref. [120] and ii) using the acceptance-independent method outlined in Ref. [116]. The results of the MC simulations were divided into bins with respect to theΛ scattering angle. In each bin, the polarisation P Y and spin correlations C ij were reconstructed. The resulting polarisation distribution is shown in panel a) of Fig. 7 with acceptance corrections and in panel b) with the acceptance-independent method. The polarisation distributions extracted with the two independent methods agree with each other as well as with the input functions.
In the same way, spin observables of the Ξ − hyperons were studied at both 4.6 GeV/c and 7.0 GeV/c. The number of signal events were 7.2·10 4 and 6.7·10 4 , respectively, samples that can be collected within a few days during Phase One. The resulting polarisation distributions as a function of cos θ Ξ obtained at each energy are shown in Fig. 8. The singlet fractions were calculated from the spin correlations and are shown in Fig. 9. A singlet fraction of 0 means that all Ξ −Ξ + states are produced in a spin triplet state, a fraction of 1 means they are all in a singlet state, and a fraction of 0.25 means the spins are completely uncorrelated. In Ref. [108], the singlet fraction is predicted to be 0 for forward-goingΞ + and closer to 1 in the backward region. This is in contrast to the single-strange case where the singlet fraction is almost independent of the scattering angle [117]. The results of the simulations shown in Fig. 9  indicate that the uncertainties in the singlet fraction will be modest at all scattering angles, which enables a precise test of the prediction from Ref. [108].
Most systematic effects that are important in cross section measurements, e.g. trigger efficiencies and luminosity, are expected to be isotropically distributed in a near 4π experiment like PANDA. This means that their impact on angular distributions, and parameters extracted  from these, are expected to be small. Hyperon polarisation studies with BESIII (e.g. [112]) instead indicate that imperfections in the Monte Carlo description of the data, due to for example gain drift in HV supplies, may be more important. Most of these effects can however only be studied once PANDA is operational and by careful Monte Carlo modelling, they can be minimised. In the simulation studies presented here, three basic consistency tests have been performed in order to reveal eventual sensitivity to detection-and reconstruction artefacts: i) comparison between generated and reconstructed distributions that are efficiency corrected ii) comparison between extracted hyperon and antihyperon parameters iii) comparison between two different efficiency correction methods. All three tests show differences that are negligible with respect to the small statistical uncertainties.

Hyperon spectroscopy
Baryon spectroscopy has been decisive in the development of our understanding of the microscopic world, the best example being the plethora of new states discovered in the 1950's and 1960's. It was found that these states could be organised according to "the Eightfold Way", i.e. SU(3) flavour symmetry, that led to formulation of the quark model by Gell-Mann and Zweig [125]. Though successful in classifying ground-state baryons and describing some of their ground-state properties, the quark model fails to explain some features of the baryon excitation spectra. This indicates that the underlying picture is more complicated. In contemporary baryon spectroscopy, the most intriguing questions are i) Which effective degrees of freedom are adequate to describe the hadronic reaction dynamics? Which baryonic excitations are efficiently and well described in a three-quark picture and which are generated by coupledchannel effects of hadronic interactions? ii) To which extent does the excitation spectra of baryons consisting of u, d, s obey SU(3) flavour symmetry? iii) Are there exotic baryon states, e.g. pentaquarks or dibaryons?
Among the theoretical tools available to study the spectra and internal properties of baryons, lattice QCD approaches have received a lot of attention thanks to the tremendous progress over the past years. Prominent examples are the mass prediction of the double charm ground state Ξ cc baryon [126][127][128][129][130], now confirmed by LHCb [131], and accurate lattice calculations of the mass splitting of the neutron and proton [132].
However, for the excited states in the light-baryon sector differences between recent calculations [133,134] remain and unambiguous conclusions cannot yet be drawn. Other approaches to baryon excitation spectra are based on the Dyson-Schwinger framework [135], on the coupled-channel chiral Lagrangian [136][137][138][139], and on the AdS/QCD approach [140][141][142]. The next step is systematic studies of the strange sector, in particular states with double and triple strangeness. These bridge the gap between the highly relativistic light quarks and the less relativistic heavy ones.

State of the art
So far, worldwide experimental efforts in baryon spectroscopy have been focused on N * and ∆ resonances. Most of the known states have masses smaller than 2 GeV/c 2 and were discovered in πN scattering experiments. In recent years, many laboratories (JLab, ELSA, MAMI, GRAAL, Spring-8 etc) have studied these resonances in photon-induced reactions [143,144]. As a result, the data bank on nucleon and ∆ spectra has become significantly bigger and a lot has been learned. However, there are several puzzles that remain to be resolved.
One example is the so-called missing resonance problem of Constituent Quark Models (CQMs): many states that are expected from these phenomenological-driven models have not been observed experimentally. This is in contrast to the Dyson-Schwinger approach whose predictions agree almost one-to-one with the experimentally measured light baryon spectra below 2 GeV/c 2 [145,146]. This observation demonstrates the shortcomings of CQMs, thereby motivating the necessity to experimentally establish the spectra of excited baryons. For a successful campaign, an experimental approach is needed in which these states are searched for and their properties are studied using various complementary initial probes such as πN , γN , and, with PANDA,pN .
Another example of an unresolved conundrum is the level ordering: The lightest baryon, i.e the nucleon, has J P = 1 2 + and the next-to-lightest baryon is expected to be its parity partner, with J P = 1 2 − . However, this is in contrast to experimental findings where the Roper N * (1440) resonance, with J P = 1 2 + , is significantly lighter than the lightest J P = 1 2 − state, i.e. the N * (1535).
A new angle to the aforementioned puzzles can be provided by studying how they carry over to strange baryons. In the single-strange sector, the missing CQM resonance problem remains. Regarding the level-ordering, the situation is very different regarding light baryons: the parity partner of the lightest Λ hyperon is the Λ(1405) which is indeed the next-to-lightest isosinglet hyperon [147]. However, the Λ(1405) is very light, and, therefore, it has been suggested to be a molecular state, see e.g. Ref. [148,149]. The existing world data on double-and triple-strange baryons are very scarce and do not allow for the kind of systematic comparisons with theory predictions that led to progress in the light and single-strange sector. Only one excited Ξ state and no excited Ω states are considered well established within the PDG classification scheme [3]. It is also worth pointing out that even for the ground state Ξ and Ω, the parity has not been determined experimentally. Furthermore, the spin determination of the Ω is not model-independent but inferred by assumptions on the Ξ c and Ω c spin [150]. It would be very illuminating to study the features of the double-and triple-strange hyperon spectra since it enables a systematic comparison of systems containing different strangeness.

Potential for Phase One
A dedicated simulation study has been performed of thē pp → ΛK − Ξ + + c.c. reaction at a beam momentum of 4.6 GeV/c. In the following, the inclusion of the charge conjugate channel is implicit. In spectroscopy, parameters like mass, widths and Dalitz plots are essential. Therefore, the focus of this study is to estimate how well such parameters can be measured with PANDA. The simulated data sample of 4.5·10 6 events includes the Ξ(1690) ± and Ξ(1820) ± resonances, decaying into ΛK − + c.c. (each 40% of the total generated events), as well as non-resonant ΛK − Ξ + + c.c. production (20% of the generated sample).
The simulated widths of the Ξ(1690) − and Ξ(1820) − resonances were 30 MeV/c 2 and 24 MeV/c 2 , respectively, in line with the PDG [3]. The event generation was performed using EvtGen [151] with the reaction topology as illustrated in Figure 10. The angular distribution of the produced Ξ * resonance are isotropically generated since no information from experimental data exist. Further technical details related to the simulation study can be found in Ref. [152]. The analysis was performed in the same way as described in Section 5.1.2. The final state is required to contain p,p, π − , π + , K − and K + . The Λ candidates were identified by combining p and π − into a common vertex and applying a mass window criterion. The Ξ − (Ξ * ) hyperons were identified by combining Λ candidates with the remaining pions (kaons). Background was further suppressed by a decay tree fit in the same way as in Section 5.1.2. The exclusive reconstruction efficiency was found to be 5.4%. We assume app →ΛKΞ +c.c. cross section of 1 µb, where the production mainly occurs through a Ξ − Ξ * + c.c. pair and where the excited cascade could be either Ξ * (1690) or Ξ * (1820). With this assumption, the reconstruction rate is 0.2 s −1 or 18000 events per day. These cross sections have never been measured and assumed of the same order as the ground-stateΞ + Ξ − [153] that was measured by Ref. [118] to be around 1 µb.
The background was studied using a DPM sample containing 10 8 events and the signal events were weighted assuming a total cross section of 50 mb. No background events survived the selection criteria and we therefore conclude that on a 90% confidence level, the signal-tobackground is S/B > 19. The numbers are summarised in Table 2.
The reconstructed Dalitz plot and ΛK − invariant mass are shown in Figure 11. The reconstruction efficiency distribution is flat with respect to the Dalitz plot variables and the angles. This is a necessary condition in order to minimise systematic effects in the planned partial-wave analysis of this final state.
In order to evaluate the Ξ andΞ resonance parameters, the ΛK − andΛK + mass distributions have been fitted with two Voigt functions combined with a polynomial. By comparing the reconstructed ΛK − andΛK + widths to the generated ones, the mass resolution was estimated to σ M = 4.0 MeV/c 2 for the Ξ(1690) − and σ M = 6.7 MeV/c 2 for the Ξ(1820) − . The obtained fit values are shown in Table 3. In both cases, the fitted masses are in good agreement with the input values.

Impact and long-term perspective
PANDA will be a strangeness factory where many different aspects of hyperon physics can be studied. Double-and triple strange hyperons are unknown territory both when it comes to production dynamics, spin observables and spectroscopy. Long-standing questions, such as relevant degrees of freedom and quark structure, can be investigated already during the first years with reduced detector  setup and luminosity. Furthermore, the measurements in Phase One provide important milestones for the foreseen precision tests of CP conservation, that will be carried out when the design luminosity and the full PANDA setup are available in the subsequent Phase Two and Three. In the latter, copious amounts of weak, two-body hyperon decays will be recorded -several million exclusively reconstructed ΛΛ pairs every hour. This enables precise measurements of the decay asymmetry parameters. In the absence of CP violation, the asymmetry parameters of a hyperon have the same magnitude but the opposite sign of those of the antihyperon, e.g. α = −ᾱ. Differences in the decay asymmetry therefore indicate violation of CP symmetry. Thepp →Ȳ Y reaction provides a clean test of CP violation, since the initial state is a CP eigenstate and no mixing between the baryon and antibaryon is expected to occur. Since hyperons and antihyperons can be produced and detected at the same rate and in very large amounts, the prospects are excellent for ground-breaking symmetry tests that could help us to understand the matterantimatter asymmetry of the Universe. With the PS185 experiment at LEAR, the CP symmetry inpp → ΛΛ was studied and a sensitivity of 10 −2 was reached [111]. This measurement was based upon ∼ 10 5 reconstructed events which statistics can be obtained in Phase One of PANDA within only a few hours of beam time. The perspective of such an asymmetry measurement in Phase Three with a ten times higher luminosity is, therefore, very promising. In addition, Phase Three opens up the possibility to study also single-charm hyperons. A systematic comparison between the strange and the charm sector will be an important step towards a coherent understanding of nonperturbative QCD at different scales.

Charm and exotics
The original constituent quark model (CQM) describes mesons and baryons. In CQM, mesons are described as quark-antiquark states (qq) interacting through a potential. One of the motivations for this description was the non-observation of mesons with strangeness or charge larger than unity, neither had states been observed with other spin and parity combinations than those consistent with fermion-antifermion pairs. However, QCD allows for any colour-neutral combination of strongly interacting quarks and gluons and therefore, CQM-based models can be extended to incorporate the dynamics of glueballs, hybrids and multiquarks. These states are often referred to as QCD exotics.
Glueballs (gg or ggg) are formed due to the selfcoupling of the colour-charged gluons. This feature of the strong interaction is of particular interest since the glueball mass has no contribution from the Higgs mechanism. Instead, it is completely dynamically generated by the strong interaction. Most glueballs predicted by QCD or phenomenological models have the same quantum numbers as mesons and hence they can mix. As a consequence, it is a challenge to unambiguously determine the glueball fraction of an observed hadronic state.
In addition to glueballs, there are meson-or baryonlike states for which QCD admits a gluonic component called hybrids, e.g. qqg. Hybrids can, in addition to the spin-parity combinations allowed for regular hadrons, also have spin-exotic quantum numbers. To establish the existence of hybrids experimentally, the decomposition of quantum numbers requires sophisticated partial-wave analysis (PWA) tools and large data samples.
Also other colourless combinations of multiquark resonances are allowed within QCD. The study of multiquarks has experienced tremendous progress during the last decade. Examples of multiquark states are tetraquarks (qqqq) or pentaquarks (qqqqq). However, many open questions remain, in particular about the internal structure of the observed states. Precision measurements of various resonance properties are needed, as well as ab-initio theoretical predictions, in order to reach deeper insights about the structure of multiquark states.
The search for exotic hadrons is being carried out at several energy scales, from the light u and d scale to the bottom quark scale. A fundamental question to be answered concerns the relevant degrees of freedom -should excited light hadrons be described in terms of quarks and gluons, or are various dynamical effects, e.g. at meson pair thresholds, more important? In the light quark sector, many resonances are broad and overlap in mass. This means that they mix if they have the same quantum numbers. The advantage of the light sector is that the production cross sections are generally large, allowing for large data samples to be collected within a short time. This is an advantage when determining spin and parity through partial-wave analyses.
The physics of hidden-charm states, such as charmonium, is expected to be very different due to the higher mass of the charm quark (m c 1.2 GeV/c 2 > Λ QCD ). The strong coupling constant in this region is α s ≈ 0.3, corresponding to an energy scale below the region in which perturbation theory starts to break down. At these energies, quark and gluon degrees of freedom become relevant. The velocity of the charm quark is relatively small, (v/c) 2 ∼0.3. Systems with charm can be partly described in a non-relativistic framework with relativistic effects added perturbatively, such as spin-spin and spin-orbit coupling [155]. One of the interesting questions is how large are the relativistic corrections actually. The structure of a separated energy scale (m c m c v/c m c (v/c) 2 ) makes heavy-quark systems, such as charmonium, ideal probes to study the transition between perturbative and nonperturbative regimes [156,157].
Meson-like systems composed of heavy and light constituent quarks, such as open-charm states, are complementary to that of hidden-charm meson-like states. Also here, various striking experimental observations have been made in the past [159,160] pointing to the possible existence of narrow resonances that do not fit into the conventional heavy-light meson pattern. A recent example is the intriguing observation of LHCb, speculating the existence of an open-charm tetraquark with a mass around 2.9 GeV/c 2 [154]. Besides spectroscopy aspects, groundstate open-charm states decay weakly, providing access to, e.g., semi-leptonic form factors. The field of open-charm spectroscopy and electro-weak processes will become accessible in the later stages of PANDA, beyond that of Phase One. Its success depends on the completion of the vertex reconstruction capabilities of PANDA and higher luminosities for excellent statistical significance. Differential cross section measurements will be accessible in Phase One, which allows for unique studies of the production mechanism of pairs of open-charm mesons and baryons in antiproton-proton collisions. Such measurements have the potential to study the intrinsic charm content of the nucleon and, thereby, shed light on the recently predicted nonvanishing asymmetric charm-anticharm sea from lattice QCD [158].
In the following, we discuss the Phase One perspectives of the meson-like spectroscopy programme of PANDA at various mass scales, starting from the light-quark sector to the hidden-charm region.
The lightest tensor (J P C = 2 ++ ) glueball is predicted in the mass range from 2 to 2.5 GeV/c 2 [162]. The possible mixing of two nonets ( 3 P 2 and 3 F 2 ) results in five expected isoscalar states. The JETSET collaboration at LEAR has reported a tensor component in the mass range around 2.2 GeV/c 2 in thepp → φφ reaction [184]. However, due to the limited size of the data sample, no firm conclusions could be drawn.
In the vicinity of meson-pair production thresholds, narrow meson-like excitations can appear. Prominent examples in the light quark sector are the a 0 (980) and the f 0 (980) scalar mesons. These states are strongly attracted by the KK threshold and are believed to have a large KK component. The narrow vector meson φ(2170), discovered by BaBar [185], is particularly interesting in this context. It does not fit into the qq model, it is comparatively narrow (≈ 83 MeV) and the mass is close to the φf 0 (980) threshold. It is debated whether the φ(2170) is an ss tetraquark or hybrid state [186]. Close to the K * K threshold, the COMPASS collaboration discovered a relatively narrow (Γ ≈ 153 MeV) axial-vector meson, the a 1 (1420) [187]. It has been interpreted as the isospin partner of the established f 1 (1420) [188]. The latter can be attributed a molecular-type KKπ component [189], opening up for a possibility that also the a 1 (1420) is a molecular-type state. The first coupled-channel calculation related to a potential axial-vector molecule state originates from [190]. There are further interpretations proposed such a triangle singularity from rescattering of the a 1 (1260) [191]. There has been significant progress in recent years in lattice scattering and coupled channel calculations including in the light meson sector [192,193].

Potential for Phase One
In the search for exotic hadrons, the gluon-rich environment and the access to allqq-like quantum numbers in formation, give PANDA a unique advantage compared to e + e − experiments. Furthermore, states with non-qq quantum numbers can be accessed in production.
The reactionpp → φφ is considered suitable for tensor glueball searches, since the production via intermediate conventional qq states is OZI suppressed in contrast to production via an intermediate glueball. Already during Phase One of PANDA, we will collect data samples of this reaction that are two orders of magnitude larger than achieved by previous experiments. A potential tensor component would reveal itself in energy scans and amplitude analyses.
The f 1 (1420) can be identified through the decay to KKπ and studied at a centre-of-mass energy of about 2.25 GeV inpp → π + π − + KKπ andpp → π 0 /η + KKπ. In the latter cases, the amplitude analysis is simpler since only one recoil (π 0 or η) is involved. The a 1 (1420) can be accessed in 3π combinations from thepp → π + π − π + π − reaction. The COMPASS analysis shows that very large samples are required [187]. Here, PANDA will profit from the large expected production cross sections inpp annihilations. The cross sections for pion modes are in the order of mb, while reactions involving kaons are in the order of 100 µb. The observed intensity of f 1 (1420) in pp → K + K 0 π − π + π − is about 1% [199]. This makes the prospects excellent for studying the f 1 (1420) as well as searching for the a 1 (1420) inpp annihilations already during Phase One of PANDA.
Furthermore, insights on the nature of the φ(2170) will be obtained by studying other production mechanisms and hitherto unmeasured decay patterns. At PANDA, the φ(2170) will be accessible in reactions involving π 0 , η, or π + π − recoils at centre-of-mass energies of about 2.6 GeV. In a similar way, searches for hybrid candidates such as π 1 (1400), π 1 (1600) and π 1 (2015) can be performed.

State of the art
In 2003, the discovery of a signal in the J/ψπ + π − channel near the D 0D * 0 threshold completely changed our understanding of the charmonium spectra [200]. Up to this point, the quark model originally published in 1978 [201] had been very successful in describing all observed states. However, the new signal, established the state denoted χ c1 (3872) or X(3872), turned out to have properties at odds with the quark model. After 2003, many more states in the charmonium and bottomonium mass range were discovered. While all states below the lowest Swave open charm threshold behave in accordance with the quark model, the states above fit neither in mass nor in other properties. This family of exotic states is now referred to as the XY Z states. Arguably the most prominent states, besides the aforementioned χ c1 (3872), are the vector-meson states Y (4260) [202] and Y (4360) [203] as well as the charged states Z(4430) [204], Z c (3900) [205], Z c (4020) [206], Z cs (3985) [207] in the charmonium sector and the charged states Z b (10610) and Z b (10650) [208] in the bottomonium sector. The most prevalent interpretations of these states are hybrid mesons (quark states with an active gluon degree of freedom), compact tetraquarks (bound systems of diquarks and anti-diquarks), hadroquarkonia (a compact heavy quarkonium surrounded by a light quark cloud) and hadronic molecules (bound systems of two mesons; when located very near the relevant S-wave threshold these can be very extended). Recent reviews of various calculations can be found in Refs. [9][10][11]15]. In particular the Z states -charged states decaying into final states that contain both a heavy quark and its antiquark -have received a lot of attention since they must contain at least four quarks [209]. In addition, lattice calculations predict a supermultiplet of hybrid mesons including exotic quantum numbers with a similar pattern in both charmonium and bottomonium [210].
As of today there is no consensus which one of the mentioned models explains the properties of the XY Z states best. Clearly more experimental information is needed to make progress. The two most pressing issues are: -Where are the spin partner states of the observed XY Z states? Their location contains valuable information about the most prominent component of the states, since different assumptions lead to different effects of spin symmetry violation [211]. PANDA is well prepared to hunt for those spin partner states, since the production mechanism is not constrained to certain quantum numbers. -What is the line shape of the near threshold states?
This allows one to especially investigate the role of the two-meson component in a given state, since a strongly coupled continuum necessarily leaves an imprint in the line shapes [15]. Moreover, a virtual state cannot have a prominent compact component [212].
PANDA can provide a significant contribution to answer these questions, in particular the second one, already in Phase One. Precision measurements of line-shape parameters of resonances provide crucial information that sheds light on their internal structure. The determination of these parameters for narrow states is particularly challenging and requires a facility with sufficient resolution to reach the necessary sensitivity.
In the following, we illustrate this by discussing the capability of PANDA to perform resonance energy scans, using the famous χ c1 (3872) state with J P C = 1 ++ as a benchmark. The χ c1 (3872) has a small natural width; until recently the 90% C.L. upper limit was estimated to be 1.2 MeV [213]. A new measurement from the LHCb data are compatible with an absolute Breit-Wigner decay width of Γ = 1.39 ± 0.24 ± 0.10 MeV for the χ c1 (3872) [214]. However, a Flatté-like line shape model where the state is described by a resonance pole with a Full-Width-at-Half-Maximum of about 220 keV is equally probable. The result from LHCb emphasises the need for precision lineshape measurements with significantly better mass resolution than offered by experiments that rely on the detector resolution, typically around a few MeV. Only experiments like PANDA, where these resonances are accessible in formation, offer a direct and thus model-independent measurement of the line-shape.
The analysis presented in the following is meant as a demonstration of the precision capabilities of PANDA, but the technique can be applied to extract key properties of other resonances as well.

Potential for Phase One
PANDA offers a unique possibility to reach sub-MeV resolution exploiting the cooled antiproton beam from the Table 4: Summary of parameter settings in the simulation study [215]. All parameters are defined in the text.  [50,70,100,130,180,250, 500] keV Line-shape parameter E f − [10.0, 9.5, 9.0, 8.8, 8.3, 8.0, 7.5, 7.0] MeV HESR. This has been demonstrated by a feasibility study of the χ c1 (3872) line-shape measurement, to be carried out in a future energy scan designed to precisely measure absolute decay widths and line shapes [215]. The χ c1 (3872), as well as all other non-exotic J P C combinations, can be created in formation inpp annihilation.
The details of the PANDA feasibility study can be found in Ref. [215]. In this paper, we focus on the conditions expected for Phase One. This implies an HESR beam momentum spread (beam energy resolution) of dp/p = 5 · 10 −5 (dE CMS = 83.9 keV) and an integrated luminosity of L = 1170 (day · nb) −1 .
In our study, we quantify i) the sensitivity of an absolute measurement of the natural decay width Γ 0 ii) the capability to distinguish two scenarios: a loosely bound (D 0 -D * 0 ) molecular state and a virtual scattering state.
Both scenarios have been studied under the assumption that PANDA will collect data in 40 energy points during 2x40 days of beam time, i.e. two days per energy point, with the Phase One operation conditions. This is considered a reasonable amount of time to allocate for this kind of measurements, especially since data for other purposes, e.g. hyperon-antihyperon pair production, can be collected in parallel.
The parameter Γ 0 is determined by fitting a Voigt function, i.e. a convolution of a Breit-Wigner with a natural decay width Γ 0 and a Gaussian with a standard deviation σ Beam , accounting for the beam momentum distribution.
The molecular line shape differs significantly from that of a less sophisticated Breit-Wigner-like resonance shape. It depends on the given decay channel (here J/ψπ + π − ) and on the dynamic Flatté parameter E f [216,217] (or the equivalent inverse scattering length, γ in [218]), that model a bound or virtual state.
For each of the six different input signal cross-sections, σ S = (150,100,75,50,30,20) nb, the full procedure of simulation, PDF generation and Breit-Wigner/Molecule line shape fitting has been carried out, employing a maximumlikelihood method.
The resulting sensitivities in terms of the relative uncertainty ∆Γ meas /Γ meas of the measured decay width are summarised for the Breit-Wigner case in Fig. 12. The corresponding sensitivity for the molecule case is parameterised in terms of the misidentification probability P mis = N mis−id /N MC . The P mis as a function of the input Flatté parameter E f,0 is shown in Fig. 13.
The available computing resources result in limited samples of DPM [41] background. This, in combination with an efficient background suppression of the order B,gen ≈ 1 · 10 −10 , results in a very small number of surviving background events which introduces an uncertainty. The impact is estimated by scaling up the number of background events determined from the 90% confidence level upper limit, according to [223]. The uncertainty due to non-resonant background from pp → J/ψρ 0 was determined in a similar way. The effect on the sensitivity is represented by bracket markers in Figs. 12 and 13.
A more compact representation of the results extracted from Figs. 12 and 13 is shown in Fig. 14 for the Breit-Wigner scenario (left panel) and the molecule scenario (right panel). In the BW case, the minimum Γ 0 is defined by the minimum width, for which a 3σ sensitivity is achieved in an absolute decay width measurement. This corresponds to a relative uncertainty ∆Γ meas /Γ meas of 33 %. In the left panel of Fig. 14, the 3σ sensitivity is shown as a function of the input σp p→X,max . Trendlines for inter-and extrapolation are added using an empirical analytical function.
In the molecule case, the capability of distinguishing a bound state from a virtual state is quantified in terms of the Flatté parameter difference ∆E f := |E f,0 − E f,th |, where E f,0 is the input Flatté parameter and E f,th is the threshold energy separating a bound from a virtual state. The difference can be extracted from Fig. 13 at P mis = 10 %, assuming E f,th = −8.5651 MeV according to  Figure 12: Sensitivity to the absolute Breit-Wigner width, parameterised in terms of the relative uncertainties ∆Γmeas/Γmeas, shown as a function of the input decay width Γ0 of a narrow resonance for six different input signal cross-sections σS. All results are extracted for the Phase One HESR running mode. The inner error bars represent the statistical uncertainties and the outer the systematic ones. The bracket markers indicate the corresponding numbers for the case of DPM [41] and non-resonant background upscaling according to [223], ignoring statistical and systematic errors.
Ref. [216,217]. The results are shown as a function of the input cross section σp p→X,max in the right panel of Fig. 14.
As expected, the larger the cross section, the better the performance in resolving small ∆E f . The achievable sensitivities have been calculated for one out of the six input cross sections, σ S = 50 nb, in line with the experimental upper limit on pp → χ c1 (3872) production provided by the LHCb experiment. For values of the natural decay width larger than Γ 0 = 110 keV a 3σ relative error ∆Γ meas /Γ meas better than 33 %, is achieved already in Phase One with 80 days of dedicated beam time for one resonance scan measurement. The nature of the state -bound or virtual -can be correctly determined with a probability of 90 % provided for ∆E f ≈ 700 keV. The presented work serves as an example, but the same approach will be applied to narrow resonances in general, achieving sub-MeV resolutions.

Impact and long-term perspectives
The planned Phase One line-shape measurement of the χ c1 (3872) and other states with J P C = 1 −− can reveal the intriguing nature of hadronic states. Once the nature of these states has been understood, it might lead to further insights on the relevant degrees-of-freedom that give rise to hadronisation at different energy scales.
In addition, PANDA has excellent discovery potential for hitherto unknown, exotic meson-like states thanks to the gluon-rich environment provided bypp annihilations as well as the access to allqq-like quantum numbers in formation. In particular, this opens up for extensive searches for spin partners of the XY Z states. Discoveries and measurements of the properties of spin partners provide valuable insights on the prominent components, since different assumptions lead to different effects of spin symmetry violation [211].
In later phases of PANDA, when the design luminosity is reached, studies of hadrons with open charm will commence. The structure and dynamics of these systems, composed of heavy and light constituent quarks, are complementary to that of hidden-charm meson-like states. The decay of the lowest lying states occurs primarily via weak processes, providing experimental access to the semileptonic form factors and the CKM parameters. Moreover, spectroscopy studies of the excited states can provide new  Figure 13: Sensitivity to theD * D molecule scenario, parameterised in terms of the mis-identification probability Pmis, shown as a function of the input Flatté parameter E f,0 of the χc1(3872) for six different input signal cross-section σS. All results are extracted for the Phase One HESR running mode. The inner error bars represent the statistical uncertainty and the outer the systematic ones. The bracket markers indicate the corresponding numbers for the case of DPM [41] and non-resonant background upscaling according to [223], ignoring statistical and systematic uncertainties.  insights in the non-perturbative QCD domain that are not accessible in the hidden-charm sector. This opens the possibility to search for exotic open-charm states. Hence, PANDA can build upon the BABAR and CLEO discoveries of the narrow exotic candidates D * s0 (2317) [159] and D s1 (2460) [160], respectively. PANDA has the potential to measure the width of the D * s0 (2317) with a resolution in the order of 0.1 MeV via an energy scan near the threshold of the associated D ± s D * s0 (2317) ∓ production [224] and to search for other higher order excitations of open-charm states. This is particularly important since the width is sensitive to a possible molecular component of the state [15,212,225,226].

Hadrons in nuclei
Hadron reactions with nuclear targets provide a great opportunity to study how nuclear forces emerge from QCD. In particular, these reactions offer an angle to the onset of colour transparency at intermediate energies, the shortdistance structure of the nuclear medium, and the effects of the nuclear potential on hadron properties. Two important aspects make antiproton probes unique in this regard: -The kinematic threshold for the production of heavy mesons (e.g. charmonia, D, D * ) and antibaryons is accessible at small beam momenta. -The existence of two-body annihilation channels at large momentum transfer.
Close to threshold, the produced particles are rather slow in the laboratory frame. Since the coherence lengths are small compared to the internucleon distance, these particles interact with the nuclear residue as ordinary hadrons.
The probability for such multiple interactions is quantified by the nuclear transparency T (A) and is given by the ratio of the cross section of an exclusive nuclear process with the corresponding elementary (nucleon) reaction. The antiproton beam gives access to hadron channels that are difficult to study with other probes at low momenta, for example J/Ψ . Slow particles are influenced by the nuclear mean field potentials. Antiprotons are particularly suitable for implanting low-momentum antibaryons or mesons into the nuclear environment, where resulting effects of the nuclear potential on their masses and decay widths can be studied. Nuclear potentials are crucial to gain valuable insights into neutron stars [227].
At higher beam momenta, the factorisation theorem mentioned in Section 4 becomes valid, splitting the reaction into a hard, pQCD calculable part and a soft part described by GPDs. This relies on the assumption that soft gluonic exchanges between the incoming and outgoing quark configurations are suppressed, which in turn is only possible if these configurations are colour neutral and have transverse sizes substantially smaller than the normal hadron size. While well-established at large momentum transfer, it is still an open question at which scale this phenomenon, known as colour transparency (CT), sets in.
Interactions at large momentum transfers also probe the short-distance (≤ 1.2 fm) structure of the nuclear medium itself. In this region, effects from non-perturbative QCD discussed in sections 4 to 6 come into play in the dynamics of the nuclear repulsive core, a rather unexplored territory [228], which is expected to have an effect on cold, dense nuclear matter such as neutron stars.

State of the art
Nuclei made of protons and neutrons have been studied for more than a century. Hypernuclei, where one of the nucleons is replaced by a hyperon, and hyperatoms, where a hyperon is attached to a nucleus in an atomic orbit, have been explored since more than six decades. As a result, valuable information about the nuclear potentials of Λ and Σ − hyperons has been obtained [229].
It was recently pointed out in Ref. [230] that inmedium interactions of antibaryons may cause compressional effects and may thus provide additional information on the nuclear equation of state [231]. The data for antibaryons in nuclei are however rather scarce. So far, only antiprotons have been subjected to experimental studies. The antiproton optical potential has been addressed in studies of elasticpA scattering at KEK [232] and LEAR [233,234]. The fits to the angular distributions of the scattered antiprotons, indicate that the potential has a shallow attractive real part Re(V opt ) = −(0 − 70) MeV and a deep imaginary part Im(V opt ) = −(70 − 150) MeV in the centre of a nucleus. This is in contrast to results from the analysis of X-ray transitions in antiprotonic atoms and of radiochemical data. Here, the real part turned out to be much deeper, Re(V opt ) = −110 MeV, whereas the imaginary part was found to be Im(V opt ) = −160 MeV [235]. However, the calculations of thepA elastic scattering as well as those of antiprotonic atoms, are sensitive to thē p potential at the nuclear periphery. The production ofp in pA and AA collisions, on the other hand, is sensitive to the antiproton potential deeply inside the nuclei and seems to favour Re(V opt ) = −(100 − 250) MeV at normal nuclear density as predicted by microscopic transport calculations [236][237][238]. Antiproton absorption cross sections on nuclei as well as the π + and proton momentum spectra produced inp annihilation nuclei at LEAR calculated within the Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) model [239] are consistent with Re(V opt ) −150 MeV, i.e. about a factor of four weaker than expected from naive G-parity relations.
In Ref. [240] it has been suggested that this discrepancy is a consequence of the missing energy dependence of the proton-nucleus optical potential in conventional relativistic mean-field models. The energy-and momentum dependence required for such an effect can be recovered by extending the relativistic hadrodynamic Lagrangian by non-linear derivative interactions [240][241][242], hence mimicking many-body forces [243]. Since hyperons and antihyperons play an important role in the interpretation of high-energy heavy-ion collisions and dense hadronic systems, it needs to be investigated how these concepts carry over to the strangeness sector. However, antihyperons annihilate quickly in nuclei and conventional spectroscopic studies are therefore challenging or even unfeasible. Instead, quantitative information about the potentials can be obtained from exclusive antihyperon-hyperon production inpA annihilations close to threshold. However, so far no such experimental data exist on nuclear potentials of antihyperons.

Potential for Phase One
In the absence of feasible spectroscopic methods, schematic calculations performed in Refs. [244][245][246] indicate that the transverse momentum asymmetry where p T (Y /Ȳ )) is the transverse momentum of the hyperon/antihyperon, is sensitive to the depth of the antihyperon potential. Other observables of interest are polarisation and coplanarity. As concluded in Section 5.1.2, a unique feature of antiproton interactions within the PANDA energy range is the large production cross sections of hyperon-antihyperon pairs. However, due to the strong absorption of antibaryons in nuclei, the exclusive production rate of antihyperon-hyperon pairs is expected to be smaller in antiproton-nucleus collisions compared to antiprotonproton interactions.
Realistic calculations for the Phase One feasibility have been performed using the GiBUU transport model [247]. Here we show recent results obtained with GiBUU, release 2017, which incorporates, inter alia, updates in the kaon dynamics and an improved parametrisations of the hyperon-nucleon (S = -1) collision channels at low hyperon momenta with respect to the previously used release 1.5 [246,248]. Non-linear derivative interactions were not included. A simple scaling factor ξ p = 0.22 was applied for the antiproton potential to ensure a Schrödinger equivalent antiproton potential of about 150 MeV at saturation density [239]. Since no experimental information exists so far for antihyperons in nuclei, G-parity symmetry was adopted as a starting point. The calculations were carried out for different values of the antihyperon scaling factor ξ Y . The calculations were performed for the following cases: -ΛΛ pair production inp+ 20 Ne interactions at p beam = 1.52 GeV/c. A beam momentum of 1.64 GeV/c is also used for the study of the pp → ΛΛ which will serve as a point of reference. At the lower beam momentum of 1.52 GeV/c, the production of Σ is strongly suppressed, hence reducing experimental ambiguities. The resulting distributions of transverse asymmetry α T as a function of the longitudinal asymmetry α L , defined in the same way but with T → L, are shown in Figs. 15 (ΛΛ) and 16 (ΛΣ − andΞ + Ξ − ). ForΛΛ, we observe a remarkable sensitivity of α T to the potential at negative values of α L (Fig. 15), and it is clear that secondary effects do not wipe out the dependence. The large α T sensitivity as well as the negative shift in α T are linked to the substantial Λ transverse momentum smearing due to secondary scattering.
In order to estimate the expected event rate we assume an interaction rate of 10 6 s −1 , 20% beam loss in the HESR due to the complex target and a reconstruction efficiency of 10%, which is slightly smaller than that of the elementarypp →ΛΛ presented in Tab. 1. With these assumptions we can obtain 2 (1)ΛΛ per second for p beam = 1.64(1.52) GeV/c. One day of data taking with 90% effective run time at 1.64 GeV/c will yield 15·10 4 events, corresponding to a sample size two times as large as the one presented in Fig. 15. One week of data taking would also enable measurements of polarisation and coplanarity.
For the results presented in the right panel of Fig. 16, about 12000 Ξ − Ξ + pairs were generated for each value of the scaling factor ξ Ξ + . With the Phase One luminosity and a Ξ − Ξ + reconstruction efficiency of 5% (slightly smaller than that of the elementarypp →Ξ + Ξ − presented in Tab. 2), this requires a running time of about two months. The studies proposed here will benefit from measurements of the reference reaction pp → Y Y . However, as discussed in Section 5.1, such measurements already constitute an important part of the hyperon production programme and can, thanks to the predicted large production rate, be completed in a very short time. The results from our calculations illustrate that even with rather conservative assumptions about luminosity, PANDA can provide unique and relevant information on the behaviour of antihyperons in nuclei already during Phase One.

Impact and long-term perspectives
Already in Section 5.1, it was concluded that PANDA will be a strangeness factory. In combination with nuclear targets, this opens up unique possibilities for pioneering studies of the nuclear antihyperon potentials already during Phase One. In the future, when the luminosity is increased, a unique programme for double-and possibly triple strange hyperatom-and hypernuclear studies will follow [248].

State of the art
Colour Transparency (CT) has mainly been studied in the high-energy regime, e.g. at Fermilab and HERA [249]. At intermediate energies, some evidence was found by the CLAS collaboration for an onset of CT in exclusive meson production with electron probes at momentum transfers of a few GeV [250,251].
Two-body hadron-nucleus reactions are also sensitive to short-range nucleon-nucleon correlations [252]. These have been studied experimentally for example in twonucleon knockout reactions with proton beams at BNL [253,254] and with electron beams at JLab [255,256]. It was found that inside ground-state nuclei, the short- range nucleon-nucleon interaction can give rise to correlated nucleon pairs with large relative momenta but small centre-of-mass momenta, called short-range correlated (SRC) pairs.

Potential for Phase One
Despite describing different physics phenomena, CT and SRC can be studied with similar probes and momentum regimes and with similar methods. Reactions with antiproton probes have the advantage that they give access to mesons that are unlikely to be produced with electron beams, for example kaons.
To establish the onset of CT in the intermediate energy regime indicated by CLAS, studies of e.g. exclusive meson production inpp andpA have been suggested [257,258]. At large momentum transfer, a qq pair is more probable to be in a small-size configuration than a qqq triplet due to combinatorics. Therefore, two-meson annihilation channels provide a very promising search-ground for such studies. It is noteworthy that the main feature of the nuclear target, i.e. the possibility of initial-and final state interactions with spectator nucleons, can be explored already for the deuteron. The wave function of the deuteron is relatively well-known which allows for more robust theoretical predictions. The simplest opportunity to study CT is the d(p, π − π 0 )p process at large momentum transfer in the elementarypn → π − π 0 reaction [259]. The "golden" channel for nuclear transparency is considered to be A(p, J/ψ)(A−1) * . During Phase One, it will be difficult to study charmonium production for heavy nuclei due to the limited luminosity, but studies of the integrated cross section with a deuteron target may be started. Calculations of exclusive charmonium production d(p, J/ψ)n [260] predict a quite large cross section of ∼ 5 nb at the quasifree peak (p lab = 4.07 GeV/c).
The same two-body antiproton reactions can be used to study the decay of a short-range correlation after removal of one nucleon, for examplep + A → h 1 + h 2 + N back + (A − 2) * , where N back refers to a backward-going nucleon [254]. In these reactions, it is possible to test the validity of factorisation of the cross section into the elementary cross section, the decay function, and the absorption factor using different final states. Such tests in combination with analogous studies at JLab would contribute to detailed understanding of the dynamics of interactions with short-range correlations and high density fluctuations in nuclei.
In SRC studies, antiprotons give access to correlated pp and pn pairs without the necessity of identifying and determining the momentum of an outgoing neutron. Instead, a struck neutron can be identified by reconstructing processes likepn → π − π 0 orpn → π + π − π − in the PANDA detector. The wave function of the SRC pair may include the contribution of non-hadronic degrees of freedom. The simplest case is again provided by the deuteron wave function which may include the ∆ − ∆ component predicted by meson-exchange model calculations [261] as well as in the quark model [262]. The presence of the ∆ ++ − ∆ − configuration in the deuteron may be tested in the exclusive reactionpd → π − π − ∆ ++ [263]. In the PANDA momentum range, the signal process due to the antiproton annihilation on the valence ∆ − dominates over two-step background processes. This is valid in a broad kinematic range of the produced ∆ ++ also for ∆ − ∆ probabilities in the deuteron as low as ∼ 0.3%.

Impact and long-term perspective
At large beam momenta, PANDA can contribute with studies of colour transparency and short-range correlated nucleon-nucleon pairs, and offers access to final states which are difficult or unfeasible to study with electron or proton beams.
The larger luminosities of the later stages of PANDA will allow for more extensive studies of charmonium production A(p, J/ψ)(A − 1) * reactions, both for deuterium target and beyond. Exclusive studies of differential cross sections and J/ψ and ψ (2S) transparency ratios shed further light on colour transparency, as discussed in detail in Refs. [258,264].
The J/ψN absorption cross section is of particular interest for studies of Quark-Gluon Plasma in heavy-ion collisions [265].

Summary and conclusions
The Standard Model of particle physics gives an accurate description of phenomena that occur at very high energies exploiting the basic interactions among quarks and gluons. However, describing why and how these quarks and gluons form hadrons remains an open question. The most prominent examples are the building blocks of visible matter, i.e. protons and the neutrons. Furthermore, it is a challenge to describe quantitatively how the effective forces between these composite objects emerge from first principles: how do protons and neutrons form atomic nuclei, and how do these form the macroscopic objects of our universe, for example neutron stars?
A central theme in strong interaction phenomena is the non-Abelian nature of QCD, i.e. the self-coupling of the force carriers. Self-coupling is present in all non-Abelian theories such as gravity, but hadrons are so far the only objects for which these effects can be studied in a controlled way in the laboratory.
The PANDA experiment will provide an advanced and multi-faceted facility for studies of different aspects of the strong interaction. PANDA will utilise a beam of antiprotons: a unique and highly versatile probe for hadron physics. The beam energy provided by the HESR storage ring is optimised to shed light on the very regime where quarks form hadrons. Combined with a near 4π multipurpose detector, PANDA will offer the broadest hadron physics programme of any existing or planned experiment in the world.
The PANDA physics programme will benefit from the recent theoretical developments (lattice QCD, effective field theory, Dyson-Schwinger approach, AdS/CFT, etc.), but also provide guidance from data to the construction of new theoretical and phenomenological tools, as well as refinements of the existing ones. The close collaboration between theory and experiment will hence be mutually beneficial and has potential to give new insights in the dynamics of non-linear interacting systems on a quantum scale.
In this paper, we have discussed the potential of PANDA during the first phase of data collection, Phase One, when the luminosity will be ≈20 times lower than the FAIR design value and the experimental setup will be slightly reduced. The four main physics domains of PANDA -nucleon structure, strangeness physics, charm and exotics, and hadrons in nuclei -have been discussed within the context of Phase One. Highlights have been outlined and the potential for PANDA to push the frontiers beyond state of the art was demonstrated for selected examples. PANDA will be a key experiment that can investigate the nucleon structure, perform line-shape measurements of non-vector charmonium-like states, study multistrange hyperons at a large scale and antihyperons in nuclei. It offers complementary approaches to topics like time-like form factors, light hadron spectroscopy and colour transparency. In later phases of the PANDA experiment, the full setup and the design luminosity enable an even wider programme that also includes open-charm production, triple-strange hyperon physics, hyperatom and hypernuclear physics and searches for physics beyond the Standard Model e.g. through hyperon decays.