Nuclear magnetic resonance chemical shifts and quadrupole couplings for different hydrogen-bonding cases occurring in liquid water:: A computational study

Nuclear magnetic resonance (NMR) parameters are determined theoretically for the oxygen and hydrogen/deuterium nuclei of differently hydrogen-bonded water molecules in liquid water at 300 K. The parameters are the chemical shift, the shielding anisotropy, the asymmetry parameter of shielding, the nuclear quadrupole coupling constant, and the asymmetry parameter of the nuclear quadrupole coupling. We sample instantaneous configurations from a Car-Parrinello molecular dynamics simulation and feed nuclear coordinates into a quantum chemical program for the calculation of NMR parameters using density-functional theory with the three-parameter hybrid exchange-correlation (B3LYP) functional. In the subsequent analysis, molecules are divided into groups according to the number of hydrogen bonds they possess, and the full average NMR tensors are calculated separately for each group. The classification of the hydrogen-bonding cases is performed using a simple distance-based criterion. The analysis reveals in detail how the NMR tensors evolve as the environment changes gradually from gas to liquid upon increasing the number of hydrogen bonds to the molecule of interest. Liquid-state distributions of the instantaneous values of the NMR properties show a wide range of values for each hydrogen-bonding species with significant overlap between the different cases. Our study shows how local changes in the environment, along with classical thermal averaging, affect the NMR parameters in liquid water. For example, a broken or alternatively extra hydrogen bond induces major changes in the NMR tensors, and the effect is more pronounced for hydrogen or deuterium than for oxygen. The data sheds light on the usefulness of NMR experiments in investigating the local coordination of liquid water.