Electroactive Conjugated Polymer / Magnetic Functional Reduced Graphene Oxide for Highly Capacitive Pseudocapacitors: Electrosynthesis, Physioelectrochemical and DFT Investigation

Article information

J. Electrochem. Sci. Technol. 2018;9(4):301-307
1Department of Chemistry, Faculty of science, University of Qom, Qom, Iran
2Department of Chemistry, Payame Noor University, Iran
3Department of chemistry, Amirkabir University of Technology, Tehran, Iran
Received 2018 July 07; Accepted 2018 July 31.

Abstract

The current study fabricated magnetic functional reduced graphene oxide (MFRGO) by relying on FeCl4- magnetic anion confined to cationic 1-methyl imidazolium. Furthermore, for improving the electrochemical performance of conductive polymer, hybrid poly ortho aminophenol (POAP)/ MFRGO films have then been fabricated by POAP electropolymerization in the presence of MFRGO nanorods as active electrodes for electrochemical supercapacitors. Surface and electrochemical analyses have been used for characterization of MFRGO and POAP/ MFRGO composite films. Different electrochemical methods including galvanostatic charge discharge experiments, cyclic voltammetry and electrochemical impedance spectroscopy have been applied to study the system performance. Prepared composite film exhibited a significantly high specific capacity, high rate capability and excellent cycling stability (capacitance retention of ~91% even after 1000 cycles). These results suggest that electrosynthesized composite films are a promising electrode material for energy storage applications in high-performance pseudocapacitors.

1. Introduction

The outspreading and development of next generation energy storage devices, supercapacitors (SCs) also called ultracapacitor, for varied applications such as portable electronic devices has received significant attentions in last decades as they are able to deliver high power density and show fast charge/discharge (C/D) capability beside of long life cycle. To improve the performance of supercapacitors and obtain higher electrochemical achievements a long range of materials have been used to consider their potential as electrode. Graphene and reduced graphene oxide (rGO) and carbon nano tubes (CNTs) and derivatives are promising materials to employ in SCs due to high specific capacitance delivery and surface area based on electrical double layer capacitance (EDLC)[1-5]. Functionalized GO (FGO) nanosheets have recently been developed as attractive fillers for electrode and membrane modifications. The functional groups of FGO could ease the dispersion of FGO in organic solvents. Additionally, the functional groups on the FGO structure can enhance the water retention and proton conductivity of the membranes.

Conductive polymers and metal oxide and their composites are other materials with high potential to develop supercapacitors based on pseudocapacitance. Many investigations have substantially carried out to consider the effects of hybrid supercapacitors based on EDLC and pseudocapacitance because this combination showed a synergistic effect and is an appropriate method to enhance electrochemical features [6-15]. Recently ionic liquids (ILs) have been employed increasingly in energy storage devices to improve their performances and great attractions due to the high thermal, chemical and electrochemical stability, low toxicity with good ionic conductivity, introduced them as an ideal material for SC applications. Moreover ionic liquids have been enthusiastically applied as electrolytes because of their wide potential window [16-25]. In addition, ionic liquids have also low viscosities with high conductivities advantages made them to employ in electrochemical devices for fast charge-discharge conditions. Pyrrolidinium and imidazolium are two main applied ionic liquids in electrochemical investigations [26-29]. Pyrrolidinium compared to imidazolium indicated larger electrochemical stability window (ESW) favorable for SCs operating at higher voltage or contributing mainly to the energy density. While Imidazolium actually provides higher ion conductivity with lower viscosity [30-33].

In this study, based on our recent achievements and successful results regarding to synthesis and application of new poly orthoaminophenol (POAP) composite films [34-46], we synthesized new magnetic functionalized reduced graphene oxide (MFRGO) (Fig. 1) and its nanocomposite with a POAP. The synthesized reduced graphene oxide and its nanocomposite were characterized by surface analyses and different electrochemical techniques, respectively. Furthermore, the potential of this novel nanocomposite as an efficient electrode material in electrochemical pseudocapacitors is demonstrated. The POAP as a conductive polymer when conjugated by functionalized graphene oxide, showed a synergistic effect being responsible for high supercapacitive performance.

Fig. 1.

Molecular structure of POAP and magnetic functional reduced graphene oxide (MFRGO). (Reprinted from [48] with permission from Elsevier Publications).

2. Experimental

All the chemical materials used in this work, obtained from Merck Chemical Co., were of analytical grade and used without further purification. Double distilled water was used throughout the experiments. All electrochemical experiments were carried out by a Potentiostat/galvanostat (Ivium V21508, Vertex). A conventional three electrode cell with an Ag/AgCl reference electrode (Argental, 3 M KCl) was used in order to carry out the electropolymerization of the POAP. A platinum wire and a carbon paste electrode was used as the counter and working electrodes respectively. Morphological investigations of the polymeric films were carried out by using SEM analysis.

For preparation of MFRGO, firstly, GO was synthesized according to the modified Hummers method from oxidation of natural graphite [47]. Modification of GO and then preparation of MRGO were carried out to as described in our already published paper [35,48]. POAP/MFRGO composites were prepared by electropolymerization in an acidic solution of monomer (ortho aminophenol) on the surface of the MFRGO modified working electrode. Electropoly-merizations were conducted by 40 consecutive cycles at the sweep rate of 50 mV·s−1 in the potentials between −0.2 to 0.9 V.

3. Results and Discussion

Fig. 2 shows the SEM images used to determine the morphology of the samples. Paper-like structures were the most common observations in all images and can be recognized as MFRGO. SEM image of POAP/MFRGO composite shows porous structure and presence of the MFRGO in the electrosynthesized composite film.

Fig. 2.

SEM image of FGO, MFRGO(Reprinted from [35 and 38] with permission from Elsevier Publications) and POAP/MFRGO composite film.

Fig. 3 gives CV curves of polymer and composite electrodes between −0.6V and 0.6V vs. Ag/AgCl in 0.1M HClO4 aqueous electrolyte solution. It can be seen that the shapes of the CV curves are more or less rectangular within the measured potential window. It is pointed out that the CV curves for the composite not only present a large background current, but also reveal obvious peak currents. The synergetic effect resulting from the interactions of POAP and MFRGO nmay affect the shape of CV curves. The CV of POAP/ MFRGO electrode shows the incorporation of MFRGO in POAP matrix not only increase the capacitance of composite states, but also save it’s ideal capacitive behavior. Shape of CV curves represents the pseudocapacitance nature of the electrode.

Fig. 3.

Cyclic voltammograms of POAP and POAP/MFRGO electrodes in 0.1M HClO4 at the sweep rate of 100 mV·s−1.

The surface morphology of POAP/ MFRGO film was studied by using the fractal concept. According to obtained results there is a power dependence between the peak current (Ipc) in cyclic voltammograms and the corresponding potential sweep rate (ν). Thus, the fractal parameter can be obtained easily by plotting the peak current against the sweep rate in log-log scale [35]. Based on this information, Fig. 4 represents cyclic voltammograms of POAP/MFRGO films that were recorded in different potential sweep rates in the range of 10-400 mVs−1. Fractal dimension has been obtained from relationship between the anodic peak current and potential sweep rate in log-log scale. The slope of the line gives the fractal dimensions of 2.60 for electrosynthesized film. The presented values for fractal dimension confirm porous structure of the electrosynthesized composite film.

Fig. 4.

Cyclic voltammograms of POAP/ MFRGO films in different potential sweep rates in the range of 10-400 mV·s−1.

To further evaluate the potential applications of the electrosynthesized nanocomposite as an electrode material for electrochemical supercapacitors, galvanostatic charging and discharging measurements were carried out between −0.5 and 0.5 V in 0.1 M HClO4 at 0.005 mA (Fig. 5a). All the charge-discharge curves of the nanocomposite, regardless of the current were symmetrical, indicating excellent electrochemical reversibility. From the discharge curves, the specific capacitance (Cs) was calculated according to the following equation [35]:

Fig. 5.

(a) Galvanostatic charge and discharge measurements of POAP/MFRGO electrode in 0.1M HClO4 solution at 0.005mA and (b) during 100 consecutive charge-discharge of composite film at 1mA.

Where I is the current loaded, m is reactive material mass, V is the potential change during discharge process and t is the discharge time. By substituting the obtained values in equation 1, the SC of POAP/MFRGO electrode was found to be 247 F·g−1. Furthermore, stability test results shows, using MFRGO in conductive polymer caused an excellent retention in stability percentage of composite electrode suggesting the good stability toward consecutive cycles (1000 cycles).

Electrochemical impedance spectroscopy (EIS) was carried out to investigate the electrode kinetics and other properties of as-prepared material [49-62]. The EIS data in Fig. 6 can be fitted by a solution resistance Rs, a charge-transfer resistance Rct, a pseudocapacitance (CPE2) from redox process of materials and a constant phase element (CPE1) to simulate the double-layer capacitance. Presented Nyquist plots shows that, the nanocomposites have faster electron transport in the bulk-film and charge transfer in the parallel POAP film/solution interface and MFRGO /solution interface. This fact may suggest that the MFRGO has an obvious improvement effect, which makes the composites have more active sites for faradic reactions and a larger specific capacitance than pure POAP [35].

Fig. 6.

Nyquist plots for POAP/MFRGO electrode.

At the end, electronic properties of the GO and MFRGO are investigated by DFT calculation. Nanoelectrochemical science has enabled us to produce new molecular structures, (such as phenylevinylene, thiophene and organic field-effect transistor) with desired properties, and thus can be used to design intelligent nano-electrochemical devices, using quantum mechanical concepts/methodology [63-65]. For example, details of the local charge and energy transfer mechanisms, in these systems/devices should be studied. For example, based on the quantum mechanical Methodology (DFT-method), the local charge and energy transfer in the graphen (G) and graphen oxide (GO) molecular systems are studied. In addition, geometry optimization and calculation of the structural and electronic/vibrational properties (such as energy gap between the frontier orbitals, HOMO and LUMO, electron density, spin density, virial force, local electrostatic force and vibrational frequency) of the G and GO molecular systems have been carried out using B3LYP/6-31G level of theory, using density functional theory (DFT) and within harmonic oscillator (HO) approximation. Sample of results obtained here are shown in Figs. 7-9 (all results not reported here for brevity).

Fig. 7.

Energies of the molecular orbitals (a), and the IR spectra (b) of the G molecular system, calculated at DFT-B3LYP/ 6-31G level of theory.

Fig. 9.

The contour maps of the local spin density (a), intramolecular virial force (b), and the HOMO (c)/LUMO (d) molecular orbital diagrams of the GO molecular system, are calculated at DFT-B3LYP/6-31G level of theory.

Fig. 8.

The contour maps of the local electron density (a), intramolecular virial force (b), and the HOMO (c)/LUMO (d) molecular orbital diagrams of the G molecular system, are calculated at DFT-B3LYP/6-31G level of theory.

Analysis of these results show that the value of the parallel and perpendicular charge/energy transfers in both G and GO molecular systems are considerable, which reflect intra-electrochemical phenomena in these systems. It can thus be predicted that the these intramolecular electrochemical response features depend on the vibrational and electronical structure/properties of these molecule systems and their variation with the external field or external bias voltage intensity. Also, it can be predicted that the reduced GO system device have a higher thermochemical performance than the G system when applied in a real electrochemical circuits. This higher performance can be attributed to the more extended π-conjugated system of the reduced GO system, induced by the applied external bias voltage, and thus the contributions of electrons/phonons (vibrational degrees of freedom) of the metallic Fe atoms, and π-conjugated electrons of pyridine groups in the intramolecular charge/energy transfer in this molecular system.

4. Conclusions

In this work, we successfully synthesized MFRGO and POAP/MFRGO composite material through electrochemical method. Multiple measurement methods were employed to study the performances of the material. Importantly, CV, CP and EIS tests show that POAP/ MFRGO composite material has better properties than POAP without MFRGO, suggesting it can be used as supercapacitor electrode material with excellent specific capacitance and ultrahigh specific power, which indicates this material is a promising electrode material used in high power applications.

References

[1]. Hou J, Shao Y, Ellis MW, Moore RB, Yi B. Phys Chem Chem Phys 2011;13(34):15384–15402. 10.1039/c1cp21915d.
[2]. Mishra AK, Ramaprabhu S. J. Phys. Chem 2011;115(29):14006–14013. 10.1021/jp2086736.
[3]. Yu Z, McInnis M, Calderon J, Zhai L, Seal S, Thomas J. Nano Energy 2015;11:611–620. 10.1016/j.nanoen.2014.11.030.
[4]. Aghazadeh M, Ganjali MR, Norouzi P. J. Mater. Sci - Mater. Electron 2016;27(7):7707–7714. 10.1007/s10854-016-4757-1.
[5]. Aghazadeh M, Ganjali MR, Norouzi P. Thin Solid films 2017;634:24–32. 10.1016/j.tsf.2017.05.008.
[6]. Liu C, Yu Z, Neff D, Zhamu A, Jang BZ. Nano Lett 2010;10(12):4863–4868. 10.1021/nl102661q.
[7]. Dey RS, Hjuler HA, Chi Q. J Mater Chem A 2015;3(12):6324–6329. 10.1039/C5TA01112D.
[8]. Frackowiak E., Jurewicz K., Delpeux S., Béguin F.. J. Power Sources 2001;97:822–825.
[9]. Portet C., Taberna P. L., Simon P., Flahaut E.. J. Power Sources 2005;139(1-2):371–378. 10.1016/j.jpowsour.2004.07.015.
[10]. Reddy Ravinder N., Reddy Ramana G.. J. Power Sources 2006;156(2):700–704. 10.1016/j.jpowsour.2005.05.071.
[11]. Reddy Ravinder N., Reddy Ramana G.. J. Power Sources 2004;132(1-2):315–320. 10.1016/j.jpowsour.2003.12.054.
[12]. Kim TY, Lee HW, Stoller M, Dreyer DR, Bielawski CW, Ruoff RS. ACS Nano 2011;5(1):436–442. 10.1021/nn101968p.
[13]. Aghazadeh M, Ganjali MR. J. Mater. Sci - Mater. Electron 2017;28(11):8144–8154. 10.1007/s10854-017-6521-6.
[14]. Aghazadeh M, Ganjali MR. J. Mater. Sci 2018;53(1):295–308. 10.1007/s10853-017-1514-7.
[15]. Aghazadeh M, Bahrami-Samani A, Gharailou D, Ghannadi Maragheh M, Ganjali MR, Norouzi P. J. Mater. Sci - Mater. Electron 2017;27(11):11192–11200.
[16]. Bhise S.C., Awale D.V., Vadiyar M.M., Patil S.K., Kokare B.N., Kolekar S.S.. J. Solid State Electrochem 2017;21(9):2585–2591. 10.1007/s10008-016-3490-2.
[17]. Liew C.-W., Ramesh S., Arof A.K.. Materials & Design 2016;92:829–835. 10.1016/j.matdes.2015.12.115.
[18]. Izmailova M.Y., Rychagov A.Y., Den’shchikov K.K., Vol’fkovich Y.M., Lozinskaya E.I., Shaplov A.S.. Russ. J. Electrochem 2009;45(8):949–950. 10.1134/S1023193509080205.
[19]. Ekka D., Roy M.N.. Ionics 2014;20(4):495–505. 10.1007/s11581-013-1003-1.
[20]. Salanne M.. Top. Curr. Chem 2017;375(3):63. 10.1007/s41061-017-0150-7.
[21]. Patil S.K., Vadiyar M.M., Bhise S.C., Patil S.A., Awale D.V., Ghorpade U.V., Kim J.H., Ghule A.V., Kolekar S.S.. J. Mater. Sci - Mater. Electron 2017;28(16):11738–11748. 10.1007/s10854-017-6978-3.
[22]. Boruń A., Bald A.. Ionics 2016;22(6):859–867. 10.1007/s11581-015-1613-x.
[23]. Qu S., Sun Y., Li J.. Ionics 2017;23(6):1607–1611. 10.1007/s11581-017-2114-x.
[24]. Tripathi M., Tripathi S.K.. Ionics 2017;23(10):2735–2746. 10.1007/s11581-017-2051-8.
[25]. Liew C.-W., Ramesh S., Arof A.K.. Int. J. Hydrogen. Energ 2014;39(6):2953–2963. 10.1016/j.ijhydene.2013.06.061.
[26]. Liu G., Ma Y., Hou X., huang Y., Chen J., Zhan G., Li C.. Ionics 2015;21(9):2567–2574. 10.1007/s11581-015-1455-6.
[27]. Lakshminarayana G., Tripathi V.S., Tiwari I., Nogami M.. Ionics 2010;16(5):385–395. 10.1007/s11581-010-0436-z.
[28]. Xu P., Gui H.-g., Ding Y.-s.. Ionics 2013;19(11):1579–1585. 10.1007/s11581-013-0901-6.
[29]. Huang L., Yao X., Yuan L., Yao B., Gao X., Wan J., Zhou P., Xu M., Wu J., Yu H., Hu Z., Li T., Li Y., Zhou J.. Energy Storage materials 2018;12:191–196. 10.1016/j.ensm.2017.12.016.
[30]. Kowsari E., Ehsani A., Dashti Najafi M., Seifvand N., Heidari A.A.. Ionics 2018;
[31]. Randström S., Appetecchi G.B., Lagergren C., Moreno A., Passerini S.. Electrochim. Acta 2007;53(4):1837–1842. 10.1016/j.electacta.2007.08.029.
[32]. Lin R., Taberna P.-L., Fantini S., Presser V., Pérez C.R., Malbosc F., Rupesinghe N.L., Teo K.B., Gogotsi Y., Simon P.. J. Phys. Chem. Lett. 2011;2(19):2396–2401. 10.1021/jz201065t.
[33]. Arbizzani C., Biso M., Cericola D., Lazzari M., Soavi F., Mastragostino M.. J. Power Sources 2008;185(2):1575–1579. 10.1016/j.jpowsour.2008.09.016.
[34]. Ehsani A., Khodayari J., Hadi M., Mohammad Shiri H., Mostaanzadeh H.. Ionics 2017;23(1):131–138. 10.1007/s11581-016-1811-1.
[35]. Ehsani A., Mohammad Shiri H., Kowsari E., Safari R., Shabani Shayeh J., Barbary M.. J. Colloid Interface. Sci 2017;490:695–702. 10.1016/j.jcis.2016.12.003.
[36]. Ehsani A., Mohammad Shiri H., Kowsari E., Safari R., Torabian J., Kazemi S.. J. Colloid interface. Sci 2016;478:181–187. 10.1016/j.jcis.2016.06.013.
[37]. Mohammad Shiri H., Ehsani A.. J. Colloid Interface. Sci 2016;5:91062–91068.
[38]. Mohammad Shiri H., Ehsani A.. J. Colloid Interface. Sci 2016;484:70–76. 10.1016/j.jcis.2016.08.075.
[39]. Ehsani A., Kowsari E., Boorboor Ajdari F., safari R., Mohammad Shiri H.. J. Colloid Interface. Sci 2018;512:151–157. 10.1016/j.jcis.2017.10.046.
[40]. Ehsani A., Mohammad Shiri H., Kowsari E., Safari R., Shabani Shayeh J., Barbary M.. J. Colloid Interface. Sci 2017;490:695–702. 10.1016/j.jcis.2016.12.003.
[41]. Mohammad Shiri H., Ehsani A.. J. Colloid interface. Sci 2016;473:126–131. 10.1016/j.jcis.2016.03.065.
[42]. Mohammad Shiri H., Ehsani A.. Bull. Chem. Soc. Jpn 2016;89(10):1201–1206. 10.1246/bcsj.20160082.
[43]. Naseri M., Fotouhi L., Ehsani A., Mohammad Shiri H.. J. Colloid interface. Sci 2016;484:308–313. 10.1016/j.jcis.2016.08.071.
[44]. Naseri M., Fotouhi L., Ehsani A., Dehghanpour S.. J. Colloid interface. Sci 2016;484:314–319. 10.1016/j.jcis.2016.09.001.
[45]. Ehsani A., Mahjani M.G., Bordbar M., Moshrefi R.. Synth. Met 2013;165:51–55. 10.1016/j.synthmet.2013.01.004.
[46]. Shabani Shayeh J., Sadeghinia M., Omid Ranaei Siadat S., Ehsani A., Rezaei M., Omidi M.. J. Colloid Interface Sci 2017;496:401–406. 10.1016/j.jcis.2017.02.010.
[47]. Hummers WS Jr, Offeman RE.. J. Am. Chem. Soc 1958;80(6):1339. 10.1021/ja01539a017.
[48]. Kowsari E., Mohammadi M.. Compos. Sci. Technol. 2016;126:106–114. 10.1016/j.compscitech.2016.02.019.
[49]. Ehsani A.. Prog. Org. Coat 2015;78:133–139. 10.1016/j.porgcoat.2014.09.015.
[50]. Ehsani A., Mahjani M.G., Moshrefi R., Mostaanzadeh H., Shayeh J.S.. RSC Advances 2014;4:20031–20037. 10.1039/C4RA01029A.
[51]. Mohammad Shiri H., Ehsani A., Jalali Khales M.. J. Colloid interface. Sc 2017;505:940–946. 10.1016/j.jcis.2017.06.086.
[52]. Hosseini M., Fotouhi L., Ehsani A., Naseri M.. J. Colloid interface. Sci 2017;505:213–219. 10.1016/j.jcis.2017.05.097.
[53]. Aljourani J., Raeissi K., Golozar M.A.. Corros. Sci. 2009;5(8):1836–1843.
[54]. Özcan M., Dehri İ., Erbil M.. Appl. Surf. Sci. 2004;236(1-4):155–164. 10.1016/j.apsusc.2004.04.017.
[55]. Bouklah M., Hammouti B., Aouniti A., Benkaddour M., Bouyanzer A.. Appl. Surf. Sci 2006;252(18):6236–6242. 10.1016/j.apsusc.2005.08.026.
[56]. Elkadi L., Mernari B., Traisnel M., Bentiss F., Lagrenée M.. Corros. Sci 2000;42(4):703–719. 10.1016/S0010-938X(99)00101-8.
[57]. Tebbji K., Hammouti B., Oudda H., Ramdani A., Benkadour M.. Appl. Surf. Sci 2005;252(5):1378–1385. 10.1016/j.apsusc.2005.02.097.
[58]. Noor E. A.. Mater. Chem. Phys 2011;131:160–169. 10.1016/j.matchemphys.2011.08.001.
[59]. Bentiss F., Lebrini M., Lagreńee M., Traisnel M., Elfarouk A., Vezin H.. Electrochim. Acta. 2007;52:6865–6872. 10.1016/j.electacta.2007.04.111.
[60]. Saha S.K., Dutta A., Ghosh P., Sukul D., Banerjee P.. Phys. Chem. Chem. Phys 2015;17:5679–5690. 10.1039/C4CP05614K.
[61]. Fotouhi L., Fathali N., Ehsani A.. Int. J. Hydrogen. Energ 2018;43(14):6987–6996. 10.1016/j.ijhydene.2018.02.123.
[62]. Naseri M., Fotouhi L., Ehsani A.. J. Electrochem. Sci. Technol 2018;9(1):28–36.
[63]. Iniewski K.. Nanoelectronics: Nanowires, Molecular Electronics,and Nanodevices McGraw-Hill; 2010.
[64]. Feringa B. L.. Molecular Switches Wiley. Weinheim: 2007.
[65]. Matta C.F., Boyd R.J.. Quantum Biochemistry Wiley. Weinheim: 2010.

Article information Continued

Fig. 1.

Molecular structure of POAP and magnetic functional reduced graphene oxide (MFRGO). (Reprinted from [48] with permission from Elsevier Publications).

Fig. 2.

SEM image of FGO, MFRGO(Reprinted from [35 and 38] with permission from Elsevier Publications) and POAP/MFRGO composite film.

Fig. 3.

Cyclic voltammograms of POAP and POAP/MFRGO electrodes in 0.1M HClO4 at the sweep rate of 100 mV·s−1.

Fig. 4.

Cyclic voltammograms of POAP/ MFRGO films in different potential sweep rates in the range of 10-400 mV·s−1.

Fig. 5.

(a) Galvanostatic charge and discharge measurements of POAP/MFRGO electrode in 0.1M HClO4 solution at 0.005mA and (b) during 100 consecutive charge-discharge of composite film at 1mA.

Fig. 6.

Nyquist plots for POAP/MFRGO electrode.

Fig. 7.

Energies of the molecular orbitals (a), and the IR spectra (b) of the G molecular system, calculated at DFT-B3LYP/ 6-31G level of theory.

Fig. 8.

The contour maps of the local electron density (a), intramolecular virial force (b), and the HOMO (c)/LUMO (d) molecular orbital diagrams of the G molecular system, are calculated at DFT-B3LYP/6-31G level of theory.

Fig. 9.

The contour maps of the local spin density (a), intramolecular virial force (b), and the HOMO (c)/LUMO (d) molecular orbital diagrams of the GO molecular system, are calculated at DFT-B3LYP/6-31G level of theory.