Probing thermal nature of matter formed at RHIC via fluctuations

Arguments for thermalization of the QCD matter created in high-energy nuclear collisions has dominantly come from the agreement of the measured yields of produced hadrons with those from statistical thermal models. Ideally for a thermalized system, in addition to mean, the higher orders of the moments of the multiplicity distribution of produced particles should also show agreement with thermal models. In this respect, simultaneously studying the moments of the event-by-event distributions of conserved quantities like net-baryon, net-strangeness and net-charge number is best suited. We present a systematic study of comparing the results from a thermal hadron resonance gas (HRG) model with data on higher moments of net-proton, net-kaon and net-charge distributions measured at RHIC beam energy scan program. The experimental acceptances in terms of rapidity and transverse momentum are used in the model calculations which also include resonance decay. For the first time, the HRG model results are found to explain the measurements up to third order of moment with a common temperature and baryonic chemical potential. These calculations have tested the thermal nature of produced net-particle distributions up to third order, thereby providing evidence for thermalization of the QCD matter formed in such high energy heavy-ion collisions.