We make use of a model chemistry within Density Functional Theory (DFT) recently presented, which is called CHIH-DFT, to calculate the molecular structure of 3,4-diphenyl-1,2,5-thiadiazole 1-oxide (TSO), as well to predict its infrared (IR), ultraviolet (UV-Vis) and fluorescence (Fluo) spectra, the dipole moment and polarizability, the free energy of solvation in different solvents as an indication of solubility, and the chemical reactivity parameters that arise from Conceptual DFT. The calculated values are compared with the available experimental data for this molecule as a means of validation of our proposed model chemistry.

1. Introduction

Since the dawn of history, man has been fascinated by the Sun, the provider of the light and warmth that sustains life on Earth. In preindustrial times, our major sources of energy wood, wind and water power derived from solar energy. The subsequent discovery and massive exploitation of fossil fuels laid down in the Earth's crust by early aeons of photosynthetic activity have conditioned the developed world to be dependent on convenient, readily available energy [1].

Despite much effort, semiconductor photovoltaic devices made with traditional inorganic semiconductors have remained sufficiently expensive that their uses are confined to a number of niches. Much effort is currently directed towards the use of thin-film semiconductors, in place of silicon wafers, since the direct fabrication of thin devices on substrates offers the prospect of lower manufacturing costs, particularly for larger area applications. An alternative approach is the use of organic, molecular semiconductors, which can be processed over large areas at relatively low temperatures, either by vacuum sublimation of molecular materials, or, preferably, by processing from solution of film-forming materials such as polymers [2].

The objective of this work is to perform a detailed calculation of the molecular structure of 3,4-diphenyl-1,2,5-thiadiazole 1-oxide (TSO), as well as to predict its infrared (IR), ultraviolet (UV-Vis), fluorescence (Fluo) spectra, some physical properties related to solubility and the chemical reactivity, by using a model chemistry within Density Functional Theory (DFT) [3] specially tailored to study heterocyclic systems [4] and to validate the calculated results by comparison with the available experimental data related to the molecular structure for this molecule. This compound has several desirable characteristics related to its use in Organic Photovoltaics and OLEDS: (i) it contains the 1,2,5-thiadiazole 1-oxide group, with molecular parameters similar to the thiophene 1-oxide derivatives that have been shown to be useful as nanostructured molecular materials; (ii) the -conjugated derivatives are generally efficient fluorophores, and as such, useful for the fabrication of nanobiosensors; (iii) it can be used as an attractive building block for Organic Molecular Materials.

2. Theory and Computational Details

All computational studies were performed with the Gaussian 03W [5] series of programs with density functional methods as implemented in the computational package. The equilibrium geometries of the molecules were determined by means of the gradient technique. The force constants and vibrational frequencies were determined by the FREQ calculations on the stationary points obtained after the optimization to check if there were true minima. Density functionals used in this study are a modification of those incorporated in the Gaussian 03W computational package [5]. To this end, we have defined a new model chemistry that we have found empirically that works well with heterocyclic molecules. The model was found by trial and error determination of the molecular structure of several heterocycles and a comparison with the experimental results. The implementation is a slight different version of the PBE0 hybrid density functional [6, 7]. In the PBE0 (or PBE1PBE one-parameter) functional, there is only one coefficient which is theoretically adjusted to 0.25, reflecting the mix of Hartree-Fock or exact exchange and the DFT exchange which is represented by the PBE density functional [8]. The correlation part is also represented by the PBE correlation functional [8] with coefficient equal to one. Our proposed density functional model, which we have called PBEg, is the same as PBE0, but with the mixing coefficient which adopts different values depending on the number of heteroatoms in the studied molecule, or, in turn, of its molecular structure. The value of can be calculated through the following empirical formula: , where FHA is the first heteroatom chosen as the one less electronegative, FV is valence factor which represents the oxidation state of the FHA (i.e., 1 for the first oxidation state, two for the second, and so on), and AHA are the number of additional heteroatoms besides the FHA. FHA is always equal to one. This formula has been empirically found by studying a large number of heterocyclic compounds with different number of heteroatoms and oxidation states. For example, for the 1,2,5-thiadiazole molecule, the coefficient will be , thus implying that 22% of HF exchange will be mixed with 78% of PBE exchange. Similarly, for 1,2,5-thiadiazole 1,1-dioxide, the coefficient will be , reflecting 50% of HF exchange and 50% of PBE exchange. For those cases in which the molecule has several heterocyclic rings, we will calculate the coefficient by averaging the coefficient for each heterocycle. As it is evident, the bounds to the empirical formula are 0 and 1 meaning no mixing of HF exchange in the first case and a purely HF model in the last case. The reasons for the dependence of the mix of HF and DFT exchange with the number of heteroatoms and their oxidation state are not apparent, but the validity of the empirical formula have been computationally assessed for several different molecular systems [4, 920]. Based on these results it is clear that our choice of this particular functional has been done because its performance is superior to other methods both in accuracy and computational cost.

A model chemistry is an unbiased, uniquely defined, and uniformly applicable theoretical model for predicting the properties of chemical systems. A model chemistry generally consists of the combination of a theoretical method wit a basis set [2130]. In order to define our model chemistry, we have to couple the proposed density functional with one or more basis sets. In this way, the new model chemistry that we have called CHIH-DFT can be represented by the expression CHIH = PBEg/basis sets. There are three different CHIH-DFT model chemistries: CHIH(small) that uses the 3-21G* basis set for geometry optimizations and frequency calculations, and the CBSB2** basis set for the calculation of the electronic properties; CHIH(medium) that uses the CBSB2** basis set for geometry optimizations and frequency calculations and the CBSB4 basis set for the electronic properties; and CHIH(large) which uses the CBSB7 basis set for geometry optimizations and frequency calculations and the CBSB1 basis set for the electronic properties. In this way, by considering a compromise between accuracy and CPU time, the CHIH(large) will be used for small heterocyclic molecules, the CHIH(medium) for medium-sized molecules and the CHIH(small) for large heterocyclic molecules. We understand that a molecule is small when the number of atoms within the cycle is five or six. When the number of atoms in the cycle is more than six, we considered it a medium-sized molecule. If there is more than one cycle, or if there are a large number of side-chain radical substituents, we say that the molecule is large. The CHIH-DFT model chemistry has been already validated for the case of antiparasitic drugs [911], flavonoids [1214], organic corrosion inhibitors [1517], alkaloids [18, 19] and antitubercular drugs [20]. Thus, in our case for the 3,4-diphenyl-1,2,5-thiadiazole 1-oxide molecule (TSO), the value of the mixing will be , implying that 36% of Hartree-Fock (HF) exchange is mixed with 64% of DFT exchange. The correlation part is taken as the PBE correlation functional with coefficient equal to one [8].

For the purpose of comparison, the molecular structure of the molecule has been computed using different model chemistries. Additional density functionals used in this study, are LSDA or SVWN [31, 32], BLYP [33, 34], B3LYP [34, 35], B3PW91 [35, 36], PBEPBE [8], PBE0 [7], TPSSTPSS [37] and TPSSh [38]. A suitable description of these density functionals is provided in some of the most important Computational Chemistry recent books [25, 27, 39, 40].

The calculation of the ultraviolet (UV-Vis) and fluorescence (Fluo) spectra of 3,4-diphenyl-1,2,5-thiadiazole 1-oxide (TSO) has been performed by solving the time dependent Kohn-Sham equations according to the method implemented in Gaussian 03W [39, 4143]. The equations have been solved for 10 excited states.

The infrared (IR), ultraviolet (UV-Vis), and fluorescence (Fluo) spectra were calculated and visualized using the Swizard program [44]. In all cases the displayed spectra show the calculated frequencies and absorption or emission wavelengths.

3. Results and Discussion

The results for the equilibrium conformation of the neutral molecule of 3,4-diphenyl-1,2,5-thiadiazole 1-oxide (or TSO, for short) calculated with the CHIH(small) model chemistry are presented in Figure 1 through a representation of the molecular structure showing the atomic labelling and numbering as well as the interatomic bond lengths and angles. The results have been compared with the experimental X-ray crystallography determination of the molecular structure of TSO [45]. The agreement is generally good: the standard error of the differences between the experimental and the calculated bond lengths and bond angles being very low. Although not shown in Figure 1, an analysis of the torsional angles reveals that the heterocycle is almost planar. The good agreement with the experimental data is more noticeable for the heterocyclic ring, because the model have been parameterized for this kind of systems [4].

As an additional comparison, we have calculated the optimized structure of the studied molecule, using several different common density functionals (local, nonlocal GGA, hybrid, meta-GGA and meta-GGA hybrid) coupled with the CBSB2** basis set. The molecular structure obtained from each calculation was compared with the experimental results by aligning both structures and calculating the Root Mean-Square Deviation (RMSD). The results can be summarized as follows: CHIH(small) = 0.026, PBE0 = 0.033, TPSSh = 0.049, B3LYP = 0.054, LSDA = 0.066, TPSS = 0.053, PBEPBE = 0.053, BLYP = 0.087. It can be concluded that the CHIH(small) model chemistry performs better than other hybrids functionals and even better than a recently proposed hybrid meta-GGA functional. The behavior of CHIH(small) and is similar and better than local and nonhybrids functionals. This is remarkable if we consider that this model chemistry is more efficient from a computational point of view.

The molecular dipole moment is perhaps the simplest experimental measure of charge distribution in a molecule. The accuracy of the overall distribution of electrons in a molecule is hard to quantify, since it involves all the multipoles. From the present calculations, the total energy, the total dipole moment and the isotropic polarizability of the ground state with the CHIH(small) model chemistry are 1121.415 au, 6.9576 Debye, and 66.94 Bohr for the TSO molecule. These results could be of interest as an indication of the solubility and chemical reactivity of the studied molecule, not only for it synthesis but for the potential polymerization as a step in its use in organic electronics and photovoltaics.

The infrared spectrum (IR) for the TSO molecule calculated with the CHIH(small) model chemistry is displayed in Figure 2. To the best of our knowledge, the experimental spectrum has not been reported in the literature. The vibrational band assignments have been done using the ChemCraft for Windows molecular visualization program. The symmetric vibration of the SO group has been found at 1225 cm. The symmetric stretching of the C=N double bonds is located at 1595 cm. At 1319 cm there is a signal belonging to heterocyclic C–C symmetric stretching, while at 752 cm there is a strong peak that corresponds to the vibrations of all the aromatic ring hydrogens. These are by far the most intense signals in the IR spectrum. There are two signals at 549 and 599 cm that belong to the asymmetric and symmetric stretching of the S–N bonds. At 835 cm there is a signal related to the vibrations of only three of the hydrogens of the benzenic rings. The asymmetric and symmetric out-of-plane vibrations of the aromatic ring carbons can be observed at 1070 and 1136 cm, respectively. The peak at 999 cm belongs to the vibrations of the C-N-S of the thiadiazole moiety. Finally, the symmetric and asymmetric in-plane vibrations of the aromatic ring carbons are observed at 1549 and 1675 cm.

As the experimental IR spectrum has not been reported, the scaling factor could not be calculated. However, from our experience based on previous calculations with similar molecules [4], it is possible to say that the scaling factor will be between 0.97 and 1. Thus, it is be expected that the model chemistry used in this work could reproduce the experimental spectrum of the TSO molecule with a certain degree of accuracy.

The ultraviolet spectrum (UV-Vis) of the TSO molecule calculated with the CHIH(small) model chemistry is displayed in Figure 3. The results for the energy of the excitations and the oscillator strengths computed for the ten low-lying singlet excited states are presented in Table 1. In order to characterize the excited state transitions presented in Table 1, we performed an analysis of all the molecular orbitals involved, taking into consideration that orbital 66 is the HOMO and orbital 67 is the LUMO. For example, in alpha molecular orbital 66 (occupied, ) main atomic orbitals are: N22 (Px on atom O2) , N26 (Px on atom O2) . In alpha molecular orbital 67 (unoccupied, ) main atomic orbitals are: N24 (Pz on atom O2) , N28 (Pz on atom O2) . This means that the HOMO has a character, while the LUMO has a character. On the basis of the above analysis of transitions and MOs, it can be said that the absorption bands of TSO in the visible and UV regions are typical - transitions.

As the HOMO-LUMO transition takes place in the ultraviolet region, close to but out of the visible zone, it can be predicted that this molecule will be colorless or slightly colored. In order to calculate the fluorescence spectrum of the studied molecule, the first singlet excited state was optimized by using the CIS/3-21G(d) model chemistry. The results of a TDDFT calculation at the CHIH(small) model chemistry on the optimized singlet excited state are not shown but reported as the fluorescent spectrum, where the HOMO-LUMO transition takes place at 613 nm, being this signal very weak. This means that the studied molecule will present fluorescence in the visible region. Thus, it could be of potential application for the fabrication of OLEDs.

The free energy of solvation G(solv) of the molecule have been calculated for TSO by resorting to the CHIH(small) model chemistry coupled with the Integral-Equation-Formalism of the Polarized Continuum Model (IEF-PCM) for different solvents as implemented in Gaussian 03. The solubility of a molecule will depend on several kinetic and thermodynamic factors. However, it can be said that the magnitude and the sign of G(solv) could be a good approximation as an index of solubility. In this way, a negative sign and a large magnitude will be an indication of increased solubility. The results of these calculations for the studied molecule can be summarized as follows: Chloroform = 4.56 Kcal/mol, Cyclohexane = 2.75 kcal/mol, Ethanol = 10.30 Kcal/mol, THF = 5.18 Kcal/mol, Toluene = 2.33 Kcal/mol, and Water = 5.20 Kcal/mol. These values could be an indication that the studied molecule will be mostly soluble in ethanol, and this can be related to the results obtained for the dipole moment and polarizability.

The HOMO and LUMO of TSO calculated with the CHIH(small) model chemistry are displayed in Figure 4. The HOMO density corresponds to the nitrogens, the oxygen and some carbon atoms of the nonheterocyclic rings, while the LUMO density is over the whole molecule. This can give us an idea of the reactivity of the molecule.

The reactive sites can be identified through an analysis of the total and orbital densities. The representation of the calculated HOMO and LUMO densities in Figure 4 show that the electrophilic attack would occur preferentially at the nitrogens and the nucleophilic attack at the heterocyclic ring carbons.

The condensed Fukui functions can also be employed to determine the reactivity of each atom in the molecule. The corresponding condensed functions are given by (for nucleophilic attack), (for electrophilic attack), and (for radical attack), where is the gross charge of atom in the molecule.

It is possible to evaluate condensed Fukui functions from single-points calculations directly, without resorting to additional calculations involving the systems with and electrons: with being the LCAO coefficients and the overlap matrix. The condensed Fukui functions are normalized, thus and .

The results from the calculation of the condensed Fukui functions for nucleophilic, electrophilic, and radical attack have been obtained by resorting to the AOMix molecular analysis program [46] and are presented in Table 2.

The sites for electrophilic attack will be those atoms bearing a negative charge and where the Fukui function is a maximum. These values confirm that the site for the electrophilic attack is the oxygen atom. The site for potential nucleophilic attack would depend on the values of on the atoms with a positive charge density. The results form Table 1 show that the site for nucleophilic attack will be the nitrogen atoms and also the heterocyclic carbon atoms. Finally, the sites for radical attack, governed by the values of will be the oxygen atom, and perhaps the sulfur atom. The benzenic rings are not very reactive in this context.

4. Conclusions

In this work, our model chemistry within DFT (the CHIH-DFT chemistry model) has been presented and the methodology has been applied to the study of a molecule which is potentially useful for OLEDs and photovoltaic devices. The molecular structure for TSO has been determined by using the CHIH(small) model chemistry. A comparison has been made with the results from the experimental X-ray crystallography for this molecule. The agreement is generally very good. It is worth noting that the results obtained with the (faster) CHIH(small) model chemistry are very similar to those obtained using the PBE0 functional with a CBSB7 basis set, thus implying that equally accurate results could be obtained in a reasonable time in the study of heterocycles. A general comparison with the results from the optimization of the molecular structure using several different density functionals have also been performed and the conclusion is that the CHIH-DFT is a good alternative to other model chemistries based on hybrids GGA and meta-GGA functionals for the prediction of molecular structures.

The shape of the frontier orbitals of this molecule were displayed as well as some electronic parameters like the total energy, the dipole moment and the polarizability. Moreover, the calculation of the condensed Fukui functions can be useful to obtain an idea of the chemical reactivity of the studied molecule. This is of special importance if we consider that the studied molecule could be the starting point for a polymeric material of potential application in organic electronics and photovoltaics, and that the chemical reactivity indices could assist the experimentalists in planing and guiding the synthetic and polymerization processes.

The infrared (IR) and ultraviolet (UV-Vis) and fluorescence (Fluo) spectra for TSO have been predicted according to the CHIH(small) model chemistry, and an assignment of the principal peaks have been achieved. The shape of the UV-Vis spectra and the maximum absorption wavelength belonging to this molecule have been presented and this have been compared with the fluorescence spectra. The results showed that the studied molecule could be an useful starting point for the preparation of materials of use in organic electronic and photovoltaics.

The free energy of solvation G(solv) of the molecule have been calculated for TSO by resorting to the CHIH(small) model chemistry coupled with the Integral-Equation-Formalism of the Polarized Continuum Model (IEF-PCM) for different solvents and the results gave an indication of ethanol as the solvent in which this molecule could be potentially soluble.

The CHIH-DFT model chemistry appears to be a useful tool for the study of the molecular structure and electronic properties of heterocycles, and further applications to several molecular systems are being pursued in our laboratory.


This work has been partially supported by Consejo Nacional de Ciencia y Tecnologa (CONACYT, Mexico) and by Fondo Mixto del Estado de Baja California (FOMIX-BC) through Project 69363. The author is a researcher of CONACYT and CIMAV.