How does relativity affect magnetically induced currents?

Under terrestrial conditions magnetic fields are omnipresent and even in environments free of external magnetic fields nuclear magnetic moments are present. The electronic probability current density vector field (jB) is the fundamental property for the description of the magnetic response of a molecule to an external magnetic field B. Various programs for the computation of numerical approximations to jB have been developed, and exemplary studies have been reviewed some years ago. Other magnetic response properties like the magnetic susceptibility, induced magnetic multipole-moments or even chemical shieldings can directly be calculated from jB. 5 jB can be determined experimentally, as it was shown in two conceptional studies using polarized neutron diffraction. Two of us (SK, MR) have developed in the framework of the program ReSpect methods for a fully relativistic four-component treatment (FR) capable of calculating jB in polyatomic molecules. In this current version of ReSpect also a physically consistent determination of the influence of solely the relativistic spin–orbit coupling (SOC) contribution to jB is possible. This is one of the first investigations where particular SOC contributions to FR molecular currents are reported. The only other work where a FR method was employed to analyze the SOC contribution to jB involved the Sternheim approximation to the diamagnetic contribution of the response current. Later Sulzer and coworkers overcame this approximation by using the simple magnetic balance approach but they have not regarded the specific nature of the SOC contributions they yielded. In this work we use a computationally efficient solution based on the restricted magnetic balance framework. We have chosen gold(I) hydride (AuH) and mercury(II) bishydride (HgH2) as subjects of this study due to their simplicity and since they are in the non-relativistic (NR) framework single reference singlets. Not taking SOC into account (computational details are provided in the ESI†), AuH sustains one large gold centered diatropic (inductively weakening the external field B) current loop (see Fig. 1b) with a distinct extension directed to the proton. At the FR level (Fig. 1a) the current topology remains unchanged, however, it appears that the current loops contain considerably more curvature around the atomic cores and most notably a kind of current cusp appears at the position of the proton, which without regarding SOC effects (non-SOC) is passed over by a smooth current. Such a cusp structure in molecular ring currents, to the best of our knowledge is previously unobserved. The current difference plot (FR minus non-SOC, Fig. 1c) shows that this cusp originates solely from SOC, causing a strongly localized ring current around the proton. It can be easily seen that a ring current of suitable orientation when added to a smooth current field with little curvature yields a total current of such a cusp shape. However, it is required that the strength of the ring current contribution in the area of the cusp is approximately constant and equal in strength to the smooth current field contribution. Apparently, these conditions are met in the case of AuH. In a homogeneous magnetic field of 1 T strength the total current of the SOC induced vortex around the hydrogen atom integrates to 1.6 nA (see the ESI† for computational details). This additional diatropic vortex causes an SOC induced high field chemical shift of the proton of about 19 ppm compared to the non-SOC calculation (FR shielding: 54 ppm, non-SOC shielding: 35 ppm). The corresponding difference (FR minus non-SOC) of only the z components of the shielding tensor, which can be directly a Chemistry of Materials, Paris-Lodron-University Salzburg, Hellbrunner Str. 34, A-5020 Salzburg, Austria. E-mail: Raphael.Berger@sbg.ac.at b Department of Chemistry University of Tromsø, N-9037 Tromsø, Norway. E-mail: stanislav.komorovsky@uit.no † Electronic supplementary information (ESI) available: Computational details. See DOI: 10.1039/c5cc05732a ‡ Current address: Laboratory for Instruction in Swedish, University of Helsinki, A. I. Virtanens Plats 1, FI-00014 University of Helsinki, Finland. Received 10th July 2015, Accepted 22nd July 2015

Under terrestrial conditions magnetic fields are omnipresent and even in environments free of external magnetic fields nuclear magnetic moments are present.The electronic probability current density vector field (j B ) is the fundamental property for the description of the magnetic response of a molecule to an external magnetic field B. Various programs for the computation of numerical approximations to j B have been developed, [1][2][3] and exemplary studies have been reviewed some years ago. 4Other magnetic response properties like the magnetic susceptibility, induced magnetic multipole-moments or even chemical shieldings can directly be calculated from j B . 5 j B can be determined experimentally, as it was shown in two conceptional studies using polarized neutron diffraction. 6,7wo of us (SK, MR) have developed in the framework of the program ReSpect 8 methods for a fully relativistic four-component treatment (FR) capable of calculating j B in polyatomic molecules.In this current version of ReSpect also a physically consistent determination of the influence of solely the relativistic spin-orbit coupling (SOC) contribution to j B is possible.This is one of the first investigations where particular SOC contributions to FR molecular currents are reported.The only other work where a FR method was employed to analyze the SOC contribution to j B involved the Sternheim approximation 9 to the diamagnetic contribution of the response current. 10Later Sulzer and coworkers overcame this approximation by using the simple magnetic balance approach but they have not regarded the specific nature of the SOC contributions they yielded. 11In this work we use a computationally efficient solution based on the restricted magnetic balance framework. 12We have chosen gold(I) hydride (AuH) and mercury(II) bishydride (HgH 2 ) as subjects of this study due to their simplicity and since they are in the non-relativistic (NR) framework single reference singlets.
Not taking SOC into account (computational details are provided in the ESI †), AuH sustains one large gold centered diatropic (inductively weakening the external field B) current loop (see Fig. 1b) with a distinct extension directed to the proton.At the FR level (Fig. 1a) the current topology remains unchanged, however, it appears that the current loops contain considerably more curvature around the atomic cores and most notably a kind of current cusp appears at the position of the proton, which without regarding SOC effects (non-SOC) is passed over by a smooth current.Such a cusp structure in molecular ring currents, to the best of our knowledge is previously unobserved.The current difference plot (FR minus non-SOC, Fig. 1c) shows that this cusp originates solely from SOC, causing a strongly localized ring current around the proton.It can be easily seen that a ring current of suitable orientation when added to a smooth current field with little curvature yields a total current of such a cusp shape.However, it is required that the strength of the ring current contribution in the area of the cusp is approximately constant and equal in strength to the smooth current field contribution.Apparently, these conditions are met in the case of AuH.In a homogeneous magnetic field of 1 T strength the total current of the SOC induced vortex around the hydrogen atom integrates to 1.6 nA (see the ESI † for computational details).This additional diatropic vortex causes an SOC induced high field chemical shift of the proton of about À19 ppm compared to the non-SOC calculation (FR shielding: 54 ppm, non-SOC shielding: 35 ppm).The corresponding difference (FR minus non-SOC) of only the z components of the shielding tensor, which can be directly (without taking into account directional averaging) related to the integrated and plotted j B field, is À30 ppm.To quantify the additional curvature of the current field in the region around the hydrogen atom one can calculate the radius of a circular loop of an infinitely thin wire, causing a field of À30 ppm of 1 T (3.0 Â 10 À5 T) at the center of the loop (hydrogen atom position) and sustaining a current of 1.6 nA using the Biot-Savart law.This yields a loop radius of 0.35 Å or a curvature (reciprocal radius) of 2.9 Å À1 .
At the molecular currents of HgH 2 at the non-SOC level of theory (Fig. 1e), we find very similar conditions to those in AuH with the exception that the extension of the currents to the proton positions is slightly more pronounced than in AuH.Going to the full relativistic level of theory we again see that more curvature is introduced into the streamlines (Fig. 1d).Unlike in AuH here a drop-shaped deformation of the currents, directed to the center of the molecule, is induced by SOC effects.Again its origins can be readily seen in the difference plot in Fig. 1f (FR minus non-SOC), showing small SOC induced localized current loops on the proton positions.In contrast to the situation in AuH these additional local current loops in HgH 2 are paratropic (inductively enhancing the external field B), leading to drop-shaped perturbations rather than the cusp-shaped ones seen in AuH.An integration (see the ESI † for computational details) of these SOC induced local current vortices yields a local paramagnetic contribution of À1.5 nA (in a magnetic field of 1 T).At the SOC level of theory the isotropic 1 H chemical shieldings (10 ppm) are found to be down field shifted by +17 ppm as compared to their non-SOC equivalents (27 ppm).An analogous consideration of the curvature as for AuH yields a loop radius of 0.41 Å or a curvature of 2.4 Å À1 .
For an explanation of these findings the following can be considered: j B in a single charged spin-free particle in the NR framework contains a sum of two terms which are conventionally called diamagnetic and paramagnetic terms, together these might be called 'Pauli-current' § (j Pauli B ).In the case of a manyparticle system described in an effective single particle picture like HF or Kohn-Sham DFT with the total spin equal to zero, like for example a closed shell molecule, j B reduces to the sum of Pauli-current terms of each particle.In the NR framework there are simple MO selection rules allowing for a prediction and a chemical interpretation of these currents, 13,14 and the physics of the NR contribution to the molecular currents appears to be well understood.For systems with non-zero spin, and again in the NR case, j B is the sum of the Pauli-current terms plus a 'spin-current' term j spin : Hereby j spin is the curl of the density of the magnetic moment l spin arising from spin: The density of the magnetic moment arising from spin, l spin , is directly proportional to the spin-density ¶ distribution q spin of the system: This means that l spin (or q spin ) is represented by a current field j spin and both are related by eqn (2).In classical electrodynamics (using j instead of j spin , and B instead of l spin ) eqn ( 2) is Ampe `res law.It relates currents to the magnetic field they induce.It can be shown that in the NR case, eqn (1) also holds for a system of n-electrons with non-zero total spin in an effective one-particle picture (summing up one-particle We are here interested in a fully relativistic description including spin-orbit effects of currents in a magnetic field j FR B for ''closed shell'' molecules (i.e.singlet states in the NR picture).Treating SOC and B as a perturbation and using an effective one-particle description, it can be shown 15 that up to second order terms, the total current where j spin B,SOC is a spin-current which is induced into the closed shell, spin-zero system by the presence of both the magnetic field and SOC, while j B,scalar is the scalar relativistic remainder of the induced current density.Note that first order SO effects contribute exclusively via j spin to the response current density, while scalar relativistic effects contribute via j spin and j Pauli like terms (in the sense of the Gordon decomposition of the fourcomponent current 16 ).Due to the short range of the SOC effect the perturbations of the singlet wave function in the first instance are spatially in close proximity to the heavy nuclei causing the SOC effect, but when such perturbations occur within a covalent bond between the heavy atom and a light atom the perturbation extends also to a spatial region around the bound light atom.8][19] In the presence of an external magnetic field B and under consideration of SOC a non-zero and non-constant spin-density is induced into the wave function.The spin-current that directly corresponds to this induced spin-density is, in very good approximation, the only spin-orbit contribution to the currents (see the ESI † for details of term partitioning) (where we have used the same definitions as for eqn (4) and j non-SOC B = j Pauli B + j B,scalar ).Thus studying pure SOC influences on magnetically induced currents in 'closed-shell' molecules can be reduced to studying the SOC induced spin-density in the magnetic field q spin B,SOC .The SOC coupling mechanism explains why the effects are confined to a surrounding of the nuclei and why they can extend to covalently bond lighter atoms as we have observed in AuH and HgH 2 .
We have developed computational methods to calculate four-component relativistic magnetically induced molecular ring currents.Case studies on AuH and HgH 2 illustrate that spin-orbit coupling influences the shape of induced currents leading to intriguing curvature and strongly localized ring current contributions to hydrogen atoms bound to gold and mercury.In this way we have revealed the ring current picture of the HALA effect.SOC induced currents in the magnetic field can be traced back directly to an induced spin-density via a relation resembling Ampe `res law.In general a non-zero and non-constant spin-density arises in 'closed shell' molecules due to the presence of both a magnetic field and SOC.
Notes and references § Alternative definitions of the term 'Pauli current' are in use.This one might be justified by the notion that it is the current terms which arise already from the Pauli equation. 20,21Which is defined as q spin = h/2 hri, integrating over spin coordinates.

Fig. 1
Fig. 1 Streamline plots in AuH and HgH 2 .B is set perpendicular to the plot plane, all atoms are in the plane.Diatropic currents (weakening B) show clockwise and paratropic currents (enhancing B) have anti-clockwise orientation.The line shading intensity is proportional to the current strengths.(a) AuH, fully relativistic; (b) AuH, without spin-orbit coupling; (c) AuH, difference plot from (a) and (b); (d) HgH 2 , fully relativistic; (e) HgH 2 , without spin-orbit coupling; (f) HgH 2 , difference plot from (d) and (e).