Toward the comprehensive calculations on the relationship between 1 H, 13 C, 31 P chemical shifts, 2 JPH , and the bonding structure of different phosphoryl benzamides.

A comprehensive investigation was performed on 1 H, 13 C, and 31 P nuclear magnetic resonance (NMR) chemical shifts (CSs) of phosphoryl benzamide derivatives (C6 H5 C(O)NHP(O)R1 R2 ), (R1 , R2  = aziridine [L1 ], azetidine [L2 ], pyrrolidine [L3 ], piperidine [L4 ], azepane [L5 ], 4-methylpiperidine [L6 ], propane-2-amine [L7 ], and 2-methylpropane-2-amine [L8 ]) by the gauge-independent atomic orbital method (GIAO) to find the most accordant level of theory with the experimental values. To achieve this goal, all the structures were optimized using the B3LYP, BP86, PBE1PBE, M06-2X, MPWB1K, and MP2 methods with 6-31+G* basis set. Computed structural parameters demonstrate that BP86 has the best agreement to the experimental values between the other methods. The def2-TZVP and aug-cc-pVDZ basis sets were also employed to inspect the effect of different types of basis sets with higher polarization and diffuse functions. The correlation between the empirical and computational values attests that 6-31+G* basis set is the optimum case regarding minimization of the costs and results. The comparison between calculated and experimental CSs at all mentioned combinations illustrated that in accordance with structural results, the best level of theory in CSs is also BP86/6-31+G*. Besides, 2 JPH values were computed with an acceptable agreement to experimental data at the optimum level of theory. The dependency between 2 JPH and the bonding structure of studied ligands was also scrutinized by the Natural Bond Orbital (NBO) analysis that interprets the relationship between the electronic properties and 2 JPH values.