Abstract
The gauge‐including atomic orbital (GIAO) method for the gauge‐invariant calculation of nuclear magnetic shielding constants has been implemented at the coupled‐cluster singles and doubles (CCSD) level. A brief description of the theory and its computational requirements is provided. Finally, the GIAO–CCSD method is applied to calculate nuclear shielding constants of N2O with two different basis sets, the larger of which contains 153 contracted Gaussian functions.