Next Article in Journal
Advances in Nanomaterials for Photovoltaic Applications
Previous Article in Journal
Criterion for Selective Area Growth of III-V Nanowires
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design of A-D-A-Type Organic Third-Order Nonlinear Optical Materials Based on Benzodithiophene: A DFT Study

1
Gansu Provincial Engineering Research Center for Organic Semiconductor Materials and Application Technology, School of Materials Science and Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China
2
Key Laboratory of Optoelectronic Technology and Intelligent Control of Ministry Education, Lanzhou Jiaotong University, Lanzhou 730070, China
3
School of Chemical Engineering, Northwest Minzu University, Lanzhou 730030, China
4
Key Laboratory of Physical Electronics and Devices of the Ministry of Education & Shaanxi Key Lab of Photonic Technique for Information, School of Electronic Science and Engineering, Faculty of Electronic and Information Engineering, Xi’an Jiaotong University, Xi’an 710049, China
*
Authors to whom correspondence should be addressed.
Nanomaterials 2022, 12(20), 3700; https://doi.org/10.3390/nano12203700
Submission received: 14 September 2022 / Revised: 14 October 2022 / Accepted: 18 October 2022 / Published: 21 October 2022

Abstract

:
The acceptor-donor-acceptor (A-D-A) type conjugated organic molecule has been widely applied in the organic optoelectronics field. A total of Nine compounds (19) were designed under the A-D-A framework, with the electron donor benzodithiophene as the core and dicyanomethylene as the acceptor moiety, modifying the benzodithiophene with the phenyl, naphthyl, and difluorinated phenyl groups. The conjugation length can be changed by introducing a thiophene π-conjugated bridge. The geometric structures, electronic structure, excited state properties, aromaticity, and the static- and frequency-dependent second hyperpolarizabilities were investigated by employing high-precision density functional theory (DFT) calculations with an aug-cc-pVDZ basis set. As a result, the three compounds with the longest conjugation length exhibit a smaller energy gap (Egap), larger UV-vis absorption coefficient, and response range, which are the three strongest third-order nonlinear optical (NLO) response properties in this work. This work systematically explored the connection between molecular structure and NLO response, which provides a rational design strategy for high-performance organic NLO materials.

Graphical Abstract

1. Introduction

With the rapid development of nonlinear optics, nonlinear optical (NLO) materials have been widely applied in the field of fiber optics [1], data transformation [2], photonic lasers, and data storage in the area of wireless communication [3]. At present, the design and synthesis of high-performance NLO materials in organic and inorganic systems still meet major challenges [4]. In recent decades, many materials exhibited good optical nonlinearity, including organic materials [5,6,7,8], perovskite quantum dots [9], nanoclusters [10,11,12,13,14], organic-inorganic hybrid materials [1,15,16], etc. Among them, organic compounds that are easy to synthesize, with the advance of lower dielectric coefficients, have attracted wide attention [7].
Compared with traditional inorganic materials, organic chromophore compounds with excellent rapid natural response rates, a high absorption coefficient, tunable energy gaps, and molecular structures are easily adjustable by molecular designing, which is considered promising for NLO materials [17]. Organic-photovoltaic (OPV) materials include electronic donors and acceptors that, based on the donor-acceptor-type (D-A) and donor-π-acceptor-type (D-π-A) molecular structures, have been widely applied in OPV devices and achieved exciting power conversion efficiency. Because of the highly π-conjugated D-A/D-π-A molecular structure, the intramolecular charge transfer (ICT) would be easier and greater for obtaining efficient NLO characteristics [7]. The properties of substituent groups and the conjugated length of the molecules can significantly influence NLO characteristics.
In this work, benzodithiophene was selected as the core, which is generally a key structure of the donor unit in D-A/D-π-A polymer electronic donors in OPV [18,19,20,21,22,23,24]. For the NLO materials, the centrosymmetry is a basic structural feature, and the molecular structure of the D-A-D type has symmetry itself. For organic molecules, because the D-A-D molecular structure has a significant ICT effect, it is easy to achieve better photo response characteristics, and the ICT of molecules can be adjusted by changing the electronic structure of the side chains and the conjugate length of the molecules so that nine molecules were designed with a terminal dicyanomethylene acceptor as the end-group, and thiophene as a conjugate extension, meaning that the electron-deficient terminal groups could be conducive to enhancing the molecule’s absorption coefficient. Among the nine compounds, as displayed in Figure 1, compounds 13, 46, and 79 were designed by extending the molecule’s conjugated length by increasing the thiophene rings. The compounds 13, 46, and 79 modified the benzodithiophene with phenyl, naphthyl, and difluorinated phenyl groups, respectively. Based on this modification strategy, the influence of the conjugated length of the molecule and the molecular electronic structure changes on linear optical and third-order NLO properties were deeply explored for compounds 19. The molecular electronic structure, energy gap, second hyperpolarizability (γ), UV−vis spectra, excited hole-electron character, and anisotropy of the induced current density (AICD) are comprehensively investigated by density functional theory (DFT) calculations. The calculated results indicate that compounds 3, 6, and 9 exhibited marvelous third-order NLO properties. With the compound’s conjugation length increased, the Egap is smaller, and the UV-vis response range is wider, with the linear response intensity and the third-order NLO response being stronger. On account of this, this NLO material design strategy will help to develop new high-performance organic NLO materials.

2. Computational Methods

All the DFT calculations were carried out using the Gaussian 16 package [25], and the Visual Molecular Dynamics (VMD) software was employed to visualize the calculated results [26]. All the geometries were optimized under the B3LYP-D3(BJ)/def2SVP level [27,28]. After geometry optimization, we further confirmed no imaginary frequency by vibrational frequency calculation. The excited states of the nine geometries were calculated by the time-dependent density functional theory (TD-DFT) using the CAM-B3LYP method with a def2SVP basis set [29]. In order to obtain proficiency and accuracy, static and frequency-dependent third-order NLO hyperpolarizabilities (γ) were calculated using the CAM-B3LYP method by using the high accuracy aug-cc-pVDZ basis set [30,31,32,33,34,35,36,37,38,39]. The molecule’s electronic structure, hyperpolarizability, excited state properties, and aromaticity were analyzed by Multiwfn 3.8(dev) code [40].
The second hyperpolarizability density ρ x x x ( 3 ) of nine molecules is contrasted between the dipole moment μ(F) after Taylor expansion and electron density, as follows [7,11].
E F = E 0 μ 0 F 1 2 α F 2 1 6 β F 3 1 24 β F 4
μ F = E F = μ 0 + α F + 1 2 β F 2 + 1 2 β F 3 +
ρ r , F = ρ 0 r + ρ 1 r F + 1 2 ρ 2 r F 2 + 1 6 ρ 3 r F 3 +
where μ0 is the permanent dipole moment and F is the uniform external electric field.
γ x x x = ρ x x x 3 r x d x
ρ x x x 3 = ρ 2 F x 2 ρ F x + 2 ρ F x ρ 2 F x 2 ( F x ) 3
The i components of γ are defined as
γ i = 1 15 j γ i j j i + γ i j i j + γ i i j j i , j = x , y , z
The total magnitude of γ is measured as
γ t o t = γ x 2 + γ y 2 + γ z 2

3. Results and Discussion

All the geometries were optimized in the gas phase under the B3LYP-D3(BJ)/def2SVP level. As shown in Figure 2, Figures S1 and S2, all the molecules are rigid and planar along the main skeleton. The HOMOs (highest occupied molecular orbitals), LUMOs (lowest unoccupied molecular orbitals), and the energy gaps (Egap) of the nine molecules are calculated as displayed in Figure 2. The HOMO energy of compound 1 was −6.46 eV, with a LUMO energy of −4.62 eV and an Egap of 1.84 eV. Similarly, the HOMO and LUMO energy of compound 2 were −5.62 eV and −4.59 eV, respectively, with an Egap of 1.03 eV. The remaining compounds 3, 4, 5, 6, 7, 8, and 9 had a HOMO energy of −5.16, −6.36, −5.58, −5.14, −6.70, −5.77, and −5.27 eV, respectively. While their LUMO energy was found to be −4.61, −4.57, −4.56, −4.59, −4.87, −4.75, and −4.73 eV, respectively, with an Egap of 0.55, 1.79, 1.02, 0.55, 1.83, 1.02, and 0.54 eV, respectively. The order of Egap for the compounds was: compound 1 < compound 2 < compound 3, compound 4 < compound 5 < compound 6, compound 7 < compound 8 < compound 9. It is clearly shown that with an increase in the molecule’s conjugated length, the Egap decreased. The TDOS (total density of states) plots were calculated to show the population evaluation of each orbital, which reflects the number of molecular orbitals (MOs) in the unit energy interval at the corresponding energy.
There is an important connection between the molecule charge distribution and linear molecular polarizability (α) [41]. The molecular polarizability can be influenced by the applied external electric field. The linear isotropic (αiso) and anisotropic (αaniso) polarizabilities, as well as the individual tensor components of all the compounds, were calculated at the CAM-B3LYP/aug-cc-PVDZ level, as given in Table 1. The αxx shows larger values among the tensor components, which suggests the polarizability mainly occurred along the x-axis. Compounds 3, 6, and 9 show larger polarizability values of 7979.27, 8018.49, and 8023.75 a.u., respectively. They are much larger than the other remaining compounds. The compounds exhibited a trend of increasing linear polarizabilities, which corresponds well with the trend of increasing the conjugated length of the molecule, as 1 < 2 < 3, 4 < 5 < 6, 7 < 8 < 9. As displayed in Figure 3, the αiso and αaniso increased with the increments in π-conjugation length. For all the compounds, compound 1 exhibits the smallest value of αiso and αaniso, which were 556.12 and 638.61 a.u., respectively. However, compound 9 exhibited the largest αiso and αaniso of 3035.44 and 7489.37 a.u., respectively. The large difference between the αiso and αaniso values for all the compounds suggests that the direction of the response of these materials to external electric fields is significant selectively.
The second-order hyperpolarizability (γ) is considered to be an important index to measure nonlinear optical coefficients. The γ values were calculated under the CAM-B3LYP/aug-cc-PVDZ level, which is shown in Figure 4a and is summarized in Table S1. The results of static third-order NLO coefficients (γtot) of these nine compounds show that when the backbone has the same conjugate length, the γtot of those molecules with difluorinated phenyl groups is stronger than that of phenyl- and naphthyl-modified side groups. As shown in Figure 4b, the results suggest that, for those compounds in which the Egap is smaller, the polarizability is larger for their smaller excitation energy. To further investigate the NLO properties of the nine compounds, the second hyperpolarizabilities γ (−2ω; ω, ω, 0) were calculated, as summarized in Table 2. The sequence of the γ (−2ω; ω, ω, 0) values were γ (∞ nm) < γ (1064 nm) < γ (532 nm) for compounds 1, 2, 4, 5, 7, and 8, and for compounds 3, 6 and 9, the γ (−2ω; ω, ω, 0) values were γ (532 nm) < γ (∞ nm) < γ (1064 nm). Which are higher than the reported bis(dicyanomethylene) end-functionalized quinoidal oligothiophenes [7].
The anisotropy of the third-order NLO polarizabilities can be visualized directly, perceived through the unit sphere representation [11,42,43], as shown in Figure 5. The change of the molecular dipole moment can be measured based on the unit sphere representation due to the length and direction of the arrows, representing the strength and direction of the dipole moment, which can show the direction and intensity of the γ. Furthermore, the hyperpolarizability density was calculated using the numerical differentiation method, and the contribution of hyperpolarizability in different regions can be clearly presented. Based on the calculated hyperpolarizability density results, as shown in Figure 6, the gray and blue-black regions denote positive and negative spatial contributions, with superposition effects for all the compounds, and the hyperpolarizability density can clearly show the magnitude of hyperpolarizability and the main contribution regions. All the distributions of X p x x x ( 3 ) ( r ) are centrosymmetric for the nine centrosymmetric molecules [42,43]. Moreover, compared with Figure 5 and Figure 6, it can be seen that the length of the π-conjugate has a significant positive effect on the third-order NLO response.
The molecular surface electrostatic potentials of the compounds were further calculated to understand the characteristics of charge distribution on the molecule surface. As shown in Figure 7a, the red and blue colors represent the relative negative and positive potential surfaces, respectively, and the terminal dicyanomethylene groups show a maximum negative potential. The maximum positive potential region appears on the central benzodithiophene core or π-bridge (thiophene moiety). In order to further confirm the molecule aromaticity, the anisotropy of the induced current density (AICD) plot of the compounds was further calculated, as shown in Figure 7b. The results indicated that all the compounds would generate ring current under a magnetic field, π-electrons delocalized within the whole molecule, meaning they are aromatic.
As displayed in Figure 8, the absorption spectra of the compounds were obtained by TDDFT calculations at the CAM-B3LYP/def2SVP level. We first obtained the excited state energy and oscillator strength, and the excited state energy was converted into nm units. They were then expanded into curves and added together, and the spectra were broadened by a Gaussian function with a full width at half maximum (FWHM) of 0.6667 eV. The transition energy E (eV), oscillator strength (f), maximum absorption wavelength λmax (nm), and molecular orbit transition contributions of the investigated compounds (19) are presented in Table 3. The calculated λmax of the compounds (19) is attributed to the electron excited from the ground state to the first excited state. The simulated compounds absorption spectra ranged from visible light and extended into the infrared region. The simulated compounds spectra scale shows a redshift tendency, and the oscillator strength is simultaneously enhanced with the increasing π-conjugate length of the molecule. Meanwhile, in order to further confirm the reliability of the simulated spectra, all the molecules have been simulated using the wB97XD method to calculate the UV-Vis spectra, which was in good agreement with CAM-B3LYP methods, as shown in Figure S3.
In order to deeply understand the charge transfer characteristics of the molecules in their excited states, we calculated the transition density matrix (TDM) plots, as shown in Figure 9. The TDM plots can be more intuitive in reflecting the information of the ICT as well as the specific contributions of the electrons, holes, and electron-hole overlaps. Here, the red frame region represents the central benzodithiophene moiety, and the green frame region represents the thiophene, phenyl, naphthyl, and difluorinated phenyl or C(CN)2 groups. Based on the calculated TDM results, all the molecules exhibited significant ICT behavior. Except for compounds 1, 3, and 5, the electron-hole pairs of other compounds were mainly localized along with the diagonal element. The results suggest the electron transition involves the whole molecule, which is beneficial for intramolecular charge transfer. Compounds 1, 3, and 5 with the structure of the benzodithiophene are substituted with the phenyl, naphthyl, and difluorinated phenyl groups, and the benzodithiophene is directly combined with the dicyanomethylene end-group. For the small, conjugated distance, for compounds 1, 3, and 5, the twist structure between benzodithiophene and the side substituents (phenyl, naphthyl, and difluorinated phenyl groups) had lost the delocalized electric to a great degree so that the structure-property of compounds 1, 3, and 5 further enhanced the electron-hole transition coherence and did not benefit from ICT, as displayed in Figure 9a,d,g. Among all the compounds, compounds 3, 6, and 9 show the most remarkable ICT effects, followed by the largest molecule conjugation length and smaller Egap. Moreover, the off-diagonal part of the TDM also exhibited the probability of electron transformation, which contributed to the molecule’s charge transfer.

4. Conclusions

In this contribution, nine molecules based on benzodithiophene as the core- and dicyanomethylene as the end-group were designed, with their molecule conjugation length changed by increasing the number of thiophene moieties. The molecular electronic structure, TDOS, molecular surface electrostatic potential, anisotropy of the AICD, UV-vis absorption spectra, and TDMs were calculated by DFT and TDDFT to deeply investigate the compounds. The calculated results demonstrate that compounds 3, 6, and 9, with long conjugated lengths, exhibited high third-order NLO responses for their smallest Egap and large UV-vis absorption coefficients. The TDM results also proved that an increase in conjugated length could promote ICT. Therefore, the enhancement of nonlinear optical response properties is mainly caused by an increase in molecule conjugated length. At the same time, the electronic structure of the side chain groups also affects NLO, and modifying those groups with strong electronegativity will enhance ICT and further enhance the NLO. Therefore, to design high-performance NLO molecular materials with an A-D-A structure, researchers should consider increasing the conjugated length of the molecule and introducing electron-withdrawing groups to the donor or acceptor units. The object of our research was to abstract molecules into isolated systems, which has a certain reference significance for the performance of the actual molecules. This work provides a strategy to design and synthesize new organic third-order NLO materials with highly responsive characteristics.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano12203700/s1; Figure S1: The front view of the optimized molecules; Figure S2: The side view of the optimized molecules; Table S1: Static third-order NLO coefficients (γ) of the nine molecules; Figure S3: The UV-vis spectra of the nine compounds calculated under the wB97XD/def2-SVP level.

Author Contributions

Investigation, P.G.; Conceptualization, J.T. and Z.L.; methodology, P.G. and J.T.; formal analysis, L.A.; data curation, Z.L. and L.A.; writing—original draft preparation, P.G.; writing—review and editing, J.T. and Z.L.; visualization, X.L.; supervision, J.L.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (62064006) Innovation Fund Project of Colleges and Universities in Gansu (2021B-063) and the Open Fund of Gansu Province Organic Semiconductor Materials and Technology Research Center (Lanzhou Jiaotong University) (GOTC2022-02).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yi, J.; Du, L.; Li, J.; Yang, L.L.; Hu, L.Y.; Huang, S.H.; Dong, Y.C.; Miao, L.L.; Wen, S.C.; Mochalin, V.N.; et al. Unleashing the potential of Ti2CTx MXene as a pulse modulator for mid-infrared fiber lasers. 2D Mater. 2019, 6, 045038–045047. [Google Scholar] [CrossRef]
  2. Pattipaka, S.; Joseph, A.; Bharti, G.P.; Raju, K.C.J.; Khare, A.; Pamu, D. Thickness-dependent microwave dielectric and nonlinear optical properties of Bi0.5Na0.5TiO3 thin films. Appl. Surf. Sci. 2019, 488, 391–403. [Google Scholar] [CrossRef]
  3. Ropagnol, X.; Kovacs, Z.; Gilicze, B.; Zhuldybina, M.; Blanchard, F.; Garcia-Rosas, C.M.; Szatmari, S.; Foldes, I.B.; Ozaki, T. Intense sub-terahertz radiation from wide-bandgap semiconductor based large-aperture photoconductive antennas pumped by UV lasers. New J. Phys. 2019, 21, 113042–113050. [Google Scholar] [CrossRef] [Green Version]
  4. Kuis, R.; Gougousi, T.; Basaldua, I.; Burkins, P.; Kropp, J.A.; Johnson, A.M. Engineering of large third-order nonlinearities in atomic layer deposition grown nitrogen-enriched TiO2. ACS Photonics 2019, 6, 2966–2973. [Google Scholar] [CrossRef]
  5. Uthayakumar, M.; Jeyakumari, A.P.; Dhandapani, A.; Shinde, V.; Arivanandhan, M. Synthesis, experimental and computational spectroscopic investigations of third-order nonlinear optical material (E)-N′-(benzo[d][1,3]dioxol-5-ylmethylene)benzo-hydrazide. J. Phys. D. Appl. Phys. 2019, 52, 395102–395113. [Google Scholar] [CrossRef]
  6. Dong, D.Y.; Li, Q.Z.; Hou, W.L.; Zhang, H.Q. Synthesis, nonlinear optical, magnetic and electrical properties of ultra-stable open-shell pancake bonding linked perylene diimide anion radicals π-oligomer. J. Mol. Struct. 2020, 1199, 127002–127007. [Google Scholar] [CrossRef]
  7. Bibi, A.; Muhammad, S.; UrRehman, S.; Bibi, S.; Bashir, S.; Ayub, K.; Adnan, M.; Khalid, M. Chemically modified quinoidal oligothiophenes for enhanced linear and third-order nonlinear optical properties. ACS Omega 2021, 6, 24602–24613. [Google Scholar] [CrossRef]
  8. Chang, H.J.; Bondar, M.V.; Liu, T.; Liu, X.; Singh, S.; Belfield, K.D.; Sheely, A.; Masunov, A.E.; Hagan, D.J.; van Stryland, E.W. Electronic nature of neutral and charged two-photon absorbing squaraines for fluorescence bioimaging application. ACS Omega 2019, 4, 14669–14679. [Google Scholar] [CrossRef]
  9. Ketavath, R.; Kattur, N.K.; Ghugal, S.G.; Kolli, H.K.; Swetha, T.; Soma, V.R.; Murali, B. Deciphering the ultrafast nonlinear optical properties and dynamics of pristine and Ni-doped CsPbBr3 colloidal two-dimensional nanocrystals. J. Phys. Chem. Lett. 2019, 10, 5577–5584. [Google Scholar] [CrossRef]
  10. Mbarak, H.; Kodeary, A.K.; Hamidi, S.M.; Mohajarani, E.; Zaatar, Y. Control of nonlinear refractive index of AuNPs doped with nematic liquid crystal under external electric field. Optik 2019, 198, 163299–163304. [Google Scholar] [CrossRef]
  11. Wang, Z.R.; Zhao, G.X.; Yan, W.Y.; Wu, K.C.; Wang, F.; Li, Q.H.; Zhang, J. Tin metal cluster compounds as new third-order nonlinear optical materials by computational study. J. Phys. Chem. Lett. 2021, 12, 7537–7544. [Google Scholar] [CrossRef]
  12. Bhattacharya, D.; Ghoshal, D.; Bhattacharya, S.; Mondal, D.; Paul, B.K.; Bose, N.; Datta, P.K.; Das, S.; Basu, M. Third-order optical nonlinearity of the CuCo0.5Ti0.5O2 nanostructure under 120 fs laser irradiation. Appl. Opt. 2019, 58, 9163–9171. [Google Scholar] [CrossRef]
  13. Sujatha, R.A.; Flower, N.A.L.; Vinitha, G.; Sharath, R.A.; Rahulan, K.M. Structural and non-linear optical response of Er3+ doped SrMoO4 nanostructures. Appl. Surf. Sci. 2019, 490, 260–265. [Google Scholar] [CrossRef]
  14. Jagadale, T.C.; Murali, D.S.; Chu, S.W. Nonlinear absorption and scattering of a single plasmonic nanostructure characterized by x-scan technique. Beilstein J. Nanotechnol. 2019, 10, 2182–2191. [Google Scholar] [CrossRef] [Green Version]
  15. Cocca, L.H.Z.; Sciuti, L.F.; Menezes, L.B.; Kohler, M.H.; Bevilacqua, A.C.; Piquini, P.C.; Iglesias, B.A.; de Boni, L. Excited-state investigations of meso-mono-substituted-(amino-ferrocenyl) porphyrins: Experimental and theoretical approaches. J. Photoch. Photobio. A 2019, 384, 112048–112058. [Google Scholar] [CrossRef]
  16. Liu, Z.W.; Gan, F.; Dong, N.N.; Zhang, B.; Wang, J.; Chen, Y. Fabrication and nonlinear optical characterization of fluorinated zinc phthalocyanine covalently modified black phosphorus/PMMA films using the nanosecond Z-scan technique. J. Mater. Chem. C 2019, 7, 10789–10794. [Google Scholar] [CrossRef]
  17. Shelton, A.H.; Valandro, S.R.; Price, R.S.; Dubinina, G.G.; Abboud, K.A.; Wicks, G.; Rebane, A.; Younus, M.; Schanze, K.S. Stereochemical effects on platinum acetylide two-photon chromophores. J. Phys. Chem. A 2019, 123, 9382–9393. [Google Scholar] [CrossRef]
  18. Tong, J.; Huang, Y.; Liu, W.; Shou, M.; An, L.; Jiang, X.; Guo, P.; Han, Y.; Liang, Z.; Li, J.; et al. Optimized molecular aggregation and photophysical process synergistically promoted photovoltaic performance in low regularity benzo[c][1,2,5]thiadiazole-based medium bandgap copolymers via modulating π bridges. J. Mater. Chem. C 2022. [Google Scholar] [CrossRef]
  19. Xue, L.; Liu, X.; Wang, Q.; Yang, M.; Du, S.; Yang, C.; Tong, J.; Xia, Y.; Li, J. Improved Performance of Organic Solar Cells by Utilizing Green Non-Halogen Additive to Modulate Active-Layer Morphology. Energy Technol. 2022, 10, 2200504. [Google Scholar] [CrossRef]
  20. Ren, Y.; Liu, X.; Li, H.; Qin, J.; Du, S.; Lu, X.; Tong, J.; Yang, C.; Li, J. Utilizing non-conjugated small-molecular tetrasodium iminodisuccinateas electron transport layer enabled improving efficiency of organic solar cells. Opt. Mater. 2022, 129, 112520. [Google Scholar] [CrossRef]
  21. Liu, X.; Liang, Z.; Du, S.; Niu, X.; Tong, J.; Yang, C.; Lu, X.; Bao, X.; Yan, L.; Li, J.; et al. Two compatible acceptors as an alloy model with a halogen-free solvent for efficient ternary polymer solar cells. ACS Appl. Mater. Interfaces 2022, 14, 9386–9397. [Google Scholar] [CrossRef]
  22. Yang, C.; Zhang, C.; Chen, C.; Ren, Y.; Shen, H.; Tong, J.; Du, S.; Xia, Y.; Li, J. A New Alcohol-Soluble Polymer PFN-ID as Cathode Interlayer to Optimize Performance of Conventional Polymer Solar Cells by Increasing Electron Mobility. Energy Technol. 2022, 10, 2200199. [Google Scholar] [CrossRef]
  23. Liang, Z.; Tong, J.; Li, H.; Wang, Y.; Wang, N.; Li, J.; Yang, C.; Xia, Y. The comprehensive utilization of the synergistic effect of fullerene and non-fullerene acceptors to achieve highly efficient polymer solar cells. J. Mater. Chem. A 2019, 7, 15841–15850. [Google Scholar] [CrossRef]
  24. Liang, Z.; Yan, L.; Si, J.; Gong, P.; Li, X.; Liu, D.; Li, J.; Hou, X. Rational design and characterization of symmetry-breaking organic semiconductors in polymer solar cells: A theory insight of the asymmetric advantage. Materials 2021, 14, 6723. [Google Scholar] [CrossRef]
  25. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16 Rev. C.01; Gaussian Inc.: Wallingford, UK, 2019. [Google Scholar]
  26. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  27. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104–154123. [Google Scholar] [CrossRef] [Green Version]
  28. Weigend, F. Accurate Coulomb-fitting basis sets for H to Rn. Phys. Chem. Chem. Phys. 2006, 8, 1057–1065. [Google Scholar] [CrossRef]
  29. Yanai, T.; Tew, D.P.; Handy, N.C. A new hybrid exchange-correlation functional using the Coulomb-attenuating method (CAM-B3LYP). Chem. Phys. Lett. 2004, 393, 51–57. [Google Scholar] [CrossRef] [Green Version]
  30. Hill, J.G.; Peterson, K.A. Gaussian basis sets for use in correlated molecular calculations. XI. Pseudopotential-based and all-electron relativistic basis sets for alkali metal (K-Fr) and alkaline earth (Ca-Ra) elements. J. Chem. Phys. 2017, 147, 244106–244117. [Google Scholar] [CrossRef] [Green Version]
  31. Prascher, B.P.; Woon, D.E.; Peterson, K.A.; Dunning, T.H.; Wilson, A.K. Gaussian basis sets for use in correlated molecular calculations. VII. Valence, core-valence, and scalar relativistic basis sets for Li, Be, Na, and Mg. Theor. Chem. Acc. 2011, 128, 69–82. [Google Scholar] [CrossRef]
  32. Dunning, T.H.; Peterson, K.A.; Wilson, A.K. Gaussian basis sets for use in correlated molecular calculations. X. The atoms aluminum through argon revisited. J. Chem. Phys. 2001, 114, 9244–9253. [Google Scholar] [CrossRef] [Green Version]
  33. Van Mourik, T.; Dunning, T.H. Gaussian basis sets for use in correlated molecular calculations. VIII. Standard and augmented sextuple zeta correlation consistent basis sets for aluminum through argon. Int. J. Quantum. Chem. 2000, 76, 205–221. [Google Scholar] [CrossRef]
  34. Wilson, A.K.; Woon, D.E.; Peterson, K.A.; Dunning, T.H. Gaussian basis sets for use in correlated molecular calculations. IX. The atoms gallium through krypton. J. Chem. Phys. 1999, 110, 7667–7676. [Google Scholar] [CrossRef] [Green Version]
  35. Wilson, A.K.; van Mourik, T.; Dunning, T.H. Gaussian basis sets for use in correlated molecular calculations 6. Sextuple zeta correlation consistent basis sets for boron through neon. J. Mol. Struc. Theochem. 1996, 388, 339–349. [Google Scholar] [CrossRef]
  36. Woon, D.E.; Dunning, T.H. Gaussian-basis sets for use in correlated molecular calculations 5. core-valence basis-sets for boron through neon. J. Chem. Phys. 1995, 103, 4572–4585. [Google Scholar] [CrossRef] [Green Version]
  37. Woon, D.E.; Dunning, T.H. Gaussian-basis sets for use in correlated molecular calculations 4. calculation of static electrical response properties. J. Chem. Phys. 1994, 100, 2975–2988. [Google Scholar] [CrossRef] [Green Version]
  38. Woon, D.E.; Dunning, T.H. Gaussian-basis sets for use in correlated molecular calculations. III. The atoms aluminum through argon. J. Chem. Phys. 1993, 98, 1358–1371. [Google Scholar] [CrossRef] [Green Version]
  39. Dunning, T.H. Gaussian-basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 1989, 90, 1007–1023. [Google Scholar] [CrossRef]
  40. Lu, T.; Chen, F.W. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef]
  41. Shinde, S.S.; Sreenath, M.C.; Chitrambalam, S.; Joe, I.H.; Sekar, N. Non-linear optical properties of disperse blue 354 and disperse blue183 by DFT and Z-scan technique. Polycycl. Aromat. Comp. 2021, 41, 1531–1548. [Google Scholar] [CrossRef]
  42. Liu, Z.; Wang, X.; Lu, T.; Yuan, A.; Yan, X. Potential optical molecular switch: Lithium@cyclo[18]carbon complex transforming between two stable configurations. Carbon 2022, 187, 78–85. [Google Scholar] [CrossRef]
  43. Liu, Z.; Lu, T.; Yuan, A.; Wang, X.; Chen, Q.; Yan, X. Remarkable size effect on photophysical and nonlinear optical properties of all-carboatomic rings, cyclo[18]carbon and its analogues. Chem. Asian J. 2021, 16, 2267–2271. [Google Scholar] [CrossRef]
Figure 1. Chemical structures of the nine compounds (19).
Figure 1. Chemical structures of the nine compounds (19).
Nanomaterials 12 03700 g001
Figure 2. DOS curves of the nine compounds. (ai) are compounds (19).
Figure 2. DOS curves of the nine compounds. (ai) are compounds (19).
Nanomaterials 12 03700 g002
Figure 3. The linear isotropic (αiso) and anisotropic (αaniso) polarizability.
Figure 3. The linear isotropic (αiso) and anisotropic (αaniso) polarizability.
Nanomaterials 12 03700 g003
Figure 4. (a) Static third-order NLO coefficients (γ) of the compounds (19); (b) the relationship between the static third-order NLO polarizabilities and Egap.
Figure 4. (a) Static third-order NLO coefficients (γ) of the compounds (19); (b) the relationship between the static third-order NLO polarizabilities and Egap.
Nanomaterials 12 03700 g004
Figure 5. (ah) the anisotropy of the hyperpolarizability of the nine compounds (19); the blue and red represent arrow difference in length.
Figure 5. (ah) the anisotropy of the hyperpolarizability of the nine compounds (19); the blue and red represent arrow difference in length.
Nanomaterials 12 03700 g005
Figure 6. (ai) are the isosurfaces of X p x x x ( 3 ) ( r ) (isovalue = 1100) of the nine compounds (19).
Figure 6. (ai) are the isosurfaces of X p x x x ( 3 ) ( r ) (isovalue = 1100) of the nine compounds (19).
Nanomaterials 12 03700 g006
Figure 7. (a) The molecular surface electrostatic potentials of all the compounds (19), red and blue, represent the negative and positive potential, respectively; (b) AICD plots (isovalue: 0.03 a.u.) of the compounds (19).
Figure 7. (a) The molecular surface electrostatic potentials of all the compounds (19), red and blue, represent the negative and positive potential, respectively; (b) AICD plots (isovalue: 0.03 a.u.) of the compounds (19).
Nanomaterials 12 03700 g007
Figure 8. The UV−vis spectra of the nine compounds.
Figure 8. The UV−vis spectra of the nine compounds.
Nanomaterials 12 03700 g008
Figure 9. (ai) TDM plots of the compounds (19).
Figure 9. (ai) TDM plots of the compounds (19).
Nanomaterials 12 03700 g009
Table 1. Calculated α (a.u.) and the individual components.
Table 1. Calculated α (a.u.) and the individual components.
αComp 1 (a.u.)Comp 2 (a.u.)Comp 3 (a.u.)Comp 4 (a.u.)Comp 5 (a.u.)Comp 6 (a.u.)Comp 7 (a.u.)Comp 8 (a.u.)Comp 9 (a.u.)
αxx952.082996.137979.27810.353028.708018.49575.883007.868023.75
αxy−34.435.86−69.60180.20292.75190.93186.8592.3415.61
αyy489.14598.37726.34921.301834.517928.69865.15601.55726.26
αxz0.000020−0.00680.02210−0.0000300.004020.036−0.000020−0.00940.023
αyz−0.000030−0.0059−0.00531−0.0000200.007150.0100.000010−0.017−0.0052
αzz227.12290.63354.587308.309371.897437.18231.75293.634356.31
αiso556.121295.053020.06679.9861411.73128.12557.591301.023035.40
αaniso638.612565.537446.75646.1232510.117355.33637.492579.087489.40
Table 2. Calculated static and dynamic second hyperpolarizability (γ).
Table 2. Calculated static and dynamic second hyperpolarizability (γ).
Incident LightSecond Hyperpolarizability (γ)
Comp 1 (a.u.)Comp 2 (a.u.)Comp 3 (a.u.)Comp 4 (a.u.)Comp 5 (a.u.)Comp 6 (a.u.)Comp 7 (a.u.)Comp 8 (a.u.)Comp 9 (a.u.)
∞ (static)1.94 × 1051.85 × 1064.18 × 1073.29 × 1052.40 × 1064.22 × 1072.45 × 1051.84 × 1064.02 × 107
532 nm1.18 × 10106.70 × 1071.70 × 1071.02 × 1099.42 × 1072.59 × 1071.66 × 1091.68 × 1082.77 × 107
1064 nm3.01 × 1073.66 × 1075.81 × 1081.01 × 1073.78 × 1072.73 × 1093.13 × 1073.83 × 1071.41 × 109
Table 3. Excitation energy (E), oscillator strengths (f), wavelengths (λmax), and MO contributions of the compounds.
Table 3. Excitation energy (E), oscillator strengths (f), wavelengths (λmax), and MO contributions of the compounds.
Compoundλmax (nm)E(eV)fMO Contributions a
1529.982.341.35H → L (98.07%)
2801.551.542.56H → L (97.40%)
31187.251.043.34H → L (98.03%)
4544.822.271.12H → L (98.00%)
5806.191.542.51H → L (97.50%)
61215.761.023.29H → L (97.90%)
7527.572.351.39H → L (98.18%)
8802.121.552.58H → L (97.50%)
91218.151.013.32H → L (98.00%)
a H = HOMO, L = LUMO.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gong, P.; An, L.; Tong, J.; Liu, X.; Liang, Z.; Li, J. Design of A-D-A-Type Organic Third-Order Nonlinear Optical Materials Based on Benzodithiophene: A DFT Study. Nanomaterials 2022, 12, 3700. https://doi.org/10.3390/nano12203700

AMA Style

Gong P, An L, Tong J, Liu X, Liang Z, Li J. Design of A-D-A-Type Organic Third-Order Nonlinear Optical Materials Based on Benzodithiophene: A DFT Study. Nanomaterials. 2022; 12(20):3700. https://doi.org/10.3390/nano12203700

Chicago/Turabian Style

Gong, Pingping, Lili An, Junfeng Tong, Xinpeng Liu, Zezhou Liang, and Jianfeng Li. 2022. "Design of A-D-A-Type Organic Third-Order Nonlinear Optical Materials Based on Benzodithiophene: A DFT Study" Nanomaterials 12, no. 20: 3700. https://doi.org/10.3390/nano12203700

APA Style

Gong, P., An, L., Tong, J., Liu, X., Liang, Z., & Li, J. (2022). Design of A-D-A-Type Organic Third-Order Nonlinear Optical Materials Based on Benzodithiophene: A DFT Study. Nanomaterials, 12(20), 3700. https://doi.org/10.3390/nano12203700

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop