Direct Photon Production in High-Energy Heavy Ion Collisions within the Integrated Hydrokinetic Model

The results on description of direct photon yields, transverse momentum spectra, and flow harmonics, measured in ultrarelativistic heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) for different collision centrality classes, analyzed within the Integrated Hydrokinetic Model (iHKM) are reviewed. The iHKM simulation results, corresponding to the two opposite approaches to the matter evolution treatment at the final stage of the system’s expansion within the model, namely, the chemically equilibrated and the chemically frozen evolution, are compared. The so-called “direct photon puzzle” is addressed, and its possible solution, suggesting the account for additional photon emission at confinement, is considered.


Introduction
The analysis of the experimental data on high-energy nucleus-nucleus collisions allows one to conclude that a new form of strongly interacting matter, extremely hot and dense one, is formed in these processes. It is believed that quarks and gluons in this new state of matter are not confined in hadrons, so that such a state is called Quark-Gluon Plasma (QGP). The experimental and theoretical study of QGP constitutes a major task for the present research programs in the field of high-density energy physics.
In particular, QGP is the only source of exotic particles, containing strange, or even heavier charm and bottom, quarks. These quarks are not included in protons and neutrons that make up usual matter we normally see around us. They can be produced only in deconfined state of matter created in relativistic A+A collisions, when high kinetic energy of colliding nuclei transforms, inter alia, into the mass of newly formed heavy quarks. The analysis of exotic particle yields, their ratios and other relevant observables measured in the experiments allow one to investigate the deep properties of QGP. On the other hand, the study of peculiarities of the Quark-Gluon matter evolution process and its eventual breakup into particles, as well as the development of realistic models of these processes facilitates the description and the interpretation of data related to the heavy flavor production.
The Integrated Hydrokinetic Model (iHKM) of high-energy heavy-ion collisions [1,2] repeatedly showed good results in simultaneous description/prediction of a variety of hadron bulk observables for different colliding nuclei at different relativistic energies [2][3][4][5][6][7]. These results particularly included the successful description of strange mesons (K, K * (892), φ(1020)) and strange baryons (Λ, Ω, Ξ) production, as well as kaon and pion femtoscopy. Thus, one can say that more or less realistic picture of relativistic A+A collision was eventually constructed within the iHKM model. The aspect which was missing in this picture was the description of experimental data on photon production and solving the so-called "direct photon puzzle".
The latter task should be considered quite important, since photon radiation produced in relativistic heavy-ion collisions is a unique probe for exploring the properties of the created system and the dynamics of its evolution, since photons are emitted at all the stages of the matter expansion and do not strongly interact with it, thus providing the undistorted information about the system's state at the moment of photon emission for each class of the produced photons [8][9][10][11].
All the photons radiated during the matter evolution are usually divided into decay photons, i.e., those coming from hadron decays (e.g., from π 0 , ω, and η meson decays) [12], and direct photons, including all the other photons apart from the decay ones. The direct photons, in turn, are subdivided into prompt photons, i.e., the photons emitted at the very initial stage of the two nuclei collision as a result of hard partonic scatterings and jet fragmentations, pre-equilibrium photons, coming from the hot guark-gluon matter at the early stage of its evolution, while it has not yet reached a nearly thermal and locallyequilibrated state (this includes also the glasma phase [13][14][15][16][17]), and thermal photons, emitted from the expanding nearly thermalized continuous Quark-Gluon medium [18], as well as from the hadron-resonance gas formed at the late stage of the collision, when continuous QGP gets finally transformed into particles [19]. One can also consider such direct photon sources as jet-photon conversion and jet-induced bremsstrahlung [20,21], involving the interaction of hard partons with the thermal QGP medium.
The thermal direct photons are of particular interest, since they can give us an insight into spatio-temporal structure, collective flow and temperature of the evolving system (however, extraction of this information is somewhat complicated due to superposition of the contributions from different temporal stages of the matter evolution in the resulting spectrum). The thermal photons are dominant in the low-p T part of the spectrum, where they form a nearly exponential p T distribution. As for the prompt photons, they can help us test the correctness of perturbative Quantum ChromoDynamics (pQCD) calculations and investigate the parton distribution functions of the colliding nuclei. The prompt photons play a major role in high-p T region of the spectrum, demonstrating a power law behavior.
Some time ago the experimental results on direct-photon production obtained for Au+Au collisions at the top RHIC energy, √ s NN = 200 GeV, and for Pb+Pb collisions at the LHC energy √ s NN = 2.76 TeV were presented in the articles by the PHENIX Collaboration [22][23][24][25] and the ALICE Collaboration [26][27][28], respectively. The results for both energies demonstrated a similar distinctive feature, namely, the unexpectedly large yields of thermal photons combined with the large elliptic flow v 2 . The first attempts to describe such a peculiarity of the data within hydrodynamic/hybrid models did not succeed: the models either described the measured p T spectra, but underestimated v 2 , or they gave large enough v 2 , however the corresponding photon yields were too low [16,[29][30][31][32]. This motivated new intensive research activity aiming to solve the arisen "direct photon flow puzzle" [13,[33][34][35].
In particular, additional mechanisms of direct photon production were proposed to be taken into account. For example, the photon yields can be increased due to the conformal anomaly in the presence of a very strong magnetic field created in non-central collisions [36,37] or due to synchrotron radiation [38,39]. One more mechanism, allowing to enhance direct-photon anisotropy is "magnetic bremsstrahlung-like radiation" of quarks in the collective color field ensuring confinement [40]. The other possible extra sources of photon radiation, associated with the confinement of QGP medium, are described in [41][42][43][44]. However, although these mechanisms look quite realistic from a theoretical point of view, accounting for them has not lead to fairly clear phenomenological predictions.
Recently, the analysis of direct-photon production in heavy-ion collisions at the mentioned RHIC and LHC energies, was also performed within the integrated hydrokinetic model [45,46]. The calculations of thermal photon emission within iHKM were carried out in the two different approaches to the description of the matter evolution at the posthydrodynamic stage-the chemically equilibrated and the chemically frozen ones. A new possible contribution to the photon spectra connected with the hadronization process was also proposed aiming to improve the observables description and, if possible, to resolve the direct photon puzzle. In this review paper we summarize and re-examine the results obtained in the original works [45,46].

Prompt Photons
The prompt photons are radiated at the earliest stage of the nucleus-nucleus interaction. They are mainly produced in QCD Compton scattering and qq annihilation processes and, additionally, in processes of QCD jets fragmentation. Naturally, the corresponding γ spectra are calculated within the perturbative QCD. In order to obtain prompt photon spectra for A+A collision case, one can take into account the experimental results presented in References [47,48], which indicate that the prompt photon momentum spectra scale with the number of binary nucleon-nucleon collisions, N coll . Thus, to calculate the spectra for Au+Au and Pb+Pb collisions one can rescale the pQCD result for p + p collision case with the corresponding binary collision number N coll , which can be estimated using a Monte Carlo Glauber approach [49].
Generally, the prompt photon production cross-section in p + p collisions can be expressed as follows where the summation is performed over all the different processes of partonic interactions, producing prompt photons, f i and f j are parton distribution functions (PDFs), D γ k is the parton-to-photon fragmentation function and dσ is the corresponding partonic process cross-section. The latter is calculated based on the perturbative expansion with respect to the strong coupling constant, α s (Q ren ).
As one can see, in order to calculate the cross-section dσ pp , one should specify the three QCD scales, Q f act , Q ren , and Q f rag , entering (1). In our consideration the values of all the three scales are fixed to Q = 0.5 p T , since it is commonly known that the smaller is the constant of proportionality between Q and p T , the better data description can be reached.
To calculate the prompt γ transverse momentum spectra in case of proton-proton collisions we utilize the JETPHOX code [50]. At high photon p T (perturbative region) one can use the optimal relation Q = 0.5 p T . However, one can face the problem, when applying the mentioned parametrization for relatively small photon momenta, p T < Q 0 /λ GeV/c, where λ is the proportionality coefficient between Q and p T , and Q 0 ≈ 1.5 GeV/c is the energy-momentum scale usually considered as the limit of pQCD applicability. Fortunately, it can be shown that changing the λ value results simply in the renormalization of the photon spectrum [13]. So, the low-p T part of the spectrum can be obtained by rescaling the spectrum, calculated at larger λ value, e.g., at Q = 4 p T , in accordance with the high-p T results for the case of Q = 0.5 p T .
In our study, the EPS09 parton distribution function [51] and the BFG II fragmentation function [52] were used in tabular form for prompt photon spectra calculations.

Pre-Equilibrium and Thermal Photons
The pre-thermal and thermal photon contributions are treated in our consideration based on the matter evolution description, obtained within the iHKM model. The model simulates the full process of the evolution of the system formed in a high-energy nuclear collision, which can be divided into several stages, each described using a relevant approach (see our previous papers, e.g., [1,2,46] for a detailed model description).
The initial state at the proper time τ 0 , serving a starting point for the subsequent pre-equilibrium energy-momentum tensor dynamics within a relaxation-time approximation [1,53], is modeled with the help of GLISSANDO code [54][55][56]. It gives us the initial energy-density transverse profile (r T ), containing contributions from the binary collision and the wounded nucleon models, which is then scaled by the factor 0 -iHKM parameter, characterizing the initial energy density in the center of the system. One more iHKM parameter, α defines the fraction of binary collisions in the resulting (r T ) distribution. During the pre-thermal expansion system gradually equilibrates, so that before the thermalization time τ th ≈ 1 fm/c, it becomes nearly thermal and thus can be further described within the viscous hydrodynamics approximation.
At the next, hydrodynamic stage the system expansion is described using the Israel-Stewart formalism with the Laine-Schroeder equation of state (EoS) for the Quark-Gluon matter [57]. We also use the minimal possible shear viscosity to entropy density ratio η/s = 0.08 ≈ 1/4π in our simulations. So, the system gradually cools down, and when the temperature reaches the particlization value, T p = 165 MeV, one performs sharp switching to the description in terms of particles based on the Cooper-Frye prescription (the particlization hypersurface is constructed using the Cornelius routine [58]). Finally, the generated particles are forwarded to the UrQMD hadron cascade [59,60], which describes the afterburner stage of the system's evolution.
The last sentence corresponds to the original iHKM model used for the description of hadronic observables. However, the description of photon emission at the final stage of matter evolution presently cannot be explicitly implemented in UrQMD, so that we use two alternative approaches of γ radiation treatment at the post-hydrodynamic stage.
The first variant is to continue the hydrodynamic evolution of the hadron matter in a chemically equilibrated way until the photon radiation becomes negligible, e.g., until the temperature drops to the values lower than T ≈ 100 MeV. In this case, one can determine the hydrodynamic velocities u µ (x) and the temperatures T(x) at the final stage of collision straightforwardly, so that the thermal photon production from the expanding hadronic medium can be easily calculated based on the well-known results for the emission of photons from a thermal hadronic system in rest. Such a method has a noticeable disadvantage, namely, the fact that at low enough temperatures (still relevant to the evolution of considered hadronic system) the chemical equilibrium gets definitely broken.
Another way of dealing with the afterburner stage suggests using the original version of hydrokinetic model-HKM [61]. The latter features chemically frozen matter evolution with continuous particlization, in course of which the local chemical equilibrium, as well as the local thermal one, are violated. The chemical composition of the expanding system does not change due to inelastic reactions, which are all switched off, except for the resonance decays. Such a picture of the late stage of the collision corresponds to the widespread conception of chemical freeze-out. It requires the introduction of separate chemical potential for each sort of hadron, thus complicating the simulation of the system's evolution as compared to the chemical equilibrium case, but still allowing to calculate all the quantities (collective velocities, temperatures, chemical potentials) necessary for estimation of the photon production. The latter, in principle, should be somewhat different from the case of chemically equilibrated evolution.
In our analysis we apply both ways of describing the evolution of hadronic system in order to compare them and get the estimation for the lower and the upper limits of the calculated direct-photon spectra and v 2 coefficients resulting from the use of such rough approximations.
In what follows we will describe how we account for the thermal photons contribution to the resulting direct-photon spectrum. However, in our study we account for the prethermal photons as well (despite the fact that, according to our estimates, their contribution to the resulting γ spectrum is only about 2% of the total thermal photon contribution), and the way we do it is quite similar to that we use for the thermal photons. The total energy-momentum tensor for the iHKM pre-thermal stage is assumed to be a sum of two parts: the hydrodynamic tensor T µν hydro (x), multiplied by the weight (1 − P (τ)), and the free-streaming tensor T µν f ree (x), multiplied by P (τ) (see [1,2] for details). Accordingly, we treat pre-thermal photons as thermal ones, produced solely by the hydrodynamic component of the total energy-momentum tensor at the pre-thermal phase of the matter evolution, i.e., the corresponding spectrum contribution is calculated the same way as for the hydrodynamics regime, but in addition it is scaled by the factor (1 − P (τ)).
The momentum spectra of the thermal photons, coming from the QGP phase, can be described well already in the leading order of strong coupling constant α s , see Reference [62].
Hence, we can use the corresponding formulas for the spectra calculations in our analysis within iHKM. The QGP thermal photon emission originates mainly from the leading order 2 → 2 processes (however, as it is demonstrated in [63,64], the higher-order processes can also bring noticeable contribution of parametrically the same order to the resulting spectrum of thermal photons). The photon radiation rate can be expressed by the following relativistically invariant formula where k is the photon momentum and u is the hydrodynamic velocity in the point x.
Accordingly, the quantity ν e (k · u) can be interpreted as the spontaneous emission rate of the photons with the momentum k. In Reference [62] the rate ν e (k · u) was evaluated by summarizing contributions from different partonic 2 → 2 processes. In addition, the radiation from inelastic qq annihilation and near-collinear bremsstrahlung processes, as well as the effects of Landau-Pomeranchuk-Migdal suppression were taken into account. As a result, the total spectrum can be written as follows: where ∆ 4 V = τ∆τ∆x∆y∆η is the volume of a 3D computation grid cell and the summation runs over all the grid cells, where T > 165 MeV, thus comprising the pre-thermal and the hydrodynamic iHKM stages. It is worth noting, that in the Formula (3) for the thermal photon spectra we neglect the corrections to the continuous Quark-Gluon matter PDFs, arising due to viscosity, since at the hydrodynamic stage of the system's expansion we assume the shear viscosity to entropy density ratio, η/s, to have the minimal possible value. At the post-hydrodynamic (hadronic) stage of the matter evolution the same Expression (3) can be applied to estimate the thermal photon emission, however the quantity ν e (k · u) should now correspond to the radiation from hadron-resonance gas. This radiation includes many different sources that should be taken into account:

•
The photons radiating from a meson gas. Here we base our consideration on the Reference [65]. We account for the hadronic reactions with π, K and ρ mesons and K * resonances, where photons are produced. The corresponding emission rates are taken from [65] with one additional prescription for the chemically frozen evolution regime: in this case the contribution from each hadronic reaction has to be multiplied by , where µ i (τ, x) represent the individual chemical potentials of hadrons participating in the reaction. • The photons coming from in-medium ρ mesons. Here we follow the results of the Reference [66]. The emission rate for these photons incorporates baryon (antibaryon) effects and in our analysis depends on the local baryon chemical potential µ B (τ, x) (for chemical equilibrium expansion), and, in addition, on individual chemical potentials of hadrons involved in reaction (for chemically frozen evolution). • The photons emitted in reactions with the three mesons π, ρ and ω [19], such as πρ → ωγ, πω → ργ, and ρω → πγ. Again, for the chemically frozen regime the corresponding contributions should be multiplied by exp ∑ µ i (τ,x)

Hadronization Photons
In the work [45] an approach was developed allowing to account for the extra portion of soft photons, which presumably should be emitted from the system during the process of Quark-Gluon matter confinement (in our study within iHKM it corresponds to an extra γ emission at the particlization phase of the system's evolution). Such additional electromagnetic radiation could take place, e.g., via the mechanism of "magnetic bremsstrahlung-like" (synchrotron) radiation, coming from the interaction of quarks, which escape from QGP, with the confining collective color field [40]. Other possible mechanisms of additional photon emission are considered in the papers [41][42][43][44]67]. Despite some of these mechanisms can improve calculated direct-photon spectra significantly, the corresponding models are still phenomenological and incomplete, since the considered processes are non-perturbative and non-equilibrium. In order to avoid uncertainties associated with the use of different theoretical models and their different calibration, in our approach we aim to implement the mechanism of confinement emission in a simplest phenomenological way possible.
Let us introduce the emission function of extra photon radiation at the particlization (hadronization) stage, G hadr (t, r, p), as follows: Then let us also consider the hypersurface of maximal emission for the photons possessing 4-momentum p. For each given p such a hypersurface σ can be constructed as a set of corresponding space-time points (t σ , r), where t σ = t σ (r, p). Following the recipe utilized in [68,69], one can switch to new space coordinates x = r + p p 0 t σ (r, p), including the saddle point t σ (r, p). The emission function G hadr (t, x, p) after that can be approximately presented (using saddle point method) as G hadr ≈ F(t, x, p) exp(− (t−t σ ) 2 2D 2 ), where F should be a function, depending smoothly on time t. The application of such an approximation gives one the following spectrum: . (5) If we then assume that the considered hypersurface of maximal emission coincides for soft photons with the common hypersurface of hadronization (in iHKM we use the isotherm T h = 165 MeV as such hadronization hypersurface), the previous equation can be rewritten as follows: where the chemical potentials at T = T h are assumed to be vanishing in case of the LHC high-energy collisions, and the Heaviside function θ(x) is introduced to cut off possible negative contributions to the momentum spectra from non-space-like pieces of the particlization hypersurface. Note, that the quantity D c (representing the emission temporal width in the obtained Expression (6) for momentum spectra) is assumed to depend on collision centrality c, reflecting the fact that the created fireball expands and hadronizes faster for non-central collisions due to higher gradients in transverse direction, caused by smaller transverse size of the system in this case. Let us further assume, following the work [44] and our articles [45,46], that the photon emission function, approximated in our consideration by the product of the two functions F · D c , can be presented in the following thermal-like form: with the factor g = 2.
The quantity γ hadr in the above formula characterizes the QGP confinement process itself (its basic properties), while the parameter d c ∝ D characterizes its temporal width and thus can be possibly different for different centrality classes. However, in our analysis, carried out in the References [45,46], we combined the product of these two parameters, γ hadr and d c , into a single parameter β ≡ d c γ hadr and used the same its value to describe the hadronization emission of photons for all the considered centrality classes at given collision energy. The value β = 0.04 was used in the case of Au+Au collisions at the top RHIC energy [46] and the value β = 0.02 was used even earlier for 2.76A TeV Pb+Pb collisions at the LHC [45]. Note, however, that the photon momentum spectra measured by the STAR Collaboration at RHIC [70] can also be described in iHKM with the "LHC value", β = 0.02, and re-tuning to β = 0.04 is necessary to describe at the same time the data on elliptic and triangular flows, measured by the PHENIX Collaboration [71].
The possible reason for the observed distinction between the β values used for RHIC and for LHC can be connected with the corresponding difference between the maximal intensities of the hadronization process at the two colliders. The hadronization intensity most likely should be smaller at the LHC, because of much higher system's expansion rate in this case, which would impede the quark coalescence. The detailed study of the interplay between different factors affecting the β parameter value is, however, outside the scope of current paper and can be carried out in further works.

Results and Discussion
In the papers [45,46] to describe the photon production at RHIC and the LHC we utilized the same model parameters tuning which previously ensured the successful description of hadron observables at the same collision energies [5,6]. As a result of analysis carried out in [45,46] it was found that the description of direct-photon p T spectra, elliptic and triangular flow p T -dependencies, v 2 (p T ) and v 3 (p T ), can be significantly improved, for both RHIC and the LHC energies, if one includes, along with the contributions from prompt, pre-thermal and thermal (QGP and hadron gas) photon emission, an additional portion of γ radiation, corresponding to the confinement process, which can be called the "hadronization photons" (see Figure 1, demonstrating the RHIC and the LHC photon spectra calculated within the iHKM for the central events (c = 0-20%); the spectra obtained with and without accounting for the hadronization emission (HE) are shown). Comparing the model spectra for the two simulation modes with the experimental points presented by the PHENIX [25] and the ALICE [28] Collaborations, we can definitely say that for both collision energies the inclusion of the HE component makes the simulation results much closer to the measured data.
In iHKM the magnitude of the HE contribution is regulated by the only one free parameter β, which characterizes the intensity of the photon emission during the process of confinement. In order to fix the value of the β parameter, one should at first perform fitting of the direct-photon spectrum for the case of most central events. Having determined in such a way the value of the HE parameter, one can further use it for all the centrality classes for the description of transverse momentum photon spectra, together with elliptic and triangular flows. Our analysis gives the value β = 0.02 for the LHC Pb+Pb collisions at the energy 2.76A TeV and the value β = 0.04 for the RHIC Au+Au collisions at 200A GeV.
As for the dependence of the simulation results on the variant of the matter evolution treatment at the afterburner stage of the collision (chemically equilibrated vs. chemically frozen evolution), the photon p T spectra corresponding to these two approaches are compared in Figure 2. From the plot one can conclude, that both calculation regimes give very close results, although the chemically frozen evolution leads to hardly noticeable lowering of the total spectra. The same concerns the results for elliptic flow coefficients v 2 .  So, despite our assumption that changing the approach to the description of the system's post-hydrodynamic expansion should in theory lead to changes in the resulting photon spectra, we see that on practice the difference is fairly insignificant. One can explain such a result by mutual compensation of particular changes brought to the model simulation by switching to another matter evolution treatment. Specifically, the emergence of additional hadron chemical potentials, which accompanies our passage to the chemically frozen evolution scenario, could be compensated by a simultaneous lowering of the temperatures, associated with certain energy density values. The latter lowering, in turn, could happen due to interplay between the changes in the system's expansion rate and in its chemical composition, taking place when one switches from the chemically equilibrated to the chemically frozen regime (see [61,68]).
An interesting detail concerning the dependency of the total direct-photon yield on the mean charged-particle multiplicity dN ch dη (see Table 1) seems to be worth noting: in our model calculations we obtained that the total yield of soft direct photons is approximately proportional to dN ch dη 1.25 [46]. This feature of the iHKM results is in agreement with the results of experimental measurements of the PHENIX Collaboration [72]. From Table 1 one can also see that for all the centralities the magnitude of the contribution from the hadronization emission photons to the total yield is comparable with that of the thermal photon contribution (which includes the emission from both QGP and hadron matter phases), while the contribution of prompt photons is about two-times smaller. The thermal and hadronization photons, however, dominate only in a relatively-low-momentum region, p T < 3 GeV/c, while for high p T > 3 − 4 GeV/c prompt photon emission gives the main contribution to the total spectrum (see Figure 1 from [45] and Figure 4 from [46]). One can also see our results on direct-photon v 2 (p T ) and v 3 (p T ) dependencies, obtained in iHKM for the Au+Au collisions at the top RHIC energy, c = 20-40%, in Figure 3. The experimental points reported by the PHENIX Collaboration [71] are shown as well for comparison purposes. Again, we see that the hadronization photons play an important role in bettering the description of the experimental data within the model. This, in particular, should imply that the expanding matter situated in a thin space-time layer enclosing the hadronization hypersurface is strongly involved into collective motion. Comparing the results on flow harmonics for different collision centralities (see Figures 6-11 from the paper [46]) one can also conclude that the relative contribution of the hadronization emission to v n values increases when one passes from central events to peripheral ones.
In Figure 4 we also show the results on v 2 (p T ) dependency for the LHC energy 2.76A TeV. In addition to the curves presented for RHIC we include here the elliptic flow coefficients calculated in iHKM for the thermal photons only (excluding the contributions from prompt and hadronization photons). The corresponding curve goes even higher and closer to the experimental data than the curve for the full iHKM calculation with the HE contribution. Such a peculiarity can be explained if one notices that the prompt photon radiation is momentum isotropic, so that its inclusion into resulting v 2 coefficients suppresses the anisotropy, specific for the thermal photon momentum spectra. This effect manifests itself the stronger, the more peripheral collision events we consider [45].

Conclusions
The previously obtained results on direct-photon yields, momentum spectra, v 2 and v 3 flow harmonics obtained within the integrated hydrokinetic model of relativistic nucleusnucleus collisions for the top RHIC energy 200A GeV and the LHC energy 2.76A TeV are summarized and reviewed.
The comparison of the chemically equilibrated and the chemically frozen descriptions of the matter evolution at the post-hadronization stage shows that the chemically frozen expansion treatment leads to a slightly lower resulting p T spectra and v n coefficients of direct photons.
The main conclusion that can be drawn from the performed analysis is that the problem of simultaneous theoretical description of the direct-photon spectra and flow harmonics, which appear to have unexpectedly large values in the experimental data, referred to as the "direct photon flow puzzle", probably can be solved if one accounts for the extra contribution to the photon production connected with the process of confinement (hadronization) in the considered strongly interacting system. The most simple phenomenological accounting for this additional portion of radiation is sufficient to significantly improve the agreement between the model and experimental results.
A satisfactory photon data description was obtained within the iHKM using the same basic model parameters tuning, which was used earlier for the description of hadronic observables at the considered collision energies. The elucidated dynamics creates a reliable basis for the study of other phenomena and also mechanisms of the production of exotic particles, e.g., those containing heavy flavors.

Funding:
The research was carried out within the NAS of Ukraine Targeted research program "Collaboration in advanced international projects on high-energy physics and nuclear physics", Agreement 7-2022 between the NAS of Ukraine and BITP of the NAS of Ukraine.