Journal of Physical Chemistry A, Vol.107, No.38, 7629-7642, 2003
Accurate C-13 and N-15 chemical shift and N-14 quadrupolar coupling constant calculations in amino acid crystals: Zwitterionic, hydrogen-bonded systems
EIM (embedded ion method), cluster, combined EIM/cluster, and isolated molecule C-13 and N-15 chemical shielding and quadrupolar coupling constant (QCC) calculations at the B3LYP level with D95**, D95++**., 6-311G**, and 6-311+G** basis sets were done on the amino acids L-alanine, L-asparagine monohydrate, and L-histidine monohydrate monohydrochloride and on the two polymorphs alpha and gamma glycine. The intermolecular interactions that are present in the amino acid crystals are accounted for in the EIM calculations by a finite array of point charges calculated from Ewald lattice sums and in the cluster calculations by a shell of neighboring molecules or molecular fragments. The combined EIM/cluster calculations utilize a cluster of molecules inside an EIM point charge array. The theoretical C-13 and N-15 principal shielding values for the amino acids studied are compared to the experimental principal shift values. In addition, theoretical CN bond orientations in the chemical shift principal axis system (PAS) are compared to the experimental orientations obtained from C-13-N-14 dipolar couplings. The theoretical QCC at the nitrogen positions are compared to experimental N-14 QCC principal values reported in the literature. The carbon and nitrogen theoretical chemical shielding, the C-N orientations, and the QCCs from the ab initio calculations show improved agreement with the experimental values when the intermolecular interactions are accounted for by EIM or cluster calculations. The EIM C-13 shielding calculations are found to give better agreement with the experimental values than cluster C-13 shielding calculations. However, to achieve good agreement between the theoretical N-14 QCC and the N-15 principal shielding values with the respective experimental values, both intermolecular electrostatic and covalent interactions have to be included explicitly in the EIM/cluster calculations.