Balazic Katja, Stare Jernej, Mavri Janez
National Institute of Chemistry, Hajdrihova 19, 1000, Ljubljana, Slovenia.
J Chem Inf Model. 2007 May-Jun;47(3):832-9. doi: 10.1021/ci600491p. Epub 2007 May 1.
In this article we studied the strong intramolecularly hydrogen-bonded system 4-methoxypicolinic acid N-oxide. The potential energy surface V = V(rOH,rOO) and the corresponding dipole moment function were calculated using the DFT B3LYP/6-31+G(d,p) level of approximation. The time-independent vibrational Schrödinger equation was solved using a rectangular grid basis set and shifted Gaussian basis set. The vibrational spectrum and metric parameters were also calculated. Effects of deuteration were considered. The calculated vibrational spectra were compared with the experimental spectra. The vibrational transition corresponding to asymmetric O-H stretching that occurs at about 1400 cm-1 compares well with the experimentally assigned O-H asymmetric stretching band centered at 1380 cm-1. The corresponding asymmetric O-D stretching band was predicted to be at 1154 cm-1, while the experimental O-D band was not assigned due to its very low intensity. Several overtones and hot transitions of significant intensities were located in the vicinity of the fundamental O-H stretching frequency, effectively broadening the infrared absorption attributed to the O-H stretching mode. This is in a good agreement with the observed broad protonic absorptions found in the infrared spectra of the title compound and its analogs. We have shown that the Gaussian basis set is the method of choice for a two-dimensional vibrational problem that requires several hundreds of vibrational basis functions and when high accuracy of the eigenvalues is required or when extending the calculations to more vibrational degrees of freedom. We have also demonstrated that for a large number of basis functions the Gramm-Schmidt orthogonalization procedure outperforms symmetric and canonical orthogonalization schemes.
在本文中,我们研究了强分子内氢键体系4-甲氧基吡啶-N-氧化物。使用DFT B3LYP/6-31+G(d,p)近似水平计算了势能面V = V(rOH,rOO)和相应的偶极矩函数。使用矩形网格基组和移位高斯基组求解了与时间无关的振动薛定谔方程。还计算了振动光谱和度量参数。考虑了氘代的影响。将计算得到的振动光谱与实验光谱进行了比较。在约1400 cm-1处出现的对应于不对称O-H伸缩的振动跃迁与实验指定的以1380 cm-1为中心的O-H不对称伸缩带吻合良好。预测相应的不对称O-D伸缩带位于1154 cm-1处,而实验O-D带由于强度非常低未被指定。在基本O-H伸缩频率附近发现了几个具有显著强度的泛音和热跃迁,有效地拓宽了归因于O-H伸缩模式的红外吸收。这与在标题化合物及其类似物的红外光谱中观察到的宽质子吸收非常一致。我们已经表明,对于需要数百个振动基函数的二维振动问题,当需要高精度的本征值或当将计算扩展到更多振动自由度时,高斯基组是首选方法。我们还证明,对于大量基函数,Gramm-Schmidt正交化过程优于对称和规范正交化方案。