Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

One-pot synthesis and property study on thieno[3,2-b]furan compounds

Weimin Maa, Jiawei Huanga, Chao Lia, Yueren Jianga, Baolin Li*a, Ting Qi*a and Xiaozhang Zhu*b
aSchool of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, P. R. China. E-mail: libl@ucas.ac.cn; qiting@ucas.ac.cn
bBeijing National Laboratory for Molecular Sciences, CAS Key Laboratory of Organic Solids, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China. E-mail: xzzhu@iccas.ac.cn

Received 30th January 2019 , Accepted 25th February 2019

First published on 1st March 2019


Abstract

Based on the regioselective intermolecular Suzuki coupling and subsequent intramolecular Ullmann C–O coupling reactions, one-pot synthesis of benzo[4,5]thieno[3,2-b]benzofurans (BTBFs) was developed after optimization of the reaction conditions including catalysts, solvents, bases, ligands and reaction times. The one-pot reaction, with only 2 mol% Pd(PPh3)4 and 2 mol% copper(I) thiophene-2-carboxylate (CuTc) as the catalysts, K3PO4·3H2O as the base and tert-butanol as the solvent, afforded moderate to good yields (up to 70%) for a variety of substrates.


Introduction

Furan-fused ring compounds exhibit high fluorescence quantum efficiency and high carrier mobility, and thus could be widely used as luminescent materials or hole transporting materials in optoelectronic devices.1,2 For instance, benzodifurans (BDFs)3,4 and naphthodifurans (NDFs)5,6 play an important role in the development of high-performance organic light-emitting diodes (OLEDs) and organic field-effect transistors (OFETs).2 Therefore, more and more attention has been paid to the design and synthesis of this type of compounds.7–15

Over the years, some synthetic methods including transition metal catalyzed intramolecular C–O coupling16–18 and bis-C/H-activated ring closure of diaryl ethers19,20 as well as ring closure of ortho-hydroxyarylalkynes21–27 have been reported for the synthesis of furan-fused ring compounds. These methods are usually high-yielding and applicable to a variety of substrates. Nevertheless, they also suffer from some drawbacks, including the usage of expensive palladium catalysts, cumbersome steps for the preparation of the substrates, and so on. Therefore, these methods are not satisfactory from a synthetic chemistry viewpoint, and new cost-effective synthetic methods for furan-fused ring compounds are urgently needed.

The search for alternative safer, cleaner and environmentally friendly synthetic procedures is one of the major activities in synthetic chemistry. To this end, the reduction of wastes together with the use of renewable feedstocks, environmentally friendly reagents and inexpensive catalysts are important pathways to access more sustainable processes.28,29 In this sense, one-pot procedure involving catalytic multistep reactions allows the decrease of energy consuming steps such as separation and purification of intermediates.30 These transformations known as tandem, domino or cascade reactions31,32 have become an important research area in organic chemistry33–35 since they can reduce raw materials and energy consumption and improve atom economy.

As shown in Scheme 1a, thieno[3,2-b]benzofurans have been synthesized by palladium-catalyzed double C–H activation from diaryl ethers.19,20 However, the reactions were performed under rather hash condition (110 or 150 °C), and the precursor diaryl ethers need to be synthesized by a multiple-step reaction. Recently, our group has reported a new strategy for the synthesis of benzo[4,5]thieno[3,2-b]benzofurans (BTBFs) through intramolecular C–O coupling by using cheap copper catalyst (Scheme 1b) under mild condition (90 °C).8 We further noticed that both the Suzuki coupling and the subsequent Ullmann reaction require similar reaction conditions: transition metal catalysts, inorganic bases, and the reaction temperature at around 90 °C. In order to further simplify the synthetic route, we thus design one-pot method30 to synthesize the target product BTBFs. Fortunately, we successfully achieved this goal after screening the reaction conditions. The yields are moderate to good (up to 70%). Meanwhile, this method shows good functional group compatibility, halogenated compounds are also suitable substrates. On the other hand, we studied the photophysical properties of this type of compounds, representing the first example of thieno[3,2-b]furans.


image file: c9ra00796b-s1.tif
Scheme 1 Synthetic strategies toward thieno[3,2-b]benzofurans.

Results and discussion

Initially, all the reactants and catalysts were added at the same time, during the screening of reaction conditions for the synthesis of BTBF (3a). Unfortunately, we didn't detect the target product. Next, we use the method of adding catalysts in two batches without any other treatment of the intermediate product (i.e. Suzuki coupling product). For the first batch, 1a, 2a and the base were mixed in the corresponding solvent and the first catalyst Pd(OAc)2 or Pd(PPh3)4 was added, then the mixture was heated at 90 °C for 24 hours. For the second batch, the second catalyst copper salt was added and heated for another 4 hours. At first, we tried Pd(OAc)2 as the catalyst for Suzuki coupling and CuI as the catalyst for Ullmann reaction, and found that no target product was detected (Table 1, entries 1 and 2). Considering that the most commonly used catalyst for Suzuki coupling is Pd(PPh3)4, we switched to use it as the catalyst in the following condition screening (Table 1, entries 3–24). When the reaction was carried out with 5 mol% CuOAc, CuI, Cu(OAc)2 or copper(I) thiophene-2-carboxylate (CuTc) as the catalyst for Ullmann reaction in N-methyl pyrrolidone (NMP) (Table 1, entries 5–8) at 90 °C, less than 10% yields of desired product were obtained. Nonetheless, we found that CuTc exhibited the highest catalytic activity. Then the effect of different solvents was evaluated (Table 1, entries 8–10), we found that the three solvents including NMP, 1,4-dioxane and tert-butanol afforded essentially the same yields (6.4–6.9%). Considering that tert-butanol is more economical and less toxic than the other two solvents, we chose tert-butanol as the solvent in the following condition screening. For comparison, we noticed that toxic solvents including toluene, N,N-dimethylformamide (DMF) or propionic acid were used in previous reports.8,20
Table 1 Screening of reaction conditionsa

image file: c9ra00796b-u1.tif

Entry Catalyst 1 Catalyst 2 Base 18-crown-6/3 Å M.S.b Solvents GC yield
a Reaction conditions: the mixture of 1a (140 mg, 0.48 mmol), 1b (55 mg, 0.4 mmol), the first catalyst Pd(PPh3)4 and the base in the solvent (5.0 mL) under N2 was stirred at 90 °C for 24 hours. Then the second copper catalyst was added and heated for another 4 hours. 1,1,2,2-Tetrachloroethane was used as an internal standard to determine the GC yield.36b 3 Å molecular sieves were activated under vacuum condition at 250 °C for 5 hours before being used as water absorbent.c Different amount of water was added as additive.d The yield in parentheses is the isolated yield.e 2-Hydroxyphenylboronic acid (2a) was replaced by 2-hydroxyphenylboronic acid pinacol ester.
1 5 mol% Pd(OAc)2 5 mol% CuI K2CO3 —/— 1,4-Dioxane 0
2 5 mol% Pd(OAc)2 5 mol% CuI K2CO3 —/— DMF 0
3 5 mol% Pd(PPh3)4 5 mol% CuI K2CO3 —/— 1,4-Dioxane 1.1%
4 5 mol% Pd(PPh3)4 5 mol% CuI K2CO3 —/— DMF 3.7%
5 5 mol% Pd(PPh3)4 5 mol% CuI K2CO3 —/— NMP 4.6%
6 5 mol% Pd(PPh3)4 5 mol% CuOAc K2CO3 —/— NMP 5.6%
7 5 mol% Pd(PPh3)4 5 mol% Cu(OAc)2 K2CO3 —/— NMP 5.1%
8 5 mol% Pd(PPh3)4 5 mol% CuTc K2CO3 —/— NMP 6.8%
9 5 mol% Pd(PPh3)4 5 mol% CuTc K2CO3 —/— 1,4-Dioxane 6.9%
10 5 mol% Pd(PPh3)4 5 mol% CuTc K2CO3 —/— t-BuOH 6.4%
11 5 mol% Pd(PPh3)4 5 mol% CuTc Cs2CO3 —/— t-BuOH 10%
12 5 mol% Pd(PPh3)4 5 mol% CuTc K3PO4 —/— t-BuOH 29%
13 5 mol% Pd(PPh3)4 5 mol% CuTc K3PO4 —/— t-BuOH/10 μL H2Oc 42%
14 5 mol% Pd(PPh3)4 5 mol% CuTc K3PO4 —/— t-BuOH/0.1 mL H2Oc 35%
15 5 mol% Pd(PPh3)4 5 mol% CuTc K3PO4 —/— t-BuOH/1 mL H2Oc 25%
16 5 mol% Pd(PPh3)4 5 mol% CuTc K3PO4·3H2O —/— t-BuOH 33%
17 5 mol% Pd(PPh3)4 5 mol% CuTc K3PO4·3H2O 5 mol%/— t-BuOH 36%
18 5 mol% Pd(PPh3)4 5 mol% CuTc K3PO4·3H2O —/200 mg t-BuOH 43%
19 5 mol% Pd(PPh3)4 5 mol% CuTc K3PO4·3H2O 5 mol%/200 mg t-BuOH 49%
20 2 mol% Pd(PPh3)4 5 mol% CuTc K3PO4·3H2O 5 mol%/200 mg t-BuOH 56%
21 0.5 mol% Pd(PPh3)4 5 mol% CuTc K3PO4·3H2O 5 mol%/200 mg t-BuOH 46%
22 5 mol% Pd(PPh3)4 2 mol% CuTc K3PO4·3H2O 5 mol%/200 mg t-BuOH 43%
23 2 mol% Pd(PPh3)4 2 mol% CuTc K3PO4·3H2O 5 mol%/200 mg t-BuOH 68% (62%)d
24 2 mol% Pd(PPh3)4 2 mol% CuTc K3PO4·3H2O 5 mol%/200 mg t-BuOH 77% (70%)d,e


We next screened a series of bases (Table 1, entries 10–12) and found that K3PO4 showed the best results (29% yield, Table 1, entry 12). After the screening of solvents and bases, we also explored the effect of the amount of water on reaction yield (Table 1, entries 13–15). To our delight, significant increase of the yield was observed when the reaction was carried out in 5 mL tert-butanol with 10 μL water (Table 1, entry 13). However, addition of excess amount of water reduces the reaction yield (Table 1, entry 15). This could be explained as followed: addition of small amount of water could increase the reaction rate of Suzuki coupling reaction since water and base are required to activate the boronic acid,37,38 and it doesn't do harm to the Ullmann reaction; while addition of large amount of water significantly reduces the catalytic activity of copper catalyst, thus does harm to the Ullmann reaction heavily and the intermediate remained when the reaction mixture was worked up. In order to simplify the experimental operation, we used K3PO4·3H2O as the source of base and water to avoid the need of addition of base and water separately (Table 1, entries 16–24).

On the other hand, we found that addition of 18-crown-6 as an additive to the reaction system can increase the solubility of K3PO4·3H2O37,39 in tert-butanol thus increase the reaction rate and improve the reaction yield (36%, Table 1, entry 17). It is well known that coordinated water in the base favors the previous Suzuki coupling,37 however, water will also coordinate with the copper catalyst in the subsequent Ullmann reaction and reduce its catalytic activity. Therefore, we introduced activated 3 Å molecular sieves as water absorbent into the reaction system after the Suzuki coupling has been carried out, and thus reaction yield has been further improved to 43% (Table 1, entry 18). By further screening the amount of two catalysts (Table 1, entries 19–23), we obtained the highest GC yield of 68% (isolated yield 62%, Table 1, entry 23). It is noted that the amount of two catalysts could be reduced to 2 mol%, which is significantly lower than our previous report [5 mol% Pd(PPh3)4, 10 mol% CuI]8 and the others' reports.19,20 Strikingly, when we used 2-hydroxyphenyl-boronic acid pinacol ester instead of 2-hydroxybenzeneboronic acid (Table 1, entry 24), GC yield can be further increased to 77% (isolated yield 70%), which is an improvement over the previously reported yield (67% total yield via two steps).

With the optimal reaction condition in hand, we next explored reaction scope. The reaction proved to be robust and versatile and allowed the synthesis of a variety of BTBFs as shown in Table 2. It worked well for various substrates bearing either electron-withdrawing or electron-donating groups on the aromatic rings, such as F, Cl, Me, C8H17, or t-Bu, affording the desired products in moderate to good yields (Table 2, 3a–3i). The yields are not high, it is understandable considering that both the Suzuki coupling and the subsequent Ullmann reaction are usually not very high yielding. To our delight, for the known compounds (3a–3b), we achieved higher yield than previous reports (70% vs. 67%, 61% for 3a; 60% vs. 46% for 3b).8,20 Furthermore, the π-extended product BTNF (3j) could be obtained through this approach in moderate yield (35%). This can extend the conjugation of this type of compounds and change their photophysical properties.

Table 2 Substrates scopea

image file: c9ra00796b-u2.tif

a Reactions were conducted on 1 (0.48 mmol), 2 (0.4 mmol) and K3PO4·3H2O (319 mg, 1.2 mmol), and all the yields are isolated ones.b Total yields via two steps in our previous report.8c Isolated yield reported by Kanai group.20
image file: c9ra00796b-u3.tif


The UV-Vis absorption and fluorescence spectra of BTBF (3a), F-BTBF (3c), F,Cl-BTBF (3d), M-BTBF (3e), and BTNF (3j) in dilute ethanol solutions were investigated. As shown in Fig. 1, the absorption spectra of these compounds exhibits fine structures, and the maximum absorption peaks at 318, 320, 319, 319 and 355 nm respectively could be assigned to π–π* transitions of the conjugated backbones containing sulfur and oxygen heteroatoms.40,41 We can see that either halogen atoms on the F-BTBF and F,Cl-BTBF molecules or methyl group in the M-BTBF molecule have neglectable effect on the maximum absorption. But for BTNF, the introduction of naphthalene ring extends the conjugated system when compared to BTBFs, leading to a significant red shift of absorption peak.


image file: c9ra00796b-f1.tif
Fig. 1 UV-Vis absorption spectra in ethanol (6 × 10−5 mol L−1).

As shown in Fig. 2, the maximum emission peaks of these compounds are located at 339, 342, 347, 340, 378 nm, respectively. The fluorescence quantum yields (Φf) of the compounds were determined to be 0.01, 0.03, 0.04, 0.10, 0.72 respectively in dilute ethanol solutions, with 9,10-diphenylanthracene (Φf = 0.95)42 as the reference standard. To be noted, BTNF exhibits significantly red-shifted fluorescence peak at 378 nm and much increased fluorescence quantum yield due to its extended π system provided by the rigid naphthalene ring. The introduction of rigid naphthalene ring may reduce the nonradiative decay rate, resulting in high fluorescence quantum yield. Thus, it is expected to be a good luminescent material. To our best knowledge, this is the first example of study on photophysical properties of thieno[3,2-b]furans.


image file: c9ra00796b-f2.tif
Fig. 2 Fluorescence spectra in ethanol (2 × 10−7 mol L−1).

To understand the frontier molecular orbitals and spectroscopic properties of these five compounds, the computation of molecular orbital geometries were performed by density functional theory (DFT) at the B3LYP/6-31G level of theory.43,44 In Fig. S2, we can see that the highest occupied molecular orbitals (HOMOs) and lowest occupied molecular orbitals (LUMOs) of these compounds are homogeneously delocalized over the entire backbone. The HOMO levels (EHOMOs) of BTBF, F-BTBF, F,Cl-BTBF, M-BTBF and BTNF are calculated to be −5.71, −5.94, −6.22, −5.66, −5.45 eV, respectively, and the LUMO levels (ELUMOs) are −1.29, −1.54, −1.82, −1.24, −1.64 eV, respectively (Fig. S3).

Compared with BTBF, F-BTBF and F,Cl-BTBF show much deeper HOMO and LUMO energy levels resulting from the existence of electron-withdrawing halogen atoms (Fig. S3). On the contrary, the introduction of an electron-donating methyl group slightly raises the HOMO and LUMO energy levels of M-BTBF. Obviously, the tendency of the calculated energy gaps of these compounds are in agreement with those of the energy gaps calculated from the UV absorption onset (Table 3 and Fig. S3).

Table 3 Photophysical properties of BTBFs and BTNF
Compounds λabs,max (nm) Ega (ev) λFL (nm) Stoke shift (nm) Φfb
a The optical band gap Eg is calculated by the formula Eg = 1240/λonset. λonset represents the largest edge of the UV-Visible absorption spectrum.b 9,10-Diphenylanthracene (Φf = 0.95) as a reference in ethanol.42
BTBF (3a) 318 3.79 339 21 0.01
F-BTBF (3c) 320 3.73 342 22 0.03
F,Cl-BTBF (3d) 319 3.75 347 28 0.04
M-BTBF (3e) 319 3.77 340 21 0.10
BTNF (3j) 355 3.41 378 23 0.72


Conclusions

In summary, we have developed a one-pot method for the construction of ladder-type BTBFs. One-pot method can not only eliminate the requirement of separation of the intermediate product, but also reduce the use of various solvents, in line with the concept of green chemistry. In addition, low loading of palladium and copper catalysts and mild reaction temperature make this strategy more economic and attractive. On the other hand, photophysical property study showed that π-extended BTNF exhibits high quantum efficiency. We believe the simple synthesis of BTBFs in one-pot method makes it more convenient to prepare these compounds and will promote their application as semiconductor materials.

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

We thank the National Natural Science Foundation of China (No. 21402194), CAS Center of Excellence in Molecular Science, Chinese Academy of Sciences and Beijing National Laboratory for Molecular Sciences (BNLMS20140158) for financial support, and the foundation supported by the University of Chinese Academy of Sciences.

Notes and references

  1. B. Li, Chin. J. Org. Chem., 2015, 35, 2487–2506 CrossRef CAS .
  2. H. Tsuji and E. Nakamura, Acc. Chem. Res., 2017, 50, 396–406 CrossRef CAS PubMed .
  3. H. Tsuji, C. Mitsui, L. Ilies, Y. Sato and E. Nakamura, J. Am. Chem. Soc., 2007, 129, 11902–11903 CrossRef CAS PubMed .
  4. H. Tsuji, C. Mitsui, Y. Sato and E. Nakamura, Adv. Mater., 2009, 21, 3776–3779 CrossRef CAS .
  5. C. Mitsui, J. Soeda, K. Miwa, H. Tsuji, J. Takeya and E. Nakamura, J. Am. Chem. Soc., 2012, 134, 5448–5451 CrossRef CAS PubMed .
  6. K. Niimi, H. Mori, E. Miyazaki, I. Osaka, H. Kakizoe, K. Takimiya and C. Adachi, Chem. Commun., 2012, 48, 5892–5894 RSC .
  7. S. Agasti, A. Dey and D. Maiti, Chem. Commun., 2017, 53, 6544–6556 RSC .
  8. D. Chen, D. Yuan, C. Zhang, H. Wu, J. Zhang, B. Li and X. Zhu, J. Org. Chem., 2017, 82, 10920–10927 CrossRef CAS PubMed .
  9. M. Heravi, V. Zadsirjan, H. Hamidi, A. Tabar and H. Parvin, RSC Adv., 2017, 7, 24470–24521 RSC .
  10. A. Abu-Hashem, H. Hussein, A. Aly and M. Gouda, Synth. Commun., 2014, 44, 2285–2312 CrossRef CAS .
  11. H. Miyabe, Molecules, 2015, 20, 12558–12575 CrossRef CAS PubMed .
  12. X. Duan, J. Feng and Z. Zhang, Synthesis, 2010, 515–519 CrossRef CAS .
  13. M. Heravi and V. Zadsirjan, Adv. Heterocycl. Chem., 2015, 117, 261–376 CrossRef CAS .
  14. D. Chen, J. Li, W. Ma, B. Li, Y. Zhen, X. Zhu, W. Hu, H. Tsuji and E. Nakamura, Asian J. Org. Chem., 2018, 7, 2228–2232 CrossRef CAS .
  15. E. Nakamura, H. Tsuji and L. Ilies, Synlett, 2014, 25, 2099–2110 CrossRef .
  16. K. Kawaguchi, K. Nakano and K. Nozaki, J. Org. Chem., 2007, 72, 5119–5228 CrossRef CAS PubMed .
  17. S. Newman, V. Aureggi, C. Bryan and M. Lautens, Chem. Commun., 2009, 5236–5238 RSC .
  18. W. Chen, Y. Zhang, L. Zhang, M. Wang and L. Wang, Chem. Commun., 2011, 47, 10476–10478 RSC .
  19. H. Kaida, T. Satoh, K. Hirano and M. Miura, Chem. Lett., 2015, 44, 1125–1127 CrossRef CAS .
  20. K. Saito, P. K. Chikkade, M. Kanai and Y. Kuninobu, Chem.–Eur. J., 2015, 21, 8365–8368 CrossRef CAS PubMed .
  21. K. Nakanishi, D. Fukatsu, K. Takaishi, T. Tsuji, K. Uenaka, K. Kuramochi, T. Kawabata and K. Tsubaki, J. Am. Chem. Soc., 2014, 136, 7101–7109 CrossRef CAS PubMed .
  22. Y. Liang, S. Tang, X. Zhang, L. Mao, Y. Xie and J. Li, Org. Lett., 2006, 8, 3017–3020 CrossRef CAS PubMed .
  23. A. Varela-Fernandez, C. Gonzalez-Rodrıguez, J. Varela, L. Castedo and C. Saa, Org. Lett., 2009, 11, 5350–5353 CrossRef CAS PubMed .
  24. A. Arcadi, S. Cacchi, M. DelRosario, G. Fabrizi and F. Marinelli, J. Org. Chem., 1996, 61, 9280–9288 CrossRef CAS .
  25. A. Furstner and P. Davies, J. Am. Chem. Soc., 2005, 127, 15024–15025 CrossRef PubMed .
  26. K. Anderson, T. Ikawa, R. Tundel and S. Buchwald, J. Am. Chem. Soc., 2006, 128, 10694–10695 CrossRef CAS PubMed .
  27. R. N. Nair, P. J. Lee, A. L. Rheingold and D. B. Grotjahn, Chem.–Eur. J., 2010, 16, 7992–7995 CrossRef CAS PubMed .
  28. I. T. Horvath and P. T. Anastas, Chem. Rev., 2007, 107, 2167–2168 CrossRef CAS PubMed .
  29. I. T. Horvath and P. T. Anastas, Chem. Rev., 2007, 107, 2169–2173 CrossRef CAS PubMed .
  30. A. C. Maria, J. Climent and S. Iborra, Chem. Rev., 2011, 111, 1072–1133 CrossRef PubMed .
  31. D. E. Fogg and E. N. Santos, Coord. Chem. Rev., 2004, 248, 2365–2379 CrossRef CAS .
  32. G. Poli and G. Giambastiani, J. Org. Chem., 2002, 67, 9456–9459 CrossRef CAS PubMed .
  33. L. F. Tietze, Chem. Rev., 1996, 96, 115–136 CrossRef CAS PubMed .
  34. J. M. Lee, Y. Na, H. Han and S. Chang, Chem. Soc. Rev., 2004, 33, 302–312 RSC .
  35. J. C. Wasilke, S. J. Obrey, R. T. Baker and G. C. Bazan, Chem. Rev., 2005, 105, 1001–1020 CrossRef CAS PubMed .
  36. F. Mo and G. Dong, Science, 2014, 345, 68–72 CrossRef CAS PubMed .
  37. N. Miyaura and A. Suzuki, Chem. Rev., 1995, 95, 2457–2483 CrossRef CAS .
  38. G. B. Smith, G. C. Dezeny, D. L. Hughes, A. O. King and T. R. Verhoeven, J. Org. Chem., 1994, 59, 8151–8156 CrossRef CAS .
  39. K. Rubina, E. Abele, P. Arsenyan, M. Fleisher, J. Popelis, A. Gaukhman and E. Lukevics, Tetrahedron, 2003, 59, 7603–7607 CrossRef CAS .
  40. Y. Zhang, F. Peng, R. He, L. Ying, W. Yang and Y. Cao, New J. Chem., 2018, 42, 2750–2757 RSC .
  41. Y. Fogel, M. Kastler, Z. Wang, D. Andrienko, G. Bodwell and K. Mullen, J. Am. Chem. Soc., 2007, 129, 11743–11749 CrossRef CAS PubMed .
  42. J. V. Morris, M. A. Mahaney and J. R. Huber, J. Phys. Chem., 1976, 80, 969–974 CrossRef CAS .
  43. C. Lee, W. Yang and R. G. Parr, Phys. Rev., 1988, 37, 785–789 CAS .
  44. M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci, G. A. Petersson, H. Nakatsuji, M. Caricato, H. X. Li, P. Hratchian, A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Hada, M. Ehara and K. Toyota, Gaussian 09, revision A.01, Gaussian, Inc., Wallingford, CT, 2009 Search PubMed .

Footnote

Electronic supplementary information (ESI) available. See DOI: 10.1039/c9ra00796b

This journal is © The Royal Society of Chemistry 2019