Charged particles in a magnetosphere are spontaneously attracted to a planet while increasing their kinetic energy via inward diffusion process. A constraint on particles' micro-scale adiabatic invariants restricts the class of motions available to the system, giving rise to a proper frame on which particle diffusion occurs. We investigate the inward diffusion process by numerical simulation of particles on constrained phase space. The results reveal the emergence of inhomogeneous density gradient and anisotropic heating, which is consistent with spacecraft observations, experimental observations, and the recently formulated diffusion model on the constrained phase space.