Molecular Design and Geometry Structuresand Stability for PyrroleSubstitutes, DFT study as Organic Solar cell system (One Anchoring System)
The theoretical electron properties ofPyrroleSubstitutes( S1, S2, S3, and S4 ) were carried out by using quantum chemical calculations. The optimized structures were obtained by the Density Functional Theory DFT/B3LYP level of theory using the basis set DVZ(d). The optimized structures of compounds have the global minimum energy. It was found that the dipole moment of compound (S2) have high values compared with the Compounds (S1,S3,S4). Global descriptors such as the MO energies of HOMO, LUMO levels and ”E, were determined and used to identify the differences in the stability and reactivity of compounds. In general, the calculated values lead to the conclusion that on the one hand the stability of the compounds are S4> S3> S2> S1. On the other hand, the theoreticalstudyofnovelacceptor-donor organicmaterialsbasedon these compounds hasbeeninvestigated.Differentelectronsidegroupswereintroducedtoinvestigate theireffectsonthe electronicstructure;HOMO,LUMO andenergy gap. Thestructural and electronic study as shown in this paper in hand forthesecompoundscouldhelptodesignmoreefficientfunctionalphotovoltaic organic materials. the adsorption energy between the compounds and the TiO2 surface was studied. The most suitable adsorption configurations high value of adsorption energy was compound S2.
Alkyl Chain Substituted, PyrroleSubstitutes, Structural and Electronic Properties, DFT/B3LYP, Band Gap
Aromatic heterocyclic compounds based on pyrrole or pyrdinium core including carboxylate derivatives, Schiff-bases and their metal complexes, are materials of much interest to researcher due to their wide applications in medicine[1-2], catalysis, luminescence, optical devices and clean energy[5-6]. The”-conjugatedmoleculesbasedonheterocyclicringandcharacterizedbyanextended ”- conjugationisofbroadinterestinmanyareasofresearch,fortheirapplicationsinorganicelectronics, particularly, inthelastdecade,optoelectronicdevicetechnology[7,8],suchaslight-emitting diode(LEDS), thinfilmtransistors(TFTs)andlow-costsolarcellsResearch in organic solar cells has attracted considerable interests in the past decade because of their potential applications as alternative light harvesting devices other than the conventional silicon-based solar cells [12-13]. However, the performance and lifetime most organic solar cells there are not satisfactory. In general, organic materials exhibit low carrier mobilities, low absorbance in near-IR region, short exciton diffusion length, and poor long-term stability. In order to improve the performance of organic photovoltaic devices, design of new photosensitizing/charge transport materials and fabrication of devices with improved architecture are essential. So,designing andsynthesizing molecules involves connecting donor (D), ”-bridge and acceptor (A) groups as D””””A system to create highly polarized withmore interestingpropertiesplayacrucialroleintechnology. Itisimportant tounderstandthenatureoftherelationshipbetweenthemolecularstructureandtheelectronicpropertiesto provideguidelinesfor thedevelopmentofnewmaterials.Nowadays much interest is devoted to the prediction and estimation of physicochemical properties of molecules, and materials by using computational. Computational chemistry (also known as molecular modeling ) is the application of computer-based models to the simulation of chemical is the processes and the computation of chemical properties and to use this for predication and understanding how electrons, atoms, and molecules interact . Density Functional Theory (DFT) is one major method of computational chemistry which has been accepted by the ab initio quantum chemistry community as a cost effective general procedure for studying physical properties of the molecules. This is because it is based on total electron density rather than on wave functions[14-19]. The work in hand attempts to study and predication of structure and electronic properties, relative stabilitiesof structural modifications were investigated for the D””””A system with an blend (alkyl chain substituted by C, O, N, S ) as donor and pyrrole as ” ”’spacer (” ”’conjugated) while carboxylic ”’ as acceptor ( an anchoring group) as imitation for organic solar cell system [20-23]. As strategy for designing solar cell analogues with reduced band gaps and whether to be materials for organic solar cell, by performing Density Functional Theory (DFT)/B3LYB level of theory using the basis set DZV(d).
Modeling and Geometry Optimization
The quantum chemical calculations were performed for 4 compounds using the PC Gamess (Firefly) .program . optimal geometry of all compounds under study computed using density functional theory with double zeta valence DZV(d) basis set and Becke’s 3-parameter exchange functional with Lee-Yang-Parr correlation energy functional(B3LYP) [25-26]. All calculations were performed on the Pentium (R)4/IPM-PC- CPU 3.00GHz, 2.00GB . The structures of compounds are shown below in Figure 1:
Figure- 1: Molecular Structure of Compounds
Results and Discussion
The geometry optimized structures ( S1, S2, S3 and S4 ) are visualized in Figure (2) below while the selected parameters of their structural data are summarized in Table 1. For all molecules,pyrrolesubstitutes (( S1, S2, S3 and S4)), optimized geometric structures, dipolemoments,total energies, the HOMO, LUMO and energy band gap (”E) calculation have been investigated after total optimization by B3LYP/DZV.
Total Energy: The structure optimization (minimum energy) obtained by the DFT(B3LYP) presented in Table 1.Thetotalenergy determinestheoccurrenceor non-occurrenceof chemicalreactionsandstereospecific paths in intra- and intermolecular processes. The total energy ofthe system composedof the internal,potential, andkineticenergy. the total energy (absolute values) for compounds are S4> S3> S2> S1.
DipoleMoment: The dipole moment (” in Debye) is another important electronic parameter used to describe the polarity of the molecule. Thisparameterthathelpsintheunderstanding of interaction between atoms in thesame ordifferent molecules dipolemomentincreases with the increaseinelectronegativity ofatoms. Also is related to the distribution of electrons in a molecule. The dipole moment (” in Debye) for compounds are S2 > S4> S1> S3. Chemical reactivityusuallyincreaseswiththe increasein dipolemoment as well as attractive for the interaction with other systems and to form complexes, from Table 1 that, the dipole moment has maximum values for compound (S2) compared with compounds (S1, S3, S4). The high dipole moment may make the compound (S2) attract other systems to interact, to form complexes, and to indicate highly polar molecules.[28-29].
Figure- 1: DFT/B3LYP ”’DZV(d) Calculated Optimized Structures of the Possible Predication for
the Compounds in Gas Phase.
HOMO-LUMOEnergyGap:The electron distribution for the highest occupied molecular orbital (HOMO), the lowest unoccupied molecular orbital (LUMO) and HOMO and LUMO energy gaps (”E) for compounds (S1, S2, S3, and S4) are depicted in figure3.The egien values of LUMO and HOMO and their energy gap reflect the chemical activity of the molecule. LUMO as an electron acceptor represents the ability to obtain an electron, while HOMO as an electron donor represents the ability to donate an electron. The smaller the LUMO, HOMO and energy gaps are the easier is it for the HOMO electrons to be excited; the higher the HOMO energies, the easier it is for HOMO to donate electrons; the lower the LUMO energies are the easier is it for LUMO to accept electrons. So from the energies band gaps it results that the stabilities of compounds are (S4> S3> S2> S1). Compounds (R1&R2) with a small HOMO-LUMO gap and A large HOMO”’LUMO gap implies a high kinetic stability and low chemical reactivity, because it is energetically unfavorable to add electrons to a high-lying LUMO or to extract electrons from a low-lying HOMO. Generally it is said that low Egap leads to red shift [30-32].
Table- 1:TheTotal Energy,Dipole moment(”), EHOMO,ELUMO,Energy Gap(eV), of (S1,S2,S3 and S4) Obtained by B3LYP/ DZV(d).
.Compound T.E(a.u) EHOMO eV ELUMO eV D.M ”E
S1 -713.34501 -5.578378 -0.307491 1.622811 5.27088
S2 -729.37607 -5.657291 -0.329260 2.074038 5.32803
S3 -749.25059 -5.502185 -0.321096 1.576536 5.18108
S3 -1072.2131 -5.382454 -0.394568 2.011611 4.98788
Here, we studied the photovoltaic properties of the all studied compounds as donor blended with TiO which is used as an acceptor in solar cell devices. To elucidate the parameters that influence the photovoltaic efficiency toward better understanding of the structure”’property relationships. The presented study of structural, electronic and optical properties for these compounds could help to design more efficient functional photovoltaic organic materials. The conjugated linker in a D””””A system acts as both a component for light harvesting and a channel for charge transport. A good conjugated linker should promote the absorption of light over a wide region, and at the same time, facilitate charge transfer. The pyrrole unit is an appealing class of conjugated linker used in photovoltaic cells. It provides effective conjugation without affecting of the stability. In this model in hand present itself as D””””A system with an blend (alkyl chain substituted by C, O, N and S ) as donor and pyrrole as ” ”’spacer (” ”’conjugated) while carboxylic ”’ as acceptor ( an anchoring group)as well as TiO2 as surface of the conduction band, as application of intermolecular charge transfer (ICT). The ICT property of a D””””A dye is strongly dependent on the electron-donating ability of D, the electron-withdrawing ability of A, as well as the electronic characteristic of the ”-conjugated bridge. It can be tuned through chemical modification of each component. The HOMO and the LUMO energy levels of the donor and acceptor components are very important factors to determine whether effective charge transfer will happen between donor and acceptor. In these compounds, the HOMO is mainly localized on the pyrrole and electron donating ligands, and the LUMO is mainly distributed over the carboxylic (Figure 4). The distribution of LUMO and HOMO is favorable for e”’cient electron injection into TiO2 via chemical bonding. Our strategy for designing solar cell analogues with reduced band gaps was based upon the observation that the HOMO in all compounds has slight change from compound S1 to S4, with high HOMO on the sulfur heteroatom, whilst the LUMO has minor slight changed on contribution on the heteroatom. We thereby reasoned that replacement of sulfur by another atom (C, O, N) have a good influence on the HOMO energy level whichwould primarily result from the small difference in electronegativity between sulfur and another atom (C, O and N) .
The band gap of S4 is much smaller than that of the other substituted compounds. This may be attributed to the number of electron attracting sulfur atom which put in side groups, and also the resonance effect tends to con”’ne the ”-electrons within the pyrrole ring and hence to prevent their delocalization along the whole conjugated chain and leads to the destabilisation of both the HOMO level and LUMO levels with decrease in the energy gap, another concepts to the small energy gap is the relative distance between the donor alkyl chain substituted, (”-pyrrole) and (carboxylic )acceptor which responsible for the strongest localized contribution to the population transfer is the LUMO orbital of the donor[33-35]. Table1, showthechange oftheelectron-donoragreateffectontheHOMOandLUMOlevels.
Table- 2:TheVoc (ev) and ” gap of (S1,S2,S3 and S4)
PCBM C60 PCBM C70 PCBM C78-D3
Compound Voc ” Voc ” Voc ”
S1 3.0925 3.3925 2.93250 3.2325 3.3925 3.6925
S2 3.0707 3.3707 2.91074 3.2107 3.3707 3.6707
S3 3.0789 3.3789 2.91890 3.2189 3.3789 3.6789
S3 3.0054 3.3054 2.84543 3.1454 3.3054 3.6054
Open-Circuit Voltage (Voc).
Generally energy gap more than 0.2 eV between the LUMO of the organicmolecules and the conduction band of the TiO2 is necessary for effective electron injection .As well as theexperimentphenomenon wasquiteconsistentwithprevious literature, whichreported thattheincreaseof theHOMOlevelsmay suggestanegativeeffectonorganicsolarcellperformanceduetothebroadergapbetweentheHOMOlevelof theorganicmoleculesand theLUMOlevelofPCBM (acceptor was reported to be between”’3.47 eV, PCBM C70 -3.54 eV, PCBM C78-D3 -4.0 eV As it’s known, the most efficient organic solar cells are based on the bulk heterojunction (BHJ) structure of the blend of organic material donors and fullerene derivative acceptors (PCBM) (opencircuitvoltage,Voc) is the maximum voltage that can be generated by the device.Although Voc is an important criterion with which we can estimate the efficiency of composite (donor/PCBM), we must take into account the difference between the LUMO of conjugated-” compound donor and those of (PCBM) acceptor noticed ” gap. This is an indicator for optimization of the open circuit voltage (Voc), which should be maximal for efficiency of solar cells.Maximizing the open-circuit voltage in a low-bandgap compound is one of the critical factors towards enabling high-ef”’ciency solar cells so that from table 1, the Maximizing the open-circuit voltage was the molecule S1 blended with C78-D3. which mean to increase of power conversion efficiency for solar cell. ThetheoreticalvaluesofVocwere calculatedfromthe followingexpression:
Themaximumopencircuitvoltage oftheBulkHeteroJunction solarcellisrelatedtothe differencebetweentheHOMOoftheelectrondonorandtheLUMOoftheelectronacceptor,takingintoaccountthe energylostduringthe photo-chargegeneration. knowingthatinorganicsolarcells,theVocisconsideredlinearlyasafunctionofthe HOMOlevelandthedonorleveloftheLUMOoftheacceptor. The obtained values of Voc of the studied molecules calculated according to the equation (1) (table 2) these values are sufficient for a possible efficient electron injection. Therefore, all the studied can be used as sensitizers because electron injection process from the exited molecule to the condition band of PCBM and subsequent regeneration is possible in sensitized solar cells. It can also be found that, the HOMO and LUMO energies of the studied compounds are slightly different. This implies that different structures play key roles on electronic properties and the effect of slight structural variations, especially the effect of the motifs branched to the molecule on the HOMO and LUMO energies is clearly seen[39-40].
On the other hand the molecules of one anchoring system S1, S2, S3 and S4 (-CH2-CO-OH-terminated) can be desorption on surface TiO2 (layers) by mono anchoring and this lead to increase a sufficient crosslinking or coordination of the carboxylic groups molecules with TiO2 surface to formation one substrate-CH2-CO-O-TiO2 surface and more stability. From energy gap value table 1, it can see the compound S4 have small energy gap more stability more coordination of the carboxylic groups molecules with TiO2 surface and this also obviously from geometry optimization. The variation of Voc of the compound understudy as a function of HOMO level is shown in figure 3. It can be seen that in general Voc of compounds varies inversely with HOMO level of the electron donor material.
Figure- 3: The Open-circuit voltage distribution versus the HOMO energy levels
On the other hand the bond length between the Donor moieties and ”-spacer is very important for loading and transfer charge between the donor and acceptor group in the device solar cells. The bond lengths and angles for compounds are listed in Table 2; As shown, there are slight changes in the bond lengths and angles of compounds, like the changes of the bond length and angles in the ring indicate the presence of -conjugation which causes of the electrostatic attraction between atoms. As well as when the replacing the carbon atom by atoms oxygen and nitrogen in the bond (1-2), this lead to short the bond while when replacing carbon by silver lead to increase in the bond distance ( be week), also this happen to the bond (1-18).Also the effect of substitution on the ring leads to redistribution of electron cloud in the ring of the compounds (S1, S2, S3, and S4 )[42-45]. These results show that the connection between donor moieties (D) via ”-bridge (pyrrole) is crucial which enhance the (ICT) character .
Figure- 3: Electronic Distribution in HOMO and LUMO Orbitals.
we study the relationship between adsorption geometry and electron injection properties of one of the most successful metal-organic, on the TiO2(111) surface(the most commonly used DSSC semiconductor). We systematically construct all possible adsorption configurations of the (S1, S2, S3, and S4 ) molecule (bidentate) on this surface by connected to TiO2 through the carboxylic acid.. we find that a large number of adsorption configurations are possible ”’ more than ten structures, with one( anchoring group carboxylic groups (bidentate )adsorbed on the TiO2 (111) surface.
In this computational work the Monte Carlo simulation techniques is used to find the preferential adsorption sites on the TiO2 surface by finding the low-energy adsorption sites on the TiO2 surface(111) Molecular mechanics (force field) tools are used to investigate the simulated adsorption system.
Geometry optimization was achieved using COMPASS force ”’eld and the Smart minimize method by high-convergence criteria. This was followed by modelling the molecular electronic structures, including the distribution of frontier molecular orbitals and Fukui indices to establish the active sites as well as thelocal reactivity of the molecules. This computational study aims to find low-energy adsorption sites to investigate the preferential adsorption of molecule on TiO2 surface aiming to find a relation between the effect of its molecular structure and its compounds efficiency.
Table 2: Selected Structural Parameters of The Optimized Compounds, Bond Distance(A”) and Bond Angles(”) Obtained By B3LYP/DZV(d).
Bond/Angle S1 S2 S3n S4
(1-2) 1.533 1.430 1.431 1.830
(1-18) 1.494 1.368 1.431 1.762
(18-19) 1.380 1.431 1.469 1.428
(18-21) 1.425 1.385 1.431 1.383
(19-20) 1.386 1.381 1.383 1.384
(20-22) 1.381 1.378 1.382 1.378
(22-21) 1.381 1.387 1.381 1.456
(2-1-18) 111.5 113.7 113.3 101.0
(1-18-19) 125.8 127.5 128.3 127.5
(1-18-21) 127.0 124.3 124.6 125.2
Figure 3. Most suitable configuration for adsorption in bridged bidentate(one anchoring groups) of compounds on TiO2 surface substrate obtained.
As seen in Table 3, The parameters presented in including total energy, in kcal mol”’1, of the substrate”’adsorbate configuration. The total energy is defined as the sum of the energies of the adsorbate components, the rigid adsorption energy and the deformation energy. In this study, the substrate energy (TiO2 surface) is taken as zero. In addition, adsorption energy in kcal mol”’1, reports energy released when the relaxed adsorbate on the substrate. The adsorption energy is defined as the sum of the rigid adsorption energy and the deformation energy for the adsorbate 2 components. The rigid adsorption energy reports the energy, in kcal mol”’1, released when the unrelaxedadsorbate components (i.e., before the geometry optimization step) are adsorbed on the substrate. The deformation energy reports the energy, in kcal mol , released when the adsorbed adsorbate components are relaxed on the substrate surface. Table3 shows also (dE /dNi), which reports the energy, in kcal mol”’1 ads , of substrate”’adsorbate configurations where one of the adsorbate components has been removed[46-50].
Table 3: The output and descriptors calculated by the Monte Carlo simulation of compoundsconfirmations on TiO2 (111) surface.
Co. Total Energy
Adsorption Energy Rigid Adsorption Energy Deformation
S1 -53.64188 -2.5181e+03 -46.4778 -2.4719e+03 -2.5183e+03 2.319
S2 -48.70808 -2.5772e+03 -54.5584 -2.5227e+03 -2.5772e+03 3.637
S3 -58.88136 -2.5555e+03 -65.2632 -2.4905e+03 -2.5558e+03 3.994
S4 -64.74342 -2.5373e+03 -64.2526 -2.4730e+03 -2.5373e+03
As seen in table 3, the adsorption energy. Figure3, Table 3 shownthe most suitable adsorption configurations for the oxygen compounds on TiO2 . It can be seen from Table 3 that the values of adsorption energy for compounds are S2 > S3> S4> S1, the compounds and TiO2 surface are determined by the van der waals forces while the electrostatic energies are as the sum of the energies of the adsorbate components and the adsorption energy. The van der Waals interactions are dominant in adsorption of lone-pair electrons of oxygen on metal surfaces. Indicate that the adsorption energies between the oxygen in carboxylic functional group in compounds and TiO2 surface are determined by the van der Waals forces. The electronegativity and the tow lone-pair of electron in oxygen lead to gethigh value of adsorption energy for compound S2 (containing oxygen group) is related to the distribution of electrons in a molecule and indicate that the oxygen compounds are likely to adsorb on the TiO2 surface and form stable films to design organic solar cell system.
In this study, the quantum chemical investigation of the geometries and electronic properties of pyrrolesubstitutes(S1, S2, S3, S4 and S5), optimized geometric structures, dipolemoments,total energies, the HOMO, LUMO and also the energy band gap (”E) calculation have been investigated after total optimization by B3LYP/DZV(d). This study is a theoretical analysis of the geometries and optoelectronic properties which displays the effect of substituted groups on the structural and optoelectronic properties of these materials and leads to the possibility to suggest these materials for organic solar cells application. The concluding remarks are:
1-The results of the optimized structures for all studied compounds show that they have similar conformations. We found that the modification of several groups does change the geometric parameters.
(bond lengths and angles)of compounds, indicate the presence of -conjugation which causes of the electrostatic attraction between atoms. And this leads to redistribution of electron cloud in the ring of the compounds (S1, S2, S3, and S4 )
2- The calculated frontier orbital energy HOMO and LUMO and energy gap showed that the energy gap of the studied molecules differ slightly from 5.27088 to 4.98788 depending on the different structures. Depend on the energies band gaps, the stabilities of compounds are (S4>S3> S2> S1). Compounds (S4) with a small HOMO-LUMO gap and A large HOMO”’LUMO gap implies a high kinetic stability and low chemical reactivity.
3-The donor moieties improve the electronic properties of the studied molecules by reducing the energy gap, making them more conductive systems. also improves the charge transfer along the molecular chain. So that the best reducing of band gap (favorable state) was compound S4.
4- The best values of Voc are indicated for the studied compounds blended with C60, C70 or C78-D3,and higher value are given for molecule S1 blended with C78-D3. which mean to increase of power conversion efficiency for solar cell.
5-Values of Voc of the studied molecules calculated range from 2.842 eV to 3.392 eV these values are sufficient for a possible efficient electron injection. Therefore, all the studied can be used as sensitizers because electron injection process.
This calculation procedure can be used as a model system for understanding the relationships between electronic properties and molecular structure and also can be employed to explore their suitability in electroluminescent devices and in related applications. Finally, the procedures of theoretical calculations can be employed to predict the electronic properties on the other compounds, and further to design novel materials for organic solar cell.
6- Monte Carlo simulation studies help to find the most stable conformation and adsorption sites for a broad range of materials, The results indicated that the studied compounds could adsorb onto TiO2 surface by the van der Waals force, and The high value of adsorption energy for compound S2.
Weissleder.R, and Pittet .M. J., Nature, 452, 580, 2008.
Jun.Y.w, Lee .J. H. and Cheon J. W., Angew. Chem., Int. Ed., 47, 5122, 2008.
Melnik.M, Koman. M., Hudecov. D., Moncol. J., Dudov .B., Glowiak .T., Mrozinski. J., and
Holloway .C. E., Inorg. Chem. Acta, 308, 2000.
Smith C. S. and Mann K. R., Chem. Mater., 21, 5042”’5049, 2009.
Kanoo.P, Sambhu .R., and Maji.T, Inorg. Chem., 50, 400 402, 2011.
Khurana. J.M., ”ORGANIC CHEMISTRY, Chemistry of Heterocyclic Compounds”University of
Delhi Delhi -110007.2006.
Nguyen. VC,and Potje-Kamloth. K,ThinSolidFilms,338,142”’148.1999.
Gill.RE.,Malliaras. GG,Wildeman.J, and Hadziioannou.G,Adv.Mater, 6,132-135.1994.
Wang GM.,Qian. SX ,Xu. JH.,Wang. WJ,X Liu. X.,Lu. XZ.,and Li FM.,
 B.C. Thompson, J.M.J. Frechet, Angew. Chem., Int. Ed. 47, 58.(2008).
 S. G”nes, H. Neugebauer, N.S. Sariciftci, Chem. Rev. 107 , 1324.(2007).
 Miller, R. D., Lee, V. Y. and Moylan, C. R., Chem. Mater, 6, 1023-1032 (1994).
 Breitung, E. M., Shu, C.-F. and McMahon, R. J., J. Am. Chem. Soc., 122, 1154-1160.(2000).
 Facchetti, A., Beverina, L., van der Boom, M. E., Dutta, E. G., Pagani, G. A. and Marks, T. J., J.
Am. Chem. Soc., 128,2142-2153 (2006).
 Davies, J. A., Elangovan, A., Sullivan, P. A., Olbricht, B. C., Bale, D. H., Ewy, T. R., Isborn, C. M.,
Eichinger,B. E., Robinson, B. H., Reid, P. J., Li, X. and Dalton, L. R., J. Am. Chem. Soc., 130,
10565- 10575 (2008) .
 Abbotto, A., Beverina, L., Manfredi, N., Pagani, G. A., Archetti, G., Kuball, H. G., Wittenburg, C.,
Heck, J. and Holtmann, J., Chem.-A Eur. J., 15, 6175-6185 (2009) .
 Pi”etrickHudhomme., EPJ Photovoltaics 4, 40401 (2013)., www.epj-pv.org.,
 Sandeep B. Mane and Chen-Hsiung Hung., J. Name., 00, 1”’ 3 | 1-13., 2013., DOI:
Corneliu I. Oprea , PetrePanait , FanicaCimpoesu , MarilenaFerbinteanu and Mihai A. G”r”u.,
Materials., 6, 2372-2392, 2013., doi:10.3390/ma6062372
 S. Jungsuttiwong, R. Tarsang, S. Pansay, T. Yakhantip, V. Promarak, T. Sudyoadsuk, T.Kaewin, S.
Saengsuwan, and S. Namuangrak., International Scholarly and Scientific Research & Innovation
Junfeng .Y and Feng Z., J. Mater. Chem., 21, 9406., 2011. DOI: 10.1039/c1jm10274e.
AA Granovsky. PC Gamess, version 7.1.F (Firefly), Moscow State University, Moscow, Russia,
 Lee CT, Yang WT, Parr RG. Physical Review B.;37:785”’789.1988.
Becke AD. Journal of Chemical Physics. 1993;98:1372”’1377.1993.
 Fukui, K.Science.218,747.1982.
Ebenso EE; Isabirye DA and Eddy NO, Int. J. Mol. Sci.,11, 2473-2498.2010.
Obot. I. B,and Johnson. A,S. Comput. Chem. 43,6658.2012.
 Socrates, G.,andWiktor, Z. J PhysChem, 123A, 12690-12697.2009.
 Ortiz, R, P., Osuna, R, M., Delgado, M., Casado, J., Jenekhe, S, A., Hernandez, V., and Navarrete,
L. J. T. Int. J. Quantum Chem. 104, 635-644.2005.
 Abdullah G. Al-Sehemi, Ahmad Irfan, Abdullah M. Asiri, Yousry Ahmed Ammar., Bull. Chem.
Soc. Ethiop, 29(1), 137-148.2015. DOI: http://dx.doi.org/10.4314/bcse.v29i1.13
 Krishna .F, Warwick J. B, Christopher J. F and Paul C. D., Int. J. Mol. Sci. 13,
17019-17047, 2012: doi:10.3390/ijms131217019
Peng Qin., Ph.D. Thesis, KTH Chemical Science and Engineering, RoyalInstitute of Technology,
SE-100 44 Stockholm, Sweden.2010.
Wenjie Fan and Wei”’qiao Deng.,Commun. Comput. Chem.1, No.2, pp. 152”’170.,2013. doi: 10.
 Angelis, F.D.; Fntacci, S.; Selloni, A. Nanotechnology, 19, 424002.2008.
Sadiki A.Y, Bouzzine. S.M., Hamidi.M, Bejjit. L, Haddad.M, Serein-Spirau. F, L”re- Porte, J. P.,
Sotiropoulos. J. M. and Bouachrine.M, Journal of Computational Methods inMolecular Design,
4 (2):10-18, 2014.
 L. Zhang, Q. Zhang, H. Ren, H. Yan, J. Zhang, H. Zhang, J. Gu, J. Sol. Energ. Materials & Sol.
M. Bourass , A. TouimiBenjelloun , M. Benzakour , M. Mcharfi , M. Hamidi , S.M.Bouzzine, F.
Serein-Spirau , T. Jarrosson, J. P. L”re-Porte, J. M. Sotiropoulos ,M. Bouachrine.,J. Mater.
Environ. Sci.6 (6) 1542-1553.2015.
A. Adad, R. Hmammouchi, TaharLakhlifi and M. Bouachrine., J. Chem. Pharm. Res.,
 Jianhui H, Mi-Hyae P, Shaoqing Z, Yan Y, Li-Min C, Juo-Hao L, and Yang Y.,Macromolecules
, 41, 6012-6018, 2008.
Ram B. A, Gao-Fong C, , Manoj R. Z, Ching-Fa Y, Eric Wei-Guang D
and Chen-Hsiung H. Chem. Asian J. 2013, 8, 2144 ”’ 2153., DOI: 10.1002/asia.201300328
Cai-Rong Z, Li-Heng H, Jian-Wu Z, Neng-Zhi J, Yu-Lin S, Hai-Min Z, Yu-Hong C, Zi-Jiang L and
Ji-Jun G.,Computational and Theoretical Chemistry., 1039, 62-70,.2014.
Jesus Baldenebro-L, Jos” Castorena-G , Norma Flores-H, Jorge Almaral-S and Daniel Glossman-
M., Int. J. Mol. Sci. 2012, 13, 4418-4432; doi:10.3390/ijms13044418
Nurcan S. T, Fatma B, A. Sezai S., Journal of Molecular Structure: THEOCHEM 857,95”’
Enrico Ronca, MariachiaraPastore, and Filippo De Angelis., Leonardo Belpassi,., Francesco
Tarantelli.,Energy Environ. Sci., 6., 183.2013.
 P. Persson,., M. J. Lundqvist,R. Ernstorfer,W. A. Goddard III,and F. Willig., J. Chem. Theory
Comput., 2, 441-451, 2006.
 Fernando Mendizabal, Alfredo Lop”z, Ramiro Arratia-P”rez, Gerald Zapata-Torres.,
Computational and Theoretical Chemistry., 1070,117-125, 2015.
Corneliu I. Oprea , PetrePanait , FanicaCimpoesu , MarilenaFerbinteanu and Mihai A. G”r”u.,
Materials, 6, 2372-2392, 2013; doi:10.3390/ma6062372
Szymon G and Marek S.,Int. J. Mol. Sci.14, 2946-2966, 2013; doi:10.3390/ijms14022946.