Gauge invariant shielding constants calculations of ?H, ?C, and ??N were calculated for acetonitrile in the gas and liquid phases. Dierent basis sets as well as dierent ab initio and DFT methods were tested to select a time-ecient level of theory with reasonable accuracy. The eect of nuclear motion on the shielding constants was also explored. To investigate solvent eects on the shielding constants of acetonitrile, dierent clusters were extracted from molecular dynamics simulations. Convergence to the experimental values varied for the dierent clusters. The geometry of the central molecule in a cluster played an important factor in reaching convergence. / Master of Science