Magnetically Induced Ring Currents from DFT-GIMIC Calculations

Aug 20, 2014 - Department of Applied Chemistry, Graduate School of Engineering, Tohoku University, Aoba-yama 6-6-07, Sendai 980-8579, Japan...
0 downloads 0 Views 1MB Size
Subscriber access provided by NATIONAL SUN YAT SEN UNIV

Article

The Antiaromatic Character Of 16#-electron Octaethylporphyrins: Magnetically Induced Ring Currents From DFT-GIMIC Calculations Heike Fliegl, Fabio Pichierri, and Dage Sundholm J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/jp5067549 • Publication Date (Web): 20 Aug 2014 Downloaded from http://pubs.acs.org on August 22, 2014

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry A is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

The Antiaromatic Character Of 16π-electron Octaethylporphyrins: Magnetically Induced Ring Currents From DFT-GIMIC Calculations Heike Fliegl,∗,† Fabio Pichierri,∗,‡ and Dage Sundholm∗,¶ Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, P.O.Box 1033 Blindern, 0315 Oslo, Norway, Department of Applied Chemistry, Graduate School of Engineering, Tohoku University, Aoba-yama 6-6-07, Sendai, 980-8579 Japan , and University of Helsinki, Department of Chemistry, P.O. Box 55 (A.I. Virtanens plats 1), FIN-00014 University of Helsinki, Finland. E-mail: [email protected]; [email protected]; [email protected] Phone: +358 (2)941 50176 . Fax: +358 (2)941 50279



To whom correspondence should be addressed Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, P.O.Box 1033 Blindern, 0315 Oslo, Norway ‡ Department of Applied Chemistry, Graduate School of Engineering, Tohoku University, Aoba-yama 6-6-07, Sendai, 980-8579 Japan ¶ University of Helsinki, Department of Chemistry, P.O. Box 55 (A.I. Virtanens plats 1), FIN-00014 University of Helsinki, Finland. †

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract The magnetically induced current-density susceptibility also called current density has been calculated for a recently synthesized octaethyl porphyrin (OEP) zinc(II) dication with formally 16π electrons. Numerical integration of the current density passing selected chemical bonds yields the current pathway around the porphyrinoid ring and the strength of the ring current. The current strengths show that the OEP-Zn(II) dication is strongly antiaromatic as also concluded experimentally. The calculation of the ring-current pathway shows that all 24 π electrons participate in the transport of the ring current, since the current splits into inner and outer branches of practically equal strengths at the four pyrrolic rings. The corresponding neutral octaethyl porphyrinoid without Zn and inner hydrogens is found to be antiaromatic sustaining a paratropic ring current along the inner pathway with 16 π electrons. The neutral OEP-Zn(II) molecule with formally 18 π electrons is found to be almost as aromatic as free-base porphyrin. However, also in this case all 26 π electrons contribute to the ring current as for free-base porphyrin. A comparison of calculated and measured 1 H NMR chemical shifts is presented. The current strength susceptibility under experimental conditions has been estimated by assuming a linear relation between experimental shielding constants and calculated current strengths.

Keywords Magnetically induced current densities, London orbitals, aromaticity, nuclear magnetic resonance shieldings, porphyrinoids

1

Introduction

Synthetic porphyrin chemistry is a very active research field because the electronic properties of the connected conjugated rings can be adjusted to yield porphyrinoids that may

2 ACS Paragon Plus Environment

Page 2 of 23

Page 3 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

be employed in technological applications. 1–4 A variety of modifications of the basic freebase porphyrin macroring have been devised to obtain desired molecular properties that can be fine-tuned by introducing substituents, forming porphyrin metal complexes, replacing pyrrole rings with other aromatic rings, expanding the porphyrin macrocycle by increasing the number of pyrrole rings, replacing the methine bridge with longer bridging units, compressing the porphyrin by removing methine bridges, etc. The electronic properties of the porphyrinoid macrocycle are to large extent determined by the number of π electrons participating in the conjugation network. 5–8 Experimental and computational studies have shown that the aromatic properties of porphyrinoids follow almost exceptionlessly the generalized H¨ uckel’s rule for aromaticity, 9–22 whereas the dominating conjugation pathway sustaining the magnetically induced ring-current can be reliably obtained only by performing explicit calculations of the current densities. 23–29 Recently, Hiramatsu et al. 30 synthesized a octaethyl porphyrin (OEP) zinc(II) dication without meso substituents that has formally 16π electrons by oxidizing the corresponding 18π porphyrin Zn(II) complex with AgSbF6 and I2 , see Figure 1. Measurements of the 1 H NMR and UV-Vis spectra as well as calculations of nucleus independent chemical shifts (NICS) suggested that the OEP-Zn(II) dication is antiaromatic, which is also expected when assuming that H¨ uckel’s rule for aromaticity holds. The detected a large blue shift of the UV-Vis spectrum and the missing Q band generally indicate that the porphyrinoid is antiaromatic. 31–35 For meso-unsubstituted porphyrinoids, the degree of aromaticity can also be estimated from the position of the 1 H NMR signal of the meso hydrogens. 30 Judging from those experimental data, Hiramatsu et al. claimed that the OEP-Zn(II) dication is the strongest antiaromatic porphyrin that has been reported so far. They also proposed that the reason for the extraordinary strong paratropicity of the porphyrinoid with formally 16π electrons is the absence of the meso substituents allowing a planar structure of the porphyrin macroring.

Yamamoto et al. have reported a nonplanar 16 π Zn porphyrinoid with eight

isobutyl groups at the β carbons and four phenyl groups in the meso positions in addition to

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

an axial Cl ligand coordinated to the metal center. 36 However, the porphyrin ring of axially coordinated 18 π OEP-Zn(II) molecules lacking meso substituents is found to be practically planar. 37 In this work, we computationally study the magnetically induced current density for OEP and OEP-Zn(II) species with 16π and 18π electrons, with the aim to determine the degree of antiaromaticity according to the magnetic criterion as well as the conjugation pathway of the ring-current flow. The article is organized as follows. In Section 2, we present the employed computational methods. The molecular structures of the studied molecules including the nomenclature are described in Section 3. The calculated magnetically induced current densities, current strengths, and current pathways are discussed in Section 4. Calculated and measured 1 H NMR chemical shifts are compared in Section 5. In Section 6, we summarize the obtained results and present our conclusions.

2

Computational methods

The molecular structures were optimized at the density functional theory (DFT) level using Becke’s three-parameter functional (B3LYP) 38,39 as implemented in Turbomole version 6.5. 40 The Karlsruhe triple-ζ quality basis sets augmented with polarization functions (def2TZVP) 41 were employed in the calculations. The initial structures of the different possible conformers were preoptimized with Gaussian 09 at the B3LYP level using double-ζ polarization quality (DGDZVP) basis sets. 42–44 The nuclear magnetic shieldings were calculated at the B3LYP level using split-valence polarization (def2-SVP) basis sets as well as at the B3LYP/def2-TZVP level using the Mpshift module of Turbomole. 45,46 The NMR chemical shifts were obtained from the calculated magnetic shielding constants by using the tetramethylsilane (TMS) value of 31.93 ppm as reference, which has been obtained at the same level of theory. 47 The magnetically induced current densities were calculated at the B3LYP/def2-SVP level

4 ACS Paragon Plus Environment

Page 4 of 23

Page 5 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

using the gimic method. 48–50 Since the use of gauge including atomic orbitals (GIAO) leads to a fast basis set convergence of the current densities, the def2-SVP basis set is found to be sufficient for calculating the current densities. gimic is an independent program that uses the atomic orbital density matrix as well as the corresponding first-order magnetically perturbed density matrices from the nuclear magnetic shielding calculations and basis-set information as input data. 48,49 Ring-current strengths and current pathways were determined by numerical integration of the current-density susceptibility (in nA/T) passing through cut planes perpendicularly to selected bonds of the investigated octaethyl porphyrins. The current-density susceptibility is for simplicity also called current density. The molecular pictures of the current pathways were obtained with jmol 51 and gimp. 52

3

Molecular structures and nomenclature

Each pyrrole ring has two ethyl groups in the β positions. Each of them can be oriented either upwards (u) or downwards (d) with respect to the porphyrin plane. The X-ray structure showed that the ethyl conformation of the neutral OEP-Zn(II) complex is uu-ud-dd-du (c1). 37 We calculated the relative energies of the uu-uu-uu-uu (2c), dd-uu-dd-uu (3c), and du-dudu-du (4c) conformers with respect to 1c. Conformer 2c has all its ethyl groups oriented upwards (or downwards), whereas 3c has every second pair of ethyl groups alternatingly oriented downwards and upwards. For 4c, the ethyl groups of each pair on the pyrrole ring are oriented in opposite directions. Thus, the relative energy of the conformers is 4c < 1c < 3c < 2c. Structure optimizations at the B3LYP/DGDZVP level showed that 4c is the most stable conformer, whereas 2c is highest in energy lying 0.42 kcal/mol above 4c. Conformer 1c, which is found in the molecular crystal, is 0.29 kcal/mol above 4c. Conformer 3c lies 0.36 kcal/mol above 4c. Thus, the energy differences between the conformers are very small implying that small perturbations such as packing effects in the crystal or the employed computational

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

level might shift the energetic order of the conformers. The final molecular structures were optimized at the B3LYP/def2-TZVP level with the preoptimized 4c conformer as initial structure. The obtained structures were employed in the calculations of the NMR chemical shieldings and current-density susceptibilities, which were calculated for the neutral OEP-Zn(II) complex, the OEP-Zn(II) dication, the neutral metal-free octaethyl porphyrinoid, and the metal-free OEP dianion. The neutral OEP-Zn(II) complex and the metal-free OEP dianion have formally 18π electrons, whereas the OEPZn(II) dication and the neutral metal-free octaethyl porphyrinoid without inner hydrogens have formally 16π electrons. Calculations were also performed on the OEP-Zn(II) dication with phenyl substituents in the meso positions to assess how the meso substituents influence the molecular structure and the current strengths. The porphyrin ring of phenyl-substituted OEP-Zn(II)2+ is almost planar with deviations of less than 12 degrees in the dihedral angles of the porphyrin ring. The Cartesian coordinates of the optimized molecular structures are given in the Supporting Information.

4

Current-density calculations

The calculated current pathways and current strengths for the octaethyl porphyrinoids are shown in Figure 2. The main current pathways are indicated with the pink lines. The integrated current strengths passing selected bonds are also given in the figure. Figure 2a shows that the neutral OEP-Zn(II) complex with formally 18π electrons is aromatic sustaining a ring current of 25.2 nA/T around the porphyrin macroring, which is of the same size as the ring-current strength of 27.2 nA/T that has previously been obtained for free-base porphyrin. 19 Thus, the 18π OEP-Zn(II) complex is almost as aromatic as free-base porphyrin. At the pyrrole rings, the current is rather evenly split into a current of 11.9 nA/T along the inner route and 13.3 nA/T passes the β carbons.

6 ACS Paragon Plus Environment

Page 6 of 23

Page 7 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Current density calculations show that the OEP-Zn(II) dication sustain a very strong paratropic ring-current of -51.7 nA/T around the porphyrin macroring. Figure 2b shows that the ring current flow of the OEP-Zn(II) dication splits at the pyrrolic rings into an inner and outer route the current strengths of which are -29.1 nA/T and -21.1 nA/T, respectively. Thus, the current is rather evenly split into the two branches. The current strength along the inner route is 8 nA/T stronger than for the one via the β carbons. The calculated ring-current strength shows that the OEP-Zn(II) dication is strongly antiaromatic as also suggested by Hiramatsu et al. 30 For both neutral OEP-Zn(II) and the corresponding dication, almost no ring current passes through the metal-nitrogen bond, which is expected because radial currents are classically forbidden. The fact that the current flow splits into two almost equally strong branches at the pyrrolic rings shows that all 24 π electrons of the dication participate in the current transport, whereas for the neutral complex, 26 π electrons sustain the ring current. Previous studies on free-base porphyrin and Mg porphyrin showed that there is almost no difference in the current pathways. 53 The same trend is expected to hold for Zn porphyrins. The calculated ring-current pathways show that the proposed 16π-electron pathway for the OEP-Zn(II) dication is not correct. For the neutral OEP-Zn(II) complex, the current does not flow along the classical 18π-electron pathway either as one has traditionally believed and which is still proposed. 29 The role of the Zn(II) dication at the center of the porphyrin is also investigated. The two zinc-free species shown in Figures 2c and 2d are obtained by removing Zn2+ from the OEP-Zn(II) dication and the neutral complex, respectively. The ring-current strength of the neutral metal-free porphyrinoid is -26.5 nA/T indicating that it is also antiaromatic as the parent OEP-Zn(II) dication. Its ring current strength is though only half the current strength of the dication. The current pathway shown in Figure 2c shows that the paratropic ring current of the metal-free porphyrinoid flows along the inner route at the pyrrolic rings. Thus, in this case only 16 π electrons sustain the ring current around the porphyrin macroring.

7 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The OEP dianion is on the other hand aromatic sustaining a ring current of 21.1 nA/T. As seen in Figure 2d, all 26 π electrons of the dianion participate in the current transport, since the current strength along the inner route is 13.6 nA/T and 8.2 nA/T flows via the β carbons. The current pattern of the anion is reminiscent of the average current pathway at the pyrrolic rings with and without inner hydrogens for free-base porphyrin. 19 However, the total current strength is 6 nA/T weaker in the anion than for free-base porphyrin. Substitution of phenyl groups to the meso carbon atom did not significantly change the current pathway around the porphyrinoid macroring for the OEP-Zn(II) dication, see Figure 3. Phenyl-substituted OEP-Zn(II)2+ is also antiaromatic according to the magnetic criterion. The total net current strength of -44.6 nA/T is 13% weaker than for the unsubstituted porphyrinoid. Mulliken population analysis indicates that the perpendicular phenyl groups pull electrons from the porphyrin macroring due to the mesomeric effect (-M) resulting in a somewhat weaker antiaromaticity. The phenyl-substituted OEP-Zn(II) dication is also a very antiaromatic molecule. However, the experimental detection of the degree of antiaromaticity is more problematic for phenyl-substituted porphyrins than for the corresponding molecules with hydrogens at the meso carbons, because of the lacking hydrogens at the meso carbons as 1 H NMR probes. Since the ring currents of the aromatic phenyl rings also shift the 1 H NMR signals of the hydrogens of the neighboring ethyl groups, the degree of aromaticity of the porphyrin ring might be difficult to assess from 1 H NMR chemical shifts.

5

1

H NMR chemical shifts

The calculated 1 H NMR chemical shifts of the meso hydrogens are compared to experimental values in Table 1. The calculated nuclear magnetic shieldings of all atoms are reported in the Supporting Information. The calculated and measured 1 H NMR chemical shifts for neutral OEP-Zn(II) agree well with a discrepancy of only 0.55 ppm. The rather small difference between calculated and measured 1 H NMR chemical shifts is due to vibrational effects,

8 ACS Paragon Plus Environment

Page 8 of 23

Page 9 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

environmental effects and the employed computational level. For the OEP-Zn(II)2+ complex, the deviation of 4.55 ppm between the calculated and measured values of -1.00 ppm and 3.55 ppm,respectively, is much larger. Thus, the comparison of calculated and measured 1 H NMR chemical shifts for the hydrogens at the meso carbons indicates that the calculated ring-current strength is larger than the experimental one. Solvent and counterion effects are omitted in the calculations. The environmental effects might be larger for the dication than for neutral porphyrins, because of the strong Coulomb interactions between the charged species in the solution. The surrounding solvent molecules and counterions most likely decrease the paratropic ring current as the electron-deficient dication is stabilized by accepting a partial charge from them. However, when neglecting environmental and vibrational effects and assigning the difference between calculated and measured 1 H NMR chemical shifts to induced ring currents, the ring-current strength of the porphyrin under experimental conditions can be estimated. Linear interpolation of the calculated ring-current strengths and 1 H NMR chemical shifts in combination with the measured 1

H NMR chemical shifts yield ring-current strengths of 21.6 nA/T for neutral OEP-Zn(II)

and -21.7 nA/T for the dication. The 1 H NMR chemical shifts of the hydrogens at the meso carbons and the ring current strengths are compared in Table 1. The ring currents of the porphyrinoids also affect the proton shifts of the ethyl groups. The calculations show that the average 1 H NMR chemical shifts of the two hydrogens at the carbons connected to the β carbon of the pyrrolic rings are shifted by -4.66 ppm when oxidizing OEP-Zn(II) to OEP-Zn(II)2+ . The methyl hydrogens of the ethyl groups are more distant from the ring current leading to an average ring-current shift of -2.13 ppm. Experimentally the average shift of the 1 H NMR signals of the ethyl groups due to oxidation are -2.67 ppm and -1.31 ppm. By comparing the calculated and measured 1 H NMR chemical shifts of the hydrogens of the ethyl groups, the ring-current strength under experimental conditions can be estimated. An analogous interpolation procedure as used for the hydrogens at the meso positions suggests that the ring-current strength of OEP-Zn(II) is 18.0-19.1

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

nA/T. For the OEP-Zn(II)2+ dication, ring-current strengths of -25.1 nA/T and -29.1 nA/T were estimated from the ethyl hydrogens using the same procedure. The calculated and measured 1 H NMR chemical shifts of the ethyl groups as well as calculated and estimated ring-current strengths are given in Table 2. The uncertainty in the ring-current strength estimated from the ethyl hydrogens is larger than the ones estimated from the hydrogens at the meso carbons, because the ring-current shift is smaller for the ethyl hydrogens than for the hydrogens at the meso carbons. The peripheral hydrogens of the ethyl groups can also be expected to be more exposed to the environment than the hydrogens at the meso carbons. Anyway, one can conclude that the OEP-Zn(II)2+ dication is strongly antiaromatic sustaining a paratropic ring current of about -22 nA/T under the experimental conditions.

6

Summary and conclusions

The magnetically induced current-density susceptibility (current density) has been calculated for a recently synthesized octaethyl porphyrin (OEP) zinc(II) dication with formally 16π. The recent experimental study suggested that the OEP-Zn(II) dication is strongly antiaromatic. 30 In this work, integrations of the current density passing selected bonds show how the ring current flow around the porphyrin macroring. The isolated OEP-Zn(II) dication is found to be very antiaromatic sustaining a paratropic ring current of -51.7 nA/T. Calculated ring-current strengths and the comparison of calculated and measured 1 H NMR chemical shifts indicate that the molecular environment such as the effect of the counter anions and solvent molecules make the molecule less antiaromatic under experimental conditions. However, the computational study confirm that the OEP-Zn(II) dication dissolved in CD2 Cl2 is indeed antiaromatic sustaining a paratropic ring current of about -22 nA/T. The calculations also show that all 24 π electrons participate in the aromatic pathway, because the strengths of the currents via the inner and outer routes at the pyrrolic rings are almost equal. Calculations on the neutral octaethyl porphyrinoid showed that it is also antiaromatic

10 ACS Paragon Plus Environment

Page 10 of 23

Page 11 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

sustaining a paratropic ring current. However, its paratropic ring current is weaker and flows mainly along the inner 16 π route. Calculations on the neutral OEP-Zn(II) complex and the OEP dianion yielded diatropic ring currents showing that the porphyrinoids with formally 18 π electrons are aromatic. Since the ring current takes both the inner and outer routes at the pyrrolic rings, all 26 π electrons contribute to the ring current around the porphyrin macroring. We recently performed a similar study on antiaromatic core-modified isophlorins that were synthesized by Reddy and Anand in 2008. 10,22 The same kind of interpolation procedure as used here to estimate the degree of antiaromaticity under experimental conditions yielded paratropic ring-current strengths of -24.8 nA/T and -30.5 nA/T for the thiophene-furan and furan based isophlorins, respectively. Thus, the strength of the paratropic ring current of the core-modified isophlorins are slightly larger than what we obtained for the OEP-Zn(II) dication dissolved in CD2 Cl2 in this work. The OEP-Zn(II)2+ dication might be the most antiaromatic porphyrin. However, when porphyrinoids like the isophlorins are considered to belong to the class of porphyrins, they seem to be somewhat more antiaromatic than OEP-Zn(II)2+ , since the isophlorins sustain a larger paratropic ring current according to our calculations.

Acknowledgement This research has been supported by the Academy of Finland through its Computational Science Research Programme (LASTU). We thank Magnus Ehrnrooth foundation for financial support. CSC – the Finnish IT Center for Science – is acknowledged for computer time. H. F. thanks for support by the Norwegian Research Council through the CoE Centre for Theoretical and Computational Chemistry (Grant No. 179568/V30). This work has received support from the Norwegian Supercomputing Program (NOTUR) through a grant of computer time (Grant No. NN4654K). F.P. thanks the Department of Applied Chemistry of the Graduate School of Engineering (Tohoku University) for financial support 11 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Supporting Information Available Cartesian coordinates and nuclear magnetic shieldings of the studied molecules. This material is available free of charge via the Internet at http://pubs.acs.org/.

References (1) Mohnani, S.; Bonifazi, D. Supramolecular Architectures Of Porphyrins On Surfaces: The Structural Evolution From 1D To 2D To 3D To Devices. Coord. Chem. Rev. 2010, 254, 2342–2362. (2) Jurow, M.; Schuckman, A. E.; Batteas, J. D.; Drain, C. M. Porphyrins As Molecular Electronic Components Of Functional Devices. Coord. Chem. Rev. 2010, 254, 2297– 2310. (3) Cho, W. J.; Cho, Y.; Min, S. K.; Kim, W. Y.; Kim, K. S. Chromium Porphyrin Arrays As Spintronic Devices. J. Am. Chem. Soc. 2011, 133, 9364–9369. (4) Hasobe, T. Porphyrin-Based Supramolecular Nanoarchitectures For Solar Energy Conversion. J. Phys. Chem. Lett. 2013, 4, 1771–1780. (5) H¨ uckel, E. Grundz¨ uge der Theorie unges¨attigter und aromatischer Verbindungen; Verlag Chemie: Berlin, 1938. (6) von Eggers Doering, W.; Detert, F. L. Cycloheptatrienylium Oxide. J. Am. Chem. Soc. 1951, 73, 876–877. (7) Heilbronner, E. H¨ uckel Molecular Orbitals Of M¨obius-type Conformations Of Annulenes. Tetrahedron Letters 1964, 5, 1923–1928. (8) Rappaport, S. M.; Rzepa, H. S. Intrinsically Chiral Aromaticity. Rules Incorporating Linking Number, Twist, And Writhe For Higher-Twist M¨obius Annulenes. J. Am. Chem. Soc. 2008, 130, 7613–7619. 12 ACS Paragon Plus Environment

Page 12 of 23

Page 13 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(9) Liu, C.; Shen, D.-M.; Chen, Q.-Y. Synthesis And Reactions Of 20 π-Electron βTetrakis(trifluoromethyl)-meso-tetraphenylporphyrins. J. Am. Chem. Soc. 2007, 129, 5814–5815. (10) Reddy, J. S.; Anand, V. G. Planar Meso Pentafluorophenyl Core Modified Isophlorins. J. Am. Chem. Soc. 2008, 130, 3718–3719. (11) Reddy, J. S.; Anand, V. G. Aromatic Expanded Isophlorins: Stable 30π Annulene Analogues With Diverse Structural Features. J. Am. Chem. Soc. 2008, 131, 15433– 15439. (12) Chang, Y.; Chen, H.; Zhou, Z.; Zhang, Y.; Sch¨ utt, C.; Herges, R.; Shen, Z. A 20πElectron Heteroporphyrin Containing A Thienopyrrole Unit. Angew. Chem. Int. Ed. 2012, 51, 12801–12805. (13) Mitsushige, Y.; Yamaguchi, S.; Lee, B. S.; Sung, Y. M.; Kuhri, S.; Schierl, C. A.; Guldi, D. M.; Kim, D.; Matsuo, Y. Synthesis Of Thieno-Bridged Porphyrins: Changing The Antiaromatic Contribution By The Direction Of The Thiophene Ring. J. Am. Chem. Soc. 2012, 134, 16540–16543. (14) Tanaka, Y.; Saito, S.; Mori, S.; Aratani, N.; Shinokubo, H.; Shibata, N.; Higuchi, Y.; Yoon, Z. S.; Kim, K. S.; Noh, S. B. et al. Metalation Of Expanded Porphyrins: A Chemical Trigger Used To Produce Molecular Twisting And M¨obius Aromaticity. Angew. Chem. Int. Ed. 2008, 47, 681–684. (15) Yoon, Z. S.; Osuka, A.; Kim, D. M¨obius Aromaticity And Antiaromaticity In Expanded Porphyrins. Nature Chemistry 2009, 1, 113–122. (16) Shin, J.-Y.; Kim, K. S.; Yoon, M.-C.; Lim, J. M.; Yoon, Z. S.; Osuka, A.; Kim, D. Aromaticity And Photophysical Properties Of Various Topology-controlled Expanded Porphyrins. Chem. Soc. Rev. 2010, 39, 2751–2767.

13 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 23

(17) Higashino, T.; Lim, J. M.; Miura, T.; Saito, S.; Shin, J.-Y.; Kim, D.; Osuka, A. M¨obius Antiaromatic Bisphosphorous Complexes Of [30]Hexaphyrins. Angew. Chem. Int. Ed. 2010, 49, 1–6. (18) Osuka, A.; Saito, S. Expanded Porphyrins And Aromaticity. Chem. Commun. 2011, 47, 4330–4339. (19) H. Fliegl and D. Sundholm, Aromatic Pathways Of Porphins, Chlorins And Bacteriochlorins. J. Org. Chem. 2012, 77, 3408–3414. (20) H.

Fliegl,

and

D.

N.

¨ Ozcan,

Sundholm,

R.

Aromatic

Mera-Adasme, Pathways

F.

In

Pichierri,

Thieno-bridged

J.

Jus´elius

Porphyrins:

Understanding The Influence Of The Direction Of The Thiophene Ring On The Aromatic Character. Mol. Phys. 2013, 111, 1364–1372. (21) R. R. Valiev, H. Fliegl, D. Sundholm, New Insights Into Magnetically Induced Current Pathways And Optical Properties Of Isophlorins. J. Phys. Chem. A 2013, 117, 9062– 9068. (22) R. R. Valiev, H. Fliegl, D. Sundholm, The Aromatic Character Of Thienopyrrolemodified 20π-Electron Porphyrinoids. Phys. Chem. Chem. Phys. 2014, 16, 11010– 11016. (23) Jux, N. The Porphyrin Twist: H¨ uckel And M¨obius Aromaticity. Angew. Chem. Int. Ed. 2008, 47, 2543–2546. (24) Ajami, D.; Hess, K.; K¨ohler, F.; N¨ather, C.; Oeckler, O.; Simon, A.; Yamamoto, C.; Okamoto, Y.; Herges, R. Synthesis And Properties Of The First M¨obius Annulenes. Chem. Eur. J. 2006, 12, 5434–5445. (25) Stepie´ n, M.; Latos-Gra˙zy´ nski, L.; Sprutta, N.; Chwalisz, P.; Szterenberg, L. Expanded

14 ACS Paragon Plus Environment

Page 15 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Porphyrin With A Split Personality: A H¨ uckel-M¨obius Aromaticity Switch. Angew. Chem. Int. Ed. 2007, 46, 7869–7873. (26) Sankar, J.; Mori, S.; Saito, S.; Rath, H.; Suzuki, M.; Inokuma, Y.; Shinokubo, H.; Kim, K. S.; Yoon, Z. S.; Shin, J. Y. et al. Unambiguous Identification Of M¨obius Aromaticity For Meso-Aryl-Substituted [28]Hexaphyrins(1.1.1.1.1.1). J. Am. Chem. Soc. 2008, 130, 13568–13579. (27) Fliegl, H.; Sundholm, D.; Taubert, S.; Pichierri, F. Aromatic Pathways In Twisted Hexaphyrins. J. Phys. Chem. A 2010, 114, 7153–7161. (28) Fliegl, H.; Sundholm, D.; Pichierri, F. Aromatic Pathways In Mono- And Bisphosphorous Singly M¨obius Twisted [28] And [30]Hexaphyrins. Phys. Chem. Chem. Phys. 2011, 13, 20659–20665. (29) Lash, T. D. Origin Of Aromatic Character In Porphyrinoid Systems. J. Porphyrins Phthalocyanines 2011, 15, 1093–1115. (30) S. Hiramatsu, S. Sugawara, S. Kojima and Y. Yamamoto, Synthesis Of The Most Antiaromatic 16π Porphyrin: An Octaethylporphyrin Zinc(II) Complex With No MesoSubstituents. J. Porphyrins Phalocyanines 2013, 17, 1–5. (31) Yamamoto, Y.; Yamamoto, A.; Furuta, S.-y.; Horie, M.; Kodama, M.; Sato, W.; Akiba, K.-Y.; Tsuzuki, S.; Uchimaru, T.; Hashizume, D. et al. Synthesis And Structure Of 16 π Octaalkyltetraphenylporphyrins. J. Am. Chem. Soc. 2005, 127, 14540–14541. (32) Cissell, J. A.; Vaid, T. P.; Yap, G. P. A. The Doubly Oxidized, Antiaromatic Tetraphenylporphyrin Complex [Li(TPP)][BF4]. Org. Lett. 2006, 8, 2401–2404. (33) Sugawara, S.; Hirata, Y.; Kojima, S.; Yamamoto, Y.; Miyazaki, E.; Takimiya, K.; Matsukawa, S.; Hashizume, D.; Mack, J.; Kobayashi, N. et al. Synthesis, Characterization,

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

And Spectroscopic Analysis Of Antiaromatic Benzofused Metalloporphyrins. Chem. Eur. J. 2012, 18, 3566–3581. (34) Kakui, T.; Sugawara, S.; Hirata, Y.; Kojima, S.; Yamamoto, Y. Anti-Aromatic 16π PorphyrinMetal Complexes With Meso-Alkyl Substituents. Chem. Eur. J. 2011, 17, 7768–7771. (35) Sugawara, S.; Kodama, M.; Hirata, Y.; Kojima, S.; Yamamoto, Y. Synthesis And Characterization Of The Most Distorted 16π Porphyrin: 16π Octaisopropyltetraphenylporphyrin (OiPTPP). J. Porphyrins Phthalocyanines 2011, 15, 1326–1334 . (36) Yamamoto, Y.; Hirata, Y.; Kodama, M.; Yamaguchi, T.; Matsukawa, S.; Akiba, K.-y.; Hashizume, D.; Iwasaki, F.; Muranaka, A.; Uchiyama, M. et al. Synthesis, Reactions, And Electronic Properties Of 16 -Electron Octaisobutyltetraphenylporphyrin. J. Am. Chem. Soc. 2010, 132, 12627–12638. (37) A. Ozarowski, H. M. Lee, A. L. Balch, Crystal Environments Probed By EPR Spectroscopy. Variations In The EPR Spectra Of CoII (octaethylporphyrin) Doped In Crystalline Diamagnetic Hosts And A Reassessment Of The Electronic Structure Of FourCoordinate Cobalt(II). J. Am. Chem. Soc. 2003, 125, 12606–12614 . (38) Becke, A. D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. (39) Lee, C.; Yang, W.; Parr, R. G. Development Of The Colle-Salvetti Correlation-energy Formula Into a Functional Of The Electron Density. Phys. Rev. B 1988, 37, 785–789. (40) Ahlrichs, R.; B¨ar, M.; H¨aser, M.; Horn, H.; K¨olmel, C. Electronic Structure Calculations On Workstation Computers: The Program System TURBOMOLE. Chem. Phys. Letters 1989, 162, 165–169, current version: see http://www.turbomole.com.

16 ACS Paragon Plus Environment

Page 16 of 23

Page 17 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(41) Weigend, F.; Ahlrichs, R. Balanced Basis Sets Of Split Valence, Triple Zeta Valence And Quadruple Zeta Valence Quality For H To Rn: Design And Assessment Of Accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297–3305. (42) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A. et al. Gaussian 09, Revision D.01. 2013; Gaussian, Inc., Wallingford CT, see www.gaussian.com. (43) N. Godbout, D. R. Salahub, J. Andzelm, and E. Wimmer, Optimization Of Gaussiantype Basis Sets For Local Spin Density Functional Calculations. Part I. Boron Through Neon, Optimization Technique And Validation. Can. J. Chem. 1992, 70, 560–571. (44) C. Sosa, J. Andzelm, B. C. Elkin, E. Wimmer, K. D. Dobbs, and D. A. Dixon, A Local Density Functional Study Of The Structure And Vibrational Frequencies Of Molecular Transition-Metal Compounds. J. Phys. Chem. 1992, 96, 6630–6636. (45) Sch¨afer, A.; Horn, H.; Ahlrichs, R. Fully Optimized Contracted Gaussian-Basis Sets For Atoms Li To Kr. J. Chem. Phys. 1992, 97, 2571–2577. (46) H¨aser, M.; Ahlrichs, R.; Baron, H. P.; Weis, P.; Horn, H. Direct Computation Of 2nd-Order SCF Properties Of Large Molecules On Workstation Computers With An Application To Large Carbon Clusters. Theoret. Chim. Acta 1992, 83, 455–470. (47) Taubert, S.; Konschin, H.; Sundholm, D. Computational Studies Of 13C NMR Chemical Shifts Of Saccharides. Phys. Chem. Chem. Phys. 2005, 7, 2561–2569. (48) Jus´elius, J.; Sundholm, D.; Gauss, J. Calculation Of Current Densities Using GaugeIncluding Atomic Orbitals. J. Chem. Phys. 2004, 121, 3952–3963. (49) Taubert, S.; Sundholm, D.; Jus´elius, J. Calculation Of Spin-current Densities Using Gauge-including Atomic Orbitals. J. Chem. Phys. 2011, 134, 054123.

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 23

(50) Fliegl, H.; Taubert, S.; Lehtonen, O.; Sundholm, D. The Gauge Including Magnetically Induced Current Method. Phys. Chem. Chem. Phys. 2011, 13, 20500–20518. (51) jmol:

An

Open-source

Java

Viewer

For

Chemical

Structures

In

3D,

http://www.jmol.org. (52) gimp: GNU Image Manipulation Program, http://www.gimp.org . (53) Jus´elius, J.; Sundholm, D. The Aromatic Character Of Magnesium Porphyrins. J. Org. Chem. 2000, 65, 5233–5237 .

18 ACS Paragon Plus Environment

Page 19 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1: Octaethyl Zn(II) porphyrin (Zn(II)OEP) and the Zn(II)OEP2+ dication.

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 23

(a)

(b)

(c)

(d)

Figure 2: Calculated current pathways (pink) and current strengths for selected bonds of (a) neutral octaethyl porphyrin zinc(II) (OEP-Zn(II)), (b) the OEP-Zn(II) dication, (c) the neutral metal-free octaethyl porphyrinoid, and (d) the metal-free OEP dianion.

20 ACS Paragon Plus Environment

Page 21 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 3: Calculated current strengths for selected bonds of the OEP-Zn(II) dication with phenyl substituents at the meso carbon atoms.

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table 1: Calculated chemical shifts (δ in ppm) for the meso hydrogens as compared to experimental data. The difference in the chemical shifts (∆δ in ppm) between the oxidized and reduced forms are also given. The calculated and estimated experimental ring-current strengths (I in nA/T) are also given. Molecule OEP-Zn(II) OEP-Zn(II)2+ OEP2− OEP

1

Calculated H NMR I 10.66 25.2 -1.00 -51.7 9.24 21.3 3.07 -25.2

∆δ – -11.66 – -6.17

1

Experimental H NMR Ia ∆δ(Exp.) 10.11 21.6 – 3.55 -21.7 -6.56 – – – – – –

a

Estimated from calculated ring-current strengths and calculated and measured 1 H NMR chemical shifts for the meso hydrogens. Table 2: The average of the calculated chemical shifts (δ in ppm) for the hydrogens of the ethyl groups of the OEP-Zn(II) molecules as compared to experimental data. The difference in the chemical shifts (∆δ in ppm) between the oxidized and reduced forms are also given. The calculated and estimated experimental ring-current strengths I are given in nA/T. Calculated Experimental 1 1 Molecule H NMR I ∆δ H NMR Ia ∆δ(Exp.) OEP-Zn(II) ethyl 4.34 25.2 – 3.97 19.1 – 2+ OEP-Zn(II) ethyl -0.32 -51.7 -4.66 1.30 -25.0 -2.67 OEP-Zn(II) phenyl 1.98 25.2 – 1.78 18.0 – 2+ OEP-Zn(II) phenyl -0.16 -51.7 -2.13 0.47 -29.1 -1.31 a 1 Estimated from calculated ring-current strengths and calculated and measured H NMR chemical shifts for the corresponding hydrogen of the ethyl group.

22 ACS Paragon Plus Environment

Page 22 of 23

Page 23 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Graphical TOC Entry

Investigation of the change in aromaticity of octaethyl Zn(II) porphyrin upon oxidation.

23 ACS Paragon Plus Environment