Nuclear quantum effects in calculated NMR shieldings of benzene; a Feynman path integral study

The Feynman path integral Monte Carlo approach has been coupled to the gauge including atomic orbital formalism in order to analyse the absolute magnetic shieldings of the benzene nuclei under the conditions of thermal equilibrium. The Hamiltonian employed in the derivation of ensemble averaged NMR quantities is of the Hartree-Fock type. The basis set used is of 6–31G quality. The spatial delocalization of the atoms leads to a deshielding of both types of benzene nuclei relative to the shieldings experienced at the minimum of the potential energy surface. This deshielding has to be traced back to bond length elongations in thermal equilibrium. The influence of the nuclear fluctuations on the NMR parameters of benzene is quantum driven up to temperatures of 400 K; classical fluctuations are of minor importance in this low-temperature window.     

Ground state and excited state properties of ethylene isomers studied by a combined Feynman path integral-ab initio approach

Ground state and excited state properties of ethylene, C₂H₄, and several ethylene isomers have been studied by Feynman path integral Monte Carlo (PIMC) simulations. The PIMC treatment of the atomic nuclei has been combined with different electronic Hamiltonians in order to analyse the influence of the nuclear degrees of freedom on electronic quantities. Electronic expectation values at the minimum of the potential energy surface (PES) have been compared with PIMC based ensemble averaged values. Ensemble averaged quantities have been derived by Hamiltonians of the ab initio type and a tight-binding (TB) one-electron model. The combined influence of anharmonicities in the interatomic potential and the quantum fluctuations of the atomic nuclei lead to ensemble averaged bondlengths r _g which are significantly larger than the parameters r _e, at the minimum of the PES. The implications of this bond length elongation for the electronic properties of ethylene are discussed. The occupied canonical molecular orbitals (CMOs) of ethylene are destabilized under the influence of the nuclear degrees of freedom while virtual CMOs are stabilized. These shifts of one-electron energies suggest a comparison of electronic excitation energies at the minimum of the PES with PIMC based ensemble averages. The quantum fluctuations of the nuclei cause a strong redistribution in the intensities of electronic transitions. Transitions, which are dipole allowed in the planar D_2h geometry of ethylene, lose intensity under the influence of nuclear quantum effects, and vice versa for electronic excitations that are dipole forbidden under D_2h symmetry. This ‘vibrational borrowing’ is enhanced with decreasing atomic masses. The Feynman centroid density has been used to calculate the anharmonic vibrational wavenumbers of C₂H₄ and C₂D₄. The results of the present PIMC simulations have been employed to emphasize general problems of electronic structure calculations based on a single nuclear configuration (i.e. the configuration at the minimum of the PES).