Models of differentially rotating protoneutron stars are calculated, using realistic equations of state of dense hot matter. Various conditions within the stellar interior, corresponding to different stages of protoneutron star evolution, are considered. Families of differentially rotating models of a given baryon mass are calculated, using a two-parameter formula describing the angular velocity profile within a rotating protoneutron star. Maximum angular momentum, which can be accommodated by a protoneutron star is calculated, for various thermal conditions in stellar interior, for a baryon mass of $1.5 M_\odot$. In the case of a thermally homogeneous (isentropic or isothermal) neutrino-opaque interior this maximum angular momentum turns out to be somewhat higher than that of a cold neutron star of the same baryon mass, rotating uniformly at the mass shedding angular velocity. However, if the protoneutron star has a thermal structure characteristic of initial state, with a low entropy (unshocked) core, and a high entropy (shocked) outer half of baryon mass, the maximum angular momentum is significantly lower. This leads to a minimum period of uniform rotation of cold neutron stars of baryon mass $\sim 1.5~M_\odot$, formed directly (i.e. without a subsequent significant accretion of mass) from protoneutron stars with shocked envelope, of about 1.7 ms and strengthens the hypothesis that millisecond pulsars are accretion accelerated neutron stars.

Rapid uniform rotation of newborn neutron stars (protoneutron stars) is studied for a range of internal temperatures and entropies per baryon predicted by the existing numerical simulations. Calculations are performed using general relativistic equations of hydrostatic equilibrium of rotating, axially symmetric stars. Stability of rotating configurations with respect to mass shedding and the axially symmetric perturbations is studied. Numerical calculations are performed for a realistic dense matter equation of state, under various assumptions concerning neutron star interior (large trapped lepton number, no trapped lepton number, isentropic, isothermal). For configurations with baryon mass well below the maximum one for the non-rotating models, the mass shedding limit depends quite sensitively on the position of the ``neutrinosphere''. The absolute upper limit on rotation frequency is, to a good approximation, obtained for the maximum baryon mass of rotating configurations. Empirical formula for the maximum rotation frequency of uniformly rotating protoneutron stars is shown to be quite precise; it actually coincides with that used for cold neutron stars. Evolutionary sequences at fixed baryon mass and angular momentum, which correspond to evolution of protoneutron stars into cold neutron stars are studied, and resulting constraints on the maximum rotation frequency of solitary pulsars are discussed.