A new set of functions, which form a basis of the massive nonpower expansion for physical observables, is presented in the framework of the analytic approach to QCD at the four-loop level. The effects due to the π meson mass are taken into account by employing the dispersion relation for the Adler function. The nonvanishing pion mass substantially modifies the functional expansion at low energies. Specifically, the spacelike functions are affected by the mass of the π meson in the infrared domain below few GeV, whereas the timelike functions acquire characteristic plateaulike behavior below the two-pion threshold. At the same time, all the appealing features of the massless nonpower expansion persist in the considered case of the nonvanishing pion mass.The renormalization group (RG) method plays a key role in the framework of the Quantum Field Theory (QFT) and its applications. Indeed, one is able to handle reliably the strong interaction processes at high energies by employing this method together with perturbative calculations. However, such perturbative solutions to the RG equation possess unphysical singularities in the infrared domain, a fact that contradicts the general principles of the local QFT, and significantly complicates the theoretical description and interpretation of the intermediate-and low-energy experimental data. Nevertheless, an effective way to overcome these difficulties is to complement the perturbative results with a proper nonperturbative insight into the infrared hadron dynamics.One of the sources of the nonperturbative information about the strong interaction processes is the dispersion relations. The idea of employing the latter together with perturbation theory forms the underlying concept of the so-called analytic approach to QFT, which was first proposed in the framework of Quantum Electrodynamics [1]. Recently, this approach has been extended to Quantum Chromodynamics (QCD) [2] and applied to the "analytization" of the perturbative power series for the QCD observables [3,4,5]. The term analytization means the restoring of the correct analytic properties in the kinematic variable of a quantity under consideration by making use of the Källén-Lehmann integral representation (positive q 2 corresponds to a spacelike momentum transfer hereinafter)with the spectral function defined by the initial (perturbative) expression for the quantity at hand:(2) However, there are several ways to embody the analyticity requirement into the RG formalism, that eventually has given rise to different models for the analytic running coupling.Thus, in the original model due to Shirkov and Solovtsov [2] the analyticity requirement (1) is imposed on the perturbative running coupling itself. At the one-loop level this leads towhereas at the higher loop levels the integral representation of the Källén-Lehmann type α(q 2 ) = 4πholds for this invariant charge. Ultimately, the prescription [2] results in the infrared finite limiting value for the running coupling (see papers [3,4,5] and references therein for t...