Sky-blue Thermally Activated Delayed Fluorescence Material Employing Diphenylethyne Acceptor for Organic Light-Emitting Diode Application

Similar documents
Motooka, Nishi, Fukuoka , Japan.

Department of Chemistry, School of life Science and Technology, Jinan University, Guangzhou , China b

Electronic Supplementary Information (ESI) for Chem. Commun.

Group 13 BN dehydrocoupling reagents, similar to transition metal catalysts but with unique reactivity. Part A: NMR Studies

Electronic Supplementary Information. Electron Mobility for All-Polymer Solar Cells

Supporting Information

A Computational Model for the Dimerization of Allene: Supporting Information

Supporting Information. 4-Pyridylnitrene and 2-pyrazinylcarbene

Scalable synthesis of quaterrylene: solution-phase

Supporting Information

Supporting Information For. metal-free methods for preparation of 2-acylbenzothiazoles and. dialkyl benzothiazole-2-yl phosphonates

Supporting Information. {RuNO} 6 vs. Co-Ligand Oxidation: Two Non-Innocent Groups in One Ruthenium Nitrosyl Complex

Supporting Information. O-Acetyl Side-chains in Saccharides: NMR J-Couplings and Statistical Models for Acetate Ester Conformational Analysis

The Activation of Carboxylic Acids via Self Assembly Asymmetric Organocatalysis: A Combined Experimental and Computational Investigation

Non-Radiative Decay Paths in Rhodamines: New. Theoretical Insights

Experimental Evidence for Non-Canonical Thymine Cation Radicals in the Gas Phase

Supporting Information

Supporting Information

CICECO, Departamento de Química, Universidade de Aveiro, Aveiro, Portugal

Two-Dimensional Carbon Compounds Derived from Graphyne with Chemical Properties Superior to Those of Graphene

SUPPORTING INFORMATION

Detailed Syntheses. K 5H[Co II W 12O 40] 15H 2O (1). The synthesis was adapted from published methods. [1] Sodium tungstate

Supporting Information. Detection of Leucine Aminopeptidase Activity in Serum Using. Surface-Enhanced Raman Spectroscopy

Supporting information

Driving forces for the self-assembly of graphene oxide on organic monolayers

Supplementary Information

Phosphine Oxide Jointed Electron Transporters for Reducing Interfacial

Highly sensitive detection of low-level water contents in. organic solvents and cyanide in aqueous media using novel. solvatochromic AIEE fluorophores

SUPPORTING INFORMATION. Ammonia-Borane Dehydrogenation Promoted by a Pincer-Square- Planar Rhodium(I)-Monohydride: A Stepwise Hydrogen Transfer

Supplementary Information

Supplementary Information:

Supporting Material Biomimetic zinc chlorin poly(4-vinylpyridine) assemblies: doping level dependent emission-absorption regimes

SUPPLEMENTARY INFORMATION

DFT and TDDFT calculation of lead chalcogenide clusters up to (PbX) 32

BINOPtimal: A Web Tool for Optimal Chiral Phosphoric Acid Catalyst Selection

Electronic relaxation dynamics of PCDA PDA studied by transient absorption spectroscopy

Cationic Polycyclization of Ynamides: Building up Molecular Complexity

Supporting Information

Observation of Guanidine-Carbon Dioxide Complexation in Solution and Its Role in Reaction of Carbon Dioxide and Propargylamines

SUPPLEMENTARY MATERIAL. Nitrogen Containing Ionic Liquids: Biodegradation Studies and

Supporting Information

Supporting Information for. Acceptor-Type Hydroxide Graphite Intercalation Compounds. Electrochemically Formed in High Ionic Strength Solutions

Ni II Tetrahydronorcorroles: Antiaromatic Porphyrinoids with Saturated Pyrrole Units

Supplementary Material: 2D-IR Spectroscopy of an AHA Labelled Photoswitchable PDZ2 Domain

SUPPORTING INFORMATION

Synthesis, electrochemical and theoretical studies of. V-Shaped donor-acceptor hexaazatriphenylene. (HAT) derivatives for second harmonic generation

Electronic Supporting Information

Elucidating the structure of light absorbing styrene. carbocation species formed within zeolites SUPPORTING INFORMATION

Ring expansion reactions of electron-rich boron-containing heterocycles. Supporting Information Contents

Dinuclear Nickel Complexes in Five States of Oxidation Using a Redox-Active Ligand. Supporting Information

Supporting Information

Metal-Free Hydrogenation Catalysis of Polycyclic Aromatic Hydrocarbons

A dominant homolytic O-Cl bond cleavage with low-spin triplet-state Fe(IV)=O formed is revealed in the mechanism of heme-dependent chlorite dismutase

Supplementary information

3,4-Ethylenedioxythiophene (EDOT) and 3,4- Ethylenedioxyselenophene (EDOS): Synthesis and Reactivity of

Organic photodiodes constructed from a single radial heterojunction. microwire

Electronic Supplementary Information (ESI)

Electronic Supplementary Information

Planar Pentacoordinate Carbon in CAl 5 + : A Global Minimum

Nucleophilicity Evaluation for Primary and Secondary Amines

Supporting Information. Synthesis, resolution and absolute configuration of chiral 4,4 -bipyridines

Table S1: DFT computed energies of high spin (HS) and broken symmetry (BS) state, <S 2 > values for different radical systems. HS BS1 BS2 BS3 < S 2 >

A Redox-Fluorescent Molecular Switch Based on a. Heterobimetallic Ir(III) Complex with a Ferrocenyl. Azaheterocycle as Ancillary Ligand.

Mechanism of Hydrogen Evolution in Cu(bztpen)-Catalysed Water Reduction: A DFT Study

Decomposition!of!Malonic!Anhydrides. Charles L. Perrin,* Agnes Flach, and Marlon N. Manalo SUPPORTING INFORMATION

Out-Basicity of 1,8-bis(dimethylamino)naphthalene: The experimental and theoretical challenge

Supporting Information

ELECTRONIC SUPPLEMENTARY INFORMATION

Supporting Information. A Polymer Acceptor with An Optimal LUMO Energy Level. for All-Polymer Solar Cells

Supporting Information

A Key n π* Interaction in N-Acyl Homoserine Lactones

Supporting Information. for. Silylation of Iron-Bound Carbon Monoxide. Affords a Terminal Fe Carbyne

A Fluorescent Heteroditopic Hemicryptophane Cage for the Selective Recognition of Choline Phosphate

Supporting Information

Ethynylene-linked Figure-eight. Octaphyrin( ): Synthesis and. Characterization of Its Two Oxidation States

High-Efficiency Organic Light-Emitting Diodes Based on Sublimable Cationic Iridium(III) Complexes with Sterically Hindered Spacers

INTERNATIONAL JOURNAL OF RESEARCH IN COMPUTER APPLICATIONS AND ROBOTICS ISSN

A mechanistic study supports a two-step mechanism for peptide bond formation on the ribosome

Ionic liquid electrolytes for reversible magnesium electrochemistry

Methionine Ligand selectively promotes monofunctional adducts between Trans-EE platinum anticancer drug and Guanine DNA base

Electronic Supplementary Information

Speciation of Adsorbed Phosphate at Gold Electrodes: A Combined. Surface-Enhanced Infrared Absorption Spectroscopy and DFT Study

1,4-Benzene-Bridged Covalent Hybrid of Triarylamine and Cyclometalated Ruthenium: A New Type of Organic-Inorganic Mixed-Valent System

Supporting Information. Reactive Simulations-based Model for the Chemistry behind Condensed Phase Ignition in RDX crystals from Hot Spots

Supplementary Information. Institut für Anorganische Chemie, Julius-Maximilians Universität Würzburg, Am Hubland, D Würzburg, Germany.

Unveiling the role of hot charge-transfer states. in molecular aggregates via nonadiabatic. dynamics

Electronic Supplementary Information for:

Supporting Information

Spin-Orbit Coupling, Spin-Spin, NMR and Electric Properties of Hydrogen Chloride Molecule

Electronic Supplementary Information

Supporting Information

Supporting Information

Synthesis and Electronic Properties of Length-Defined 9,10-Anthrylene Butadiynylene Oligomers

Truong Ba Tai, Long Van Duong, Hung Tan Pham, Dang Thi Tuyet Mai and Minh Tho Nguyen*

Supporting Information for. Chemoselective Asymmetric Dearomative [3+2] Cycloaddition Reactions of Purines with Aminocyclopropanes

Supporting Information. Copyright Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim, 2008

Supporting Information Computational Part

Supplementary Figure 1. Energy diagram of constitutional isomers 1a, 5a and 2a.

Interfacial Hydration Dynamics in Cationic Micelles Using 2D-IR and NMR

Supporting Information

Transcription:

Electronic Supplementary Material (ESI) for Journal of Materials Chemistry C. This journal is The Royal Society of Chemistry 2017 Supporting Information Sky-blue Thermally Activated Delayed Fluorescence Material Employing Diphenylethyne Acceptor for Organic Light-Emitting Diode Application Zuozheng He, Xinyi Cai, Zhiheng Wang, Yunchuan Li, Zhida Xu, Kunkun Liu, Dongcheng Chen, and Shi-Jian Su* State Key Laboratory of Luminescent Materials and Devices and Institute of Polymer Optoelectronic Materials and Devices, South China University of Technology, Guangzhou 510640, China The authors contribute equally to this work. *Corresponding author. E-mail: mssjsu@scut.edu.cn

Synthesis. All solvents and reagents were used as received from commercial suppliers without any further purification. Synthetic routes of the compounds are shown in Scheme S1, and they were further purified by repeated temperature gradient vacuum sublimation. Br I Si(CH 3 ) 3 PdCl 2 (PPh 3 ) 2, CuI DBU, H 2 O Br Br Br Br Pd(OAc) 2, (t-bu) 3 P H or H O O O Scheme S1. Synthetic routes of the investigated compounds. Synthesis of 1,2-bis(4-bromophenyl)ethyne (DPE-2Br): A solution of 1-bromo-4-iodobenzene (3 g, 10 mmol) was dissolved in toluene (50 ml) in a three necked 100 ml flask. Bis(triphenylphosphine)palladium(II) chloride (PdCl 2 (PPh 3 ) 2, 421 mg, 6 mol%), CuI (190 mg, 10 mol%), 1,8-diazabicyclo[5.4.0]undec-7-ene (DBU) (8.8 ml, 60 mmol) and (trimethylsilyl)acetylene (0.69 ml, 5 mmol) were added under 2 atmosphere, followed by the addition of H 2 O (72 μl, 0.4 eq.). The reaction mixture was stirred at room temperature in the dark overnight. Upon completion of the reaction, HCl was poured into the crude mixture and then it was extracted with chloroform. The organic layer was dried over anhydrous a 2 SO 4, and the solvent removed under reduced pressure. The crude product was purified by silica gel column chromatography using petroleum ether (PE) as eluent, to give a yellow powder in 60% yield (2 g). 1 H MR (500 MHz, chloroform-d) δ 7.51 7.46 (m, 4H), 7.40 7.35 (m, 4H). Synthesis of 1,2-bis(4-(9,9-dimethylacridin-10(9H)-yl)phenyl)ethyne (DPE-DDMAc): (1,2-Bis(4-bromophenyl)ethyne (1 g, 3 mmol ), 9,9-dimethyl-9,10-dihydroacridine (1.4 g, 6.6 mmol), and K 2 CO 3 (2.4 g, 18 mmol) were added into a 100 ml three neck flask with 60 ml toluene in 2 atmosphere. After degassing for 15 min, acetic acid palladium (II) (67.2 mg, 0.3 mmol) and tritert-butylphosphine (1 ml, 1 mmol) were added. Subsequently, the mixture was stirred and refluxed overnight, and then cooled to room temperature. After removing the solvent in vacuum, the mixture was partitioned between DCM and water. The combined organic layers were washed with brine, dried over Mg 2 SO 4 and purified by column chromatography on silica gel (eluent: DCM/PE=3/1), affording a white solid in 80% yield (1.4 g). 1 H MR (500 MHz, chloroform-d) δ 7.87 7.79 (m, 4H), 7.47 (dd, J = 7.6, 1.7 Hz, 4H), 7.40 7.33 (m, 4H), 7.04 6.89 (m, 8H), 6.31 (dd, J = 8.2, 1.3 Hz, 4H), 1.70 (s, 12H). 13 C MR (126 MHz, chloroform-d) δ 141.44, 140.68, 140.66, 134.19, 131.77, 131.74, 131.56, 131.26, 131.19, 130.15, 126.45, 126.41, 125.30, 125.21, 125.14, 123.02, 120.82, 120.78, 120.69, 114.02, 113.89, 89.74, 36.00, 31.22, 31.02, 30.99. MALDI-TOF MS (mass m/z): calc 592.2878, found 592.2883 [M+]. Synthesis of 1,2-bis(4-(10H-phenoxazin-10-yl)phenyl)ethyne (DPE-DPXZ): (1,2-Bis(4-bromophenyl)ethyne (1 g, 3 mmol ), 10H-phenoxazine (1.2 g, 6.6 mmol), and K 2 CO 3 (2.4 g, 18 mmol) were added into a 100 ml three neck flask with 60 ml toluene in 2 atmosphere. After

degassing for 15 min, acetic acid palladium (II) (67.2 mg, 0.3 mmol) and tri-tert-butylphosphine (1 ml, 1 mmol) were added. Subsequently, the mixture was stirred and refluxed overnight, and then cooled to room temperature. After removing the solvent in vacuum, the mixture was partitioned between DCM and water. The combined organic layers were washed with brine, dried over Mg 2 SO 4 and purified by column chromatography on silica gel (eluent: DCM/PE=3/1), affording a white solid in 82% yield (1.3 g). 1 H MR (500 MHz, chloroform-d) δ 7.83 7.73 (m, 4H), 7.42 7.33 (m, 4H), 6.82 6.46 (m, 12H), 6.02 5.90 (m, 4H). 13 C MR (126 MHz, chloroform-d) δ 143.96, 139.15, 134.34, 133.99, 131.10, 123.27, 123.22, 121.62, 115.55, 113.24, 89.69. MALDI-TOF MS (mass m/z): calc 540.1871, found 540.1843 [M+]. Chemical Characterization. Fig. S1. 1 H MR spectrum of compound DPE-2Br in deuterated CDCl 3.

Fig. S2. 1 H MR spectrum of compound DPE-DDMAc in deuterated CDCl 3. Fig. S3. 13 C MR spectrum of compound DPE-DDMAc in deuterated CDCl 3.

Fig. S4. 1 H MR spectrum of compound DPE-DPXZ in deuterated CDCl 3. Fig. S5. 13 C MR spectrum of compound DPE-DPXZ in deuterated CDCl 3.

Fig. S6. MALDI-TOF mass spectrum of compound DPE-DDMAc. Fig. S7. MALDI-TOF mass spectrum of compound DPE-DPXZ.

Fig. S8. Schematic diagram of sp hybrid orbitals. Computation. Computational details: All of the simulations were performed using the Gaussian 09_B01 program package. 1-3 Detailed calculation approach affording E 0-0 ( 1 CT), E 0-0 ( 3 CT) and E 0-0 ( 3 LE) are Equations S1 to S3. E 0 0 ( 1 CT) = E VA (S 1, OHF) 0.24 E 0 0 ( 3 CT) = E 0 0 (S 1 ) [ E VA(S 1, OHF) E VA (S 1, OHF) E VA (S 1, B3LYP) E VA(T 1, B3LYP)] E 0 0 ( 3 LE) = E VA (T 1 ) [ E VA(S 1, OHF) E VA (S 1, B3LYP)] 0.09 (S3) (S1) (S2) Vertical Absorption Energy (ev) 4 DPE-DPXZ (E VA (S 1 )) DPE-DPXZ (E VA (T 1 )) 3 2 DPE-DDMAc (E VA (S 1 )) DPE-DDMAc (E VA (T 1 )) DMAc (E VA (T 1 )) PXZ (E VA (T 1 )) DPE (E VA (T 1 )) 10 100 OHF (%) Fig. S9. Vertical absorption energy depending on HF% curves of the final compounds DPE-DPXZ and DPE-DDMAc as well as those of the intermediates (DMAc, PXZ and DPE) are shown for comparison. Fig. S10. Front views of horizontal oriented backbone (-1807.84419541 a.u.) and vertical oriented backbone (-1807.84288910 a.u.) DPE-DDMAc optimized in gas phase in B3LYP/6-31G* level. We tentatively took the horizontal oriented backbone one as the initial guess for the subsequent

simulation for its ground state energy is lower than the vertical oriented backbone one. Fig. S11. Electrostatic potential (ESP) distribution of two conformers of DPE-DDMAc. The positive electrostatic sites mainly localize on the DPE moiety for both molecules which induced moderate electron-withdrawing ability of the DPE unit. LUMO -1.71 LUMO+1-4.94 HOMO HOMO-1 Fig. S12. Frontier molecular orbital distributions of DPE-DDMAc calculated on B3LYP/6-31G* level in gas phase. Fig. S13. Front views of horizontal oriented backbone (-1722.38936220 a.u.) and vertical oriented backbone (-1722.38783975 a.u.) DPE-DPXZ optimized in gas phase in B3LYP/6-31G* level. We tentatively took the horizontal oriented backbone one as the initial guess for the subsequent simulation for its ground state energy is lower than the vertical oriented backbone one. Fig. S14. Electrostatic potential (ESP) distribution of two conformers of DPE-DPXZ. The positive electrostatic sites mainly localize on the DPE moiety for both molecules which induced moderate electron-withdrawing ability of the DPE unit.

-1.84 LUMO LUMO+1-4.72 HOMO-1 HOMO Fig. S15. Frontier molecular orbital distributions of DPE-DPXZ calculated on B3LYP/6-31G* level in gas phase. α DPE-DDMAc (α) DPE-DPXZ (β) Energy (a.u.) β O O -200-160 -120-80 -40 0 Dihedral angle ( ) Fig. S16. Potential surface scan of DPE-DDMAc and DPE-DPXZ in gas phase in B3LYP/6-31G* level. The scanning dihedral angles are α for DPE-DDMAc and β for DPE-DPXZ, respectively. Fig. S17. The molecular conformation of DPE-DDMAc in single crystal structures. Fig. S18. Molecular packing of DPE-DDMAc in single crystal structures.

Table S1. Photophysical properties of the investigated compounds. a) data measured in toluene solution (0.01 mm); b) data measured in neat film state; c) data measured in doped films (10 wt% in DPEPO); d) calculated results from TD-DFT simulation. ε/mm Materials λ abs /nm λ PL /nm QY/% c) τ F /ns τ TADF /μs cm -1 k F k ISC k IC k TADF ( 10 7 s -1 ) ( 10 7 s -1 ) ( 10 7 s -1 ) ( 10 3 s -1 ) E 0-0 (T 1 ) E 0-0 (S 1 ) (ev) ΔE ST /ev (ev) 23.0/25.0 DPE-DDMAc 355 a) /363 b) 3.60 a) 428 a) /415 b) /433 c) (±3.0) 6.70 c) / 3.73 c) 11.2 c) / 3.16 a), 1 CT 2.61 a), 3 LE 0.55 a) /0.44 d) 19.0/30.0 DPE-DPXZ 381 a) /390 b) 3.75 a) 468 a) /475 b) /475 c) (±3.0) 5.58 c) 542.7 c) 1.79 c) 11.9 c) 4.18 c) 0.55 c) 2.92 a), 1 CT 2.52 a), 3 LE 0.40 a) /0.18 d) DMAc 336 a) 3.69 a), 1 LE 3.12 a), 3 LE 0.57 a) PXZ 336 a) 3.69 a), 1 LE 3.13 a), 3 LE 0.56 a) DPE 323 a) 3.84 a), 3 LE 2.72 a), 3 LE 1.12 a)

Photo-physical Properties. a Molar extinction coefficient (mm cm -1 ) 40 30 20 10 0 300 400 500 600 Wavelength (nm) DPE-DDMAc DPE-DPXZ phos_dpe-ddmac phos_dpe-dpxz phos_dpe 1.5 1.2 0.9 0.6 0.3 0.0 PL intensity (a.u.) Absorption intensity (a.u.) 1.5 1.2 0.9 0.6 0.3 DPE-DPXZ_neat flim 10 wt.% DPE-DPXZ in DPEPO DPE-DDMAc_neat flim 10 wt.% DPE-DDMAc in DPEPO 0.0 0.0 200 300 400 500 600 700 Wavelength (nm) Fig. S19. a, UV-vis, PL and phosphorescence (at 77 K) spectra of the investigated compounds in toluene solution (0.01 mm). b, UV-vis absorption and photoluminescence spectra of neat films and doped films (10 wt% in DPEPO) of compounds DPE-DPXZ and DPE-DDMAc were shown, respectively. b 1.5 1.2 0.9 0.6 0.3 PL intensity (a.u.) 10 4 10 wt% DPE-DDMAc in DPEPO_vacuum Peak counts 10 3 10 2 10 1 0 40 80 120 160 200 Time (ns) 10 5 10 4 DPE-DPXZ neat film_air DPE-DPXZ neat film_vacuum Peak counts 10 3 10 2 10 1 0 40 80 120 160 200 Time (ns) Fig. S20. Transient decay spectra of the investigated compounds measured in air and vacuum condition at room temperature using a fluorescence lifetime spectrometer (C11367-01, Hamamatsu Photonics).

Thermal Properties. 110 100 DPE-DDMAc DPE-DPXZ 90 383 394 Mass (%) 80 70 60 50 T g =122 Exothermic 40 50 100 150 200 250 30 150 300 450 600 750 Temperature ( C) Fig. S21. Thermal properties of DPE-DDMAc and DPE-DPXZ were shown. The corresponding thermal decomposition temperatures (referring to 5% mass loss) were over 380, which were high enough for vacuum evaporation process. For DPE-DPXZ, the glass transition temperature of 122 was observed. Cyclic Voltammetry. 0.455 Ferrocene DPE-DDMAc DPE-DPXZ Current 0.539 0.889 0.679 0.0 0.5 1.0 1.5 Potential (V vs Ag/AgCl) Fig. S22. Cyclic voltammograms of two compounds measured in dichloromethane: acetonitrile (3:1) solution using n-bu 4 PF 6 as the supporting electrolyte. The detailed evaluation of ionization potential (IP) and electron affinity (EA) potentials were listed in Table 1.

OLED Characterization. Device Fabrication and Characterization: Glass substrates pre-coated with a 95-nm-thin layer of indium tin oxide (ITO) with a sheet resistance of 10 Ω per square were thoroughly cleaned in ultrasonic bath of acetone, isopropyl alcohol, detergent, deionized water, and isopropyl alcohol and treated with O 2 plasma for 20 min in sequence. Organic layers were deposited onto the ITO coated glass substrates by thermal evaporation under high vacuum (<5 10 4 Pa). Cathode, consisting of a 1 nm-thin layer of LiF followed by a 100 nm thin Al layer, was patterned using a shadow mask with an array of 3 mm 3 mm openings. Deposition rates are 1 2 Å s -1 for organic materials, 0.1 Å s -1 for LiF, and 6 Å s -1 for Al, respectively. Electroluminescence (EL) spectra were recorded by an optical analyzer, Photo Research PR705. The current density and luminance versus driving voltage characteristics were measured by Keithley 2420 and Konica Minolta chromameter CS-200. EQE was calculated from the luminance, current density, and EL spectrum, assuming a Lambertian distribution. Table S2. Summary of performances of OLED devices I, II and III. Device EBL V O * (V) Max. Efficiency CE max EQE max (cd/a) (%) EQE 100 * (%) L max (cd m -2 ) CIE (x, y) (@ 1mA cm 2 ) I TAPC 3.8 2.72 2.89 1.14 2599 (0.271, 0.434) II TAPC 3.3 7.27 3.52 2.14 5995 (0.176, 0.282) III DAcDB 3.4 22.8 10.8 2.58 718 (0.180, 0.314) * V O and EQE 100 are obtained at 1 and 100 cd m -2, respectively.

References 1. 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. akatsuji, M. Caricato, X. Li, H. P. Hratchian, A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. akajima, Y. Honda, O. Kitao, H. akai, T. Vreven, J. A. Montgomery Jr, J. E. Peralta, F. Ogliaro, M. J. Bearpark, J. Heyd, E.. Brothers, K.. Kudin, V.. Staroverov, R. Kobayashi, J. ormand, K. Raghavachari, A. P. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi,. Rega,. J. Millam, M. Klene, J. E. Knox, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, R. L. Martin, K. Morokuma, V. G. Zakrzewski, G. A. Voth, P. Salvador, J. J. Dannenberg, S. Dapprich, A. D. Daniels, Ö. Farkas, J. B. Foresman, J. V. Ortiz, J. Cioslowski and D. J. Fox, Gaussian 09, Revision B.01, 2009, Gaussian, Inc., Wallingford CT, U.S.A. 2. S. Huang, Q. Zhang, Y. Shiota, T. akagawa, K. Kuwabara, K. Yoshizawa and C. Adachi, J. Chem. Theory Comput., 2013, 9, 3872-3877. 3. Q. Zhang, B. Li, S. Huang, H. omura, H. Tanaka and C. Adachi, at. Photon., 2014, 8, 326-332.