We discuss how the non-equilibrium process of pion production within the Zubarev approach of the non-equilibrium statistical operator leads to a theoretical foundation for the appearance of a non-equilibrium pion chemical potential for the pion distribution function for which there is experimental evidence in experiments at the CERN LHC.