The Kwiecinski equations for the QCD evolution of the unintegrated parton
distributions in the transverse-coordinate space (b) are analyzed with the help
of the Mellin-transform method. The equations are solved numerically in the
general case, as well as in a small-b expansion which converges fast for b
Lambda_QCD sufficiently small. We also discuss the asymptotic limit of large bQ
and show that the distributions generated by the evolution decrease with b
according to a power law. Numerical results are presented for the pion
distributions with a simple valence-like initial condition at the low scale,
following from chiral large-N_c quark models. We use two models: the Spectral
Quark Model and the Nambu--Jona-Lasinio model. Formal aspects of the equations,
such as the analytic form of the b-dependent anomalous dimensions, their
analytic structure, as well as the limits of unintegrated parton densities at x
-> 0, x -> 1, and at large b, are discussed in detail. The effect of spreading
of the transverse momentum with the increasing scale is confirmed, with
growing asymptotically as Q^2 alpha(Q^2). Approximate formulas for
for each parton species is given, which may be used in practical
applications.Comment: 18 pages, 6 figures, RevTe