The pre-processed and “cleaned” functional data was registered to the individual's structural scan and standard space images using FNIRT, then optimized using boundary-based-registration approach [28], and finally spatially smoothed using an isotropic Gaussian kernel of 6 mm full width at half maximum (FWHM).