Three-body decay functions in space-like parton branches are implemented to evaluate transverse-momentum-dependent (TMD) parton distribution functions in the next-to-leading logarithmic (NLL) order of quantum chromodynamics (QCD). Interference contributions due to the next-to-leading-order terms are taken into account for the evaluation of the transverse momenta in initial state parton radiations. Some properties of the decay functions are also
examined. As an example, the calculated results are compared with those evaluated by an algorithm proposed in [M. A. Kimber, A. D. Martin, and M. G. Ryskin, Eur. Phys. J. C 12, 655 (2000)], [M. A. Kimber, A. D. Martin, andM.G.Ryskin, Phys.Rev.D63, 11402 (2001)], [G.Watt, A. D. Martin, and M. G. Ryskin, Eur. Phys. J. C 31, 73 (2003)], and [A. D. Martin, M. G. Ryskin, and G.Watt, Eur. Phys. J. C 66, 167 (2010)], in which the TMD parton distributions are defined based
on the kt -factorization method with angular ordering conditions due to interference effects.