Research Article | Open Access

Andrea Alparone, "A Quantum Chemical Study on Structures and Electronic (Hyper)polarizabilities of 2,2′-Biselenophene Rotamers", *International Scholarly Research Notices*, vol. 2013, Article ID 251635, 5 pages, 2013. https://doi.org/10.1155/2013/251635

# A Quantum Chemical Study on Structures and Electronic (Hyper)polarizabilities of 2,2′-Biselenophene Rotamers

**Academic Editor:**P. O. Westlund

#### Abstract

Geometries, IR and Raman spectra, nucleus independent chemical shifts, and static electronic (hyper)polarizabilities of the equilibrium conformations of 2,2′-biselenophene were determined in vacuum using density functional theory (DFT) computations. At the DFT-PBE0/6-31G+pdd′ level the *antigauche* structure characterized by the dihedral angle of 157° is the global minimum, whereas the *syngauche* rotamer ( = 40°) lies ca. 0.7 kcal·mol^{−1} above the *antigauche* form. The structural and spectroscopic properties as well as the electronic polarizability of the *antigauche* are similar to those of the *syngauche* structure. On the other hand, the dipole moments and first-order hyperpolarizabilities are strongly influenced by the conformational characteristics, increasing by ca. a factor of five when passing from the *antigauche* to the *syngauche* form.

#### 1. Introduction

Oligomers and polymers based on five-membered heterocycles have received great attention as promising conductive and nonlinear optical (NLO) materials [1–3]. Although the major efforts have been principally directed towards polyfurans, polypyrroles, and especially polythiophenes, recent studies have been also dedicated to the properties of polyselenophenes [4, 5]. Selenophene, owing to the intrinsic effects of the heavy atom, is considered an interesting building-block for the design of NLO devices [6–10]. Physicochemical properties of -conjugated oligomers and polymers are usually influenced by the twisting degree of the backbone as well as by the extension of the electron delocalization [11, 12]. However, structural and electronic properties of large oligomers and polymers can be extrapolated by using data of smaller oligomers [11, 12].

Differently from the monomer, little is known about the physicochemical properties of the smallest oligoselenophenes. Structures and torsional potentials of 2,2′-biselenophene, 2,2′:5′,2′′-terselenophene, and 2,2′:5′,2′′:5′′,2′′′-quaterselenophene have been previously investigated by using *ab initio* and density functional theory (DFT) methods [13, 14]. On the basis of the most recent theoretical results obtained in vacuum, 2,2′-biselenophene is predicted to exist in two nonplanar minimum-energy conformations, characterized by dihedral angles of ca. 150° (*antigauche*) and 40° (*syngauche*) [14]. The *antigauche* is the lowest-energy conformer, with the *syngauche* being predicted to lie above the *antigauche* structure by less than 1 kcal/mol [13, 14]. The torsional potentials of 2,2′-biselenophene for the 0°–360° rotation around the bond are described by flat four-well potentials, showing a high degree of conformational flexibility [13].

The present work reports important physicochemical properties of the most stable rotamers of 2,2′-biselenophene by using computational methods. The structures of the investigated conformations are displayed in Figure 1. Specifically, we determined the geometries, vibrational spectra, nucleus independent chemical shifts, dipole moments, and static electronic polarizabilities and first-order hyperpolarizabilities in the gas phase. To the best of our knowledge experimental and theoretical hyperpolarizabilities of 2,2′-biselenophenes are not available so far, whereas there are a number of studies available for the monomer selenophene [6–10]. Electronic first-order hyperpolarizabilities and the related NLO properties are often affected by conformational characteristics and could be helpful tools for isomeric identification [15–20].

#### 2. Computational Details

The calculations were performed using the GAUSSIAN 03 program [21]. The geometries of the *antigauche* and *syngauche* forms of 2,2′-selenophene (Figure 1) were optimized using the PBE0 DFT method [22], together with the basis set previously developed by Kamada et al. [23] for structural and (hyper)polarizability computations. It consists of 6-31G basis set for hydrogen atoms, 6-31G+pdd for carbon atoms (, , and ), and LANL1DZ+pdd for selenium atoms (, , and ), which also includes effective core-potential functions. This basis is here denoted as 6-31+pdd′. The vibrational analysis performed at the PBE0/6-31+pdd′ level under the harmonic approximation shows that all the studied structures are equilibrium geometries (no imaginary frequencies). The static electronic polarizabilities () and first-order hyperpolarizabilities () were computed at the PBE0/6-31G+pdd′ level through the coupled-perturbed HF (CP-HF) theory [24]. The following orientationally invariant quantities were employed:

#### 3. Results and Discussion

The *antigauche* form is the global minimum-energy structure, characterized by the dihedral angle of 157°. The *syngauche* form ( = 40°) lies above the *antigauche* conformation by 0.68 kcal·mol^{−1}, in good agreement with previous *ab initio* and DFT results [13, 14], especially with the B3LYP/6-311++G** computations ( = 0.73 kcal·mol^{−1}) [14]. Thus on the basis of the present results, in the gas phase the *antigauche* is the prevailing structure (ca. 76%), with the *syngauche* being however an important conformation (ca. 24%). Note that, as can be appreciated from the data reported in Figure 1, the geometries of the *gauche* rotamers are very similar to each other. In particular, the bond lengths and angles of the *antigauche* differ, respectively, by no more than 0.003 Å () and 1.3° () from the corresponding values of the *syngauche* form. The structural similarity is also confirmed by the spectroscopic results. For both the rotamers we determined the nucleus independent chemical shift perpendicular to the plane of the selenophene ring at the distance of 1 Å [] [25, 26]. This parameter is closely connected to aromaticity, a negative value denoting an aromatic ring [25, 26]. The PBE0/6-31+pdd′ values for the *antigauche* and *syngauche* 2,2′-biselenophenes are very close to each other, being calculated to be −17.6 and −17.3 ppm, respectively. Note that, when passing from the monomer to the dimer, the ring aromaticity decreases by ca. 20%.

The above results on the calculated geometries are also corroborated by the IR and Raman spectra obtained at the PBE0/6-31+pdd′ level under the harmonic approximation. The simulated vibrational spectra in the 500–1800 cm^{−1} range are displayed in Figure 2. The complete sets of wavenumber values (), IR intensities (), and Raman activities () of the *antigauche* and *syngauche* forms of 2,2′-biselenophene are given in Tables 1 and 2. The most active Raman peak of the *antigauche* form placed at 1526 cm^{−1} ( Å^{4}/amu) is principally attributed to the C=C + C–C bonds stretching vibration with the nonnegligible contribution of in-plane C–H bending deformations (mode no. 9, Table 1). A graphical representation of the atom vector displacements involved in this mode is given in Figure 2(a). The corresponding transition for the *syngauche* conformation is blue-shifted by only 1 cm^{−1} (Table 2). Additionally, the IR spectrum of both rotamers is dominated by an isolated band located at 673 cm^{−1} for the *antigauche* ( km/mol, mode no. 28) and at 682 cm^{−1} for the *syngauche *( km/mol, mode no. 27). As can be seen in Figure 2(b), this transition is assigned to a pure out-of-plane bending deformation.

| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

The calculations were carried out at the PBE0/6-31 + pdd′ level. |

| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

The calculations were carried out at the PBE0/6-31 + pdd′ level. |

**(a)**

**(b)**

Beside to the structural and spectroscopic properties, now we explore the electric properties of the *gauche* rotamers. As a test case, we computed the electric properties also for the monomer, for which some experimental and high-level theoretical estimates are available from the literature. The experimental gas-phase dipole moment () of selenophene at 0.39 D [27] is overestimated by the present calculations ( D), whereas the observed value of 74.90 a.u. [28] is reasonably reproduced by the PBE0/6-31+pdd′ level ( a.u., −6%). Unfortunately, the experimental first-order hyperpolarizability of selenophene is not available to date. However, the PB0/6-31G+pdd′ value of selenophene of 98.2 a.u. is in excellent agreement with the datum previously predicted by the CCSD(T)/6-31G+pdd′ computations at 98.4 a.u. (with a deviation of −0.2%) [7].

Similarly to the monomer, both 2,2′-biselenophenes show a rather low polarity. However, due to the mutual arrangement of the selenophene moieties, when going from the *antigauche* to the *syngauche* form, the value increases by ca. a factor of five ( D versus D). A somewhat different behaviour occurs for the electronic polarizabilities, with the (*syngauche*) value of 150.76 a.u. being very close to the corresponding datum of the *antigauche* structure (152.74 a.u., +1%). By contrast, both the magnitude and direction of the vector (see Figure 1) are strongly influenced by the conformation, in line with the results obtained for the dipole moments. Specifically, the vector for the *antigauche *form is directed along the -axis, whereas for the *syngauche* rotamer it is aligned along the -axis. More importantly, the present PB0/6-31G+pdd′ calculations predict (*syngauche*) > (*antigauche*) by ca. a factor of five.

#### 4. Conclusions

The present results indicate that the calculated electronic polarizability and structural and spectroscopic properties of 2,2′-biselenophenes are little dependent on the conformation. Differently, both the magnitude and direction of the dipole moment and first-order hyperpolarizability are significantly influenced by the structural features. Therefore the *gauche* rotamers could be unambiguously discriminated on the basis of their rather different NLO signals.

#### Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

#### References

- R. D. McCullough, “The chemistry of conducting polythiophenes,”
*Advanced Materials*, vol. 10, no. 2, pp. 93–116, 1998. View at: Google Scholar - P. N. Prasad and D. J. Williams,
*Introduction to Nonlinear Optical Effects in Molecules and Polymers*, Wiley, New York, NY, USA, 1991. - J. Zyss,
*Molecular Nonlinear Optics: Materials, Physics and Devices*, Academic Press, London, UK, 1994. - A. Patra and M. Bendikov, “Polyselenophenes,”
*Journal of Materials Chemistry*, vol. 20, no. 3, pp. 422–433, 2010. View at: Publisher Site | Google Scholar - D. Bhattacharyya and K. K. Gleason, “Low band gap conformal polyselenophene thin films by oxidative chemical vapor deposition,”
*Journal of Materials Chemistry*, vol. 22, no. 2, pp. 405–410, 2012. View at: Publisher Site | Google Scholar - K. Kamada, M. Ueda, T. Sakaguchi, K. Ohta, and T. Fukumi, “Femtosecond optical Kerr study of heavy atom effects on the third-order optical nonlinearity of thiophene homologues: purely electronic contribution,”
*Chemical Physics Letters*, vol. 263, no. 1-2, pp. 215–222, 1996. View at: Google Scholar - K. Kamada, M. Ueda, H. Nagao et al., “Molecular design for organic nonlinear optics: polarizability and hyperpolarizabilities of furan homologues investigated by ab initio molecular orbital method,”
*Journal of Physical Chemistry A*, vol. 104, no. 20, pp. 4723–4734, 2000. View at: Publisher Site | Google Scholar - S. Millefiori and A. Alparone, “(Hyper)polarizability of chalcogenophenes C4H4X (X = O, S, Se, Te) conventional ab initio and density functional theory study,”
*Journal of Molecular Structure: THEOCHEM*, vol. 431, no. 1-2, pp. 59–78, 1998. View at: Google Scholar - S. Millefiori and A. Alparone, “Theoretical determination of the vibrational and electronic (hyper)polarizabilities of C
_{4}H_{4}X (X = O, S, Se, Te) heterocycles,”*Physical Chemistry Chemical Physics*, vol. 2, no. 11, pp. 2495–2501, 2000. View at: Publisher Site | Google Scholar - B. Jansik, B. Schimmelpfennig, P. Norman, P. Macak, H. Ågren, and K. Ohta, “Relativistic effects on linear and non-linear polarizabilities of the furan homologues,”
*Journal of Molecular Structure*, vol. 633, no. 2-3, pp. 237–246, 2003. View at: Publisher Site | Google Scholar - J. L. Brédas, G. B. Street, B. Thémans, and J. M. André, “Organic polymers based on aromatic rings (polyparaphenylene, polypyrrole, polythiophene): evolution of the electronic properties as a function of the torsion angle between adjacent rings,”
*Journal of Chemical Physics*, vol. 83, no. 3, pp. 1323–1329, 1985. View at: Google Scholar - V. Hernandez, C. Castiglioni, M. Del Zoppo, and G. Zerbi, “Confinement potential and
*π*-electron delocalization in polyconjugated organic materials,”*Physical Review B*, vol. 50, no. 14, pp. 9815–9823, 1994. View at: Publisher Site | Google Scholar - S. Millefiori and A. Alparone, “Theoretical investigation of the structure and conformational behaviour of small selenophene oligomers,”
*Synthetic Metals*, vol. 95, no. 3, pp. 217–224, 1998. View at: Google Scholar - G. Sánchez-Sanz, I. Alkorta, and J. Elguero, “A theoretical study of the conformation of 2,${2}^{\prime}$-bifuran, 2,${2}^{\prime}$-bithiophene, 2,${2}^{\prime}$-bitellurophene and mixed derivatives: chalcogen-chalcogen interactions or dipole-dipole effects?”
*Computational and Theoretical Chemistry*, vol. 974, no. 1–3, pp. 37–42, 2011. View at: Publisher Site | Google Scholar - G. J. M. Velders, J.-M. Gillet, P. J. Becker, and D. Feil, “Electron density analysis of nonlinear optical materials. An ab initio study of different conformations of benzene derivatives,”
*Journal of Physical Chemistry*, vol. 95, no. 22, pp. 8601–8608, 1991. View at: Google Scholar - J. Lipiński and W. Bartkowiak, “Conformation and solvent dependence of the first and second molecular hyperpolarizabilities of charge-transfer chromophores. Quantum-chemical calculations,”
*Chemical Physics*, vol. 245, no. 1-3, pp. 263–276, 1999. View at: Google Scholar - J. S. Salafsky, “Second-harmonic generation as a probe of conformational change in molecules,”
*Chemical Physics Letters*, vol. 381, no. 5-6, pp. 705–709, 2003. View at: Publisher Site | Google Scholar - A. Plaquet, B. Champagne, F. Castet et al., “Theoretical investigation of the dynamic first hyperpolarizability of DHA-VHF molecular switches,”
*New Journal of Chemistry*, vol. 33, no. 6, pp. 349–356, 2009. View at: Publisher Site | Google Scholar - Z.-M. Su, H.-L. Xu, Z.-R. Li, S. Muhammad, F. L. Gu, and K. Harigaya, “Knot-isomers of möbius cyclacene: how does the number of knots influence the structure and first hyperpolarizability,”
*Journal of Physical Chemistry C*, vol. 113, no. 34, pp. 15380–15383, 2009. View at: Publisher Site | Google Scholar - A. Alparone, “Comparative study of CCSD(T) and DFT methods: electronic (hyper)polarizabilities of glycine,”
*Chemical Physics Letters*, vol. 514, no. 1–3, pp. 21–25, 2011. View at: Publisher Site | Google Scholar - M. J. Frisch, G. W. Trucks, H. B. Schlegel et al.,
*GAUSSIAN 03, Revision E.01. Gaussian*, J. A. Pople. - C. Adamo and V. Barone, “Toward reliable density functional methods without adjustable parameters: the PBE0 model,”
*Journal of Chemical Physics*, vol. 110, no. 13, pp. 6158–6170, 1999. View at: Google Scholar - K. Kamada, M. Ueda, H. Nagao et al., “Effect of heavy atom on the second hyperpolarizability of tetrahydrofuran homologs investigated by ab initio molecular orbital method,”
*International Journal Quantum Chemistry*, vol. 70, pp. 737–743, 1998. View at: Google Scholar - H. Sekino and R. J. Bartlett, “Frequency dependent nonlinear optical properties of molecules,”
*Journal of Chemical Physics*, vol. 85, no. 2, pp. 976–989, 1986. View at: Google Scholar - Z. Chen, C. S. Wannere, C. Corminboeuf, R. Puchta, and P. von Ragué Schleyer, “Nucleus-independent chemical shifts (NICS) as an aromaticity criterion,”
*Chemical Reviews*, vol. 105, no. 10, pp. 3842–3888, 2005. View at: Publisher Site | Google Scholar - N. S. Mills and K. B. Llagostera, “Summation of nucleus independent chemical shifts as a measure of aromaticity,”
*Journal of Organic Chemistry*, vol. 72, no. 24, pp. 9163–9169, 2007. View at: Publisher Site | Google Scholar - F. Fringuelli, G. Marino, and A. Taticchi, “Tellurophene and related compounds,”
*Advances in Heterocyclic Chemistry*, vol. 21, pp. 119–173, 1977. View at: Publisher Site | Google Scholar - G. R. Dennis, I. R. Gentle, G. L. D. Ritchie, and C. G. Andrieu, “Field-gradient-induced birefringence in dilute solutions of furan, thiophen and selenophen in cyclohexane,”
*Journal of the Chemical Society, Faraday Transactions 2*, vol. 79, no. 4, pp. 539–545, 1983. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2013 Andrea Alparone. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.