Accurate predictions for molecular magnetisabilities and rotational g tensors using a simple and efficient DFT approach

Efficient calculation of the rotational g tensor from auxiliary density functional theory

The computation of the rotational g tensor with the recently developed auxiliary density functional theory (ADFT) gauge including atomic orbital (GIAO) methodology is presented. For the rotational g tensor, the calculation of the magnetizability tensor represents the most demanding computational task. With the ADFT-GIAO methodology, the CPU time for the magnetizability tensor calculation can be dramatically reduced. Therefore, it seems most desirable to employ the ADFT-GIAO methodology also for the computation of the rotational g tensor. In this work, the quality of rotational g tensors obtained with the ADFT-GIAO methodology is compared with available experimental data as well as with other theoretical results at the Hartree-Fock and coupled-cluster level of theory. It is found that the agreement between the ADFT-GIAO results and the experiment is good. Furthermore, we also show that the ADFT-GIAO g tensor calculation is applicable to large systems like carbon nanotube models containing hundreds of atom and thousands of basis functions.

On equilibrium structures of the water molecule

Equilibrium structures are fundamental entities in molecular sciences. They can be inferred from experimental data by complicated inverse procedures which often rely on several assumptions, including the Born-Oppenheimer approximation. Theory provides a direct route to equilibrium geometries. A recent high-quality ab initio semiglobal adiabatic potential-energy surface (PES) of the electronic ground state of water, reported by Polyansky et al. [ ibid. 299, 539 (2003)] and called CVRQD here, is analyzed in this respect. The equilibrium geometries resulting from this direct route are deemed to be of higher accuracy than those that can be determined by analyzing experimental data. Detailed investigation of the effect of the breakdown of the Born-Oppenheimer approximation suggests that the concept of an isotope-independent equilibrium structure holds to about 3 x 10(-5) A and 0.02 degrees for water. The mass-independent [Born-Oppenheimer (BO)] equilibrium bond length and bond angle on the ground electronic state PES of water is r(e) (BO)=0.957 82 A and theta e (BO)=104.48(5) degrees , respectively. The related mass-dependent (adiabatic) equilibrium bond length and bond angle of H2 (16)O is r(e) (ad)=0.957 85 A and theta e (ad)=104.50(0) degrees , respectively, while those of D2 (16)O are r(e) (ad)=0.957 83 A and theta e (ad)=104.49(0) degrees . Pure ab initio prediction of J=1 and 2 rotational levels on the vibrational ground state by the CVRQD PESs is accurate to better than 0.002 cm(-1) for all isotopologs of water considered. Elaborate adjustment of the CVRQD PESs to reproduce all observed rovibrational transitions to better than 0.05 cm(-1) (or the lower ones to better than 0.0035 cm(-1)) does not result in noticeable changes in the adiabatic equilibrium structure parameters. The expectation values of the ground vibrational state rotational constants of the water isotopologs, computed in the Eckart frame using the CVRQD PESs and atomic masses, deviate from the experimentally measured ones only marginally, especially for A0 and B0. The small residual deviations in the effective rotational constants are due to centrifugal distortion, electronic, and non-Born-Oppenheimer effects. The spectroscopic (nonadiabatic) equilibrium structural parameters of H2 16O, obtained from experimentally determined A'0 and B'0 rotational constants corrected empirically to obtain equilibrium rotational constants, are r(e) (sp)=0.957 77 A and theta e (sp)=104.48 degrees .