Hi Teo and Matt, Thanks for pointing out the problem. I'll read more carefully into those references. The reason I didn't use the periodic and decoupling recoupling scheme is that the QM region should be isolated and should interact with its images. Regards, Shihao