Theoretical investigation of the static (dynamic) polarizability and second hyperpolarizability of DAAD quadrupolar push–pull molecules. A comparison among HF (TD-HF), DFT (TD-B3LYP), and MP2 (TD-MP2) methods (original) (raw)
Related papers
The Journal of chemical physics, 2008
The donor/acceptor (D/A) substituted π-conjugated organic molecules possess extremely fast nonlinear optical (NLO) response time that is purely electronic in origin. This makes them promising candidates for optoelectronic applications. In the present study, we utilized four hybrid density functionals (B3LYP, B97-2, PBE0, BMK), Hartree–Fock, and second order Møller–Plesset correlation energy correction, truncated at second-order (MP2) methods with different basis sets to estimate molecular first hyperpolarizability (β) of D/A-substituted benzenes and stilbenes (D = OMe, OH, NMe2, NH2; A = NO2, CN). The results of density functional theory (DFT) calculations are compared to those of MP2 method and to the experimental data. We addressed the following questions: (1) the accurate techniques to compare calculated results to each other and to experiment, (2) the choice of the basis set, (3) the effect of molecular planarity, and (4) the choice of the method. Comparison of the absolute values of hyperpolarizabilities obtained computationally and experimentally is complicated by the ambiguities in conventions and reference values used by different experimental groups. A much more tangible way is to compare the ratios of β’s for two (or more) given molecules of interest that were calculated at the same level of theory and measured at the same laboratory using the same conventions and reference values. Coincidentally, it is the relative hyperpolarizabilities rather than absolute ones that are of importance in the rational molecular design of effective NLO materials. This design includes prediction of the most promising candidates from particular homologous series, which are to be synthesized and used for further investigation. In order to accomplish this goal, semiquantitative level of accuracy is usually sufficient. Augmentation of the basis set with polarization and diffuse functions changes β by 20%; however, further extension of the basis set does not have significant effect. Thus, we recommend 6-31+G* basis set. We also show that the use of planar geometry constraints for the molecules, which can somewhat deviate from planarity in the gas phase, leads to sufficient accuracy (with an error less than 10%) of predicted values. For all the molecules studied, MP2 values are in better agreement with experiment, while DFT hybrid methods overestimate β values. BMK functional gives the best agreement with experiment, with systematic overestimation close to the factor of 1.4. We propose to use the scaled BMK results for prediction of molecular hyperpolarizability at semiquantitative level of accuracy.
Conformational dependence of the second hyperpolarizability of quadrupolar molecules
Journal of Molecular Structure-theochem, 2009
Hatree–Fock calculations at ab initio and semiempirical levels were carried out for the averaged polarizability 〈α〉 and second hyperpolarizability 〈γ〉 of two pairs of quadrupolar isomers with different donor and acceptor groups. These properties were correlated with the antibonding/bonding π occupation number (π∗/π ratio). It was found that isomers with extended π systems had low π∗/π ratios and high 〈α〉 and 〈γ〉 values, while low 〈α〉 and 〈γ〉 values were obtained for isomers with large π∗/π ratios and no extended π system. The PM3 and PM6 〈α〉 values were found to be in excellent agreement with the HF/6-31+G(d,p) ones. The PM3 values for 〈γ〉 were significantly larger than those calculated by HF/6-31+G(d,p), with an average PM3/HF ratio of 1.43. The PM6 results were noticeably better with a ratio of 0.85. The calculation of 〈α〉 and 〈γ〉 at MP2/6-31+G(d,p) level for representative isomers showed that the contribution of the electron correlation to their values was small and that the HF/6-31+G(d,p) method provides reliable values at much lower computational cost.
2016
A series of divinyl-pyrene derivatives of the form D-vinyl-pyrene-vinyl-A, in which D corresponds to an electron donor group and A to an electron acceptor group, were studied in this work. The first purpose was to determine the optimal HF % exchange as incorporated in a range of hybrid functionals (M06HF, M062X, M06L, CAM-B3LYP, PBE0, BMK, and B3LYP) capable to produce, reliably and as close as possible to those obtained from MP2 calculations, NLO parameters and, in particular, first-order static hyperpolarizabilities. The CAM-B3LYP functional was revealed to be the most suitable one. The pair N(CH 3) 2 /NO 2 was then determined as the most efficient pair of groups in producing appreciable NLO responses. The effect of the substitution position on the pyrene moiety was also investigated, whereby aligning the two substituents involving the D and A groups in the direction of the dipole moment as in the (1,6 DVP) derivatives was shown to be most favorable for increasing the NLO parameters. ■ INTRODUCTION Because of their numerous applications in the areas of photonics and optoelectronics, 1−3 various types of NLO active materials have been developed over the recent years. These materials are used in technical applications like lasers, 4 optical data storage devices, biological imaging, 5 signal processing, and in electro-optical modulation as well. 6,7 These materials range from being inorganic to organic to organometallic. 8−11 In order to identify additional and more efficient materials for further NLO applications, many studies were conducted on organic materials, including, for instance, the polycyclic aromatic hydrocarbons (PAH) or, polynuclear aromatic hydrocarbons. 12−15 These latter are a series of hydrocarbons in which carbon atoms are arranged in multiple aromatic rings sharing one or more C−C bond. The resulting structures are molecules wherein all carbon and hydrogen atoms are lying in the same plane. Accordingly, their physical and chemical characteristics make them eligible as π-conjugated bridges 16 or chromophores in push−pull systems, 17 which contain a donor (D) and an acceptor (A) group end-capping the π-conjugated backbone. 18−20 Thereby, an enhancement of the intramolecular charge transfer (ICT) is obtained, 21−23 which results in a strong NLO response. The latter is closely connected to several factors such as the energy gap (ΔE H-L = ε LUMO − ε HOMO), the occupation quotient π*/π, the central bond length, the bond-length alternation (BLA), the delocalization energy E(2), the electronic transitions parameters, the polarizability α, and, finally, the most important factors: the static first order hyper-polarizability β 0 and the dipole moment μ 0 , 24−26 which have been at the center of much interest from theoretical chemists 27−34 for the past couple of decades. In this field, pure and hybrid DFT methods have been demonstrated to fail in calculating reliable hyperpolarizabilities, 35−38 while post-HF methods reproduce hyperpolarizabilities with a good accuracy 39 and are considered in the present work as a reference. The wrong asymptotic behavior of DFT methods 40 has been shown to be at the origin of the resulting inaccurate values. 41,42 Therefore, enhancing the incorporated HF exchange (X) fraction in hybrid DFT fuctionals as well as long-range corrections become very important in order to improve results, and many papers have appeared dealing with this goal. 43,44,37,38 Despite the popularity of the Minnesota M06 functionals, 45,46 their performance in calculating NLO responses of push−pull molecules has not been extensively tested. In this work, we focus on calculating NLO responses of various divinylpyrene derivatives that we consider as sufficiently large systems to test to what extent DFT methods can produce satisfactory results. Actually, as a special case of a PAH, we have pyrene, which consists of four fused benzene rings 47,48 (Figure 1). Since some years, pyrene and its derivatives have attracted considerable attention in the field of NLO because of their remarkable ICT. Actually, they have been frequently used as materials in organic light emitting diodes (OLED), in organic field effect transistors (OFET), and in organic photovoltaic cells. 49,50 So far, some 1-monosubstituted pyrene derivatives have been studied using the DFT-B3LYP/6-31g (d,p) level of theory. 51 It was demonstrated that the pyrene structure is strongly affected by the nature of the substituent. These variations are then obviously reflected in appreciable changes in the dipole moments values. Another work also using the B3LYP/6-31g (d,p) and dealing with some 1,3,6,8 tetra-hologeno-substituted