Second Harmonic Generation, Electrooptical Pockels Effect, and Static First-Order Hyperpolarizabilities of 2,2[variant prime]-Bithiophene Conformers: An HF, MP2, and DFT Theoretical Investigation

Abstract

The static and dynamic electronic (hyper)polarizabilities of the equilibrium conformations of 2,2[variant prime]-bithiophene (anti-gauche and syn-gauche) were computed in the gas phase. The calculations were carried out using Hartree-Fock (HF), Møller-Plesset second-order perturbation theory (MP2), and density functional theory methods. The properties were evaluated for the second harmonic generation (SHG), and electrooptical Pockels effect (EOPE) nonlinear optical processes at the typical λ=1064 nm of the Nd:YAG laser. The anti-gauche form characterized by the S-[subscript]C2[/subscript] -[subscript]C[superscript]2[variant prime][/superscript] [/subscript] -S dihedral angle of 137° (MP2/6-311G**) is the global minimum on the potential energy surface, whereas the syn-gauche rotamer (S-[subscript]C2[/subscript] -[subscript]C[superscript]2[variant prime][/superscript] [/subscript] -S = 48°, MP2/6-311G**) lies ca. 0.5 kcal/mol above the anti-gauche form. The structural properties of the gauche structures are rather similar to each other. The MP2 electron correlation effects are dramatic for the first-order hyperpolarizabilities of the 2,2[variant prime]-bithiophenes, decreasing the HF values by ca. a factor of three. When passing from the anti-gauche to the syn-gauche conformer, the static and frequency-dependent first-order hyperpolarizabilities increase by ca. a factor of two. Differently, the electronic polarizabilities and second-order hyperpolarizabilities of these rotamers are rather close to each other. The syn-gauche structure could be discriminated from the anti-gauche one through its much more intense SHG and EOPE signals.

(ProQuest: ... denotes non-US-ASCII text omitted.)

Academic Editor:Miquel Solà

1, Department of Chemistry, University of Catania, Viale A. Doria 6, 95125 Catania, Italy

Received 29 September 2013; Accepted 30 October 2013

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.

1. Introduction

Thiophene-based oligomers and polymers are an interesting class of π -conjugated materials for the development and construction of conductive and nonlinear optical (NLO) devices [1-5]. Electronic properties of π -conjugated polymeric systems are significantly affected by the twisting degree of the backbone and extension of the electron delocalization although molecular structure and physicochemical properties of extended oligomers and polymers are usually modelled through smaller oligomeric chains [6, 7].

The C_α -C_{α '} bonded bithiophene oligomer, 2,2[variant prime]-bithiophene (Figure 1), is the principal building block of polythiophene chains, extensively characterized by experimental and theoretical studies. In the solid state, 2,2[variant prime]-bithiophene predominantly exists as a planar anti -structure together with a nonnegligible fraction of planar syn -conformation (ca. 15%) [8]. A slightly different picture occurs in the gas phase: two nonplanar minimum-energy structures are observed to coexist (Figure 1), which are characterized by S-C₂ -C_2' -S dihedral angles of 148° ± 3° (anti-gauche ) and 36° ± 5° (syn-gauche ) [9]. On the basis of electron diffraction data [9] and experimental fluorescence spectra [10], the anti-gauche is the minimum-energy form on the potential energy surface of the 2,2[variant prime]-bithiophene, being more stable than the syn-gauche rotamer by 1.16±0.13 kcal/mol (enthalpy difference extrapolated from the dependence of the relative abundance in the 58-130°C range of temperatures) [10] and 0.18 kcal/mol at 100°C (energy difference) [9]. On the theoretical side, there are many contributions in the literature about relative stabilities and torsional potentials of 2,2[variant prime]-bithiophene conformations [11-15]. In agreement with experiment, correlated ab initio and density functional theory (DFT) levels concordantly predict both the gauche structures as stationary points on the PESs of 2,2[variant prime]-bithiophene, the anti-gauche being predicted to be the global minimum [11-15]. The torsional potentials for the rotation around the C₂ -C_2' bond are characterized by flat double-well potentials, allowing a high degree of conformational flexibility to oligothiophenes [11-15]. In addition, the influence of the torsional potential on the electronic polarizabilities of 2,2[variant prime]-bithiophene rotamers was analyzed by Lukes et al. [14] using HF/aug-cc-pVDZ and MP2/aug-cc-pVDZ computations. However, the average electronic polarizabilities of 2,2[variant prime]-bithiophene are little affected by the conformation, varying within 5-6% [14]. Therefore, this response electric property is little informative for identification of conformations. On the other hand, electronic first-order hyperpolarizabilities (β ), the related NLO electrooptical Pockels effect (EOPE), and second harmonic generation (SHG) properties are usually much more influenced by structural characteristics than electronic polarizabilities, being of potential utility for discrimination of conformers [16-19].

Figure 1: Structures of thiophene and 2,2[variant prime]-bithiophene rotamers. Colours: white (hydrogen), grey (carbon), and orange (sulphur).

[figure omitted; refer to PDF]

In the present investigation, we have calculated the static and dynamic electronic (hyper)polarizabilities of the most stable conformations of 2,2[variant prime]-bithiophene in the gas phase, aiming to identify physicochemical properties useful to discriminate the different rotamers. The frequency-dependent first-order hyperpolarizabilities were obtained for the SHG and EOPE NLO processes at the typical experimental wavelength of the Nd:YAG laser (1064 nm). Neither experimental nor theoretical electronic β values of 2,2[variant prime]-bithiophenes are known so far, whereas some computational studies were previously carried out on the first-order hyperpolarizabilities of the monomer of thiophene [20-25].

2. Computational Details

The geometries of the anti-gauche and syn-gauche 2,2[variant prime]-bithiophene rotamers (Figure 1) were optimized using ab initio Hartree-Fock (HF) and Møller-Plesset second-order perturbation theories (MP2) with the 6-311G** basis set. The vibrational analysis obtained under the harmonic approximation confirmed that all the investigated structures are stationary points (no imaginary frequencies). The calculations were also carried out on the thiophene molecule for comparison.

For thiophene and 2,2[variant prime]-bithiophenes, we computed the dipole moments (μ ), static and dynamic electronic polarizabilities (α ), first-order hyperpolarizabilities (β ), and second-order hyperpolarizabilities (γ ). To this purpose, the HF and MP2 levels as well as a series of hybrid DFT methods such as B3LYP [26, 27], PBE0 [28], BH&HLYP [26], and B97-1 [29] were employed. Additionally, we investigated the performances of the long-range corrected ω B97X-D functional [30], recently employed with success for response electric property calculations [31-34]. The present computations were thoroughly carried out using the polarised and diffuse Sadlej's POL basis set [35]. There are many indications in the literature showing that this basis set is adequate for predicting response electric properties of organic compounds [31, 36-39]. However, for thiophene as a test case, we also performed (hyper)polarizability calculations using the larger correlation-consistent Dunning's triple-zeta basis set (aug-cc-pVTZ) [40]. At the HF level, the static α and β values were determined analytically by means of the coupled-perturbed HF (CP-HF) theory [41, 42], whereas the MP2 and DFT α and β data were obtained through a finite-field (FF) numerical scheme illustrated in detail by Kurtz and coworkers [43]. For the FF computations, we used a field strength amplitude of 0.005 a.u.. The accuracy of the numerical procedure was verified at the HF level by comparing the FF-HF and CP-HF (hyper)polarizability values. The static γ values were determined at the HF/POL level by means of the FF procedure. Frequency-dependent polarizabilities [α(-ω;ω) ] and first-order hyperpolarizabilities were computed for the 2,2[variant prime]-bithiophenes by using the CP-HF procedure at the characteristic Nd:YAG laser wavelength (λ ) of 1064 nm (...ω=0.04282 a.u.). Specifically, the SHG [β(-2ω;ω,ω) ] and EOPE [β(-ω;ω,0) ] NLO processes were investigated. At this λ value, resonance enhancement effects for the SHG phenomenon are expected to be rather negligible, since the 2...ω value of 0.08564 a.u. (2.33 eV) is rather far from the lowest-energy absorption of planar 2,2[variant prime]-bithiophene which is observed in gas at 3.86 eV [44, 45].

As commonly used in the literature, the calculated physicochemical properties are here expressed as dipole moment (μ ), average polarizability (Y9;αYA; ), polarizability anisotropy (Δα ), first-order hyperpolarizability aligned along the μ direction (_βμ ), and average second-order hyperpolarizability (Y9;γYA; ), which are orientationally invariant quantities: [figure omitted; refer to PDF] where _βi (i=x,y,z) is given by _βi =(1/3)_∑j=x,y,z (_βijj +_βjij +_βjji ) [figure omitted; refer to PDF] For the (hyper)polarizabilities, atomic units are used throughout the work. Conversion factors to the SI are 1 a.u. of α (e^2a02Eh-1 ) = 1.648778 × 10^-41 C² m² J^-2 ; 1 a.u. of β (e^3a03Eh-2 ) = 3.206361 × 10^-53 C³ m³ J^-2 ; 1 a.u. of γ (e^4a04Eh-3 ) = 6.235377 × 10^-65 C⁴ m⁴ J^-3 . All calculations were performed using the GAUSSIAN 09 program [46].

3. Results and Discussion

3.1. Geometries and Energetics of Thiophene and 2,2[variant prime]-Bithiophene Rotamers

The structural parameters of the equilibrium conformations of 2,2[variant prime]-bithiophene (anti-gauche and syn-gauche ) as well as those of the monomer obtained in the gas phase at the MP2/6-311G** level are displayed in Figure 2. On the basis of the MP2/6-311G** calculations, the anti-gauche form is the global minimum, while the syn-gauche form is less stable by 0.51 kcal·mol^-1 , in excellent agreement with previous correlated ab initio results [11-15] and in particular with the highest CCSD(T)/6-31G** and CCSD(T)/cc-pVDZ levels (0.49 kcal·mol^-1 ) [13]. Thus following the MP2/6-311G** results, the anti-gauche rotamer is the prevailing form of 2,2[variant prime]-bithiophene in vacuum (ca. 70%), the syn-gauche being, however, an important conformation (ca. 30%).

Figure 2: Gas phase MP2/6-311G** geometrical parameters of thiophene and 2,2[variant prime]-bithiophene anti-gauche and syn-gauche rotamers. The experimental data reported in parentheses are taken from [47] (thiophene) and [9] (2,2[variant prime]-bithiophene). The MP2/6-311G** (experimental) S-C₂ -C_2' -S dihedral angles are 137° (148°, see [9]) and 48° for the anti-gauche and syn-gauche conformation, respectively.

[figure omitted; refer to PDF]

The MP2/6-311G** geometry of thiophene is in good agreement with the experimental gas phase microwave structure [47]. In particular, the experimental C-S, C-C, and C=C bond lengths are well reproduced by the present computations. As can be appreciated from the data reported in Figure 2, the agreement between the observed [9] and calculated geometries for the 2,2[variant prime]-bithiophene is slightly less satisfactorily. In particular, the present computations underestimate the experimental typical S-C₂ -C_2' -S dihedral angle (148°) by 11° and the C-C interring bond length (1.46 Å) by ca. 0.01 Å. However, the MP2/6-311G** S-C₂ -C_2' -S dihedral angles of 137° (anti-gauche rotamer) and 48° (syn-gauche rotamer) are in reasonable agreement with the CCSD(T)/6-31G* estimates (142° and 44°, resp.) [13]. Some additional observations can be made: (i) the equilibrium geometry of the thiophene ring shows only little differences when passing from the monomer to the dimer; (ii) the geometries of the gauche conformations are very similar to each other. Indeed, the bond lengths and angles of the anti-gauche rotamer differ, respectively, by no more than 0.01 Å and 1° from the data of the syn-gauche form.

3.2. Effects of the Basis Set and Theoretical Level on the Electronic (Hyper)polarizabilities of Thiophene

As widely documented in the literature, accurate estimates of electronic (hyper)polarizabilities need polarized and diffuse basis sets as well as introduction of electron correlation contributions [48-52]. In particular for the series of thiophene (C₄ H₄ X, X = O, S, Se, Te) [23] and pyrrole (C₄ H₄ YH, Y = N, P, As, Sb, Bi) [53] homologues, the correlated MP2 method reproduces reasonably well the response electric properties obtained using high-level CCSD(T) and MP4-SDTQ calculations. In the current study, we explored the effects of the basis set and computational method on the μ , α , and β values of thiophene as a case test, for which some experimental and high-level theoretical data are known. Specifically, we compared the POL and the aug-cc-pVTZ basis sets. The POL basis set consists of [3s2p] functions for hydrogen atom, [5s3p2d] for carbon and [7s5p2d] for sulphur, giving a total of 164 basis functions on thiophene, which are ca. half of the functions for the aug-cc-pVTZ basis set ([4s3p2d] for H, [5s4p3d2f] for C, and [6s5p3d2f] for S). The calculations were performed at the HF level and the data are collected in Table 1. The results show that, when passing from the HF/POL to HF/aug-cc-pVTZ level, only marginal variations are observed on the calculated properties. Indeed, the μ , Y9;αYA; , Δα , and β_μ values vary by 1.43, 0.22, 0.04, and 10.41%, respectively. It is worth mentioning that the above results are in line with previous (hyper)polarizability studies using the two basis sets [22, 31, 54]. However, it is important to notice that for thiophene the HF/aug-cc-pVTZ hyperpolarizability computations require a CPU demand by one order of magnitude greater than that for the HF/POL calculations. Thus, the use of the POL basis set can be considered a valid compromise between accuracy and computational cost and will be entirely adopted for the subsequent MP2 and DFT calculations on thiophene and on the 2,2[variant prime]-bithiophenes.

Table 1: Dipole moments (D), static electronic polarizabilities (a.u.), and first-order hyperpolarizabilities (a.u.) of thiophene^a .

	HFb	MP2	BH and HLYP	PBE0	B3LYP	B97-1	ω B97X-D	MP4-SDTQc	Exp.
μ z	0.69 (0.70)	0.42	0.53	0.49	0.46	0.48	0.52	0.70	0.55d
α x x	68.61 (68.69)	71.49	7019	71.04	72.18	71.86	70.63	68.32
α y y	44.37 (44.51)	45.17	43.81	43.77	44.48	44.35	44.25	44.61
α z z	75.81 (76.00)	76.84	76.43	76.95	78.02	77.68	76.79	75.50
Y9; α Ya;	62.93 (63.07)	64.50	63.48	63.92	64.89	64.62	63.89	62.81	60.71e ; 64.90f ; 65.17g ; 66.12h
Δ α	28.52 (28.54)	29.36	29.99	30.66	31.03	30.83	29.94	27.99	23.25e ; 31.90f ; 21.56g
β x x z	-5.2 (-1.7)	-1.7	-1.2	1.7	2.2	2.9	-0.8	-5.1
β y y z	6.2 (7.3)	9.4	9.9	12.8	12.6	12.3	9.7	10.9
β z z z	40.4 (40.5)	4.5	28.3	26.8	25.1	26.8	30.8	7.7
β μ	41.3 (46.1)	12.2	37.0	41.3	39.9	42.0	39.7	13.5

^a The calculations were performed with the POL basis set on the MP2/6-311G^** geometries.

^b The values in parentheses refer to the HF/aug-cc-pVTZ computations.

^c Values taken from [23]. Basis set: 6-31G (3d) + pd.

^d Reference [55].

^e Reference [56].

^f Reference [57].

^g Reference [58].

^h Reference [59].

The effects of electron correlation as evaluated at the MP2/POL level are rather significant for the dipole moment but negligible for the polarizabilities, decreasing the HF/POL μ datum by ca. 0.3 D (-40%) and increasing the Y9;αYA; value by ca. 2.5% and the Δα value by ca. 2.9%. All the DFT methods give similar μ and α values to each other, being in reasonable agreement with the MP2/POL data. The observed gas phase dipole moment of 0.55 D [55] is underestimated by the MP2 and DFT calculations (by 0.03-0.13 D), the smallest deviation being obtained by the ω B97X-D functional. The experimental Y9;αYA; (Δα ) data, comprised between 61 and 66 a.u. (23 and 32 a.u.) [56-59], are reasonably well reproduced by all the present theoretical methods, including the HF level. More importantly, the introduction of electron correlation contributions is crucial for the first-order hyperpolarizability. In fact, when passing from the HF/POL to MP2/POL level, the β_μ value reduces by ca. a factor of three. The effect is especially conspicuous for the _βzzz component, which decreases by one order of magnitude. Note that our HF versus MP2 comparison for thiophene agrees with those previously obtained with other basis sets [22, 23]. Unfortunately, the experimental first-order hyperpolarizability of thiophene is not available so far. However, it is of interest mentioning that the present MP2/POL β_μ (thiophene) value is in good agreement with the datum previously predicted by the high-level MP4-SDTQ calculations (β_μ = 13.5 a.u.) [23], showing a difference of 1.3 a.u. (-9.6%). On the other hand, similar to the HF/POL computations, all the DFT methods overestimate the MP2/POL β_μ value by a factor between 3.0 and 3.4, the BH&HLYP functional giving the closest value. Note that the failure of the traditional DFT methods in the prediction of the electronic (hyper)polarizabilities especially of π -conjugated compounds is well known and has been exhaustively illustrated by Champagne and co-workers [60]. Quite surprisingly, it is worth noting that the use of the long-range corrected ω B97X-D functional does not improve significantly the performances obtained using the conventional functionals.

3.3. Static and Frequency-Dependent (Hyper)polarizabilities of 2,2[variant prime]-Bithiophene Rotamers

Table 2 lists the dipole moments and static (hyper)polarizabilities of the anti-gauche and syn-gauche 2,2[variant prime]-bithiophene forms calculated at the HF/POL, MP2/POL, and BH&HLYP/POL levels. By analogy to the monomer, both the 2,2[variant prime]-bithiophene forms exhibit a somewhat low polarity. Nevertheless, when passing from the syn-gauche to the anti-gauche conformer, the μ value decreases by ca. a factor of two, due to the mutual disposition of the monomeric thiophene rings. Note that, although the syn-gauche conformation reveals a rotated structure, its μ value is slightly less than 2 × μ (thiophene). On the other hand, in the case of the polarizabilities, the variations between the studied rotamers are significantly minor. In fact, the static MP2/POL Y9;αYA; value for the syn-gauche form is only 0.51 a.u. smaller than the corresponding value for the anti-gauche conformation (-0.4%). For both the dimers, _αxx is the largest component, recovering ca. 50% of the total polarizability (_αxx + _αyy + _αzz ). The MP2/POL Y9;αYA; data overestimate the experimental datum obtained in tetrahydrofuran solution by ca. 10 a.u. (+7.7%) [59]; however, they are in good agreement with those previously computed at the MP2/aug-cc-pVDZ level for conformations in the 0-180° S-C₂ -C_2' -S dihedral angle range [14]. Note that the present Δα values are slightly more influenced by the structure, increasing by ca. 3 a.u. (+3.6%) when passing from the syn-gauche to the anti-gauche rotamer. By analogy to the results found for the monomer, in comparison to the MP2/POL data, the HF/POL and BH&HLYP/POL methods furnish good Y9;αYA; and Δα estimates (within 2.5-4.4% and 0.4-4.6%, resp.). Similarly to the electronic polarizabilities, the Y9;γYA; values are little dependent on the conformation, decreasing by ca. 5% when going from the anti-gauche to the syn-gauche form. In addition, as for the calculated polarizabilities, the dominant γ component lies along the x -axis for both the gauche structures, amounting to ca. 50% of the total second-order hyperpolarizability (_γxxxx + _γyyyy + _γzzzz + _2γxxyy + _2γxxzz + _2γyyzz ).

Table 2: Dipole moments (D), static electronic polarizabilities (a.u.), first-order hyperpolarizabilities (a.u.), and second-order hyperpolarizabilities (a.u.) of 2,2[variant prime]-bithiophene rotamers^a .

	anti-gauche	syn-gauche
	HF	MP2	BH and HLYP	HF	MP2	BH and HLYP
μ	0.52	0.31	0.40	1.21	0.72	0.93
α x x	178.53	184.76	186.63	177.29	183.13	184.95
α y y	84.85	86.72	84.30	125.80	127.96	126.58
α z z	127.99	129.81	128.68	86.73	88.65	86.27
Y9; α Ya;	130.46	133.76	133.20	129.94	133.25	132.60
Δ α	81.71	85.28	89.22	78.70	82.29	85.98
β x x y	2.7	- 3.1	7.0	0.0	0.0	0.0
β y y y	29.9	30.3	32.2	0.0	0.0	0.0
β z z y	16.1	- 3.3	7.0	0.0	0.0	0.0
β x x z	0.0	0.0	0.0	45.7	31.5	59.0
β y y z	0.0	0.0	0.0	22.5	- 4.0	27.0
β z z z	0.0	0.0	0.0	40.8	16.3	21.9
β μ	48.7	23.9	46.2	109.0	43.8	107.9
γ x x x x	102128			90352
γ y y y y	21824			22464
γ z z z z	21392			21216
γ x x y y	13312			12960
γ x x z z	6432			7744
γ y y z z	7536			6928
Y9; γ Ya; b	39981			37859

^a The calculations were performed with the POL basis set on the MP2/6-311G** geometries.

^b The HF/POL Y9;γYA; value of the thiophene monomer is calculated to be 14850 a.u..

The dispersion effects (Table 3) evaluated at ...ω=0.04282 and 0.08564 a.u. increase the static Y9;αYA; (Δα ) values of both the conformations, respectively, by ca. 2 and 7% (4 and 18%). It is of interest noting that the Y9;αYA; and Y9;γYA; data of the 2,2[variant prime]-bithiophenes are greater than twice the corresponding values for the monomer, suggesting a some degree of interring π -conjugation in the dimers.

Table 3: Static and frequency-dependent (...ω=0.04282 a.u.) electronic polarizabilities (a.u.) and first-order hyperpolarizabilities (a.u.) of 2,2[variant prime]-bithiophene rotamers^a .

	anti-gauche	syn-gauche
Y9; α Ya; (0;0 )	130.46	129.94
Y9; α Ya; (-ω;ω )	132.57 (140.15)	131.97 (139.14)
Δ α (0;0 )	81.71	78.70
Δ α (-ω;ω )	84.68 (96.33)	81.43 (92.48)
β μ (0;0 )	48.7	109.0
β μ (-ω;ω;0 )	51.3	116.1
β μ (-2ω;ω;ω )	54.7	126.8

^a The calculations were performed at the HF/POL level on the MP2/6-311G** geometries. The values in parentheses refer to the frequency-dependent polarizabilities at ...ω=0.08564 a.u..

Differently from the calculated polarizabilities and second-order hyperpolarizabilities, both the magnitude (Table 1) and direction (Figure 1) of the _βμ vector are noticeably affected by the structural features, almost following the behaviour found for the dipole moments. In fact, for the anti-gauche form, the _βμ property which is aligned along the y -axis is predicted to be smaller than that for the syn-gauche rotamer (with _βμ aligned along the z -axis) by ca. a factor of two. This result is almost independent from the theoretical level and is observed for the static first-order hyperpolarizabilities as well as for the NLO SHG and EOPE processes. As for the monomer, the introduction of the MP2 electron correlation contributions is negative for _βμ , reducing the HF/POL data by 50-60%, the largest effect being found for the syn-gauche rotamer. Similar to the HF behavior, the BH&HLYP functional overestimates the MP2/POL _βμ values by 95-165%. It is of interest to note that, at the HF/POL level when passing from the monomer to the syn-gauche 2,2[variant prime]-bithiophene, the static _βμ value increases by a factor of 2.6, although the dimer exhibits a nonplanar arrangement. Note that the corresponding monomer[arrow right]dimer increases obtained at the BH&HLYP/POL and MP2/POL levels are greater, being calculated to be 2.9 and 3.6, respectively.

Not surprisingly, for both the rotamers, _βμ (SHG) > _βμ (EOPE) > _βμ (static) (Table 3). The dispersion effects estimated at ...ω=0.04282 a.u. increase the static _βμ values by 5-6% for the EOPE and 12-16% for the SHG NLO phenomenon, the largest percentages being predicted for the syn-gauche conformation.

4. Conclusions

The dipole moments and static and frequency-dependent electronic (hyper)polarizabilities of the anti-gauche and syn-gauche minimum-energy conformations of 2,2[variant prime]-bithiophene were studied in the gas phase using ab initio HF, MP2, and DFT methods. The NLO properties for the SHG and EOPE phenomena were explored at λ=1064 nm. The effects of electron correlation at MP2 level are remarkable especially for the first-order hyperpolarizabilities, reducing the HF data by 50-60%. The DFT methods, although furnishing good performance for the dipole moments and polarizabilities, significantly overestimate the MP2 first-order hyperpolarizabilities. The polarizabilities and second-order hyperpolarizabilities are little influenced by the conformation. By contrast, both the magnitude and direction of the dipole moment and first-order hyperpolarizabilities are strongly affected by the structural characteristics, the magnitudes increasing when going from the anti-gauche to the syn-gauche form by ca. a factor of two. On the basis of the present findings, the 2,2[variant prime]-bithiophene rotamers might be identified through experimental SHG and EOPE NLO measurements.

Conflict of Interests

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

References

[1] J. Roncali, "Conjugated poly(thiophenes): synthesis, functionalization, and applications," Chemical Reviews , vol. 92, no. 4, pp. 711-738, 1992.

[2] J. Zyss Molecular Nonlinear Optics: Materials, Physics and Devices , Academic Press, London, UK, 1994.

[3] R. D. McCullough, "The chemistry of conducting polythiophenes," Advanced Materials , vol. 10, no. 2, pp. 93-116, 1998.

[4] D. Fichou Handbook of Oligo- and Polythiophenes , Wiley-VCH, Weinheim, Germany, 1999.

[5] A. Mishra, C.-Q. Ma, P. Bäuerle, "Functional oligothiophenes: molecular design for multidimensional nanoarchitectures and their applications," Chemical Reviews , vol. 109, no. 3, pp. 1141-1176, 2009.

[6] V. Hernandez, C. Castiglioni, M. Del Zoppo, G. Zerbi, "Confinement potential and π -electron delocalization in polyconjugated organic materials," Physical Review B , vol. 50, no. 14, pp. 9815-9823, 1994.

[7] J. L. Brédas, G. B. Street, B. Thémans, 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," The Journal of Chemical Physics , vol. 83, no. 3, pp. 1323-1329, 1985.

[8] P. A. Chaloner, S. R. Gunatunga, P. B. Hitchcock, "Redetermination of 2, 2[variant prime]-bithiophene," Acta Crystallographica C , vol. 50, pp. 1941-1942, 1994.

[9] S. Samdal, E. J. Samuelsen, H. V. Volden, "Molecular conformation of 2,2[variant prime]-bithiophene determined by gas phase electron diffraction and ab initio calculations," Synthetic Metals , vol. 59, no. 2, pp. 259-265, 1993.

[10] J. E. Chadwick, B. E. Kohler, "Optical spectra of isolated s-cis- and s-trans-bithiophene: torsional potential in the ground and excited states," Journal of Physical Chemistry , vol. 98, no. 14, pp. 3631-3637, 1994.

[11] A. Karpfen, C. H. Choi, M. Kertesz, "Single-bond torsional potentials in conjugated systems: a comparison of ab initio and density functional results," Journal of Physical Chemistry A , vol. 101, no. 40, pp. 7426-7433, 1997.

[12] S. Millefiori, A. Alparone, A. Millefiori, "Conformational properties of thiophene oligomers," Journal of Heterocyclic Chemistry , vol. 37, no. 4, pp. 847-853, 2000.

[13] G. Raos, A. Famulari, V. Marcon, "Computational reinvestigation of the bithiophene torsion potential," Chemical Physics Letters , vol. 379, no. 3-4, pp. 364-372, 2003.

[14] V. Lukes, M. Breza, S. Biskupic, "Structure and electronic properties of bithiophenes. I. Torsional dependence," Journal of Molecular Structure , vol. 618, no. 1-2, pp. 93-100, 2002.

[15] M. A. V. Ribeiro Da Silva, A. F. L. O. M. Santos, J. R. B. Gomes, M. V. Roux, M. Temprado, P. Jiménez, R. Notario, "Thermochemistry of bithiophenes and thienyl radicals. A calorimetric and computational study," Journal of Physical Chemistry A , vol. 113, no. 41, pp. 11042-11050, 2009.

[16] 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.

[17] G. J. M. Velders, J.-M. Gillet, P. J. Becker, 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.

[18] E. Hendrickx, K. Clays, A. Persoons, C. Dehu, J. L. Brédas, "The bacteriorhodopsin chromophore retinal and derivatives: an experimental and theoretical investigation of the second-order optical properties," Journal of the American Chemical Society , vol. 117, no. 12, pp. 3547-3555, 1995.

[19] J. Lipinski, 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.

[20] V. Keshari, W. M. K. P. Wijekoon, P. N. Prasad, S. P. Karna, "Hyperpolarizabilities of organic molecules: Ab initio time-dependent coupled perturbed Hartree--fock--roothaan studies of basic heterocyclic structures," Journal of Physical Chemistry , vol. 99, no. 22, pp. 9045-9050, 1995.

[21] S. Millefiori, A. Alparone, "(Hyper)polarizability of chalcogenophenes C₄ H₄ X (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.

[22] S. Millefiori, A. Alparone, "Theoretical determination of the vibrational and electronic (hyper)polarizabilities of C₄ H₄ X (X = O, S, Se, Te) heterocycles," Physical Chemistry Chemical Physics , vol. 2, no. 11, pp. 2495-2501, 2000.

[23] K. Kamada, M. Ueda, H. Nagao, K. Tawa, T. Sugino, Y. Shmizu, K. Ohta, "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.

[24] K. Jug, S. Chiodo, P. Calaminici, A. Avramopoulos, M. G. Papadopoulos, "Electronic and vibrational polarizabilities and hyperpolarizabilities of azoles: a comparative study of the structure-polarization relationship," Journal of Physical Chemistry A , vol. 107, no. 20, pp. 4172-4183, 2003.

[25] B. Jansik, B. Schimmelpfennig, P. Norman, P. Macak, H. Ågren, 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.

[26] A. D. Becke, "A new mixing of Hartree-Fock and local density-functional theories," The Journal of Chemical Physics , vol. 98, no. 2, pp. 1372-1377, 1993.

[27] C. Lee, W. Yang, R. G. Parr, "Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density," Physical Review B , vol. 37, no. 2, pp. 785-789, 1988.

[28] J. P. Perdew, K. Burke, M. Ernzerhof, "Generalized gradient approximation made simple," Physical Review Letters , vol. 77, no. 18, pp. 3865-3868, 1996.

[29] F. A. Hamprecht, A. J. Cohen, D. J. Tozer, N. C. Handy, "Development and assessment of new exchange-correlation functionals," Journal of Chemical Physics , vol. 109, no. 15, pp. 6264-6271, 1998.

[30] J.-D. Chai, M. Head-Gordon, "Long-range corrected hybrid density functionals with damped atom-atom dispersion corrections," Physical Chemistry Chemical Physics , vol. 10, no. 44, pp. 6615-6620, 2008.

[31] 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.

[32] A. Alparone, "Evolution of electric dipole (hyper)polarizabilities of β -strand polyglycine single chains: an ab initio and DFT theoretical study," Journal of Physical Chemistry A , vol. 117, pp. 5184-5194, 2013.

[33] M. Alipour, A. Mohajeri, "Assessing the performance of density functional theory for the dynamic polarizabilities of amino acids: treatment of correlation and role of exact exchange," International Journal of Quantum Chemistry , vol. 113, pp. 1803-1811, 2013.

[34] A. Alparone, "Theoretical study on the static and dynamic first-order hyperpolarizabilities of adenine tautomers," Molecular Physics , 2013.

[35] A. J. Sadlej, "Medium-size polarized basis sets for high-level correlated calculations of molecular electric properties," Collection of Czechoslovak Chemical Communications , vol. 53, pp. 1995-2016, 1988.

[36] U. Eckart, M. P. Fülscher, L. Serrano-Andrés, A. J. Sadlej, "Static electric properties of conjugated cyclic ketones and thioketones," Journal of Chemical Physics , vol. 113, no. 15, pp. 6235-6244, 2000.

[37] T. Pluta, A. J. Sadlej, "Electric properties of urea and thiourea," Journal of Chemical Physics , vol. 114, no. 1, pp. 136-146, 2001.

[38] U. Eckart, V. E. Ingamells, M. G. Papadopoulos, A. J. Sadlej, "Vibrational effects on electric properties of cyclopropenone and cyclopropenethione," Journal of Chemical Physics , vol. 114, no. 2, pp. 735-745, 2001.

[39] A. Alparone, S. Millefiori, "Gas and solution phase electronic and vibrational (hyper)polarizabilities in the series formaldehyde, formamide and urea: CCSD(T) and DFT theoretical study," Chemical Physics Letters , vol. 416, no. 4-6, pp. 282-288, 2005.

[40] D. E. Woon, T. H. Dunning Jr., "Gaussian basis sets for use in correlated molecular calculations. IV. Calculation of static electrical response properties," The Journal of Chemical Physics , vol. 100, no. 4, pp. 2975-2988, 1994.

[41] H. Sekino, R. J. Bartlett, "Frequency dependent nonlinear optical properties of molecules," The Journal of Chemical Physics , vol. 85, no. 2, pp. 976-989, 1986.

[42] S. P. Karna, M. Dupuis, "Frequency dependent nonlinear optical properties of molecules: formulation and implementation in the HONDO program," Journal of Computational Chemistry , vol. 12, pp. 487-504, 1991.

[43] H. A. Kurtz, J. J. P. Stewart, K. M. Dieter, "Calculation of the nonlinear optical properties of molecules," Journal of Computational Chemistry , vol. 11, pp. 82-87, 1990.

[44] D. Birnbaum, B. E. Kohler, "Lowest energy excited singlet state of 2,2[variant prime]:5[variant prime],2[variant prime]- terthiophene, an oligomer of polythiophene," The Journal of Chemical Physics , vol. 90, no. 7, pp. 3506-3510, 1989.

[45] S. Siegert, F. Vogeler, C. M. Marian, R. Weinkauf, "Throwing light on dark states of α -oligothiophenes of chain lengths 2 to 6: radical anion photoelectron spectroscopy and excited-state theory," Physical Chemistry Chemical Physics , vol. 13, no. 21, pp. 10350-10363, 2011.

[46] M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci, G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P. Hratchian, A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, J. J. A. Montgomery, J. E. Peralta, F. Ogliaro, M. Bearpark, J. J. Heyd, E. Brothers, K. N. Kudin, V. N. Staroverov, R. Kobayashi, J. Normand, K. Raghavachari, A. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, R. L. Martin, K. Morokuma, V. G. Zakrzewski, G. A. Voth, P. Salvador, J. J. Dannenberg, S. Dapprich, A. D. Daniels, O. Farkas, J. B. Foresman, J. V. Ortiz, J. Cioslowski, D. J. Fox GAUSSIAN 09, Revision A. 02 , Gaussian, Wallingford, Conn, USA, 2009.

[47] B. Bak, D. Christensen, L. Hansen-Nygaard, J. Rastrup-Andersen, "The structure of thiophene," Journal of Molecular Spectroscopy , vol. 7, no. 1-6, pp. 58-63, 1961.

[48] J. E. Rice, R. D. Amos, S. M. Colwell, N. C. Handy, J. Sanz, "Frequency dependent hyperpolarizabilities with application to formaldehyde and methyl fluoride," The Journal of Chemical Physics , vol. 93, no. 12, pp. 8828-8839, 1990.

[49] H. Sekino, R. J. Bartlett, "Molecular hyperpolarizabilities," Journal of Chemical Physics , vol. 98, no. 4, pp. 3022-3037, 1993.

[50] F. Sim, S. Chin, M. Dupuis, J. E. Rice, "Electron correlation effects in hyperpolarizabilities of p-nitroaniline," Journal of Physical Chemistry , vol. 97, no. 6, pp. 1158-1163, 1993.

[51] G. Maroulis, "Hyperpolarizability of H₂ O revisited: accurate estimate of the basis set limit and the size of electron correlation effects," Chemical Physics Letters , vol. 289, no. 3-4, pp. 403-411, 1998.

[52] D. Xenides, G. Maroulis, "Basis set and electron correlation effects on the first and second static hyperpolarizability of SO₂ ," Chemical Physics Letters , vol. 319, no. 5-6, pp. 618-624, 2000.

[53] A. Alparone, H. Reis, M. G. Papadopoulos, "Theoretical investigation of the (hyper)polarizabilities of pyrrole Homologues C₄ H₄ XH (X = N, P, As, Sb, Bi). A coupled-cluster and density functional theory study," Journal of Physical Chemistry A , vol. 110, no. 17, pp. 5909-5918, 2006.

[54] A. Alparone, "Dipole (hyper)polarizabilities of fluorinated benzenes: an ab initio investigation," Journal of Fluorine Chemistry , vol. 144, pp. 94-101, 2012.

[55] T. Ogata, K. Kozima, "Microwave spectrum, barrier height to internal rotation of methyl group of 3-methylthiophene, and dipole moments of 3-methylthiophene and thiophene," Journal of Molecular Spectroscopy , vol. 42, no. 1, pp. 38-46, 1972.

[56] C. G. Le Fèvre, R. J. W. Le Fèvre, B. Purnachandra Rao, M. R. Smith, "Molecular polarisability. Ellipsoids of polarisability for certain fundamental heterocycles," Journal of the Chemical Society , pp. 1188-1192, 1959.

[57] M. H. Coonan, I. E. Craven, M. R. Hesling, G. L. D. Ritchie, M. A. Spackman, "Anisotropic molecular polarizabilities, dipole moments, and quadrupole moments of (CH₂ )₂ X, (CH₃ )₂ X, and C₄ H₄ X (X = O, S, Se). Comparison of experimental results and ab initio calculations," Journal of Physical Chemistry , vol. 96, no. 18, pp. 7301-7307, 1992.

[58] G. R. Dennis, I. R. Gentle, G. L. D. Ritchie, 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.

[59] M.-T. Zhao, B. P. Singh, P. N. Prasad, "A systematic study of polarizability and microscopic third-order optical nonlinearity in thiophene oligomers," The Journal of Chemical Physics , vol. 89, no. 9, pp. 5535-5541, 1988.

[60] B. Champagne, E. A. Perpète, S. J. A. Van Gisbergen, E.-J. Baerends, J. G. Snijders, C. Soubra-Ghaoui, K. A. Robins, B. Kirtman, "Assessment of conventional density functional schemes for computing the polarizabilities and hyperpolarizabilities of conjugated oligomers: an ab initio investigation of polyacetylene chains," Journal of Chemical Physics , vol. 109, no. 23, pp. 10489-10498, 1998.