NLO Productions of ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document} and KS0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^0_{\mathrm{S}}$$\end{document} with a global extraction of the jet transport parameter in heavy-ion collisions

In this work, we pave the way to calculate the productions of ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document} and KS0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^0_{\mathrm{S}}$$\end{document} mesons with large pT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_T$$\end{document} in p+p and A+A collisions both at RHIC and LHC. The fragmentation functions (FFs) of the ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document} meson in vacuum at next-to-leading order (NLO) are obtained by evolving the NLO DGLAP evolution equations with rescaled ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document} FFs at initial scale Q02=1.5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q_0^2=1.5$$\end{document} GeV2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {GeV}^2$$\end{document} from a broken SU(3) model, and the FFs of KS0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^0_{\mathrm{S}}$$\end{document} in vacuum are taken from AKK08 parametrization directly. Within the framework of the NLO pQCD improved parton model, we arrive at good descriptions of the experimental data on ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document} and KS0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^0_{\mathrm{S}}$$\end{document} in p+p both at RHIC and LHC. With the higher-twist approach, to take into account jet quenching effect by medium-modified FFs, nuclear modification factors for ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document} meson and KS0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^0_{\mathrm{S}}$$\end{document} meson both at RHIC and LHC are presented with different sets of jet transport coefficients q^0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{q}}_0$$\end{document}. Then we make a global extraction of q^0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{q}}_0$$\end{document} both at RHIC and LHC by confronting our model calculations with all available data on six identified mesons: π0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi ^0$$\end{document}, η\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document}, ρ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho ^0$$\end{document}, ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document}, ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document}, and KS0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^0_{\mathrm{S}}$$\end{document}. The minimum value of total χ2/d.o.f\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi ^2/d.o.f$$\end{document} for productions of these mesons gives the best value of q^0=0.5GeV2/fm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{q}}_0=0.5\mathrm ~GeV^2/fm$$\end{document} for Au+Au collisions with sNN=200\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_{\mathrm{NN}}}=200$$\end{document} GeV at RHIC, and q^0=1.2GeV2/fm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{q}}_0=1.2\mathrm ~GeV^2/fm$$\end{document} for Pb+Pb collisions with sNN=2.76\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_{\mathrm{NN}}}=2.76$$\end{document} TeV at LHC, respectively, with the QGP spacetime evolution given by the event-by-event viscous hydrodynamics model IEBE-VISHNU. With these global extracted values of q^0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{q}}_0$$\end{document}, nuclear modification factors of π0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi ^0$$\end{document}, η\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document}, ρ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho ^0$$\end{document}, ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document}, ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document}, and KS0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^0_{\mathrm{S}}$$\end{document} in A+A collisions are presented, and predictions of yield ratios such as ω/π0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega /\pi ^0$$\end{document} and KS0/π0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^0_{\mathrm{S}}/\pi ^0$$\end{document} at the high-pT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_T$$\end{document} regime in heavy-ion collisions both at RHIC and LHC are provided.


Introduction
The jet quenching effect describes energy dissipation of an energetic parton when it traverses through the hot and dense QCD medium, which is produced shortly after high energy nuclear collisions [1]. The single hadron production suppression at the highp T regime when compared to scaled p+p data is one of the primary perturbative probes to study the properties of this de-coupled quark and gluon QCD matter [2]. π 0 is the best measured final-state hadron, and its nuclear modification factor R AA as a function of transverse momentum p T is interpreted as the consequence of the jet quenching effect and it did help us to constrain the strength of the jet-medium interaction and the properties of the QCD medium [3][4][5][6][7][8][9]. In 2013, the Jet Collaboration summarized different jet quenching theoretical frameworks and compared the jet transport parameter extracted by different energy loss models using the same hydro description of QCD medium [10]. Meanwhile, the R AA for different final-state identified hadrons have been measured, and their production suppressions as well as patterns of their yield ratios have been observed [4,5,[11][12][13][14][15][16]. It is of interest and it is a challenge to describe the cross sections of leading hadrons of different types and their yield ratios with each other in heavy-ion collisions (HIC) both at RHIC and LHC with a unified model of jet quenching, which should shed light on the flavor dependence of the parton energy loss and the intrinsic properties of the identified hadron productions in p+p and A+A reactions [17][18][19][20][21][22][23][24].
We have achieved a better understanding of the suppression patterns of different mesons by conducting calculations and by analysis of the R AA and particle ratios for other identified hadron productions [25][26][27]. We find that understanding of the suppression pattern of the leading hadron requires one to take into account all three factors [25]: the initial hard jet spectrum, the energy loss mechanism and the parton fragmentation functions in vacuum. The flavor dependence of the energy loss will result in a decrease of the fraction of gluon fragmenting contribution in the nuclear-nuclear collision and an increase of the fraction of the quark fragmenting contribution. The productions of final-state mesons, like π 0 , ρ 0 , η [25][26][27], are dominated by the quark fragmentation contribution at the highp T regime in p+p collisions. The energy loss effect will only enhance the domination of the quark fragmenting contribution in A+A collisions. Therefore ρ 0 /π 0 and η/π 0 in p+p and A+A collisions coincide at largep T , since they are only determined by the ratios of the FFs in vacuum. But the production of φ meson is dominated by the gluon fragmenting contribution in p+p collision, therefore we can observe a separation of φ/π 0 in p+p and A+A. It is of great interest to gain better understanding of the suppression pattern of different final-state hadrons, and furthermore to provide systematical pQCD predictions for the current experimental measurement of the identified hadrons. In this letter, we mainly investigate one of the low-mass vector mesons (ω) and one of the pseudoscalar mesons (K 0 S ), together with π 0 , η, ρ 0 and φ to achieve these goals.
In this article, we investigate the production of two other mesons, ω and K 0 S , with large p T in A+A collisions, which has never been computed so far, to the best of our knowledge. The ω meson is constituted of a similar valence quark of π 0 with larger mass 782.65 MeV and spin 1. Kaons are a group of lightest mesons, carrying strangeness components, K 0 S ( ds+sd 2 ) is one type of kaon, consisting of s-quark, d-quark and their corresponding anti-quarks. The productions of these two mesons have been measured in p+p and A+A collisions both at RHIC and LHC, but a theoretical description is lacking. Within the NLO pQCD improved parton model, we calculate ω and K 0 S yields at the highp T regime in heavy-ion collisions, by employing medium-modification fragmentation functions (FFs) due to gluon radiation in the hot/dense QCD medium in the higher-twist approach of jet quenching [22,23,[28][29][30], the same approach as in our calculations on the productions of η, ρ 0 and φ mesons in HIC [25][26][27].
It is noted in the previous studies that, for consistency, we implemented the Hirano hydro description [31,32] to describe the spacetime evolution of the QGP fireball. In this work, we utilize a state-of-the-art, event-by-event (2+1)-D viscous hydrodynamics model (IEBE-VISHNU) [33] to give the spacetime evolution information of the hot and dense medium. Due to the changes of the medium description, the strength of the jet-medium interaction characterized by the jet transport parameterq 0 should be re-extracted. Taking advantage of the systematical study of the nuclear modification factors with respect to p T and the large amount of experimental data of π 0 , η, ρ 0 and φ both at RHIC and LHC, with two more mesons ω and K 0 S calculated in this article, we can make a global extraction of the jet transport parameterq 0 with all the available experimental data of these six identified mesons R AA .
The article is organized as follows. We first present the theoretical framework of computing single hadron cross sections in p+p collision in Sect. 2, and we give the p+p baseline in investigating the in-medium modification of these productions. In Sect. 3, we discuss the inclusive hadron production in A+A collisions with the medium-modified FFs, based on the higher-twist approach of parton energy loss, and we present the nuclear modification factors R AA for ω and K 0 S both at RHIC and LHC. Section 4 shows the global extraction of the jet transport coefficientq 0 both at RHIC and LHC by confronting our model calculations with all available data on six identified mesons: π 0 , η, ρ 0 , φ, ω, and K 0 S . In Sect. 5, we make predictions on the identified hadron yield ratios ω/π 0 and K 0 S /π 0 in p+p and A+A collisions, and we compare them with experimental data if applicable. We give a brief summary in Sect. 6.

Large p T yield of ω and K 0 S meson in p+p
We start with the productions of ω and K 0 S mesons at NLO in p+p collisions. In the framework of the pQCD improved parton model, the inclusive hadron production in p+p collisions can be given by the convolution of three parts, the parton distribution functions (PDFs) in a proton, the hard partonic scattering cross section, denoted dσ/dt (up to the order of α 3 s ), and the parton fragmentation functions (FFs) to the final-state hadron D q(g)→h (z h , Q 2 ). One may convolute PDFs and partonic cross section dσ/dt into the initial hard (parton-) jet spectrum F q,g ( p T ), and we have Here, the quark and gluon fragmenting contributions are written separately to facilitate future discussions. A next-toleading order (NLO) Monte Carlo code has been employed to calculate leading hadron productions in the p+p collision [36], and the CT14 parametrization of PDFs for the free proton [37] has been implemented. In the framework, as long as the parametrization of the specific final-state hadron FFs in vacuum is available, one can predict its production yield at large p T in elementary p+p collisions in principle.
In order to make the NLO calculation of ω and K 0 S mesons in p+p, the NLO parton FFs of these two mesons are needed. The parton FFs of the K 0 S meson at NLO can be found in the AKK08 parametrization [38], while there is not such a kind of global parametrization for ω FFs in vacuum. So we have to rely on theoretical models and utilize ω parton FFs at a starting scale Q 2 0 = 1.5 GeV 2 provided by a broken SU(3) model [39,40]. We note that our preceding investigations on leading ρ 0 and φ production [26,27] have benefited from this broken SU(3) model. In this model, the independent parton FFs of different flavor are reduced into three independent functions named valence function, V (x, Q 2 0 ), sea function, γ (x, Q 2 0 ) and gluon function, D g (x, Q 2 0 ), by considering the SU(3) flavor symmetry with a symmetry breaking parameter and also addressing the isospin and charge conjugation invariance of the vector mesons. We can write these three functions into a standard polynomial at the starting low energy scale of Q 2 0 = 1.5 GeV 2 as where the full set of parameters (a, b, c, d, e) defined in V , γ , D g and a few additional parameters defined for each vector meson such as the strangeness suppression factor λ, the vector mixing angle θ , the sea suppression factor f ω sea , f u 1 (ω) and f ω g have been determined in the broken SU(3) model in Refs. [39,40], and we multiply the parameter a i by 3 to make the best fit to the yield of the ω meson in p+p collisions. All these parameters for fitting the initial FFs of ω meson are listed in Table 1. To obtain a NLO parton FF for ω at any energy scale Q, we employ the numerical NLO DGLAP evolution program provided in Ref. [41] with the initial parton FFs starting scale Q 2 0 = 1.5 GeV 2 as input. We demonstrate in Fig. 1 the NLO DGLAP evolved FFs of the ω meson as functions of z h at fixed Q 2 = 1.5 GeV 2 and Q 2 = 100 GeV 2 on the left and as functions of Q at fixed z h = 0.4 and z h = 0.6 on the right. We find in the typical fraction region ( [27]. Due to the fact that in parton FFs the gluon contribution exceeds the quark contribution, an overwhelming advantage of the gluon fragmenting contribution of the final-state ω production in p+p collision is expected in the competition with the quark fragmenting contribution. In Fig. 2, we plot the NLO FFs of K 0 S from the AKK08 parametrization in the same manner as the ω meson in Fig. 1. We find the strange quark FF therefore competition of the gluon and quark fragmenting contributions is expected in the K 0 S production in p+p collisions based on the pattern discovered in the previous study of the π 0 , η, ρ 0 and φ mesons.
With the NLO FFs in vacuum of both the ω and the K 0 S mesons, we are able to confront the numerical results of the ω and K 0 S meson productions in p+p collisions up to NLO with the available experimental data both at RHIC and LHC shown in Figs. 3 and 4. Note that, in our simulations, the hard scales such as factorization scale, renormalization scale and fragmentation scale are chosen to be the same and proportional to the p T of the final-state hadron. We see that the numerical results match well with the experimental data on ω meson spectra both at RHIC, √ s NN = 200 GeV, and LHC, √ s NN = 7 TeV, and the hard scales, μ = 0.5 p T , as shown in Fig. 3. We also demonstrate in Fig. 4 the good agreement of the theoretical calculations of the leading K 0 S meson with the STAR data and the ALICE data, where the hard scales are fixed to be μ = 1.0 p T to give the best fit.

Large p T yield of ω and K 0 S meson in A+A
To calculate the ω and K 0 S productions in A+A collisions, the jet quenching effect should be included. In this article we utilize the higher-twist approach of parton energy loss, which relates the parton energy loss due to its multiple scattering in QCD medium and medium-induced gluon radiation to twist-four processes, and this shows that these processes give rise to additional terms in QCD evolution equations and lead to the effectively medium-modified fragmentation functions; thus the partonic energy loss effect can be taken into account by replacing the FFs in vacuum with the effec-  Fig. 4 The theoretical results of the K 0 s meson production in p+p collisions at √ s = 200 GeV RHIC compared with STAR data [35] and at √ s = 2.76 TeV LHC compared with ALICE data [16] tively medium-modified FFs [28][29][30]42,43]. By assuming a thermal ensemble of quasi-particle states in the hot/dense medium and neglecting the multiple particle correlations inside the medium, the effective medium-modified FFs are given as [22,23,[25][26][27]: The contribution of the medium-induced gluon radiation is attributed to the medium-modified splitting functions represented by Δγ q→qg and Δγ q→gq . In this formula, we assume that the energy loss of the jet propagating through the medium is totally carried away by the radiative gluons, then the convolution of these energy loss kernels Δγ q→qg (Δγ q→gq ) with the (DGLAP) evolved vacuum FFs at any scale D h q,g ( z h z , Q 2 ) implies the assumption that the fast parton first loses its energy in the medium and then fragments into final-state hadrons in vacuum.
The medium-modified splitting functions depend on the properties of the local medium, which cannot be determined directly by the theoretical calculation itself. The jet transport parameterq, which is defined as the average squared transverse momentum broadening per unit length, is therefore introduced to profile the dependency of the local medium properties in these energy loss kernels [10]. The jet transport parameterq can be phenomenologically assumed to be proportional to the local parton density in the hot/dense medium to adopt a further medium description: where ρ med is the parton (quark and gluon) density in the medium at a given temperature, q 0 is the initial jet transport parameter at the center of the bulk medium at initial time τ 0 , p μ is the four-momentum of the jet and u μ is the four-flow velocity of the probed local medium in the collision frame. In this article we employ the state-of-art, event-by-event viscous hydrodynamics description IEBE-VISHNU [33] to give the spacetime evolution information of the hot/dense medium, such as temperature, energy density, four-velocity of the local medium at any evolution time and position, and the formation time of QGP, τ 0 = 0.6 fm. Therefore the only undetermined parameter in the calculation isq 0 , representing the strength of the jet-medium interaction, which can be fixed by fitting the data of leading hadrons (usually π meson) suppression in A+A collisions [10,[25][26][27].
When averaging the medium-modified FFs over the initial production position and the jet propagation direction, we are able to directly replace the vacuum FFs in the p+p formalism with these averaged medium-modified FFs ; therefore the formalism of the calculated cross section of the single hadron productions in HIC would take the form of Here, N AB coll (b) = d 2 rt A (r )t B (|b − r|) is the number of binary nucleon-nucleon collisions at a certain impact parameter b in A+B collisions; it can be calculated by using the Glauber model [44]. f a/A (x a , μ 2 ) represents the effective PDFs inside a nucleus. In our calculations, we employed EPPS16 NLO nuclear PDFs [45] to include the initial-state cold nuclear matter effects on single hadron productions [46].
The nuclear modification factor R AA as a function of p T is calculated as cross sections in A + A collisions divided by the ones in p+p collision, scaled by the averaged number of binary nucleon-nucleon collisions with the chosen impact parameter b [3]: The theoretical results of R AA at various values ofq 0 = 0.4 − 0.7 GeV 2 /fm for both ω and K s mesons are presented in Figs. 5 and 6.
When obtaining a suitable value ofq 0 by comparing the theoretical calculations of R AA for ω and K s with the corresponding data, there are two caveats. First, since the data of ω and K 0 S meson at large p T in A+A collisions are rather limited and come with a large uncertainty, it is difficult to make a good constraint onq 0 with these data. Second, there is the caveat that we employ the IEBE-VISHNU hydrodynamics model to describe the spacetime evolution of the fireball, which gives different information of physics quantities, such as temperature and density, from those provided by other hydro models, such as the Hirano hydro description [31,32]. Therefore, we could not take advantage of the extracted value ofq 0 in Refs. [10,22,23,[25][26][27], where the Hirano hydro description has been utilized.   [4,5], ρ 0 [5], K 0 S [5], η [11], φ [12], and ω [15] With these two cautions in mind, and realizing that with our model we are now ready to make a systematic study of six types of identified mesons such as π 0 , η, ρ 0 , φ, ω, and K 0 S in heavy-ion collisions, it will be of great interest to make a global extraction of the jet transport coefficientq 0 both at RHIC and LHC by confronting our model calculations (with the IEBE-VISHNU hydro model) against all available data on six identified mesons: π 0 , η, ρ 0 , φ, ω, and K 0 S , and then make precise calculations on the nuclear modification factors of these mesons including ω, and K 0 S , and the yield ratios of these six mesons in HIC.

Global extraction ofq 0 with R AA for six identified mesons
We perform a systematic calculation of the R AA of 6 identified mesons (π 0 , ρ 0 , η, φ, ω, K 0 S ) to compare with all their available experimental data both at RHIC and LHC in Figs. 7 and 8.
In order to extract the best value of the initial jet transport parameterq 0 , we perform a χ 2 fit to compare the theoretical results with differentq 0 and the available experimental data of all the different final hadrons. We have In the above equations, D i represents the experimental grids and T i is our theoretical prediction at input parameter a. σ 2 i means the ith systematic and statistical experimental errors. We show in the top panel of Fig. 9 the derived χ 2 averaged by the number of the compared data points for dif-  Fig. 9 we plot the curve χ 2 /d.o.f as a function ofq 0 both at RHIC and LHC, where the minimum of the curve χ 2 /d.o.f with respect toq 0 presents the best fit of theory with data. We then observe that at RHIC the minimum point of χ 2 /d.o.f of all six identified mesons gives the best value ofq 0 = 0.5(+ 0.15/− 0.05) GeV 2 /fm. It is noted that the production of π 0 will carry the largest weight due to its more abundant data with relatively smaller error bars. We also derive the best value ofq 0 in the Pb+Pb collisions at LHC √ s NN = 2.76 TeV to beq 0 = 1.2(+ 0.25/ − 0.15) GeV 2 /fm, though the curve χ 2 /d.o.f at LHC is much flatter that at the RHIC. So theoretical results withq 0 = 1.1 − 1.4 GeV 2 /fm should all give decent descriptions as regards the data at LHC.
It is noted that in our current model the extracted values of jet transport coefficientq 0 both at RHIC and LHC are smaller than that by the JET Collaboration in Ref. [10] and the ones in our preceding calculations [25][26][27]. These differences come mainly from the different hydro models utilized in the studies. We have checked [47] that if we employ the Hirano hydro description in the calculation, the extracted values of q 0 from our global fitting will be consistent with those in Refs. [10,[25][26][27], though the model in this article has the potential to give more precise extraction of the jet transport coefficients when more data of identified hadrons in A+A collisions become available in the near future.
5 Particle ratios of ω/π 0 and K 0 s /π 0 in A+A With the global extracted value ofq 0 discussed in Sect. 4, we are able to further investigate the particle ratio of ω and K 0 S both at RHIC and LHC. We first calculate ω/π 0 ratio as a function of p T and show the results in p+p and Au+Au at RHIC in the left panel of Fig. 10, where the PHENIX experimental data on the ω/π 0 ratio in p+p are also illustrated. An enhancement of the ratio in A+A relative to that in p+p is found in the smallp T region, whereas we have a small suppression in the highp T regime. We also predict the ω/π 0 ratio as a function of p T in p+p and Pb+Pb collisions with √ s N N = 2.76 TeV at LHC. In Fig. 10 we do not see the overlapping of the curves ω/π 0 in p+p and A+A at the highp T regime, as we have observed for the ratios η/π 0 [25] and ρ 0 /π 0 [26] in p+p and A+A.
To understand this feature deeper, we plot the gluon and quark (fragmentation) contribution fractions to ω (and π 0 ) in p+p and Au+Au at RHIC in Fig. 11. One sees in the p+p collision at RHIC that the production of ω is dominated by gluon fragmentation. In the A+A collisions, the gluon suffers more energy loss than the light quark due to its larger color factor, which decreases the gluon contribution fraction and increases the light quark contribution fraction, but the On the other hand, π 0 meson production is light quark fragmentation dominant in p+p, and the jet quenching effect further enhances this dominance of light quark fragmentation to π 0 in A+A. Therefore, we see that the yield ratio of ω/π 0 in A+A should be suppressed relative to that in p+p due to the energy loss effect, and it separates from the one in p+p even at very high p T . We also compute the K 0 S /π 0 ratio as a function of p T both at RHIC and LHC in Fig. 12. We find the curves in A+A and in p+p are approaching to each other with p T increasing, and an obvious coincidence of these two curves is seen at LHC. We show the gluon and quark contribution fractions to the K 0 S meson (and the π 0 meson) as functions of p T in Fig. 13. We find in p+p collisions that the productions of both K 0 S and π 0 at very large transverse momenta are dominated by quark fragmentation. In A+A collisions, the gluon contribution shall be further suppressed because the gluon generally loses more energy. Thus, both in p+p and A+A collisions, the ratio K 0 S /π 0 should be largely determined by the ratio of quark FFs for K 0 S (D K 0 S q (z h , Q 2 )) to quark FFs for π 0 (D π 0 q (z h , Q 2 )) at very high p T , where these FFs vary slowly with the momentum fraction z h . This is very similar to the case of η 0 /π 0 at high p T [25]. Even though in A+A collisions the jet quenching effect can shift z h of quark FFs, if the quark FFs have a rather weak dependence on z h and p T , we can see that at the very highp T regime the curves for K 0 S /π 0 in A+A and p+p come close to each other, and they even coincide at LHC.

Summary
In summary, we obtain the NLO FFs of the ω meson in vacuum by evolving the rescaled ω FFs from a broken SU(3) model at a starting scale Q 2 0 = 1.5 GeV 2 , and we directly employ NLO K 0 S FFs in vacuum from the AKK08 parameterizations; the numerical simulation of productions of both ω and K 0 S matches well with the experimental data in p+p reactions. With the IEBE-VISHNU hydro profile of the QCD medium, we calculate the nuclear modification factors of the ω and K 0 S mesons as well as π 0 , η, φ, ρ 0 in A+A collisions both at RHIC and LHC, including the jet quenching effect in a higher-twist approach. The global extraction of the jet transport parameterq 0 is made, with comparison of the theoretical calculation and the experimental data of all six identified mesons: π 0 , ρ 0 , η, φ, ω, K 0 S . Furthermore, we predict the yield ratios of ω/π 0 both at RHIC and LHC, and a fairly good agreement of the theoretical results and experimental data is found at RHIC. Theoretical predictions of K 0 S /π 0 ratios as functions of p T at RHIC and LHC are also presented.

A: Input parameters for ω in broken SU(3) model
See Table 1   Table 1 The parameters for the input fragmentation functions of ω at the initial scale Q 2 = 1.5 GeV 2 based on broken SU(3) model in our calculation ω V a 2.16