The structures of a series of beryllium containing complexes have been optimized at the B3LYP/6-31G(d) level and their Be-9 magnetic shielding values have been determined using B3LYP/6-311G+g(2d,p) and the gauge-including atomic orbital (GIAO) method. The calculated chemical shifts are in excellent agreement with experimental values. The performance of a variety of NMR methods (SGO, IGAIM, CSGT) were also examined but were found to be inferior to the GIAO method at the chosen level of theory employed. The theoretical method has been utilized to predict the beryllium chemical shifts of structurally characterized complexes for which no measured Be-9 NMR spectrum exists, and to investigate a literature complex with an unusual Be-9 NMR chemical shift. A new standard for beryllium NMR in nonaqueous solvents has been suggested.