Upgrading the Baryonic Matter at the Nuclotron Experiment at NICA for Studies of Dense Nuclear Matter

The Nuclotron at the Joint Institute for Nuclear Research in Dubna can deliver gold beams with kinetic energies between 2 and 4.5 A GeV. In heavy-ion collisions at these energies, it is expected that the nuclear fireball will be compressed by up to approximately four times the saturation density. This offers the opportunity to study the high-density equation-of-state (EOS) of nuclear matter in the laboratory, which is needed for our understanding of the structure of neutron stars and the dynamics of neutron star mergers. The Baryonic Matter at the Nuclotron (BM@N) experiment will be upgraded to perform multi-differential measurements of hadrons including (multi-) strange hyperons, which are promising probes of the high-density EOS, and of new phases of quantum chromodynamic (QCD) matter. The layout of the upgraded BM@N experiment and the results of feasibility studies are presented.


Introduction
The experimental exploration of the high-density nuclear matter equation of state is an important goal for contemporary and future research, both in the laboratory and for astronomical observations. Heavy-ion beams provided by the future accelerator centers, the Facility for Antiproton and Ion Research (FAIR) in Darmstadt, Germany [1], and the Nuclotron-based Ion Collider Facility (NICA) at the Joint Institute for Nuclear Research (JINR) in Russia [2], offer unique possibilities to produce baryonic matter at high densities. The fixed-target Compressed Baryonic Matter (CBM) experiment at the FAIR [3] and the Multi-Purpose Detector (MPD) at the NICA collider [4] are designed to measure diagnostic probes, which will provide new information on the quantum chromodynamic (QCD) phase diagram at large baryon-chemical potentials, including the high-density equation-of-state (EOS). As an intermediate step, the Baryonic Matter at the Nuclotron (BM@N) fixed-target experiment has been installed at the Nuclotron at the JINR [5]. The current version of the BM@N experiment is well suited to perform experiments with beams of light ions. To measure Au + Au collisions at beam energies up to 4.5 A GeV and with reaction rates up to 20 kHz, the experiment will be upgraded with an evacuated beam pipe, with a hybrid tracking system comprising of 4 stations of double-sided microstrip silicon sensors and 7 planes of two-coordinate gaseous electron multiplier (GEM) detectors (already partly existing), and with a zero degree calorimeter for event characterization. The silicon tracking detector system is presently under development in close cooperation with the silicon detector group of the CBM collaboration. This article will briefly outline the envisaged physics program of the upgraded BM@N experiment and the first results of the physics performance simulations.
The Nuclotron beam energies are well suited to produce dense nuclear matter in the collision of heavy nuclei. This is illustrated in Figure 1, which depicts the central density in the reaction volume of a central Au + Au collision at a beam energy of 5 A GeV as a function of time, as predicted by several transport models and a hydro-dynamical calculation [6]. According to these calculations, the nuclear fireball will be compressed to about 5 times the saturation density ρ o at a beam kinetic energy of 5 A GeV. Therefore, experiments using the Nuclotron beam allow to study the EOS at neutron star core densities. Depending on the EOS, the central densities of neutron stars with about two solar masses vary between 4 and 8 ρ o [7]. The central densities increase for neutron star models with smaller maximum mass, i.e., for softer equations of state.
Particles 2019, 2 FOR PEER REVIEW 2 comprising of 4 stations of double-sided microstrip silicon sensors and 7 planes of two-coordinate gaseous electron multiplier (GEM) detectors (already partly existing), and with a zero degree calorimeter for event characterization. The silicon tracking detector system is presently under development in close cooperation with the silicon detector group of the CBM collaboration. This article will briefly outline the envisaged physics program of the upgraded BM@N experiment and the first results of the physics performance simulations. The Nuclotron beam energies are well suited to produce dense nuclear matter in the collision of heavy nuclei. This is illustrated in Figure 1, which depicts the central density in the reaction volume of a central Au + Au collision at a beam energy of 5 A GeV as a function of time, as predicted by several transport models and a hydro-dynamical calculation [6]. According to these calculations, the nuclear fireball will be compressed to about 5 times the saturation density ρ0 at a beam kinetic energy of 5 A GeV. Therefore, experiments using the Nuclotron beam allow to study the EOS at neutron star core densities. Depending on the EOS, the central densities of neutron stars with about two solar masses vary between 4 and 8 ρo [7]. The central densities increase for neutron star models with smaller maximum mass, i.e., for softer equations of state.

The High-Density Nuclear-Matter Equation-of-State
A fundamental ingredient for our understanding of nuclear matter is its equation-of-state (EOS), which determines the evolution of supernova explosions, the size and radius of neutron stars, and the dynamics of neutron star mergers. The EOS describes the relationship between density, pressure, volume, temperature, energy, and isospin asymmetry, and can be expressed as the energy per nucleon as a function of density: with the isospin asymmetry parameter δ = (ρn -ρp)/ρ. Symmetric matter is stable around saturation density ρ0 with a binding energy of E/A(ρo) = −16 MeV at saturation density ρo, the slope δ(E/A)(ρo)/δρ = 0, and the curvature Knm = 9ρ 2 δ 2 (E/A)/δρ 2 with Knm the nuclear incompressibility. The EOS for symmetric nuclear matter is shown in Figure 2 for soft and hard Skyrme forces, and from microscopic ab-initio calculations [8].
Experiments at GSI have investigated the EOS of symmetric matter up to about 2 ρo by measuring the excitation of K + meson production in Au + Au and C + C collisions at sub-threshold beam energies from 0.8 and 1.5 A GeV [9]. Sub-threshold kaon production proceeds via multiple collisions involving pions and Delta resonances. At high densities, which are reached if the EOS is soft, multiple collisions occur more frequently, and the kaon yield is enhanced. The measured data can be reproduced by transport model calculations when taking into account momentum-dependent interactions and a soft EOS with a nuclear incompressibility of about 200 MeV [8,10]. The data on the elliptic flow of protons, deuterons, tritons, and 3 He in Au + Au collisions at beam kinetic energies

The High-Density Nuclear-Matter Equation-of-State
A fundamental ingredient for our understanding of nuclear matter is its equation-of-state (EOS), which determines the evolution of supernova explosions, the size and radius of neutron stars, and the dynamics of neutron star mergers. The EOS describes the relationship between density, pressure, volume, temperature, energy, and isospin asymmetry, and can be expressed as the energy per nucleon as a function of density: with the isospin asymmetry parameter δ = (ρ n -ρ p )/ρ. Symmetric matter is stable around saturation density ρ o with a binding energy of E/A(ρ o ) = −16 MeV at saturation density ρ o , the slope δ(E/A)(ρ o )/δρ = 0, and the curvature K nm = 9ρ 2 δ 2 (E/A)/δρ 2 with K nm the nuclear incompressibility. The EOS for symmetric nuclear matter is shown in Figure 2 for soft and hard Skyrme forces, and from microscopic ab-initio calculations [8].
between 0.4 and 1.5 A GeV have been explained by the isospin-dependent quantum molecular dynamics (IQMD) transport calculations when taking into account momentum-dependent interactions and assuming a soft EOS for symmetric matter, with nuclear incompressibility of Knm = 190 ± 30 MeV [11]. The EOS with hard Skyrme forces, i.e., with nuclear incompressibility of Knm = 380 MeV, as shown in Figure 2 is in contradiction with these results. The BM@N experiment will explore the EOS of symmetric matter up to about 4 times the saturation density by measuring the elliptic flow of protons, which is driven by the pressure gradient in the fireball and, therefore, it is sensitive to the EOS. Pioneering experiments at the Alternating Gradient Synchrotron (AGS) in Brookhaven have already measured both the directed and the elliptic flow of protons in Au + Au collisions at beam kinetic energies between 2 and 11 A GeV [12]. However, the interpretation of these flow data using a transport model was not conclusive, in particular in the Nuclotron energy range. The data on the directed flow were compatible with a soft EOS, while elliptic flow data indicate a hard EOS [13]. Therefore, the aim of the BM@N experiment is to repeat these measurements with improved statistics and precision.
Another very promising observable that is sensitive to the high-density EOS of symmetric matter is the excitation function of multi-strange (anti-)hyperons. According to transport models, Ξ and Ω hyperons are produced in sequential collisions involving kaons and lambdas and, therefore, are sensitive to the density of the fireball [14,15]. This sensitivity is expected to increase at beam energies close to or even below the production threshold. For example, the production of a Xi hyperon in a proton-proton collision pp → Ξ − K + K + p requires a minimum proton energy of Ethr = 3.7 GeV. Nevertheless, Xi hyperons can be produced in heavy-collisions via strangeness exchange reaction at energies above the lambda threshold of 1.6 GeV via lambda-lambda collisions Λ 0 Λ 0 → Ξ − p. Even Omega hyperons can be produced, when the Xi collides with another lambda Λ 0 Ξ − → Ω − n or with a kaon K − Ξ − → Ω − π − . Therefore, the measurement of multi-strange hyperons will be a central goal of the BM@N research program at the Nuclotron.

Exploring Mixed Phases of QCD Matter
According to chiral mean field model calculations, a transition from nuclear to quark matter is expected during the evolution of a neutron star merger with a total mass of 2.8 solar masses, when a density between 3 and 4 ρ0 and at a temperature of about 50 MeV is reached [16]. Such densities, although at higher temperatures, can be created in Au + Au collisions at Nuclotron energies. A promising observable, which is expected to be sensitive to the onset of such a phase transition, is the yield of multi-strange hyperons. In heavy-ion collisions at ultra-relativistic energies, the yield of Experiments at GSI have investigated the EOS of symmetric matter up to about 2 ρ o by measuring the excitation of K + meson production in Au + Au and C + C collisions at sub-threshold beam energies from 0.8 and 1.5 A GeV [9]. Sub-threshold kaon production proceeds via multiple collisions involving pions and Delta resonances. At high densities, which are reached if the EOS is soft, multiple collisions occur more frequently, and the kaon yield is enhanced. The measured data can be reproduced by transport model calculations when taking into account momentum-dependent interactions and a soft EOS with a nuclear incompressibility of about 200 MeV [8,10]. The data on the elliptic flow of protons, deuterons, tritons, and 3 He in Au + Au collisions at beam kinetic energies between 0.4 and 1.5 A GeV have been explained by the isospin-dependent quantum molecular dynamics (IQMD) transport calculations when taking into account momentum-dependent interactions and assuming a soft EOS for symmetric matter, with nuclear incompressibility of K nm = 190 ± 30 MeV [11]. The EOS with hard Skyrme forces, i.e., with nuclear incompressibility of K nm = 380 MeV, as shown in Figure 2 is in contradiction with these results.
The BM@N experiment will explore the EOS of symmetric matter up to about 4 times the saturation density by measuring the elliptic flow of protons, which is driven by the pressure gradient in the fireball and, therefore, it is sensitive to the EOS. Pioneering experiments at the Alternating Gradient Synchrotron (AGS) in Brookhaven have already measured both the directed and the elliptic flow of protons in Au + Au collisions at beam kinetic energies between 2 and 11 A GeV [12]. However, the interpretation of these flow data using a transport model was not conclusive, in particular in the Nuclotron energy range. The data on the directed flow were compatible with a soft EOS, while elliptic flow data indicate a hard EOS [13]. Therefore, the aim of the BM@N experiment is to repeat these measurements with improved statistics and precision.
Another very promising observable that is sensitive to the high-density EOS of symmetric matter is the excitation function of multi-strange (anti-)hyperons. According to transport models, Ξ and Ω hyperons are produced in sequential collisions involving kaons and lambdas and, therefore, are sensitive to the density of the fireball [14,15]. This sensitivity is expected to increase at beam energies close to or even below the production threshold. For example, the production of a Xi hyperon in a proton-proton collision pp → Ξ − K + K + p requires a minimum proton energy of E thr = 3.7 GeV. Nevertheless, Xi hyperons can be produced in heavy-collisions via strangeness exchange reaction at energies above the lambda threshold of 1.6 GeV via lambda-lambda collisions Λ 0 Λ 0 → Ξ − p. Even Omega hyperons can be produced, when the Xi collides with another lambda Λ 0 Ξ − → Ω − n or with a kaon K − Ξ − → Ω − π − . Therefore, the measurement of multi-strange hyperons will be a central goal of the BM@N research program at the Nuclotron.

Exploring Mixed Phases of QCD Matter
According to chiral mean field model calculations, a transition from nuclear to quark matter is expected during the evolution of a neutron star merger with a total mass of 2.8 solar masses, when a density between 3 and 4 ρ o and at a temperature of about 50 MeV is reached [16]. Such densities, although at higher temperatures, can be created in Au + Au collisions at Nuclotron energies. A promising observable, which is expected to be sensitive to the onset of such a phase transition, is the yield of multi-strange hyperons. In heavy-ion collisions at ultra-relativistic energies, the yield of (anti-) hadrons including multi-strange (anti-) hyperons and light nuclei can be perfectly described by the statistical hadronization model, which assumes chemical equilibrium [17]. In view of the fact, that the hyperon-nucleon scattering cross-section is small, this observation was interpreted as an indication that the system had undergone a transition from a partonic phase to the hadronic final state, with the equilibration of multi-strange hyperons being driven by multi-body collisions in the high-particle density regime near the phase boundary [18]. Agreement of the hyperon yield with thermal model calculations was also found at a beam kinetic energy of 40 A GeV in Pb + Pb collisions at the SPS [19]. The particle yields measured in Ar + KCl collisions at a beam kinetic energy of 1.76 A GeV can also be explained by the thermal model, except for Ξ − hyperons. The yield of Ξ − hyperons exceeds the thermal model prediction by about a factor 24 ± 9 [20], indicating that these particles are far off the chemical equilibrium. Therefore, the measurement of multi-strange hyperons in Au + Au collisions at different Nuclotron beam energies will explore the onset of equilibration of multi-strange hyperons at high net-baryon densities.

The BM@N Experiment
The layout of the proposed BM@N configuration is shown in Figure 3. The tracks of charged particles will be determined with two detector systems: 4 stations of double-sided microstrip silicon sensors downstream the target, and a set of 7 planes of two-coordinate GEM (gaseous electron multiplier) detectors mounted downstream of the silicon stations. Both the silicon tracking system (STS) and the GEM stations will be operated in the magnetic field (maximum value 1.2 T) of a large aperture dipole magnet with a gap of 1 m. The silicon detector stations are currently under development for the CBM experiment [21]. The tracking system behind the magnet actually consists of 2 large drift chambers (DCH), which will be replaced by cathode strip chambers for track measurement in Au + Au collisions. For particle identification, the time-of-flight will be measured with 3 walls of resistive plate chambers (mRPC) with strip readout. The T0 detector is positioned around the target and will provide the start signal for the time-of-flight (TOF) measurements together with the trigger signal for the data acquisition. The total energy of the projectile spectator fragments will be measured by a zero degree calorimeter (ZDC) to determine the centrality of the collision and the orientation of the reaction plane. The existing BM@N detector systems are described in [22].

Radiation Calculations
The current version of the BM@N experiment is used to measure reactions with low-intensity beams of light ions, which traverse the setup in the air without any beam pipe. For future measurements with a high-intensity gold beam, the setup has to be upgraded with an evacuated beam pipe. Such a beam pipe has been designed, consisting of carbon fiber tubes of 1 mm wall thickness with increasing diameter and several kinks to follow the deflection of the beam in the magnetic field. This design has been validated by FLUKA [23] calculations assuming a gold beam with a kinetic energy of 4.5 A GeV, a profile with a width of σ = 3.5 mm, a divergence of 1 mrad, and an intensity of 2 × 10 6 ions/s, hitting a gold target with a thickness of 250 µm located at Z = 0. The resulting distribution of the gold ions projected in the X-Z plane is shown in Figure 4. The calculations demonstrate, that the curvature and the increasing diameter of the proposed beam pipe are well suited to guide the high-intensity beam and to minimize radiation outside the pipe. The FLUKA calculation takes into account the divergence of the beam, the multiple scattering of the gold ions in the target, and the reactions of the beam in the material of the beam pipe. The particles produced in the target including the delta-electrons are not shown in Figure 4.
(STS) and the GEM stations will be operated in the magnetic field (maximum value 1.2 T) of a large aperture dipole magnet with a gap of 1 m. The silicon detector stations are currently under development for the CBM experiment [21]. The tracking system behind the magnet actually consists of 2 large drift chambers (DCH), which will be replaced by cathode strip chambers for track measurement in Au + Au collisions. For particle identification, the time-of-flight will be measured with 3 walls of resistive plate chambers (mRPC) with strip readout. The T0 detector is positioned around the target and will provide the start signal for the time-of-flight (TOF) measurements together with the trigger signal for the data acquisition. The total energy of the projectile spectator fragments will be measured by a zero degree calorimeter (ZDC) to determine the centrality of the collision and the orientation of the reaction plane. The existing BM@N detector systems are described in [22].  Particles 2019, 2 FOR PEER REVIEW 5 The current version of the BM@N experiment is used to measure reactions with low-intensity beams of light ions, which traverse the setup in the air without any beam pipe. For future measurements with a high-intensity gold beam, the setup has to be upgraded with an evacuated beam pipe. Such a beam pipe has been designed, consisting of carbon fiber tubes of 1 mm wall thickness with increasing diameter and several kinks to follow the deflection of the beam in the magnetic field. This design has been validated by FLUKA [23] calculations assuming a gold beam with a kinetic energy of 4.5 A GeV, a profile with a width of σ = 3.5 mm, a divergence of 1 mrad, and an intensity of 2 × 10 6 ions/s, hitting a gold target with a thickness of 250 µm located at Z = 0. The resulting distribution of the gold ions projected in the X-Z plane is shown in Figure 4. The calculations demonstrate, that the curvature and the increasing diameter of the proposed beam pipe are well suited to guide the high-intensity beam and to minimize radiation outside the pipe. The FLUKA calculation takes into account the divergence of the beam, the multiple scattering of the gold ions in the target, and the reactions of the beam in the material of the beam pipe. The particles produced in the target including the delta-electrons are not shown in Figure 4.

Radiation Calculations
In the region of the silicon stations, the ionizing dose reaches values of about 10 Gy after 2 months of beam on target. Assuming a total running time of the BM@N experiment of 5 years and 4 months per year, a dose of about 100 Gy is reached, resulting in mild damage to the sensors in the inner part of the stations. The equivalent neutron fluence in the silicon stations reaches values below 10 10 neq/cm 2 after 2 months, corresponding to a lifetime fluence of 10 11 neq/cm 2 , which is well within the radiation tolerance of the sensors.
In the region of the GEM stations, the ionizing dose reaches values of about 1 Gy after 2 months of beam on target, corresponding to a lifetime dose of 10 Gy. The equivalent neutron fluence in the GEM stations also reaches values below 10 10 neq/cm 2 after 2 months, corresponding to a lifetime fluence of 10 11 neq/cm 2 . Both values can be tolerated by the GEM detectors.

Physics Performance Simulations of the BM@N Hybrid Tracking System
As a first step, the lambda hyperon reconstruction performance of the BM@N hybrid tracking system has been studied based on 1000 central Au + Au collisions with a beam kinetic energy of 4 A GeV produced with the DCM-QGSM event generator [24]. For track reconstruction, the latest version of the CBM Cellular Automaton (CA) tracking code was used [25]. The BM@N hybrid tracking system comprises of 4 silicon stations and 6 GEM stations. The silicon stations are located at 30, 50, 70, and 90 cm downstream of the target, followed by GEM detectors at 120 cm up to 270 cm, with a gap of 30 cm between the stations.
The reconstruction efficiency for primary tracks as a function of momentum is presented in the left panel of Figure 5. The primary tracks, which create only hits in the 4 silicon stations, can be reconstructed with an efficiency of better than 90% above a momentum of about 0.6 GeV/c, whereas the track reconstruction efficiency for the hybrid system reaches a maximum of about 90% between 1 and 2.5 GeV/c, and then drops to about 80% at 6 GeV/c. The reason for the lower efficiency of the Silicon + GEM system is the low granularity of the GEMs, which leads to a large number of clone hits being misinterpreted as real hits. The momentum resolution for primary particles as a function of momentum is shown in the right panel of Figure 5 for the hybrid tracker. Only for momenta below p = 0.5 GeV/c, the momentum resolution is worse than Δp/p = 0.006. In the region of the silicon stations, the ionizing dose reaches values of about 10 Gy after 2 months of beam on target. Assuming a total running time of the BM@N experiment of 5 years and 4 months per year, a dose of about 100 Gy is reached, resulting in mild damage to the sensors in the inner part of the stations. The equivalent neutron fluence in the silicon stations reaches values below 10 10 n eq /cm 2 after 2 months, corresponding to a lifetime fluence of 10 11 n eq /cm 2 , which is well within the radiation tolerance of the sensors.
In the region of the GEM stations, the ionizing dose reaches values of about 1 Gy after 2 months of beam on target, corresponding to a lifetime dose of 10 Gy. The equivalent neutron fluence in the GEM stations also reaches values below 10 10 n eq /cm 2 after 2 months, corresponding to a lifetime fluence of 10 11 n eq /cm 2 . Both values can be tolerated by the GEM detectors.

Physics Performance Simulations of the BM@N Hybrid Tracking System
As a first step, the lambda hyperon reconstruction performance of the BM@N hybrid tracking system has been studied based on 1000 central Au + Au collisions with a beam kinetic energy of 4 A GeV produced with the DCM-QGSM event generator [24]. For track reconstruction, the latest version Particles 2019, 2 486 of the CBM Cellular Automaton (CA) tracking code was used [25]. The BM@N hybrid tracking system comprises of 4 silicon stations and 6 GEM stations. The silicon stations are located at 30, 50, 70, and 90 cm downstream of the target, followed by GEM detectors at 120 cm up to 270 cm, with a gap of 30 cm between the stations.
The reconstruction efficiency for primary tracks as a function of momentum is presented in the left panel of Figure 5. The primary tracks, which create only hits in the 4 silicon stations, can be reconstructed with an efficiency of better than 90% above a momentum of about 0.6 GeV/c, whereas the track reconstruction efficiency for the hybrid system reaches a maximum of about 90% between 1 and 2.5 GeV/c, and then drops to about 80% at 6 GeV/c. The reason for the lower efficiency of the Silicon + GEM system is the low granularity of the GEMs, which leads to a large number of clone hits being misinterpreted as real hits. The momentum resolution for primary particles as a function of momentum is shown in the right panel of Figure 5 for the hybrid tracker. Only for momenta below p = 0.5 GeV/c, the momentum resolution is worse than ∆p/p = 0.006. Lambda hyperons are reconstructed by the cellular automaton (CA) track finder by identifying the vertex of the lambda decay into two charged particles (Λ → p + π − ), assuming that the decay products are either a pion or a proton, and vice versa. It is worthwhile to note, that no particle identification of the lambda decay products has been applied. The resulting proton-pion invariant mass spectrum is shown in Figure 6. The mass resolution of the lambdas is σ = 1.1 MeV/c 2 , the signalto-background ratio is S/B = 5.2.  Lambda hyperons are reconstructed by the cellular automaton (CA) track finder by identifying the vertex of the lambda decay into two charged particles (Λ → p + π − ), assuming that the decay products are either a pion or a proton, and vice versa. It is worthwhile to note, that no particle identification of the lambda decay products has been applied. The resulting proton-pion invariant mass spectrum is shown in Figure 6. The mass resolution of the lambdas is σ = 1.1 MeV/c 2 , the signal-to-background ratio is S/B = 5.2.
The phase space distributions of lambdas and their decay products emitted in central Au + Au collisions at a beam kinetic energy of 4 A GeV as generated with the DCM-QGSM code are shown in the upper row of Figure 7 (lambdas left, protons center, pions right). The proton and pion identification has been performed based on Monte Carlo information. The lower row depicts the distribution of lambdas (left), decay protons (center), and decay pions (right) reconstructed by the hybrid tracking system. According to Figure 7, about 3800 lambdas are produced in 1000 central Au + Au collisions at 4 A GeV kinetic beam energy. The number of the decay products correspond to the branching ratio of the lambda decay Λ → p + π − , which is BR = 0.64. The lambda reconstruction efficiency is slightly above 10% for the hybrid detector system. Lambda hyperons are reconstructed by the cellular automaton (CA) track finder by identifying the vertex of the lambda decay into two charged particles (Λ → p + π − ), assuming that the decay products are either a pion or a proton, and vice versa. It is worthwhile to note, that no particle identification of the lambda decay products has been applied. The resulting proton-pion invariant mass spectrum is shown in Figure 6. The mass resolution of the lambdas is σ = 1.1 MeV/c 2 , the signalto-background ratio is S/B = 5.2.  Figure 7, about 3800 lambdas are produced in 1000 central Au + Au collisions at 4 A GeV kinetic beam energy. The number of the decay products correspond to the branching ratio of the lambda decay Λ → p + π − , which is BR = 0.64. The lambda reconstruction efficiency is slightly above 10% for the hybrid detector system.

Expected Particle Yields
The ultimate goal of the BM@N research program devoted to the nuclear EOS would be to measure the excitation functions for the production of Λ, Ξ − , and Ω − hyperons in central Au + Au and central C + C collisions for beam kinetic energies of 2, 3, and 4 A GeV. The multiplicities for Λ and Ξ − hyperons per central Au + Au collision as predicted by a statistical hadronization model [26] are listed in Table 1 for the 3 bombarding energies, together with the expected yields recorded per week. The collision rate is assumed to be 10 kHz. The yield estimate takes into account the reconstruction efficiencies (10% for Λ, and 1% for Ξ − ) and the respective branching ratios. According to these numbers, no Ξ − hyperon can be reconstructed in 1000 central Au + Au collisions, which have been simulated so far. In order to demonstrate the feasibility of Ξ − or even Ω − hyperon measurements, simulations with much better statistics have to be performed.

Expected Particle Yields
The ultimate goal of the BM@N research program devoted to the nuclear EOS would be to measure the excitation functions for the production of Λ, Ξ − , and Ω − hyperons in central Au + Au and central C + C collisions for beam kinetic energies of 2, 3, and 4 A GeV. The multiplicities for Λ and Ξ − hyperons per central Au + Au collision as predicted by a statistical hadronization model [26] are listed in Table 1 for the 3 bombarding energies, together with the expected yields recorded per week. The collision rate is assumed to be 10 kHz. The yield estimate takes into account the reconstruction efficiencies (10% for Λ, and 1% for Ξ − ) and the respective branching ratios. According to these numbers, no Ξ − hyperon can be reconstructed in 1000 central Au + Au collisions, which have been simulated so far. In order to demonstrate the feasibility of Ξ − or even Ω − hyperon measurements, simulations with much better statistics have to be performed.  [26], and the yields measured per week taking into account the reconstruction efficiencies (10% for Λ, and 1% for Ξ − ) and the branching ratios (BR = 0.64 for Λ → p + π − and also for Ξ − → Λ + π − → p + π − + π − ). The kinetic beam energies are indicated, the collision rate is 10 kHz. According to statistical model predictions, also the measurement of hypernuclei should be possible at Nuclotron beam energies [27]. The multiplicity of 3 Λ H hypernuclei per central Au + Au collisions at 4 A GeV is predicted to be 10 -2 , and the multiplicity of both 4 Λ H and 4 Λ He hypernuclei is 10 −3 . The decay topology of 3 Λ H → 3 He + π − and 4 Λ H → 4 He + π − is similar to the lambda decay. Assuming the hypernuclei reconstruction has the same branching ratio as for the lambdas and a reconstruction efficiency of 1%, a measured yield of about 3 × 10 4 3 Λ H and 3 × 10 3 4 Λ H per week of running at a collision rate of 1 kHz should be feasible. The reconstruction of the decay 4 Λ He → 3 He + p + π − should be also possible. These measurements will allow one to determine the lifetime of these hypernuclei with high precision. It should be noted, that the yields mentioned above have been estimated without taking into account a duty factor, neither for the accelerator nor for the experiment

Conclusions and Outlook
The upgrade of the BM@N experiment with an evacuated beam pipe, a hybrid tracking system comprising of 4 stations of double-sided microstrip silicon sensors and 7 planes of two-coordinate GEM detectors (already partly existing), and a zero degree calorimeter for event characterization, will allow one to investigate Au + Au collisions at Nuclotron energies with interaction rates up to several kHz. These measurements open the possibility, to explore high-density EOS for symmetric matter, and to scout the location of a possible mixed phase of QCD matter. Both the detector developments and the feasibility studies are works in progress. The next simulation steps include particle identification via time-of-flight (TOF), and a substantial increase of the event statistics to be able to demonstrate the feasibility to reconstruct multi-strange hyperons and light hypernuclei. It is worthwhile to note that up to date no multi-strange hyperons have been measured in heavy collision systems in the Nuclotron beam energy range. The goal of the BM@N experiment is to provide the first data on these promising probes of dense QCD matter.