<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">CC</journal-id><journal-title-group><journal-title>Computational Chemistry</journal-title></journal-title-group><issn pub-type="epub">2332-5968</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/cc.2020.84005</article-id><article-id pub-id-type="publisher-id">CC-103869</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Chemistry&amp;Materials Science</subject></subj-group></article-categories><title-group><article-title>
 
 
  Electronic Absorption Spectra and Third-Order Nonlinear Optical Property of Dinaphtho[2,3-b:2’,3’-d]Thiophene-5,7,12,13- Tetraone (DNTTRA) and Its Phenyldiazenyl Derivatives: DFT Calculations
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ziran</surname><given-names>Chen</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Yujin</surname><given-names>Zhang</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Zhanrong</surname><given-names>He</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Yonghua</surname><given-names>Guan</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Yuan</surname><given-names>Li</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Hongping</surname><given-names>Li</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff3"><addr-line>Institute for Energy Research, Jiangsu University, Zhenjiang, China</addr-line></aff><aff id="aff2"><addr-line>Department of Electronic and Electrical Engineering, Sichuan Vocational and Technical College, Suining, China</addr-line></aff><aff id="aff1"><addr-line>Department of Architecture and Environment Engineering, Sichuan Vocational and Technical College, Suining, China</addr-line></aff><pub-date pub-type="epub"><day>30</day><month>10</month><year>2020</year></pub-date><volume>08</volume><issue>04</issue><fpage>43</fpage><lpage>60</lpage><history><date date-type="received"><day>10,</day>	<month>October</month>	<year>2020</year></date><date date-type="rev-recd"><day>27,</day>	<month>October</month>	<year>2020</year>	</date><date date-type="accepted"><day>30,</day>	<month>October</month>	<year>2020</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Third-order nonlinear optical (NLO) materials have broad application prospects in high-density data storage, optical computer, modern laser technology, and other high-tech industries. The structures and frequencies of Dinaphtho[2,3-b:2’,3’-d]thiophene-5,7,12,13-tetraone (DNTTRA) and its 36 derivatives containing azobenzene were calculated by using density functional theory B3LYP and M06-2X methods at 6-311++g(d, p) level, respectively. Besides, the atomic charges of natural bond orbitals (NBO) were analyzed. The frontier orbitals and electron absorption spectra of A-G5 molecule were calculated by TD-DFT (TD-B3LYP/6-311++g(d, p) and TD-M06-2X/6-311++g(d, p)). The NLO properties were calculated by effective finite field FF method and self-compiled program. The results show that 36 molecules of these six series are D-π-A-π-D structures. The third-order NLO coefficients γ (second-order hyperpolarizability) of the D series molecules are the largest among the six series, reaching 10
  <sup>7</sup> atomic units (10
  <sup>-33</sup> esu) of order of magnitude, showing good third-order NLO properties. Last, the third-order NLO properties of the azobenzene ring can be improved by introducing strong electron donor groups (e.g. -N(CH
  <sub>3</sub>)
  <sub>2</sub> or -NHCH
  <sub>3</sub>) in the azobenzene ring, so that the third-order NLO materials with good performance can be obtained.
 
</p></abstract><kwd-group><kwd>Azobenzene</kwd><kwd> Dinaphtho[2</kwd><kwd>3-b:2’</kwd><kwd>3’-d]Thiophene-5</kwd><kwd>7</kwd><kwd>12</kwd><kwd>13-Tetraone</kwd><kwd> Density Functional Theory</kwd><kwd> Electronic Absorption Spectra</kwd><kwd> Third-Order Nonlinear Optical Property</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Third-order nonlinear optical (NLO) materials have the characteristics of fast nonlinear optical response, large nonlinear polarizability, wide response band, high optical damage domain value, excellent flexibility, good chemical and thermal stability, and are easy to modify molecular structure. This kind of material has broad application prospects in high-density data storage, optical computer, modern laser technology, and other high-tech industries [<xref ref-type="bibr" rid="scirp.103869-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref4">4</xref>]. Therefore, many chemists are attracted to design and synthesize new compounds with strong nonlinear optical effects. A large number of research results show that the third-order nonlinearity of materials is closely related to the type of molecular structure [<xref ref-type="bibr" rid="scirp.103869-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref6">6</xref>]. For example, the linear and nonlinear polarizabilities of polyene and polymethine dyes have been shown to be highly correlated with the bond length alternation between adjacent C-C bonds in the conjugated chain. At present, the research focus of nonlinear optical materials is mainly focused on the construction of organic molecules containing electron donor (D), electron acceptor (A) and π-conjugation [<xref ref-type="bibr" rid="scirp.103869-ref7">7</xref>]. For example, D-π-A, A-π-D-π-A, D-π-A-π-D and other organic molecules [<xref ref-type="bibr" rid="scirp.103869-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref9">9</xref>], focus on increasing the third-order NLO coefficients (second-order hyperpolarizability), so as to improve the optical properties of organic materials [<xref ref-type="bibr" rid="scirp.103869-ref10">10</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref11">11</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref12">12</xref>].</p><p>Quinone heterocyclic compounds are a kind of conjugated organic functional materials [<xref ref-type="bibr" rid="scirp.103869-ref13">13</xref>]. The intramolecular charge transfer compound synthesized by using it as the electron donor has become the preferred material for organic superconductors [<xref ref-type="bibr" rid="scirp.103869-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref15">15</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref16">16</xref>]. It is widely used in organic semiconductors, airport devices, organic light-emitting diodes, photovoltaic cells and other fields. At present, anthraquinone and naphthoquinone derivatives have been widely studied in the third-order nonlinear optical materials. The nonlinear optical effect of this kind of molecule mainly depends on its molecular configuration. The third-order hyperpolarizability of this kind of molecule increases obviously with the growth of its conjugated chain. At the same time, the sulfur-containing heterocycle in the conjugated heterocycle structure contributes significantly to its third-order hyperpolarizability due to its high hole mobility values and air stability, so it has attracted the attention of chemical researchers [<xref ref-type="bibr" rid="scirp.103869-ref17">17</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref18">18</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref19">19</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref21">21</xref>]. For example, Dibenzo [b,i] thiazene-5,7,12,14-tetraketone was synthesized from 2,3-dichloro-1,4-naphthoquinone and sodium sulfide by Gao’s research group [<xref ref-type="bibr" rid="scirp.103869-ref13">13</xref>].</p><p>Azo aromatic organic compounds are a kind of nonlinear optical materials with good properties. Because of its photosensitivity and photoisomerism, it has a great potential application in nonlinear fields such as optical information storage [<xref ref-type="bibr" rid="scirp.103869-ref22">22</xref>] - [<xref ref-type="bibr" rid="scirp.103869-ref27">27</xref>]. Many researchers have studied its third-order nonlinear optical properties. For example, Guo et al. [<xref ref-type="bibr" rid="scirp.103869-ref28">28</xref>] study the third-order nonlinearity of materials by increasing π electron delocalization; Kang et al. [<xref ref-type="bibr" rid="scirp.103869-ref29">29</xref>] study the influence of different substituents on the nonlinear optical properties of azobenzene derivatives; Wang et al. [<xref ref-type="bibr" rid="scirp.103869-ref30">30</xref>] report that the introduction of other organic molecules on azobenzene molecules will enhance the nonlinear optical properties of azobenzene.</p><p>In recent years, there have been experimental and theoretical reports on the structure and NLO properties of quinone heterocycles containing diazobenzene [<xref ref-type="bibr" rid="scirp.103869-ref31">31</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref32">32</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref33">33</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref34">34</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref35">35</xref>]. DNTTRA has good rigid conjugate plane and good delocalization. The large π conjugated parent structure would interact with the azo group after conjugated bridging the electron donor of phenyl, which can further enhance the delocalization of the electron, and thus can be designed as an organic electron transport material. These structural molecules belong to D-π-A-π-D type organic molecules, and their chemical modification can improve the nonlinear optical properties of the molecules. Therefore, in this work, we calculated the molecular structure and third-order NLO of DNTTRA and 36 derivatives as shown in Figure1 and FigureS1. The effects of azobenzene on the third-order NLO properties of DNTTRA molecules at 3.9 sites, 2.9 sites, 2.10 sites, 1.8 sites, 1.11 sites and 4.8 sites were also discussed. Secondly, we further studied the introduction of -OH, -OCH<sub>3</sub>, -NH<sub>2</sub>, -NHCH<sub>3</sub>, -N(CH<sub>3</sub>)<sub>2</sub> into the para position of azobenzene ring. It provides a theoretical reference for the further design and synthesis of the third-order nonlinear optical materials containing DNTTRA with excellent properties.</p></sec><sec id="s2"><title>2. Computational Details</title><p>The finite field (FF) method is an effective method for calculating the NLO properties of molecules [<xref ref-type="bibr" rid="scirp.103869-ref35">35</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref36">36</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref37">37</xref>]. In this paper, the structures of A-G5 molecule were optimized by using Gaussian 16 program at B3LYP/6-311++g(d, p) and M06-2X/6-311++g(d, p) computational level, which has been widely used in the NLO predictions [<xref ref-type="bibr" rid="scirp.103869-ref38">38</xref>] [<xref ref-type="bibr" rid="scirp.103869-ref39">39</xref>]. Based on the optimized structures, the electronic absorption spectra of these molecules were calculated by TD-M06-2X/6-311++g (d, p) method and the first superpolarizability (β<sub>&#181;</sub>) in the dipole moment direction were calculated.</p><p>Then, the first hyperpolarizabilities (β<sub>&#181;</sub>) of 37 molecules under different applied electric fields (0.001 - 0.004 a.u.) are calculated by using the finite field method and self-made program. Compared with the first hyperpolarizability (β<sub>&#181;</sub>) in the direction of dipole moment when there is no applied electric field, the value of applied electric field corresponding to the closest value of the two is found. Then, under the external electric field, the components of the third-order nonlinear optical properties obtained by calculation are the credible values of the calculated molecules. Finally, the average values of the third-order NLO properties are calculated by using the following formulas (γ) [<xref ref-type="bibr" rid="scirp.103869-ref35">35</xref>]:</p><p>γ = ( γ xxxx + γ yyyy + γ zzzz + 2 γ xxyy + 2 γ xxzz + 2 γ yyzz ) / 5 (1)</p><p>where the γ<sub>xxxx</sub>, γ<sub>yyyy</sub>, γ<sub>zzzz</sub>, γ<sub>xxyy</sub>, γ<sub>xxzz</sub>, γ<sub>yyzz</sub> are the fourth-order tensor components. The center of thiophene ring is chosen as the coordinate origin when calculating the third-order NLO properties. The XY plane is a conjugate plane, the Z axis is perpendicular to the molecular plane, and the dipole moment of the molecule is roughly the same as the X axis.</p></sec><sec id="s3"><title>3. Results and Discussion</title><sec id="s3_1"><title>3.1. Structures</title><p>The geometric optimization and frequency calculation of A-G5 molecule were carried out at the theoretical level of B3LYP/6-311++g (d, p) and M06-2X/6-311++g (d, p). Among the 37 molecules, A is C2 point group and the rest are C1 point group. Among the 37 molecules calculated, the length of C-C single bond ranges from 0.143 to 0.150 nm; whereas the single bond length of C-N ranges from 0.137 to 0.145 nm [<xref ref-type="bibr" rid="scirp.103869-ref9">9</xref>]. These structure parameters are similar with our previous works [<xref ref-type="bibr" rid="scirp.103869-ref37">37</xref>].</p><p>Among the 37 molecular structures, the dihedral angle of naphthalene and thiophene is between 2.5˚ and 4.8˚, and the dihedral angle of naphthalene and azobenzene is between 1.4˚ and 1.8˚, shown that the molecule is a quasi-planar molecule. Compared with B, C and D, the molecular planarity is poor. The degree of conjugation in conjugated system is consistent with the degree of overlap of π-electron orbital, that is, the degree of overlap of π-electron orbital decreases when the degree of conjugation decreases, thus the charge transfer in molecule cannot be realized effectively, and the NLO effect of molecule is weakened. It can be predicted that the derivatives of azobenzene at 1.8 sites, 1.11 sites and 4.8 sites of DNTTRA molecule have less NLO properties than those of azobenzene at 2.9 sites, 3.9 sites, and 2.10 sites respectively.</p></sec><sec id="s3_2"><title>3.2. Charge Analysis</title><p>In order to explore the correlation between charge transfer and NLO property of organic NLO materials, The charge of 11 organic molecules, such as B, C, D, E, F, G, D1, D2, D3, D4, D5, were calculated by natural bond orbital (NBO) method. The natural population analysis (NPA) charge values of each part of the molecule were listed in <xref ref-type="table" rid="table1"><xref ref-type="table" rid="table">Table </xref>1</xref>. Among them, I and III are benzene ring or phenyl derivative units, and II are conjugated bridges contained DNTTRA and two azo groups.</p><p><xref ref-type="table" rid="table1"><xref ref-type="table" rid="table">Table </xref>1</xref> shows that all the II structural units of 11 molecules are negatively charged. It is known that the conjugated bridge units of DNTTRA and two azo groups are electron acceptors. Compared with the parent D molecule, when R is -OH, -OCH<sub>3</sub>, -NH<sub>2</sub>, -NHCH<sub>3</sub>, -N(CH<sub>3</sub>)<sub>2</sub> and other groups, the negative charge of the II structure unit increases, but does not change its negative charge characteristics. I and III structural units are positively charged and play the role of electron donor in molecules, which will significantly infulence the third-order NLO properties. Compared with the parent D molecule, when R is -OH, -OCH<sub>3</sub>, -NH<sub>2</sub>, -NHCH<sub>3</sub>, -N(CH<sub>3</sub>)<sub>2</sub> and other groups respectively, the positive charge increases, but does not change its positive charge characteristics. Therefore, these molecules can be regarded as D-π-A-π-D structures. With the enhancement of electron-donating ability of substituents introduced by para-position of benzene ring, the electron delocalization increases, which predicts that the third-order NLO properties of molecules will also be enhanced.</p></sec><sec id="s3_3"><title>3.3. Frontier Orbitals and Electron Absorption Spectra</title><p>The molecular gap value is closely related to the intramolecular charge transfer [<xref ref-type="bibr" rid="scirp.103869-ref37">37</xref>]. In order to better understand the nature of molecular charge transfer, the B3LYP and M06-2X methods were used to calculate the energy gap of A-G5 molecules. The calculation results are shown in <xref ref-type="fig" rid="fig2">Figure 2</xref> and TableS1.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1"><xref ref-type="table" rid="table">Table </xref>1</xref></label><caption><title> NPA Charge distribution of 11 molecules</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Comp.</th><th align="center" valign="middle" >B</th><th align="center" valign="middle" >C</th><th align="center" valign="middle" >D</th><th align="center" valign="middle" >E</th><th align="center" valign="middle" >F</th><th align="center" valign="middle" >G</th><th align="center" valign="middle" >D1</th><th align="center" valign="middle" >D2</th><th align="center" valign="middle" >D3</th><th align="center" valign="middle" >D4</th><th align="center" valign="middle" >D5</th></tr></thead><tr><td align="center" valign="middle" >Q(I)/e</td><td align="center" valign="middle" >0.216</td><td align="center" valign="middle" >0.216</td><td align="center" valign="middle" >0.215</td><td align="center" valign="middle" >0.212</td><td align="center" valign="middle" >0.213</td><td align="center" valign="middle" >0.218</td><td align="center" valign="middle" >0.241</td><td align="center" valign="middle" >0.245</td><td align="center" valign="middle" >0.263</td><td align="center" valign="middle" >0.274</td><td align="center" valign="middle" >0.280</td></tr><tr><td align="center" valign="middle" >Q(II)/e</td><td align="center" valign="middle" >−0.390</td><td align="center" valign="middle" >−0.390</td><td align="center" valign="middle" >−0.393</td><td align="center" valign="middle" >−0.418</td><td align="center" valign="middle" >−0.418</td><td align="center" valign="middle" >−0.375</td><td align="center" valign="middle" >−0.437</td><td align="center" valign="middle" >−0.447</td><td align="center" valign="middle" >−0.483</td><td align="center" valign="middle" >−0.507</td><td align="center" valign="middle" >−0.516</td></tr><tr><td align="center" valign="middle" >Q(III)/e</td><td align="center" valign="middle" >0.216</td><td align="center" valign="middle" >0.216</td><td align="center" valign="middle" >0.215</td><td align="center" valign="middle" >0.212</td><td align="center" valign="middle" >0.211</td><td align="center" valign="middle" >0.218</td><td align="center" valign="middle" >0.241</td><td align="center" valign="middle" >0.246</td><td align="center" valign="middle" >0.264</td><td align="center" valign="middle" >0.275</td><td align="center" valign="middle" >0.280</td></tr></tbody></table></table-wrap><p>It can be seen from <xref ref-type="fig" rid="fig2">Figure 2</xref> that the trend of energy gap calculated by the two methods is exactly the same. The energy gap of a molecule is the largest. B, C, D, E, F and G molecules all decrease in the order of A-G-E-F-B-C-D, that is to say, the D molecular energy gap obtained by introducing azobenzene into the 2.10 sites of A molecule is the minimum. Taking B, C, D, E, F and G as the parent molecules, introducing -OH, -OCH<sub>3</sub>, -NH<sub>2</sub>, -NHCH<sub>3</sub>, -N(CH<sub>3</sub>)<sub>2</sub> and other groups into the azobenzene para position respectively, the change trend of the energy gap values of the six series of derivative molecules is exactly the same, both of which decrease with the enhancement of electron donation ability. Among the six series, D series has the lowest energy gap. It is predicted that the series of molecules obtained by introducing azobenzene into the 2.10 sites of A molecule are most favorable for electronic transition.</p><p>Molecular orbitals play an important role in charge transfer. For the calculated six series B, C, D, E, F and G, the HOMO and LUMO electron orbital density distributions are similar. Therefore, <xref ref-type="fig" rid="fig3">Figure 3</xref> only shows the frontier orbital eigenvalues and frontier molecular orbital diagram of 12 molecules such as A, B, C, D, E, F, G and D1-D5. It can be seen from <xref ref-type="fig" rid="fig3">Figure 3</xref> that among the seven molecules A, B, C, D, E, F and G, although the increase of the eigenvalues of HOMO of D molecule is not the most, the decrease of the eigenvalues of LUMO is the largest, resulting in the minimum of the energy gap value. In the D1-D5 series, the eigenvalues of HOMO and LUMO orbitals increased in turn, but the eigenvalues of HOMO increased more, resulting in the decrease of energy gap value in turn. Based on this, it can be preliminarily predicted that with the enhancement of the molecular electron donation ability, the more obvious the intramolecular charge transfer will be. In addition to the F molecule, the HOMO electron orbital is mainly distributed on the azobenzene chain on both sides of DNTTRA. The lone pair electrons on the O and N atoms on the -OH, -OCH<sub>3</sub>, -NH<sub>2</sub>, -NHCH<sub>3</sub>, -N(CH<sub>3</sub>)<sub>2</sub> groups introduced by the azobenzene pair contribute to the HOMO electron orbital, all of which belong to the bonding π orbitals. The LUMO electron orbital of the molecule is mainly distributed on the DNTTRA</p><p>conjugate surface, and they are π* anti-bonding orbits. When electrons transit from HOMO to LUMO, the electron orbital concentrates from azobenzene side chain to DNTTRA. The results show that 18 molecules of B, C, D, E, F, G series have better intramolecular charge transfer effect.</p><p>The electronic absorption spectra of 12 molecules (A, B, C, D, E, F, G and D1-D5) are calculated by TD-M06-2X method, as shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>, <xref ref-type="fig" rid="fig5">Figure 5</xref> and TableS2. It can be seen from <xref ref-type="fig" rid="fig4">Figure 4</xref> that among the seven molecules A, B, C, D, E, F and G, the wavelength of the lowest energy absorption peak of molecule A is 314 nm. B, C, D, E, F and G, are terminated with azobenzene at different positions on both sides of A molecule, and the conjugation range of the molecule increases, resulting in the red shift of the wavelength of the lowest energy absorption peak of the molecule. Especially, the wavelength of the</p><p>absorption peak of D molecule substituted at 2.10 position is the largest (347 nm), which is nearly 100 nm compared with that of the parent A molecule. It is shown that the introduction of azobenzene at the 2.10 position of the parent a molecule is the most favorable for charge transfer, which is consistent with the conclusion of energy gap value analysis. <xref ref-type="fig" rid="fig5">Figure 5</xref> shows that with the D molecule as the parent, D1, D2, D3, D4 and D5 molecules obtained by introducing different electron donating groups into the para position of azobenzene ring, the wavelength of the lowest energy absorption peak increases in turn (significant red shift) with the enhancement of the electron donating ability of the substituent. This is mainly due to the difference in the nature of electronic transition. The absorption peaks of D1 originate from the HOMO-1→LUMO + 1, HOMO-2→LUMO + 2 and HOMO-4→LUMO transitions of S<sub>0</sub>→S<sub>7</sub>; the</p><p>absorption peaks of D2 originate from the HOMO-1→LUMO + 1, HOMO-3→LUMO + 2 and HOMO-5→LUMO transitions of S<sub>0</sub>→S<sub>5</sub>. The absorption peaks of D4 and D5 come from the mixed electronic transitions of S<sub>0</sub>→S<sub>3</sub>, HOMO-1→LUMO + 1 and HOMO-4→LUMO. To sum up, azobenzene was introduced into the two sides of DNTTRA, the introduction of azobenzene (D molecule) at the 2.10 sites is the most favorable for intramolecular charge transfer. After the azobenzene is terminated to the strong electron group (such as -NHCH<sub>3</sub> or -N(CH<sub>3</sub>)<sub>2</sub>, the electron transition energy is relatively reduced, the lowest energy absorption peak is significantly red shifted, and the electron is easy to be excited. Therefore, it can be predicted that the azobenzene has strong third-order NLO properties.</p></sec><sec id="s3_4"><title>3.4. Third-Order NLO Properties</title><p>Using the finite field FF method, an electric field ranging from 0.001 to 0.004 a.u. which applies in both of the negative and positive X, Y, Z directions (the change range is set to 0.00025 to 0.0005 a.u.), it is explored that the first hyperpolarizability (β<sub>&#181;</sub>) of each molecule under applied electric field is closest to the corresponding value of applied electric field. In addition, the non-zero components and average values of the third-order polarizability of 37 molecules calculated under electric field are shown in <xref ref-type="table" rid="table2"><xref ref-type="table" rid="table">Table </xref>2</xref>.</p><p><xref ref-type="table" rid="table2"><xref ref-type="table" rid="table">Table </xref>2</xref> shows that the γ<sub>xxxx</sub> of the third-order polarizability of 37 molecules is the largest, indicating that the third-order NLO properties of DNTTRA derivatives are mainly due to charge transfer in the X-axis direction. Compared with molecule A, the third-order NLO coefficients (γ) increased significantly when azobenzene was introduced at 3.9 sites, 2.9 sites, 2.10 sites, 1.8 sites, 1.11 sites, and 4.8 sites respectively. The third-order NLO coefficient (γ) of azobenzene is the largest when azobenzene is introduced at 2.10 sites. The third-order NLO coefficients (γ) of azobenzene at 3.9 sites, 2.9 sites and 2.10 sites are much larger than those of azobenzene at 1.8 sites, 1.11 sites and 4.8 sites. The results show that azobenzene at 1.8 sites, 1.11 sites and 4.8 sites is not as good as azobenzene at 3.9 sites, 2.9 sites and 2.10 sites, which is consistent with the structural coplanarity analysis of the preceding molecules. The best substitution of azobenzene for A molecule is at 2.10 sites. Among the six series of substituted derivatives such as B, C, D, E, F and G, the third-order NLO coefficient (γ) of the three series derivatives such as B, C and D is obviously larger than that of the three series derivatives such as E, F and G. The third-order NLO coefficients of D series derivatives are the largest in the six series of substituted derivatives, reaching 10<sup>7</sup> order of magnitude atomic units (10<sup>−</sup><sup>33</sup> esu), showing good third order NLO properties. The third-order NLO coefficients (γ) of each series of derivatives have the same change rule, that is, they increase in turn according to the substitution of OH, -OCH<sub>3</sub>, -NH<sub>2</sub>, -NHCH<sub>3</sub>, -N(CH<sub>3</sub>)<sub>2</sub>. This is mainly due to the induction effect and conjugation effect of the introduced donor electrons. -N(CH<sub>3</sub>)<sub>2</sub> and-NHCH<sub>3</sub> have strong electron-donating induction and conjugation effects, while -OH and -OCH<sub>3</sub> have relatively weak induction and conjugation effects.</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2"><xref ref-type="table" rid="table">Table </xref>2</xref></label><caption><title> The third-order NLO properties of 37 molecules (in 10 <sup>5</sup> a .u.)</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Comp.</th><th align="center" valign="middle" >γ<sub>xxxx</sub></th><th align="center" valign="middle" >γ<sub>xxyy</sub></th><th align="center" valign="middle" >γ<sub>xxzz</sub></th><th align="center" valign="middle" >γ<sub>yyyy</sub></th><th align="center" valign="middle" >γ<sub>yyzz</sub></th><th align="center" valign="middle" >γ<sub>zzzz</sub></th><th align="center" valign="middle" >γ</th></tr></thead><tr><td align="center" valign="middle" >A</td><td align="center" valign="middle" >3.39</td><td align="center" valign="middle" >0.30</td><td align="center" valign="middle" >0.17</td><td align="center" valign="middle" >0.26</td><td align="center" valign="middle" >0.05</td><td align="center" valign="middle" >0.14</td><td align="center" valign="middle" >0.97</td></tr><tr><td align="center" valign="middle" >B</td><td align="center" valign="middle" >26.31</td><td align="center" valign="middle" >2.89</td><td align="center" valign="middle" >0.27</td><td align="center" valign="middle" >0.40</td><td align="center" valign="middle" >0.13</td><td align="center" valign="middle" >0.27</td><td align="center" valign="middle" >7.11</td></tr><tr><td align="center" valign="middle" >C</td><td align="center" valign="middle" >29.13</td><td align="center" valign="middle" >3.31</td><td align="center" valign="middle" >0.12</td><td align="center" valign="middle" >1.31</td><td align="center" valign="middle" >0.13</td><td align="center" valign="middle" >0.19</td><td align="center" valign="middle" >7.55</td></tr><tr><td align="center" valign="middle" >D</td><td align="center" valign="middle" >21.59</td><td align="center" valign="middle" >7.33</td><td align="center" valign="middle" >0.27</td><td align="center" valign="middle" >3.99</td><td align="center" valign="middle" >0.40</td><td align="center" valign="middle" >0.67</td><td align="center" valign="middle" >8.45</td></tr><tr><td align="center" valign="middle" >E</td><td align="center" valign="middle" >9.63</td><td align="center" valign="middle" >2.22</td><td align="center" valign="middle" >0.49</td><td align="center" valign="middle" >2.72</td><td align="center" valign="middle" >0.25</td><td align="center" valign="middle" >0.74</td><td align="center" valign="middle" >3.80</td></tr><tr><td align="center" valign="middle" >F</td><td align="center" valign="middle" >11.87</td><td align="center" valign="middle" >3.13</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >2.50</td><td align="center" valign="middle" >0.63</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >5.13</td></tr><tr><td align="center" valign="middle" >G</td><td align="center" valign="middle" >5.68</td><td align="center" valign="middle" >0.99</td><td align="center" valign="middle" >0.25</td><td align="center" valign="middle" >3.46</td><td align="center" valign="middle" >0.25</td><td align="center" valign="middle" >0.37</td><td align="center" valign="middle" >2.49</td></tr><tr><td align="center" valign="middle" >B1</td><td align="center" valign="middle" >51.87</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >2.50</td><td align="center" valign="middle" >0.63</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >3.75</td><td align="center" valign="middle" >13.25</td></tr><tr><td align="center" valign="middle" >B2</td><td align="center" valign="middle" >58.75</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >3.13</td><td align="center" valign="middle" >1.88</td><td align="center" valign="middle" >1.88</td><td align="center" valign="middle" >14.50</td></tr><tr><td align="center" valign="middle" >B3</td><td align="center" valign="middle" >72.50</td><td align="center" valign="middle" >0.63</td><td align="center" valign="middle" >−1.87</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >−1.25</td><td align="center" valign="middle" >−3.75</td><td align="center" valign="middle" >12.75</td></tr><tr><td align="center" valign="middle" >B4</td><td align="center" valign="middle" >120.63</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >1.53</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >−6.25</td><td align="center" valign="middle" >24.75</td></tr><tr><td align="center" valign="middle" >B5</td><td align="center" valign="middle" >96.79</td><td align="center" valign="middle" >3.95</td><td align="center" valign="middle" >3.95</td><td align="center" valign="middle" >3.95</td><td align="center" valign="middle" >3.95</td><td align="center" valign="middle" >7.90</td><td align="center" valign="middle" >26.47</td></tr><tr><td align="center" valign="middle" >C1</td><td align="center" valign="middle" >43.29</td><td align="center" valign="middle" >6.25</td><td align="center" valign="middle" >0.63</td><td align="center" valign="middle" >0.50</td><td align="center" valign="middle" >1.88</td><td align="center" valign="middle" >6.25</td><td align="center" valign="middle" >13.37</td></tr><tr><td align="center" valign="middle" >C2</td><td align="center" valign="middle" >52.49</td><td align="center" valign="middle" >6.25</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >2.50</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >−1.25</td><td align="center" valign="middle" >13.25</td></tr><tr><td align="center" valign="middle" >C3</td><td align="center" valign="middle" >65.10</td><td align="center" valign="middle" >6.88</td><td align="center" valign="middle" >−0.63</td><td align="center" valign="middle" >0.63</td><td align="center" valign="middle" >−0.63</td><td align="center" valign="middle" >−0.63</td><td align="center" valign="middle" >15.25</td></tr><tr><td align="center" valign="middle" >C4</td><td align="center" valign="middle" >85.01</td><td align="center" valign="middle" >8.75</td><td align="center" valign="middle" >−0.63</td><td align="center" valign="middle" >0.63</td><td align="center" valign="middle" >−1.25</td><td align="center" valign="middle" >−1.25</td><td align="center" valign="middle" >19.63</td></tr><tr><td align="center" valign="middle" >C5</td><td align="center" valign="middle" >106.25</td><td align="center" valign="middle" >100.01</td><td align="center" valign="middle" >−2.50</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >−1.25</td><td align="center" valign="middle" >−2.50</td><td align="center" valign="middle" >23.50</td></tr><tr><td align="center" valign="middle" >D1</td><td align="center" valign="middle" >29.85</td><td align="center" valign="middle" >12.79</td><td align="center" valign="middle" >3.20</td><td align="center" valign="middle" >11.73</td><td align="center" valign="middle" >3.52</td><td align="center" valign="middle" >3.20</td><td align="center" valign="middle" >16.76</td></tr><tr><td align="center" valign="middle" >D2</td><td align="center" valign="middle" >32.50</td><td align="center" valign="middle" >12.50</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >6.88</td><td align="center" valign="middle" >2.51</td><td align="center" valign="middle" >3.75</td><td align="center" valign="middle" >15.13</td></tr><tr><td align="center" valign="middle" >D3</td><td align="center" valign="middle" >35.56</td><td align="center" valign="middle" >15.80</td><td align="center" valign="middle" >−1.97</td><td align="center" valign="middle" >13.83</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >1.98</td><td align="center" valign="middle" >15.81</td></tr><tr><td align="center" valign="middle" >D4</td><td align="center" valign="middle" >57.34</td><td align="center" valign="middle" >24.58</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >20.48</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >25.40</td></tr><tr><td align="center" valign="middle" >D5</td><td align="center" valign="middle" >95.80</td><td align="center" valign="middle" >30.62</td><td align="center" valign="middle" >0.12</td><td align="center" valign="middle" >17.07</td><td align="center" valign="middle" >0.12</td><td align="center" valign="middle" >0.37</td><td align="center" valign="middle" >34.40</td></tr><tr><td align="center" valign="middle" >E1</td><td align="center" valign="middle" >10.00</td><td align="center" valign="middle" >3.75</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >6.88</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >5.63</td></tr><tr><td align="center" valign="middle" >E2</td><td align="center" valign="middle" >16.25</td><td align="center" valign="middle" >5.01</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >5.01</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >6.75</td></tr><tr><td align="center" valign="middle" >E3</td><td align="center" valign="middle" >11.85</td><td align="center" valign="middle" >3.96</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >1.98</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >−1.98</td><td align="center" valign="middle" >3.95</td></tr><tr><td align="center" valign="middle" >E4</td><td align="center" valign="middle" >18.13</td><td align="center" valign="middle" >6.41</td><td align="center" valign="middle" >−1.07</td><td align="center" valign="middle" >5.33</td><td align="center" valign="middle" >−3.19</td><td align="center" valign="middle" >−2.13</td><td align="center" valign="middle" >5.12</td></tr><tr><td align="center" valign="middle" >E5</td><td align="center" valign="middle" >21.73</td><td align="center" valign="middle" >11.85</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >13.83</td><td align="center" valign="middle" >1.98</td><td align="center" valign="middle" >1.98</td><td align="center" valign="middle" >13.04</td></tr><tr><td align="center" valign="middle" >F1</td><td align="center" valign="middle" >12.79</td><td align="center" valign="middle" >2.13</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >4.26</td><td align="center" valign="middle" >1.06</td><td align="center" valign="middle" >−1.06</td><td align="center" valign="middle" >4.48</td></tr><tr><td align="center" valign="middle" >F2</td><td align="center" valign="middle" >14.99</td><td align="center" valign="middle" >4.37</td><td align="center" valign="middle" >0.62</td><td align="center" valign="middle" >5.62</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >6.88</td></tr><tr><td align="center" valign="middle" >F3</td><td align="center" valign="middle" >16.89</td><td align="center" valign="middle" >2.82</td><td align="center" valign="middle" >−1.02</td><td align="center" valign="middle" >3.58</td><td align="center" valign="middle" >−0.77</td><td align="center" valign="middle" >−1.28</td><td align="center" valign="middle" >4.25</td></tr><tr><td align="center" valign="middle" >F4</td><td align="center" valign="middle" >20.26</td><td align="center" valign="middle" >2.13</td><td align="center" valign="middle" >−1.06</td><td align="center" valign="middle" >4.26</td><td align="center" valign="middle" >−1.07</td><td align="center" valign="middle" >−1.07</td><td align="center" valign="middle" >4.69</td></tr><tr><td align="center" valign="middle" >F5</td><td align="center" valign="middle" >30.00</td><td align="center" valign="middle" >5.62</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >4.99</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >10.01</td></tr><tr><td align="center" valign="middle" >G1</td><td align="center" valign="middle" >6.87</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >−0.63</td><td align="center" valign="middle" >0.63</td><td align="center" valign="middle" >2.01</td></tr><tr><td align="center" valign="middle" >G2</td><td align="center" valign="middle" >8.97</td><td align="center" valign="middle" >1.95</td><td align="center" valign="middle" >1.17</td><td align="center" valign="middle" >3.90</td><td align="center" valign="middle" >1.17</td><td align="center" valign="middle" >1.17</td><td align="center" valign="middle" >4.53</td></tr><tr><td align="center" valign="middle" >G3</td><td align="center" valign="middle" >8.12</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >0.00</td><td align="center" valign="middle" >3.12</td><td align="center" valign="middle" >0.63</td><td align="center" valign="middle" >1.25</td><td align="center" valign="middle" >3.25</td></tr><tr><td align="center" valign="middle" >G4</td><td align="center" valign="middle" >11.25</td><td align="center" valign="middle" >0.63</td><td align="center" valign="middle" >−0.63</td><td align="center" valign="middle" >1.87</td><td align="center" valign="middle" >−1.25</td><td align="center" valign="middle" >−1.25</td><td align="center" valign="middle" >1.63</td></tr><tr><td align="center" valign="middle" >G5</td><td align="center" valign="middle" >13.75</td><td align="center" valign="middle" >2.50</td><td align="center" valign="middle" >−1.25</td><td align="center" valign="middle" >5.00</td><td align="center" valign="middle" >0.63</td><td align="center" valign="middle" >−1.25</td><td align="center" valign="middle" >4.25</td></tr></tbody></table></table-wrap><p>Therefore, the third-order NLO coefficient (γ) of N(CH<sub>3</sub>)<sub>2</sub> and -NHCH<sub>3</sub> mainly depends on the electron-donating ability of azobenzene para-substituents. Azobenzene was introduced at the 2.10 sites of DNTTRA, and strong donor electron groups -N(CH<sub>3</sub>)<sub>2</sub> and -NHCH<sub>3</sub> were introduced at azobenzene para-position. These molecules can be designed as third-order NLO materials with good properties.</p></sec></sec><sec id="s4"><title>4. Conclusion</title><p>The structures and frequencies of 37 azobenzene derivatives such as DNTTRA (A molecule) and B-G5 were calculated by M06-2X methods at 6-311++g(d, p) computational level. The natural bond orbital (NBO) method was used to calculate the charge distribution of 11 organic molecules such as B, C, D, E, F, G, D1, D2, D3, D4 and D5. The results show that the 36 molecules of six series are D-π-A-π-D structure. The planarity of the three series of molecules obtained by introducing azobenzene at 2.10 sites, 3.9 sites, and 2.9 sites of A molecule is better than that of substituted at 1.8 sites, 1.11 sites, and 4.8 sites. The γ value of the third-order NLO coefficient of the molecule also increases obviously, reaching 10<sup>7</sup> order of magnitude atomic units (10<sup>−</sup><sup>33</sup> esu), showing good third-order NLO properties. Compared with 3.9 sites, 2.9 sites, 1.8 sites, 1.11 sites, and 4.8 sites substituted azobenzene (D) terminated at the 2.10 sites of A molecule, the minimum energy absorption peak wavelength of B, C, E, F and G molecule increased significantly. The γ value of the third-order NLO coefficient of the molecule increased by 1.1 - 3.4 times, indicating that the introduction of azobenzene at 2.10 sites of molecule A is better than that at 3.9 sites, 2.9 sites, 1.8 sites, 1.11 sites, and 4.8 sites. The third-order NLO properties of the azobenzene ring can be improved by terminating the strong donor electron groups (e.g. -NHCH<sub>3</sub> or -N(CH<sub>3</sub>)<sub>2</sub>) in the azobenzene ring, so that the third-order NLO materials with good properties can be obtained.</p></sec><sec id="s5"><title>Acknowledgements</title><p>This work is financially supported by the Project of Science and Technology Department of Sichuan Province (No. 2015GZ0343); the Project of Sichuan Provincial Department of Education (Nos. 15ZA0346, 17ZA0346); the National Natural Science Foundation of China (No. 21808092), A Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions. This research work is supported by the high performance computing platform of Jiangsu University.</p></sec><sec id="s6"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s7"><title>Cite this paper</title><p>Chen, Z.R., Zhang, Y.J., He, Z.R., Guan, Y.H., Li, Y. and Li, H.P. (2020) Electronic Absorption Spectra and Third-Order Nonlinear Optical Property of Dinaphtho[2,3-b:2’,3’-d]Thiophene-5,7,12,13-Tetraone (DNTTRA) and Its Phenyldiazenyl Derivatives: DFT Calculations. Computational Chemistry, 8, 43-60. https://doi.org/10.4236/cc.2020.84005</p></sec><sec id="s8"><title>Supporting Information</title><table-wrap id="table3" ><label><xref ref-type="table" rid="table">Table </xref>S1</label><caption><title> Eigenvalues of frontier orbitals for 37 molecules (in eV∙mol<sup>−</sup><sup>1</sup>)</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Comp.</th><th align="center" valign="middle" ></th><th align="center" valign="middle" >E1</th><th align="center" valign="middle" >E2</th><th align="center" valign="middle" >E3</th><th align="center" valign="middle" >E4</th><th align="center" valign="middle" >E5</th><th align="center" valign="middle" >F1</th><th align="center" valign="middle" >F2</th><th align="center" valign="middle" >F3</th><th align="center" valign="middle" >F4</th><th align="center" valign="middle" >F5</th></tr></thead><tr><td align="center" valign="middle"  rowspan="3"  >B3LYP</td><td align="center" valign="middle" >E<sub>(HOMO)</sub></td><td align="center" valign="middle" >−6.474</td><td align="center" valign="middle" >−6.355</td><td align="center" valign="middle" >−5.964</td><td align="center" valign="middle" >−5.731</td><td align="center" valign="middle" >−5.620</td><td align="center" valign="middle" >−6.410</td><td align="center" valign="middle" >−6.277</td><td align="center" valign="middle" >−5.881</td><td align="center" valign="middle" >−5.644</td><td align="center" valign="middle" >−5.533</td></tr><tr><td align="center" valign="middle" >E<sub>(LUMO</sub><sub>)</sub></td><td align="center" valign="middle" >−3.512</td><td align="center" valign="middle" >−3.457</td><td align="center" valign="middle" >−3.345</td><td align="center" valign="middle" >−3.274</td><td align="center" valign="middle" >−3.246</td><td align="center" valign="middle" >−3.539</td><td align="center" valign="middle" >−3.486</td><td align="center" valign="middle" >−3.381</td><td align="center" valign="middle" >−3.312</td><td align="center" valign="middle" >−3.285</td></tr><tr><td align="center" valign="middle" >E<sub>g</sub></td><td align="center" valign="middle" >2.962</td><td align="center" valign="middle" >2.898</td><td align="center" valign="middle" >2.619</td><td align="center" valign="middle" >2.457</td><td align="center" valign="middle" >2.374</td><td align="center" valign="middle" >2.871</td><td align="center" valign="middle" >2.791</td><td align="center" valign="middle" >2.500</td><td align="center" valign="middle" >2.332</td><td align="center" valign="middle" >2.248</td></tr><tr><td align="center" valign="middle"  rowspan="3"  >M06-2X</td><td align="center" valign="middle" >E<sub>(HOMO</sub><sub>)</sub></td><td align="center" valign="middle" >−7.725</td><td align="center" valign="middle" >−7.599</td><td align="center" valign="middle" >−7.163</td><td align="center" valign="middle" >−6.906</td><td align="center" valign="middle" >−6.744</td><td align="center" valign="middle" >−7.644</td><td align="center" valign="middle" >−7.519</td><td align="center" valign="middle" >−7.080</td><td align="center" valign="middle" >−6.819</td><td align="center" valign="middle" >−6.661</td></tr><tr><td align="center" valign="middle" >E<sub>(LUMO</sub><sub>)</sub></td><td align="center" valign="middle" >−2.607</td><td align="center" valign="middle" >−2.563</td><td align="center" valign="middle" >−2.479</td><td align="center" valign="middle" >−2.423</td><td align="center" valign="middle" >−2.390</td><td align="center" valign="middle" >−2.634</td><td align="center" valign="middle" >−2.591</td><td align="center" valign="middle" >−2.513</td><td align="center" valign="middle" >−2.458</td><td align="center" valign="middle" >−2.426</td></tr><tr><td align="center" valign="middle" >E<sub>g</sub></td><td align="center" valign="middle" >5.118</td><td align="center" valign="middle" >5.036</td><td align="center" valign="middle" >4.684</td><td align="center" valign="middle" >4.483</td><td align="center" valign="middle" >4.354</td><td align="center" valign="middle" >5.010</td><td align="center" valign="middle" >4.928</td><td align="center" valign="middle" >4.567</td><td align="center" valign="middle" >4.361</td><td align="center" valign="middle" >4.235</td></tr><tr><td align="center" valign="middle" >Comp.</td><td align="center" valign="middle" ></td><td align="center" valign="middle" >G1</td><td align="center" valign="middle" >G2</td><td align="center" valign="middle" >G3</td><td align="center" valign="middle" >G4</td><td align="center" valign="middle" >G5</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle"  rowspan="3"  >B3LYP</td><td align="center" valign="middle" >E<sub>(HOMO)</sub></td><td align="center" valign="middle" >−6.504</td><td align="center" valign="middle" >−6.376</td><td align="center" valign="middle" >−6.001</td><td align="center" valign="middle" >−5.798</td><td align="center" valign="middle" >−5.678</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >E<sub>(LUMO</sub><sub>)</sub></td><td align="center" valign="middle" >−3.623</td><td align="center" valign="middle" >−3.576</td><td align="center" valign="middle" >−3.464</td><td align="center" valign="middle" >−3.397</td><td align="center" valign="middle" >−3.371</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >E<sub>g</sub></td><td align="center" valign="middle" >2.881</td><td align="center" valign="middle" >2.800</td><td align="center" valign="middle" >2.537</td><td align="center" valign="middle" >2.401</td><td align="center" valign="middle" >2.307</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle"  rowspan="3"  >M06-2X</td><td align="center" valign="middle" >E<sub>(HOMO</sub><sub>)</sub></td><td align="center" valign="middle" >−7.754</td><td align="center" valign="middle" >−7.619</td><td align="center" valign="middle" >−7.208</td><td align="center" valign="middle" >−6.941</td><td align="center" valign="middle" >−6.797</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >E<sub>(LUMO</sub><sub>)</sub></td><td align="center" valign="middle" >−2.710</td><td align="center" valign="middle" >−2.671</td><td align="center" valign="middle" >−2.592</td><td align="center" valign="middle" >−2.556</td><td align="center" valign="middle" >−2.529</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >E<sub>g</sub></td><td align="center" valign="middle" >5.044</td><td align="center" valign="middle" >4.948</td><td align="center" valign="middle" >4.616</td><td align="center" valign="middle" >4.385</td><td align="center" valign="middle" >4.268</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr></tbody></table></table-wrap><table-wrap id="table4" ><label><xref ref-type="table" rid="table">Table </xref>S2</label><caption><title> The data of electronic absorption spectrum for E, F and G</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Comp.</th><th align="center" valign="middle" >Electronic transition</th><th align="center" valign="middle" >λ/nm</th><th align="center" valign="middle" >E/ev</th><th align="center" valign="middle" >f</th><th align="center" valign="middle" >Major contribution</th></tr></thead><tr><td align="center" valign="middle"  rowspan="6"  >E</td><td align="center" valign="middle"  rowspan="3"  >S<sub>0</sub>  S<sub>1</sub></td><td align="center" valign="middle"  rowspan="3"  >488.4</td><td align="center" valign="middle"  rowspan="3"  >2.54</td><td align="center" valign="middle"  rowspan="3"  >0.0011</td><td align="center" valign="middle" >HOMO-3  LUMO + 2 (34%)</td></tr><tr><td align="center" valign="middle" >HOMO-2  LUMO + 3 (29%)</td></tr><tr><td align="center" valign="middle" >HOMO-2  LUMO (10%)</td></tr><tr><td align="center" valign="middle"  rowspan="3"  >S<sub>0</sub>  S<sub>7</sub></td><td align="center" valign="middle"  rowspan="3"  >325.5</td><td align="center" valign="middle"  rowspan="3"  >3.81</td><td align="center" valign="middle"  rowspan="3"  >0.5763</td><td align="center" valign="middle" >HOMO-1  LUMO (39%)</td></tr><tr><td align="center" valign="middle" >HOMO  LUMO + 1 (16%)</td></tr><tr><td align="center" valign="middle" >HOMO-2  LUMO (12%)</td></tr><tr><td align="center" valign="middle"  rowspan="6"  >F</td><td align="center" valign="middle"  rowspan="3"  >S<sub>0</sub>  S<sub>1</sub></td><td align="center" valign="middle"  rowspan="3"  >488.2</td><td align="center" valign="middle"  rowspan="3"  >2.54</td><td align="center" valign="middle"  rowspan="3"  >0.0006</td><td align="center" valign="middle" >HOMO-3  LUMO + 2 (46%)</td></tr><tr><td align="center" valign="middle" >HOMO-3  LUMO + 3 (17%)</td></tr><tr><td align="center" valign="middle" >HOMO-3  LUMO (11%)</td></tr><tr><td align="center" valign="middle"  rowspan="3"  >S<sub>0</sub>  S<sub>7</sub></td><td align="center" valign="middle"  rowspan="3"  >324.2</td><td align="center" valign="middle"  rowspan="3"  >3.82</td><td align="center" valign="middle"  rowspan="3"  >0.4107</td><td align="center" valign="middle" >HOMO-1  LUMO (42%)</td></tr><tr><td align="center" valign="middle" >HOMO-1  LUMO + 1 (12%)</td></tr><tr><td align="center" valign="middle" >HOMO-3  LUMO (13%)</td></tr><tr><td align="center" valign="middle"  rowspan="7"  >G</td><td align="center" valign="middle"  rowspan="3"  >S<sub>0</sub>  S<sub>1</sub></td><td align="center" valign="middle"  rowspan="3"  >486.8</td><td align="center" valign="middle"  rowspan="3"  >2.55</td><td align="center" valign="middle"  rowspan="3"  >0.0001</td><td align="center" valign="middle" >HOMO-3  LUMO + 2 (38%)</td></tr><tr><td align="center" valign="middle" >HOMO-2  LUMO + 3 (32%)</td></tr><tr><td align="center" valign="middle" >HOMO-2  LUMO + 1 (7%)</td></tr><tr><td align="center" valign="middle"  rowspan="4"  >S<sub>0</sub>  S<sub>15</sub></td><td align="center" valign="middle"  rowspan="4"  >286.6</td><td align="center" valign="middle"  rowspan="4"  >4.33</td><td align="center" valign="middle"  rowspan="4"  >0.4945</td><td align="center" valign="middle" >HOMO-7  LUMO (25%)</td></tr><tr><td align="center" valign="middle" >HOMO  LUMO + 2 (18%)</td></tr><tr><td align="center" valign="middle" >HOMO-1  LUMO + 3 (16%)</td></tr><tr><td align="center" valign="middle" >HOMO  LUMO (13%)</td></tr></tbody></table></table-wrap></sec></body><back><ref-list><title>References</title><ref id="scirp.103869-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Getmanenko, Y.A., Hales, J.M., Balu, M., Fu, J., Zojer, E., Kwon, O., Mendez, J., Thayumanavan, S., Walker, G., Zhang, Q., Bunge, S.D., Bredas, J.L., Hagan, D.J., Van Stryland, E.W., Barlow, S. and Marder, S.R. (2012) Characterisation of a Dipolar Chromophore with Third-Harmonic Generation Applications in the Near-IR. Journal of Materials Chemistry, 22, 4371-4382. https://doi.org/10.1039/c2jm15599k</mixed-citation></ref><ref id="scirp.103869-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">He, G.S., Zhu, J., Baev, A., Samoc, M., Frattarelli, D.L., Watanabe, N., Facchetti, A., Agren, H., Marks, T.J. and Prasad, P.N. (2011) Twisted pi-System Chromophores for All-Optical Switching. Journal of the American Chemical Society, 133, 6675-6680. https://doi.org/10.1021/ja1113112</mixed-citation></ref><ref id="scirp.103869-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Li, Z.A., Wu, W., Ye, C., Qin, J. and Li, Z. (2009) Two Types of Nonlinear Optical Polyurethanes Containing the Same Isolation Groups: Syntheses, Optical Properties, and Influence of Binding Mode. The Journal of Physical Chemistry B, 113, 14943-14949. https://doi.org/10.1021/jp907135f</mixed-citation></ref><ref id="scirp.103869-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Mohbiya, D.R. and Sekar, N. (2018) Electronic Structure and Spectral Properties of Indole Based Fluorescent Styryl Dyes: Comprehensive Study on Linear and Non-Linear Optical Properties by DFT/TDDFT Method. Computational and Theoretical Chemistry, 1139, 90-101. https://doi.org/10.1016/j.comptc.2018.07.015</mixed-citation></ref><ref id="scirp.103869-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Hales, J.M., Matichak, J., Barlow, S., Ohira, S., Yesudas, K., Bredas, J.L., Perry, J.W. and Marder, S.R. (2010) Design of Polymethine Dyes with Large Third-Order Optical Nonlinearities and Loss Figures of Merit. Science, 327, 1485-1488. https://doi.org/10.1126/science.1185117</mixed-citation></ref><ref id="scirp.103869-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Scarpaci, A., Nantalaksaku, A., Hales, J.M., Matichak, J.D., Barlow, S., Rumi, M., Perry, J.W. and Marder, S.R. (2012) Effects of Dendronization on the Linear and Third-Order Nonlinear Optical Properties of Bis(thiopyrylium) Polymethine Dyes in Solution and the Solid State. Chemistry of Materials, 24, 1606-1618. https://doi.org/10.1021/cm3002139</mixed-citation></ref><ref id="scirp.103869-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Shkir, M., AlFaify, S., Arora, M., Ganesh, V., Abbas, H. and Yahia, I.S. (2018) A First Principles Study of Key Electronic, Optical, Second and Third Order Nonlinear Optical Properties of 3-(4-chlorophenyl)-1-(pyridin-3-yl) prop-2-en-1-one: A Novel D-π-A Type Chalcone Derivative. Journal of Computational Electronics, 17, 9-20. https://doi.org/10.1007/s10825-017-1050-3</mixed-citation></ref><ref id="scirp.103869-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Lin, J., Sa, R., Zhang, M. and Wu, K. (2015) Exploring Second-Order Nonlinear Optical Properties and Switching Ability of a Series of Dithienylethene-Containing, Cyclometalated Platinum Complexes: A Theoretical Investigation. The Journal of Physical Chemistry A, 119, 8174-8181. https://doi.org/10.1021/acs.jpca.5b03456</mixed-citation></ref><ref id="scirp.103869-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Bondu, F., Quertinmont, J., Rodriguez, V., Pozzo, J.-L., Plaquet, A., Champagne, B. and Castet, F. (2015) Second-Order Nonlinear Optical Properties of a Dithienylethene-Indolinooxazolidine Hybrid: A Joint Experimental and Theoretical Investigation. Chemistry—A European Journal, 21, 18749-18757. https://doi.org/10.1002/chem.201502728</mixed-citation></ref><ref id="scirp.103869-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Alam, M.M., Kundi, V. and Thankachan, P.P. (2016) Solvent Effects on Static Polarizability, Static First Hyperpolarizability and One- and Two-Photon Absorption Properties of Functionalized Triply Twisted M&amp;#246bius Annulenes: A DFT Study. Physical Chemistry Chemical Physics, 18, 21833-21842. https://doi.org/10.1039/C6CP02732F</mixed-citation></ref><ref id="scirp.103869-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Teran, N.B., He, G.S., Baev, A., Shi, Y., Swihart, M.T., Prasad, P.N., Marks, T.J. and Reynolds, J.R. (2016) Twisted Thiophene-Based Chromophores with Enhanced Intramolecular Charge Transfer for Cooperative Amplification of Third-Order Optical Nonlinearity. Journal of the American Chemical Society, 138, 6975-6984. https://doi.org/10.1021/jacs.5b12457</mixed-citation></ref><ref id="scirp.103869-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Yang, M., Jacquemin, D. and Champagne, B. (2002) Intramolecular Charge Transfer and First-Order Hyperpolarizability of Planar and Twisted Sesquifulvalenes. Physical Chemistry Chemical Physics, 4, 5566-5571. https://doi.org/10.1039/b207514h</mixed-citation></ref><ref id="scirp.103869-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Gao, J., Cheng, L. and Chen, X. (1999) Synthesis and Properties of Quinone Heterocyclic Organic Third Order Nonlinear Optical Materials. Chinese High Technology Letters, 2, 45-49.</mixed-citation></ref><ref id="scirp.103869-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">van Dijk, E.H., Myles, D.J., van der Veen, M.H. and Hummelen, J.C. (2006) Synthesis and Properties of an Anthraquinone-Based Redox Switch for Molecular Electronics. Organic Letters, 8, 2333-2336. https://doi.org/10.1021/ol0606278</mixed-citation></ref><ref id="scirp.103869-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Zhao, L., Wang, A.B., Wang, W.K., Yu, Z.B., Chen, S. and Yang, Y.S. (2012) Preparation and Electrochemical Performance of Aminoanthraquinone Derivative as Cathode Materials in Rechargeable Lithium Batteries. Acta Physico-Chimica Sinica, 28, 596-602.</mixed-citation></ref><ref id="scirp.103869-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">El-Aal, R.M.A., Koraiem, A.I.M. and El-Deen, N.S. (2004) Pyrazolo Quinone Heterocyclic Compounds and Metal Complex Derivatives in the Synthesis of Cyanine Dyes. Dyes and Pigments, 63, 301-314. https://doi.org/10.1016/j.dyepig.2004.03.008</mixed-citation></ref><ref id="scirp.103869-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Fujita, T., Atahan-Evrenk, S., Sawaya, N.P.D. and Aspuru-Guzik, A. (2016) Coherent Dynamics of Mixed Frenkel and Charge-Transfer Excitons in Dinaphtho[2,3- b:2’3’-f]thieno[3,2-b]-thiophene Thin Films: The Importance of Hole Delocalization. The Journal of Physical Chemistry Letters, 7, 1374-1380. https://doi.org/10.1021/acs.jpclett.6b00364</mixed-citation></ref><ref id="scirp.103869-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Guo, X., Yao, B., Jiang, G., Cheng, Y., Xie, Z., Wang, L., Ing, X. and Wang, F. (2008) Synthesis and Characterization of Color-Stable Electroluminescent Polymers: Poly(dinaphtho[1,2-a:1’,2’-g]-s-indacene)s. Journal of Polymer Science: Polymer Chemistry, 46, 4866-4878. https://doi.org/10.1002/pola.22821</mixed-citation></ref><ref id="scirp.103869-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Milvich, J., Zaki, T., Aghamohammadi, M., Roedel, R., Kraft, U., Klauk, H. and Burghartz, J.N. (2015) Flexible Low-Voltage Organic Phototransistors Based on Air-Stable Dinaphtho[2,3-b:2’,3’-f]thieno[3,2-b]thiophene (DNTT). Organic Electronics, 20, 63-68. https://doi.org/10.1016/j.orgel.2015.02.007</mixed-citation></ref><ref id="scirp.103869-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Fan, W., Yin, Z., Ma, Y., Wang, B., Chen, S., Tang, C. and Zheng, Q. (2014) Dinaphtho-s-indacene-Based Copolymers for Inverted Organic Solar Cells with High Open-Circuit Voltages. Polymer, 55, 2262-2270. https://doi.org/10.1016/j.polymer.2014.03.017</mixed-citation></ref><ref id="scirp.103869-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Ishino, Y., Miyata, K., Sugimoto, T., Watanabe, K., Matsumoto, Y., Uemura, T. and Takeya, J. (2014) Ultrafast Exciton Dynamics in Dinaphtho[2,3-b:2’3’-f]thieno[3,2- b]-thiophene Thin Films. Physical Chemistry Chemical Physics, 16, 7501-7512. https://doi.org/10.1039/c3cp54157f</mixed-citation></ref><ref id="scirp.103869-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Tabiryan, N., Hrozhyk, U. and Serak, S. (2004) Nonlinear Refraction in Photoinduced Isotropic State of Liquid Crystalline Azobenzenes. Physical Review Letters, 93, Article ID: 113901. https://doi.org/10.1103/PhysRevLett.93.113901</mixed-citation></ref><ref id="scirp.103869-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Jaunet-Lahary, T., Chantzis, A., Chen, K.J., Laurent, A.D. and Jacquemin, D. (2016) Designing Efficient Azobenzene and Azothiophene Nonlinear Optical Photochromes. The Journal of Physical Chemistry C, 118, 28831-28841. https://doi.org/10.1021/jp510581m</mixed-citation></ref><ref id="scirp.103869-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Zhao, F., Wang, C., Zeng, Y., Jin, Z. and Ma, G. (2013) Ultrafast Third-Order Nonlinear Optical Properties of an Azobenzene-Containing Ionic Liquid Crystalline Polymer. Chemical Physics Letters, 558, 100-103. https://doi.org/10.1016/j.cplett.2012.12.043</mixed-citation></ref><ref id="scirp.103869-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Virkki, M., Tuominen, O., Forni, A., Saccone, M., Metrangolo, P., Resnati, G., Kauranen, M. and Priimagi, A. (2015) Halogen Bonding Enhances Nonlinear Optical Response in Poled Supramolecular Polymers. Journal of Materials Chemistry C, 3, 3003-3006. https://doi.org/10.1039/C5TC00484E</mixed-citation></ref><ref id="scirp.103869-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Brzozowski, L. and Sargent, E.H. (2001) Azobenzenes for Photonic Network Applications: Third-Order Nonlinear Optical Properties. Journal of Materials Science: Materials in Electronics, 12, 483-489. https://doi.org/10.1023/A:1012446007088</mixed-citation></ref><ref id="scirp.103869-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Bandara, H.M. and Burdette, S.C. (2012) Photoisomerization in Different Classes of Azobenzene. Chemical Society Reviews, 41, 1809-1825. https://doi.org/10.1039/C1CS15179G</mixed-citation></ref><ref id="scirp.103869-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Li, M.M., Zhu, B.H., Ran, X., Liu, B. and Guo, L.J. (2016) Third Order Nonlinear Optical Properties of New Azobenzene Derivatives. Acta Physica Sinica, 65, Article ID: 024207.</mixed-citation></ref><ref id="scirp.103869-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Li, N., Lu, J., Li, H. and Kang, E.-T. (2011) Nonlinear Optical Properties and Memory Effects of the Azo Polymers Carrying Different Substituents. Dyes and Pigments, 88, 18-24. https://doi.org/10.1016/j.dyepig.2010.04.010</mixed-citation></ref><ref id="scirp.103869-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">Zeng, Y., Pan, Z.-H., Zhao, F.-L., Qin, M., Zhou, Y. and Wang, C.-S. (2014) Nonlinear Optical Properties of an Azobenzene Polymer. Chinese Physics B, 23, Article ID: 024212. https://doi.org/10.1088/1674-1056/23/2/024212</mixed-citation></ref><ref id="scirp.103869-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">Cai, Z., Zhou, M. and Gao, J. (2010) Studies on the Synthesis and Structure Nonlinear Optical Properties of Anthracene Two Ketones. Acta Photonica Sinica, 39, 823-828. https://doi.org/10.3788/gzxb20103905.0823</mixed-citation></ref><ref id="scirp.103869-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">Shahab, S., Filippovich, L., Sheikhi, M., Kumar, R., Dikusar, E., Yahyaei, H. and Muravsky, A. (2017) Polarization, Excited States, Trans-cis Properties and Anisotropy of Thermal and Electrical Conductivity of the 4-(phenyldiazenyl)aniline in PVA Matrix. Journal of Molecular Structure, 1141, 703-709. https://doi.org/10.1016/j.molstruc.2017.04.014</mixed-citation></ref><ref id="scirp.103869-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">Sayin, K., Kurtoglu, N., Kose, M., Karakas, D. and Kurtoglu, M. (2016) Computational and Experimental Studies of 2-(E)-hydrazinylidenemethyl-6-methoxy-4-(E)- phenyldiazenyl Phenol and Its Tautomers. Journal of Molecular Structure, 1119, 413-422. https://doi.org/10.1016/j.molstruc.2016.04.097</mixed-citation></ref><ref id="scirp.103869-ref34"><label>34</label><mixed-citation publication-type="other" xlink:type="simple">Kose, M., Kurtoglu, N., Gumussu, O., Tutak, M., McKee, V., Karakas, D. and Kurtoglu, M. (2013) Synthesis, Characterization and Antimicrobial Studies of 2-{(E)-(2- hydroxy-5-methylphenyl)imino methyl}-4-(E)-phenyldiazenyl Phenol as a Novel Azo-Azomethine Dye. Journal of Molecular Structure, 1053, 89-99. https://doi.org/10.1016/j.molstruc.2013.09.013</mixed-citation></ref><ref id="scirp.103869-ref35"><label>35</label><mixed-citation publication-type="other" xlink:type="simple">Hristova, S., Deneva, V., Pittelkow, M., Crochet, A., Kamounah, F.S., Fromm, K.M., Hansen, P.E. and Antonov, L. (2018) A Concept for Stimulated Proton Transfer in 1-(phenyldiazenyl)naphthalen-2-ols. Dyes and Pigments, 156, 91-99. https://doi.org/10.1016/j.dyepig.2018.03.070</mixed-citation></ref><ref id="scirp.103869-ref36"><label>36</label><mixed-citation publication-type="other" xlink:type="simple">Meng, Q.H., Yan, W.F., Yu, M.J. and Huang, D.Y. (2003) A Study of Third-Order Nonlinear Optical Properties for Anthraquinone Derivatives. Dyes and Pigments, 56, 145-149. https://doi.org/10.1016/S0143-7208(02)00123-7</mixed-citation></ref><ref id="scirp.103869-ref37"><label>37</label><mixed-citation publication-type="other" xlink:type="simple">Chen, Z., Li, Y., Guan, Y. and Li, H. (2019) Rational Design of the Nonlinear Optical Materials Dinaphtho[2,3-b:2’,3’-d]thiophene-5,7,12,13-tetraone (DNTTRA) and Its Phenyldiazenyl Derivatives Using First-Principles Calculations. Journal of Computational Electronics, 18, 6-15. https://doi.org/10.1007/s10825-019-01300-y</mixed-citation></ref><ref id="scirp.103869-ref38"><label>38</label><mixed-citation publication-type="other" xlink:type="simple">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., Li, X., Caricato, M., Marenich, A.V., Bloino, J., Janesko, B.G., Gomperts, R., Mennucci, B., Hratchian, H.P., Ortiz, J.V., Izmaylov, A.F., Sonnenberg, J.L., Williams-Young, D., Ding, F., Lipparini, F., Egidi, F., Goings, J., Peng, B., Petrone, A., Henderson, T., Ranasinghe, D., Zakrzewski, V.G., Gao, J., Rega, N., Zheng, G., Liang, W., Hada, M., Ehara, M., Toyota, K., Fukuda, R., Hasegawa, J., Ishida, M., Nakajima, T., Honda, Y., Kitao, O., Nakai, H., Vreven, T., Throssell, K., Montgomery, J.A., Jr., J.E.P., Ogliaro, F., Bearpark, M.J., Heyd, J.J., Brothers, E.N., Kudin, K.N., Staroverov, V.N., Keith, T.A., Kobayashi, R., Normand, J., Raghavachari, K., Rendell, A.P., Burant, J.C., Iyengar, S.S., Tomasi, J., Cossi, M., Millam, J.M., Klene, M., Adamo, C., Cammi, R., Ochterski, J.W., Martin, R.L., Morokuma, K., Farkas, O., Foresman, J. B. and Fox, D.J. (2016) Gaussian 16. Gaussian, Inc., Wallingford.</mixed-citation></ref><ref id="scirp.103869-ref39"><label>39</label><mixed-citation publication-type="other" xlink:type="simple">Li, H.P., Chang, Y.H., Zhu, W.S., Zhu, S.W., Jiang, W., Zhang, M., Zhou, Y.W., Xia, J.X. and Li, H.M. (2016) The Selectivity for Sulfur Removal from Oils: An Insight from Conceptual Density Functional Theory. AIChE Journal, 62, 2087-2100. https://doi.org/10.1002/aic.15161</mixed-citation></ref></ref-list></back></article>