In summary, we have developed a method to calculate intermolecular dispersion interaction based on a fragmentation technique. We represented the quantum fluctuations of electron densities for these fragments by anisotropic quantum Drude oscillators. It turns out that the method's success heavily depends on the type of fragmentation used. This aspect has been tested with comparing the monomer polarizability with fragment polarizability. We applied our methodology on S66X8 database to understand its efficacy as well as drawbacks. We find that the method works less accurately only for stacked dimers where the electron density is spreadout and therefore highly polarizable. In general, however, this method provides a simple yet effective approximation of dispersion interactions in large molecular systems with accuracy at least at par with popular D2 method. Particularly, we find that except 8 systems, QDO-corrected interaction energy for all other complexes lies within 1 kcal./mol of CCSD(T) energy. We also find that our correction due to QDO is maximum for dispersionbound complexes. For such systems, the QDO energy almost provides 100% binding energy.