Structural parameter-modulated nonlinear optical amplitude of acceptor – p – D – p – donor-con ﬁ gured pyrene derivatives: a DFT approach † – 14250

In the present study, organic pyrene-based derivatives were selected for NLO investigation. The reference compound named methyl ( E )-2-cyano-3-(5-(pyren-1-yl)thiophen-2-yl)-3-acrylate ( MCPTR ) was taken for the design of its derivatives, abbreviated from MCPTD1 to MCPTD8 compounds. The nonlinear optical (NLO) properties, frontier molecular orbitals (FMOs), natural bonding orbital (NBO), and UV-vis analyses of molecules ( MCPTR – MCPTD8 ) were executed by M06 level with 6-31G(d,p) basis set. The UV-vis investigation showed that all designed compounds exhibited a redshift, and the maximum wavelength was studied in MCPTD7 (832.330 nm). The HOMO – LUMO band gaps of MCPTD1 – MCPTD8 were found to be smaller as compared to those of MCPTR (3.210 eV). The global reactivity information was correlated with band gap values; MCPTD7 , having a lower band gap, exhibited smaller hardness values (0.0321 E h ) with larger softness values (15.5763 E h ). The natural bond orbital analysis (NBO) helped to elucidate the hyper conjugative interactions, along with the stability and electron-transfer process. The dipole moment ( m ), average polarizability h a i , ﬁ rst hyperpolarizability ( b vec ) and second hyperpolarizability h g i were computed for MCPTR – MCPTD8 . Consequently, all designed compounds ( MCPTD1 – MCPTD8 ) possessed greater NLO responses than the reference compound ( MCPTR ). Interestingly, MCPTD7 showed a smaller energy gap and remarkable NLO response among MCPTD1 – MCPTD8 compounds. The highest m total , h a i , b vec and h g i values for MCPTD7 were observed as 7.200, 2.40 (cid:2) 10 (cid:3) 22 esu, 2.84 (cid:2) 10 (cid:3) 27 esu and 8.6024 (cid:2) 10 7 esu, respectively. Aptitude towards the NLO material relied upon the position of di ﬀ erent groups, the conjugated system donor and acceptor regions. The high NLO response reveals the fact that this class of pyrene-based derivatives with a thiophene linker has remarkable contributions towards NLO technological applications.


Introduction
Organic compounds with donor-p-bridge-acceptor (D-p-A) conguration showed excellent nonlinear optical (NLO) properties employed in various areas of research like molecular switching, optical modulation/memory, laser, biophysics and surface interface science. [1][2][3] Organic NLO compounds have advantages over the inorganic based on their enhanced electronic penetration, larger molecular polarizability and rapid response times. 4 The p-based electronic delocalization permits organic compounds to exhibit rapid response times, resulting in promising NLO responses. 5,6 The non-centrosymmetric pconjugated spacer compounds boost the charge transfer process among electron-donors (D) and acceptors (A), accomplishing 1 st , 2 nd and 3 rd order nonlinear amplitudes. The NLO materials have demonstrated large intramolecular charge transference (ICT), transmission of donor-acceptor p-clouds via p-linkers, higher ground electronic state dipole moments, along with excited transition states. [7][8][9] In order to compute the structural and electronic characteristics, the HOMO and LUMO band gap greatly inuenced the charge transfer in these compounds. [10][11][12] The HOMO and LUMO properties of compounds can be changed by altering conjugation, and their NLO response may be improved with a reduction of the band gap. Appropriate donors, acceptors and p-linker arrangement play a signicant role in the band gap reductions for the designed organic compound. 13,14 Therefore, it is essential to establish such D-p-A-congured compounds containing promising NLO responses. DFT studies were executed to obtain a comprehensive understanding of the NBO, FMO and UV-vis analysis.
The pyrene structure has a solid p-electron delocalization and uorescent property because of its planar aromatic features. Its solid structure is based upon the two-dimensional geometry having a strong tendency toward p assembly, and the great variation in functional groups onto graphite, fullerenes 15,16 and graphene 17,18 make it attractive towards NLO applications. 19 We have designed an appropriate donor-p-spacer-acceptor of pyrene-based molecules containing an acceptor (cyano vinyl ester) allied to pyrene centered by thienyl or phenyl conjugating linkers. 19,20 These structural variations give promising NLO results towards pi spacer end capped donor-acceptor frameworks. [21][22][23][24] These compounds may receive signicant interest in their optoelectronic applications in the zone of photonic gadgets, sensors and optical exchanging. 25 These computational methods have become a great curiosity in NLO containing a D-p-A system within pyrene-based compounds. The literature is overwhelmed with various structures, including donor-acceptor, donor-p-spacer-acceptor, acceptor-p-spacerdonor-p-linker-acceptor, donor-p-linker-acceptor-p-spacerdonor, donor-acceptor-p-acceptor and donor-donor-pacceptor. 26 The extent of p-conjugation and the nature of substituents strongly affect the NLO response of the molecules. 26,27 Experimental and theoretical research studies enhanced the NLO response, which is created through participation of strong donor and acceptor groups on one of two sides of a suitable p-conjugated structure. The pyrene units embedded from both sides with electronegative substituents will lead to surprising bathochromic shis in the optical transitions. Aer halogen substitution, a decline occurs in the pyrene transition energy because of uneven stabilization of the HOMO, exhibiting reduced orbital coefficients within the ring. The more strongly electron-withdrawing cyano entity is predicted to further red-shi the transitions of pyrene. Electronegative substitution impacts the electronic properties of the parent molecule to a much greater extent than assumed. 28 The main objective towards NLO applications is the discovery of extremely active compounds with large second-order polarizabilities (b), and hyperpolarizability that relates with the intramolecular charge transference (ICT) of the compounds. The experimental and theoretical analyses revealed that excellent second order NLO properties arise because of the D-p-A system planarity, auxiliary donors, donors, and acceptor units, and the p-electron twisted model. These ndings stimulated towards the second and third order NLO responses of pyrene-based compounds that are explored in the current ndings. This is the rst detailed DFT calculations, electronic, and NLO responses of D-p-A pyrene system compounds, and these nding may provide valuable outcomes in order to acquire new pyrene-based NLO compounds. The third-order NLO computations designed a remarkable series of pyrene derivatives with representative donor-p-acceptor structures.
Computational procedure DFT calculations were executed to obtain the electronic properties, absorption spectra and NLO properties of pyrene-rooted compounds. In recent work, quantum chemical studies of MCPTR and designed compounds MCPTD1-MCPTD8 were done by M06/6-31G(d,p) level using the Gaussian 09 program package. 29 The Gaussian 09 embedded NBO 3.1 program package was used for natural bond orbital (NBO) analysis by applying the aforesaid level of DFT and basis set. Similarly, M06/6-31G(d,p) level of TD-DFT was also employed for FMO analysis and UV-vis spectral analysis. It is used in determining photophysical properties. Gauss View 5.0 was used for organizing input les. 30 Furthermore, Gauss View 5.0, 30 Avogadro 31 and Chemcra 32 were used for the interpretation of the output results. The dipole moment was calculated by eqn (1). 33 The average polarizability hai was determined with eqn (2). 34 Analysis of the Gaussian output le yielded ten hyperpolarizability tensors oriented along the x, y, and z directions: The magnitude of the total rst hyperpolarizability (b vec ) was calculated with the help of eqn (3). 34 The second hyperpolarizability was determined with the help of eqn (4). 33 hgi ¼ 1/5[g xxxx + g yyyy + g zzzz + 2(g xxxx + g yyyy + g zzzz )] (4)

Results and discussion
This NLO study was investigated using the theoretical designing of NLO-active pyrene-based derivatives with efficient conguration of D-p-D-p-A-type. The designed compounds have the following parts: acceptor (methyl cyano), p-conjugated linkers (thiophene), and donor (pyrene). The properties of the entitled chemical systems were altered by applying some structural modications, such as acceptors and p-conjugated linkers. In particular, the following computed parameters: (i) linear polarizability (a), (ii) rst hyperpolarizability (b), (iii) second hyperpolarizability (g), (iv) absorption wavelength, (v) NBO transitions and (vi) HOMO-LUMO energies were obtained to a remarkable extent.
Usually, theoretical and experimental investigations suggest that a signicant NLO response can be achieved through the unication of strong donor and acceptor groups at the peripheral ends of a p-conjugated system. 1 In the recent study, we used the original compound named methyl (E)-2-cyano-3-(5-(pyren-1yl)thiophen-2-yl)acrylate with D-p-A conguration, which is abbreviated as MCPTR. 19 The designed series of NLO-based compounds were developed by structural modeling of MCPTR by modication of the acceptor. First of all, we designed MCPTD1 from the parent molecule (MCPTR) by keeping the D-p part preserved, as taken from the parent molecule (MCPTR), and introducing a pyren-1-yl moiety as a donor and thiophen as p-spacer units before the acceptor. Furthermore, this MCPTD1 was used to design MCPTD2-MCPTD8 compounds by modulating the end acceptor (Scheme 1 and Fig. 1). Thus, the newly designed compounds consisted of IUPAC names and their abbreviations as hexyl (E)-3-(5-(6-(5-(4-(bis (9,9- Fig. 1. It is expected that this work will trigger further research studies on their synthesis and analysis of their marvelous effects in terms of procient NLO materials. Moreover, the optimized compounds are shown in Fig. 2.

Electronic structure
Frontier molecular orbitals (FMOs) contain two main orbitals; HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital). 26 FMOs help in explaining the chemical affinity and interaction of the designed molecule with other moieties. 35 It also helps us to deduce the reactive sites in any p-electron systems. 17,36,37 The FMOs energy gap (E gap ¼ E LUMO À E HOMO ) elucidates the electron transference properties, reactivity, hardness and soness of the molecule. 38 The higher value of E gap for molecules indicates more stability, less chemical affinity and greater hardness value, which means that they may cause a change in the electronic arrangement. 12 In contrast, a molecule that is so with high reactivity and low stability has a smaller E gap value. Moreover, the E gap is regarded as a basic characteristic for analyzing the NLO properties. The HOMO is located on such orbital having an electron-donating tendency, while the LUMO is considered as the orbital that is partially lled, having an electron-accepting tendency. 19 The decreased HOMO-LUMO band gap, which can be highly polarizable, can elevate the NLO response. 37,39 FMO analysis of MCPTR and MCPTD1-MCPTD8 was carried out, and the results are collected in Table 1. Table 1 shows that the parent molecule has calculated values of HOMO/LUMO at 2.557/À5.767 eV, which is in close accordance with the experimentally determined HOMO/LUMO values of À2.780/À5.730 eV. 40 This good agreement points out that the implemented DFT computations, i.e., M06 level with 6-31G(d, p) basis set is apt to investigate MCPTD1-MCPTD8. The reference compound (MCPTR) has a larger band gap value, i.e., 3.210 eV, than all of its derivatives. This energy band gap is decreased to 2.646 eV in MCPTD1 due to the incorporation of the donor as N-(9,9-dimethyl-9H-uoren-3-yl)-9,9-dimethyl-Nphenyl-9H-uoren-3-amine and pi spacer (2-phenylthiophene) units, owing to enhancement of the conjugation and resonance factors. This band gap is further lowered in MCPTD2 as the stronger acceptor (Z)-2-(2-ethylidene-3-oxo-2,3-dihydro-1Hinden-1-ylidene)malononitrile) is substituted with the former used acceptor in MCPTR and MCPTD1, thus, reducing the band Scheme 1 Sketch map of the designed compounds.
gap value up to 2.160 eV due to the enhanced capability of the push-pull mechanism. Furthermore, the above mentioned band gap was reduced in MCPTD3 due to the incorporated one cyano and two uoro moieties at the structure of acceptor as (2-(5,6-diuoro-2-methylene-3-oxo-2,3-dihydro-1H-inden-1ylidene)malononitrile), and one more uoro moiety is introduced in the vicinity of the prior attached uoro moieties at the structure of the acceptor in MCPTD4. Consequently, the electrons migrate towards the acceptor part in these molecules due to the strong accepting nature of the cyano and uoro groups, which may be reason for the further reduction in their band gap values, such as 2.074 eV for MCPTD3 and 2.012 eV for MCPTD4. Further reduction in the band gap is examined, such as 2.002 and 1.948 eV for MCPTD5 and MCPTD6, respectively. In fact, the structure of MCPTD5 is designed via the replacement of the two uoro groups with two chloro groups at the structure of the acceptor in the MCPTD3 molecule. In the same way, the structure of MCPTD6 is designed via replacement of the three uoro groups with three chloro groups at the structure of the acceptor of MCPTD4. As it is expected, the greater electronegativity of the moieties will increase the movement of the electrons in the directions of the acceptor part because of the inductive electron withdrawal (F > Cl), 41 as observed in MCPTD3 and MCPTD4. However, at the same time, the resonance effect may compete with the inductive effect. The F and Cl groups are electrondonating due to the resonance effect (Cl > F). 39,41 The chloro groups at the structure of the acceptor are found in MCPTD5 and MCPTD6, which showed less inductive electron withdrawal as compared to the uoro groups. 42 However, a reduction in the band gap was observed in both molecules in comparison that of MCPTD3 and MCPTD4, which may be due to the (a) resonance effect, (b) better orbital matching, and (c) specic geometry of the acceptor moiety. A much lower HOMO-LUMO band gap was calculated at 1.748 eV in MCPTD7. This reduction in the band gap is due to the four substituted cyano (-CN) groups at the acceptor part having a larger (-I) effect than F and Cl. These four cyano groups can withdraw more electrons in the direction of the acceptor part. Subsequently, this factor enhances the charge transfer and reduces its band gap. This is the smallest value of E gap in the group of all investigated compounds. However, compound MCPTD8 was found with a larger value of E gap than MCPTD7 because it has just two substituted cyano (-CN) groups at the acceptor part, although one extra benzene ring was found in its acceptor part. Therefore, by structural tailoring of the donor, spacer and acceptor unit of MCPTR, the HOMO-LUMO  energy gap was remarkably reduced in all investigated compounds.
In short, the HOMO-LUMO energy gap decreases in this order: . This order conrms that the designed compounds incorporating electronegative units would be a remarkable aspect to lower the E gap values, hence, signicant the response of NLO. 43 The magnitude of the band gap describes the experience of ICT from the donor to acceptor part facilitated via pi linkers, and provides information about the NLO structure related function associations. [44][45][46] The contour sides of the frontier molecular orbitals (FMOs) are used to describe the transference of charges, as displayed in Fig. 3. In the reference compound (MCPTR), the charge density for HOMO is located over (N-(cyclohexa-1,5-dien-1-yl)-N-phenyl-4-(5-(pyren-1-yl)thiophen-2yl)aniline), while the dominant portion of the acceptor part involved in the electronic distribution is methyl (E)-2-cyano-3(thiophen-2-yl)acrylate. In the designed compounds (MCPTD1-MCPTD8), the charge density for HOMO was predominantly present over (N-(cyclohexa-1,5-dien-1-yl)-Nphenyl-4-(5-(pyren-1-yl)thiophen-2-yl)aniline) and the minor charge density was located over the rst pi-linker. The LUMO was concentrated partly on the second bridge moiety, and mainly on the acceptor moiety. Thus, a signicant charge transfer in the entitled compounds was administered from the donor to an acceptor via p-linker. This assistance of the charge transference conrms that all of the investigated compounds should be procient NLO substances. The designed compounds MCPTD1-MCPTD8 are more reactive, having larger wavelengths. Therefore, they are relatively readily polarizable compared to the reference compound. 34 The other values of HOMOÀ1, LUMO+1, HOMOÀ2, and LUMO+2 are displayed in Table S10. † The pictographic display of HOMOÀ1, LUMO+1, HOMOÀ2, and LUMO+2 is shown in Fig. S1. †

UV-vis analyses
In order to gain insight into the absorption spectra for the excited states, TD-DFT calculations were executed via M06 level of theory and 6-31G(d,p) basis set combination. The TD-DFT investigation was executed with dichloromethane (DCM) as a solvent. 47 It is expected that the polar medium is involved in the stabilization of the p-p* state associated with the n-p* state by the use of an appropriate electronic level. 39 Usually, the energy of the interactions of the compound in dichloromethane is controlled by non-covalent interactions (NCIs) and polarity effects. 42 This factor denes the dipolar interactions and hydrogen bond performing a major part, while stabilizing the rst singlet electronic level of the molecules. 48 As the solvent polarity is increased, the molecule shis to red shi. The excited state is regarded as more polar than the ground state; hence, the excited state is more stabilizing as compared to the ground state in DCM. 42 Our results revealed that all molecules (reference MCPTR and designed MCPTD1-MCPTD8) display absorbance in the visible region ( Table 2). The results discussed in Table 2 shows that the largest absorption peak of MCPTR was found at 481.21 nm, which is comparable to their corresponding experimental l max of 424 nm. The designed molecules MCPTD1-MCPTD8 cover the absorption range from 547.050 nm to 832.330 nm. Generally, the l max values of the designed MCPTD1-MCPTD8 molecules were detected to be greater in measure than the value of MCPTR (Table 2). It can be seen that the l max values are greatly inuenced by the end-capped electron acceptor motifs, which in turn push the absorption spectra to move towards a red shi. [49][50][51][52] The computed maximum absorption peak for MCPTR was obtained at 481.212 nm with 2.577 eV transition energy and oscillator strength (0.619), showing 99% molecular orbital contributions from HOMO to LUMO. The introduction of the donor as N-(9,9-dimethyl-9H-uoren-3-yl)-9,9-dimethyl-N-phenyl-9H-uoren-3-amine and pi spacer (2-phenylthiophene) in the MCPTD1 molecule effectively increased the l max value to 547.050 nm with 2.266 eV transition energy. Furthermore, the l max value of MCPTD2 was observed as 664.505 nm, which is higher than that for MCPTD1. This is because of the stronger acceptor as (Z)-2-(2-ethylidene-3-oxo-2,3dihydro-1H-inden-1-ylidene)malononitrile). For MCPTD3, the maximum absorption wavelength was further enhanced and the transition energy was reduced due to incorporation of two uoro groups in the structure of the acceptor as (2-(5,6-diuoro-2methylene-3-oxo-2,3-dihydro-1H-inden-1-ylidene) malononitrile). Additionally, in MCPTD4, the absorption wavelength was further moved to red shi, and the transition energy was reduced to 1.722 eV. The alteration in the above parameters for MCPTD4 might be due to one extra uoro moiety, which was submitted in the neighborhood of the antecedent attached uoro moieties in the structure of the acceptor. The calculated maximum wavelengths for MCPTD5 and MCPTD6 were obtained at 724.286 and 743.667 nm with 1.712 and 1.667 eV transition energies showing 91% (HOMO to LUMO) and 92% (HOMOÀ1 to LUMO), respectively. Similarly, for MCPTD7 and MCPTD8, the maximum wavelengths were calculated at 832.330 and 691.600 nm with 1.490 and 1.793 eV transition energies and oscillator strengths as 0.277 and 0.488. The best efficacy for the acceptor in MCPTD7 induces the highest red shiing in accordance with the lowering of the band gap results, as discoursed above. The l max values of all molecules are found in the subsequent decreasing order: The lowest excitation energy value as 1.490 eV is disclosed in MCPTD7, owing to the inuence of the strong electron withdrawing acceptor entity and conjugation. The increasing order of the excitation energy values is as follows: Moreover, some absorption values of MCPTR-MCPTD8 are shown in S14-S22. † The enhancement in the wavelengths and lower transition energy values uncovered that the MCPTD1-MCPTD8 molecules demonstrated larger charge transport facility. Subsequently, an easy excitation may occur between HOMO and LUMO. A former HOMO-LUMO study of MCPTD1-MCPTD8 supports the UV-vis data with respect to the minimum transition energy and higher l max values, which guide them to enhance their photoelectric effects.
To sum things up, MCPTD7 has the capacity for the lower most transition energy, minimum band gap and largest l max value, which may be contemplated like an appropriate material for using its optoelectronic properties in the NLO eld.

Global reactivity parameters
The strength of FMOs (E gap ¼ E LUMO À E HOMO ) is valuable to discover the global reactivity, such as the ionization potential (IP), global soness (s), global hardness (h), global electrophilicity index (u), electron affinity (EA), 43 electronegativity (X), and the chemical potential (m). 38,[53][54][55][56] The IP and EA are calculated by utilizing the following equations.
The chemical hardness (h), chemical potential (m), electronegativity (X), global soness (s), and electrophilicity index (u) were determined by utilizing Koopmans's theorem. 57 The ionization potential expressed the electron-accepting and electron donating abilities, and this equals the energy compulsory to extract one electron from HOMO. The reactivity and stability of the compounds are also related to the chemical potential of the species. 50 The chemical potential, energy gaps, hardness, and stability of the compound are directly related, whereas the relation is inverse in the case of the reactivity. 50 Hence, the molecule with greater energy gap is considered to be a harder molecule, and has low reactivity and more kinetic stability.
The electronegativity of the substituent controls the stabilization of the molecules. The position of the molecule holding the electronegative substituent plays a role in its stability. 58 The electronegativity explicates the capability of a representative compound to attract electrons to it. The molecule manifests more stability, lesser chemical reactivity and greater hardness if the band gap becomes higher.
The examined hardness value gradually decreased through MCPTR-MCPTD8. The hardness value calculated in MCPTD1 is 0.0486 E h , which falls to a value of 0.0397 E h in MCPTD2 and further decreases to 0.0321 E h in MCPTD7. In the same manner, the chemical potential was found to decrease, as explained in Table 3. Hence, the trend of the calculated hardness values across the designed molecules was absolutely the same as the order of the decrease of E gap . The decreasing order of the hardness and chemical potential of the investigated molecules was: However, the soness is another factor that reveals the reactivity of the molecules as directly related to their polarizability. The value of the soness expanded in the designed compounds from the MCPTR soness (8.4818 E h ) to MCPTD7 (15.5763 E h ). Interestingly, MCPTD7 shows the highest soness value, and it is regarded as the most polarizable species among all designed compounds. The higher soness values indicated that all designed compounds are polarizable. Consequently, these molecules may hold potential NLO ndings. Moreover, all global reactivity information of the studied compounds is correlated with the band gap values of MCPTR-MCPTD8. The molecules with lower band gaps exhibited smaller values of hardness and chemical potential with larger values of soness.

Natural bond orbital (NBO) analysis
NBO analysis is recognized as the most accurate framework to achieve insights into the electron-transfer process, bond interaction and elucidation of the hyper conjugative interaction among the nucleophile and electrophile. 59 The study of NBO can be useful in the interpretation of the steady image formation in donor-p-acceptor architectures, charge density transitions from fully occupied, bonded or a donor to partially-lled, non-bonded type NBOs. 60 Delocalizing reactions were evaluated by a second order perturbation access. For every transition donor (i) and acceptor (j), the stabilization energy E (2) accompanied with delocalization i / j was measured by the following equation: where E (2) is the stabilization energy, F(i / j) is the diagonal, q i is the donor orbital occupancy, and 3 j and 3 i are the off-diagonal NBO Fock or Kohn-Sham medium element. 59 There is a strong interaction between i and j when E (2) has a large value, and the whole system has a greater degree of conjugation. The results show that the withdrawing/donating ability of the substituents plays a dominant role in inducing the atomic charge on the carbon atom associated with substituents. Highly electronegative groups could increase the hyperconjugation as the nitrogen, uoro and chloro group insertion, and accordingly stabilize the lone pair of the atom in the transition states. Thus, the electronegativity plays an important role in stabilization. 58 In the second-order perturbation method, the representative interactions are obtained in Table 4. The overlapping of the orbitals causes the hyper conjugative type events to occur from s / s*, p / p*, LP / s* and LP / p*. These transitions are imperative for the intramolecular charge transfer and results in system stabilization. In p / p* transitions, the charge transfer taking place within the molecule is more dominant as compared to s / s*. In MCPTR, the important p / p* electronic transition, which expressed the highest value of stabilization energy (26.91 kcal mol À1 ), is for p (C53-C54) / p* (C26-C28). In contrast, p (C56-C58) /p* (C56-C58) was of the lowest energy transition, having 0.67 kcal mol À1 energy. Moreover, the s / s* transition as s (C26-H27) / s* (C53-S55) had 7.79 kcal mol À1 as highest energy and s (C44-C47) / s* (C47-H52) having 0.51 kcal mol À1 lowest energy transitions. Furthermore, other transitions such as LP / p* had the highest value at 51.37 kcal mol À1 for the LP2 (O33) / p* (C31-O32) transition, and LP / s* had the lowest stabilization energy value of 0.54 kcal mol À1 for the LP1 (S55) / s* (C17-H18) transition. A bunch of transitions exhibiting conjugation are demonstrated in Table S23 Table S25. † In the chemical structure of MCPTD3, the p / p* electronic transition exhibiting the highest value of stabilization energy (30.08 kcal mol À1 ) was p (C27-C28) / p* (C25-C103). In contrast, the transition p (C111-N112) / p* (C109-N110) exhibited the lowest energy as 0.69 kcal mol À1 . The highest energy from the s / s* transition s (C25-H26) / s* (C27-S29) was found as 9.05 kcal mol À1 . In contrast, the lowest energy in s (C28-H31) / s* (C25-H26) was uncovered as 0.50 kcal mol À1 . Furthermore, LP1 / p* had the highest value of 28.94 kcal mol À1 in LP1 (N40) / (C119-C123), while the lowest value (0.63 kcal mol À1 ) was observed in LP1 / p* as the LP1 (S29) / (C16-H17) transition. Other transitions exhibiting established conjugation are summarized in Table S26. † In the framework of MCPTD4, the p / p* transition had the highest value of stabilization energy (30.68 kcal mol À1 ) for p (C27-C28) / p* (C25-C103). In contrast, p (C110-N111) / p* (C108-N109) was found to have the lowest energy as 0.70 kcal mol À1 . The s / s* transition s (C25-H26) / s* (C27-S29) had 9.11 kcal mol À1 as the highest energy. The s (C84-H87) / s* (C61-C84) consisted of the lowest energy, having 0.50 kcal mol À1 . The other LP / p* had the highest value at 28.90 kcal mol À1 for the transition LP1 (N40) / p* (C118-C122) and the lowest value (0.52 kcal mol À1 ) was for the LP2 (F113) / s* (C100-C101) transition. Other transitions exhibiting conjugation are shown in Table S27. † For MCPTD5, the p / p* transition that has the highest value of stabilization energy was of p (C27-C28) / p* (C25-C103) with 30.59 kcal mol À1 . The p (C110-N111) / p* (C108-N109) was of the lowest energy, having 0.69 kcal mol À1 . The s / s* transition as s (C25-H26) / s* (C27-S29) had 9.07 kcal mol À1 . The s (C118-H122) / s* (N40-C120) transition was of the lowest energy, having 0.52 kcal mol À1 . LP / p* had the high value of 28.86 kcal mol À1 for LP1 (N40) / p* (C116-C120), and the lowest value (0.52 kcal mol À1 ) was observed for the LP2 (F113) / s* (C100-C101) transition. Other transitions exhibiting conjugation can be seen in Table S28. † In MCPTD6, the p / p* transition having the highest value of stabilization energy at 31.15 kcal mol À1 was for the p (C27-C28) / p* (C25-C103) transition, and the p (C49-C50) / p* (C41-C43) was of the lowest energy, having 0.70 kcal mol À1 . The s / s* transition as s (C25-H26) / s* (C27-S29) had 9.22 kcal mol À1 as the highest energy. The s (C118-H122) / s* (N40-C120) had the lowest energy at 0.52 kcal mol À1 . The LP / p* transition had the highest value at 28.79 kcal mol À1 for LP1 (N40) / p* (C116-C120), and the lowest value (0.57 kcal mol À1 ) was observed for the LP2 (Cl125) / s* (C97-C98) transition. Other transitions exhibiting conjugation are listed in Table  S29. † For MCPTD7, the p / p* transition having the highest value of stabilization energy was for p (C27-C28) / p* (C25-C103) as 32.49 kcal mol À1 , whereas p (C110-N111) / p* (C108-N109) consisted of the lowest energy as 0.68 kcal mol À1 . The s / s* transition: s (C25-H26) / s* (C27-S29) had 9.23 kcal mol À1 as the highest value, and s (C28-H31) / s* (C25-H26) comprised the lowest energy at 0.50 kcal mol À1 . Moreover, LP / p* had the highest value at 29.41 kcal mol À1 for LP2 (S29) / p* (C30-C32), while the lowest value (3.08 kcal mol À1 ) was observed for the LP1 (N40) / s* (C49-C50) transition. The remaining transitions exhibiting conjugation are listed in Table S30. † For MCPTD8, the p / p* transition contained the highest value of stabilization energy at 29.53 kcal mol À1 for p (C27-C28) / p* (C25-C103). However, the lowest energy (0.60 kcal mol À1 ) was uncovered for p (C49-C50) / p* (C41-C43). Furthermore, the s / s* transition, such as s (C25-H26) / s* (C27-S29), comprised the highest energy at 9.12 kcal mol À1 . The s (C30-C32) / s* (C27-S29) was found to have the lowest energy at 0.51 kcal mol À1 . Furthermore, LP / p* has the highest value at 26.58 kcal mol À1 due to LP2 (S29) / p* (C30-C32), while the lowest value observed in LP1 (N40) / s* (C49-C50) was the 3.08 kcal mol À1 transition. Many transitions exhibiting conjugation are listed in Table S31. † Among all the entitled compounds, MCPTD7 had the highest stability due to the extended hyper-conjugation with 32.49 kcal mol À1 . All of the entitled compounds exhibited more stability in comparison with the reference compound. The overall stability order is as follows: Therefore, NBO analysis of these compounds entails that the extended hyper conjugation and robust intramolecular transference of the charge play a remarkable role in stabilizing these compounds, and signify the charge transfer properties that are important for NLO properties.

Nonlinear optical (NLO) properties
The electronic characteristics of the designed compounds were measured for the optical analyses, such as linear response: polarizability (a) and nonlinear responses: rst hyper polarizability (b) and second hyper polarizability (g). 61 An applied external electric eld interacts with the electronic densities of the extended conjugated systems. Consequently, it changes their dipole moments. 62 The delocalization strength of the electrons in any system can be enhanced using substituents (electron donor and acceptor) at reasonable places of the conjugated system. 43,63 This enhancement in the asymmetric distribution of the electrons can be a reason for the promising nonlinear optical phenomena of molecules. 62 In compounds MCPTD1-MCPTD8, the resonance of the p-electrons between the phenyl rings has been increased, and the C-C p-bond is not positioned. Therefore, it takes part in the extended conjugation of the two or more phenylic rings. 64 So, the delocalization increased from MCPTD1-MCPTD8 through the addition of the chlorine, uorine and cyano groups as a substituent, which dris electrons marvelously. Thus, introducing a pi-spacer to link different appropriate acceptor parts is effective to decrease the E gap , and thus to increase polarizabilities. 62 Moreover, the polarity has a signicant impact on the dipole moment (m) by enhancing their NLO values. 43 Dipole moments are produced from electronegativity differences. The larger the difference in the electronegativity, the larger the dipole moment. 65 The dipole moment (m) can be seen by the results mentioned in Table 5. The dipole moment (m), average polarizability hai, second-order polarizability (b vec ), and second-order hyperpolarizability hgi values of the entire compound are calculated in esu units, 66 and are listed in Table 5, while their major contributing tensors are discussed in detail in S11-S13. † The reference molecule (MCPTR) had the lowest average polarizability. The strong acceptor moiety with stronger electron-withdrawing groups caused the highest value of average polarizability. Therefore, MCPTD7 had the predominant value of the dipole moment as 7.200 a.u. The decreasing order of the dipole moment values was observed for the designed compounds with the reference as MCPTD7 > MCPTD5 > MCPTD6 > MCPTD3 > MCPTD4 > MCPTD8 > MCPTD2 > MCPTR > MCPTD1. The addition of the electronegative atom pulled the electron density to that atom; in this way, it rendered the other parts of the molecule as a result of this increased linear polarizability. 67 The trend for the average linear polarizability values of the designed compounds with the reference was observed as: Similar to the dipole moment, the highest value of the average polarizability was observed for the designed compound MCPTD7 due to the presence of four cyano groups, which enhanced the electron density towards the A unit and made it more electronwithdrawing. Hence, a remarkable increase in its value was observed, i.e., 2.40 Â 10 À22 esu. Furthermore, a systematic relationship was obtained between the molecular structures and b values. The b parameter usually increased with the strength of the acceptor substituents attached with MCPTR, like the uoro, chloro and cyano groups, which contribute to the molecular nonlinearity. Moreover, the contribution of the extension of the conjugated system to b is dominant too as the substitutions take place. Thus, the b vec value for all designed compounds decreased in the following order: has been shown to possess the highest value of b vec as 2.84 Â 10 À27 esu compared to the other designed compounds. However, MCPTD1 and MCPTD2 showed minimal b vec values, owing to the poor ICT. The introduction of a uorine (electron acceptor) in MCPTD3 and MCPTD4 enhanced the accepting power on the opposite end. MCPTD3 has an effectively high NLO value due to the strong accepting tendency of the two uorine atoms attached to it. There is the phenyl ring as a donor moiety and three atoms of the uoro (-F) moiety competing on the opposite end in the case of MCPTD4. Thus, the b vec value was boosted from 1.66 Â 10 À27 esu of MCPTD3 to 1.86 Â 10 À27 esu in MCPTD4. The installation of a chlorine electron acceptor onto the phenyl ring in MCPTD5 and MCPTD6 enhances the accepting power again. MCPTD5 has a high NLO value because of the strong accepting ability of the chlorine. There is the phenyl ring as a donor moiety and three atoms of the chloro (-Cl) moiety competing on the other end. Thus, the b vec value is shied from 1.92 Â 10 À27 esu in MCPTD5 to 2.02 Â 10 À27 esu in MCPTD6. MCPTD7 has the highest NLO value of all because it consists of four cyano (-CN) groups. In these molecules, the maximum and minimum NLO responses can be ascribed to the prociency of the charge transfer from the donor to acceptor through their corresponding p-linked interactions. Moreover, upsurges in the hyperpolarizability values in the molecules discussed here arise in association with the delocalization of the p-electrons. This delocalization decreases the HOMO-LUMO energy difference. The calculated values of the average polarizability and second-order polarizability of the investigated molecules are signicantly greater than that for urea (b vec (urea) ¼ 0.372 Â 10 À30 esu), a standard molecule for the analysis of the NLO response. 68 The hyperpolarizability values of MCPTR, MCPTD1, MCPTD2, MCPTD3, MCPTD4, MCPTD5, MCPTD6, MCPTD7, and MCPTD8 were determined to be 557, 1709, 3551, 4459, 5004, 5158, 5445, 7632 and 4474 times larger than the second-order polarizability value of urea, respectively.
The rst hyperpolarizability values of all designed compounds were uncovered with the opposite order of the E gap between the HOMO-LUMO orbitals. The abovementioned factors affect the second hyper polarizability hgi amplitudes for MCPTD1-MCPTD8 too. 69 The maximum hgi amplitude (8.6024 Â 10 7 esu) was also found in MCPTD7. Overall, the decreasing trend of all of the studied compounds was observed as: It may be deduced from the aforementioned discussion that the electronwithdrawing nature of the acceptors with the pi-conjugation and thienyl moieties played a pivotal role and yielded remarkable NLO amplitudes. 69

Conclusion
A series of structurally related molecules (MCPTD1-MCPTD8) were designed from the pyrene-based synthesized compound (MCPTR). The inuence of different acceptors of MCPTD1-MCPTD8 was carefully explored for the NLO properties. These results revealed the fact that the acceptors have promising inuence over the donor-p-acceptor structures, and strongly tuned the entire features of the designed compounds over the synthesized compound. MCPTD1-MCPTD8 showed the highest wavelength with low transitional energy, a broader extent of orbital formation, and oscillating strength values compared to MCPTR. The larger red shied spectrum value (l max ¼ 832.330 nm) was found for MCPTD7. The FMO-based investigation indicates that HOMOs are located over the donor moiety and its nearby p-linker, while LUMOs are settled on the acceptor in all of the molecules. The HOMO-LUMO energy-gap was reduced from 3.210 eV to À2.088 eV for MCPTR-MCPTD8, respectively. Interestingly, the HOMO-LUMO energy-gap of MCPTD7 was decreased by up to 1.748 eV. Contrarily, the NLO values have been enhanced from 2.07 Â 10 À28 esu to 1.66 Â 10 À27 esu for MCPTR-MCPTD8. Surprisingly, the hai, b vec and hgi values were found to be 2.40 Â 10 À22 , 2.84 Â 10 À27 and 8.6024 Â 10 7 esu, respectively, for MCPTD7, which was greater than the reference (MCPTR) and MCPTD1-MCPTD8. Experimental researchers may be attracted to our designed compounds owing to the impressive NLO results. We could enlist the promising abilities of our designed compounds for establishing the better and responsive NLO. The larger NBO and NLO response in MCPTD7 was produced as a consequence of the delocalization of the psystem. Delocalization lessens the HOMO-LUMO band gap and ensures stability of the molecules. Thus, it can exhibit excellent material properties in photoresponsive compounds.

Conflicts of interest
There are no conicts to declare. Table 5 The computed dipole moment (m), polarizability hai, second-order polarizabilities (b), and second-order hyper polarizability hgi of the entitled compounds (MCPTR-MCPTD8) in (a.u)