In order to mimic the Brownian particle in liquid, molecular dynamics calculations of dilute solutions of spherical fullerene molecules with various sizes in liquid Ar were performed. To establish the scaling equation for self-diffusion coefficient, D, of the fullerenes, the dependence of D was examined on the mass ratio of solute to solvent and on the energy-parameter ratio used in the Lennard-Jones potentials. The dependence on the energy-parameter ratio remains up to C_{540}, whereas D rapidly becomes independent of the mass ratio with increasing mass of the solute. The product of the scaling equations obtained for the D of the solute and for shear viscosity, η_{sv}, for the solvent gives a relation which replaces the Stokes-Einstein relation based on the hydrodynamics. The present expression does not need both the boundary conditions and the hydrodynamic particle size, but instead the energy-parameter ratio, packing fraction of solvent, and bare size of solute. From the viewpoint of the tackiness at the boundary, the cage correlation function around the diffusing particle was examined; it was found that the decay time of the function depends mainly on the the energy-parameter ratio. Therefore, the energy-parameter ratio accounts for the main part of both the boundary conditions and the hydrodynamic particle size in the Stokes equation, which have so far been ill-defined in any molecular theories.