Due to increase in the environmental pollution and limited resource of fossil fuel, the development of solar energy harvesting is considered as a promising technology to substitute the fossil fuel, eliminate environmental pollution and sustain development of human societies.
DSSCs have attracted substantial attention, not only because of their environmental friendliness, low manufacturing cost, simple fabrication process and tunable aesthetic color, but also due to their lower dependence on the angle of incident light and higher efficiency under low irradiance compared to traditional Si-based solar cells  - . DSSC is composed of two sandwiched electrodes, working and counter electrode which is made of transparent conducting glass mainly composed of fluorine doped tin oxide (FTO) or indium doped tin oxide. Working electrode is made of thin layer of semiconducting materials such as titanium dioxide (TiO2), zinc oxide (ZnO), tin dioxide (SnO2) coated with dyes on top of FTO  . Counter electrode is made of FTO with thin layer of platinum or carbon . Electrolyte is placed between counter electrode and photoelectrode to mediate electrons . Among of these components of DSSC, the sensitizer (dye) is a crucial element which plays a significant role in light capturing, electron injection and dye regeneration which impact on power conversion efficiency (PCE) as well as the stability of the devices.
However DSSCs are emerging as potential candidates for next generation solar cells owing to low-cost production process  but they face challenges regarding sensitizer such as high price of some dyes, dye degradation, and toxicity of some dyes . Since the invention of DSSCs in the year of 1991 by Brian O’Regan and Michael Grätzel  , ruthenium complexes   , metal-free organic dyes sensitizers      and natural dyes    have been investigated.
Although the Incident Photon to Charge Carrier Efficiency (IPCE) of DSSCs based on ruthenium dyes has reached 12%    under simulated AM 1.5 irradiation (100 mW・cm−2), ruthenium is rare, expensive, and toxic in nature that limits its potentially wide application. Despite the highest efficiency of 13%  has been reported for DSSCs based on free metal dyes (organic dyes), the high price hampers their extensive use. DSSCs based on natural dyes have poor photoelectrochemical properties but the natural dye sensitizers are cheap and environmentally friendly.
To overcome challenges facing DSSCs, theoretical and experimental studies of physical and chemical properties of dyes are required as they contribute to understanding the relationship between structure, properties of dyes and performance of a solar cell and hence facilitate design and synthesis of efficient dyes for DSSCs application   - . In selecting and designing an efficient dye, electron donating part, unit to adjust the absorption spectrum and the electron-acceptor part of the dye  must be taken into design consideration. Mainly the dye to be an efficient sensitizer should be designed with a target of eliminating dye aggregation, minimizing dye recombination and maximizing the driving force for electron injection from the dye to the semiconductor and driving force for dye regeneration, broaden up the absorption spectra which finally improve the overall efficiency of the DSSCs   . Combination of dyes with different absorption properties has been investigated intending to enhance light harvesting efficiency as well as general performance of DSSCs  - .
In this study, the effect of combination of natural dye lawsone and synthetic metal free indoline D131 dye for the purpose of utilization in DSSCs has been investigated. The computations of vibrational and electronic absorption spectra of the indoline D131, lawsone and designed complex molecule as well as mechanisms of electron transitions and energy levels alignment are included in the theoretical part. The UV-Vis spectra of the individual dyes solutions and their mixture are considered in the experimental part.
2. Materials and Methods
2.1. Computational Details
The initial coordinates of indoline and lawsone dyes were taken from ChemSpider database  and NIST  , respectively. The geometrical parameters were optimized using density function theory DFT/B3LYP5 with the different basis sets, 3-21G, 6-31G and 6-31G(d,p). The optimized coordinates were used as inputs for computation of vibrational and electronic spectra.
The complex dye molecule C45H32N2O4 was designed through modeling of esterification reaction between indoline D131 and lawsone dyes. The structure of the designed dye was optimized using DFT/B3LYP5 starting from the lowest STO-3G through 3-21G, 6-31G to highest 6-31G(d,p) basis set. The energy, Gibbs free energy and enthalpy of the complex formation reaction
were computed at different basis sets with the corresponding optimized coordinates. The energy of the reaction was calculated as a difference between the total energies of the products and reactants:
In quantum chemical computations, the total energies of species correspond to temperature 0 K and normal pressure. The enthalpy of the reaction was found using the ∆rE and zero point vibration energy (ZPVE) Δrε
where h is the Plank’s constant, c is the speed of light in free space, and are the sums of the vibrational frequencies of the products and reactants, respectively. The Gibbs free energy of the reaction was obtained by the following equations:
where and are the enthalpy and entropy of the reaction at temperature T. The thermodynamic functions, entropies and enthalpy increments , of the participants of the reaction have been calculated with OpenThermo software .
The vertical excitation and electronic absorption spectra of the lawsone dye, indoline D131 dye and the designed complex molecules both in vacuum and solvent (chloroform) were computed at the TD-DFT B3LYP5/6-31G(d,p) level of theory. The solvent effects were accounted by employing the polarized continuum model (PCM). All quantum chemical computations were performed using the Firefly QC package  which is partially based on the GAMESS (US) source code . The geometrical structures, vibrational and electronic spectra were visualized and analyzed using the Chemcraft  and MacMolPlt software .
2.2. Experimental Details
The indoline D131 dye with a purity of 97% and lawsone dye with purity 97% were bought in powder form from Sigma-Aldrich Chemical Company and were used without further purification for solution preparation at room temperature. A single solvent both for D131 and lawsone was needed. Chloroform was found to be appropriate solvent   ; it was bought from Sigma-Aldrich Chemical supplier.
The solutions of 1.15 × 10−4 M for lawsone and 2.95 × 10−5 M for indoline D131 were prepared and used for the UV-Vis spectra measurements. The mass 8 mg of lawsone dissolved in 40 ml of chloroform to have 1.15 × 10−3 M and then 1 ml of the solution was added to 9 ml of chloroform resulting in 1.15 × 10−4 M. For the case of indoline D131 dye, 6 mg of the sample was dissolved in 40 ml of chloroform resulting in 2.95 × 10−4 M, then 2 ml of this solution was added to 18 ml of chloroform to get 2.95 × 10−5 M indoline D131 solution. The solution of the dyes mixture was prepared by combination of the indoline D131 and lawsone dye solutions in the ratio of 1:1. The UV-Vis spectra of the samples were recorded in the region 190 - 820 nm on a single beam 2800 UV-Vis spectrophotometer (Hitachi U2000) with quartz cuvettes of 10 mm path length.
3. Results and Discussion
3.1. Molecular Structure of Individual Dyes and Complex Molecule
The equilibrium geometrical structures of lawsone, indoline D131 and the designed complex molecules computed at B3LYP5/6-31G(d,p) are given in Figure 1.
Figure 1. The optimized geometrical structure of molecules: (a) lawsone, (b) indoline D131, (c) complex. The visualization of molecular structures was performed by using the ChemCraft software .
The lawsone molecule (2-hydroxy-1, 4-naphthoquinone molecule) contains two six-membered rings, two ketone groups and one hydroxyl group (Figure 1(a)). Indoline dye D131 is composed of 4-(2, 2 diphenylethenyl) phenyl) substituent acting as a donor and π conjugated to the indoline ring, cynoacrylic acid group as an acceptor  (Figure 1(b)).
In the complex molecule, the coupling of the molecules occurs via a new chemical bond C1-O2-C37 (Figure 1(c)). The formation of this bond is accompanied by detachment of the H1-atom from cynoacrylic acid group of indoline D131 dye and hydroxyl O1-H6 from lawsone dye to release water molecule. Schematic representation of H2O molecule detachment is shown below.
The selected geometrical parameters of the individual dyes are compared with those in the complex molecule in Table 1 and Table 2; the differences between respective parameters are shown in the far right columns. The bond C1-O2-C37
Table 1. Comparison of selected geometrical parameters of lawsone and complex molecule computed at B3LYP5/6-31G(d,p).
Table 2. Comparison of selected geometrical parameters of indoline D131 dye and complex molecule computed at B3LYP5/6-31G(d,p).
in the complex is specified with parameters R(C1-O2) = 1.397 Å, R(O2-C37) = 1.370 Å and ÐC1-O2-C37 = 117.0˚; which are comparable with those of the ester C-O-C linkage, 1.43 Å   and 119˚ . The most noticeable change is observed in the vicinity of the new chemical bond: +0.026 Å (elongation) for C37-O2 bond in the lawsone moiety and −0.022 Å (shortening) for C1=O1 in the indoline part. The biggest change is seen in the bond angles, +3.6˚ (O2-C37-C36) and −4.7˚ (O2-C37-C38) in the lawsone part, and in the dihedral angles, 3.9˚ (C17-C16-C6-C9) and +4.2˚ (C10-C11-C1-O1) in the indoline part. The difference in other respective parameters is not significant; it does not exceed 0.009 Å for bond lengths and 1.5˚ for bond angles. The parameters of the far parts of the complex molecule remain almost the same as in the individual molecules.
The thermodynamic characteristics of the complex formation reaction (1): energy ∆rE, zero point vibration energy ∆rε, enthalpies ∆rH˚(0) and ∆rH˚(298), Gibbs free energy ∆rG˚(298) and entropy ∆rS˚(298) computed at different basis sets are given in Table 3.
As is seen, the extension of the basis set from 3-21G to 6-31G results in the decrease of the value of ∆rE by 15.3 kJ・mol−1 but the polarization functions addition does not influence much on this value. The values of enthalpy and Gibbs free energy follow the trend of the energy ∆rE as the ZPVE values are not much sensitive to the basis set change.
The values of enthalpies ∆rH˚ are positive that implies the reaction (1) is endothermic. The entropy ∆rS˚(298) is negative hence depicting an increase of the order of the system. The Gibbs free energy ∆rG˚(298) being positive shows non-spontaneity of the complex formation through direct combination of the dye molecules.
3.2. Vibrational Spectra of the Species
The IR spectra of the individual dyes and complex molecule computed at B3LYP5/6-31G(d,p) are shown in Figures 2(a)-(c). The absence of imaginary frequencies proves that the geometrical structures of the species correspond to minima on the potential energy surfaces.
Table 3. The thermodynamic characteristics of the complex formation reaction (1): energy ∆rE, zero point vibration energy ∆rε, enthalpies ∆rH˚(0) and ∆rH˚(298), Gibbs free energy ∆rG˚(298) and entropy ∆rS˚(298).
Figure 2. The IR spectra of molecules: (a) lawsone, (b) indoline D131, (c) complex.
The most intensive modes in the spectrum of lawsone are assigned to vibrations 1680 cm−1 (C=C stretching), 1280 cm−1 (C-C and C-O stretching), 1740 and 1770 cm−1 (C=O stretching), 1180 cm−1 (C-C-H and C-O-H bending), 990 cm−1 (C-C-C bending), 520 cm−1 (C-O-H out of plane bend), and 3800 cm−1 (O-H stretching); and those in indoline are 1530 - 1550 cm−1 (C-C and C-N stretching in cycles), 1390 cm−1 (C-C and C-O stretching near carboxylic group), 1810 cm−1 (C=O stretching), 1190 cm−1 (C-C stretching and C-C-C bending near carboxylic group), 1080 cm−1 (C-C-O-H bending), 3770 cm−1 (O-H stretching). All these frequencies are close to typical vibrational wavenumbers of the respective functional groups.
The IR spectrum of the complex molecule apparently contains the vibrational modes of the lawsone part, e.g., 1763 cm−1 (C=O stretching), 1696 cm−1 (C=C stretching), and indoline part, 1820 cm−1 (C=O stretching) and 1530 - 1550 cm−1 (C-C and C-N stretching), but these modes slightly differ from those of the individual dyes.
The uniqueness of the spectrum in comparison to individual dyes is the appearance of very intensive bands which relate to the new formed C37-O2-C1 bond, these modes at 1203 and 1231 cm−1 are assigned to the C-O stretching motion; in addition the peak at 1023 cm−1 is observed which corresponds to stretching of C-C and C-O bonds in vicinity of the new bond. Evidently the complex molecule IR spectrum does not have O-H vibration mode as the hydroxyl group is removed through the esterification reaction.
3.3. Electronic Spectra of the Species
The electronic absorption spectra for lawsone, indoline D131 and the designed complex are computed both for vacuum and chloroform solution at the TD-DFT B3LYP5/6-31G(d,p); the PCM is used to account the solvent effect.
The vertical excitation energies, wavelengths, oscillator strengths and electronic configurations of the transitions are listed in Table 4. The singlet?singlet excitations with nonzero oscillator strengths (f > 0.001) are taken into account. For the lawsone molecule the most probable transitions in vacuum are to the fifth and eighth excited states and attributed to H-4 → L and H-1 → L + 1, while in chloroform solution to the fifth and seventh excited states assigned to H-3 → L and H → L + 1 transitions. For the indoline both in vacuum and solution the highest oscillator strength is observed for the electron transition H → L, to the first excited state. For the complex molecule the most probable transition occurs to the third excited state and assigned to H → L + 1 transition, both for vacuum and solution. Worth to note that in spectra of lawsone the band of the maximum absorption wavelength λmax is in the ultraviolet region while for the D131 and designed complex molecule the values of λmax correspond to visible region. Moreover for the complex, the additional peak is appeared at 655 (vacuum) and 814 nm (chloroform) but of low probability. From vacuum to chloroform solution, the red shift of λmax is observed for all three species that is mostly common
Table 4. The excitation energies Eex, absorption wavelengths λ, oscillator strengths f and electronic transition configurations of the individual dyes and designed complex molecule C45H32N2O4 computed at the TD-DFT B3LYP5/6-31G(d,p).
effect mentioned elsewhere e.g.   .
The theoretical electronic absorption spectra are compared to experimental UV-Vis spectra of dyes in chloroform solution in Figures 3-5.
The theoretical UV-Vis spectra of lawsone (Figure 3(a) and Figure 3(b)) depict a red shift of maximum wavelength from 279 (vacuum) to 283 nm (chloroform). The experimental UV-Vis spectrum of lawsone in chloroform solution shows the highest absorbance band in UV region with absorption peak at 292 nm and lower absorbance band at 338 nm. Thus both theoretical and experimental results are in accordance and indicate maximal peaks in UV region with slight differences in the peaks position.
In the theoretical UV-Vis spectrum of the indoline D131 dye, the peak of maximum wavelength and highest intensity is observed at 447 (vacuum) and 476 nm (chloroform) hence the red shift is seen from vacuum to solution. The experimental spectrum of the indoline solution (Figure 4(c)) has broad bands with maxima located at 294, 334, 354 and 450, 488 nm. Thus the experimental results are comparable with the theoretical and also show similarity with data reported in  .
In the theoretical UV-Vis spectrum of the designed complex molecule (Figure 5(a) and Figure 5(b)), the highest peak is at 450 (vacuum) and 488 nm (chloroform), thus the red shift is also seen. In the experimental spectrum of the mixture (Figure 5(c)), the absorption bands lie in the UV-Vis region with maximum absorption peak at 500 nm. Two peaks in visible region (450 and 500 nm) apparently correspond to those in spectrum of the individual indoline D131 dye (450 and 488 nm) while other peaks correspond to both individual dyes. Generally despite the slight dissimilarities, both the peak intensities and wavelengths in the spectrum of the mixture refer to those of the individual dyes.
Figure 3. UV-Vis spectrum of lawsone: (a) computed for vacuum; (b) computed for chloroform solution; (c) experimental for chloroform solution.
Figure 4. UV-Vis spectrum of indoline D131: (a) computed for vacuum; (b) computed for chloroform solution; (c) experimental for chloroform solution.
Figure 5. UV-Vis spectrum: (a) of complex molecule computed for vacuum; (b) of complex molecule computed for chloroform solution; (c) measured experimentally for the mixture of lawsone and indoline in chloroform solution.
3.4. Analysis of Frontier Molecular Orbitals and Energy Levels Alignment
Analysis of molecular orbitals can provide useful information on electronic structure of molecules and mechanisms of electron transitions and chemical reactions    . However for understanding the mechanism of electron transitions the probabilities of the transitions must be taken into account, therefore the analysis of the whole electronic spectrum is required.
The frontier and adjacent molecular orbitals which mostly participate in the transitions are selected according to the data listed in Table 4 and shown in Figures 6-8 as computed for chloroform solution.
For the lawsone molecule, the shape of orbitals depicts that H-3 and HOMO are delocalized over both rings, benzenoid and quinoid, whereas LUMO and L + 1 are centered mainly on the quinoid ring. This implies that π → π* transitions involve electron localization to the quinoid ring of lawsone.
In the D131, the shape of the HOMO indicates that the electron density is delocalized on the extended donor unit, while the LUMO is localized on the acceptor groups (cyanoacrylic acid). The HOMO → LUMO transition contributes mostly to the first excited state transition (Table 4). For the H-1 MO, electron density is delocalized over the entire molecule except the cyclopentane group.
For the designed complex molecule, the HOMO has the electron density only on the indoline part, mainly on the benzene rings, indole unit and cynoacrylic group adjacent to the linkage with lawsone. In the LUMO, the electron density is located mostly on the lawsone part with a small contribution from the linkage. In the L+1 MO, the electrons are concentrated on the indole unit and cynoacrylic part, while in both end parts and cyclopentane no any electron density is observed. The shape of the H-1 MO resembles that of HOMO since electrons are located in the indoline moiety and absent on the lawsone part.
In the most probabale transitions, H → L + 1 and H-1 → L + 1, the excitons formed are coupled via the new C-O-C bond between the indoline and lawsone moieties. Thus the excitons’ size in the complex is enlarged compared to individual molecules. Therefore we suggest that recombination processes might be
Figure 6. Frontier and adjacent MOs of lawsone molecule: (a) HOMO and LUMO; (b) H-3 MO; (c) L + 1 MO.
Figure 7. Frontier and adjacent MOs of indoline D131 molecule: (a) HOMO and LUMO; (b) H-3 MO.
hindered in the complex that in turn favours sensitizing efficiency of the dye.
Proper energy level alignment between semiconductor, sensitizer and electrolyte is among of the basic requirements for DSSCs  . The simplest approach is to sketch the frontier MOs energy levels together with conduction and valence bands of a semiconductor and redox potential of electrolyte. Nevertheless more strict approach is to implement the excitation energies obtained by TD-DFT and calculate the excited state oxidation potentials (ESOPs)   . The values of ESOP were found as sum of the ground state occupied MOs energies and respective excitation energies, e.g. ESOP = ε(HOMO) + Eex. The
Figure 8. Frontier and adjacent MOs of complex molecule: (a) HOMO and LUMO; (b) L + 1 MO; (c) H-1 MO.
calculated energies of frontier and adjacent MOs, excitation energies, energy gaps between relevant MOs and ESOPs are shown in Table 5. As is seen the energy gaps for each dye are systematically bigger compared to respective excitation energies; the maximum difference between Eg and Eex is observed for lawsone dye, 0.60 eV and 0.74 eV for vacuum and chloroform solution, respectively. Thus the ESOPs are lower compared to respective unoccupied orbitals.
The energy levels of the relevant occupied MOs and ESOPs of the molecules under the study are shown in Figure 9; the results are compared with the valence and conduction bands of semiconductors TiO2    , ZnO   and redox potential level of the electrolyte . According to photosensitizer requirements, HOMO and LUMO levels of a dye must match conduction band (CB) of semiconductor and redox potential of electrolyte; the HOMO level of the
Table 5. The energies of molecular orbitals ε(MO), excitation energies Eex, energy gaps Eg, and excited state oxidation potentials ESOPs; all values in eV.
Figure 9. Energy level diagram of the relevant occupied MOs, transition energies and ESOPs of lawsone, indoline D131 and designed complex molecule computed for chloroform solution. The conduction and valence bands of TiO2    and ZnO   semiconductors and redox potential of the electrolyte  are given for comparison.
sensitizer has to be lower than the redox potential for efficient dye regeneration, while the LUMO must be slightly higher than the CB edge so that the electron injection to semiconductor is energetically favorable .
In our case instead of LUMO we consider the ESOP energy levels (Figure 9). The HOMO and H-1 levels of lawsone are far below the redox potential of the electrolyte, whereas the ESOPs are above the CB edge of TiO2 but only ESOP = −1.6 eV is above the CB of ZnO. The HOMO and H-1 levels of the indoline and designed molecule are below the electrolyte redox potential, while the ESOPs are above the CB edges of both semiconductors. Therefore regarding the proper energy levels alignment, the designed complex, lawsone and indoline satisfy the requirement to be used as a photosensitizer for TiO2 or ZnO semiconductors. The designed complex is predicted to be better sensitizer as it possesses smaller excitation energies (3.20, 1.52 and 2.54 eV) compared to the indoline (3.28 and 2.60 eV) and lawsone (4.38, 3.27 and 5.39 eV). Furthermore, the transition from HOMO to the first excited state favours sensitization of TiO2 as the respective ESOP = −3.63 eV is close to the CB edge (−4.0 eV). The transitions from H-1 and from HOMO to the second excited state are of high probability to be involved in sensitizing the ZnO semiconductor, as the two corresponding ESOP levels (−2.69 and −2.61 eV) are slightly above the CB of ZnO (−2.96 eV).
The complex molecule was theoretically designed via combination of natural lawsone dye and synthetic metal free indoline D131 dye. The equilibrium geometry structures, electronic and IR spectra of the individual dyes and complex molecule were computed using DFT and TD-DFT methods; also the UV-Vis spectra of the individual dyes’ solutions and their mixture were measured and analyzed. In the complex molecule, the formation of the link C-O-C between the lawsone and indoline moieties is accompanied with slight change of the geometrical parameters of the moieties, but clearly exhibited in the IR spectrum of the complex as the most intensive peaks originate from the new formed bond. Based on thermodynamic approach, the complex formation reaction is found to be endothermic and non-spontaneous. Compared to the individual dyes, in the theoretical optical spectrum of the designed molecule the highest peaks are red shifted. As for the experimental UV-Vis spectrum of the mixture, both the peak intensities and wavelengths, in general refer to those of the individual dyes. The electron density distribution over frontier and adjacent MOs and most probable electronic transitions have been scrutinized. The electron excitations in the complex bring to lager size of the excitons that might hinder electron-hole recombination. Also the complex possesses better energy level alignment with the conduction band of the TiO2 and ZnO semiconductors and redox level of the electrolyte. Therefore we suggest the designed complex is a promising sensitizer for efficient DSSCs with better PCE compared to that of the individual dyes.
 Yuan, H., Wang, W., Xu, D., Xu, Q., Xie, J., Chen, X., Zhang, T., Xiong, C., He, Y., Zhang, Y., Liu, Y. and Shen, H. (2018) Outdoor Testing and Ageing of Dye-Sensitized Solar Cells for Building Integrated Photovoltaics. Solar Energy, 165, 233-239.
 Abdalla, J., Huang, Y.-W., Yu, Q.-J., Wang, J.-Z., Wang, J.-N., Yu, C.-L., Gao, S.-Y., Jiao, S.-J., Wang, D.-B. and Alarabi, A. (2017) TiCl4 Surface-Treated SnO2 Photoanodes for Self-Powered UV Photodetectors and Dye-Sensitized Solar Cells. Materials Technology, 32, 443-450.
 Arbab, A.A., Peerzada, M.H., Sahito, I.A. and Jeong, S.H. (2017) A Complete Carbon Counter Electrode for high Performance Quasi Solid State Dye Sensitized Solar Cell. Journal of Power Sources, 343, 412-423.
 Nishio, Y., Yamaguchi, T., Yamguchi, T., Nishio, K. and Hayase, S. (2016) Transparent Conductive Oxide-Less Dye-Sensitized Solar Cells (TCO-Less DSSC) with Titanium Nitride Compact Layer on Back Contact Ti Metal Mesh. Journal of Applied Electrochemistry, 46, 551-557.
 Shalini, S., Balasundaraprabhu, R., Kumar, T.S., Prabavathy, N., Senthilarasu, S. and Prasanna, S. (2016) Status and Outlook of Sensitizers/Dyes Used in Dye Sensitized Solar Cells (DSSC): A Review. International Journal of Energy Research, 40, 1303-1320.
 Lee, C.W., Lu, H.P., Lan, C.M., Huang, Y.L., Liang, Y.R., Yen, W.N., Liu, Y.C., Lin, Y.S., Diau, E.W.G. and Yeh, C.Y. (2009) Novel Zinc Porphyrin Sensitizers for Dye-Sensitized Solar Cells: Synthesis and Spectral, Electrochemical, and Photovoltaic Properties. Chemistry—A European Journal, 15, 1403-1412.
 Hara, K., Sato, T., Katoh, R., Furube, A., Ohga, Y., Shinpo, A., Suga, S., Sayama, K., Sugihara, H. and Arakawa, H. (2003) Molecular Design of Coumarin Dyes for Efficient Dye-Sensitized Solar Cells. The Journal of Physical Chemistry B, 107, 597-606.
 Santhanamoorthi, N., Lo, C.-M. and Jiang, J.-C. (2013) Molecular Design of Porphyrins for Dye-Sensitized Solar Cells: A DFT/TDDFT Study. The Journal of Physical Chemistry Letters, 4, 524-530.
 Wu, Y., Marszalek, M., Zakeeruddin, S.M., Zhang, Q., Tian, H., Gr?tzel, M. and Zhu, W. (2012) High-Conversion-Efficiency Organic Dye-Sensitized Solar Cells: Molecular Engineering on D-A-π-A Featured Organic Indoline Dyes. Energy & Environmental Science, 5, 8261-8272.
 Wongcharee, K., Meeyoo, V. and Chavadej, S. (2007) Dye-Sensitized Solar Cell Using Natural Dyes Extracted from Rosella and Blue Pea Flowers. Solar Energy Materials and Solar Cells, 91, 566-571.
 Chang, H., Wu, H., Chen, T., Huang, K., Jwo, C. and Lo, Y. (2010) Dye-Sensitized Solar Cell Using Natural Dyes Extracted from Spinach and Ipomoea. Journal of Alloys and Compounds, 495, 606-610.
 Chiba, Y., Islam, A., Watanabe, Y., Komiya, R., Koide, N. and Han, L. (2006) Dye-Sensitized Solar Cells with Conversion Efficiency of 11.1 %. Japanese Journal of Applied Physics, 45, L638-L640.
 Mathew, S., Yella, A., Gao, P., Humphry-Baker, R., Curchod, B.F., Ashari-Astani, N., Tavernelli, I., Rothlisberger, U., Nazeeruddin, M.K. and Gratzel, M. (2014) Dye-Sensitized Solar Cells with 13% Efficiency Achieved through the Molecular Engineering of Porphyrin Sensitizers. Nature Chemistry, 6, 242-247.
 Xu, J., Zhang, H., Liang, G., Wang, L., Weilin, X., Cui, W. and Zengchang, L. (2010) DFT Studies on the Electronic Structures of Indoline Dyes for Dye-Sensitized Solar Cells. Journal of the Serbian Chemical Society, 75, 259-269.
 Han, M.W., Ekanayake, P., Ming, L.C. and Yoong, V.N. (2015) DFT/TD-DFT Studies on the Lawsone (Henna) as a Photosensitizer for Dye-Sensitized Solar Cells. Applied Mechanics and Materials, 789-790, 56-60.
 Sun, C., Li, Y., Song, P. and Ma, F. (2016) An Experimental and Theoretical Investigation of the Electronic Structures and Photoelectrical Properties of Ethyl Red and Carminic Acid for DSSC Application. Materials, 9, 1-22.
 Le Bahers, T., Pauporté, T., Scalmani, G., Adamo, C. and Ciofini, I. (2009) A TD-DFT Investigation of Ground and Excited State Properties in Indoline Dyes Used for Dye-Sensitized Solar Cells. Physical Chemistry Chemical Physics, 11, 11276-11284.
 Keerthi, A., Chua, M.H., Timothy Chan, T.Y., Liu, Y., Wang, Q. and Valiyaveettil, S. (2016) Synthesis of Multi-Donor Dyes and Influence of Molecular Design on Dye-Sensitized Solar Cells. The Royal Society of Chemistry, 6, 51807-51815.
 Ren, P., Zhang, Y., Luo, Z., Song, P. and Li, Y. (2017) Theoretical and Experimental Study on Spectra, Electronic Structure and Photoelectric Properties of Three Nature Dyes Used for Solar Cells. Journal of Molecular Liquids, 247, 193-206.
 Morita, S., Ikegami, M., Wei, T.C. and Miyasaka, T. (2014) Quantum Conversion Enhancement with TiOx Compact Layers for ITO-Plastic-Film-Based Low-Tempe-rature-Processed Dye-Sensitized Photoelectrodes. European Journal of Chemical Physics and Physical Chemistry, 15, 1190-1193.
 Shibayama, N., Inoue, Y., Abe, M., Kajiyama, S., Ozawa, H., Miura, H. and Arakawa, H. (2015) Novel Near-Infrared Carboxylated 1, 3-Indandione Sensitizers for Highly Efficient Flexible Dye-Sensitized Solar Cells. Chemical Communications, 51, 12795-12798.
 Sun, X., Wang, Y., Li, X., Agren, H., Zhu, W., Tian, H. and Xie, Y. (2014) Cosensitizers for Simultaneous Filling up of both Absorption Valleys of Porphyrins: A Novel Approach for Developing Efficient Panchromatic Dye-Sensitized Solar Cells. The Royal Society of Chemistry, 50, 1-4.
 Msangi, A., Pogrebnoi, A. and Pogrebnaya, T. (2018) Combination of Natural Dye (Crocetin) and Synthetic Dye (Indoline D205) for DSSCs Application. International Journal of Computational and Theoretical Chemistry, 6, 1-13.
 Sharma, G.D., Angaridis, P.A., Pipou, S., Zervaki, G.E., Nikolaou, V., Misra, R. and Coutsolelos, A.G. (2015) Efficient Co-Sensitization of Dye-Sensitized Solar Cells by Novel Porphyrin/Triazine Dye and Tertiary Aryl-Amine Organic Dye. Elsevier, 25, 295-307.
 Nguyen, L.H., Mulmudi, H.K., Sabba, D., Kulkarni, S.A., Batabyal, S.K., Nonomura, K., Gratzel, M. and Mhaisalkar, S.G. (2012) A Selective Co-Sensitization Approach to Increase Photon Conversion Efficiency and Electron Lifetime in Dye-Sensitized Solar Cells. Physical Chemistry Chemical Physics, 14, 16182-16186.
 Xie, Y., Wu, L., Han, L. and Gao, J. (2017) Co-Sensitization by Triarylamine Dyes for Improved Dye-Sensitized Solar Cells. Physical Status Solidi Applications and Materials Science, 214, 1-6.
 Pamain, A., Pogrebnaya, T. and King'ondu, C.K. (2014) Natural Dyes for Solar Cell Application: UV-Visible Spectra and Outdoor Photovoltaic Performance. Research Journal in Engineering and Applied Sciences, 3, 332-336.
 Lemmon, E., McLinden, M., Friend, D., Linstrom, P. and Mallard, W. (2011) NIST Chemistry WebBook, Nist Standard Reference Database Number 69. National Institute of Standards and Technology, Gaithersburg.
 Schmidt, M.W., Baldridge, K.K., Boatz, J.A., Elbert, S.T., Gordon, M.S., Jensen, J.H., Koseki, S., Matsunaga, N., Nguyen, K.A. and Su, S. (1993) General Atomic and Molecular Electronic Structure System. Journal of Computational Chemistry, 14, 1347-1363.
 Mahkam, M., Kafshboran, H.R. and Nabati, M. (2014) Synthesis and Characterization of Novel Colored Polymers Based on Lawsone Natural Compound. Designed Monomers and Polymers, 17, 784-794.
 Dentani, T., Kubota, Y., Funabiki, K., Jin, J., Yoshida, T., Minoura, H., Miura, H. and Matsui, M. (2009) Novel Thiophene-Conjugated Indoline Dyes for Zinc Oxide Solar Cells. New Journal of Chemistry, 33, 93-101.
 Allen, F.H. and Kirby, A.J. (1984) Bond Length and Reactivity. Variable Length of the Carbon-Oxygen Single Bond. Journal of the American Chemical Society, 106, 6197-6200.
 Hrmova, M., Varghese, J.N., De Gori, R., Smith, B.J., Driguez, H. and Fincher, G.B. (2001) Catalytic Mechanisms and Reaction Intermediates along the Hydrolytic Pathway of a Plant β-D-Glucan Glucohydrolase. Structure, 9, 1005-1016.
 Makuraza, J., Pogrebnaya, T. and Pogrebnoi, A. (2015) Vibrational and Electronic Spectra of Natural Dyes Constituents for Solar Cell Application: DFT and TDDFT Study. International Journal of Materials Science and Applications, 4, 314-324.
 Liu, Z., Wang, X.-F., Wang, Z., Ojima, H., Hong, Z., Tian, W. and Kido, J. (2013) Indoline-Based Donor Molecule for Efficient Co-Evaporated Organic Photovoltaics. Organic Electronics, 14, 2210-2215.
 Lin, H., Zhu, S.-G., Chen, P.-Y., Li, K., Li, H.-Z. and Peng, X.-H. (2013) DFT Investigation of a High Energy Density Polynitro Compound, 2,2’-Bis(trinitromethyl)-5,5’-azo-1,2,3,4-tetrazole. Central European Journal of Energetic Materials, 10, 325-338.
 Wu, C.J., Yang, L.H., Fried, L.E., Quenneville, J. and Martinez, T.J. (2003) Electronic Structure of Solid 1,3,5-triamino-2,4,6-trinitrobenzene under Uniaxial Compression: Possible Role of Pressure-Induced Metallization in Energetic Materials. The American Physical Society, 67, 1-7.
 Wei, T., Zhu, W., Zhang, X., Li, Y.-F. and Xiao, H. (2009) Molecular Design of 1, 2, 4, 5-Tetrazine-Based High-Energy Density Materials. The Journal of Physical Chemistry A, 113, 9404-9412.
 Badders, N., Wei, C., Aldeeb, A., Rogers, W. and Mannan, M. (2006) Predicting the Impact Sensitivities of Polynitro Compounds Using Quantum Chemical Descriptors. Journal of Energetic Materials, 24, 17-33.
 Lungu, J., Oprea, C., Dumbrava, A., Enache, I., Georgescu, A., Radulescu, C., Ionita, I., Cimpoca, G. and Girtu, M. (2010) Heterocyclic Azodyes as Pigments for Dye Sensitized Solar Cells: A Combined Experimental and Theoretical Study. Journal of Optoelectronics and Advanced Materials, 12, 1969-1975.
 Oprea, C.I., Dumbrav?, A., Enache, I., Lungu, J., Georgescu, A., Moscalu, F., Oprea, C. and G?r?u, M.A. (2011) Role of Energy Level Alignment in Solar Cells Sensitized with a Metal-Free Organic Dye: A Combined Experimental and Theoretical Approach. Physical Status Solidi Applications and Materials Science, 208, 2467-2477.
 Oprea, C.I., Frecu?, B., Minaev, B.F. and G?r?u, M.A. (2011) DFT Study of Electronic Structure and Optical Properties of Some Ru-and Rh-Based Complexes for Dye-Sensitized Solar Cells. Molecular Physics, 109, 2511-2523.
 De Angelis, F., Fantacci, S. and Selloni, A. (2008) Alignment of the Dye’s Molecular Levels with the TiO2 Band Edges in Dye-Sensitized Solar Cells: A DFT-TDDFT Study. Nanotechnology, 19, 1-7.
 Xu, Y. and Schoonen, M.A.A. (2000) The Absolute Energy Positions of Conduction and Valence Bands of Selected Semiconducting Minerals. American Mineralogist, 85, 543-556.