DFT studies on nonlinear optical properties of neutral nest-shaped heterothiometallic [MOS3Py5Cu3X] (M=Mo, W; X=F, Cl, Br, I) clusters.