Abstract

The study concerns molecular design and theoretical study of novel prospective organic sensitizers based on D–A–π–A concept. A modification of the concept (DA–π–A) and new ways for design of organic sensitizers are proposed. The study shows that π-block plays a key role for the wavelength shift in UV-Vis spectra. An existence of element of higher periods in the π-block provides a noticeable bathochromic shift. The thienyl bridge in A-block provides a red shift with respect to phenyl bridge for 1,2,5oxadiazolo3,4-cpyridines while the phenyl bridge provides a red shift with respect to thienyl bridge for 1,2,5selenadiazolo3,4-cpyridines. The maximal wavelengths are observed for DA-blocks with aliphatic/alicyclic fragments we proposed. An influence of enthalpy of formation and polarity of DA-blocks on bathochromic shift is shown. The quantitative dependencies of wavelengths on the enthalpies and dipole moments are determined. The maximal oscillator strengths are observed for the planar DA-block. Therefore, the conditions for design of compounds possessing the maximal wavelengths and oscillator strengths are different: the latter presupposes the planar structure of DA-block while the first presupposes an existence of nonplanar aliphatic/alicyclic fragments. Thus, the further design of organic sensitizers should include some combination of these concepts.

1. Introduction

Chemistry of chalcogen-containing compounds is currently one of the most rapidly developing areas of heterocyclic chemistry. It has been established that a number of structures at the boundary of organic and inorganic chemistry, characterized by a high ratio of heteroatoms (sulfur, selenium, and nitrogen) relatively to carbon, have useful properties as materials with interesting optical and electrically conductive properties such as high biological activity [14]. For example, some of them containing 1,2,5-chalcogenadiazole ring are prospective high performance and low cost components of dye-sensitized solar cells [57], anticancer and anti-HIV1 agents [811], which have drawn great interest of academic and industrial specialists. Therefore, a design of new 1,2,5-chalcogenadiazole containing compounds that determine properties useful for engineering and medicine is of a great importance. Thus, the objective of the paper is a theoretical study of the previously unknown heterocyclic systems containing 1,2,5-chalcogenadiazole ring which may be of interest for specific materials with useful properties for components of small-molecule organic solar cells (SMOSCs). The mainstream for the design of chalcogen-containing heterocyclic systems SMOSCs for these aims is a construction of “donor–p bridge–acceptor” (D–p–A) configuration due to its convenient modulation of the intramolecular charge-transfer nature. Recently, it has been shown [7] that incorporation of additional electron-acceptors (such as benzothiadiazole, benzotriazole, quinoxaline, phthalimide, diketopyrrolopyrrole, thienopyrazine, thiazole, triazine, cyanovinyl, and cyano- and fluorosubstituted phenyl) into the p bridge (termed the D–A–p–A) configuration displays several advantages such as regulation of the molecular energy levels, bathochromic shift in UV-Vis spectrum, and distinct improvement of photovoltaic performance and stability. A group from East China University of Science and Technology has proposed a novel “D–A–π–A” [7] concept for designing novel organic sensitizers, in which several kinds of electron-withdrawing units are incorporated into the π bridge to tailor molecular structures and optimize energy levels. They have systematically demonstrated that the incorporated electron-withdrawing additional acceptor can be treated as an “electron trap,” showing several distinguished merits such as (i) essentially facilitating the electron transfer from the donor to the acceptor/anchor; (ii) conveniently tailoring the solar cell performance with a facile structural modification on the additional acceptor; (iii) improving circuit photovoltage with the nitrogen-containing heterocyclic group; (iv) conveniently tuning the molecular energy gap and modulating the response of the light-harvesting range with the new resulting absorption band; and (v) most importantly being capable of greatly improving the sensitizer photostability. Organic sensitizers containing additional electron-withdrawing units (or D–A–π–A dyes) are reviewed with specific concern on the relationship between molecular structures and absorption and energy levels as well as photovoltaic performances. From the quantum viewpoint the spectral and electronic properties are dependent on the energy and distribution of the highest occupied (HOMO) and the lowest unoccupied (LUMO) molecular orbital such as being dependent on the threshold between their energies. Therefore, more specifically, our study is focused on fused [1,2,5]thiadiazolo[3,4-c]pyridine and [1,2,5]selenadiazolo[3,4-c]pyridine derivatives incorporated into D–A–π–A configuration and their comparison with the similar [1,2,5]oxadiazolo[3,4-c]pyridine derivatives using DFT B3LYP calculations at the 6/311G(d,p) level of theory for the further design of new compounds with perspective properties. The specific interest to [1,2,5]thiadiazolo[3,4-c]pyridine and [1,2,5]selenadiazolo[3,4-c]pyridine derivatives is concerned with the fact that these derivatives condensed with strong electronegative rings such as pyridine or pyridazine may be used in organic solar cells due to their ability to convert light into cheap electricity.

Taking into account D–A–π–A concept it is necessary to select the building blocks, that is, suitable donor groups (D), acceptor groups (A), and π-fragments (π) for the design of new prospective molecules as organic sensitizers. In our study we use both traditional building blocks reviewed previously by us and by other authors and extended list of donor groups that includes some new previously unexplored fragments. Based on our recent investigations and the works of other researchers [7, 1221] where three variants of π-fragments (π-block), [1,2,5]thiadiazolo[3,4-c]pyridine, [1,2,5]selenadiazolo[3,4-c]pyridine, and [1,2,5]oxadiazolo[3,4-c]pyridine, showed good characteristics, they were selected for the study. Cyanoacrylic acid as an electron-acceptor, bridged with a phenyl or thienyl fused π-block, was selected as acceptors (A-block) since it is traditionally used in that role in a number of studies (e.g., [7, 22]). Frequently, substituted amines are used as donor groups and here we are also not original in their use but our list of donor groups (D-block) includes some previously unexplored fragments. Moreover, some substituents represent fused donor and acceptor fragments at the same time and thus it is possible to invent slightly another DA–π–A concept for the design of featured organic sensitizers (where DA represents a fused indissoluble agglomerate comprising an electron donor and electron-acceptor group). Also, in order to distinguish the acceptors at the right and left side of π-block we have introduced a somewhat other designation D–A–π–A′ and D–A–π–A′. All these building blocks are shown in Scheme 1. Therefore, totally 15 building blocks named from a to o are used in this work. These building blocks allow a combinatorial design of 84 different molecules listed in Table 1.

3. The Analysis of Characteristics of the Designed Organic Sensitizers

The geometry optimization along with conformational search has been performed for each molecule listed in Table 1 at the DFT B3LYP 6-311G(d,p) level of theory. Then, the HOMO and LUMO energies, excitation energies, oscillator strengths, and wavelength for the minimal excitations have been calculated using TD SCF at the DFT B3LYP 6-311G(d,p). All these characteristics are represented in Table 1. The quantum computations were performed using GAMESS software, release May 2013 R1 [23, 24].

Statistical analysis of the calculated characteristics shows that π-block plays a key role for the quantum and spectral characteristics of the designed compounds. So, the mean value of HOMO energy of [1,2,5]oxadiazolo[3,4-c]pyridine derivatives is substantially lower than that for [1,2,5]thiadiazolo[3,4-c]pyridine and [1,2,5]selenadiazolo[3,4-c]pyridine derivatives (Table 2) while the HOMO energies of [1,2,5]thiadiazolo[3,4-c]pyridine and [1,2,5]selenadiazolo[3,4-c]pyridine derivatives are statistically comparable. So, it is not surprising that the presence of an element of higher periods (S, Se) in the heterocyclic π-block reduces the HOMO level which in its turn facilitates a bathochromic shift in the UV-Vis spectra. The LUMO energies behave unusually in a sequence O, S, and Se: their mean values for [1,2,5]oxadiazolo[3,4-c]pyridine and [1,2,5]selenadiazolo[3,4-c]pyridine derivatives are statistically comparable but substantially lower than that for [1,2,5]thiadiazolo[3,4-c]pyridine derivatives (Table 2) that should provide a red shift in the UV-Vis spectra for O and Se containing π-blocks. The intervals of HOMO and LUMO energies for [1,2,5]oxadiazolo[3,4-c]pyridine, [1,2,5]thiadiazolo[3,4-c]pyridine, and [1,2,5]selenadiazolo[3,4-c]pyridine derivatives are illustrated in Figure 1. Thus, a substitution of oxygen by sulfur in a π-block yields in an increase of both HOMO and LUMO levels while the substitution by selenium leads to an increase of only HOMO energy at an approximately constant LUMO energy. Therefore, a reduction of the excitation energy and correspondingly a bathochromic shift for [1,2,5]selenadiazolo[3,4-c]pyridine derivatives summarizing the findings on the HOMO and LUMO energies can be predicted certainly.

Actually, the gaps between HOMO and LUMO and correspondingly the excitation energies for the selenium containing compounds are significantly less than that for other molecules. The LUMO and HOMO difference is less on average by 0.14–0.15 eV and the excitation energy – by 0.23 eV. At the same time, the [1,2,5]oxadiazolo[3,4-c]pyridine and [1,2,5]thiadiazolo[3,4-c]pyridine derivatives are statistically comparable in these characteristics (Table 2). Obviously, the mean value of wavelength for selenium containing compounds is greater than that for others. The bathochromic shift for selenium containing compounds relative to others is approximately 26 nm. [1,2,5]Selenadiazolo[3,4-c]pyridine derivative 51 in Table 1 (admn) possesses the maximal value of the wavelength 548 nm among the designed molecules. The minimal value of the wavelength 417 nm is observed for the molecule 47 (bcln in Table 1).

The analysis of the influence of the acceptors of A′-block shows that there are no overall effects on the whole set of the designed molecules but the effects are found separately for [1,2,5]oxadiazolo[3,4-c]pyridine and [1,2,5]selenadiazolo[3,4-c]pyridine derivatives. For [1,2,5]oxadiazolo[3,4-c]pyridine derivatives the A′-block has no significant influence on the HOMO energy, but the compounds with thienyl bridge in A′-block (o fragment) possess less value of LUMO energy than the compounds with phenyl bridge (n fragment) on average by 0.15 eV that reduces the gap between LUMO and HOMO and excitation energy and correspondingly provides a bathochromic shift in UV-Vis spectra. Actually, the mean values of the gap between LUMO and HOMO for the compounds with thienyl bridge in A′-block (o fragment) are lesser than that for the compounds with phenyl bridge (n fragment) on average by 0.14 eV, the excitation energy is lesser on average by 0.12 eV, and the wavelength is greater on average by 23 nm. This mean value for the [1,2,5]oxadiazolo[3,4-c]pyridine with thienyl bridge is 500.9 ± 9.8 nm which is comparable with the [1,2,5]selenadiazolo[3,4-c]pyridine derivatives (Table 2). [1,2,5]Oxadiazolo[3,4-c]pyridine 82 in Table 1 (gko) possesses the maximal value of wavelength 532 nm among the designed molecules.

There is no influence of A′-block observed on the quantum and spectral characteristics of [1,2,5]thiadiazolo[3,4-c]pyridine derivatives.

For [1,2,5]selenadiazolo[3,4-c]pyridine derivatives the A′-block has no significant influence on the HOMO and LUMO energies, but in contrast to [1,2,5]oxadiazolo[3,4-c]pyridine derivatives the compounds with thienyl bridge in A′-block (o fragment) possess greater gap between LUMO and HOMO than the compounds with phenyl bridge (n fragment) on average by 0.09 eV that provides an increase of excitation energy and correspondingly provides a hypsochromic shift in UV-Vis spectra. Actually, the mean values of the gap between LUMO and HOMO for the compounds with thienyl bridge in A′-block (o fragment) are greater than that for the compounds with phenyl bridge (n fragment) on average by 0.09 eV, the excitation energy is greater on average by 0.12 eV, and the wavelength is lesser on average by 16 nm. The mean value for the [1,2,5]selenadiazolo[3,4-c]pyridine derivatives with phenyl bridge is 524.6 ± 8.0 nm.

No significant dependency of the HOMO and LUMO energies, the excitation energy, and the wavelength on the quantum derived characteristics of D–A- (DA-) block and on the substitution type in the π-block (h, i, and j or k, l, and m in Scheme 1) is observed, so the calculation of more than 500 descriptors was performed using ChemoSophia on-line software (http://www.chemosophia.com/) [25] both for whole molecules (molecules 1–84 in Table 1) and for D–A- (DA-) blocks (building blocks ag in Scheme 1) in radical state.

Analysis of the descriptors showed that the HOMO energy has relationships with enthalpy of formation () of D–A- (DA-) blocks calculated within ChemoSophia Elastic Model [25] but separately for [1,2,5]oxadiazolo[3,4-c]pyridine derivatives and for [1,2,5]thiadiazolo[3,4-c]pyridine/[1,2,5]selenadiazolo[3,4-c]pyridine derivatives. The enthalpies of formation for D–A- (DA-) blocks are represented in Table 3. The relationships are shown in Figure 2 and can be described using the following polynomial equation:where is HOMO energy and , , and are coefficients.

The coefficients for [1,2,5]oxadiazolo[3,4-c]pyridine derivatives are as follows: ; ; and . The correlation coefficient ; the standard deviation .

The coefficients for [1,2,5]thiadiazolo[3,4-c]pyridine and [1,2,5]selenadiazolo[3,4-c]pyridine derivatives are as follows: ; ; and . The correlation coefficient ; the standard deviation .

Thus the values of and coefficients for [1,2,5]oxadiazolo[3,4-c]pyridine, [1,2,5]thiadiazolo[3,4-c]pyridine, and [1,2,5]selenadiazolo[3,4-c]pyridine derivatives coincide. So, the only difference in equations is the free terms which means the same form of the relationships but the curve for [1,2,5]thiadiazolo[3,4-c]pyridine and [1,2,5]selenadiazolo[3,4-c]pyridine derivatives lies above the curve for [1,2,5]oxadiazolo[3,4-c]pyridine derivatives approximately by 0.15 eV. Figure 2 and (1) show that the reduction of leads to increase of HOMO energy and in turn may provide a reduction of excitation energy and a bathochromic shift. Table 3 shows that the lower values are the characteristic feature of DA-blocks (e, f, and g) we proposed. These building blocks include alkyl and alicyclic components in their structure. Thus, the DA-blocks (e, f, and g) containing organic sensitizers should be more prospective than earlier proposed D–A-blocks (ac, bc, ad, and bd). Actually, the DA containing compounds possess less mean gap between LUMO and HOMO on average 0.12 eV, less mean excitation energy on average 0.07 eV, and greater wavelength on average 13 nm than D–A containing compounds.

Additionally, correlations between HOMO energy and dipole moment of D–A- (DA-) blocks calculated within ChemoSophia Elastic Model are observed. The greater polarity of the blocks provides the greater HOMO energy (the correlation coefficients 0.90, 0.88, and 0.88 for [1,2,5]oxadiazolo[3,4-c]pyridine, [1,2,5]thiadiazolo[3,4-c]pyridine, and [1,2,5]selenadiazolo[3,4-c]pyridine derivatives, correspondingly). The mean value of dipole moment of DA containing compounds (8.9 Debye) is greater than that for D–A containing compounds on average by 3.6 Debye.

An analysis of the oscillator strengths shows the following: [1,2,5]oxadiazolo[3,4-c]pyridine derivative 31 in Table 1 (bdho) possesses the maximal value of the oscillator strength 1.233 among the designed molecules. The minimal value of the oscillator strength 0.146 nm is observed for the molecule 72 (admo in Table 1). Generally, the largest oscillator strengths are observed for molecules with completely planar D–A-block (bc and bd fragments). The mean value of the oscillator strengths for the molecules with bc or bd fragments is 1.02 ± 0.12 while for other molecules the mean value of the oscillator strengths is 0.42 ± 0.15. Moreover, the diapasons of oscillator strengths for bc and bd containing molecules and for all other molecules do not intersect. The minimal value of oscillator strength among bc and bd containing molecules is 0.7999 (molecule 10 bdhn in Table 1) while maximal value of oscillator strength among all other molecules is 0.7386 (molecule 64, acko, in Table 1). The planar structure of the D–A-block of bc and bd containing molecules provides in turn the greater planarity and the minimal deviation from π-plane of the π-block for the whole bc and bd containing molecules. So, the standard deviation of atoms from π-plane of the π-block for the whole bc and bd containing molecules is 4.94 ± 0.14 Å while for the other molecules the standard deviation is 6.73 ± 0.42 Å. Additionally the D-block b contains the minimal number of electron-poor atoms (hydrogens) among all D- and DA-blocks which provides stronger electron donor property. Therefore, it is necessary to find more planar D–A-blocks containing minimal number of hydrogens for the design of novel prospective organic sensitizers with the maximal oscillator strength.

4. Conclusions

Thus, a molecular design and theoretical study of novel prospective organic sensitizers based on D–A–π–A concept have been fulfilled. A modification of the concept (DA–π–A) and new ways for molecular design of organic sensitizers have been proposed. As in previous researches, it has been shown that the π-block plays a key role for the wavelength shift in UV-Vis spectra; also we have shown that an existence of an element of higher periods of the Periodic System in the π-block provides a noticeable bathochromic shift. This study determines that the thienyl bridge in the A-block provides a red shift with respect to the phenyl bridge for [1,2,5]oxadiazolo[3,4-c]pyridine derivatives and vice versa the phenyl bridge in the A-block provides a red shift with respect to the thienyl bridge for [1,2,5]selenadiazolo[3,4-c]pyridine derivatives. No significant influence of the bridge in A-block was observed for [1,2,5]thiadiazolo[3,4-c]pyridine derivatives. The maximal values of wavelengths have been observed for DA-blocks with aliphatic and alicyclic fragments we proposed. It has been shown in the study that the enthalpy of formation reduction and increase of polarity of D–A- (DA-) blocks yields in bathochromic shift in UV-Vis spectra. The quantitative dependencies of wavelengths on the enthalpies of formation and dipole moments are determined. It has been shown that the maximal values of the oscillator strength are observed for the planar D and A components in D–A- (DA-) block. Therefore, it has been shown in this study that the conditions for the design of compounds possessing the maximal wavelength and the maximal oscillator strength are different: the latter presupposes the planar structure of D–A- (DA-) block while the first presupposes an existence of nonplanar aliphatic and/or alicyclic fragments in D–A- (DA-) block. Thus, the next steps for the molecular design of novel prospective organic sensitizers should include some combination of these concepts.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work was financially supported by the Russian Science Foundation (Grant no. 15-13-10022).

Supplementary Materials

The Supplementary Material includes all atomic coordinates and description of all valent bonds for all computed molecules in SDF format.

  1. Supplementary Material