Band gap modulation in polythiophene and polypyrrole-based systems.

In this paper, the structural and electronic properties of polythiophene and polyprrrole-based systems have been investigated using first-principles calculations both in periodic and oligomer forms. Of particular interest is the band gap modulation through substitutions and bilayer formation. Specifically, S has been substituted by Se and Te in polythiophene, leading to polyseleophene and polytellurophene, respectively, and N has been substituted by P and As in polypyrrole. The values obtained of the binding energy suggest that all the systems studied can be realized experimentally. Stacking (bilayer formation) of pure polythiophene, polypyrrole and their derivatives leads to linear suppression of the band gap or HOMO-LUMO gap as a function of the stacking. Mixed bilayers, including one formed from polythiophene on top of polypyrrole, have also been considered. Overall, a wide range of band gaps can be achieved through substitutions and stacking. Hybrid (B3LYP) calculations also suggest the same trend in the band gap as PBE calculations. Trends in the binding energy are similar for both periodic and molecular calculations. In addition, Γ-point phonon calculations were performed in order to check the stability of selected systems.

In 1976, π-n class="Chemical">conducting polymers were discovered by Heeger, MacDiarmid, Shirakawa and co-workers. This discovery opened the way for comprehensive investigations and understanding of various aspects of the π-conducting polymers from both physical and chemical points of views1. The π-conjugated polymers possess alternating single and double bonds of carbon atoms leading to one π electron for each of the carbon atoms23456. The carbon atoms are sp2p hybridized and overlap along the polymer chain, inducing electron delocalization in the polymer chain78. The electronic delocalization opens up a path for charge transport along the polymer chain, and as a consequence, the polymer chain can have conducting or semiconducting nature91011121314151617181920212223. Organic π-conjugated polymers have been investigated widely in experiment and theory. These polymers are of great interest because of their potential to be integrated into electronic and optical devices12425262728293031323334353637383940 as well as sensors4414243444546.

It hn class="Chemical">as already been reporpan> class="Chemical">ted that the electronic properties of such polymers can be engineered/tuned without much efforts either by chemical modification or atomic/molecular scale doping846474849. n class="Chemical">Polyacetylene can be doped at room temperature by various dopants over a wide range of concentrations such that the structural and electronic properties can be tuned150. Doping has been widely utilized for various conjugated polymers. It was found that dopants play a crucial role in changing their properties compared to those of conventional semiconductors46. The doping can be achieved either by charge transfer or by application of an external electric field. By applying these doping techniques, the electronic and optical properties of conducting polymers can be engineered widely, including the transition from semiconductor to metal or further to insulator depending upon the dopant concentration48515253. This induces an outstanding opportunity for switching between conducting and insulating properties, which indeed paves the way for applications in optoelectronics, such as organic polymer based transistors, photoresistances, light-emitting diodes, and polythiophene based organic solar cells54.

There have been n class="Chemical">several efforts to address the electronic structure of n class="Chemical">polythiophene derivatives using first-principles calculations535556575859606162 Substitution in polythiophenes, such as replacement of H by CH3, NH2, NO2, and Cl, has been studied. It was found that the electronic structure can easily be tuned62. The optical properties of fully conjugated cyclo[n]thiophenes have been studied in the experiment and theory and it has been demonstrated that the band gap and optical properties depend on structure of the thiophene molecules6364. In addition, polythiophene and its derivatives are expected to be utilized in electronic memory devices. These materials have been in focus for the last two decades5462. Injection of electrons from the highest occupied molecular orbital (HOMO) to the lowest unoccupied molecular orbital (LUMO) is achieved without difficulties in this class of materials65. Recently, synthesis of polyselenophene/polytellurophene (replacement/substitution of the S atoms by Se/Te) has been achieved6166676869. The material can be used as an active layer in field-effect transistors68. Moreover, in a recent experiment, it has been reported that such systems are promising for CO2 capture70.

Using a n class="Chemical">conjugated n class="Chemical">polymer, a well-established barrier to Li-ion drift has been achieved experimentally using polythiophene and polypyrrole, which in fact provides a mechanism for creating high-performance memory devices71. For such a system, the conjugated polymer barrier layer made of polypyrrole remains conducting at the interface between the metal oxide and Li+ ion doped polythiophene. Moreover, the structural and electronic properties of polythiophene doped by Li/Cl in various stacking schemes in both periodic and oligomer forms have been investigated53. It has been suggested that these types of study will be useful for understanding the structural properties and the tunability of the electronic states, which should be an important step toward construction of polythiophene-based electronic devices. In this paper, we investigate the structural and electronic properties of S substituted by Se and Te in polythiophene and N substituted by P and As in polypyrrole using first-principle density functional theory (DFT) based calculations. It is believed that the comprehensive analysis of our investigation will be very useful for understanding the structural properties and the band gap or HOMO-LUMO gap of the systems under study, which provides a basis for constructing polythiophene and polypyrrole-based electronic devices. Moreover, in case of the periodic calculations, we also performed hybrid (B3LYP) calculations in order to determine the effect of the hybrid functional on the band gap. Furthermore, Γ-point phonon calculations were performed to confirm the stability of selected systems.

Methods

Periodic calculations were performed using DFT within the generalized gradient approximation in the Perdew, Burke, and Ernzerhof (n class="Gene">PBE) parametrizationpan>72 pan> class="Chemical">as implemented in the Quantum-ESPRESSO code73. The effects of spin-orbit coupling (spin-orbit relativistic effects with PBE pseudopotential) were also taken into account. For the systems with stacking configurations, the van der Waals interactions (DFT-D) were included in order to obtain correct dispersion74. A relatively high plane wave cutoff energy of 800 eV with a larger Monkhorst-Pack 32 × 1 × 1 k-mesh were employed. The considered supercell contains four monomers of thiophene as well as pyrrole and has lattice constants of a = 15.59 Å and b = 15.00 Å. A vacuum layer of 20 Å was used in order to nullify any spurious interactions due to the periodic boundary conditions. The atomic positions were optimized very well with all the forces converged to 0.0005 eV/Å. The following systems were considered with a vacuum along the y- and z-directions (see Tables 1, 2, 3 and 4) (i) monolayer/bilayer of pure polythiophene, (ii) monolayer/bilayer of polythiophene substituted by Se (polyselenophene), (iii) monolayer/bilayer of polythiophene substituted by Te (polytellurophene), (iv) monolayer/bilayer of pure polypyrrole, (vi) monolayer/bilayer of polypyrrole substituted by P, (vi) monolayer/bilayer of polypyrrole substituted by As, and (vii) bilayer made of polythiophene and polypyrrole. In addition, various mixed bilayers were considered as summarized in Tables 1, 2, 3 and 4. All of the above systems were also considered for the molecular case. In addition, the B3LYP hybrid functional75 was employed in order to examine the effect of the hybrid functional on the band gap. Note that the B3LYP calculations were performed only for the periodic cases and not extended to the molecular cases due to limitations in computational resources. We believe that the trends in the band gap would be the same in molecular and periodic calculations.

Table 1

Structural parameters, periodic calculations: selected bond lengths (in Å) and bond angles (in °).

System	X–H	C–C	C–X	C-X-C	C-C-C	X-C-C
Polythiophene derivatives, monolayers
S	—	1.38–1.41	1.73	93	113–128	109–119
Se	—	1.37–1.40	1.88	85	114–129	112–117
Te	—	1.35–1.38	2.08	76	115–127	116–119
Bilayer (identical layers) polythiophene with S, Se, and Te
S-S	—	1.38–1.43	1.74	92	113–129	109–121
Se-Se	—	1.39–1.44	1.88	85	113–128	112–117
Te-Te	—	1.37–1.45	2.07	76	116–127	116–119
Bilayer (non-identical layers) polythiophene with S, Se, and Te
S-Se	—	1.37–1.39	1.73–1.88	85–92	114–129	110–121
S-Te	—	1.36–1.38	1.74–2.06	76–92	113–128	111–120
Se-Te	—	1.36–1.38	1.88–2.07	76–85	114–128	112–116
Polythiophene with alternative S replaced by Se and Te
S-Se-Te1	—	1.38–1.40	1.72–2.08	79–91	117–129	112–120
S-Se-Te2	—	1.39–1.41	1.71–2.07	79–92	115–128	111–119
Monolayer polypyrrole with N, P, and As
N	1.01	1.40–1.55	1.44	113	110–128	102–128
P	1.43	1.37–1.42	1.82	88	114–126	110–122
As	1.55	1.36–1.42	1.99	85	117–126	110–123
Bilayer (identical layers) polypyrrole with N, P, and As
N-N	1.01	1.42–1.55	1.45	112	110–128	102–129
P-P	1.54	1.37–1.42	1.99	84	117–126	109–129
As-As	1.54	1.38–1.42	1.99	84	117–125	109–123
Bilayer (non-identical layers) polypyrrole with N, P, and As
N-P	1.01–1.43	1.40–1.55	1.45–1.82	89–113	110–128	102–128
N-As	1.02–1.55	1.36–1.42	1.45–1.48	84–113	110–128	102–128
P-As	1.43–1.54	1.37–1.42	1.83–1.99	84–89	114–126	109–123
Bilayer composed of pure polythiophene and polypyrrole
S-P	1.02	1.39–1.42	1.44–1.74	93–113	110–129	110–128

Table 2

Periodic calculations: the calculated values of the valence and conduction bandwidths (in eV), band gap (in eV), band gap in valence and conduction bands at the Γ point (in eV), binding energy (in eV), interlayer spacing (in Å), and buckling (in Å), where X = S, Se, Te, N, P, and As and S-Se-Te1 and S-Se-Te2 correspond to single layer and bilayers.

System	VBW	CBW	Egap (PBE)	Egap (B3LYP)			Ebinding	dint	Buckling
Polythiophene derivatives, monolayers
S	0.74	0.63	1.12	2.08	—	—	−0.483	—	0.10
Se	0.76	0.69	0.68	1.53	—	—	−0.462	—	0.12
Te	0.74	0.81	0.26	0.97	—	—	−0.455	—	0.17
Bilayer (identical layers) polythiophene with S, Se, and Te
S-S	1.14	0.62	0.83	1.81	0.14	0.21	−0484	3.33–3.43	0.11
Se-Se	0.75	0.70	0.51	1.33	0.16	0.24	−0.463	3.36–3.44	0.12
Te-Te	0.79	0.73	0.22	0.81	0.31	0.40	−0.457	3.36–3.42	0.16
Bilayer (non-identical layers) polythiophene with S, Se, and Te
S-Se	0.77	0.63	0.64	1.51	0.31	0.29	−0.476	3.42–3.62	0.46
S-Te	0.75	0.73	0.44	1.29	0.49	0.44	−0.469	3.39–3.65	0.72
Se-Te	0.78	0.72	0.36	1.19	0.32	0.40	−0.462	3.39–3.68	0.84
Polythiophene with alternative S replaced by Se and Te
S-Se-Te1	0.61	0.63	0.70	1.54	—	—	−0.498	—	0.12
S-Se-Te2	0.77	0.65	0.64	1.48	0.17	0.14	−0.496	3.48–3.75	0.74
Monolayer polypyrrole with N, P, and As
N	0.61	0.47	1.88	3.01	—	—	−0.492	—	0.11
P	0.88	0.79	0.34	1.27	—	—	−0.461	—	0.17
As	0.67	0.72	0.61	1.58	—	—	−0.448	—	1.91
Bilayer (identical layers) polypyrrole with N, P, and As
N-N	0.71	0.61	0.86	2.58	0.13	0.26	−0.494	3.18–3.68	1.52
P-P	0.74	0.69	0.51	1.49	0.15	0.22	−0.464	1.73–2.69	1.87
As-As	0.77	0.73	0.22	1.08	0.30	0.40	−0.450	1.73–2.72	1.89
Bilayer (non-identical layers) polypyrrole with N, P, and As
N-P	0.27	0.70	0.46	1.81	0.35	—	−0.491	2.31–2.67	0.81
N-As	0.81	0.70	0.29	1.17	0.24	0.12	−0.462	1.51–3.98	0.99–2.24
P-As	0.57	0.64	0.54	1.88	0.08	—	−0.453	1.74–2.01	0.51–1.64
Bilayers composed of pure polythiophene and polypyrrole
S-P	0.60	0.67	0.72	1.72	0.41	1.03	−0.491	2.29–3.76	1.69

Table 3

Structural parameters, molecular calculations: selected bond lengths (in Å) and bond angles (in °).

System	X–H	C–C	C–X	C-X-C	C-C-C	X-C-C
Polythiophene derivatives, monolayers
S	—	1.40–1.44	1.75	92	113–129	110–120
Se	—	1.40–1.42	1.91	87	115–129	111–120
Te	—	1.38–1.41	2.08	82	118–128	108–122
Bilayer (identical layers) polythiophene with S, Se, and Te
S-S	—	1.39–1.44	1.75	91	114–128	111–120
Se-Se	—	1.38–1.43	1.90	87	115–129	109–120
Te-Te	—	1.38–1.42	2.12	82	118–128	109–122
Bilayer (non-identical layers) polythiophene with S, Se, and Te
S-Se	—	1.38–1.42	1.73–1.91	87–91	113–129	110–121
S-Te	—	1.38–1.43	1.75–2.12	81–92	113–128	111–121
Se-Te	—	1.38–1.43	1.91–2.12	81–87	116–129	109–122
Polythiophene with alternative S replaced by Se and Te
S-Se-Te1	—	1.41–1.43	1.73–1.12	82–83	114–128	109–122
S-Se-Te2	—	1.40–1.42	1.74–2.11	81–84	113–127	111–123
Monolayer polypyrrole with N, P, and As
N	1.01	1.42–1.44	1.38	111	108–132	107–122
P	1.44	1.36–1.43	1.85	90	114–127	110–124
As	1.54	1.37–1.43	1.99	86	117–127	109–124
Bilayer (identical layers) polypyrrole with N, P, and As
N-N	1.01	1.39–1.44	1.38	110	108–131	106–122
P-P	1.43	1.38–1.42	1.85	90	115–127	108–124
As-As	1.54	1.39–1.41	1.99	86	117–127	108–124
Bilayer (non-identical layers) polypyrrole with N, P, and As
N-P	1.02–1.73	1.39–1.44	1.38–1.85	91–110	115–127	106–124
N-As	1.01–1.54	1.39–1.44	1.96–1.38	85–110	108–131	107–1.21
P-As	1.43–1.54	1.39–1.42	1.85–1.99	86–90	114–126	108–123
Bilayer composed of pure polythiophene and polypyrrole
S-P	1.01	1.37–1.43	1.38–1.73	92–110	107–131	107–121

Table 4

Molecular calculations: The calculated values of the highest occupied (HOMO) and lowest unoccupied (LUMO) molecular orbitals (in eV), HOMO-LUMO (in eV), binding energy (in eV), interlayer spacing (in Å), buckling (in Å), bond lengths (in Å), and bond angles (in °).

System	HOMO	LUMO	HOMO-LUMO	Ebinding	dint	Buckling
Polythiophene derivatives, monolayers
S	−4.547	−3.034	1.51	−0.596	—	0.00
Se	−4.535	−3.226	1.31	−0.582	—	0.00
Te	−4.516	−3.377	1.14	−0.569	—	0.00
Bilaye Bilayer (identical layers) polythiophene with S, Se, and Te
S-S	−4.273	−3.112	1.16	−0.601	3.41–3.53	0.25
Se-Se	−4.335	−3.261	1.07	−0.587	3.44–3.61	0.35
Te-Te	−4.277	−3.381	0.89	−0.575	3.39–3.52	0.46
Bilayer Bilayer (non-identical layers) polythiophene with S, Se, and Te
S-Se	−4.308	−3.118	1.19	−0.594	3.33–3.49	0.15–0.89
S-Te	−4.380	−3.279	1.10	−0.588	3.31–3.50	0.40–0.79
Se-Te	−4.345	−3.303	1.04	−0.577	3.41–3.58	0.42–0.65
Polythiophene with alternative S is replaced by Se and Te
S-Se-Te1	−4.533	−3.203	1.23	−0.582	—	0.00
S-Se-Te2	−4.511	−3.136	1.37	−0.579	3.42–3.51	0.23–0.27
Monolayer polypyrrole with N, P, and As
N	−3.936	−1.687	2.25	−0.609	—	0.09–0.12
P	−4.464	−3.570	0.89	−0.562	—	0.25–1.12
As	−4.429	−3.617	0.81	−0.548	—	0.27–1.14
Bilayer (identical layers) polypyrrole with N, P, and As
N-N	−3.945	−1.745	2.20	−0.610	0.96–5.53	1.72–3.93
P-P	−4.454	−3.591	0.86	−0.599	0.75–3.57	0.60–0.71
As-As	−4.366	−3.605	0.76	−0.588	3.02–3.44	0.41–0.48
Bilayer (non-identical layers) polypyrrole with N, P, and As
N-P	−3.995	−3.580	0.41	−0.586	2.17–2.73	0.87–3.38
N-As	−3.992	−3.592	0.40	−0.582	1.45–3.36	0.38–2.32
P-As	−4.410	−3.682	0.73	−0.578	2.43–2.89	0.66–3.21
Bilayer composed of pure polythiophene and polypyrrole
S-P	−3.931	−3.057	0.87	−0.589	3.55–3.75	0.76–0.98

In the molecular calculations (hexamers), the sysn class="Chemical">tems were optimized usinpan>g DFT-D (n class="Gene">PBE)72 as implemented in the Amsterdam Density Functional (ADF) code767778. Previously, we have shown that this setup provides reasonable results while comparing to the corresponding experimental findings53. Note that experimentally and theoretically, oligomers of six or less thiophene rings are found to be sufficient to reasonably model such systems79. All the calculations were performed using double zeta Slater-type orbital basis sets (ZORA-DZ). Relativistic effects including spin-orbit effects were included (ZORA with spin-orbit)808182, they cannot be neglected due to larger values of the spin-orbit coupling and scalar relativistic effects especially for Se and Te atoms. All the geometries were optimized until the energy was converged to 0.0003 eV. A default integration accuracy parameter of 4.0 was used for geometry optimizations and single point calculations for all systems under study.

Results and Discussion

Polythiophene and its derivatives (Periodic approach)

Vibrational properties

Vibrational Raman and infrared frequencies are really un class="Chemical">seful to provide inpan>formationpan> onpan> a molecular or atomic level. The vibrationpan>al frequencies can easily be detected in experiment and the frequencies can vary between different materials as well as due to distortions in the structures arising from doping or defects838485. Therefore, first of all, we performed a series of phonon calculations at the Γ-point phonon for polythiophene and its derivatives in order to check the structural stability. Indeed, only stable structures are presented in this paper as listed in Table 1. Each atom has three degrees of vibrational freedom, each monolayer of polythiophene and its derivatives have a total of 28 atoms in the system, thus, the total number of phonon modes in these structures is 84, see Fig. 1(a), whereas the number of phonon modes becomes 168 for identical and non-identical bilayers (28 atoms in each layer), see Fig. 1(b,c). The highest phonon frequency for a monolayer of polythiophene is found to be 3195 cm−1 at the Γ-point, which agrees well with an experimental observation83. Note that no imaginary phonon frequencies are observed at the Γ-point for all the system under study, which confirms that the systems are stable. The highest phonon frequency for monolayers of polyselenophene and polytellurophene is 3124 cm−1 and 3077 cm−1 at the Γ-point, respectively. Our obtained values show good agreement with a previous theoretical prediction86. It is clear that the phonon frequencies are decreasing while moving from polythiophene to polytellurophene (S to Te), which is called phonon softening87. The phonon softening occurs due to the the weakening of C-Se and C-Te bond lengths as compared to C-S bond lengths. However, the phonon frequencies become harder in case of identical as well as non-identical bilayers as compared to their monolayer counterparts because the interlayer interactions (out-of-plane vibrations) play a pivotal role. The obtained value of the phonon frequencies for the bilayers are found to be 3207 cm−1, 3178 cm−1, 3198 cm−1, 3210 cm−1, 3260 cm−1, and 3310 cm−1 for S-S, Se-Se, Te-Te, S-Se, Se-Te, and S-Te (see Fig. 1(b,c), respectively.

Figure 1

Γ-point phonon modes of polythiophene derivatives (a) monolayers (S, Se, and Te), (b) identical bilayers (S-S, Se-Se, Te-Te), and (c) non-identical bilayers (S-Se, S-Te, and Se-Te), see the optimized structures depicted in Fig. 3(a) and left panel of Fig. 4.

Furthermore, similar to the n class="Chemical">polythiophene derivatives, we also performed a n class="Chemical">series of phonon calculation at the Γ-point for n class="Chemical">polypyrrole and its derivatives to check the structural stability. Only stable structures are presented in Table 1, Fig. 3(b). Each monolayer of polypyrrole and its derivatives has a total of 32 atoms in a unit cell; therefore, the total number of phonon modes in these structures is 96, see Fig. 2(a). But the number of phonon modes becomes 192 for identical and non-identical bilayers (32 atoms in each layer) as depicted in the Fig. 2(b,c). We observe a highest phonon frequency of 3559 cm−1 at the Γ-point for the monolayer of polypyrrole, which agrees well with a previously reported value88. No imaginary phonon frequencies are observed at the Γ-point for all the system under consideration, which confirms that the structures presented in the paper are stable. The highest phonon frequency for P and As amounts to 3103 cm−1 and 3074 cm−1 at the Γ-point. It should be noted that we do not have any previously published data to be compared with our results for these systems. Similarly to the polythiophene, it is clear that the phonon frequencies are decreasing while moving from N to P to As; the weaker bonding is responsible for the softening of the phonon as we move from N to P to As, see Fig. 2(a). However, the phonon frequency becomes softer/harder in case of identical/non-identical bilayers as compared to their monolayer counterparts because the interlayer interactions (out-of-plane vibrations) play a crucial role. In contrast to polythiophene bilayers, the phonon frequencies in case of identical bilayers become softer due to the fact of a strong deformation in the atomic structure in both of the layers, while non-identical bilayers follow the same trend as the polythiophene derivatives (non-identical bilayers). The calculated values of the highest phonon frequencies for bilayers are found to be 3094 cm−1, 2828 cm−1, 2640 cm−1, 3533 cm−1, 3195 cm−1, and 3541 cm−1 for N-N, P-P, As-As, N-P, P-As, and As-N (see Fig. 2(b,c)), respectively.

Figure 2

Γ-point phonon modes of polypyrrole derivatives (a) monolayers (N, P, and As), (b) identical bilayers (N-N, P-P, and As-As), and (c) non-identical bilayers (N-P, N-As, and P-As, see the optimized structures presented in Fig. 3(b) and right panel of Fig. 4.

Figure 3

Top view of the structures (periodic approach) of (a) polythiophene and (b) polypyrrole based materials, where yellow, blue, black, and gray balls represent C, H, S or Se or Te, and N or P or As atoms, respectively.

Structural arrangements

The structure of monolayer pan> class="Chemical">polythiophene is shown in Fig. 3(a). The polyselenophene and n class="Chemical">polytellurophene have very similar structures. Side views of the optimized structures of various bilayers are presented in Fig. 4. The considered supercell of the monolayer has 16 carbon (C) atoms, 8 Hydrogen (H) atoms, and 4 Sulphur (S) or Selenium (Se) or Tellurium (Te) atoms. The structural parameters, which consist of various bond lengths and bond angles, are summarized in Table 1. The values obtained of bond lengths and bond angles agree well with the available experimentally and theoretically obtained values536689. Polythiophene, containing S atoms in the pentagonal rings of the polymer backbone, has been well studied in both experiment and theory, whereas the effect of Se or Te instead of S is not well understood. Recently, replacement/substitution of S by Se has been successfully realized experimentally, resulting in polyselenophene536190. It has been found that such a substitution is crucial in order to control the electronic properties of the polythiophene based polymer, especially the band gap.

Figure 4

Optimized bilayer structures, periodic calculations.

Yellow, blue, black, green, and gray balls represent C, H, S/As, Se/P, and Te/N atoms, respectively.

For pristine n class="Chemical">polythiophene, the buckling, which is defined pan> class="Chemical">as the perpendicular out-of-plane distance between atoms in the same layer, is found to be 0.10 Å (Table 2), C–C bond lengths are found to be 1.38–1.41 Å, and the C–S bond length amounts to 1.73 Å. While the n class="Chemical">C-S-C, C-C-C, and S-C-C bond angles are found to be 93°, 113–128°, and 109–119°, respectively. The calculated structural parameters for polyselenophene are found to be similar to those of polythiophene except for the C–Se bond length, which grows as compared to the C–S bond length in polythiophene and is found to be 1.88 Å. At the same time, the C–C bond length shrinks slightly, in good agreement with the experimentally obtained value61.

Our calculated values for the structural parameters of polytellurophene are listed in Table 1 as well. The buckling is found to be 0.17 Å, slightly increased as compared to polythiophene (0.10 Å) and polyselenophene (0.12 Å). The C–C bond lengths shrink to 1.35–1.38 Å as compared to the C–C bond lengths in polythiophene and polyselenophene, while the C–Te bond length grows to be 2.08 Å. All other parameters are similar to those of the previous cases. It should be noted that C–C bond lengths decrease from S to Te and the C–X (X = S, Se, Te) bond length grows from Se to Te. Note that the exact structural modifications at finite temperature can be seen from molecular dynamics studies, which indeed is close to the experimental observations, as recently shown by Salvador and co-workers91.

In addition, we have studied bilayers of n class="Chemical">polythiophene, pan> class="Chemical">polyselenophene, and polytellurophene in two scheme (i) both layers are identical (S-S, Se-Se, and Te-Te) and (ii) the two layers are not identical (S-Se, S-Te, and Se-Te). The side views of the optimized structures are depicted in Fig. 4 (left side). The interlayer spacings for scheme (i) are found to be 3.33–3.43 Å, 3.36–3.44 Å, and 3.36–3.42 Å, for polythiophene, polyselenophene, and polytellurophene, respectively. The buckling is found to increase from polythiophene to polytellurophene and the C–C bond lengths are found to decrease from polythiophene to polytellurophene, respectively. Whereas the C–X bond length is found to be increasing from polythiophene to polytellurophene, almost identical to the monolayers. The other structural parameters vary in a similar pattern as in case of the corresponding monolayers, compare Table 1.

The interlayer spacinpan>g amounts to 3.42–3.62 Å, 3.39–3.65 Å, anpan>d 3.39–3.68 Å, respectively, inpan> case of scheme (ii). Importantly, the buckling is quite large as compared to that of scheme (i); the probable ren class="Chemical">ason for the increment in the buckling is the interaction between the X atoms, e.g. the S atom of one layer with the Se atom of another overlying layer. This is not the case for scheme (i) because both the layers are identical. In addition, C–S, C–Se, and C–Te bond lengths of 1.73–1.88 Å, 1.74–2.06 Å, and 1.88–2.07 Å are found for S-Se, S-Te, and Se-Te bilayers, respectively, again this is different from scheme (i).

Furthermore, n class="Chemical">polymers that n class="Chemical">contain alternating S, Se and Te atoms in consecutive rings have also been considered, either in a monolayer or bilayer form (referred to S-Se-Te1 and S-Se-Te2 in Tables 1, 2, 3 and 4). It is found that the buckling in the bilayer is about 6 times larger than that of the monolayer because of multi-atomic interactions, i.e. S-S, Se-Se, and Te-Te atomic level interactions (note that AA stacking has been considered). The remaining structural parameters are consistent with the other systems discussed above. Based on our calculations, such systems are promising to reduce the band gap, which will be discussed in detail below.

Electronic structure

n class="Chemical">Polythiophene and its derivatives are a cln class="Chemical">ass of materials that have been widely used as conjugated materials in electronic devices such as optoelectronics9092. Nevertheless, there are certain difficulties to integrate polythiophene in devices mainly due to the limited absorption profile. In order to improve device performance, lots of research activities have been focused on the development of new low band gap polymers made from polythiophene and its derivatives. In Fig. 5, the electronic structures of polythiophene, polyselenophene, and polytellurophene are addressed. The calculated value of the valence bandwidth (which is defined as the energy difference between the first band located right below the Fermi level at Γ-point to the M-point) is found to be almost constant at about 0.74 eV, whereas the conduction bandwidth (which is defined as the energy difference between the first band located right above the Fermi level at Γ-point to the M-point) increases significantly from 0.63 eV to 0.69 eV and 0.81 eV, respectively for polythiophene, polyselenophene, and polytellurophene, see Table 2. Concurrently, the value obtained of the band gap is reduced from 1.12 eV to 0.68 eV to 0.26 eV for polythiophene, polyselenophene, and polytellurophene, respectively, compare Fig. 5(a–c). Whereas the experimental value of the band gap for polythiophene was found to be 2.0 eV, which indeed agrees well with our calculated band gap in the framework of the hybrid (B3LYP) calculations of 2.08 eV. Using the same functional, the obtained values of the band gap are 1.53/0.97 eV for polyselenophene/polytellurophene. It has already been established that the band gaps obtained from the hybrid (B3LYP) calculations reproduce the experimental band gaps very well61. Note that the trend in the band gap using hybrid (B3LYP) calculations (reduction in the band gap) from polythiophene to polytellurophene and their derivatives follows the same trend as the PBE calculations, see Table 2. The values obtained of the band gap for polythiophene and polyselenophene agree well with previous reports628993, however, there is no previous experimentally/theoretically obtained value of the band gap for polytellurophene to be compared. Note that replacing the S atom in polythiophene by heavier atoms such as Se or Te results in lowering the band gap. Previous studies indicated that the HOMO is not contributed by the S atom; on the other hand the LUMO has significant contributions from the S atom9094. Thus, replacement of S by Se or Te should have a large contribution to the LUMO. Because of the lower ionization potential of the Se or Te atoms, their incorporation into polythiophene should lead to a lowered LUMO energy level, and hence, provide for a smaller energy gap as compared to its S counterpart.

Figure 5

The electronic structure of single layers of (a) polythiophene, (b) polyselenophene, and (c) polytellurophene.

Importantly, lowering the band gap along with increasinpan>g the pan> class="Chemical">conduction bandwidth make the material a good candidate to be n-doped, and it should possess a large n-type conductivity as compared with the systems with larger band gap95. Moreover, the conductivity, ionization potential, and electron affinity should be largest for polytellurophene amongst the three materials. Similarly, the effective masses of holes and electrons in the valence as well as conduction bands should follow the same trend as conductivity, ionization potential, and electron affinity. The calculated value of the binding energy [E = (E − E − E − E)/N, where E, E, E, and E are the total energy of the supercell, total energy of isolated C atoms, total energy of isolated S or Se or Te atoms, and total energy of isolated H atoms, while N is the total number of atoms in the system] decreases slightly from polythiophene (−0.483 eV) to polytellurophene (−0.455 eV). Note that, essentially, all the three systems should be easily realized experimentally because of only slight differences in the binding energy. Likewise, the nature of the bonding is expected to be quite similar in the whole network/backbone. Indeed, polytellurophene and polyselenophene have been synthesized34383968969798.

The electronic band stn class="Chemical">ructures of bilayer pan> class="Chemical">polythiophene, polyselenophene, and polytellurophene, where both the layers are identical, are addressed in Fig. 6. The calculated value of the valence bandwidth decreases from polythiophene to polytellurophene, while the conduction bandwidth increases, similar to the monolayer counterparts, see Table 2. The band gap is reduced substantially as compared to the monolayer counterparts; the reduction in the band gap can be understood as a result of the interaction between the two identical layers53. The obtained value of the band gap using the hybrid (B3LYP) functional becomes larger than the value obtained from the PBE functional, as expected, and, importantly, the band gap trend is similar as compared to the PBE functional, see Table 2. A number of synthetic routes have been proposed experimentally to produce various low band gap polymers using bulk or bilayer schemes, including donor-π-acceptor conjugated polymers, for example99, as well as polythiophenes grafted with side chains100 copolymerized with fluorene moieties101. Such low band gap materials are promising for harvesting solar photons102. Therefore, our results point to new forms of conducting polymers that would be useful in solar cells. Moreover, a sizeable band gap is found for two immediate valence and conduction bands; the values of the band gap in the valence as well as conduction bands edges are found to be increasing from polythiophene to polytellurophene, see the data summarized in Table 2. Importantly, there are two possible optical transitions during light absorption and exciton formation for bilayer systems. These are π − π transitions. Such transitions have been observed experimentally for doped π-conjugated polymers and supported by theoretical calculations for doped polythiophene2853. The calculated values of the E are found to be similar to those of the monolayer counterparts, which means that the bonding of the atoms is quite similar for monolayer and bilayer.

Figure 6

The electronic structure of bilayers of (a) polythiophene, (b) polyselenophene, and (c) polytellurophene.

However, a significant variation in the valence and n class="Chemical">conductionpan> banpan>d gap is founpan>d for mixed bilayers of pan> class="Chemical">polythiophene and polyselenophene or polythiophene and polytellurophene or polyselenophene and polytellurophene (compare Fig. 7(a–c)) because of the interaction between the two different layers. This can also be understood from the strong modification in the structural geometry (Fig. 4) as compared to the identical bilayers discussed above. On the other hand, the other parameters including bandwidths, band gap and E follow similar trends to those of the identical bilayers discussed above.

Figure 7

The electronic structure of mixed bilayer systems: (a) polythiophene and polyselenophene, (b) polythiophene and polytellurophene, and (c) polyselenophene and polytellurophene.

In addition, we have investigated the electronpan>ic properties of the pan> class="Chemical">S-Se-Te1 and S-Se-Te2 systems, where the composition of these systems has been discussed above in the structural arrangements section. These systems have the highest values of the E among the systems under consideration. Hence, it is expected that such systems should be energetically very stable, and the constituent atoms should be bonded strongly in the pentagonal backbone. Thus, these systems could be utilized for device application, for example solar cells102. All the calculated parameters are summarized in Tables 1 and 2. The electronic band structure is addressed in Fig. 8(a,b). The value obtained of the band gap is 0.70 eV (1.54 eV using B3LYP functional) for the monolayer and 0.64 eV (1.48 eV using B3LYP functional) for the bilayer with valence and conduction band gaps of 0.17 eV and 0.14 eV for S-Se-Te2. Because of the interlayer interaction, we observe two bands immediately below and above the Fermi level, making the system promising for optical transitions28.

Figure 8

The electronic structure of (a) S-Se-Te1 (b) S-Se-Te2.

Molecular calculations

The structures under pan> class="Chemical">consideration for the molecular calculations are addressed in Fig. 9 along with their n class="Chemical">corresponding HOMO and LUMO orbitals. Whereas, the HOMO and LUMO energies, HOMO-LUMO gap, E, interlayer distance, buckling, bond distances, and bond angles for all the systems are summarized in Tables 3 and 4. Note that there are significant differences in the calculated values of the structural parameters, energies, and HOMO-LUMO gap between periodic and molecular calculation, however both approaches maintain internal consistency. These differences in the parameters can be understood due to the fact that in periodic calculations, the systems are fixed along the y- or z-directions, while in molecular calculations the structures are subjected to modifications in all three directions. Thus, in our opinion, molecular calculations provide slightly better results for comparing with the real experimental situation, especially regarding geometries, whereas, periodic structures are also important because electronic devices has been synthesized using crystalline structures103. Note that, in case of mixed bilayers, the second layer is shifted as compared to the first layer, see Fig. 9. This indicates that in case of scheme (ii) (see section 3.1.1) a slipped ABA stacking pattern is energetically stable as compared to AAA non-slipped stacking pattern in case of scheme (i). In a recent experiment, both types of stacking patterns have been observed for two-dimensional π-conjugated polymers104. The fact is that the unshifted layers are energetically less favourable, which means that this arrangement is only a local minimum on the potential energy surface for these particular systems.

Figure 9

The HOMO and LUMO of the systems under consideration, where yellow, blue, black, green, and gray balls represent C, H, S, Se, and Te atoms, respectively.

The calculan class="Chemical">ted energy of the HOMO is slightly inpan>creasing (decreasing in magnitude) and that of LUMO decreasing (becoming more negative) from polythiophene to polytellurophene and a similar trend is achieved for the HOMO-LUMO gap, see Table 4. The values obtained of the HOMO-LUMO gap are higher (1.51 eV to 1.14 eV for polythiophene to polytellurophene) then those of the periodic approach due to the reasons as discussed above. The E is found to be significantly higher as compared to the periodic approach; it amounts to −0.596 eV, −0.582 eV, and −0.569 eV for polythiophene, polyselenophene, and polytellurophene, respectively. This clearly indicates that the constituent atoms are more strongly bonded with each other in a molecular form as compared to the crystalline form (periodic arrangement along x- and y-directions). The C–C bond lengths are again found to be shrinking upon substitution of the S atom in polythiophene by Se or Te. The C–S, C–Se, and C–Te bond lengths are found to be 1.75 Å, 1.91 Å, and 2.08 Å, respectively. The bond length between C and the heteroatom increases with increasing atomic number due to the larger ionic size of the heteroatoms, in good agreement with previous reports87105106107.

The C–X–C bond angle is found to decren class="Chemical">ase with inpan>creasing size of the heteroatoms (S, Se, and Te) and becomes 92° for polythiophene and 82° for polytellurophene. The other angles are found to be quite similar to those of the periodic approach, compare Tables 1 and 3. Essentially, the obtained parameters are more or less following similar trends as their periodic counterparts. The values obtained of the band gap for polythiophene, polyselenophene, and polytellurophene agree well with previous experimental and theoretical reports53616895.

In case of bilayers composed of two identical layers of polythiophene, polyselenophene or polytellurophene, the averaged value of the interlayer separation amounts to 3.41 Å to 3.53 Å, see Table 3. We observed variable values of the interlayer separation because significant buckling is induced in both the layers due to the interlayer interactions. This agrees well with experimental observation in case of π stacked conjugated polymers108109. The buckling is found to be 0.25 Å to 0.46 Å for polythiophene to polytellurophene, respectively; such a buckling has already been observed for similar systems53. The values obtained of the band gap are reduced significantly (1.16 eV for S, 1.07 eV for Se, and 0.89 eV for Te atoms in the backbone) as compared to those of their monolayer counterparts because of additional molecular orbitals that have developed close to the Fermi level. Similar behavior has been obtained previously for bilayer polythiophene with different chain lengths53.

However, in the case of mixed bilayer (see Fig. 4), variable values of the interlayer spacing, buckling and bond lengths between C atoms and the heteroatoms are found. The reason for obtaining variable values is the interaction between the two different heteroatoms layers. The values obtained of the various parameters are listed in Table 3. The HOMO-LUMO gap is slightly higher then those of identical bilayers (Table 4). The variable value of the band gap which is close to 1.0 eV is promising for optical devices operated in low band gap regimes. Here, we obtain significant and predictable changes in the band gap of the polymers either by substituting heavier heteroatoms in the conjugated backbone or by creating bilayers. The significant lowering of the band gap (HOMO-LUMO gap) is due to the reduced aromaticity of a selenophene ring compared to a thiophene ring. The lower aromaticity in the main ring of the system lowers the HOMO-LUMO gap due to increased contributions of the quinoid structure of the main chain61. A HOMO-LUMO gap lower then that of the monolayers and slightly higher then that of either of the bilayers is achieved for monolayer as well as bilayers of the systems where adjacent S atoms are substitutes by Se and Te atoms, the systems are also called S-Se-Te1 and S-Se-Te2, see Table 4. In summary, the calculations show that polythiophene-based systems can have wide ranges of the HOMO-LUMO gap (band gap), and they could be utilized for various electronic as well as optical devices accordingly. Therefore, synthesis of these systems would be useful.

Polypyrrole and its derivatives (Periodic approach)

The n class="Chemical">polypyrrole-based structures are presented in Fig. 3(b) and right panel of Fig. 4, where X refers to N, P, and As atoms. The supercell of the pristine monolayer of polypyrrole consists of 16 carbon (C) atoms, 12 Hydrogen (H) atoms, and 4 Nitrogen (N) or Phosphorus (P) or Arsenic (As) atoms. The calculated values of the structural parameters using periodic and molecular approaches are listed in Tables 1 and 3, respectively. For polypyrrole, the values obtained of the bond angles agree well with the available report95. In case of polypyrrole, N atom is connected to each of the pentagonal rings of the backbone and one H atom, which is well explored110. Experimentally, it has been suggested that the N atom can be replaced by a P or As atom68. Therefore, we have extended our study by replacing N by P and As, such that these materials could be utilized in electronic devices such as optoelectronics.

For pristine n class="Chemical">polypyrrole, the buckling is found to be 0.11 Å with a N–H bonpan>d lenpan>gth of 1.01 Å. The C–C bonpan>d lenpan>gth is slightly lonpan>ger thanpan> that of the pan> class="Chemical">polythiophene-based structures. The C–N bond length is 1.44 Å, the C-X-C bond angle is 113°, which is much larger than the C-X-C bond angle in the polythiophene counterpart, while the other bond angles are similar to those of polythiophene-based structures, see Table 1. Moreover, buckling, X–H bond lengths, and C–X bond lengths are increasing when replacing N by P and As, see Table 1.

We also exn class="Chemical">tended our study to two difpan> class="Chemical">ferent types of bilayers systems, (i) both the layers are identical, i. e. n class="Disease">N-N, P-P, and As-As, and (ii) the two layers are not identical i.e. N-P, N-As, and P-As, see Fig. 4. Variable interlayer distances of 3.18–3.68 Å, 1.73–2.69 Å, and 1.73–2.72 Å are obtained for the identical bilayer systems N-N, P-P, and As-As with a strong buckling, respectively. This indicates that there are strong structural modifications upon creating bilayers. In case of P and As monolayers or P-P and As-As bilayers, the P-H and As-H bonds point out-of-plane in order to get a stable structure, while the N-H bond points in-plane in case of N monolayer or bilayer of N-N. This is the reason why the interlayer distance in case of P-P and As-As bilayers is significantly smaller (1.73–2.72 Å) than that of the N-N bilayer (3.18–3.68 Å). Pristine polypyrrole bilayers (N-N) have been realized in experiment111. At this point, we believe that the other bilayers should also be realized experimentally, which would widen the research area. Experimentally it has been proposed that the bilayers made of polpyrrole (N-N) are promising for applications in amperometric ion sensing112. The buckling is found to be smaller and slightly less structural modifications are obtained in case of bilayers made of two non-identical layers (N-P, N-As, and P-As), see Table 1 and Fig. 4. Moreover, the electronic structure of the structures presented in the right panel of Fig. 4 are addressed in Fig. 10. Note that, in case of bilayers, the second layer is shifted as compared to the first layer, which indicates that the ABA kind of stacking is energetically favourable in these cases similar to bilayers made of polythiophene-based systems, as discussed above104.

Figure 10

The electronic structure of a single layer of pyrrole derivatives with (a) N, (b) P, and (c) As.

n class="Chemical">Polypyrrole is a wide band gap organic pan> class="Chemical">polymer with many applications, such as flexible batn class="Chemical">teries, semiconductors, large-area optical displays, and sensors113. It has a larger band gap and hence, lower conductivity in the neutral state as compared to polythiophene. It has been reported that the effective mass of holes and electrons is comparatively higher than that of other polymers such as polythiophene, polyacetylene, and polyfuran95. The calculated band structures of polypyrrole and its P and As analogues are presented in Fig. 10(a–c). The calculated values of the valence and conduction bandwidths along with the band gaps are summarized in Table 2. The value obtained of the band gap is found to be 1.88 eV, which agrees well with a previously calculated value using first-principles calculations95. The band gap shrinks significantly and becomes 0.34 eV and 0.61 eV for N atom replaced by P atom and As atom (see Fig. 10(b,c)), respectively. Moreover, the value of the band gap is found to be larger as compared to the PBE functional when using the hybrid (B3LYP) functional, as expected, see Table 2 for the comparison. This fact is true for all the calculations for polypyrrole and its derivatives. As discussed above, such a replacement is quite feasible in a real experimental situation70. A similar trend is found for the E as compared to the corresponding polythiophene-based monolayers.

The calculated electronic banpan>d stpan> class="Chemical">ructures for bilayers, where both the layers are identical, are addressed in Fig. 11. It is clear from the figure that there are two bands immediately above and below the Fermi level with a finite gap, contributed from each of the layers. The valence and conduction bandwidths are growing from the N-N to As-As bilayers, while the band gap shrinks from N-N to As-As. A similar trend is also found to the corresponding polythiophene-based bilayers discussed above. The E again follows the same trend as for the corresponding monolayers. Because of strong structural deformation, non-identical bilayers made of N-P and P-As do not provide two bands immediately above the Fermi level; rather they provide two bands immediately below the Fermi level. Whereas, in case of the N-As bilayer, two bands immediately above and below the Fermi level are obtained, which is due to less modification in the structure, compare Fig. 12(a–c). Moreover, the other parameters including bandwidths, band gap and E have similar trends to those of the bilayers made of two identical layers.

Figure 11

The electronic structure of identical bilayers of (a) N-N, (b) P-P, and (c) As-As.

Figure 12

The electronic structure of mixed polypyrrole derived bilayers, (a) N-P, (b) N-As, and (c) P-As.

The structures under pan> class="Chemical">consideration are shown in Fig. 13 with their corresponding HOMO and LUMO. The HOMO and LUMO energies, HOMO-LUMO gap, E, distance between the two layers, buckling, bond distances, and bond angles for all the structures under investigation are presented in Table 3. The electronic and structural parameters calculated using the molecular approach are slightly different from those of the periodic approach. The reason behind that is again the fact that in the periodic calculations, the systems are forced not to be moved along the y- or z-directions whereas in molecular calculations this is not the case. The values obtained of the HOMO and LUMO energies are significantly increasing in magnitude from N to P or As, which was not the case for the polythiophene-based systems, compare Table 4. The value obtained of the HOMO-LUMO gap for polypyrrole (N) is found to be 2.25 eV and decreases to as low as 0.89 eV and 0.81 eV for N replaced by P and As in polypyrrole, respectively. The HOMO-LUMO gap is larger than that calculated using the periodic approach (1.88 eV) and close to the experimentally and theoretically obtained value of the HOMO-LUMO gap of 2.8 eV95114115.

Figure 13

The HOMO and LUMO of the systems under consideration.

The E is found to be significantly larger than that obtained from the periodic approach, indicating a strong bonding of atoms in this approach. The E are −0.609 eV, −0.562 eV, and −0.548 eV for N, P, and n class="Chemical">As sysn class="Chemical">tems, respectively, see Table 2. The buckling is also irregular and higher than that of the periodic approach. Whereas, other structural parameters such as bond lengths and bond angles are found to be similar to the values obtained using the periodic approach. For bilayers (made of the two identical layers), a slight reduction in the HOMO and LUMO energies is achieved with slightly smaller HOMO-LUMO gaps of 2.20 eV, 0.89 eV, and 0.76 eV for bilayers of N-N, P-P, and As-As, respectively. Note that the HOMO is contributed by one of the layers and the LUMO is contributed by the other layer in case of N-N and P-P, while in case of As-As, the HOMO and LUMO both are contributed from both the layers, in agreement with previous findings53. Due to the very strong structural modification and corrugation, we observed variable values of the interlayer separation with significant buckling; this agrees well with experimental observation in case of π stacked conjugated polymers as discussed above108109. The trends for the other structural parameters remain similar to those obtained from periodic calculations.

Moreover, in cn class="Chemical">ase of mixed bilayers n class="Chemical">N-P, N-As, and P-As (Fig. 10), the N-P bilayer system has significantly larger LUMO energy as compared to N-N but N-As and P-As bilayers have similar trends of HOMO or LUMO to the P-P and As-As bilayers. The value obtained of the HOMO-LUMO gap is strongly reduced (0.41 eV, 0.40 eV, and 0.73 eV, for N-P, N-As, and P-As bilayers, respectively) as compared to their identical bilayer counterparts (see Table 2). The strong reduction in the HOMO-LUMO gap can be understood from the different interlayer interactions and structural modifications. Again in these bilayers, the HOMO is contributed from one layer and the LUMO is contributed from the other layer. The E is slightly smaller than that of the identical counterparts, and larger interlayer separation and smaller buckling are found. Variable X–H bonds are found due to strong structural modification. All other structural parameters are found to be following similar trends as compared to N-N, P-P, and As-As bilayers.

Bilayer of polythiophene and polypyrrole

In addition, we further studied the structural anpan>d electronic properties of mixed bilayers made of n class="Chemical">polythiophene and n class="Chemical">polypyrrole (S-P). The structural and electronic parameters using both the periodic and molecular approaches are summarized in Tables 1 and 2. Valence and conduction bandwidths of 0.60 eV and 0.67 eV, a band gap of 0.72 eV using PBE functional and 1.72 eV using hybrid (B3LYP) functional, a band gap immediately above the Fermi level of 1.03 eV, and a band gap immediately below the Fermi level of 0.41 eV are obtained. Other structural parameters are similar to those of other systems under study. A HOMO-LUMO gap of 0.87 eV is found in the molecular calculation, which is close to that obtained in the periodic calculation (0.72 eV), see Tables 1 and 2. Other structural parameters (bond lengths and bond angles) have again similar trends as other systems under study. It should be noted that such systems are promising to achieve low band gap π-conducting polymers, and could therefore be utilized in electronic devices9092.

Conclusion

Using first-principles calculations, the structural anpan>d electronpan>ic properties (the band gap modulationpan>) of pan> class="Chemical">polythiophene and polyprrrole-based systems (in periodic as well as in oligomer approaches) have been investigated. The effects in the structural and electronic properties for S substituted by Se and Te in case of polythiophene and N substituted by P and As for polypyrrole have been addressed. Furthermore, the effect of the stacking in pure polythiophene as well as in pure polypyrrole have been studied and a linear reduction of the band gap or HOMO-LUMO gap as function of the stacking was found. In addition, we also employed hybrid functional based calculations for periodic cases in order to understand the effect of the hybrid functional on the band gap. It was found that the band gap is larger using hybrid (B3LYP) calculations as expected, while the trend in the band gap is similar to that of the PBE calculations.

This study hn class="Chemical">as also been expan> class="Chemical">tended to polythiophene and n class="Chemical">polytellurophene and to polypyrrole derivatives where the N atom is substituted by P and As atoms. In addition, we performed Γ-point phonon calculation for selected systems to check the stability of these systems. The values obtained of the E suggest that all the systems under consideration are thermodynamically stable and can therefore be realized experimentally. The trends in the E as a function of stacking as well as substitution remain similar for both the periodic and molecular approaches. Moreover, the bilayer constructed from polythiophene on top of polypyrrole has also been investigated. It is believed that the comprehensive analysis of this investigation will be useful for fundamental understanding of the structural and electronic properties and the band gap or HOMO-LUMO gap, which could be very important for constructing polythiophene and polypyrrole-based electronic devices because a wide range of the band gap can be obtained.

Additional Information

How to cin class="Chemical">te this article: Kaloni, T. P. et al. Band gap modulation in n class="Chemical">polythiophene and polypyrrole-based systems. Sci. Rep. 6, 36554; doi: 10.1038/srep36554 (2016).

Publisher’s non class="Chemical">te: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.