The static and dynamic electronic (hyper)polarizabilities of the equilibrium conformations of 2,2′-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–C2–C2′–S dihedral angle of 137° (MP2/6-311G**) is the global minimum on the potential energy surface, whereas the syn-gauche rotamer (S–C2–C2′–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′-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.
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′-bithiophene (Figure 1), is the principal building block of polythiophene chains, extensively characterized by experimental and theoretical studies. In the solid state, 2,2′-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–C2–C2’–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′-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′-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′-bithiophene, the anti-gauche being predicted to be the global minimum [11–15]. The torsional potentials for the rotation around the C2–C2’ 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′-bithiophene rotamers was analyzed by Lukeš et al. [14] using HF/aug-cc-pVDZ and MP2/aug-cc-pVDZ computations. However, the average electronic polarizabilities of 2,2′-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].
Structures of thiophene and 2,2′-bithiophene rotamers. Colours: white (hydrogen), grey (carbon), and orange (sulphur).
In the present investigation, we have calculated the static and dynamic electronic (hyper)polarizabilities of the most stable conformations of 2,2′-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′-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′-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′-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′-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′-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 (〈α〉), polarizability anisotropy (Δα), first-order hyperpolarizability aligned along the μ direction (βμ), and average second-order hyperpolarizability (〈γ〉), which are orientationally invariant quantities:
(1)μ=μx2+μy2+μz2,〈α〉=13(αxx+αyy+αzz),Δα={(αxy2+αxz2+αyz212)12[(αxy2+αxz2+αyz212)(αxx-αyy)2+(αxx-αzz)2+(αyy-αzz)2hhh+6(αxy2+αxz2+αyz212)]}1/2,βμ=∑i=x,y,zβiμi|μ|,
where βi(i=x,y,z) is given by βi=(1/3)∑j=x,y,z(βijj+βjij+βjji)(2)〈γ〉=15(γxxxx+γyyyy+γzzzz+2γxxyyhh+2γxxzz+2γyyzz).
For the (hyper)polarizabilities, atomic units are used throughout the work. Conversion factors to the SI are 1 a.u. of α (e2a02Eh-1) = 1.648778 × 10−41 C2m2J−2; 1 a.u. of β (e3a03Eh-2) = 3.206361 × 10−53 C3m3J−2; 1 a.u. of γ (e4a04Eh-3) = 6.235377 × 10−65 C4m4J−3. All calculations were performed using the GAUSSIAN 09 program [46].
3. Results and Discussion3.1. Geometries and Energetics of Thiophene and 2,2′-Bithiophene Rotamers
The structural parameters of the equilibrium conformations of 2,2′-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′-bithiophene in vacuum (ca. 70%), the syn-gauche being, however, an important conformation (ca. 30%).
Gas phase MP2/6-311G** geometrical parameters of thiophene and 2,2′-bithiophene anti-gauche and syn-gauche rotamers. The experimental data reported in parentheses are taken from [47] (thiophene) and [9] (2,2′-bithiophene). The MP2/6-311G** (experimental) S–C2–C2’–S dihedral angles are 137° (148°, see [9]) and 48° for the anti-gauche and syn-gauche conformation, respectively.
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′-bithiophene is slightly less satisfactorily. In particular, the present computations underestimate the experimental typical S–C2–C2’–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–C2–C2’–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 (C4H4X, X = O, S, Se, Te) [23] and pyrrole (C4H4YH, 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 μ, 〈α〉, Δα, 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′-bithiophenes.
Dipole moments (D), static electronic polarizabilities (a.u.), and first-order hyperpolarizabilities (a.u.) of thiophenea.
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
αxx
68.61 (68.69)
71.49
7019
71.04
72.18
71.86
70.63
68.32
αyy
44.37 (44.51)
45.17
43.81
43.77
44.48
44.35
44.25
44.61
αzz
75.81 (76.00)
76.84
76.43
76.95
78.02
77.68
76.79
75.50
〈α〉
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
βxxz
−5.2 (−1.7)
−1.7
−1.2
1.7
2.2
2.9
−0.8
−5.1
βyyz
6.2 (7.3)
9.4
9.9
12.8
12.6
12.3
9.7
10.9
βzzz
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
aThe calculations were performed with the POL basis set on the MP2/6-311G** geometries.
bThe values in parentheses refer to the HF/aug-cc-pVTZ computations.
cValues taken from [23]. Basis set: 6-31G (3d) + pd.
dReference [55].
eReference [56].
fReference [57].
gReference [58].
hReference [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 〈α〉 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 〈α〉 (Δα) 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′-Bithiophene Rotamers
Table 2 lists the dipole moments and static (hyper)polarizabilities of the anti-gauche and syn-gauche 2,2′-bithiophene forms calculated at the HF/POL, MP2/POL, and BH&HLYP/POL levels. By analogy to the monomer, both the 2,2′-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 〈α〉 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 〈α〉 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–C2–C2’–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 〈α〉 and Δα estimates (within 2.5–4.4% and 0.4–4.6%, resp.). Similarly to the electronic polarizabilities, the 〈γ〉 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).
Dipole moments (D), static electronic polarizabilities (a.u.), first-order hyperpolarizabilities (a.u.), and second-order hyperpolarizabilities (a.u.) of 2,2′-bithiophene rotamersa.
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
αxx
178.53
184.76
186.63
177.29
183.13
184.95
αyy
84.85
86.72
84.30
125.80
127.96
126.58
αzz
127.99
129.81
128.68
86.73
88.65
86.27
〈α〉
130.46
133.76
133.20
129.94
133.25
132.60
Δα
81.71
85.28
89.22
78.70
82.29
85.98
βxxy
2.7
-3.1
7.0
0.0
0.0
0.0
βyyy
29.9
30.3
32.2
0.0
0.0
0.0
βzzy
16.1
-3.3
7.0
0.0
0.0
0.0
βxxz
0.0
0.0
0.0
45.7
31.5
59.0
βyyz
0.0
0.0
0.0
22.5
-4.0
27.0
βzzz
0.0
0.0
0.0
40.8
16.3
21.9
βμ
48.7
23.9
46.2
109.0
43.8
107.9
γxxxx
102128
90352
γyyyy
21824
22464
γzzzz
21392
21216
γxxyy
13312
12960
γxxzz
6432
7744
γyyzz
7536
6928
〈γ〉b
39981
37859
aThe calculations were performed with the POL basis set on the MP2/6-311G** geometries.
bThe HF/POL 〈γ〉 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 〈α〉 (Δα) values of both the conformations, respectively, by ca. 2 and 7% (4 and 18%). It is of interest noting that the 〈α〉 and 〈γ〉 data of the 2,2′-bithiophenes are greater than twice the corresponding values for the monomer, suggesting a some degree of interring π-conjugation in the dimers.
Static and frequency-dependent (ℏω=0.04282 a.u.) electronic polarizabilities (a.u.) and first-order hyperpolarizabilities (a.u.) of 2,2′-bithiophene rotamersa.
anti-gauche
syn-gauche
〈α〉(0;0)
130.46
129.94
〈α〉(-ω;ω)
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
aThe 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′-bithiophene, the static βμ value increases by a factor of 2.6, although the dimer exhibits a nonplanar arrangement. Note that the corresponding monomer→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′-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′-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.
RoncaliJ.Conjugated poly(thiophenes): synthesis, functionalization, and applications19929247117382-s2.0-4243893235ZyssJ.1994London, UKAcademic PressMcCulloughR. D.The chemistry of conducting polythiophenes1998102931162-s2.0-0031700973FichouD.1999Weinheim, GermanyWiley-VCHMishraA.MaC.-Q.BäuerleP.Functional oligothiophenes: molecular design for multidimensional nanoarchitectures and their applications20091093114111762-s2.0-6464910698610.1021/cr8004229HernandezV.CastiglioniC.Del ZoppoM.ZerbiG.Confinement potential and π-electron delocalization in polyconjugated organic materials19945014981598232-s2.0-000117492910.1103/PhysRevB.50.9815BrédasJ. L.StreetG. B.ThémansB.AndréJ. M.Organic polymers based on aromatic rings (polyparaphenylene, polypyrrole, polythiophene): evolution of the electronic properties as a function of the torsion angle between adjacent rings1985833132313292-s2.0-22844433265ChalonerP. A.GunatungaS. R.HitchcockP. B.Redetermination of 2, 2′-bithiophene19945019411942SamdalS.SamuelsenE. J.VoldenH. V.Molecular conformation of 2,2′-bithiophene determined by gas phase electron diffraction and ab initio calculations19935922592652-s2.0-0027624539ChadwickJ. E.KohlerB. E.Optical spectra of isolated s-cis- and s-trans-bithiophene: torsional potential in the ground and excited states19949814363136372-s2.0-0010773150KarpfenA.ChoiC. H.KerteszM.Single-bond torsional potentials in conjugated systems: a comparison of ab initio and density functional results199710140742674332-s2.0-0031549619MillefioriS.AlparoneA.MillefioriA.Conformational properties of thiophene oligomers20003748478532-s2.0-0033782638RaosG.FamulariA.MarconV.Computational reinvestigation of the bithiophene torsion potential20033793-43643722-s2.0-014156928710.1016/j.cplett.2003.08.060LukešV.BrezaM.BiskupičS.Structure and electronic properties of bithiophenes. I. Torsional dependence20026181-2931002-s2.0-003711182410.1016/S0166-1280(02)00429-3Ribeiro Da SilvaM. A. V.SantosA. F. L. O. M.GomesJ. R. B.RouxM. V.TempradoM.JiménezP.NotarioR.Thermochemistry of bithiophenes and thienyl radicals. A calorimetric and computational study20091134111042110502-s2.0-7794956712110.1021/jp905531cSalafskyJ. S.Second-harmonic generation as a probe of conformational change in molecules20033815-67057092-s2.0-024247536610.1016/j.cplett.2003.10.032VeldersG. J. M.GilletJ.-M.BeckerP. J.FeilD.Electron density analysis of nonlinear optical materials. An ab initio study of different conformations of benzene derivatives19919522860186082-s2.0-33751500018HendrickxE.ClaysK.PersoonsA.DehuC.BrédasJ. L.The bacteriorhodopsin chromophore retinal and derivatives: an experimental and theoretical investigation of the second-order optical properties199511712354735552-s2.0-0029141221LipińskiJ.BartkowiakW.Conformation and solvent dependence of the first and second molecular hyperpolarizabilities of charge-transfer chromophores. Quantum-chemical calculations19992451–32632762-s2.0-0039441117KeshariV.WijekoonW. M. K. P.PrasadP. N.KarnaS. P.Hyperpolarizabilities of organic molecules: Ab initio time-dependent coupled perturbed Hartree—fock—roothaan studies of basic heterocyclic structures19959922904590502-s2.0-33751155296MillefioriS.AlparoneA.(Hyper)polarizability of chalcogenophenes C4H4X (X = O, S, Se, Te) conventional ab initio and density functional theory study19984311-259782-s2.0-0002074771MillefioriS.AlparoneA.Theoretical determination of the vibrational and electronic (hyper)polarizabilities of C4H4X (X = O, S, Se, Te) heterocycles2000211249525012-s2.0-003421416710.1039/b001020kKamadaK.UedaM.NagaoH.TawaK.SuginoT.ShmizuY.OhtaK.Molecular design for organic nonlinear optics: polarizability and hyperpolarizabilities of furan homologues investigated by ab initio molecular orbital method200010420472347342-s2.0-003374138610.1021/jp993806yJugK.ChiodoS.CalaminiciP.AvramopoulosA.PapadopoulosM. G.Electronic and vibrational polarizabilities and hyperpolarizabilities of azoles: a comparative study of the structure-polarization relationship200310720417241832-s2.0-003760759510.1021/jp022403mJansikB.SchimmelpfennigB.NormanP.MacakP.ÅgrenH.OhtaK.Relativistic effects on linear and non-linear polarizabilities of the furan homologues20036332-32372462-s2.0-004247061510.1016/S0166-1280(03)00278-1BeckeA. D.A new mixing of Hartree-Fock and local density-functional theories1993982137213772-s2.0-34250817103LeeC.YangW.ParrR. G.Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density19883727857892-s2.0-034549110510.1103/PhysRevB.37.785PerdewJ. P.BurkeK.ErnzerhofM.Generalized gradient approximation made simple19967718386538682-s2.0-4243943295HamprechtF. A.CohenA. J.TozerD. J.HandyN. C.Development and assessment of new exchange-correlation functionals199810915626462712-s2.0-000165276210.1063/1.477267ChaiJ.-D.Head-GordonM.Long-range corrected hybrid density functionals with damped atom-atom dispersion corrections20081044661566202-s2.0-5584911739910.1039/b810189bAlparoneA.Comparative study of CCSD(T) and DFT methods: electronic (hyper)polarizabilities of glycine20115141–321252-s2.0-8005292279610.1016/j.cplett.2011.08.010AlparoneA.Evolution of electric dipole (hyper)polarizabilities of β-strand polyglycine single chains: an ab initio and DFT theoretical study201311751845194AlipourM.MohajeriA.Assessing the performance of density functional theory for the dynamic polarizabilities of amino acids: treatment of correlation and role of exact exchange201311318031811AlparoneA.Theoretical study on the static and dynamic first-order hyperpolarizabilities of adenine tautomers201310.1080/00268976.2013.861088SadlejA. J.Medium-size polarized basis sets for high-level correlated calculations of molecular electric properties19885319952016EckartU.FülscherM. P.Serrano-AndrésL.SadlejA. J.Static electric properties of conjugated cyclic ketones and thioketones200011315623562442-s2.0-034273525510.1063/1.1290012PlutaT.SadlejA. J.Electric properties of urea and thiourea200111411361462-s2.0-003514826910.1063/1.1328398EckartU.IngamellsV. E.PapadopoulosM. G.SadlejA. J.Vibrational effects on electric properties of cyclopropenone and cyclopropenethione200111427357452-s2.0-003512814610.1063/1.1331358AlparoneA.MillefioriS.Gas and solution phase electronic and vibrational (hyper)polarizabilities in the series formaldehyde, formamide and urea: CCSD(T) and DFT theoretical study20054164–62822882-s2.0-2844445423510.1016/j.cplett.2005.09.088WoonD. E.DunningT. H.Jr.Gaussian basis sets for use in correlated molecular calculations. IV. Calculation of static electrical response properties19941004297529882-s2.0-26444474547SekinoH.BartlettR. J.Frequency dependent nonlinear optical properties of molecules19868529769892-s2.0-36549102995KarnaS. P.DupuisM.Frequency dependent nonlinear optical properties of molecules: formulation and implementation in the HONDO program199112487504KurtzH. A.StewartJ. J. P.DieterK. M.Calculation of the nonlinear optical properties of molecules1990118287BirnbaumD.KohlerB. E.Lowest energy excited singlet state of 2,2′:5′,2′- terthiophene, an oligomer of polythiophene1989907350635102-s2.0-0013596856SiegertS.VogelerF.MarianC. M.WeinkaufR.Throwing light on dark states of α-oligothiophenes of chain lengths 2 to 6: radical anion photoelectron spectroscopy and excited-state theory2011132110350103632-s2.0-7995746889010.1039/c0cp02712jFrischM. J.TrucksG. W.SchlegelH. B.ScuseriaG. E.RobbM. A.CheesemanJ. R.ScalmaniG.BaroneV.MennucciB.PeterssonG. A.NakatsujiH.CaricatoM.LiX.HratchianH. P.IzmaylovA. F.BloinoJ.ZhengG.SonnenbergJ. L.HadaM.EharaM.ToyotaK.FukudaR.HasegawaJ.IshidaM.NakajimaT.HondaY.KitaoO.NakaiH.VrevenT.MontgomeryJ. J. A.PeraltaJ. E.OgliaroF.BearparkM.HeydJ. J.BrothersE.KudinK. N.StaroverovV. N.KobayashiR.NormandJ.RaghavachariK.RendellA.BurantJ. C.IyengarS. S.TomasiJ.CossiM.RegaN.MillamJ. M.KleneM.KnoxJ. E.CrossJ. B.BakkenV.AdamoC.JaramilloJ.GompertsR.StratmannR. E.YazyevO.AustinA. J.CammiR.PomelliC.OchterskiJ. W.MartinR. L.MorokumaK.ZakrzewskiV. G.VothG. A.SalvadorP.DannenbergJ. J.DapprichS.DanielsA. D.FarkasO.ForesmanJ. B.OrtizJ. V.CioslowskiJ.FoxD. J.2009Wallingford, Conn, USAGaussianBakB.ChristensenD.Hansen-NygaardL.Rastrup-AndersenJ.The structure of thiophene196171–658632-s2.0-0003072227RiceJ. E.AmosR. D.ColwellS. M.HandyN. C.SanzJ.Frequency dependent hyperpolarizabilities with application to formaldehyde and methyl fluoride19909312882888392-s2.0-36549092205SekinoH.BartlettR. J.Molecular hyperpolarizabilities1993984302230372-s2.0-36449009375SimF.ChinS.DupuisM.RiceJ. E.Electron correlation effects in hyperpolarizabilities of p-nitroaniline1993976115811632-s2.0-0001751863MaroulisG.Hyperpolarizability of H2O revisited: accurate estimate of the basis set limit and the size of electron correlation effects19982893-44034112-s2.0-0032510862XenidesD.MaroulisG.Basis set and electron correlation effects on the first and second static hyperpolarizability of SO220003195-66186242-s2.0-0001495775AlparoneA.ReisH.PapadopoulosM. G.Theoretical investigation of the (hyper)polarizabilities of pyrrole Homologues C4H4XH (X = N, P, As, Sb, Bi). A coupled-cluster and density functional theory study200611017590959182-s2.0-3374446908810.1021/jp0574605AlparoneA.Dipole (hyper)polarizabilities of fluorinated benzenes: an ab initio investigation201214494101OgataT.KozimaK.Microwave spectrum, barrier height to internal rotation of methyl group of 3-methylthiophene, and dipole moments of 3-methylthiophene and thiophene197242138462-s2.0-0002130737Le FèvreC. G.Le FèvreR. J. W.Purnachandra RaoB.SmithM. R.Molecular polarisability. Ellipsoids of polarisability for certain fundamental heterocycles1959118811922-s2.0-37049060808CoonanM. H.CravenI. E.HeslingM. R.RitchieG. L. D.SpackmanM. A.Anisotropic molecular polarizabilities, dipole moments, and quadrupole moments of (CH2)2X, (CH3)2X, and C4H4X (X = O, S, Se). Comparison of experimental results and ab initio calculations19929618730173072-s2.0-0000634566DennisG. R.GentleI. R.RitchieG. L. D.AndrieuC. G.Field-gradient-induced birefringence in dilute solutions of furan, thiophen and selenophen in cyclohexane19837945395452-s2.0-000081699810.1039/F29837900539ZhaoM.-T.SinghB. P.PrasadP. N.A systematic study of polarizability and microscopic third-order optical nonlinearity in thiophene oligomers1988899553555412-s2.0-0010256429ChampagneB.PerpèteE. A.Van GisbergenS. J. A.BaerendsE.-J.SnijdersJ. G.Soubra-GhaouiC.RobinsK. A.KirtmanB.Assessment of conventional density functional schemes for computing the polarizabilities and hyperpolarizabilities of conjugated oligomers: an ab initio investigation of polyacetylene chains19981092310489104982-s2.0-003852753010.1063/1.477731