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 [13]. 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 [610]. 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 [610]. Electronic first-order hyperpolarizabilities and the related NLO properties are often affected by conformational characteristics and could be helpful tools for isomeric identification [1520].

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.

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.