Research Article  Open Access
Theoretical Study on the Static (Hyper)Polarizabilities of αtBuωCNPoly(methylphenyl)silane
Abstract
The static linear and nonlinear optical properties of the σconjugated polymer αtBuωCNpoly(methylphenyl)silane (PMS) are studied at the CoupledPerturbed HartreeFock (CPHF) level with 631+G(d) basis set. The calculated results reveal that the static first hyperpolarizabilities of this system increase with the main chain length and have a good agreement with experiments. The (hyper)polarizabilities per unit cell have been extrapolated to infinite chain limit and a comparison is made to those of polysilane and polyacetylene (PA). Besides, other structural properties depending on the σconjugated Si–Si skeleton length are investigated as well. Electron correlation effect is estimated and it turns out that the MP2 static first hyperpolarizability is about times larger than the corresponding CPHF value for the polymer with .
1. Introduction
In recent years, polysilanes have attracted increasingly extensive attention due to their unique physical and chemical properties resulting from σelectrons delocalized along the silicon backbone. In this regard, they have resulted in a variety of technological applications, such as conductor and semiconductor [1], photoconductive materials [2], organic multilayer, light emitting diodes (LEDs) [3], highdensity optical data storage materials [4], electro luminescence (EL) devices [5, 6], and nonlinear optical (NLO) materials [7, 8].
The polysilanes represent a completely new class of potentially interesting nonlinear optical (NLO) materials, and in contrast to the classical πconjugated polymers, their physical, chemical, and optical properties show significant differences. Polysilanes are linear polymers of silicon, and the σelectrons of the polymer backbone are delocalized. This σconjugation gives rise to electronic properties that allow for possible applications as electroluminescent, nonlinear, optical, lithographic, and semiconductor materials. The σelectrons delocalization associated with the phenomenon, in which the silicon backbone itself is a chromophore, greatly influences the optical properties of polysilanes. Therefore, in spite of saturation of all the Si atoms on the main chain, polysilanes have strong absorption in the UV in the region of 250–400 nm, which is attributed to the σσ* transition or the σ^{3d} transition and relevant to the configuration and length of the Si–Si backbone. As NLO materials, in comparison with πconjugated carbonbackbone polymers, polysilanes own many unique advantages; for example, they have optical transparency in the visible spectrum, their UV absorption peak can be able to be shifted, they are easily to be made into films and applied to optic research, and they are of good thermal and chemical stability [9].
Many substituted polysilanes have been investigated from experimental points of view and the majority of them focus on the thirdorder optical nonlinear coefficient . Shi et al. [10] have synthesized poly(methylphenyl)silane (PMPS) and measured its thirdorder nonlinear optical susceptibility of 5.8 × esu by third harmonic generation (THG) with a pulse laser width at 1064 nm for thin films. Then, Tang and coworkers [11] have explored a novel strategy to develop polysilanes functionalized with a high density of nonlinear optical chromophores and the poled film reveals a resonant value of 19 pm/V by second harmonic generation (SHG) measurement at the frequency of 1064 nm. However, what would the situation of the secondorder nonlinear optical susceptibility of PMPS derivatives be? There is little scientific attention from the theoretical point of view.
For an asymmetric molecule, its secondorder nonlinear optical susceptibility is nonzero. In this paper, electric properties of PMS are studied in order to assess the trends with increasing Si–Si chain length for the cases , and 14. Meanwhile, the mechanism of the variation of the static (hyper)polarizability as a function of chain length is elucidated by analyzing the Si–Si bond lengths and natural population analysis (NPA) charges on Si atoms along the backbone as well as energy gaps between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO).
2. Computational Details
All of the calculations have been performed through the following procedures.(1)The conformations of the polymers were optimized for alltrans conformation, which was assumed to be the ground state geometry, at the B3LYP/631G(d) level of theory. We constrained the Si–Si skeleton in a zigzag plane and assumed syndiotactic structure throughout the polymers. The longitudinal axis was set by the line connecting the first and second to last silicon atoms of the main chain. With the optimized geometries, NPA charges of the polymers were investigated at the same level as B3LYP/631G(d).(2)MP2 method was employed for calculating shorter chains of PMS to estimate the electron correlation effects on NLO properties of PMS as it is known that conventional DFT method is unsuitable for describing the chain length dependence of the hyperpolarizability in conjugated systems [12, 13]. Thus, CPHF method has been employed to compute the static polarizability and first hyperpolarizability.
In the presence of an electric field, energy can be expressed as where is the energy without the applied electrostatic field, is the component of the field along the th Cartesian direction, and , , and may be called a component of dipole, polarizability, and the firstorder hyperpolarizability tensor, respectively. Based on the computed tensor components, we focus on , , and . Their expressions are written as follows: where
All calculations were performed with the GAUSSIAN 03 [14] program package. Net atomic charges were obtained using the NPA included in the natural bond orbital (NBO) algorithm [15, 16]. The plots of molecular orbitals were generated with the GaussView program (Gaussian, Inc., Wallingford, CT, USA) [17].
3. Results and Discussion
3.1. Conformations
The schematic illustration of optimized conformations of PMS is shown in Figure 1. The important geometrical parameters for , and 14 are listed in Table 1. As seen in Table 1, for a certain polymer molecule, the longest Si–Si bond is always in the central part of it. Meanwhile, from to , with the increase of Si main chain length of the system, one can find that the Si–Si bond length is increased. The variation tendency of the central Si–Si bond length is depicted in Figure 2. The central Si–Si bond length ranges from 2.393 to 2.411 a.u., and it is almost convergent when is greater than 10. From this tendency, one can deduce that the Si–Si bond dissociation energies are getting smaller and smaller; consequently the σσ* transitions become easier as the Si skeleton length increases. However, the bond lengths between substituents and Si atoms at both edges of the backbones (tBu–Si, Si–CN) are barely changed. Compared with polysilanes, whose Si–Si bond length is around 2.36 Å, PMS polymer exhibits relatively longer Si–Si bond lengths (2.38–2.41 Å), because its main Si–Si chain is positively charged by the influence of the larger electron affinity of phenyls.

The Si–Si–Si bond angles of the PMS polymer are ranging from 112.5 to 113.8 degree. At the two ends of the chain, the bond angles become a little bit smaller.
3.2. Charge Distribution
In order to understand the electronic properties of PMS, atomic charges were calculated by means of NPA at B3LYP/631G(d) level. Table 2 shows the atomic charges of Si atoms and the group charges of tBu and –CN of the PMS. From the data in Table 2, one can see that regardless of the chain length, the atomic charges of Si on one end of the Si–Si backbone, which is connected to the tBu group, are near 1.41 a.u., and while the atomic charges of Si on the other end, which is connected to the –CN group are near 1.25 a.u., the latter is smaller than the former. A proper explanation of this result could be given that the charge transfer from electrondonating group to electronaccepting group along the Si–Si backbone has taken place so that the atomic charges of Si connecting to the electrondonating group tBu are more positive than those of Si connecting to the electronaccepting group –CN. In the central region of the chain, there is a concentration of charge distribution and the curve of the average charges of the middle two Si atoms is shown in Figure 3. The values are ranging from 0.84 to 0.88 a.u., and it is almost convergent when is greater than 6. On the other hand, charges on the both groups of tBu and –CN are getting more and more negative with increasing; however, charges on –CN are always more negative than those of tBu and their differences of are decreased.

3.3. The Frontier Orbital Energies
The frontier orbital energies of PMS were evaluated by B3LYP method with the 631G(d) basis set and listed in Table 3. With the increase of repeating unit, the energies of HOMO are increased and the energies of LUMO are decreased; thus the energy gaps () between LUMO and HOMO are decreased. From the schematic illustration of HOMO and LUMO as shown in Figure 4, one can see that the HOMO is localized on both the main chain and phenyls and is combined of σ and π type; while the LUMO is antibonding σ* orbital. Therefore, it is expected that when the Si–Si main chain gets longer, the σσ* conjugation along the Si–Si backbone will be stronger.

(a)
(b)
3.4. Static Electronic (Hyper)Polarizabilities
3.4.1. Basis Sets Effect
Different basis sets were chosen as following: 3–21G, 631G, 631G(d), 631G(d, p), 631+G(d), 631+G(d, p), 631++G(d, p), and 6311++G(d, p), to study the effect of basis sets on dipole , polarizability , and the first hyperpolarizability . The results are listed in Table 4.

From the data, one can see that the choice of basis sets has a great influence on the results. For and , with the increase of the size of basis set, both of them become bigger. For , however, three basis sets 631+G(d), 631+G(d, p), and 631++G(d, p) provide nearly the same results. Comparing basis sets 631G and 631G(d) with 631G(d, p), we find that becomes larger when polarization functions are introduced, but the second polarization function has little effect on it. Comparing basis sets 631G(d) with 631+G(d) and 631+G(d, p) with 631++G(d, p), it is found that becomes smaller when diffuse functions are added; however, the second diffuse function also has little affect on . This implies that basis sets of H atoms contribute only a bit to the first hyperpolarizability . Based on these analyses, the basis set of 631+G(d) was chosen, which could be a good one for less time consuming and relatively reliable results.
3.4.2. Electron Correlation Effect
The first hyperpolarizabilities of PMS () were calculated by both CPHF and MP2 method, using 631+G(d) basis set to discuss the electron correlation effect. The results are listed in Table 5. From Table 5, one can see that the electron correlation effect has much significant influence on the first hyperpolarizability compared with basis sets effect. When , MP2 is 1.44 times larger than that of CPHF; when , the multiple rises to 1.79; and when , the corresponding multiple becomes 2.05. So one can estimate that the electron correlation effect on the first hyperpolarizability is about 1.5 ~ 2.0 times larger than that of CPHF for polymer .

3.4.3. Static (Hyper)Polarizabilities
The CPHF/631+G(d) static first hyperpolarizabilities of PMS (), as calculated for different chain lengths, are given in Table 6. From this table it can be seen that (i) component, which represents for the value in the longitudinal direction of the Si–Si skeleton, makes the dominant contribution to and increases much faster than the others; (ii) , , and increase as the chain length elongated; however, the amount of increase for each component is discriminating. Among them, has the biggest rate of increase, which ranges from 354.10 to 1564.33 a.u. from to ; is ranging from 77.2 to 306.0 a.u.: it first increases quickly and then slow down when increases. The ’s increase is the minimum, rising from 1.95 to 2.09 a.u.
 
; . 
The polarizability and first hyperpolarizability per unit cell of alltrans PMS were evaluated as with or . Thus or represents the or of each unit cell of –(SiCH_{3}C_{6}H_{5})–. They were extrapolated to the infinite chain limit. Assuming asymptotic limit as approaches infinity, the values of and have been fitted to where , , , and are the fitting parameters. The extrapolated value for infinite polymer is thus .
For a physical quantity and a certain fitting function, there might be different fitted limiting value when different fitting data ranges are chosen. The optimal fitting data ranges for are chosen from 4 to 14, whereas it is chosen from 6 to 14 for . The fitting parameters and extrapolated values are listed in Table 7. The corresponding fitting curves of and can be found in Figure 5. The fitting limit (coefficient of ) for and is approximately equal to 103.1 a.u. and 5.3 a.u., respectively. To make a comparison, Kirtman and Hasan [18] have calculated the linear polarizabilities of alltrans polysilane with 631 G basis set at CPHF level and it was found that the longitudinal linear polarizability per SiH_{2} unit is 65.6–65.8 a.u. In addition, Hurst et al. [19] have researched on large πconjugated organic molecule PA at the level of CPHF/631G(d, p), and the corresponding value is 60.8 per CH unit. The polarizability of –(SiCH_{3}C_{6}H_{5})– unit is greater than that of SiH_{2} unit and CH unit.

(a)
(b)
For experiments, Li et al. [20] have successfully synthesized several kinds of functionalized polysilanes based on PMPS and containing NLO chromophores in the side chain. One of the synthesized structures PII1 (see in Figure 6) has some similarity to that of PMS polymer because both of them have the –(SiCH_{3}C_{6}H_{5})– unit. The macroscopic secondorder nonlinear susceptibility d_{33} for PII1 was measured by SHG experiments and its value is 13 pm/V (31.0 × 10^{−9} esu) at 1064 nm fundamental wavelength, which is about 9 times larger than the experimental measurement for urea crystal (3.4 × 10^{−9} esu) [21]. of urea crystal as calculated by Wu et al. [22] is about 27.4 × 10^{−32} esu by TDCPHF method with 631+G basis set. In our system, when , is 2.6 × 10^{−30} esu, which is around 9 times larger than that of urea. Therefore, our calculated value should be close to the magnitude of the experimental result. If one considers the electron correlation effect, is estimated around 8.0 × 10^{−30} esu, which is more than 30 times larger than that of urea.
4. Conclusions
In this work, we present a quantumchemical theoretical analysis of the static first hyperpolarizabilities of a series of αtBuωCNpoly(methylphenyl)silanes at the level of ab initio CPHF method. The effects of basis sets and electron correlation on the first hyperpolarizability are discussed. The basis set, 631+G(d), is chosen for the CPHF calculation. The calculated results are also sensitive to electron correlation effect; the MP2 results are about 1.5 ~ 2.0 times larger than those of CPHF for polymer .
Dependences of NLO properties on polymer chain length are considered. It turns out that in the donorbridgeacceptor structure, there is a concentration of charge distribution in the central region of the chain and the atomic charge on Si at the end of acceptor is smaller than that of donor. With the increase of polymerization degree, the average values of Si–Si bond length are increased and the energy gaps between LUMO and HOMO are decreased; so the σσ* conjugation effect is gradually enhanced. When , the value reaches to 2.6 × 10^{−30} esu at CPHF level (about 8.0 × 10^{−30} esu after MP2 correction), which should be close to the magnitude of the experimental result.
Conflict of Interests
All the authors do not have conflict of interests with the mentioned commercial identities as confirmed.
Acknowledgments
The authors are grateful to the financial support from the National Science Foundation of China (no. 21073067). All the calculations were carried out on the Linux cluster systems in their laboratory at SCNU.
References
 S. Mimura, H. Naito, Y. Kanemitsu, K. Matsukawa, and H. Inoue, “Optical properties of (organic polysilane)(inorganic matrix) hybrid thin films,” Journal of Luminescence, vol. 87, pp. 715–717, 2000. View at: Publisher Site  Google Scholar
 K. Matsukawa, T. Tamai, and H. Inoue, “Photoconductive properties of polysilane copolymers with pendant siloxane groups,” Applied Physics Letters, vol. 77, no. 5, pp. 675–677, 2000. View at: Google Scholar
 A. Sharma, M. Katiyar, Deepak, S. Seki, and S. Tagawa, “Room temperature ultraviolet emission at 357 nm from polysilane based organic light emitting diode,” Applied Physics Letters, vol. 88, no. 14, Article ID 143511, 2006. View at: Publisher Site  Google Scholar
 M. Hiramoto, Y. Sakata, and M. Yokoyama, “Patterned emission in organic electroluminescent device using photodecomposition of polysilane film by UV light,” Japanese Journal of Applied Physics A, vol. 35, no. 9, pp. 4809–4812, 1996. View at: Google Scholar
 L. Sacarescu, I. Mangalagiu, M. Simionescu, G. Sacarescu, and R. Ardeleanu, “Polysilanes. A new route toward high performance EL devices,” Macromolecular Symposia, vol. 267, no. 1, pp. 123–128, 2008. View at: Publisher Site  Google Scholar
 K. Ebihara, S. Matsushita, S. Koshihara et al., “Ultraviolet luminescence and electroluminescence of polydihexylsilane,” Journal of Luminescence, vol. 72–74, pp. 43–45, 1997. View at: Google Scholar
 F. Kajzar, J. Messier, and C. Rosilio, “Nonlinear optical properties of thin films of polysilane,” Journal of Applied Physics, vol. 60, no. 9, pp. 3040–3044, 1986. View at: Publisher Site  Google Scholar
 J. C. Baumert, G. C. Bjorklund, D. H. Jundt et al., “Temperature dependence of the thirdorder nonlinear optical susceptibilities in polysilanes and polygermanes,” Applied Physics Letters, vol. 53, no. 13, pp. 1147–1149, 1988. View at: Publisher Site  Google Scholar
 D. B. Chen, G. Q. Li, and Z. H. Zhou, “Polysilanesa class of new and distinct nonlinear optical materials,” Chemical Research and Application, vol. 5, no. 2, pp. 8–14, 1993. View at: Google Scholar
 B. C. Shi, M. G. Zhang, G. W. Li, Y. Y. Li, M. Y. Wang, and Z. Li, “The synthesis of poly(methylphenyl)silane and investigation on its nonlinear optical properties,” Journal of NanJing Normal University, vol. 20, no. 3, pp. 26–29, 1997. View at: Google Scholar
 H. D. Tang, J. D. Luo, J. G. Qin, H. Kang, and C. Ye, “A novel synthetic strategy to develop secondorder nonlinear optical polysilanes for potential photorefractive effects,” Macromolecular Rapid Communications, vol. 21, no. 16, pp. 1125–1129, 2000. View at: Google Scholar
 B. Champagne, E. A. Perpète, D. Jacquemin et al., “Assessment of conventional density functional schemes for computing the dipole moment and (hyper)polarizabilities of pushpull πconjugated systems,” Journal of Physical Chemistry A, vol. 104, no. 20, pp. 4755–4763, 2000. View at: Google Scholar
 M. TorrentSucarrat, M. Solà, M. Duran, J. M. Luis, and B. Kirtman, “Basis set and electron correlation effects on ab initio electronic and vibrational nonlinear optical properties of conjugated organic molecules,” Journal of Chemical Physics, vol. 118, no. 2, pp. 711–718, 2003. View at: Publisher Site  Google Scholar
 M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., Gaussian 03, Revision E. 01, Gaussian Inc., Wallingford, Conn, USA, 2004.
 A. E. Reed and F. Weinhold, “Natural bond orbital analysis of nearHartreeFock water dimer,” The Journal of Chemical Physics, vol. 78, no. 6, pp. 4066–4073, 1983. View at: Google Scholar
 A. E. Reed, R. B. Weinstock, and F. Weinhold, “Natural population analysis,” The Journal of Chemical Physics, vol. 83, no. 2, pp. 735–746, 1985. View at: Google Scholar
 R. D. Dennington II, T. A. Keith, J. Millam, A. B. Nielsen, A. J. Holder, and J. Hiscocks, GaussView, Version 5.0, Gaussian Inc., Wallingford, Conn, USA, 2009.
 B. Kirtman and M. Hasan, “Linear and nonlinear polarizabilities of transpolysilane from ab initio oligomer calculations,” The Journal of Chemical Physics, vol. 96, no. 1, pp. 470–479, 1992. View at: Google Scholar
 G. J. B. Hurst, M. Dupuis, and E. Clementi, “Ab initio analytic polarizability, first and second hyperpolarizabilities of large conjugated organic molecules: applications to polyenes C_{4}H_{6} to C_{22}H_{24},” The Journal of Chemical Physics, vol. 89, no. 1, pp. 385–395, 1988. View at: Google Scholar
 J. Li, Z. Li, H. Tang, H. Zeng, and J. Qin, “Polysilanes with NLO chromophores as pendant groups by utilizing different synthetic strategies,” Journal of Organometallic Chemistry, vol. 685, no. 12, pp. 258–268, 2003. View at: Publisher Site  Google Scholar
 J. M. Halbout and C. L. Tang, “Properties and applications of ureaeds,” in Nonlinear Optical Properties of Organic Molecular and Crystals, D. S. Chemla and J. Zyss, Eds., Academic Press, New York, NY, USA, 1987. View at: Google Scholar
 K. C. Wu, R. J. Sa, J. Li et al., “Theoretical study on the linear and nonlinear optical polarizabilities of urea clusters,” Jiegou Huaxue, vol. 20, no. 4, pp. 323–324, 2001. View at: Google Scholar
Copyright
Copyright © 2013 Jing Li et al. 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.