Spatial Contributions to 1 H NMR Chemical Shifts of Free-Base Porphyrinoids

: A recently developed methodology for calculating, analyzing, and visualizing nuclear magnetic shielding densities is used for studying spatial contributions including ring-current contributions to 1 H nuclear magnetic resonance (NMR) chemical shifts of aromatic and anti-aromatic free-base porphyrinoids. Our approach allows a visual inspection of the spatial origin of the positive (shielding) and negative (deshielding) contributions to the nuclear magnetic shielding constants. Diatropic and paratropic current-density ﬂuxes yield both shielding and deshielding contributions implying that not merely the tropicity of the current density determines whether the contribution has a shielding or deshielding character. Instead the shielding or deshielding contribution is determined by the direction of the current-density ﬂux with respect to the studied nucleus. The present study shows the power of the present approach that can be used for determining spatial contributions to 1 H NMR magnetic shieldings. The shielding constants can be assigned to atoms by numerically integrating the MSD in atomic domains. The calculations yield the magnetically induced ring current contributions to the shielding constants of the the inner and outer hydrogen atoms of typical aromatic and anti-aromatic porphyrinoids. Here, we have focused on 1 H NMR shielding constants. However, a similar analysis can also be performed for 13 C NMR and 15 N NMR shielding constants.


Introduction
This article is dedicated to Professor Riccardo Zanasi on the occasion of his 70th birthday. He graduated from the University of Modena, Italy in 1975, studying electric and magnetic properties of molecules. He is one of the pioneers in the field of current-density calculations with his first article from 1981 on the ring-current model of the cyclopropenyl cation [1]. The "continuous transformation of the origin of the current density (CTOCD) method" from 1994 [2], which Keith and Bader called "continuous set of gauge transformations (CSGT)", was the starting point for the SYSMO and SYSMOIC program packages [3,4]. Professor Zanasi has published about 70 scientific articles on computational studies of magnetically induced current densities in molecules with "current density" in the title. He is the main developer of the SYSMO code. The topic of this article is closely related to the research interests of Professor Zanasi, since we have calculated and analyzed magnetic shielding densities (MSD) of free-base porphyrinoids. MSD is the spatial contribution to the chemical shifts in nuclear magnetic resonance (NMR) spectroscopy and calculated using the magnetically induced current density ( Figure S1).
For the aromatic benzene molecule, it is well known that the magnetically induced ring current leads to a downfield shift of the 1 H NMR signal, because the diatropic ring current induces a secondary magnetic field with field lines forming closed loops around the ring current. The field lines point in the opposite direction to the applied magnetic field inside the ring, while outside the ring they are parallel to the external field, see Figure 1. The hydrogen atoms of benzene outside the ring-current pathway are deshielded because the induced magnetic field enhances the applied one. However, a recent study of the MSD of benzene showed that the ring current leads to regions with significant shielding and deshielding contributions [5]. The ring-current contribution to the 1 H NMR chemical shift of benzene can be estimated to be about 1.7 ppm from the experimental chemical shifts of benzene (5.6 ppm) and the vinylic hydrogen atoms of cyclohexene (7.3 ppm). Figure 1. Illustration of the induced magnetic field (black arrows) when the aromatic benzene is exposed to an external magnetic field (yellow) perpendicular to the molecular ring. The magnetically induced current density circling in the classical direction is called diatropic (blue) and the one flowing in the opposite direction is paratropic (red). The induced magnetic field in the center of the ring points in the opposite direction relative to the applied field, while it deshields the nuclear magnetic moment of the hydrogen atoms outside the ring.
The origin of the 1 H NMR chemical shifts is even more complex when studying aromatic or anti-aromatic molecules bearing hydrogen atoms inside the ring-current pathway. The aromatic free-base trans-porphyrin is an example of such molecules. The accepted notion is that the 1 H NMR signal of the inner hydrogen is upfield shifted and the 1 H NMR signals of the meso and β hydrogen nuclei are downfield shifted because the direction of the induced magnetic field is in the opposite direction to the applied one inside the porphyrin ring and parallel to the external field outside the ring. The chemical shifts of the 1 H NMR signals of anti-aromatic molecules are influenced in the opposite manner compared to aromatic ones because anti-aromatic molecules sustain strong paratropic ring currents. Measurements of the 1 H NMR chemical shifts of the inner or outer hydrogen atoms of porphyrinoids are often used for assessing their aromatic character experimentally. Calculated 1 H NMR chemical shifts can be used for estimating the direction and the strength of the magnetically induced ring current and for identifying porphyrinoid tautomers [6][7][8].
The ring-current model of aromatic and anti-aromatic molecules [9,10] with a currentdensity pathway forming a closed loop in the π orbitals above and below the molecular ring might be a somewhat oversimplified explanation of the aromatic or anti-aromatic contribution to the 1 H NMR chemical shifts. Since the vector potential of the nuclear magnetic moment declines rapidly with the distance from the nucleus, the largest contributions to the 1 H NMR chemical shifts originate from the current density near the considered nucleus.
Instead of the well known nucleus-independent-chemical-shift (NICS) functions [11][12][13], we investigate the MSD calculated using the magnetically induced current density. The MSD yields detailed information about spatial contributions to the NMR chemical shifts. The MSD calculated using the Biot-Savart expression for NMR chemical shifts can be more useful for interpreting experimental NMR data than alternative approaches [14], since it is possible not only to quantify the amount of positive and negative contributions to the shielding constant of a specific nucleus but also to identify the spatial origin of the most relevant contributions.
In this work, we apply our methodology for analyzing spatial contributions to the 1 H NMR magnetic shieldings of inner and outer hydrogen atoms of free-base transporphyrin, free-base trans-isophlorin, and free-base trans-norcorrole ( Figures S2 and S3). Free-base trans-porphyrin serves as the typical aromatic porphyrin. Free-base transisophlorin proposed by Woodward represents anti-aromatic porphyrins [15,16], whereas anti-aromatic norcorroles are the smallest tetrapyrrole porphyrinoids that have been synthesized [17][18][19]. The underlying theory is briefly described in Section 2. The computational levels are presented in Section 3.1 and the employed numerical methods for calculating the shielding density are described in Section 3.2. The obtained shielding densities are discussed in Section 4. Finally, the study is summarized and concluded in Section 5.

Theory
Nuclear magnetic shieldings are usually calculated using the gradient theory of electronic structure calculations as the mixed second derivative of the electronic energy with respect to the strength of the external magnetic field and the size of the nuclear magnetic moment [20][21][22][23][24]. The components of the nuclear magnetic shielding tensor, σ, in the α and β Cartesian directions can also be calculated as the second derivative of the magnetic interaction energy expressed as a spatial integral over the scalar product of the magnetically induced current-density susceptibility (∂J B (r)/∂B β ) due to the external magnetic field (B) and the first derivative of the vector potential of the nuclear magnetic moment (m I ) with respect to the size of the nuclear magnetic moment (∂A m I (r)/∂m I α ) calculated at vanishing B and m I [14,[25][26][27], The vector potential A m I (r) of the nuclear magnetic dipole moment in international system of units (SI) is given by where R I is the position of nucleus I and µ 0 is the vacuum permeability [28]. This scheme is implemented in the gauge-including magnetically induced currents (GIMIC) method [29][30][31][32][33][34] for calculating the current-density susceptibility tensor (CDT) [35,36] as described in Section 3.2. Similarly, the nuclear magnetic shielding tensor can be calculated from the current density induced by the nuclear magnetic moment multiplied with the vector potential of the external magnetic field [14,[36][37][38]. This approach has been employed in other studies [39,40], however it has not been considered in this work. Even though the NMR shieldings obtained with the gradient theory and the integration approaches are the same, the integration method can provide information about orbital and spatial contributions to a given NMR chemical shift [36,37,41]. When a numerical representation of the current density is gauge-origin independent, the calculated NMR chemical shifts obtained with the integration approach are also independent of the gauge origin. We use the GIMIC method [29][30][31][32] for calculating current densities, however, the CTOCD approach also leads to gauge-origin independence [2,4,[42][43][44]. The spatial distribution of the shielding density provides detailed information about the origin of the individual elements of the nuclear magnetic shielding tensor, as well as the shielding constants [5,[44][45][46][47][48][49][50]. Dividing the magnetic shielding density into positive and negative parts yields the spatial origins of the shielding and deshielding contributions to the shielding tensor and the isotropic shielding constants, providing a rigorous physical basis for interpreting NMR chemical shifts.

Electronic and Molecular Structure Calculations
The molecular structures of free-base trans-porphyrin, free-base isophlorin, free-base trans-norcorrole were optimized with Turbomole [51] version 7.5 employing the B3LYP density functional [52][53][54], the def2-TZVP basis set using the resolution of the identity (RI) approximation [55,56], and the m5 quadrature grid [57]. The molecular structure of isophlorin was assumed to be planar. The planar isophlorin structure is not a minimum on the potential energy surface but a saddle point with several imaginary vibrational frequencies due to out-of-plane motions of the inner hydrogens. The optimized molecular structures are shown in Figure 2 and the Cartesian coordinates are given in the electronic supporting information (ESI). The NMR shielding constants were calculated with Turbomole at the BHandHLYP (LIBXC ID 436) [53,58,59] level of theory using gauge-including atomic orbitals (GIAO) [20,[60][61][62]. The BHandHLYP functional that has 50% Hartree-Fock exchange yielded accurate magnetizabilities in a recent benchmark study [41]. For all NMR calculations the pcseg-3 basis set has been employed which has been optimized for the calculation of NMR parameters [63]. The pcseg-3 basis sets have been obtained from the EMSL basis set exchange library [64][65][66]. For the visualizations, the smaller def2-TZVP basis set has been employed which was tested to be sufficient.

Magnetic Shielding-Density Calculations
The use of finite one-particle basis sets in the calculation of CDT introduces a gauge dependence even though the exact solution of the Schrödinger equation is gauge invariant. However, the calculated CDT can be made gauge-origin independent by using GIAOs. GIAOs also lead to a fast basis-set convergence of the current density and magnetic properties, because GIAOs are correct to first order in the magnetic response for any choice of the gauge origin, whereas ordinary basis functions are correct only to zeroth order [67]. The GIAOs are given by [20,60] where i is the imaginary unit and χ (0) µ (r) is a standard Gaussian-type basis function centered at R µ . The gauge origin is eliminated from the CDT expression when using GIAOs [29,31,32]. The CDT expression of the GIMIC method is where D is the density matrix and ∂D/∂B are the magnetically perturbed density matrices in the atomic-orbital basis, αβγ is the Levi-Civita symbol with α, β, γ ∈ {x, y, z}.h(r) denotes that the singular |r − R I | −3 denominator appearing in all terms is omitted. The non-singular magnetic interaction operators are and where R I is the position of nucleus I. All terms containing the gauge origin R O cancel in Equation (4), making the CDT and the Biot-Savart expression independent of the gauge origin. All nuclear-position terms in Equation (4) also cancel, eliminating the coordinates of the nucleus I from the CDT expression. The current density depends implicitly on the nuclear positions, since the basis functions are located at the nuclei. The Biot-Savart expression for the nuclear magnetic shielding tensor of nucleus I, σ I αβ , is then where µ 0 is the vacuum permeability and α, β, γ, δ ∈ {x, y, z}. The integrand in Equation (7) is the magnetic shielding density (MSD) containing spatial contributions to the nuclear magnetic shielding of nucleus I [5,14,25,26,44,46,[48][49][50]. Visualization of the MSD yields information about the spatial origin of the NMR chemical shift of nucleus I. NMR chemical shifts can be interpreted by plotting the positive and negative parts of the integrand separately. The visualization provides information about shielding and deshielding contributions originating from the relative direction of the current density with respect to the investigated atom I [5,32,68]. Contributions to NMR shieldings can be calculated by integrating the Biot-Savart expression in Equation (7) numerically using the CDT calculated in the integration points. We have implemented a numerical integration scheme into the GIMIC program for calculating spatial contributions to nuclear magnetic shieldings and magnetizabilities [5,34,41]. Atomic contributions can be obtained by integrating atomic domains generated by the NUMGRID library [69] using Becke's multicenter scheme [70]. The atomic domains are constructed using the Becke partitioning scheme with the iteration order k = 3 in the construction of the cut-off function [70]. Angular integration of the atomcentered domains is performed with Lebedev's angular grid [71] and the radial integration grid is constructed as suggested by Lindh et al. [72]. The density matrix, the magnetically perturbed density matrices and the basis-set data are obtained from Turbomole [51] calculations of NMR shielding constants.

Free-Base trans-Porphyrin
We have analyzed the ring-current contribution to the 1 H NMR shielding constants by plotting the spatial distribution of its isotropic shielding constant (σ iso ). The 1 H NMR shielding density of trans-porphyrin in Figure 3 shows that the inner hydrogen atoms have shielding contributions in the molecular plane from the ring current along the outer edge of the molecule [73]. However, a stronger shielding contribution originates from the innermost pathway via the four pyrrole nitrogen atoms. The largest shielding contribution appears near the studied hydrogen, because of the singular denominator of the vector potential of the magnetic moment of the inner hydrogen. The nitrogen atom also sustains a local diatropic current-density vortex extending around the N-H bond. The paratropic ring current inside the porphyrin ring and the diatropic current-density vortex of the nitrogen of the adjacent pyrrole rings cause deshielding. A deshielding area is also seen at the other inner hydrogen where the paratropic ring-current flux inside the porphyrin ring is closest to the studied inner hydrogen. The MSD of σ iso of the inner hydrogen is seen 1 Å above the molecular plane in Figure 3. The effect of the ring current is weaker there than in the molecular plane. The deshielding contributions are also very small at a distance of 1 Å from the molecular plane because the paratropic ring current vanishes there. The MSD of the meso hydrogen in the molecular plane shown in Figure 3 reveal a strong shielding contribution near the hydrogen. A deshielding contribution caused by the diatropic ring current flux in the vicinity of the ipso carbon and the nearest α carbon atoms of the adjacent pyrrole rings. The more distant part of the diatropic ring current shields the nucleus of the meso hydrogen. The paratropic ring current inside the porphyrin results in a shielding contribution inside the meso carbon. One can also see alternating shielding and deshielding contributions from the paratropic ring current inside the porphyrin ring due to its relative direction with respect to the nucleus of the studied meso hydrogen. The MSD of the meso hydrogen 1 Å above the molecular plane shown in Figure 3 has a similar pattern as in the molecular plane. However, the contributions away from the molecular plane are much smaller.
The 3D picture of the MSD of the inner hydrogen in Figure 4 shows the shielding contribution from the diatropic ring current along the outer edge of the molecule and the deshielding contribution from the paratropic ring current inside the porphyrin ring. Lots of atoms contribute significantly to the isotropic shielding constant of the inner hydrogen of 37.52 ppm, because the distances from it to many of the other atoms are relatively short. The local contribution from the inner H is only 12.85 ppm. The contributions from N ipso , C α of the same pyrrole ring, the nearest C meso , and the nearest C α atoms of the adjacent pyrrole rings are 5.28, 2.35, 1.53, and 1.16 ppm, respectively. The rest of the atomic contributions to σ iso of the inner hydrogen is in the range of ±1 ppm. Contributions from each atom to the σ iso of the hydrogen atoms are reported in the ESI and summarized in Figure 5. The 3D pictures of the MSD of the meso and β hydrogen atoms are reported in the ESI.

Free-Base Isophlorin
The spatial distribution of the isotropic MSD of free-base isophlorin in Figure 6 shows that the inner hydrogen atoms and the meso hydrogen atoms have deshielding contributions from the paratropic ring current in the molecular plane along the outer edge of the molecule. The shielding contributions originate mainly from the current density in the vicinity of the studied inner hydrogen.
The 3D contour of the ring-current contribution to σ iso of the inner hydrogen is shown in Figure 7. The ring current following the inner route is stronger than the one via the β carbons which means that its contribution to the shielding density is bigger. Furthermore, the vector potential of the nuclear magnetic moment, A m I (r) declines rapidly with increasing distance from the nucleus I.
Dividing the MSD of the inner hydrogen into atomic domains that are integrated separately shows that many atoms contribute significantly to its σ iso value of −6.28 ppm. Since the distances from the inner hydrogen to adjacent atoms are short, the vector potential of the nuclear magnetic moment of the inner hydrogen has a significant amplitude at the atomic domains of the neighboring atoms leading to large atomic contributions to σ iso of the inner hydrogen.   The local contribution from the atomic domain of the inner hydrogen is 16.46 ppm, whereas contributions from adjacent atoms are negative. The contributions from the adjacent and trans inner hydrogen atoms are −0.79 ppm and −0.57 ppm, respectively. The contributions from N ipso , C α of the same pyrrole ring, the nearest C meso , and the C α atoms of the neighboring pyrrole ring are 0.96, −2.18, −2.38, and −1.85 ppm, respectively. The rest of the atomic contributions are in the range of ±1 ppm. Contributions from each atom to the σ iso of the inner hydrogen are summarized in Figure 8 and the atomic contributions from all atoms to σ iso of the inner hydrogen are reported in the ESI. The paratropic ring current along the outer edge of the molecule in the molecular plane causes deshielding of the meso hydrogen, which has also shielding and deshielding contributions along the inner ring-current pathway at the adjacent pyrrole rings due to the relative direction of the ring current with respect to the studied nucleus. The MSD of the meso hydrogen in Figure 6 has shielding contributions near the hydrogen and the ipso carbon originating from the paratropic ring current inside the hydrogen and from the diatropic current-density flux outside the hydrogen. The MSD contributions at the molecular plane and 1 Å away from it has a similar shape but is weaker as the distance from the molecular plane increases.
Calculations of the atomic contributions to the MSD show that the shielding domain of the meso hydrogen atoms includes the nearest α carbons. Since the distances from the meso hydrogen to other atoms are longer than for the inner hydrogen, the contributions from the adjacent atoms are somewhat smaller than for the α carbons.
The largest atomic contributions to σ iso of the meso hydrogen are 18.80, 12.40, 2.70, and 1.66 ppm from the meso hydrogen, C meso , and from the two closest C α atoms. The largest atomic contributions to σ iso of the β hydrogen are 18.48, 7.06, 2.94, and 1.23 ppm from the studied β hydrogen, C ipso , the adjacent C β and C meso , respectively. Contributions from each atom to the σ iso of the hydrogen atoms are summarized in Figure 8. The atomic contributions to σ iso of the meso and β hydrogen atoms from all atoms are reported in the ESI. The 3D pictures of the MSD of the meso and β hydrogen atoms are reported in the ESI.

Free-Base trans-Norcorrole
The MSD of the inner hydrogen of free-base trans-norcorrole in Figure 9 is dominated by deshielding contributions around the outer edge of the molecule that originate from the strong paratropic ring current [74][75][76]. The ring-current contribution is strong in the molecular plane as well as in a plane 1 Å away from it, where the deshielding contributions completely dominate. Diatropic contributions appear in the molecular plane between the α carbons and the nitrogen atoms inside the pyrrole rings due to the local diatropic currentdensity flux around the nitrogen moiety, which leads to a large shielding domain near the studied inner hydrogen that extends to the plane 1 Å away from the molecular plane. The 3D contour of the MSD of the inner hydrogen in Figure 10 shows a completely dominating deshielding contribution originating from the paratropic ring current. The diatropic ring current along the outer perimeter of free-base trans-norcorrole results in shielding contributions that are larger in the vicinity of the inner hydrogen than on the remote side of the molecule because the magnetic vector potential declines rapidly at longer distances. The significant shielding domain appears near the studied inner hydrogen.
Calculations of MSD contributions from atomic domains show that many atoms contribute to σ iso of the inner hydrogen due to the short distances between it and the atoms inside norcorrole. The atomic contribution of 10.28 ppm from the studied inner hydrogen is shielding, whereas the atomic contributions from atoms in the vicinity are deshielding. The contributions from the C α atoms are in the range of −1.09 to −1.99 ppm. The contribution from the nearest meso carbon is −2.04 ppm. The contributions from the adjacent and trans nitrogen atoms are −2.19 and −1.01 ppm, respectively. The contribution from the ipso nitrogen is only −0.78 ppm. The small negative contributions from many atoms leads to a negative σ iso of −8.73 ppm even though the local contribution is of about the same size but positive. The MSD of the meso hydrogen consists of alternating shielding and deshielding contributions in the molecular plane, as well as in the plane 1 Å away from it. The main shielding areas appear in the vicinity of the studied meso hydrogen originating from the paratropic ring current inside it and the diatropic current density outside the meso hydrogen. The current density around the nitrogen moieties leads to alternating shielding and deshielding contributions, whereas the typical ring-current contributions from other parts of the molecule are not very pronounced. The main atomic contributions to the meso hydrogen are 18.56, 8.93, and 1.06 ppm from the meso hydrogen, C ipso and one of the nearest C α atoms, respectively. The atomic shielding contributions are localized to the nearest atoms because other atoms are very far away from the meso hydrogen, implying that the amplitude of the vector potential of its nuclear magnetic moment vanishes at atoms on the other side of the porphyrinoid ring.
The magnetic shieldings of the β hydrogen atoms are in the range of 26.89 to 28.16 ppm. Since the differences are small, we only discuss one type of β hydrogen atoms. The σ iso value of atom 21, which is one of the β hydrogen atoms of a pyrrole ring without an inner hydrogen, has a local contribution of 18.24 ppm. The largest atomic contributions from adjacent atoms are 3.73, 2.23, 1.25, and 1.16 ppm from C ipso , nearest C α , nearest C β that is not C ipso and the nitrogen of the same pyrrole ring, respectively. Only the nearest atoms contribute to the σ iso of the β hydrogen atoms due to the long distances to many of the other atoms.
Contributions from each atom to the σ iso of the hydrogen atoms are summarized in Figure 11. The atomic contributions to the σ iso values from all atoms are reported in the ESI. The 3D pictures of the MSD of the meso and β hydrogen atoms are reported in the ESI. Figure 11. Atomic contributions (in ppm) to the isotropic 1 H NMR shielding constants of free-base trans-norcorrole. From above: the inner hydrogen, the meso hydrogen, the β hydrogen atoms of the pyrrole ring with an inner hydrogen, and the β hydrogen atoms of the pyrrole ring without an inner hydrogen. The atoms are numbered along the x axis in the same order as in ESI.

Summary and Conclusions
In the present studies on free-base trans-porphyrin, free-base isophlorin, and free-base trans-norcorrole, we have employed a method to calculate and visualize nuclear magnetic shielding densities (MSD) that was recently implemented in the GIMIC program. The MSD of the isotropic 1 H NMR shielding constants has been visualized in the molecular plane and in a plane 1 Å away from it. We also presented 3D pictures of the MSD showing the spatial origin of magnetic shielding and deshielding contributions. We have divided the molecules into atomic domains and determined atomic contributions to the isotropic 1 H NMR shielding constants by numerically integrating each domain separately.
Visualization of the MSD in the two planes shows that the main shielding contributions to the isotropic 1 H NMR shielding constant of the inner hydrogen of free-base transporphyrin originate from the diatropic ring current along the inner pathway, whereas the contributions from the outer one is smaller due to the longer distance between the inner hydrogen and the atoms along the edge of the molecule. The MSD contribution fades out as one goes away from the molecular plane. The current density in the vicinity of the studied inner hydrogen also contributes significantly to its shielding constant. The MSD of the meso hydrogen of free-base trans-porphyrin has shielding contributions from the ring current along the inner pathway and on the outside of the meso hydrogen. The diatropic ring current inside the meso hydrogen deshields its nuclear magnetic moment. Many atomic domains contribute to isotropic 1 H NMR shielding constant of the inner hydrogen, whereas the contributions to the shielding constants of the meso and β hydrogen atoms atoms are dominated by contributions from the local atomic domains.
The MSD of the inner hydrogen of free-base isophlorin has large deshielding contributions from the paratropic ring current and the shielding domain appears mainly in the vicinity of the studied inner hydrogen. The contributions to MSD are stronger in the molecular plane than in the plane 1 Å away from it. The 3D picture of the MSD shows that the contribution from the inner pathway is larger than for the inner one at the remote pyrrole ring, whereas at the adjacent pyrrole rings the deshielding contribution mainly appears along the outer pathway. The meso hydrogen has deshielding contributions from the remote part of the paratropic ring current, whereas its contribution is shielding near the meso hydrogen. The current density on the outside of the meso hydrogen also shields its nuclear magnetic moment. Many atomic domains deshield the nuclear magnetic moment of the inner hydrogen, whereas the local contribution is shielding. The shielding contributions to the shielding constants of the meso and β hydrogen atoms are dominated by contributions from the atomic domains near the studied nucleus, whereas the deshielding contributions are, in comparison, very small.
Free-base trans-norcorrole sustains a paratropic ring current that leads to strong deshielding contributions to the MSD of the inner hydrogen. The deshielding contributions are strong in the molecular plane and in a plane 1 Å away from it. The local diatropic current density around the nitrogen moieties leads to shielding contributions that are largest at the studied inner hydrogen. The 3D picture of the MSD of the inner hydrogen shows that the molecule is dominated by deshielding contributions. Shielding contributions appear at the outer edge of the molecular ring and near the studied inner hydrogen. Integration of atomic domains shows that the largest shielding contribution appear locally, whereas the deshielding contributions from many atomic domains are almost twice as large, yielding a net deshielding contribution to the shielding constant of the inner hydrogen. The MSD of the meso hydrogen is dominated by local shielding contributions, whereas the remote part of the paratropic ring current deshields it. The local diatropic current density around the nitrogen moieties yields alternating shielding and deshielding contributions depending on the relative direction of the current-density flux with respect to the nucleus of the meso hydrogen. The shielding contributions from the atomic domains near the meso and β hydrogen atoms dominate their MSDs, whereas they have very few and small deshielding atomic contributions.
The present study shows the power of the present approach that can be used for determining spatial contributions to 1 H NMR magnetic shieldings. The shielding constants can be assigned to atoms by numerically integrating the MSD in atomic domains. The calculations yield the magnetically induced ring current contributions to the shielding constants of the the inner and outer hydrogen atoms of typical aromatic and anti-aromatic porphyrinoids. Here, we have focused on 1 H NMR shielding constants. However, a similar analysis can also be performed for 13 C NMR and 15 N NMR shielding constants.