Holographic Photon Production in Heavy Ion Collisions

The thermal-photon emission from strongly coupled gauge theories at finite temperature is calculated using holographic models for QCD in the Veneziano limit (V-QCD). The emission rates are then embedded in hydrodynamic simulations combined with prompt photons from hard scattering and the thermal photons from hadron gas to analyze the spectra and anisotropic flow of direct photons at RHIC and LHC. The results from different sources responsible for the thermal photons in QGP including the weakly coupled QGP (wQGP) from perturbative calculations, strongly coupled $\mathcal{N}=4$ super Yang-Mills (SYM) plasma (as a benchmark for reference), and Gubser's phenomenological holographic model are then compared. It is found that the direct-photon spectra are enhanced in the strongly coupled scenario compared with the ones in the wQGP, especially at high momenta. Moreover, both the elliptic flow and triangular flow of direct photons are amplified at high momenta for V-QCD and the SYM plasma. The results are further compared with experimental observations.

In perturbative quantum chromodynamics (pQCD), the complete leading-order result of photon production in thermal equilibrium incorporating 2 to 2 scattering and collinear emission has been obtained in [1,2] via the hard-thermal-loop (HTL) resummation, where the perturbative calculation is performed at high temperature with the resummation including dressed propagators and vertices (e.g. the parton self energy computed with soft external momenta as the so called HTL approximation is proportional to g 2 s T 2 .) in loop diagrams [3]. The next-to-leading-order calculation was accomplished in [4].
However, due to the strong coupling of the QGP around the deconfinement temperature, perturbative approaches may become invalid for the study of electromagnetic properties of the QGP in Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC).
Although the electric conductivity of the QGP has been recently obtained from lattice simulations [5,6], the non-perturbative computation of full spectra of photon or dilepton production in sQGP is still formidable.
The gauge/gravity duality, which connects the strongly coupled gauge theories and classical supergravity, provides an alternative way to tackle non-perturbative problems [7][8][9]. It has been widely applied to investigate the heavy ion collisions (see [10] and the references therein). In holography, the study of thermal-photon production and dilepton production from a strongly coupled N = 4 SYM plasma was initiated by [11].
Similar calculations were performed in top-down models that include quarks in the fundamental representation by adding flavor branes to the sYM theory, [12,13]. In addition, the photoemission at intermediate t'Hooft coupling was analyzed in [14,15]. Further studies incorporated strong magnetic fields [16][17][18] and pressure anisotropy [19][20][21], which are relevant to electromagnetic probes generated at the early stage of heavy ion collisions. Moreover, the photon production in the non-equilibrium state in holography were presented in [22][23][24][25].
On the other hand, recent measurements of large anisotropic flow of direct photons comparable to the hadron flow in RHIC [26] and LHC [27] lead to a new puzzle. It is generally believed that electromagnetic probes only record the local information in heavy ion collisions. Due to the weak electromagnetic coupling, they barely interact with the QGP after they are produced. Some of the recent studies indicate that the interaction between electromagnetic probes and the QGP could be enhanced compared with the case in the quantum electrodynamics (QED) plasma [28,29], whereas such an enhanced interaction should be negligible due to short lifetime of the QGP.
In general, since thermal photons inherit the flow from their sources, it is expected that the anisotropic flow of direct photons at intermediate energy around, 1 − 3 GeV, should be smaller than the flow of hadrons given that a large portion of direct photons are generated in early times with small momentum anisotropy. Theoretically, the state-of-the-art simulation incorporating the emission of prompt photons from hard interactions in early times before thermalization and of thermal photons produced in the weakly coupled QGP with HTL resummation and the rate from hadron gas along with the plasma expansion dictated by viscous hydrodynamics, under-predicts the elliptic flow compared with the experimental measurements [30,31].
On the other hand, a variety of non-conventional mechanisms have been proposed to address the puzzle [18,[32][33][34][35][36][37][38][39]. A more recent study, [31], which incorporates the IP-Glasma initial states and the influence from both shear and bulk viscosities with updated photon emission rates, improves the agreement with experimental observations [40]. In addition, the effects from late equilibration of (anti-)quarks in the QGP evolution was investigated in [41], which yield mild suppression of the spectrum but an enhancement of flow for direct photons.
Despite the validity at the RHIC and LHC temperature, the leading-order emission rate of thermal photons from the weakly coupled QGP is applied [1,2] to describe the thermal radiation from the deconfined phase in the simulations due to the difficulty in tackling photon emission from the sQGP.
In this paper we employ two bottom-up holographic models to model the sQGP, both of which break conformal invariance and match several properties of QCD at finite temperature as calculated from lattice simulations. In holography, the running of the coupling can be characterized by a bulk scalar field (dual to the operator driving the RG flow) with a corresponding potential in the gravity dual. This type of models may be regarded as effective theories of QCD in the infrared (IR) regime. Although the validity of such models may not be rigorously justified, they quantitatively agree with many properties of QCD and provide an alternative route to probe the direct-photon problem and improve our understandings for the present tension between theories and experiments.
One of the models we employ was introduced by Gubser and Nellore [42,43]. In [44], the model was utilized to fit the lattice data for the electric charge susceptibility and qualitatively describes the electric conductivity from lattice data near the deconfinement phase transition [45]. Furthermore, the thermal-photon production in such a model was analyzed and was briefly reported in [29]. Hereafter we will call this model of GN model.
The other model we will use is V-QCD. This is a more sophisticated holographic model and it is based on the improved holographic QCD formalism [46][47][48][49]. The flavor degrees of freedom are included by adding N f brane-antibrane pairs with a bulk tachyon field that describes the spontaneous breaking of chiral symmetry [50][51][52][53].
The framework of V-QCD works in the Veneziano limit and includes the backreaction of flavor on color [54]. Its phase structure at zero temperature, finite temperature and finite temperature and density was analyzed in a series of works [55][56][57][58][59][60][61][62][63]. In particular, several transport coefficients and flavor current two point functions were analysed at finite temperature and chemical potential in [60]. V-QCD contains in particular the dynamics of chiral symmetry breaking that is absent from the GN model.
The photon-emission rates from these holographic models are then convoluted with the medium evolution. Furthermore, the contributions from prompt photons and thermal photons from hadron gas are incorporated to compute both the spectra and flow of direct photons in both RHIC and LHC energies. Here, we summarize our findings in short: 1. Holographic photon production results in enhanced spectra of direct-photons especially for direct photons with high transverse (perpendicular to the beam direction in experiments) momenta p T .
2. The flow is suppressed in both holographic models compared to pQCD with HTL resummation at low p T , while V-QCD enhances the flow at high p T .
3. Holographic models improve the agreement with experimental observations in spectra, but still underestimate the flow.
4. In small collision systems, holographic models give rise to enhancements in both spectra and flow.
Based on the findings, we also present our conclusions below.
1. The hadronic contribution for thermal photons plays a central role for direct-photon flow in low p T .
2. In high p T , thermal photons from the QGP phase are responsible for the flow. In this case, the direct-photon flow can be further increased when the emission rate is more prominent at low temperature and the yield is more dominant than prompt-photon contribution for the corresponding model.
3. In small collision systems, larger spectra of thermal photons from the QGP phase simultaneously yield larger flow for direct photons.
The paper is organized as following. In section II, we introduce basic ideas about the constructions of holographic QCD from bottom-up approaches. In section III, we elaborate the holographic models we will apply for photon production. We then show the technical details for the computations of photon emission rates and compare the results from distinct models in section IV. In section V, we first expatiate the setup of the full simulation for the direct-photon production in heavy ion collisions. Then we present the spectra and elliptic flow in RHIC and LHC, in which we further make discussions and comparisons between the results and experimental observations. Finally, we make concluding remarks in section VI.

II. HOLOGRAPHIC QCD
There are several approaches to the holographic study of strongly coupled QCD. The deconfined phase of the theory at intermediate temperatures already bears some similarities with the N = 4 SYM theory. An effort to describe holographically QCD features in more detail, led to the development of phenomenological bottom-up holographic models, [46,47,64,65] and [42].
Current holographic approaches to the QCD problem with light quarks are carried in the context of the probe approximation whereby the light quarks are inserted as spectator D-branes in the limit where N f N c . The Veneziano limit consists of addressing the double limit N f,c → ∞ but fixed λ = g 2 N c and x = N f /N c in the context of holographic models or V-QCD. In the vacuum V-QCD exhibits chiral symmetry breaking for x ≤ 4, restores chiral symmetry for 4 ≤ x ≤ 5.5, [54]. Beyond the Banks-Zaks point or x > 5.5, [66], the theory becomes QED-like.
The bulk action describing the system S g is the gluon action which is the same as the improved-holographic QCD (ihQCD) action, [65]. S f is the brane-antibrane action which was introduced in [67], and was proposed as a low energy effective action for the holographic meson sector in [50][51][52][53]. The last part, S a , is the action of the CP-odd sector which contains the coupling of the flavor-singlet mesons to the axion, which comes from the closed string sector, [50,51,57].
The full form of the first two parts of the action will be presented below. Since, we currently do not study the flavor-singlet CP-odd excitations and the vacuum state is CPeven the last term is not discussed further. To study the excitation spectrum we expand the above equations to quadratic order around the vacuum solution.
A. The glue sector The glue part of QCD at large N c and strong t' Hooft coupling is described holographically by a two-derivative effective action with bulk fields that correspond to the lowest dimension operators of the theory. Those include the metric and the dilaton. According to gauge/gravity duality, the metric and the dilaton are dual to the energy-momentum and the TrF 2 operator, respectively. The following gluonic action is used to describe the dynamics is the 5-dimensional Plank mass and λ = e φ is the exponential of the bulk dilaton with the corresponding potential V g . λ is interpreted as the holographic t' Hooft coupling. The Ansatz for the finite temperature vacuum solution is where the conformal factor e 2A is identified as the energy scale in field theory and f (r) is the black hole factor. The choice of the potential V g (φ) depends on the details of each holographic model. In this work, we will consider three different models in order to compare their results for the photon production in Heavy-Ion collisions. Our models include the prototype AdS-Schwarzschild background, where V g (λ) is constant, the Improved Holographic QCD model, [48,49,65,68], and the GN model, [42]. In the GN model, the U (1) flavor current is described phenomenologically by including a bulk U (1) gauge field, with the standard kinetic term multiplied by a dilaton dependent potential.
In V-QCD, a more thorough approach is followed since the dynamics of the flavor sector (both in the confined and the deconfined phase) is described by introducing appropriate bulk fields for all the low dimension operators of the theory. Those include the quark condensatē ψ R ψ L , as well as the left and right flavor currentsψ L γ µ ψ L andψ R γ µ ψ R .
To describe the dynamics of the above operators and couple them to the glue part of the theory, we effectively consider N f pairs ofD 4 and D 4 coincident branes in the 5 dimensional bulk spacetime. The dynamics of such a a brane system in flat space-time was described in [67]. The lowest lying fields of the system include a complex scalar field T IJ (dual to the quark condensate ) that transforms in the (N f ,N f ) of the flavor group SU (N f ) L ×SU (N f ) R . It is the famous tachyon of brane-antibrane systems in flat space and it also a tachyon in the asymptotically AdS space relevant for holography. The tachyon instability translates into the fact that the quark mass operator is a relevant operator in QFT.
There are also two U (N f ) gauge fields, A L µ and A R µ , dual to the left and right flavor currents. The dynamics of both the tachyon and the gauge fields is determined by the action, which was proposed by Sen. The condensation of the tachyon describes the chiral symmetry breaking in the vacuum of the theory. In the current article, we are interested in the finite temperature, deconfined state of the model in which chiral symmetry is restored. This is described by a black hole that has dilaton hair but no tachyon hair, [55].

III. THE DIFFERENT HOLOGRAPHIC MODELS
The glue part of the action is the same as (2.3) with the dilaton set to zero and V g = 12 2 , where is the AdS radius. The background solution at finite temperature is given by ( In our discussion, we will omit the 5-dimensional sphere which is a part of the background metric because it will not play any role. The flavor part of the action describing the U (1) bulk gauge field dual to the photon current is simply

B. The Gubser and Nellore model
In the case of the GN model, the dilaton potential is chosen in order for the model to reproduce the the equation of state of the Quark Gluon Plasma near T c 1 , [44] and [29].
where γ = 0.606, b 2 = 0.703, b 4 = −0.12 and b 6 = 0.0044 and = 1. The background metric and dilaton are found numerically by integrating Einstein's equations. The flavor sector here contains also a U (1) gauge field dual to the photon which is coupled to gravity by the following action 1 Note the this model is only capable of describing the deconfined phase. At zero temperature it is different from QCD as it is not confining and does not have a mass gap.
where the potential of the coupling of the dilaton to the gauge field is taken to be, [29], where g 2 e , is a dimensionless constant which is matched to the overall amplitude of the relevant observables, and a 1 = 0.4. This is a phenomenological choice of the potential which matches the lattice results for the electric susceptibility of the QGP for T < 1.5T c . In [29], g e is chosen to saturate the electric conductivity of the SYM plasma at T = 3.5T c . In the subsequent study [69] with the inclusion of chemical potentials, a different choice of g e is considered, which results in smaller electric conductivity and photon-emission rates in the zero chemical potential.

C. V-QCD
The bulk fields corresponding to the lowest dimension meson operators include two vector gauge fields and a complex scalar, the tachyon. Their effective holographic action was proposed in [52,53,70], where the definitions of A L/R are where g M N is the induced metric on the brane-antibrane pair. The fields A L , A R as well as the tachyon fieldT are N f × N f matrices in the flavor space. A L/R . The Plank mass, which appears as an overall factor in front of both S g and S f , is fixed by requiring the pressure of the system to approach the large temperature limit of free non-interacting fermions and bosons. This fixes (M ) 3 = (1 + 7x/4)/45π 2 [55], and the AdS radius in the presence of backreacting flavors is 3 = 3 0 1 + 7x 4 . The quarks are taken to have the same mass, so the tachyon field is proportional to the unit matrix in flavor space,T = τ (r)I N f . The field τ (r) is real.
The near-boundary expansion of the tachyon field matches the UV running of t' Hooft coupling and the anomalous dimension of the quark mass operator [54]. The tachyon close to the boundary r → 0 is where the power γ is matched to the coefficients of the anomalous dimension ofqq. Here m q and qq denote the quark mass and quark condensate, respectively. In the IR, the tachyon field diverges and the tachyon potential vanishes as it is argued in [67]. As it is shown in [52,53,70] brane -antibrane condensation in confining backgrounds leads to chiral symmetry breaking. The form of the tachyon potential that we use is The rest of the potentials κ(λ) and w(λ) are taken to be independent of T and have an analytic expansion close to the boundary in terms of λ. Chiral symmetry breaking, thermodynamics and the meson spectra constrain their IR asymptotics, [57] potentials V f 0 (λ), κ(λ) and a(λ) are given by (3.14) It should be noted that the potentials in IR are remarkably close to the non-critical string theory values on flat space-time, except form the logarithmic corrections, see [56][57][58]. In our present work, we choose w(λ), in such a way that the holographic electric conductivity of the QGP in the deconfined phase.
The boundary of the bulk space-time is taken to be at r = 0 and the field-asymptotics are modified from the usual AdS asymptotics such that they reproduce the perturbative running of t' Hooft coupling. Hence in the UV, the scale factor of the metric A and the 't Hooft coupling λ are where is the AdS radius and Λ is the UV scale of the theory. The dilaton potential is chosen such that it reproduces the perturbative running of t' Hooft coupling in the UV. Even if QCD in the UV is not expected to have a dual description in terms of a gravity theory, since it is weakly coupled, the UV asymptotics of our model provide the correct UV-boundary conditions for the IR description of the system. The requirement of confinement, gapped and discrete glueball spectrum as well as linear Regge trajectories fix the IR behavior of the potential to V (λ) ∼ λ 4/3 √ ln λ as λ → ∞. Considering a simple interpolation of the asymptotic forms of the potential in the two regions we obtain where 0 is the AdS radius in case of no backreacting flavors.

D. The vector field action
To calculate correlation functions of the electric current in QGP, we focus on the vector gauge field in the bulk, which is defined as , in terms of the left and right gauge fields, [60]. We expand the action, (3.10), to quadratic order in V M In the chirally symmetric phase of the model the tachyon is zero,T = 0, hence the above potentials depend only on the dilaton. Hence, in this case the coupling function of the dilaton to the gauge field reads In the current paper, we make two different choices of the potential w(λ) and explore its consequences to the final photon spectrum and it's flow.
The two choices are motivated by the exploration of the holographic possibilities. The reason is as follows: Lattice calculations provide the electrical conductivity of QCD in a finite energy window, namely T c < T < 2T c . This constraints directly w(λ) in a finite range of couplings λ c < λ < λ cc . w(λ) for λ > λ c is not constrained from current lattice results.
The two distinct parametrizations above, although both in agreement with lattice results for λ < λ cc are different for λ > λ cc . In the second case, w 2 (λ) increases faster than w 1 (λ), for low values of λ or equivalently high temperature 2 This choice was made in order to explore the consequences of a high UV-contribution in the final photon spectrum. This will be explained in more detail in the following.

E. The vector fluctuations
We calculate the action and the equations of motion for the fluctuation of the vector gauge field in the deconfined background, where chiral symmetry is restored (τ (r) = 0), [60]. The calculation is common for the different models which we consider, since they only differ in the background solutions and the choice of the potential V e (λ). The quadratic action for the fluctuation of the vector field in the V r = 0 gauge is Without loss of generality, we may take k 3 = k = 0 and k 0 = −ω. Then, the two decoupled equations for the transverse V ⊥ i and the longitudinal gauge invariant field the vector field in the bulk is written as V ⊥ i = V 3 = ψ V , where i = 1, 2 denotes the transverse coordinates.

A. The Spectral function
The retarded correlator in momentum space is defined as where a, b are the SU (N f ) indices. The correlator of the vector and transverse axial-vector current are proportional to P µν = η µν − kµkν k 2 , G ab R µν = P µν Π ab (k 2 ). In the thermal states that we consider Lorentz symmetry is broken and only rotational symmetry is left. Then, the projector is split in transverse and longitudinal parts with respect to the spatial momentum k, P µν = P T µν + P L µν , where the transverse part is defined as and P L µν = P µν − P T µν . The Green's function reads In thermal equilibrium, the photon production is given by the light-like correlator (see e.g. [3]), with t a being SU(N f ) flavor matrices. Taking the second derivative with respect to the sources V i we find the transverse part of the retarded two-point function Π T is calculated from the bulk on-shell action and is normalized in order to be the same as as in the flavor singlet case. It is found directly from (4.28) where ψ V is the solution of (3.23) with infalling boundary condition on the horizon and ψ V ( ) = 1 at the boundary. The spectral density is defined as with the spatial indices being contracted. The spectral density can be calculated using the membrane paradigm by redefining a variable proportional to the canonical momentum of The new variable satisfies the following first order equation The incoming regularity condition on the horizon for ψ V translates to ζ (r h ) = 0, hence where the functions with a subindex h correspond to their values at the horizon r = r h . The spectral density in terms of ζ is In the numerical computation of ρ, we use the first order equation, (4.34), since the numerical errors close to the boundary are reduced significantly, [60]. The electric conductivity can be extracted from the time-like limit ω, k → 0 limit of the transverse vector spectral function as Solving (4.34) for k = 0 and ω = 0, with the boundary condition (4.35), we find Equivalently, the electric conductivity can be obtained from the light-like spectral function Therefore, the electric conductivity is proportional to photon-emission rate in low energy.

B. Electric Conductivity and Photon-Emission Rates
In Fig.1, we present the electric conductivity from distinct models in comparison with the latest result from lattice simulations, where the extra normalization factor C em = 2e 2 /3 comes from taking the trace over flavors. From holography, the conductivity rescaled by temperature for the strongly coupled SYM plasma has an analytic expression [11], which is constant in T due to conformal invariance. On the contrary, the most recent simulation for three-flavor QCD [6] shows that the electric conductivity rescaled by T monotonically decreases against T near T c based on the conformal anomaly. The GN model, which saturates the result for the SYM plasma at high T , roughly matches the lattice result around T c but starts to deviate for T > 1.5T c . On the other hand, the VQCD1 and VQCD2 approximately fit the lattice result above T c . For T > 2T c , due to the lack of lattice results, we make different setup in the models such that VQCD2 leads to more rapid growth for σ/T .
In Fig.2, we present the photon-emission rates from different models compared with the ones from pQCD at leading order based on the hard-thermal-loop resummation, which includes the 2 to 2 scatterings [71] and collinear emissions [2], where we take N c = N f = 3 and α EM = e 2 /(4π) = 1/137.
In general, similar to the finding in the SYM plasma [11], the emission rates from holography models have distinct features at low energy compared with the ones in the weakly coupled scenario. However, it is worthwhile to note that the 2 to 2 collisions result in similar shapes of emission rates compared with those from holographic models, whereas the collinear emissions lead to IR divergence. In contrast, the photon production in holographic models incorporate both mechanisms, which "might" suggests that the collinear emissions are suppressed in strong coupling. As shown in the photon production from the SYM plasma with finite t'Hooft coupling in holography [14,72], the peak of the photon emission rate shifts to small k/T when the coupling is reduced.
The finding in holography seems to be consistent with the trend found in the weakly coupled scenario, where the emission rate has blueshift when the coupling is increased [1,2,11]. Nevertheless, unlike the SYM plasma and pQCD with a fixed coupling, non-conformal models now result in the increase of the rescaled rates with temperature, which is qualitatively in accordance with the results for σ/T . As shown in Fig.2, the increase is most prominent at k/T ≈ 2. Although the positions of the peaks in k/T approximately remain unchanged with T for all holographic models and pQCD here, they should in fact shift to the infrared (IR) region in the unit of k at lower temperature.

SIONS
In this section, we explore the phenomenological significance of using different sets of the QGP photon emission rates in modeling the direct photon production in relativistic heavy-ion collisions.

A. Model Setup
In relativistic heavy-ion collisions there are many emission sources to produce photons [73,74]. The dominant contribution to those photons produced with transverse momentum p T < 4 GeV comes from prompt production and thermal radiation. In this work, we only consider these two sources for direct photons in the low p T region. Prompt photons include productions from the initial hard scattering processes (Compton scattering and qq annihilations) as well as from QCD jet fragmentation. The p T -spectra of prompt photons in p-p collisions can be computed using perturbative QCD up to next-to-leading order (NLO) [75,76].
The prompt photon production in nucleus-nucleus collisions is estimated using the number of binary collision N coll scaled photon spectra in pp collisions with cold nuclear and isospin effects included [76,77]. Here, we used the parameterization presented in Ref. [78] for the prompt sources in Pb+Pb collisions at 2.76 A TeV and in Au+Au collisions at 200 A GeV.
To compute thermal photon production from expanding relativistic heavy-ion collisions, we convolute the thermal photon emission rates with event-by-event hydrodynamic medium. We use the state-of-the-art medium evolution as described in Ref. [78,79]. The hydrodynamic simulations are tightly constrained by various hadronic observables. Event-by-event IP-Glasma initial conditions are matched to hydrodynamics at τ sw = 0.4 fm/c and then evolved with lattice QCD based equation of state, s95p-v1 [80,81]. Both shear and bulk viscous effects are included in the medium evolution.
Thermal photon radiation is considered from fluid cells whose temperatures are above 105 MeV. The different sets of photon production rates in the QGP phase, as described in the previous section, are applied to the T > T c (= 180 MeV) temperature region. For the QGP rates derived from a weakly-coupled QCD plasma, 2 to 2 scatterings [71] and collinear emissions [2] are included, which are referred as pQCD rates in previous sections. In the hadronic phase, whose temperature is between 105 and 180 MeV, we use the current most complete set of hadronic photon emission rates, which includes contribution from mesonmeson reactions [82], many-body ρ-spectral function [83,84], π − π bremsstrahlung [85,86], and emission from π − ρ − ω interactions [87].
At the current stage, shear and bulk viscous corrections are only known for a subset of photon production channels. In order to make fair comparisons among the different sets of QGP rates, we will not include any viscous corrections to photon emission rates in this work. Detailed analysis of the viscous effects on direct photon observables were presented in Refs. [31,71,78,88].
The momentum distribution of thermal photons is computed by first producing photons in the local rest of frame of every fluid cell, whose temperature T (x) is higher than the system's freeze-out temperature T freezeout . Then these photons are boosted with the corresponding fluid velocity u(x) to the lab frame, where the thermal photon emission rate is denoted as q dR γ d 3 q . The direct photon spectra is the sum of thermal and prompt photons, Because the thermal photon production rate is suppressed by ∼ α EM /α S compared to hadrons, they are rare in heavy-ion collisions. Thus, the direct photons anisotropic flow coefficients need to be evaluated by correlating with the majority of the soft charged hadrons to ensure sufficient statistics. This is known as the scalar-product method [89]. In the theoretical calculations, the scalar-product direct photon anisotropic flow coefficients, v γ n {SP}, can be computed as [73,90], where v γ n (p T ) and Ψ γ n (p T ) are the magnitude and angle of the n-order harmonic flow for direct photons at a given transverse momentum p T . They are defined as the Fourier coefficients of the photon momentum distribution, supersymmetric partners is expected to give the highest rate and this is turns out to be correct. The GN and VQCD models are non supersymmetric and have the same number of (perturbative) degrees of freedom as QCD.
Therefore, among the different holographic rates, the SYM rates give the most direct photons. Although the electric conductivities in the two VQCD models exceed the one in SYM model for T > 1.5T c , the direct photon yields from the VQCD models are smaller compared to the SYM results. This is because that most of the thermal photon radiations are coming from the phase transition region, 150 < T < 250 MeV, where the space-time volume is the largest. In this temperature region, the photon emission rates are suppressed in the VQCD models compared to the SYM rates. On the other hand, the GN model leads to smaller spectra compared to the ones for VQCD and SYM models as expected from the electric conductivity and emission rates.
In Fig. 4, individual contribution to thermal photons are directly compared among different sets of QGP photon emission rates. Compared to prompt and hadronic radiation  This p T window which is dominated by the QGP photon emission is slightly larger with the holographic rates. The thermal radiation with the SYM rates is about factor of 2 larger than the results using other rates. The VQCD2 rates give the flattest slope of the thermal photon spectra, i.e. largest mean p T . This is because there are more photons radiated from high temperature fluid cells by construction, as for λ large we took w(λ) to be larger than in other models. The thermal radiation from the SYM rates is about factor of 2 larger than the results using other rates.
As the collision energy increases from RHIC to LHC, both the thermal radiation and prompt photon productions increase. The growth of the prompt component is a little bit faster than the thermal radiation. A similar trend was found in the predictions at even higher collision energy [93]. Recently, sizable thermal radiation was found in high multiplicity light-heavy collisions in Ref. [96]. Owing to large pressure gradients, small collision systems, such as p+Pb and d+Au collisions, expand more rapidly compared to the larger Au+Au collisions. This leads to a smaller hadronic phase in these collision systems. Most of the thermal photons come from the hot QGP phase, T > 180 MeV [96]. Hence, the difference between the QGP photon emission rates should be more distinctive in these small collision systems.
In contrast to nucleus-nucleus collisions that were analyzed in the previous sections, full In this work, we will use these calibrated hydrodynamic medium to study the sensitivity of direct photon observables in small systems to the different sets of QGP photon emission rates.
In Figs. 7, direct photon spectra and their anisotropic coefficients are shown for top 0-1% p+Pb collisions at 5.02 TeV and 0-5% d+Au collisions at 200 GeV. In high-multiplicity events of small collision systems, thermal radiation can reach up to a factor of 2 of the prompt contribution. Similar to nucleus-nucleus collisions, the direct photon spectra using the emission rates that are derived from strongly coupled theory are larger than the QCD rates.
The difference is smaller in the d+Au collisions compared to p+Pb collisions at the

E. Discussion
In general, direct-photon production with holographic models is enhanced compared to the one with pQCD especially in intermediate p T , which may stem from the blue shift of thermal-photon-emission rates in strongly coupled scenarios as shown in Section IV, which is also in accordance with Fig.4. On the other hand, we find that the flow of direct photons with holographic models is generally suppressed in low p T . On the contrary, in high p T , the flow from V-QCD and the SYM plasma surpasses the one from pQCD. To analyze the hierarchy of the flow of direct photons by changing sources in the QGP phase, we further illustrate the flow of "thermal photons" generated in the QGP phase, where the photons from hadronic phase are excluded, from distinct sources.
We firstly focus on nucleus-nucleus collisions especially for the LHC collisional energy. As shown in Fig.8, except for the one from pQCD, the hierarchy of flow for "thermal photons" is similar to the one for direct photons in high p T as shown in Fig.6. Since the photon emission at high p T from the QGP phase dominates the one from hadronic phase as illustrated in Fig.4, the momentum anisotropy of the thermal photons from the QGP phase is the only cause of flow. The high p T regime thus allows us to unambiguously compare the effects on flow from different sources without the contamination from hadronic contributions.
It turns out that pQCD, SYM, and VQCD1 lead to larger flow of thermal photons from the QGP phase for both v 2 and v 3 in this region as shown in Fig.8. Since the flow is related to the weight of photon emission at different temperature (equivalently the time elapse) during the expansion of the medium, we may illustrate the ratios of the emission rates at T to the ones at T c in high p T to have an intuitive understanding of the hierarchy of flow. As shown in Fig.9, the high-p T photons are dominantly generated at high temperature 3 . The smaller ratios of emission rates in high temperature should result in larger flow due to larger weights of photons emitted at lower temperature. The hierarchy of emission-rate ratios at high temperature in Fig.9 thus approximately manifests the hierarchy of thermal-photon flow at high p T in Fig.8.
It is worthwhile to note here that the hierarchy of the spectral functions in high energies could be distinct from the hierarchy in low energies. Nevertheless, for the direct-photon flow, one finds that the one for pQCD rate is suppressed at high p T as shown in Fig.6, which contradicts the thermal-photon flow from the QGP phase. The suppression of flow for pQCD stems from the dilution by prompt photons, while the dilution for holographic models is less severe due to larger rates of holographic models at high p T . In conclusion, the direct-photon flow at high p T is governed by the competition between the dilution by prompt photons and the ratios of emission rates of thermal photons in the QGP phase at high temperature.
On the other hand, at low p T or intermediate p T , the thermal photons from hadronic phase become more significant, which play a central role for the direct-photon flow. Due to the interplay between thermal photons created in QGP and hadronic phases, it is more difficult to analyze the hierarchy of flow from distinct models. Nonetheless, the enhancement of the direct-photon spectra for holographic models at low p T is minor. The discrepancies between theoretical predictions and experimental observations in spectra and flow at low p T are more associated with underestimation of thermal photons from hadronic radiation in our study. In the simulations for small collision systems, due to the suppression of thermal photons from the confined phase as shown in (a) and (b) of Fig.7, the flow is primarily engendered by photons from the QGP phase. In addition, one finds that the thermal-photon spectra are more dominant than the prompt-photons spectrum in the 0 − 1% p+Pb collisions, whereas the thermal-photon spectra become more comparable to the prompt-photons spectrum in the 0 − 5% d+Au collisions. The direct-photon flow in the latter case is accordingly suppressed in comparison with the former case as illustrated in (c)-(f) of Fig.7. Moreover, since the spectra from holographic models are more prominent than the one from pQCD in the former case, we find that the flow for holographic models dominates the flow for pQCD at high p T .
However, in the latter case, the flow from different models are more comparable due to their approximate spectra at high p T . In intermediate p T , because of the considerable suppression of prompt photons, the flow reaches maximum in this region. At low p T , photons from hadron gas may start to affect the flow though their influence is highly suppressed compared to the case in AA collisions. On the other hand, since the spectra for different holographic models and pQCD converge, particularly in the 0 − 5% d+Au collisions, the splitting of flow is also reduced. Future comparisons with experimental observations in small collision systems could further verify the importance of the influence from thermal photons in the QGP phase on anisotropic flow.

VI. CONCLUDING REMARKS
In this work, we have evaluated the direct-photon production in heavy ion collisions from holography convoluted with the medium evolution and the inclusion of other sources. The agreement with experiments in spectra is improved compared with the previous study by using the pQCD rate. On the other hand, the deviation in flow is increased at low p T but decreased at high p T .
In small collision systems, where the experimental data of direct photons have not been available, holographic models lead to enhancements in both spectra and flow. Our findings may highlight the strong influence of thermal photons from the QGP phase on the direct-photon flow at high p T , where hadronic contributions are highly suppressed. The enhancement of flow in this region stems from the amplification of the weight of late-time emission and the amplitude of thermal-photon emission in the QGP phase.
In our study, the blue-shift due to the increase of couplings could in general lead to the latter effect, whereas the former one is in fact model dependent.
On the contrary, our study may further suggest that the hadronic contributions are responsible for the flow at low p T . In contrast to the nucleus-nucleus collisions, the dominance of QGP photons could be more pronounced in small collision systems for larger p T window.
Therefore, future measurements of direct-photon spectra and flow in small systems will be crucial to understand the electromagnetic property of the sQGP.