Triple high energy nuclear and hadron collisions - a new method to study QCD phase diagram at high baryonic densities

We propose an entirely new method to study the phase diagram of strongly interacting matter by means of scattering the two colliding beams at the fixed target. Here we present the results of simulations of the most central triple nuclear collisions with the UrQMD 3.4 model for the beam center-of-mass collision energies sNN=2.76TeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_{NN}} = 2.76\ \hbox {TeV}$$\end{document} and sNN=200GeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_{NN}} = 200\ \hbox {GeV}$$\end{document}. The main outcome of our modeling is that even at these very high collision energies the initial baryonic charge densities are about 3 times higher than the ones achieved in the ordinary binary nuclear collisions. As a result, for instance, the yields of protons and Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda $$\end{document}-hyperons are strongly enhanced in the triple nuclear collisions. The other prospective applications of this method are briefly discussed. Among them we consider the low energy collisions of three nuclei of lead, passing through an intermediate system with an electric charge of 246 units which exceeds essentially the critical value of 173 and, hence, this may be of crucial importance to study the spontaneous emission of positron-electron pairs from the vacuum. We present the convincing arguments that the triple nuclear collisions method will allow the high energy nuclear physics community to create a new frontier in the studies of the QCD phase diagram and to lift up these studies to an entirely new level.


Introduction
Investigation of the phase diagram of quantum chromodynamics (QCD) is the prime task of several ambitious experimental programs that are supplemented by rather sophisticated theoretical and numerical approaches. In the ongoing experiments on high energy nucleus-nucleus (A+A) collisions the center-of-mass energy of collisions ranges from √ s N N = 2.42 GeV (HADES), √ s N N = 3.1 GeV at RHIC BNL (Beam Energy Scan II) to √ s N N = 17.3 GeV at SPS CERN (NA61/SHINE) and further up to √ s N N = 5.02 TeV at the ALICE CERN. The major goals of these and completed experimental programs in the discussed collision energy range were to find the new state of matter, the quark-gluon plasma, to determine the collision energy thresholds of two phase transitions (PTs) expected to exist in QCD matter and to locate their (tri)critical endpoint(s). Unfortunately, after almost forty years of experimental studies of the QCD phase diagram nowadays we have only serious arguments that the quark-gluon plasma has, indeed, been created in A+A collisions, while two other goals are still far from being achieved [1] 1 .
Such a situation persists to exist due to objective reasons, the main of which are (I) the absence of PTs in a strict thermodynamic sense in A+A collisions due to the finite (or even small!) size of systems created in these reactions, (II) the inability of the lattice formulation of QCD to guide the experiments on A+A collisions and QCD phenomenology at high baryonic charge densities and (III) the lack of detailed experimental data measured with small steps in collision energy. Nevertheless, in our opinion under such hard circumstances, namely without a rigorous theory of phase transitions in finite systems, without clear and unambiguous signals of QCD PTs, the experimentalists and QCD phenomenologists did a great job and got some outstanding results and discoveries being more often guided by their enthusiasm and good luck rather than by the firm theoretical foundations. But the question is whether with the existing experimental methods and theoretical know-how the high energy nuclear physics (HENP) community will be able to accomplish its mission and to achieve the major goals during, let's say, a decade from now?
During the last 20 years, many promising theoretical results were obtained. But in our opinion on the theory side, the Functional Renormalization Group Method [2][3][4] and Dyson-Schwinger equations [5][6][7] provide the most coherent and consistent framework to study the QCD phase diagram. On the side of heavy ion collision phenomenology, two groups, using entirely different approaches to analyze the data, almost simultaneously came to similar conclusions: that the first order PT of (partial) chiral symmetry restoration to a phase of nearly massless hadrons occurs at the center-ofmass energies √ s N N 4−5 GeV [8][9][10][11][12], and that the deconfinement PT of color degrees of freedoms is either a very weak first-order PT or a second-order one with the threshold energy of A+A collisions being about √ s N N 9 − 10 GeV [8][9][10][11][12]. A conclusion that at finite values of baryonic chemical potential the chiral symmetry restoration PT occurs prior to the deconfinement one, i.e. at lower baryonic charge densities, obtained from the analysis of experimental data is in line with the well-known theoretical approaches [13][14][15]. However, the situation with the location of the (tri)critical endpoint is somewhat confusing yet [16,17], although the new estimates on its location made in Ref. [4] show that the situation is not hopeless. In fact, at present it is not even known for sure whether in QCD there exists a tricritical endpoint of these two PTs, or each PT has its own critical endpoint. Therefore, in order to accomplish its main mission, the HENP community has to do what it has not been done before. It is clear that on the theory side it is absolutely necessary to develop a rigorous theory of liquid-gas PTs in finite systems. Despite a limited progress achieved on a basis of the exact analytical solutions found for finite systems for a wide class of cluster models [18], the formulation of a rigorous theory of the liquid-gas PTs in finite systems is still in its infancy. The development of such an approach will move our experimental searches for the (tri)critical endpoint in A+A collisions to a solid ground.
On the experimental side, it is also necessary to get a reliable source of information about the equation of state of QCD matter which would be alternative and complementary to the A+A collision experiments. The importance of such an approach was understood recently and, hence, several groups independently suggested to employ the neutron stars and their coalescence to study the QCD matter equation of state at very high baryonic densities. Thus, neutron stars provide a possibility to test the region of the phase diagram from the vanishing temperature, typical for an equilibrated neutron star, up to around 70 MeV, reached during the neutron star mergers [19,20]. Moreover, the astrophysical observations, e.g. observations of two neutron stars with mass above 2M [21,22], simultaneous measurements of mass and radius of compact stars with NICER telescope [23,24], together with the tidal deformability limit obtained from the first binary neutron star merger GW170817 [25] have already constrained the equation of state at high density. Despite a big progress compact stars and on-ground A+A experiments so far can not probe the same part of the QCD phase diagram. As a step towards this direction the future experimental programs, i.e. NICA (JINR) [27] and FAIR (GSI) [26], are aimed to study the phase diagram of strongly interacting matter at highest baryon densities ever created in terrestrial laboratories.
In order to get an additional and independent information about the QCD matter equation of state at higher baryonic densities than the ones achieved in A+A collisions, here we suggest to perform the triple nuclear collisions (TNC) by bombarding the fixed target with two colliding beams. In this work, we mainly demonstrate the advantages of the TNC at very high collision energies of two beams √ s N N = 200 GeV (RHIC BNL) and √ s N N = 2760 GeV (LHC CERN), while some results for the low energies of collisions are also discussed. The work is organized as follows. The next section is devoted to the analysis of TNC at LHC and the highest RHIC energies of collisions. In Sect. 3 we discuss new abilities and new challenges of the TNC method, while our conclusions are summarized in Sect. 4.

Modeling triple nuclear collisions by UrQMD 3.4
The idea to search for the TNC happening in the fixed solid target bombarded by two colliding beams was presented in 2018 at the Conference "CERN-Ukraine cooperation: current state and prospects" [28]. The TNC require two colliding beams and, hence, we are forced to study the high energy domain first, since the collider experiments at LHC CERN and RHIC BNL were designed to rather high collision ener- gies. Moreover, the possibility to use the LHCb-like detector for the fixed target experiments at LHC energies [29][30][31] and the STAR Fixed Target Program has been discussed already [32]. Now we suggest to extend this frame to employ it for investigating the TNC.
In order to demonstrate new opportunities of the TNC method to study the QCD phase diagram, we use the Ultrarelativistic Quantum Molecular Dynamics (UrQMD) transport model [33,34] (version 3.4) to simulate such collisions. For two decades the UrQMD [33,34] is successfully used to predict and describe the A+A collision experiments [35,36]. Although the UrQMD is widely used for RHIC BNL collision energies, its usage for A+A collisions at the LHC CERN energy √ s N N = 2.76 TeV should be made with some care [35]. Nevertheless, as shown in Ref. [37] the UrQMD is able to rather accurately describe the pseudorapidity distribution of charged particles d N ch dη for the pseudorapidity values |η| ≤ 1 and |η| > 4 for 0-5% and 5-10% of most central A+A collisions (see a discussion of Fig. 6 in [37]). Based on this result we studied the Pb+Pb collisions with a zero impact parameter b = 0 at √ s N N = 2.76 TeV using the UrQMD 3.4 to model the transport of hadrons and strings [33][34][35]. This setup is justified by the fact that in the present work we would like to compare the initial baryonic charge density and the bulk properties of hadron production achieved in the usual A+A collisions and in the TNC, while the analysis of phase transformations occurring during a collision process we reserve for the future explorations.
As one can see from Table 1 the UrQMD results for b = 0 found at midrapidity are very accurate for yields of protons, are within 2σ deviation for antiprotons and negative kaons, while for positive kaons they are less accurate. The largest deviation of about 30% one finds for the pions. In our opinion, this is rather good coincidence between the data and model, since we used a vanishing impact parameters, while the ALICE data correspond to the 0-5% bin of centrality class. Note also that the deviation of calculated hadronic yields from the measured ones should not create a problem, since this is a common feature of the A+A collisions and the TNC and, hence, there is a good hope that in the ratios of their spectra and their integrals such deviations will not play an important role or will be strongly reduced. Figure 1 shows that the pseudorapidity distribution of charged particles at |η| < 1 is overestimated by the UrQMD by about 17%. Keeping these result in mind, we conclude that the UrQMD is able to reasonably well reproduce the bulk properties (amplitude and shape) of hadron η distributions of the data measured in 0-5% most central Pb+Pb collisions at √ s N N = 2.76 TeV. The quality of A+A collision description at √ s N N = 200 GeV is similar [35]. Probably, one can improve the UrQMD setup to make a better agreement with the data. However, we strongly believe that such an accuracy is sufficient for the present work which aims to demonstrate to the HENP community the new abilities of the TNC method. Now we turn to the TNC results (see Fig. 4 of Supplemental Materials for some examples and estimates for the TNC rate). Generally, simultaneous collision of three nuclei is almost impossible. In reality, the sequential TNC might happen, i.e. (A+A)+A and A+(A+A) collisions when one beam nucleus collides with a fixed target with a following creation of a fireball that interacts with the second beam nucleus after some time delay. To calculate the observables one should simulate sequential TNC with different time delays and then average over a delay time distribution. This effect should somehow smear particle distributions and reduce magnitudes of signals in comparison to the ideal TNC case. However, treatment of realistic cases is much more complicated and, hence, we leave this analysis for a future research. In this work, we will focus only on the absolutely central TNC with a simultaneous collision of two nuclei from beams with a target nucleus or, in the other words, the central simultaneous A+A+A events. Although, such collisions are idealized, apparently, they will provide an upper estimate of yields of secondary hadrons compared to the most central A+A collisions.
To adapt the UrQMD model for the simultaneous central TNC simulations we use two colliding nuclei with zero impact parameter. The initial distance between these nuclei is increased to be sufficient to fit the third nucleus between them. Afterwards the third nucleus, with the centre in the origin of coordinate system, was added between two colliding nuclei using standard UrQMD routine for nuclei initialization. also see that at mid rapidity the yields of charged pions and positive kaons sizably enhanced, while for the LHC energy this effect for all mesons is very weak. These result show one that one possible way to detect the TNC collisions is to look for a simultaneous and enhanced yields of protons and Λ-hyperons. In Fig. 3 we demonstrate the 3-to-2 nuclei enhancement factor of hadronic transversal momentum spectra in the TNC. Interestingly, there is an entirely different p T behavior of this 3-to-2 nuclei enhancement factor at RHIC and LHC energies! At RHIC energy the momentum spectra of hadrons are enhanced by a p T -independent factor (at least for p T < 2 GeV). At the same time for the LHC energy of collisions, this scale factor is p T -dependent and it decreases, if p T increases. Moreover, as one can see from the lower panel of Fig. 3 the number of slow-moving mesons, i.e. the ones with p T < 0.4 GeV is enhanced in the TNC compared to the A+A collisions, while the number of faster mesons is suppressed. In other words, it seems that the dense baryonic medium strongly decelerates the charged mesons at midrapidity or acts as a kind of density trap for mesons. This shows us that the hydrodynamic flow patterns expected at LHC energies in TNC may be highly nontrivial. Moreover, this figure shows that the effect is strongest for positive kaons and, hence, it In this case, it is possible that in some events the clusters of mesonic matter can be sufficiently large to make a condensate. However, it is also possible that one has to take into account an enhanced number of slow-moving Λ-hyperons and, hence, one can alternatively expect not a mesonic con- densate, but a formation of strange matter droplets or even the strangelets [40,41] We hope that the Hanbury-Brown-Twiss interferometry will be able to verify or disprove these hypotheses. Moreover, we believe that the conditions of the possible formation of dense and slow-moving strange matter in the TNC, i.e. appearance of density traps, should be studied further at higher and lower collision energies using more specialized transport models than the UrQMD. Perhaps, the 3-fluid hydro model [42] and its successor known as THESEUS [43] are the best candidates for such a study.
In Fig. 4 we show the time evolution of the baryonic charge density in the central cell of the volume 27 fm 3 . We analyzed the range of cell volumes from 1 fm 3 to 27 fm 3 and found that for these values the results become volume independent for the central cell volume above 6 fm 3 . Note that the time  In other words, the initial baryonic charge density in the central cell in the case of TNC is about 3 times higher than for the A+A collisions. However, the energy density achieved in the TNC is practically the same as in the A+A collisions.
To quantify the increase of baryonic chemical potential achieved in the TNC we employ the MIT Bag Model equation of state (EoS) [44] for 3 quark flavors and 3 colors. For massless quarks, its pressure as the function of system temperature T and baryonic chemical potential μ B can be written as where the vacuum pressure B vac was chosen as B 1 4 vac = 206 MeV [45]. This value of B vac is a conservative estimate which is widely accepted [46].
From the pressure (1) one can find the baryonic charge density ρ B , entropy density s, and energy density using the thermodynamic identities Equating ρ B M B and B M to the quantities found in central cell for the TNC, one can find T (t) and μ B (t) at different times t of evolution. The results for T (t) and μ B (t) are shown in Fig. 5 both for the TNC (empty symbols) and for the A+A collisions (filled symbols). In Fig. 5 we demonstrate all solutions found by this procedure, but it is apparent that one cannot trust them below the pseudocritical curve (solid curve in Fig. 5) taken from the recent lattice QCD work [47].
For the sake of comparison, in Fig. 5 we show also the parameters of chemical freeze-out in A+A central collisions found by the most advanced version of the hadron resonances gas model based on the induced surface tension EoS [48][49][50]. Within the induced surface tension approach these parameters for the collision energy √ s N N = 2.76 TeV were first found in Ref. [48], while the points for MeV for such reactions at √ s N N = 200 GeV. In our opinion these are quite reasonable estimates keeping in mind the simplicity of the MIT Bag Model and several assumptions adopted here.
Although the main purpose of Fig. 5 is to quantify the increase of baryonic chemical potentials achieved in the course of TNC compared to the A+A collisions, we would like also to point out a peculiar property, which may be a pure coincidence. From  Prediction for the √ s N N dependence of the K + π + A A A ratio for the central-simultaneous TNC based on existing data from E866/E917 [52], NA49 [53,54], STAR [55][56][57] and ALICE [38] collaborations. Scaled SPS data are represented by filled crosses when scaled AGS, RHIC, and LHC data are represented by empty symbols. The found maxima are discussed in the text It is easy to see that for the Bag Model pressure (1) the equation of isentropic expansion s B M /ρ B M B = const is equivalent to the condition μ B /T = const. Inspecting  Fig. 5, one can see that at the moment of crossing the pseudocritical curve [47] the trajectory of the central cell evolution is very close to the condition μ B /T = const, i.e. to an isentropic expansion. Apparently, this feature of the TNC should be also studied in more detail along with the mechanisms of entropy generation in such collisions.
Also, we want to mention that with the TNC method one can vary μ B only changing the central nucleus at constant collision energy. This could lead us to the new method of study QCD phase diagram with at different μ B but almost constant T values.

New abilities and new challenges of the TNC
In the preceding section, we discussed the new opportunities of the TNC method at very high collision energies. Keeping, however, in mind the fact that the RHIC can, in principle, operate in a colliding mode of two nuclei at lower energies and that the NICA accelerator at JINR [51] is a low energy collider as well, in Fig. 6 we summarize the collision energy dependence of the 3-to-2 nuclei enhancement factor for the central-simultaneous TNC. From this figure, one can deduce that, most probably, the effects of proton, Λ-hyperon and positive kaon production will be even stronger at √ s N N ∈ [10; 40] GeV. According to the results of Refs. [8][9][10][11][12] this is a region just above the A+A collision energy threshold of the deconfinement PT. Due to the importance of this energy range, we prepared the predictions for the K + π + A A A ratio collision energy dependence which is expected in the TNC. Using the 3-to-2 nuclei enhancement factors shown in Fig. 6 and the well-known experimental data for K + π + A A ratio measured in the A+A collisions, we constructed Fig. 7 according to the The main question that puzzles one is that the SPS and RHIC data extrapolation to the TNC lead to the different location of the K + π + A A A peak (see Fig. 7). This figure also demonstrates the crucial importance to have accurate experimental data for analysis and for reliable theoretical predictions. The point is that according to Refs. [8][9][10][11][12] the location of K + π + A A A peak, either at √ s N N = 7.7 GeV (prediction based on SPS data, shown by a cross inside the circle in Fig. 7) or at √ s N N = 9.2 GeV (prediction based on RHIC data, shown by a cross inside the square in Fig. 7), can lead to entirely different interpretations of the physics case: either the mass of strange quark carrier starts to change at 7.7 GeV, i.e. after the chiral symmetry restoration PT in hadronic phase and before the deconfinement, or the peak at 9.2 GeV can be considered as a new signal of deconfinement PT, or there may exist two peaks of this ratio. Hence, it would be important to measure this ratio in the TNC in order to find out an exact location of its peak, to verify our prediction, and to determine what set of A+A data, either SPS data or RHIC ones, is compatible with the TNC results. The existing theoretical models for the EoS of dense QCD phases can also make their predictions for the outcome of the TNC prior to the experiments. Then after experiments, the HENP community can easily see what EoS is the most realistic one. It is clear that the TNC method not only opens new research opportunities, but at the same time it creates new challenges for the whole row of hot topics discussed presently by the HENP community. Among the most prospective hot topics, we can name the jet tomography [58,59] in a dense and quark rich medium (compared to A+A collisions) created in the TNC, the chiral magnetic effect [60], and the chiral vortical effect [61,62] for the non-central TNC. Also, the non-central TNC may shed a new light on the polarization phenomenon of (anti)Λ-hyperons [63,64] inside the baryonic charge rich medium. At this stage of research, one can only guess about the outcome of the TNC method for resolving the puzzles of formation of nuclear clusters in nuclear collisions recently discussed in [50,[65][66][67][68][69][70][71][72].
Among the other challenges, it is necessary to mention the numerical modeling of the TNC by the other transport models and by the hydrodynamic ones. Our educated guess is that at low collision energies the TNC will finally require to use the correct cut-off formula [73] for the emission of secondary particles which, in contrast to the approximative Cooper-Frye formula, does not have the contributions of negative number of particles. Moreover, in this case for reliable modeling of the TNC at low collision energies, it will be also necessary to reformulate the equations of relativistic hydrodynamics and hydro-kinetics according to the correct equations of energy-momentum and charge conservations [75,76] which are mathematically consistent with the boundary (kinetic freeze-out) conditions. Besides, an experimental detection and theoretical modeling of the elliptic, triangular and higher order flow patterns [77,78] for the TNC may be another challenge for the HENP community.
We believe that the TNC method can help to ultimately resolve the problem of creating the systems with electric charges Q which are essentially larger than the critical one Z cr +173 units of elementary charge e in heavy ion collisions. This problems has an extended history [79], but with somewhat unclear output. Due to its fundamental aspects, namely the creation of an unstable vacuum and spontaneous emission of positron-electron pairs from it, this problem still attracts a high attention of theoreticians [80]. The TNC method has some advantages over the other ones, since the total electric charge of three lead nuclei is Z PbPbPb = 246, i.e. it is about 34% larger than the electric charge of two colliding uranium nuclei Z UU = 184. Moreover, the effects of Lorentz contraction of two colliding nuclei can essentially increase the electric field even for not very relativistic energies. This gives us a hope that in the TNC the spontaneous emission of positron-electron pairs can be essentially higher and, hence, the ratio of signal to the background, in this case, can be more favorable than in the A+A collisions.
Ab initio simulations of the electromagnetic field created in the A+A collision is a very complicated task [81]. Therefore, for the TNC we made a simplified model and compared it with the results obtained from the instantaneous nucleon density distributions found from the UrQMD 3.4 simulations.
In Fig. 8 we demonstrate the evolution of the apparent charge Q app (y, t) Z cr e = 4π E y (y) Z cr e y 2 , for y ≥ R Pb , seeing at the time t on the distance |y| along the y-axis which is perpendicular to the beam axis z. Here R Pb denotes the radius of the lead nucleus. The center of the coordinate system is located at the geometrical center of the target nucleus.
The y-component of the electric field strength E y (y, t) in the TNC was found from the velocity and spacial distributions of electric charges by the UrQMD. Knowing these distributions one can calculate the electric field strength E y (y, t) and the apparent charge (4). A simple analytical model (its results are shown in Fig. 8 by the short dashed curve) is presented in the Supplementary Materials. It is used to a quick extrapolation of the UrQMD results for the Q app to later times in order to save the CPU time. Also, it provides a simple analytical parameterization, which we found to be reasonably good. From the time dependence of Q app shown in Fig. 8 one can see that during the course of the TNC the apparent charge increases to a maximal value which for the γfactor of projectile nuclei being γ = 2.5 can be as large as max Q app (|y| = 2R Pb , t) 2.5 Z cr e 432.5e! This is an enormous increase compared to the A+A collisions. Moreover, in contrast to the ordinary static case in the TNC, the charge Q app (y, t) increases with the distance from the center of the target nucleus. This means that a strong electric field can produce the positron-electron pairs in a larger volume which, hopefully, may lead to a stronger emission of such pairs. Actually, if this is the case and the number of spontaneously emitted pairs is large, we suggest to use the methods of Hanbury-Brown-Twiss interferometry to these pairs and to study the positron-positron, positron-electron, or electron-electron correlations in order to find out the spacetime picture of their emission process in the TNC. Apparently, the collisions of three uranium nuclei with the total electric charge of 276 elementary ones would be even more interesting.
Above we considered one extreme, namely the triple collisions of heavy atomic nuclei, but it is clear that the other extreme, i.e. the triple collisions of stable hadrons ( p + p + p or p + n + p or even p + p +p ect) are also possible to perform. Of course, the proton target, i.e. gaseous hydrogen, can be easily made, while for the neutron target one better uses the deuterium target. Since the protons are composite objects their triple, but sequential collisions can give us a unique opportunity to study the interaction of partonic medium created in a collision of two protons with the third colliding proton. This may help us to find out the new features of the many-parton distribution functions inside the proton in the new regime. This is another hot topic of modern high energy physics discussed recently in [82].
Also, it would be interesting to exploit the fixed target regime in the electron-positron colliders. Lepton colliders collide fundamental particles, hence the initial state of each event is known, and, moreover, compared to hadronic colliders, the measurements with higher precision can be achieved. Using the existing electron-positron colliders and the new ones, like the Compact Linear Collider (CLIC), it will be possible to study the interaction of products of electron-positron collisions with nuclei of a target. The other applications of the TNC method are briefly discussed in the Supplemental Materials.

Conclusions and perspectives
In this work, we gave convincing arguments that the TNC method will allow the HENP community to create a new frontier in the studies of the QCD phase diagram. This new frontier will make a strategic change in the present situation from the status "still hard to win" to the status "success guaranteed". Our conclusion is based on a result that a strong enhancement of baryon yields in the TNC is found. Predictions for enhancement factors are made for various current and future collider's experiments. Some new features of electron-positron pair production in the TNC are underlined.
However, it is clear that the detection of TNC is highly non-trivial and seems unattainable at the current stage of technology development. Nevertheless, in this article, we propose the theoretical concept of this new method for studying the properties of QCD matter.