Continuing our development of a mathematical theory of stochastic microlensing, we study the random shear and expected number of random lensed images of different types. In particular, we characterize the first three leading terms in the asymptotic expression of the joint probability density function (p.d.f.) of the random shear tensor at a general point in the lens plane due to point masses in the limit of an infinite number of stars. Up to this order, the p.d.f. depends on the magnitude of the shear tensor, the optical depth, and the mean number of stars through a combination of radial position and the stars' masses. As a consequence, the p.d.f.s of the shear components are seen to converge, in the limit of an infinite number of stars, to shifted Cauchy distributions, which shows that the shear components have heavy tails in that limit. The asymptotic p.d.f. of the shear magnitude in the limit of an infinite number of stars is also presented. Extending to general random distributions of the lenses, we employ the Kac-Rice formula and Morse theory to deduce general formulas for the expected total number of images and the expected number of saddle images. We further generalize these results by considering random sources defined on a countable compact covering of the light source plane. This is done to introduce the notion of {\it global} expected number of positive parity images due to a general lensing map. Applying the result to microlensing, we calculate the asymptotic global expected number of minimum images in the limit of an infinite number of stars, where the stars are uniformly distributed. This global expectation is bounded, while the global expected number of images and the global expected number of saddle images diverge as the order of the number of stars.

Stochastic microlensing is a central tool in probing dark matter on galactic scales. From first principles, we initiate the development of a mathematical theory of stochastic microlensing. Beginning with the random time delay function and associated lensing map, we determine exact expressions for the mean and variance of these transformations. We characterize the exact p.d.f. of a normalized random time delay function at the origin, showing that it is a shifted gamma distribution, which also holds at leading order in the limit of a large number of point masses at a general point of the lens plane. For the large number of point masses limit, we also prove that the asymptotic p.d.f. of the random lensing map under a specified scaling converges to a bivariate normal distribution. We show analytically that the p.d.f. of the random scaled lensing map at leading order depends on the magnitude of the scaled bending angle due purely to point masses as well as demonstrate explicitly how this radial symmetry is broken at the next order. Interestingly, we found at leading order a formula linking the expectation and variance of the normalized random time delay function to the first Betti number of its domain. We also determine an asymptotic p.d.f. for the random bending angle vector and find an integral expression for the probability of a lens plane point being near a fixed point. Lastly, we show explicitly how the results are affected by location in the lens plane. The results of this paper are relevant to the theory of random fields and provide a platform for further generalizations as well as analytical limits for checking astrophysical studies of stochastic microlensing.