Skip to main content
  • Research article
  • Open access
  • Published:

DFT and TD-DFT calculation of new thienopyrazine-based small molecules for organic solar cells



Novel six organic donor-π-acceptor molecules (D-π-A) used for Bulk Heterojunction organic solar cells (BHJ), based on thienopyrazine were studied by density functional theory (DFT) and time-dependent DFT (TD-DFT) approaches, to shed light on how the π-conjugation order influence the performance of the solar cells. The electron acceptor group was 2-cyanoacrylic for all compounds, whereas the electron donor unit was varied and the influence was investigated.


The TD-DFT method, combined with a hybrid exchange-correlation functional using the Coulomb-attenuating method (CAM-B3LYP) in conjunction with a polarizable continuum model of salvation (PCM) together with a 6-31G(d,p) basis set, was used to predict the excitation energies, the absorption and the emission spectra of all molecules.


The trend of the calculated HOMO–LUMO gaps nicely compares with the spectral data. In addition, the estimated values of the open-circuit photovoltage (Voc) for these compounds were presented in two cases/PC60BM and/PC71BM.


The study of structural, electronics and optical properties for these compounds could help to design more efficient functional photovoltaic organic materials.


The organic bulk heterojunction solar cells (BHJ) are considered as one of the promising alternative used for renewable energy. This is attributed to their several advantages to fabricate the flexible large-area devices and also to their low cost compared to other alternatives based on inorganic materials [1, 2]. Generally, the organic BHJ solar cells based on the mixture of electron donor (material organic) and electron acceptor materials as PCBM or its derivatives and have been utilized in the aim to harvest the sunlight. Over the past few years, considerable effort has been focused on improving organic solar cells (OSC) performance to achieve power conversion efficiencies (PCE) of 10%. The following strategies have been adopted for this purpose [313]: (1) design of the new photoactive materials able to increase the efficiency of photoconversion such as fullerenes and π-conjugated semiconducting polymers; (2) use of functional layers of buffering, charge transport, optical spacing, etc., and; (3) morphological tuning of photoactive films by post-annealing, solvent drying, or processing by using additives. After many efforts, the design of the organic BHJ solar cells based on polymer semiconducting (PSCs) as an electron donating and PCBM as an electron accepting showed impressive performances in converting solar energy to electrical energy. Finally, the power conversion efficiency (PCE) was improved in the range of 7–9.2% [1421] for single layer PSCs and 10.6% [14] for tandem structured PSCs. These kinds of solar cells based on polymers have potential applications in next-generation solar cells compared to dye-sensitized solar cells (DSSC) and inorganic thin-film. On the other hand, considerable research has been directed to developing an efficient small-molecule organic used as a semiconductors and to improve their performance in the organic solar cells (OSCs), with the near-term goal of achieving a PCE comparable to that of polymer solar cells (PSCs) [2224].

Small-molecule organic semi-conductors are more suitable than polymer-based ones for mass production because the latter suffer from poor reproducibility of the average molecular weight, high dispersity, and difficulties in purification. Recently, the small molecule for organic solar cells (SMOSCs) with PCEs exceeding 6% have been reported [25] thus making solution-processed SMOSCs strong competitors to PSCs. This inspires us to develop a new low band gap for small molecules for organic solar cells application. In order to achieve high current density in SMOSCs, utilizing new donor molecules that can efficiently absorb the sunlight at the maximum solar flux region (500–900 nm) of the solar spectrum, because the energy conversion efficiency of the small molecule for organic solar cells is directly attached to the light harvesting ability of the electron donor molecules. In addition, to get high open circuit voltage (Voc), the HOMO levels of the donor molecules should be down a −5.0 eV, in which this factor is calculated by the difference between the HOMO and LUMO levels of the donor and acceptor materials, respectively. The most small molecule organic semiconductors used in solar cells have a push–pull structure comprising electron donors and acceptors in objective to enhance the intramolecular charge transfer (ICT) and the band gap becomes narrow and then, yielding higher molar absorptivity [2225]. A common strategy to enhance the power conversion efficiency of low band gap conjugated molecules as an alternating (D-A) or (D-π-A) structures because this improves the excitation charge transfer and transport [26]. Different authors described in recent studies the importance of compounds with D-π-A structure and their role in the elaboration of the organic solar cell [2729]. The organic material based on thienopyrazine has been used as a donor unit; still receive considerable attention for their exceptional optoelectronic properties [30, 31]. Knowledge about the optoelectronic properties of these new materials can help with the design of new materials with optimized properties for solar energy conversion. In our previous works [32, 33], we have reported a theoretical study of photovoltaic properties on a series of D-π-A structures of thienopyrazine derivatives as photoactive components of organic BHJ solar cells.

In order to obtain materials with more predominant capability, the development of novel structures is now being undertaken following the molecular engineering guidelines, the theoretical studies on the electronic structures of these materials have been done in order to rationalization the properties of known ones and the prediction those of unknown ones [26]. As is known, the knowledge of the HOMO and LUMO levels of the materials is crucial in studying organic solar cells. The HOMO and LUMO energy levels of the donor and of the acceptor compounds present an important factor for photovoltaic devices which determine if the charge transfer will be happen between donor and acceptor. The thienopyrazine derivatives would be much more promising for developing the panchromatic materials for photovoltaic, and thus, provide much higher efficiencies if new absorption bands could be created in the visible light region.

In this paper, we report a strategy to control the band-gap and different optoelectronics properties by using the DFT method on a series of no symmetrical branched molecules based on thienopyrazine as a central core and cyanoacrylic acid as the end group connected with different π-conjugated groups Xi, as shown in Fig. 1. We think that the presented study for these compounds listed in Fig. 1 bout their structural, electronic and optical properties could help to design more efficient functional photovoltaic organic materials, for aim to find the best material which is used as a donor electron in BHJ device in the solar cell.

Fig. 1
figure 1

Chemical structure of study compounds Pi (i = 1–6)

Computational methods

All calculations were carried out using density functional theory (DFT) with B3LYP (Becke three-parameter Lee–Yang–Parr) exchange-correlation functional [34]. 6-31G(d,p) was used as a basis set for all atoms (C, N, H, O, S). Recently, Tretiak and Magyar [35] have demonstrated that the charge transfer states can be achieved in D-π-A structure a large fraction of HF exchange is used. A newly designed, functional, the long range Coulomb-attenuating method (CAM-B3LYP) considered long-range interactions by comprising 81% of B88 and 19% of HF exchange at short-range and 35% of B88 and 65% of HF exchange at long-range [36]. Furthermore, The CAM-B3LYP has been used especially in recent work and was demonstrated its ability to predict the excitation energies and the absorption spectra of the D-π-A molecules [3740]. Therefore, in this work, TD-CAM-B3LYP method has been used to simulate the vertical excitation energy and electronic absorption spectra. It is important to take into account the solvent effect on theoretical calculations when seeking to reproduce or predict the experimental spectra with a reasonable accuracy. Polarizable continuum model (PCM) [41] has emerged in the last two decades as the most effective tools to treat bulk solvent effects for both the ground and excited states. In this work, the integral equation formalism polarizable continuum model (IEF-PCM) [42, 43] was used to calculate the excitation energy. The oscillator strengths and excited state energies were investigated using TD-DFT calculations on the fully DFT optimized geometries.

By using HOMO and LUMO energy values for a molecule, chemical potential, electronegativity and chemical hardness can be calculated as follows [44]:

$$\mu = \left( {E_{HOMO} + E_{LUMO} } \right)\;/\; 2$$

Chemical potential

$$\eta = \left( {E_{LUMO} - E_{HOMO} } \right)\;/\; 2$$

(Chemical hardness),

$$\chi =-\;(E_{HOMO} + E_{LUMO} )\;/\;2$$


all calculations were performed using the Gaussian 09 package [45].

Results and discussion

Ground state geometry

The optimized structures of all molecules obtained with the B3LYP/6-31G(d,p) level, are presented in Fig. 2.

Fig. 2
figure 2

Optimized geometries obtained by B3LYP/6-31G(d,p) of the studied molecules

Figure 2 shows the definition of torsional angles Φ1 and Φ2 between D and π-spacer A and π-spacer respectively, intramolecular charge transfer (ICT) which is represented by the π-spacer and the bridge bonds between D and π-spacer and A and π-spacer were marked as LB1 and LB2 respectively, using compound [P1] as an example (see Fig. 2). Torsional angles Φ1 and Φ2 are the deviation from coplanarity of π-spacer with the donor and acceptor and the LB1 and LB2 are the bond lengths of π-spacer from the donor and acceptor. The torsional angles (Φ1 and Φ2), and bridge lengths (LB1 and LB2) are listed in Table 1.

Table 1 Optimized selected bond lengths and bond angles of the studied molecules obtained by B3LYP/6-31G(d,p) level [the unit of bond lengths is angstroms (Å), the bond angles and dihedral angles is degree (°)]

As shown in Table 1, all calculations have been done by using DFT/B3LYP/6-31G(d,p) level. The large torsional angle Φ1 of the compounds P1, P2, P3, P4, P5 and P6 suggest that strong steric hindrance exists between the donor and π-spacer.

For P2, the dihedral angles Φ1 formed between the donor group and π-spacer is 0.78°, indicating a smaller conjugation effect compared to the other compounds where the coplanarity can be observed, but this geometry of P2 allows inhibiting the formation of π-stacked aggregation efficiently. Furthermore, the dihedral angles Φ2 of all compounds is very small (2.77, 2.95, 2.85, 2.82, 2.84 and 2.76) wich indicates that the acceptor (cyanoacrylic unit) is coplanar with π-spacer (thiophene–thienopyrazine–thiophene). In the excited state (S1), we remark that the dihedral angles Φ1 for all compounds are significantly decreased in comparison with those in the ground state (S0), except P2 and P6, Φ1 is almost similar to that of the ground state. It indicates that the nature of the S1 state of the molecular skeleton of all compounds is different from the S0 state, and the complete coplanarity in S1 state triggers the fast transfer of the photo-induced electron from S0 to S1.

The shorter value from the length of bridge bonds between π-spacer and the donor (LB1) and in another side between π-spacer and acceptor (LB2) favored the ICT within the D-π-A molecules. However, in the ground state (S0) the calculated critical bond lengths LB1 and LB2 are in the range of 1.421–1.462 Å showing especially more C=C character, except the compound P6, which enhances the π-electron delocalization and thus decreases the LB of the studied compounds and then favors intramolecular charge transfer ICT. On the other hand, upon photoexcitation to the excited state (S1), the bond lengths and torsional angles for these compounds significantly decreased in comparison with those in the ground state (S0), especially the linkage between the π-spacer and the acceptor moiety (LB2). These results indicate that the connection of acceptor group (2-cyanoacrylic acid) and the π-bridge is crucial for highly enhanced ICT character, which is important for the absorption spectra red-shift.

Electronic properties

Among electronic applications of these materials is their use as organic solar cells, we note that theoretical knowledge of the HOMO and LUMO energy levels of the components is crucial in studying organic solar cells. The HOMO and LUMO energy levels of the donor and of the acceptor components for photovoltaic devices are very important factors to determine whether the effective charge transfer will happen between donor and acceptor. The experiment showed that the HOMO and LUMO energies were obtained from an empirical formula based on the onset of the oxidation and reduction peaks measured by cyclic voltammetry. But in the theory, the HOMO and LUMO energies can be calculated by DFT calculation. However, it is noticeable that solid-state packing effects are not included in the DFT calculations, which tend to affect the HOMO and LUMO energy levels in a thin film compared to an isolated molecule as considered in the calculations. Even if these calculated energy levels are not accurate, it is possible to use them to get information by comparing similar oligomers or polymers.

The calculated frontier orbitals HOMO, LUMO and band gaps by using B3LYP/6-31G(d,p) level of six compounds (P1, P2, P3, P4, P5and P6) are listed in Table 2. The values of HOMO/LUMO energies are −5.025/−3.057 eV for P1, −5.276/−3.293 eV for P2, −5.091/−3.099 eV for P3, −5.139/−3.124 eV for P4, −5.155/−3.140 eV for P5 and −3.140/−3.159 for P6 and corresponding values of energy gaps are 1.968 eV for P1, 1.983 eV for P2, 1.992 eV for P3, 2.015 eV for P4, 2.015 eV for P5 and 2.171 eV for P6. The calculated band gap Eg of the studied model compounds increases in the following order P1 < P2 < P3 < P4 = P5 < P6. The much lower Eg of P1, P2 and P3 compared to that of P6 indicates a significant effect of intramolecular charge transfer, which would make the absorption spectra red shifted. However, the Eg values of P1, P2 and P3 are smaller than that of P6. This is clearly due to the effect of the electron-donor unit which is strong of P1, P2, and P3 than that of other compounds. All molecules present low energy gap are expected to have the most outstanding photophysical properties especially P1.

Table 2 Calculated EHOMO, ELUMO levels, energy gap (Eg), dipole moment (ρ) and other quantum parameters chemical as electronegativity (χ), chemical potential (μ) and chemical hardness (η) values of the studied compounds obtained by B3LYP/6-31G(d,p) level

Quantum chemical parameters

Generally, the molecules having a large dipole moment, possesses a strong asymmetry in the distribution of electronic charge, therefore can be more reactive and be sensitive to change its electronic structure and its electronic properties under an external electric field. Through the Table 2, we can observe that the dipole moment (ρ) of compounds P1 and P4 are greater than others compounds, therefore we can say that these compound are more reactive that other compound, indeed, these compounds are more favorite to liberate the electrons to PCBM.

On another side, we note that the PCBM has the smallest value of the chemical potential (μ = −4.9) compared to six compounds (P1, P2, P3, P4, P5, and P6) (see Table 2), this is a tendency to view the electrons to escape from compound Pi has a high chemical potential to PCBM which has a small chemical potential, therefore PCBM behaves as an acceptor of electrons and others compounds Pi behave as a donor of electrons. For the electronegativity, we remark that the PCBM has a high value of electronegativity than other compounds (P1, P2, P3, P4, P5, and P6) (Table 2), thus the PCBM is the compound that is able to attract to him the electrons from others compounds. In another hand, we remark that the PCBM compound has a high value of chemical hardness (η) in comparison with other six compounds, this indicates that the PCBM is very difficult to liberate the electrons, while the other compounds are good candidates to give electrons to the PCBM (see Table 2).

Figure 3 shows the frontier molecular orbitals for all the Six compounds (computed at B3LYP/6-31G(d,p) level). The FMOs of all six models have analogous distribution characteristics. All HOMOs show the typical aromatic features with electron delocalization for the whole conjugated molecule and are mainly localized at the donor parts and conjugated spacer, whereas the LUMOs are concentrated on the π-spacer and at the acceptor moieties (cyano acrylic unit). In another hand, the HOMO possesses an anti-bonding character between the consecutive subunits, while the LUMO of all oligomers shows a bonding character between the two adjacent fragments, so the lowest lying singlet states are corresponding to the electronic transition of π–π* type. Therefore the photoexcited electron will be transferred from donor moiety (donor of an electron) to the acceptor group during the excitation process, which is of benefit to the injection of the photoexcited electrons to the LUMO of the semiconductor (PCBM). In another side, we remark that the acceptor group (–CCNCOOH) of all compound has a considerable contribution to the LUMOs which could lead to a strong electronic coupling with PCBM surface upon photoexcitation electron and thus improve the electron injection efficiency, and subsequently enhance the short-circuit current density Jsc.

Fig. 3
figure 3

The contour plots of HOMO and LUMO orbitals of the studied compounds Pi

Photovoltaic properties

Generally, the power conversion efficiency (PCE) is the most commonly used parameter to compare the performance of various solar cells, and to describe it for any compounds, some important parameters has been evaluated such as the short-circuit current density (JSC), the open circuit voltage (VOC), the fill factor (FF), and the incident photon to current efficiency (Pinc). The power conversion efficiency (PCE) was calculated according to the following Eq. (1):

$$PCE\, = \;\frac{{J_{SC} \;V_{OC} FF}}{{P_{inc} }}$$

where the JSC is estimated by the maximum current which flows in the device under illumination when no voltage is applied, in which dependent on the morphology of the device and on the lifetime and the mobility of the charge carriers [46].

The maximum open-circuit voltage (Voc) of the BHJ is determined by the difference between the HOMO of the donor (π-conjugated molecule) and the LUMO of the acceptor, taking into account the energy lost during the photo-charge generation [47, 48]. It has been found that the VOC is not very dependent on the work functions of the electrodes [49, 50].

The theoretical values of open-circuit voltage Voc of the BHJ solar cell have been calculated from the following expression [47, 48]:

$$V_{OC} \; = \;\left| {E_{HOMO}^{Donor} } \right|\; - \;\left| {E_{LUMO}^{Acceptor} } \right|\; - \;0.3$$

where the represents the elementary charge, and the value of 0.3 V is an empirical factor. Scharber et al. [48] proposed the Eq  (2) using −4.3 eV as LUMO energy for the PC71BM.

In addition, low LUMO of the π-conjugated compounds and a high LUMO of the acceptor of the electron (PC71BM, PC60BM) increase the value of VOC, which contributes a high efficiency of the solar cells [48, 50].

The theoretical values of the open circuit voltage Voc of the studied molecules range from 1.499 to 1.804 eV in the case of PC60BM and 0.425 to 0.73 eV in the case of PC71BM (Table 3), these values are sufficient for a possible efficient electron injection into LUMO of the acceptor.

Table 3 Energy values of ELUMO (eV), EHOMO (eV), Egap (eV) and the open circuit Voltage Voc (eV) and LUMOdonor−LUMOacceptorof the studied molecules obtained by B3LYP/6-31G(d,p) level

In other side the Table 3 and the Fig. 4 show that the differences (LD − LA) of LUMO energy levels between those new designed donors (P1, P2, P3, P4, P5 and P6) and the acceptor of PC60BM is larger than 0 eV except P2. The same remark in case PC71BM, the differences (LD − LA) energy is also larger than 0 eV, which ensures efficient electron transfer from the donor to the acceptor (PC60BM, PC71BM) except P2 in case PC60BM because is more lower to 0 eV. This makes the transfer of electron from this compound (P2) to LUMO of PC60BM very difficult (LUMO of P2 is located below to LUMO of PC60BM).

Fig. 4
figure 4

Sketch of B3LYP/6-31G(d,p) calculated energies of the HOMO, LUMO level of study molecules

Therefore, all the studied molecules can be used as BHJ because the electron injection process from the excited molecule to the conduction band of PCBM and the subsequent regeneration is possible in an organic sensitized solar cell.

It is possible to assess the ideal performance donor, according to the position of its [ELUMO (donor) − ELUMO (acceptor)] energy and its band gap (Fig. 5). Theoretically, a maximum energy conversion efficiency of about 10% could be achieved for CPOs [51, 52] an oligomer having a LUMO energy level between −3.8 and −4.0 eV and a band gap between 1.2 and 1.9 eV has a theoretical power conversion efficiency between 8 and 10%. In a tandem configuration, the combination of two polymers band gap of 1.8 eV and 1.5 or 1.5 and 1.2 eV in two active layers separated to increase the effectiveness of a complete device for achieving a conversion efficiency of energy theoretical about 15%. We note that the higher power conversion efficiency could be achieved for P2 is 4 and 3% for P3.

Fig. 5
figure 5

Calculated efficiency under AM1.5G illumination for single junction devices based on composites that consist of a donor with a variable band gap and LUMO level and an acceptor with a variable LUMO level [34]

Optical properties

To understand the electronic transitions from our compounds, the quantum calculation on electronic absorption spectra in the gaseous phase and solvent (chloroform) was performed using TD-DFT/CAM-B3LYP/6–31G(d, p) level. The calculated absorption wavelengths (ʎ max), oscillator strengths (ƒ) and vertical excitation energies (E) for gaseous phase and solvent (chloroform) were carried out and listed in Table 4. The spectra show a similar profile for all compounds which present a main intense band at higher energies from 548.16 to 591.46 nm for gas phase and 574.33 to 625.38 for chloroform solution and were assigned to the ICT transitions. From Table 4, we could find that as the donor group changing, the first vertical excitation energies (E) were changed in decreasing order in both phases (gaseous and solvated): P6 > D5 > P4 > P2 > P3 > P1 showing that there is a red shift when passing from P6 to P1. We remark that the transition which has the larger oscillator strength is the most probable transition from the ground state to an excited state of all transitions, corresponding to excitation from HOMO to LUMO of gas phase and chloroform solution, This electronic absorption corresponds to the transition from the molecular orbital HOMO to the LUMO excited state, is a π–π* transition. These results indicate that all molecules have only one band in the Visible region (λabs > 400 nm) (Fig. 6) and P1 could harvest more light at the longer-wavelength which is beneficial to further increase the photo-to-electric conversion efficiency of the corresponding solar cells. So the lowest lying transition can be tuned by the different π-spacer.

Table 4 Absorption spectra data obtained by TD-DFT methods for the title compounds at CAM-B3LYP/6-31G(d,p) optimized geometries in the gas phase and in solvent phase (chloroform)
Fig. 6
figure 6

Simulated UV–visible optical absorption spectra of the title compounds with the calculated data at the TD-DFT/CAM-B3LYP/6-31G(d,p) level in chloroform solvent

In order to study the emission photoluminescence properties of the studied compounds Pi (i = 1 to 6), the TDDFT/CAM-B3LYP method was applied to the geometry of the lowest singlet excited state optimized at the CAM-B3LYP/6–31 (d, p), and the theoretical emission calculations with the strongest oscillator are presented in Table 5. The emission spectra arising from the S1 state is assigned to π* → π and LUMO → HOMO transition character for all molecules. Through analyzing the transition configuration of the fluorescence, we found that the calculated fluorescence has been just the reverse processed of the lowest lying absorption. Moreover, the observed red-shifted emission of the photoluminescence (PL) spectra when passing from P1 to P6 is in reasonable agreement with the obtained results of absorption. We can also note that relatively high values of Stocks Shift (SS) are obtained from all compounds P1 (179.64 nm), P2 (176.64), P3 (181.49 nm), P4 (178.33 nm), P5 (177.26 nm) and P6 (152.68 nm) (Table 5), this indicate that the compounds which have a weak Stocks Shift present a minimal conformational reorganization between ground state and excited state. Indeed, this stops the intermolecular transfer charge and delaying the injection phenomenon from LUMO of the compounds to LUMO of PCBM. In fact, the Stokes shift, which is defined as the difference between the absorption and emission maximums (EVA–EVE), is usually related to the bandwidths of both absorption and emission bands [53].

Table 5 Emission spectra data obtained by TD-DFT methods for the title compounds at B3LYP/6–31G(d,p) optimized geometries in chloroform solvent

Excited state lifetimes

The radiative lifetimes (in au) have been computed for spontaneous emission using the Einstein transition probabilities according to the following formula [54]:

$$\tau ={C^{3} } {\bigg/}{2(E_{Flu} )^{2} {\text{f}}}$$

where (c) is the velocity of light, E Flu is the excitation energy, and ƒ is the oscillator strength (O.S.). The computed lifetimes (τ), for the title compounds are listed in Table 5. However, an increase in lifetimes of Pi will retard the charge recombination process and enhance the efficiency of the photovoltaics cells. So, long radiative lifetimes facilitate the electron transfer upon the photoexcited electron, from LUMO of electron-donor to LUMO of electron-acceptor, thus lead to high light-emitting efficiency. The radiative lifetimes of the study compounds are from 7.61 to 7.11 ns and increases in the following order P4 < P1 < P2 < P5 < P3 < P6. This result is sufficient to obtain a high light-emitting efficiency, especially for P6.


We have used the density functional theory method to investigate the geometries and electronic properties of some thienopyrazine-derivatives in alternate donor-π-acceptor structure. The modification of chemical structures can greatly modulate and improve the electronic and optical properties of pristine studied materials. The electronic properties of new conjugated materials based on thienopyrazine and heterocyclic compounds and different acceptor moieties have been computed by using 6-31G(d,p) basis set at a density functional B3LYP level, in order to guide the synthesis of novel materials with specific electronic properties. The concluding remarks are:

The predicted band gaps by using DFT-B3LYP/6-31G(d,p) are in the range of 1.968–2.171 eV, knowing that the small band gap due to the increasing of the displacement of the electron between donor and acceptor spacer is very easy. The much lower Eg of P1, P2, and P3 compared to other compounds a significant effect of intramolecular charge transfer. However, the Eg values of P1, P2 and P3 are smaller than that of P6.

The theoretical values of the open circuit voltage Voc of the studied molecules range from 1.499 to 1.804 eV in the case of PC60BM and 0.425 to 0.73 eV in the case of PC71BM, these values are sufficient for a possible efficient electron injection. After the results, we note that all the studied molecules can be used as BHJ because the electron injection process from the excited molecule to the conduction band of PCBM and the subsequent regeneration is possible in an organic sensitized solar cell. It is concluded that We note that the higher power conversion efficiency could be achieved for P2 is 4 and 3% for P3.

The TD-DFT calculations, at least TD-CAM-B3LYP/6-31G(d,p) was used to replicate the optical transitions in order to predict the excited and emission states; the predicted result of the absorption wavelengths for P1, P2, P3, P4, P5, and P6 is 805.02, 794.65, 801.53, 793.82, 790.72 and 727.01 nm respectively.

The decreasing of the band gap of these six materials due to increasing the absorption wavelengths, then the best commands which can be used in photovoltaic cells such as donor of electronic, is one which has the small band gap and large wavelengths, thus all compounds (1–6) are appropriate to do this role.


  1. Sariciftci NS, Heeger AJ, Nalwa HS (1997) Handbook of organic conductive molecules and polymers. Wiley, New York, p 414

    Google Scholar 

  2. Chen HY, Hou J, Zhang S, Liang Y, Yang G, Yang Y, Li G (2009) Polymer solar cells with enhanced open-circuit voltage and efficiency. Nat Photonics 3(11):649–653

    Article  CAS  Google Scholar 

  3. Hoppe H, Sariciftci NS (2006) Morphology of polymer/fullerene bulk heterojunction solar cells. J Mater Chem 16(1):45–61

    Article  CAS  Google Scholar 

  4. Helgesen M, Søndergaard R, Krebs FC (2010) Advanced materials and processes for polymer solar cell devices. J Mater Chem 20(1):36–60

    Article  CAS  Google Scholar 

  5. Park SH, Roy A, Beaupre S, Cho S, Coates N, Moon JS, Heeger AJ (2009) Bulk heterojunction solar cells with internal quantum efficiency approaching 100 and percent. Nat Photonics 3(5):297–302

    Article  CAS  Google Scholar 

  6. Price SC, Stuart AC, Yang L, Zhou H, You W (2011) Fluorine substituted conjugated polymer of medium band gap yields 7 % efficiency in polymer-fullerene solar cells. J Am Chem Soc 133(12):4625–4631

    Article  CAS  Google Scholar 

  7. Zhou H, Yang L, Stuart AC, Price SC, Liu S, You W (2011) Development of fluorinated benzothiadiazole as a structural unit for a polymer solar cell of 7 % efficiency. Angew Chem 123(13):3051–3054

    Article  Google Scholar 

  8. Ma W, Yang C, Gong X, Lee K, Heeger AJ (2005) Thermally stable, efficient polymer solar cells with nanoscale control of the interpenetrating network morphology. Adv Funct Mater 15(10):1617–1622

    Article  CAS  Google Scholar 

  9. Yang C, Lee JK, Heeger AJ, Wudl F (2009) Well-defined donor–acceptor rod–coil diblock copolymers based on P3HT containing C 60: the morphology and role as a surfactant in bulk-heterojunction solar cells. J Mater Chem 19(30):5416–5423

    Article  CAS  Google Scholar 

  10. Lee K, Kim JY, Park SH, Kim SH, Cho S, Heeger AJ (2007) Air-stable polymer electronic devices. Adv Mater 19(18):2445–2449

    Article  CAS  Google Scholar 

  11. Lee JK, Coates NE, Cho S, Cho NS, Moses D, Bazan GC, Heeger AJ (2008) Efficacy of TiOx optical spacer in bulk-heterojunction solar cells processed with 1, 8-octanedithiol. Appl Phys Lett 92(24):3308

    Google Scholar 

  12. Peet J, Kim JY, Coates NE, Ma WL, Moses D, Heeger AJ, Bazan GC (2007) Efficiency enhancement in low-bandgap polymer solar cells by processing with alkane dithiols. Nat Mater 6(7):497–500

    Article  CAS  Google Scholar 

  13. Lee JK, Ma WL, Brabec CJ, Yuen J, Moon JS, Kim JY, Heeger AJ (2008) Processing additives for improved efficiency from bulk heterojunction solar cells. J Am Chem Soc 130(11):3619–3623

    Article  CAS  Google Scholar 

  14. You J, Dou L, Yoshimura K, Kato T, Ohya K, Moriarty T, Yang Y (2013) A polymer tandem solar cell with 10.6 % power conversion efficiency. Nat Commun 4:1446

    Article  Google Scholar 

  15. Chu TY, Lu J, Beaupré S, Zhang Y, Pouliot JR, Wakim S, Tao Y (2011) Bulk heterojunction solar cells using thieno [3,4-c] pyrrole-4,6-dione and dithieno [3, 2-b: 2′, 3′-d] silole copolymer with a power conversion efficiency of 7.3 %. J Am Chem Soc 133(12):4250–4253

    Article  CAS  Google Scholar 

  16. Sharma SS, Sharma GD, Mikroyannidis JA (2011) Improved power conversion efficiency of bulk heterojunction poly(3-hexylthiophene): PCBM photovoltaic devices using small molecule additive. Sol Energy Mater Sol Cells 95(4):1219–1223

    Article  CAS  Google Scholar 

  17. Son HJ, Wang W, Xu T, Liang Y, Wu Y, Li G, Yu L (2011) Synthesis of fluorinated polythienothiophene-co-benzodithiophenes and effect of fluorination on the photovoltaic properties. J Am Chem Soc 133(6):1885–1894

    Article  CAS  Google Scholar 

  18. Amb CM, Chen S, Graham KR, Subbiah J, Small CE, So F, Reynolds JR (2011) Dithienogermole as a fused electron donor in bulk heterojunction solar cells. J Am Chem Soc 133(26):10062–10065

    Article  CAS  Google Scholar 

  19. Small CE, Chen S, Subbiah J, Amb CM, Tsang SW, Lai TH, So F (2012) High-efficiency inverted dithienogermole-thienopyrrolodione-based polymer solar cells. Nat Photonics 6(2):115–120

    Article  CAS  Google Scholar 

  20. Dou L, You J, Yang J, Chen CC, He Y, Murase S, Yang Y (2012) Tandem polymer solar cells featuring a spectrally matched low-bandgap polymer. Nat Photonics 6(3):180–185

    Article  CAS  Google Scholar 

  21. He Z, Zhong C, Su S, Xu M, Wu H, Cao Y (2012) Enhanced power-conversion efficiency in polymer solar cells using an inverted device structure. Nat Photonics 6(9):591–595

    Article  Google Scholar 

  22. Roncali J (2009) Molecular bulk heterojunctions: an emerging approach to organic solar cells. Acc Chem Res 42(11):1719–1730

    Article  CAS  Google Scholar 

  23. Walker B, Kim C, Nguyen TQ (2010) Small molecule solution-processed bulk heterojunction solar cells. Chem Mater 23(3):470–482

    Article  Google Scholar 

  24. Demeter D, Rousseau T, Leriche P, Cauchy T, Po R, Roncali J (2011) Manipulation of the open-circuit voltage of organic solar cells by desymmetrization of the structure of acceptor–donor–acceptor molecules. Adv Funct Mater 21(22):4379–4387

    Article  CAS  Google Scholar 

  25. Sun Y, Welch GC, Leong WL, Takacs CJ, Bazan GC, Heeger AJ (2012) Solution-processed small-molecule solar cells with 6.7 % efficiency. Nat Mater 11:44–48

    Article  CAS  Google Scholar 

  26. Bundgaard E, Krebs FC (2007) Large-area photovoltaics based on low band gap copolymers of thiophene and benzothiadiazole or benzo-bis (thiadiazole). Sol Energy Mater Sol Cells 91(11):1019–1025

    Article  CAS  Google Scholar 

  27. Tian H, Yang X, Cong J, Chen R, Teng C, Liu J, Sun L (2010) Effect of different electron donating groups on the performance of dye-sensitized solar cells. Dyes Pigm 84(1):62–68

    Article  CAS  Google Scholar 

  28. Han H, Liang M, Tang K, Cheng X, Zong X, Sun Z, Xue S (2011) Molecular design of triarylamine dyes incorporating phenylene spacer and the influence of alkoxy substituent on the performance of dye-sensitized solar cells. J Photochem Photobiol A 225(1):8–16

    Article  CAS  Google Scholar 

  29. Kono T, Murakami TN, Nishida JI, Yoshida Y, Hara K, Yamashita Y (2012) Synthesis and photo-electrochemical properties of novel thienopyrazine and quinoxaline derivatives, and their dye-sensitized solar cell performance. Org Electron 13(12):3097–3101

    Article  CAS  Google Scholar 

  30. Campos LM, Tontcheva A, Günes S, Sonmez G, Neugebauer H, Sariciftci NS, Wudl F (2005) Extended photocurrent spectrum of a low band gap polymer in a bulk heterojunction solar cell. Chem Mater 17(16):4031–4033

    Article  CAS  Google Scholar 

  31. Nietfeld JP, Schwiderski RL, Gonnella TP, Rasmussen SC (2011) Structural effects on the electronic properties of extended fused-ring Thieno [3, 4-b] pyrazine analogues. J Org Chem 76(15):6383–6388

    Article  CAS  Google Scholar 

  32. Bourass M et al (2013) DFT theoretical investigations of p-conjugated molecules based on thienopyrazine and different acceptor moieties for organic photovoltaic cells. J Saudi Chem Soc. doi:10.1016/j.jscs.2013.01.003

    Google Scholar 

  33. Bourass M, Fitri A, Benjelloun AT, Mcharfi M, Hamidi M, Serein-Spirau F, Bouachrine M (2013) DFT and TDDFT investigations of new thienopyrazine-based dyes for solar cells: Effects of electron donor groups. Der Pharma Chemica 5(5):144–153

    CAS  Google Scholar 

  34. Becke AD (1993) Density-functional thermochemistry. III. The role of exact exchange. J Chem Phys 98(7):5648–5652

    Article  CAS  Google Scholar 

  35. Magyar RJ, Tretiak S (2007) Dependence of spurious charge-transfer excited states on orbital exchange in TDDFT: large molecules and clusters. J Chem Theory Comput 3(3):976–987

    Article  CAS  Google Scholar 

  36. Yanai T, Tew DP, Handy NC (2004) A new hybrid exchange-correlation functional using the Coulomb-attenuating method (CAM-B3LYP). Chem Phys Lett 393(1):51–57

    Article  CAS  Google Scholar 

  37. Preat J (2010) Photoinduced energy-transfer and electron-transfer processes in dye-sensitized solar cells: TDDFT insights for triphenylamine dyes. J Phys Chem C 114(39):16716–16725

    Article  CAS  Google Scholar 

  38. Camino B, De La Pierre M, Ferrari AM (2013) Photoelectrochemical properties of the CT1 dye: A DFT study. J Mol Struct 1046:116–123

    Article  CAS  Google Scholar 

  39. Irfan A, Jin R, Al-Sehemi AG, Asiri AM (2013) Quantum chemical study of the donor-bridge-acceptor triphenylamine based sensitizers. Spectrochim Acta Part A Mol Biomol Spectrosc 110:60–66

    Article  CAS  Google Scholar 

  40. Jungsuttiwong S, Tarsang R, Sudyoadsuk T, Promarak V, Khongpracha P, Namuangruk S (2013) Theoretical study on novel double donor-based dyes used in high efficient dye-sensitized solar cells: the application of TDDFT study to the electron injection process. Org Electron 14(3):711–722

    Article  CAS  Google Scholar 

  41. Tomasi J, Mennucci B, Cammi R (2005) Quantum mechanical continuum solvation models. Chem Rev 105(8):2999–3094

    Article  CAS  Google Scholar 

  42. Cossi M, Barone V (2001) Time-dependent density functional theory for molecules in liquid solutions. J Chem Phys 115(10):4708–4717

    Article  CAS  Google Scholar 

  43. Adamo C, Barone V (2000) A TDDFT study of the electronic spectrum of s-tetrazine in the gas-phase and in aqueous solution. Chem Phys Lett 330(1):152–160

    Article  CAS  Google Scholar 

  44. Pearson RG (1986) Absolute electronegativity and hardness correlated with molecular orbital theory. Proc Natl Acad Sci 83(22):8440–8441

    Article  CAS  Google Scholar 

  45. Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Montgomery JA Jr, Vreven T, Kudin KN, Burant JC, Millam JM, Iyengar SS, Tomasi J, Barone V, Mennucci B, Cossi M, Scalmani G, Rega N, Petersson GA, Nakatsuji H, Hada M, Ehara M, Toyota K, Ukuda R, Hasegawa J, Ishida M, Nakajima T, Honda Y, Kitao O, Nakai H, Klene M, Li X, Knox JE, Hratchian HP, Cross JB, Adamo C, Jaramillo J, Gomperts R, Stratmann RE, Yazyev O, Austin AJ, Cammi R, Pomelli C, Ochterski JW, Ayala PY, Morokuma K, Voth GA, Salvador P, Dannenberg JJ, Zakrzewski VG, Dapprich S, Daniels AD, Strain MC, Farkas O, Malick DK, Rabuck AD, Raghavachari K, Foresman JB, Ortiz JV, Cui Q, Baboul AG, Clifford S, Cioslowski J, Stefanov BB, Liu G, Liashenko A, Piskorz P, Komaromi I, Martin RL, Fox DJ, Keith T, Al-Laham MA, Peng CY, Anayakkara A, Challacombe M, Gill PMW, Johnson B, Chen W, Wong MW, Gonzalez C, Pople JA (2009) Gaussian 09, Revision A02. Gaussian Inc, Wallingford CT

    Google Scholar 

  46. Shaheen SE, Brabec CJ, Sariciftci NS, Padinger F, Fromherz T, Hummelen JC (2001) 2.5 % efficient organic plastic solar cells. Appl Phys Lett 78(6):841–843

    Article  CAS  Google Scholar 

  47. Wu Z, Fan B, Xue F, Adachi C, Ouyang J (2010) Organic molecules based on dithienyl-2, 1, 3-benzothiadiazole as new donor materials for solution-processed organic photovoltaic cells. Sol Energy Mater Sol Cells 94(12):2230–2237

    Article  CAS  Google Scholar 

  48. Scharber MC, Mühlbacher D, Koppe M, Denk P, Waldauf C, Heeger AJ, Brabec CJ (2006) Design rules for donors in bulk-heterojunction solar cells—towards 10 % energy-conversion efficiency. Adv Mater 18(6):789–794

    Article  CAS  Google Scholar 

  49. Brabec CJ, Cravino A, Meissner D, Sariciftci NS, Fromherz T, Rispens MT, Hummelen JC (2001) Origin of the open circuit voltage of plastic solar cells. Adv Funct Mater 11(5):374–380

    Article  CAS  Google Scholar 

  50. Frohne H, Shaheen SE, Brabec CJ, Müller DC, Sariciftci NS, Meerholz K (2002) Influence of the anodic work function on the performance of organic solar cells. ChemPhysChem 3(9):795–799

    Article  CAS  Google Scholar 

  51. Koster LJA, Mihailetchi VD, Blom PWM (2006) Bimolecular recombination in polymer/fullerene bulk heterojunction solar cells. Appl Phys Lett 88(5):052104

    Article  Google Scholar 

  52. Minnaert B, Burgelman M (2007) Efficiency potential of organic bulk heterojunction solar cells. Prog Photovoltaics Res Appl 15(8):741–748

    Article  CAS  Google Scholar 

  53. May V, Kühn O (2000) Intramolecular Electronic Transitions. Charge and Energy Transfer Dynamics in Molecular Systems, 3rd edn. Wiley, New York, p 255–307

  54. Lukeš V, Aquino A, Lischka H (2005) Theoretical study of vibrational and optical spectra of methylene-bridged oligofluorenes. J Phys Chem A 109(45):10232–10238

    Article  Google Scholar 

Download references

Authors’ contributions

MB, ATB, MB and MM done the quantum calculation, analyzed and interpreted the data of materials, analysis tools or data; wrote the paper. MH, SMB and MB proposed the studied compounds and checked the analyzed and interpreted the data of materials, analysis tools or data. All authors read and approved the final manuscript.


This work was supported by Volubilis Program (No MA/11/248), and the convention CNRST/CNRS (Project chimie 1009).

Competing interests

The authors declare that they have no competing interests.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Mohamed Bourass.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bourass, M., Benjelloun, A.T., Benzakour, M. et al. DFT and TD-DFT calculation of new thienopyrazine-based small molecules for organic solar cells. Chemistry Central Journal 10, 67 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: