  A special purpose N-body simulation has been built to understand the tidal
heating of the smallest dark matter substructures ($10^{-6}\msun$ and 0.01pc)
from the grainy potential of the Milky Way due to individual stars in the disk
and the bulge. To test the method we first run simulations of single encounters
of microhalos with an isolated star, and compare with analytical predictions of
the dark particle bound fraction as a function of impact parameter. We then
follow the orbits of a set of microhalos in a realistic flattened Milky Way
potential. We concentrate on (detectable) microhalos passing near the Sun with
a range of pericenter and apocenter. Stellar perturbers near the orbital path
of a microhalo would exert stochachstic impulses, which we apply in a Monte
Carlo fashion according to the Besancon model for the distribution of stars of
different masses and ages in our Galaxy. Also incorporated are the usual
pericenter tidal heating and disk-shocking heating. We give a detailed
diagnosis of typical microhalos and find microhalos with internal tangential
anisotropy are slightly more robust than the ones with radial anisotropy. In
addition, the dark particles generally go through of a random walk in velocity
space and diffuse out of the microhalos.
  We show that the typical destruction time scales are strongly correlated with
the stellar density averaged along a microhalo's orbit over the age of the
stellar disk. We also present the morphology of a microhalo at several epochs
which may hold the key to dark matter detections.
