Electronic Structure of Edge-Modified Graphene Quantum Dots

Jul 14, 2017 - ... Química, Universidad de la República, C.P. 11800, Montevideo, Uruguay ... Multilayer phosphorene quantum dots in an electric fiel...
0 downloads 0 Views 3MB Size
Article pubs.acs.org/JPCC

Electronic Structure of Edge-Modified Graphene Quantum Dots Interacting with Polyaniline: Vibrational and Optical Properties Dominique Mombrú,† Mariano Romero,*,† Ricardo Faccio,*,† and Á lvaro W. Mombrú*,† †

Centro NanoMat/CryssMat/Física, DETEMA, Facultad de Química, Universidad de la República, C.P. 11800, Montevideo, Uruguay S Supporting Information *

ABSTRACT: There are recent experimental reports on the preparation of edge-modified graphene quantum dots (GQDs) and its loading into conductive polymers for energy applications. Recently, we have studied the influence of the graphene quantum dots (GQDs) addition to polyaniline (PANI) conductive polymer in the electronic transport, finding important modifications in the structural and electrical properties respect to the pristine polymer. In this work, we present a mixed approach, based in firstprinciples molecular dynamics with time-dependent density functional calculations, for the description of the electronic dynamics, particularly for the determination of the optical properties. We found that edge modification of the GQDs favors the PANI/GQDs interaction, promoting very interesting electronic properties. Particularly, when hydroxyl chemical functional groups are present at the edges of the GQD a very good electron−hole separation during the photon excitation is obtained, demonstrating its potential use for optical applications in general and solar cells in particular.

1. INTRODUCTION There is a recent interest in the preparation and characterization of conductive polymer-based devices for photoresponse,1−3 sensing,4 and energy conversion5 applications. Polyaniline (PANI) is one of the most studied conductive polymer due to its easy preparation, low cost, and high conductivity.6 There are recent reports on the preparation and electrical transport properties on PANI-based nanocomposites, including the additions of different materials, particularly graphene quantum dots (GQDs).2,3 The literature shows very interesting works that use first-principles methods to understand the mechanism of growth of carbon nanostructures, see, for example, Goyenola et al.7 However, there are very few theoretical insights regarding the edge modification chemical nature of these GQDs and its electronic interactions with conductive polymers. Our research group contributed experimentally for the understanding of the correlations between the chemical structure and its electronic properties,2 although it was only for the case of aggregates of small percentage of GQD to PANI. The interplay between donor−acceptor properties are crucial, not only for the electrical transport, but also for the optical response that could be of interest for solar cells.8 Many of these properties are directly affected for the electronic structure of the systems, thus it is very relevant to perform © 2017 American Chemical Society

quantum-chemical simulations in order to get insight on these properties. However, due to the complexity and the considerable number of atoms involved in a “first-principles” calculations, it is necessary to apply mixed methodologies. That is the reason why in literature very few reports on the topic are found, considering the modeling but utilizing small size PANI monomers,9 or even performing periodic boundary calculations with small unit cells that limit the potential geometry optimization of the corresponding polymers.10−12 Here we present a very first approximation for the quantum chemical modeling of nanocomposite systems, considering a moderate number of atoms to describe the interaction between PANI-GQDs nanocomposites. Our methodology involves first a density functional theory (DFT) molecular dynamics (MD) simulation, followed by a time-dependent density functional theory (TD-DFT) calculation, which is necessary to obtain an accurate description of excited states for a better description of vertical transitions. This last procedure allows us to get right optical and electronic structure description of the system. Received: April 17, 2017 Revised: July 6, 2017 Published: July 14, 2017 16576

DOI: 10.1021/acs.jpcc.7b03604 J. Phys. Chem. C 2017, 121, 16576−16583

Article

The Journal of Physical Chemistry C

Figure 1. Structure and geometry for PANI−X with (a) X = H, (b) X = OH, (c) X = COOH, (d) X = COOHOH, and (e) PANI. In the lower panel, we present the electronic density differences maps (EDDM) for the most intense optical transition. Note: Yellow regions indicate a loss of electron density in the transition to the excited state and blue regions indicate a gain of electron density in a transition to the excited state.

2. METHODS As pointed out in the early section, our objective is to perform precisely first-principles calculations involving a considerable number of atoms. This task is limited by its computational cost but it can be addressed by combining molecular dynamics simulations, followed by well suited quantum-chemical calculations. We considered four types of edge chemical functionalization named as GQD−X, where X= H, OH, COOH, and COOHOH, which correspond to the following empirical formulas: C150H30, C160H30O10, C160H30O20, and C160H30O15, respectively. These sizes of the molecules have been demonstrated as adequate for describing graphene or graphene oxide successfully from the vibrational point of view.13 In the case of the PANI polymer, we considered the emeraldine salt ([C6H4NH]2[C6H4N]2)n with n = 7, see Figure 1. In order to get reasonably good starting geometries we perform DFT14,15 utilizing Gaussian 09,16 selecting the B3LYP17−20 xc-potential at the 6-31g(d) level of theory. Then, vibrational analysis is performed in order to confirm the stability and local minima in all the starting compounds. For the modeling of the polymeric PANI/GQDs nanocomposite, we performed molecular dynamics simulations, in the framework of DFT utilizing the SIESTA21,22 code. SIESTA adopts a linear combination of numerical atomic orbitals. We utilized double-ζ basis sets with polarization orbitals (DZP) for

all the atoms. We choose the van der Waals (vdW)-DF xcpotential as implemented by Dion et al.23 within the LMKLL reparametrization of Lee et al.24 We selected a mesh cutoff of MC = 200 Ry that ensures a dense grid for charge density of the system. With the aim to describe an isolated system, we constructed a unit cell sufficiently large for avoiding interaction among periodic images, thus choosing only the Γ-point for the Brillouin-zone sampling. After setting up the SIESTA simulation, we proceeded with the molecular dynamics simulation, selecting a NVT ensemble with a T = 300 K, and time step of 1 fs and a total time of 2 ns. The Nosé−Hoover thermostat was employed in this setup. This procedure allows us to map all the configurational space and thus provide an accurate description of this polymeric system. Finally, we proceed with time-dependent DFT (TDDFT)25,26 calculations, utilizing the same conditions for the Gaussian09 set up, selecting the first 20 singlet−singlet excitations for optical study. With this procedure, we will report on the energy interactions between GQD and PANI, and describe its potential use for optical applications.

3. RESULTS 3.1. Geometry Optimization and Vibrational Properties. As mentioned before, the basis of 150 carbon atoms for describing the graphene-like core in all the molecules have been demonstrated as adequate for the study of its vibrational 16577

DOI: 10.1021/acs.jpcc.7b03604 J. Phys. Chem. C 2017, 121, 16576−16583

Article

The Journal of Physical Chemistry C

Figure 2. (Top rows) structure; (bottom rows) electronic density differences maps (EDDM) depicting the most relevant optical transitions for PANI-GQD−X with (a) X = H, (b) X = OH, (c) X = COOH, and (d) X = COOHOH. Note: Yellow regions indicate a loss of electron density in the transition to the excited state and blue regions indicate a gain of electron density in a transition to the excited state.

Figure 3. (a) IR spectra for GQD−X with X = H, COOHOH, OH, and COOH and (b) Raman spectra for GQD−X, with X = H, COOHOH, OH, and COOH.

properties.13 Anyway, the direct interpretation of the vibrational modes seems to be complex. In one hand, this is due to large number of atoms and thus the large number of degrees of freedom involved in the system. On the other hand, the chemical substitution breaks the symmetry and the degeneracy of modes, giving places to broad bands in the IR and Raman spectra too. The first step involves the geometrical optimization of the four oxidized GQDs; the final geometries are presented in Figures 1 and 2. As mentioned before, we performed edge chemical doping, because edges show better formation energies27 due to their higher chemical reactivity when compared to “bulk” atoms. In order to confirm the local minima of the obtained geometries, after performing the full

DFT geometrical optimization for the four structures we performed vibrational analysis including infrared (IR) and Raman spectra, see Figure 3. In all cases, the graphene plane remains flat with the most important structural changes localized at the edges. In all the cases, we obtain all real frequencies and the spectra provides the main features for an easy identification of the different chemical functional groups. When the Raman spectra are analyzed, there are no major differences between the different GQDs, see Figure 3. This is a consequence to the size of the GQDs, where the effect of the edge doping seems to be diluted and the differences becomes reduced. All the Raman spectra present the most intense peaks located at 1310 ± 10 cm−1, corresponding to the graphene-like 16578

DOI: 10.1021/acs.jpcc.7b03604 J. Phys. Chem. C 2017, 121, 16576−16583

Article

The Journal of Physical Chemistry C

Figure 4. (a) Temperature and (b) total energy (Etotal) evolution during the molecular dynamics simulation for PANI/GQD−H.

Table 1. Wavelength, Oscillation Strength, and Major Contribution for PANI−H, GQD−H, and PANI-GQD−H PANI−H GQD−H PANI-GQD−H

wavelength (nm)

oscillation strength

major contribution

358.2 360.0 748.1 749.9 742.3 730.4

0.4669 0.391 1.2469 1.2429 1.3218 0.9398

HOMO−2→LUMO+5 (25%) HOMO−3→LUMO+4 (18%) HOMO→LUMO+8 (29%) HOMO−1→LUMO (47%) HOMO→LUMO+1 (47%) HOMO−1→LUMO+1 (49%) HOMO→LUMO (45%) HOMO−5→LUMO (46%) HOMO−2→LUMO+1 (44%) HOMO−5→LUMO+1 (72%)

breathing band, ascribed to the stretching ν(C−C), usually labeled as the graphene’s G-band. In the case of bending modes of X−H, minor peaks can be observed at both sides of the main peak in the range of ∼1200−1700 cm−1. The peaks close to ∼1200 cm−1 correspond to δ(C−H) bending modes. Then, the peaks close to ∼1700 cm−1 correspond to stretching plus bending ν(C−C) + δ(C−H) modes, mainly located at edge positions in this case. In the region of ∼3150−3300 cm−1, the stretching ν(C−H) modes are identified. When considering the others GQDs, the contribution of −OH and −COOH functional groups are observed. As mentioned before, the bands in the Raman spectra changes slightly due to the presence of these groups and the symmetry breaking. In particular, the signals coming from the −COOH functional groups appeared as very weak in the region of ∼1800 cm−1, corresponding to the stretching ν(CO) mode. In the region of ∼3650−3710 cm−1, the stretching ν(O−H) modes appeared for X = OH, COOH, and COOH. When the IR spectra are analyzed, the situation changes due to the nature of the IR’s selection rule, and the differences between all the GQDs are more evident. Figure 3 include the sketches of the most relevant vibrational modes of the IR spectrum. It is evident the contribution form the stretching ν(CO) mode of the −COOH groups. The contribution of the bending δ(O−H) modes are observed in the region of ∼650 cm−1, contribution to the broadening of the whole vibrational spectra, as mentioned before. Even though the vibrational differences are not so significant for a straight identification, the presence of characteristic bands, coming from the −COOH and −OH chemical functional groups, are possible to be identified. These characteristics will

become very different when addressing the optical properties, see further. 3.2. MD Simulations. The MD simulations demonstrate the good affinity between PANI and all the studied GQDs. An important distortion of the graphene plane is observed with an interesting interaction with the PANI polymer, in particular with the central region of the GQDs (see Figure 2). On this basis, we can state that there is a good affinity between both systems, justifying the study of its electronic structure, in particular describing the optical properties. For reference, we present in Figure 4 the total energy and temperature evolution during the MD simulation in SIESTA, demonstrating final stabilization of the system PANI/GQD−H. In order to study the interaction between PANI and GQD− X system we calculated the interaction energy (Eint) according to Eint = EPANI/GQD−X − EPANI − EGQD−X, where EPANI/GQD−X, EPANI and EGQD−X corresponds to the corrected total energy of PANI/GQD−X system, isolated PANI, and isolated GQD−X, respectively. Because of the localized nature of the utilized basis set, a “basis set superposition error” (BSSE), was incorporated by the implementation of total energy corrections, utilizing a floating basis set denominated as “ghost atoms” in the SIESTA code. The interaction energies correspond to −0.21, −2.66, −5.11, and −5.92 eV for X = H, OH, COOH, and COOHOH, respectively. There is not a clear tendency on the interactions energies, but what is clear is the necessity of the presence of oxidizing groups for the enhancement of the interaction between PANI and modified GQD. 3.3. Optical Properties. In order to study the optical properties of the PANI/GQD composites, we performed TDDFT calculations, describing the first 20 singlet−singlet transitions. We started this task, modeling the isolated GQDs 16579

DOI: 10.1021/acs.jpcc.7b03604 J. Phys. Chem. C 2017, 121, 16576−16583

Article

The Journal of Physical Chemistry C Table 2. Wavelength, Oscillation Strength, and Major Contribution for PANI−OH, GQD−OH, and PANI-GQD−OH PANI−OH GQD−OH PANI-GQD−OH

wavelength (nm)

oscillation strength

major contribution

371.9 383.7 921.8 910.8 1375.2 1517.0

0.1111 0.1046 0.834 0.4548 0.054 0.0284

HOMO−2→LUMO+1 (78%) HOMO→LUMO+3 (42%) HOMO−1→LUMO+3 (34%) HOMO−1→LUMO+1 (60%) HOMO→LUMO+1 (28%) HOMO−4→LUMO (55%) HOMO−8→LUMO (30%)

Table 3. Wavelength, Oscillation Strength, and Major Contribution for PANI−COOH, GQD−COOH, and PANI-GQD− COOH PANI−COOH GQD−COOH PANI-GQD−COOH

wavelength (nm)

oscillation strength

major contribution

372.2 350.9 773.2 772.0 915.4

0.7049 0.3953 1.2795 1.232 0.062

HOMO−1→LUMO+1 (73%) HOMO−1→LUMO+6 (38%) HOMO−1→LUMO (37%) HOMO→LUMO+1 (35%) HOMO−1→LUMO+1 (45%) HOMO→LUMO (26%) HOMO−5→LUMO (53%)

Table 4. Wavelength, Oscillation Strength, and Major Contribution for PANI−COOHOH, GQD−COOHOH, and PANIGQD−COOHOH PANI−COOHOH GQD−COOHOH PANI-GQD−COOHOH

wavelength (nm)

oscillation strength

major contribution

364.2 378.6 853.6 815.6 1107.1

0.3445 0.2525 0.925 0.6124 0.1225

HOMO→LUMO+7 (33%) HOMO→LUMO+6 (37%) HOMO→LUMO+7 (21%) HOMO−1→LUMO+1 (58%) HOMO→LUMO (30%) HOMO−3→LUMO (27%) HOMO−1→LUMO (24%) HOMO−7→LUMO (38%) HOMO−5→LUMO (22%)

Figure 5. UV−vis spectra for (a) PANI, GQD−H, and PANI-GQD−H; (b) PANI, GQD−OH and PANI-GQD−OH;( c) PANI, GQD−COOH and PANI-GQD−COOH; and (d) PANI, GQD−COOHOH and PANI-GQD−COOHOH.

16580

DOI: 10.1021/acs.jpcc.7b03604 J. Phys. Chem. C 2017, 121, 16576−16583

Article

The Journal of Physical Chemistry C

modified, presenting important similarities to GQD−H in difference to GQD−OH, see Figure 1. Now the main transitions are close located at 772.0 and 773.2 nm, corresponding to HOMO−1 → LUMO+1 and HOMO−1 → LUMO with oscillation strength of 1.232 and 1.280, respectively. Finally, we discuss the GQD−COOHOH system, which due to the presence OH groups shows a situation similar to GQD− OH, see Figure 1. The main transitions are located at 815.6 and 853.6 nm with the main contributions of HOMO−3 → LUMO and HOMO−1 → LUMO+1 with oscillation strength 0.612 and 0.925, respectively. In general, GQD−H and GQD−COOH present the main absorption peaks in a shorter range in wavelength and with stronger oscillation strength in comparison to GQD−OH and GQD−COOHOH. These differences arise from the nature of the frontier orbital due to the edge doping. PANI/GQDs. When considering the PANI/GQD nanocomposite, the interaction induces important changes in the geometry, see Figure 2 and the corresponding electronic structure. As it was mentioned before, the geometrical structure and the corresponding binding energy indicates that GQD/ PANI interaction is important. In the case of PANI/GQD−H, there exists a blue shift in the absorption band in comparison to GQD−H with transitions located at 730.4 and 742.3 nm, described as HOMO−5 → LUMO+1 and HOMO−5 → LUMO with oscillation strength 0.940 and 1.322, respectively. With a closer inspection of the transitions, it could be understood as a PANI → GQD−H electronic transition, as observed in Figure 2. Nevertheless, there are some internal π* → π* from the GQD−H that decrease the degree of charge separation in the excitation process. Here, it is important to mention that when comparing to pure PANI and GQDs, the oscillation strength and thus the optical absorption is drastically reduced for the nanocomposite, but this is just for comparison purposes with the corresponding isolated systems. When going to the PANI/GQD−OH case, we observe again a red shift of the absorption band in comparison to GQD−OH. The transitions are located at 1375.2 and 1517.0 nm, described as HOMO−4 → LUMO and HOMO−8 and LUMO, with oscillation strengths 0.054 and 0.0284, respectively. When analyzing the topography of the electronic transition, it is observed a clear π* → π* excitation from PANI to GQD−OH, thus constituting a very good charge separation during the transition. In the case of PANI/GQD−COOH, the optical band exhibits a red shift when compared to GQD−COOH. The main transition corresponds to 915.4 nm, mainly represented as HOMO−5 → LUMO, with an oscillation strength of 0.062. In difference to PANI/GQD−H and PANI/GQD−OH, the excitation for PANI/GQD−COOH can be ascribed as a pure internal π* → π* transition in the core of GQD−COOH. For PANI/GQD−COOHOH, there exists a red shift in the absorption band, when compared to GQD−COOHOH. The transitions are located 1107.1 nm, described as HOMO−7 → LUMO with oscillation strength of 0.123. When analyzing the topography of the excitation process, it is observed again a clear π* → π* excitation from PANI to GQD−COOHOH, similarly to PANI/GQD−OH, thus constituting a very good charge separation during the transition. Finally, an important and interesting aspect is the study of the dependence of electronic properties as a function of the

and PANI. Figure 2 presents the electronic density differences maps (EDDM) for the most intense optical transition, describing the gain and loss in the electron density as a consequence of the optical transition. Because of the organic chemical nature of the systems, all the transitions can be described as π*→ π* for the GQDs and PANI. Tables 1, 2, 3, and 4 and Figure 5 collect all the relevant information for the optical properties of the system. PANI. In the case of isolated PANI, we calculated the optical properties for this polymer in the final conformation for the four PANI/GQDs geometries. The differences in the structural conformations induce changes in the oscillation strength, but there are no relevant changes in the corresponding wavelength. For example, in the case of PANI−H the main electronic transitions are located at 358.2 and 360 nm with the most relevant contribution coming from HOMO−2 → LUMO+5 and HOMO → LUMO+8, leading to an oscillation strength of 0.467 and 0.391 respectively. In the case of PANI−OH, the main transitions are located at 371.8 and 383.7 nm with the main contributions of HOMO−2 → LUMO+1 and HOMO → LUMO+3, and the oscillation strength is reduced to 0.111 and 0.105 respectively. For PANI−COOH, the situation changes: the transitions are located at 350.9 and 372.2 nm with the main contribution from HOMO−1 → LUMO+6 and HOMO−1 → LUMO+1 but in this case the oscillation strength increased up to 0.3952 and 0.7049. Finally, for PANI-COOHOH, the main transitions are located at 364.2 and 378.6 nm, corresponding to the following main contributions, HOMO → LUMO+7 and HOMO → LUMO+6, with an oscillation strength corresponding to 0.344 and 0.253, respectively. As it can be observed, the main transition is always located in a limited range of 350 to 379 nm but with important changes in the oscillation strength that strongly depends on the chain conformation. In particular, when polymer chain becomes linear the oscillation strength increases. GQDs. Here we proceed with an analogous treatment of the electronic structure of the isolated GQDs in their modified conformations when interacting with PANI. Here again, the most relevant optical transitions can be understood as π*→ π* in general with an important role of the states coming from the central region of the GQD, that could be seen as a graphenelike region. The presence of different chemical substitutes, X = H, OH, COOH, and COOHOH, clearly changes the pure graphene electronic structure. In the case of GQD−H, the situation is the closest to the ideal graphene system with clear π*→ π* transitions, highly delocalized in the core of the graphene sheet. The main transitions are located at 748.1 and 749.9 nm, described as HOMO−1 → LUMO/LUMO+1 and HOMO−1 → LUMO+1 with oscillation strength 1.247 and 1.243, respectively. The situation changes for GQD−OH, not only due to the conformational changes of the central graphene region but mainly because of the change given by the breakage of symmetry introduced by edge doping, see Figure 1. Here, an important contribution to the frontiers orbitals is observed with states coming from the edges of the structure. That is the reason why the main transitions now are more separated, located at 910.8 and 921.8 nm in comparison to other cases. The main transitions correspond to HOMO → LUMO+1 and HOMO−1 → LUMO+1 with oscillation strength 0.455 and 0.834, respectively. In the case of GQD−COOH, the electronic π-states delocalization in the core of the graphene region is slightly 16581

DOI: 10.1021/acs.jpcc.7b03604 J. Phys. Chem. C 2017, 121, 16576−16583

Article

The Journal of Physical Chemistry C

(4) Lin, W.-P.; Liu, S.-J.; Gong, T.; Zhao, Q.; Huang, W. PolymerBased Resistive Memory Materials and Devices. Adv. Mater. 2014, 26, 570−606. (5) Abdulrazzaq, O.; Bourdo, S. E.; Woo, M.; Saini, V.; Berry, B. C.; Ghosh, A.; Biris, A. S. Comparative Aging Study of Organic Solar Cells Utilizing Polyaniline and Pedot:Pss as Hole Transport Layers. ACS Appl. Mater. Interfaces 2015, 7, 27667−27675. (6) Bhadra, S.; Khastgir, D.; Singha, N. K.; Lee, J. H. Progress in Preparation, Processing and Applications of Polyaniline. Prog. Polym. Sci. 2009, 34, 783−810. (7) Goyenola, C.; Stafströ m, S.; Schmidt, S.; Hultman, L.; Gueorguiev, G. K. Carbon Fluoride, Cfx: Structural Diversity as Predicted by First Principles. J. Phys. Chem. C 2014, 118, 6514−6521. (8) Geethalakshmi, D.; Muthukumarasamy, N.; Balasundaraprabhu, R. Csa-Doped Pani/Tio2 Hybrid Bhj Solar Cells − Material Synthesize and Device Fabrication. Mater. Sci. Semicond. Process. 2016, 51, 71−80. (9) Chen, X. P.; Jiang, J. K.; Liang, Q. H.; Yang, N.; Ye, H. Y.; Cai, M.; Shen, L.; Yang, D. G.; Ren, T. L. First-Principles Study of the Effect of Functional Groups on Polyaniline Backbone. Sci. Rep. 2015, 5, 16907. (10) Chen, X.; Sun, S.; Wang, X.; Li, F.; Xia, D. Dft Study of Polyaniline and Metal Composites as Nonprecious Metal Catalysts for Oxygen Reduction in Fuel Cells. J. Phys. Chem. C 2012, 116, 22737− 22742. (11) Duan, Y.; Liu, J.; Zhang, Y.; Wang, T. First-Principles Calculations of Graphene-Based Polyaniline Nano-Hybrids for Insight of Electromagnetic Properties and Electronic Structures. RSC Adv. 2016, 6, 73915−73923. (12) Ullah, H.; Shah, A.-u.-H. A.; Bilal, S.; Ayub, K. Dft Study of Polyaniline Nh3, Co2, and Co Gas Sensors: Comparison with Recent Experimental Data. J. Phys. Chem. C 2013, 117, 23701−23711. (13) Almeida de Mendonça, J. P.; et al. Structural and Vibrational Study of Graphene Oxide Via Coronene Based Models: Theoretical and Experimental Results. Mater. Res. Express 2016, 3, 055020. (14) Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864−B871. (15) Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133− A1138. (16) 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 B.01; Gaussian, Inc.: Wallingford, CT, 2009. (17) Becke, A. D. Density-Functional Thermochemistry. Iii. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648−5652. (18) 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: Condens. Matter Mater. Phys. 1988, 37, 785−789. (19) Stephens, P. J.; Devlin, F. J.; Chabalowski, C. F.; Frisch, M. J. Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields. J. Phys. Chem. 1994, 98, 11623−11627. (20) Vosko, S. H.; Wilk, L.; Nusair, M. Accurate Spin-Dependent Electron Liquid Correlation Energies for Local Spin Density Calculations: A Critical Analysis. Can. J. Phys. 1980, 58, 1200−1211. (21) Soler, J. M. The Siesta Method for Ab Initio Order- N Materials Simulation. J. Phys.: Condens. Matter 2002, 14, 2745. (22) Ordejón, P. Order-N Tight-Binding Methods for ElectronicStructure and Molecular Dynamics. Comput. Mater. Sci. 1998, 12, 157−191. (23) Dion, M.; Rydberg, H.; Schröder, E.; Langreth, D. C.; Lundqvist, B. I. Van Der Waals Density Functional for General Geometries. Phys. Rev. Lett. 2004, 92, 246401. (24) Lee, K.; Murray, É. D.; Kong, L.; Lundqvist, B. I.; Langreth, D. C. Higher-Accuracy Van Der Waals Density Functional. Phys. Rev. B: Condens. Matter Mater. Phys. 2010, 82, 081101.

concentration of GQD−X on PANI. But this constitutes with our methodolog, an important computational demand. In the Supporting Information, we present the results for PANI/ GQD−X (X = COOHOH) for ratios 2:1 and 3:1. What is observed is that the presence of more PANI chains generates a red shift of the π* → π* transitions when compared to the 1:1 ratio, thus contributing to the near-infrared region of the optical spectra. It is clear that the concentration effect plays a role, but what is still valid for all the studied systems is that the transitions with better oscillation strength involves π* → π* from PANI to GQD−X.

4. CONCLUSIONS We have presented a mixed methodology utilizing MD-DFT simulations followed by very accurate TD-DFT calculations. We demonstrated that the electronic interaction between PANI polymer and GQDs is enhanced when oxidizing groups are incorporated in the GQD, which further introduces structural asymmetries, which is very relevant for the optical properties. In particular, the presence of OH groups at the edges of the GQD−OH and GQD−COOHOH reduces the symmetry of the electronic topography, favoring the charge separation upon an optical excitation, thus constituting a very promising system for optical applications in general and organic solar cells in particular. The main conclusion of the work establishes the potential use of nanocomposites PANI/GQDs for the tuning of the optical applications.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.7b03604. Additional information, tables, and figures (PDF)



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected]. *E-mail: [email protected]. *E-mail: [email protected]. ORCID

Ricardo Faccio: 0000-0003-1650-7677 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors wish to thank the Uruguayan funding institutions CSIC, ANII and PEDECIBA.



REFERENCES

(1) Sreejith, S.; Hansen, R.; Joshi, H.; Kutty, R. G.; Liu, Z.; Zheng, L.; Yang, J.; Zhao, Y. Quantum Dot Decorated Aligned Carbon Nanotube Bundles for a Performance Enhanced Photoswitch. Nanoscale 2016, 8, 8547−8552. (2) Mombrú, D.; Romero, M.; Faccio, R.; Mombrú, Á . W. Tuning Electrical Transport Mechanism of Polyaniline−Graphene Oxide Quantum Dots Nanocomposites for Potential Electronic Device Applications. J. Phys. Chem. C 2016, 120, 25117−25123. (3) Mombrú, D.; Romero, M.; Faccio, R.; Mombrú, Á . W. Effect of Graphene-Oxide on the Microstructure and Charge Carrier Transport of Polyaniline Nanocomposites under Low Applied Electric Fields. J. Appl. Phys. 2017, 121, 045109. 16582

DOI: 10.1021/acs.jpcc.7b03604 J. Phys. Chem. C 2017, 121, 16576−16583

Article

The Journal of Physical Chemistry C (25) Bauernschmitt, R.; Ahlrichs, R. Treatment of Electronic Excitations within the Adiabatic Approximation of Time Dependent Density Functional Theory. Chem. Phys. Lett. 1996, 256, 454−464. (26) Runge, E.; Gross, E. K. U. Density-Functional Theory for TimeDependent Systems. Phys. Rev. Lett. 1984, 52, 997−1000. (27) Faccio, R.; Pardo, H.; Mombrú, A. W. Hybrid Compounds Based on Fullerene and Polycyclic Aromatic Hydrocarbons with Absorption in the near Infrared Region. Comput. Theor. Chem. 2013, 1018, 50−58.

16583

DOI: 10.1021/acs.jpcc.7b03604 J. Phys. Chem. C 2017, 121, 16576−16583