I calculated the NMR spectra (post-optimization) for ethylene glycol (HOCH2CH2OH) and the deuterated one (HOCD2CD2OH) in Gaussian09 using DFT (BLYP/6-311++G(d,p)). I have attached the input and output files for the deuterated case. Now I don't understand why I am getting peaks for deuterium in proton spectra. Is there any mistake in my input file ?