A method for calculating the various components of the magnetically induced current-density tensor using gauge-including atomic orbitals is described. The method is formulated in the framework of analytical derivative theory, thus enabling implementation at the Hartree-Fock self-consistent-field (HF-SCF) as well as at electron-correlated levels. First-order induced current densities have been computed up to the coupled-cluster singles and doubles level (CCSD) augmented by a perturbative treatment of triple excitations [CCSD(T)] for carbon dioxide and benzene and up to the full coupled-cluster singles, doubles, and triples (CCSDT) level in the case of ozone. The applicability of the gauge including magnetically induced current method to larger molecules is demonstrated by computing first-order current densities for porphin and hexabenzocoronene at the HF-SCF and density-functional theory level. Furthermore, a scheme for obtaining quantitative values for the induced currents in a molecule via numerical integration over the current flow is presented. For benzene, a perpendicular magnetic field induces a (field dependent) ring current of 12.8 nA T(-1) at the HF-SCF level using a triple-zeta basis set augmented with polarization functions (TZP). At the CCSD(T)/TZP level the induced current was found to be 11.4 nA T(-1). Gauge invariance and its relation to charge-current conservation is discussed.

The magnetically induced current densities for ring-shaped hydrocarbons are studied at the density functional theory (DFT) and second-order Møller-Plesset (MP2) levels using gauge-including atomic orbitals. The current densities are calculated using the gauge-including magnetically induced current approach. The calculations show that all studied hydrocarbon rings sustain strong diatropic and paratropic ring currents when exposed to an external magnetic field, regardless whether they are unsaturated or not. For nonaromatic rings, the strength of the paratropic current flowing inside the ring is as large as the diatropic one circling outside it, yielding a vanishing net ring current. For aromatic molecules, the diatropic current on the outside of the ring is much stronger than the paratropic one inside, giving rise to the net diatropic ring current that is typical for aromatic molecules. For antiaromatic molecules, the paratropic ring-current contribution inside the ring dominates. For homoaromatic molecules, the diatropic current circles at the periphery of the ring. The ring current is split at the CH(2) moiety; the main fraction of the current flow passes outside the CH(2) at the hydrogens, and some current flows inside the carbon atom. The diatropic current does not take the through-space short-cut pathway, whereas the paratropic current does take that route. Calculations of the ring-current profile show that the ring current of benzene is not transported by the pi electrons on both sides of the molecular ring. The strongest diatropic ring current flows on the outside of the ring and in the ring plane. A weaker paratropic current circles inside the ring with the largest current density in the ring plane. Due to the ring strain, small unconjugated and saturated hydrocarbon rings sustain a strong ring current which could be called ring-strain current. Nuclear magnetic shieldings calculated for 1,3,5-cycloheptatriene and homotropylium at the DFT and MP2 levels agree well with experimental values.

No abstract

The gauge-including magnetically induced current method for calculating the components of the current-density tensor using gauge-including atomic orbitals has been extended to treating open-shell molecules. The applicability of the method is demonstrated by calculations of first-order induced current densities on cyclobutadiene, Al(3), and B(3) at correlated ab initio levels of theory. For comparison, current-density calculations were also performed on the lowest closed-shell singlet state of cyclobutadiene as well on the closed-shell Al(3)(-) and B(3)(-) anions. The ring-current susceptibilities of the open-shell species are computed at the Hartree-Fock self-consistent-field, second-order Møller-Plesset perturbation theory, and coupled-cluster singles and doubles levels, whereas for the closed-shell systems also density functional theory calculations were employed. Explicit values for the current strengths caused by α and β electrons as well as the difference, representing the spin current, were obtained by numerical integration of the current-density contributions passing a plane perpendicular to the molecular ring. Comparisons of the present results to those recently obtained for the lowest triplet state of biphenyl emphasize that electron correlation effects must be considered for obtaining an accurate description of spin-current densities.

No abstract