The transverse momentum spectrum of η meson in relativistic heavy-ion collisions is studied at the Next-to-Leading Order (NLO) within the perturbative QCD, where the jet quenching effect in the QGP is incorporated with the effectively medium-modified η fragmentation functions using the higher-twist approach. We show that the theoretical simulations could give nice descriptions of PHENIX data on η meson in both p+p and central Au+Au collisions at the RHIC, and also provide numerical predictions of η spectra in central Pb + Pb collisions with √ sNN = 2.76 TeV at the LHC.The ratios of η/π 0 in p + p and in central Au + Au collisions at 200 GeV are found to overlap in a wide pT region, which matches well the measured ratio η/π 0 by PHENIX. We demonstrate that, at the asymptotic region when pT → ∞ the ratios of η/π 0 in both Au + Au and p + p are almost determined only by quark jets fragmentation and thus approach to the one in e + e − scattering; in addition, the almost identical gluon (quark) contribution fractions to η and to π result in a rather moderate variation of η/π 0 distribution at intermediate and high pT region in A + A relative to that in p+p; while a slightly higher η/π 0 at small pT in Au+Au can be observed due to larger suppression of gluon contribution fraction to π 0 as compared to the one to η. The theoretical prediction for η/π 0 at the LHC has also been presented. PACS numbers: 12.38.Mh; 25.75.-q; 13.85.NiThe strong suppression of single hadron production at large transverse momentum [1,2] has provided the convincing evidence of the jet quenching phenomena discovered in relativistic heavy-ion collisions (HIC) [3]. Extensive phenomenological investigations [4-9] and experimental measurements [10-16] on the suppression of single hadron spectra at high p T have been carried out at both the RHIC and the LHC. As the first observable of jet quenching phenomena, the yield suppression of inclusive hadrons is arguably the most thoroughly studied quantity of jet quenching, and provides an indispensible tool to extract the properties of the hot medium created in nucleus-nucleus collisions by comparing theoretical calculations with experimental measurements, such as the jet transport coefficientq [17]. The interplay between theory and experiment on the single hadron production will help constraining the longitudinal distribution of parton energy loss in hot/dense QCD medium, and better understanding the jet-medium interactions after being combined with studies of full jets which also shed light on the angular distribution of the medium-induced gluon radiation and thus constrain the transverse distribution of parton energy loss as well [18][19][20][21][22][23].So far, most of the theoretical calculations on single hadron productions in HIC focus on π meson or charged hadrons (where π also giving a predominant contribu- * tion), and there are very few studies on other identified hadrons [24][25][26][27]. We note that η meson is the second important source of decay electrons and photons just after the π 0 , and...