OPTO-ELECTRONICS AND NONLINEAR OPTICAL PROPERTIES OF ISOINDOLINE-1,3-DIONE-FULLERENE20-ISOINDOLINE-1,3-DIONE USING DENSITY FUNCTIONAL THEORY

OPTO-ELECTRÓNICA Y PROPIEDADES ÓPTICAS NO LINEALES DE ISOINDOLINA-1,3-DIONA-FULLERENO20-ISOINDOLINA-1,3-DIONA UTILIZANDO LA TEORÍA FUNCIONAL DE LA DENSIDAD

RESAN, SAMIRA; Al-ANBER, MOHANNED

SAMIRA RESAN

Department of Physics, University of Basrah, College of Science, Karmat-Ali Road, 261004, Basrah, Iraq Irak

MOHANNED Al-ANBER

mohanned.mohammed@uobasrah.edu.iq

Department of Physics, University of Basrah, College of Science, Karmat-Ali Road, 261004, Basrah, Iraq Irak

Revista de la Facultad de Ciencias

Universidad Nacional de Colombia, Colombia

ISSN: 2357-5549

Periodicity: Semestral

vol. 12, no. 2, 2023

revistafc_med@unal.edu.co

Received: 11 February 2023

Accepted: 23 April 2023

URL: http://portal.amelica.org/ameli/journal/115/11540681154068004/

DOI: https://doi.org/10.15446/rev.fac.cienc.v12n2.107224

License Creative Commons
This work is licensed under Creative Commons Attribution-NonCommercial-NoDerivs 3.0 International.

Non-profit publishing model to preserve the academic and open nature of scientific communication

ePub generated from XML JATS4R

Abstract: Linear and nonlinear properties of isoindoline-1,3-dione-fullerene-isoindoline-1,3-dione with the five phases were studied. The geometries of the isoindoline-1,3-dione-fullerene-isoindoline-1,3-dione combination were attained using density functional theory DFT/B3LYP/6-31G+(d,p). The calculated properties are static polarizability, anisotropy of polarizability and first static hyperpolarizability at equilibrium geometry. In this article, the theoretical design of NLO molecules is discussed. The results show that there is a good electronic transfer within this molecule, and have a reasonably good propensity for nonlinear optical activity.

Keywords: NLO properties, A-D-A, Fullerene, Intramolecular charge transfer, Torsional angles.

Resumen: Se estudiaron las propiedades lineales y no lineales de isoindolina-1,3-diona-fullereno-isoindolina-1,3-diona con las cinco fases. Las geometrías de la combinación de isoindolina-1,3-diona-fullereno-isoindolina-1,3-diona se obtuvieron utilizando la teoría funcional de la densidad DFT/B3LYP/6-31G+(d,p). Las propiedades calculadas son la polarizabilidad estática, la anisotropía de la polarizabilidad y la primera hiperpolarizabilidad estática en la geometría de equilibrio. En este artículo, se discuten los resultados del diseño teórico de moléculas NLO. Los resultados muestran que existe una buena transferencia electrónica dentro de esta molécula y tienen una propensión razonablemente buena a la actividad óptica no lineal.

Palabras clave: Propiedades NLO, A-D-A, Fullereno, transferencia de carga intramolecular, ángulos de torsión.

1 Introduction

The electro-optical and nonlinear optical (NLO) materials have received considerable attention due to their wide range of potential applications (Kiven et al. 2023). When the donor and acceptor groups are present at opposite ends, it leads to significant change and plays an essential role in materials D–A systems; introducing potent electron-donor groups raises the HOMO levels. At the same time, introducing potent electron-withdrawing groups lowers the LUMO levels, resulting in the narrow HOMO–LUMO energy gap (Hashemi et al. 2019). Thus, the fundamental design criterion is the selection of proper electron-donor and acceptor units to achieve desired HOMO and LUMO levels beneficial for developing organic-optoelectronic materials. Additionally, potent electron donors and acceptors increase the delocalization of π–electrons, making the material highly polarizable, which causes remarkable optical nonlinearities in such systems. Since A-D-A materials possess distinct advantages like lightweight, good film-forming properties, relatively low manufacturing cost, biocompatibility, moderate to high conductivity, and easy tunability of desired properties, they are broadly used as active components in optoelectronic/photonic devices such as OPVs, OLEDs, OFETs, and NLO (Khalid et al. 2023).

The objective of the present paper is to design new organic compounds. Where in nonlinear optics, organic materials on delocalized π-electron systems have attracted significant research interest in carbon compounds; for this reason, the fullerene (C20) was adopted as an electron-donor group. At the same time, two molecules (isoindoline-1,3-dione), which were adopted as two terminal acceptor groups (Janjua 2017), are linked on the surface of fullerene, as shown in Figure 1. Therefore, isoindoline-1,3-dione-fullerene-isoindoline-1,3-dione (IDF) was designed as an A-D-A type of conjugated system. In addition, two terminal acceptor groups were linked with different sites on the surface of fullerene, see Figure 2, so that they make a semi-angle (φ) between these two terminal acceptor groups around the fullerene, so five phases have appeared (IDF¹, IDF², IDF³, IDF⁴ and IDF⁵, respectively). The theoretical calculations in this study are focused on the detailed analysis of the effects of various rotations (θ) of one isoindoline-1,3-dione group with the fullerene, see Figure 1, for the five phases.

Figure 1
Figure 1. Shows the general molecular design of the NLO structure (isoindoline-1,3-dione-fullerene-isoindoline-1,3-dione) and the rotation angle. Source: Elaborated by the authors.

Elaborated by the authors.

Figure 2
Figure 2. Shows the optimized geometry of the design NLO structure (isoindoline-1,3-dione-fullerene-isoindoline-1,3-dione) of the five phases (IDF¹, IDF², IDF³, IDF⁴ and IDF⁵, respectively) using the B3LYB/6-31+G(d,p) method. Source: Elaborated by the authors.

Elaborated by the authors.

2 Method and Computational Procedure

The ground-state geometry optimizations for isoindoline-1,3-dione-fullerene-isoindoline-1,3-dione of the five molecular phases, see Figure 2, have been explored by performing the density functional theory (DFT), which is one of the best computation methods of study for various types of molecules (Wazzan et al., 2016; Kadadevarmath et al., 2013). The torsion potentials were obtained for the isoindoline-1,3-dione-fullerene-isoindoline-1,3-dione) (IDF) as a function of the dihedral angle between the isoindoline-1,3-dione group and fullerene. During the scan process, the whole geometrical parameters were simultaneously relaxed so that they varied between 0◦ and 180◦ in 10◦ steps (Resan et al. 2020). The electric dipole moment μ of molecules is a fundamental interest quantity in a structural molecule. While a molecule is subject to an external electric field ϵ, the molecular charge density may reset so that the dipole moment may change (Resan et al. 2020). This change can be described as the first derivative of the molecular energy E to a component of the electric field (ϵi) that gives a component of the electric dipole moment in symbols (Shkir et al. 2022):

μi = (∂E/∂ϵi) (ϵ=0) (1)

where the total dipole moment (μ):

μ = (μx ²+μy ²+μz ²) ½ (2)

The polarizability can be understood as the gradient of the induced dipole:

αij = (∂²E/∂ϵi∂ϵj) (ϵ=0) (3)

The average static polarizability <α> tensor is defined as:

<α> = 1/3 (αxx+ αyy+ αzz) (4)

where αxx, αyy, and αzz are the polarizability matrix diagonal elements (Shkir et al. 2022). The anisotropic polarizability Δα amplitudes usually defined as (Shkir et al. 2022):

Δα = (½) [(αxx-αyy)²+(αyy-αzz)²+(αzz-αxx)²+6(αxy ²+αxz ²+ αyz ²)] ½ (5)

The equations for the calculation of molecular static first-hyperpolarizability are given as follows:

βijk =(∂³E/∂ϵi∂ϵj∂ϵk) (ϵ=0) (6)

The static first-hyperpolarizability (β) is the third-rank tensor and is expressed as a 3×3×3 matrix. Conferring to Kleinman symmetry (βxyy= βyxy= βyyx= βyyz= βyzy= βzyy,…) likewise and these 27 components can be reduced to 10 (Wagnière 1986). The complete equation for calculating the magnitudes of β from Gaussian 09 output provides ten components of this matrix as βxxx; βxxy; βxyy; βyyy; βxxz; βxyz; βyyz; βxzz; βyzz; βzzz, respectively, are reported in atomic units (AU), so that;

βtot = [βx ² + βy ² + βz ²] ½ (7)

where βx= (βxxx + βxyy + βxzz), βy= (βyyy + βyzz + βyxx), βz= (βzzz + βzxx + βzyy) (Günay et al. 2020).

The first hyperpolarizability component along the direction of the dipole moment is represented by βμ which is usually defined as:

βμ= (μxβx+ μyβy+ μzβz)/μ (8)

The XY-plane hyperpolarizability (βxy-plane), which reflects the amount of βtot in the XY-plane of the molecule, is given as (Khalid et al. 2019):

βxy-plane = βxxx+ βxxy+ βxyy+ βyyy (9)

The geometrical ground state, the electric dipole moment (μ), the mean polarizability <α>, and the first-order static hyperpolarizability (β) were calculated using the Gaussian 09W program (Frisch et al. 2009). Output files were visualized via Gauss View 5 software (Dennington et al 2009). The isoindoline-1,3-dione-fullerene-isoindoline-1,3-dione phases were determined by performing B3LYP (Becke’s three-parameter hybrid functional using the LYP correlation functional)(Soscún et al. 2002; Mbala et al. 2021). With the basis set 6-31+G(d,p) (Soscún et al. 2002) where it is a split-valence basis set, which includes both diffuse (s and p-functions) for non-hydrogen atoms and added d-polarization functions on non-hydrogen atoms and p-polarization functions for hydrogen (Dikomang 2019). Generally, many researchers have extensively used DFT/ B3LYP/6-31G+(d,p) for the first principle screening of nonlinear optical molecules (Yang et al. 2016; Parol et al. 2020; Resan et al. 2020).

3 Results and discussion

In this article, we present the theoretical design results for NLO molecules. Because of the unique ability of intramolecular charge transfer between D and A, where A-D-A systems require consideration of nonlinear optical (NLO) features(Khalid et al. 2022). The design NLO structure (IDF), which contains five phases (IDF1, IDF2, IDF3, IDF4, and IDF5, respectively), must be identified to conduct our inquiry. At first, it was necessary to find the most stable form of the design NLO structures (IDF1, IDF2, IDF3, IDF4 and IDF5, respectively) before rotation, where the relative stability (ΔE) among them are very little and the five phases have same stability approximately. To add resilience to the NLO response of (IDF) throughout its design, the geometrical effect for the five phases was tested as a function of the torsional angles (Shettigar et al. 2006). Figure 3 shows the fluctuation of relative stability ΔE with torsion angles for the five phases of IDF1, IDF2, IDF3, IDF4, and IDF5, respectively. The fifth phase's stability (IDF5) shows more instability as the torsional angles increase forward to θ=70◦ and 170◦. However, the other phases showed very low instabilities compared to the results of Bahat et al. (2009) and Resan et al. (2020). Any results for IDF5 at a high torsion angle will be ignored in favour of the other four phases due to requiring high energy for rotating, which is maybe because of the steric obstacle between the acceptor group (isoindoline-1,3-dione) and fullerene (Bahat et al., 2009; Resan et al., 2020; Bahrani et al., 2022).

Figure 3
Figure 3: Shows the relative stability ∆E with the torsional angles for the design NLO (−..−..− IDF¹ , − − −− IDF² , −.−.−.− IDF³ , − − − − − IDF⁴ and _______ IDF⁵ , respectively) using the B3LYB/6−31+G(d,p) method. Source: Elaborated by the authors.

Elaborated by the authors

The dipole moment reflects the molecular charge distribution and is given as a vector in three dimensions. Therefore, it can be used as a descriptor to depict the charge movement across the molecule. The direction of the dipole moment vector in a molecule depends on the centers of positive and negative charges. For charged systems, its value depends on the choice of origin and molecular orientation. It is well known that the higher values of dipole moment, molecular polarizability, and hyperpolarizability are important for more active NLO properties. With the dihedral angle rotation, Figure 4 depicts the dipole moment (μ) and its three primary components. As the phases change from 1 to 5, the dipole moment values fluctuate due to angle rotates. Generally, the order of μ of designed structures was: IDF⁵> IDF⁴> IDF³> IDF²> IDF¹, but they have been little fluctuation with rotation compared with the results by Resan et al. (Resan et al. 2020). As for the μx component, phases IDF⁵ and IDF⁴ are higher than the others and IDF⁴ showed a switch behavior at 100°, see Figure 4. In addition, μy has high values (IDF⁴> IDF³> IDF² respectively) and contributes significantly to total dipole moment values. Also, IDF³ showed good switch behavior for μy at 50° and 150°, with the caveat that this rotation does not need much energy, as seen in Figure 2. The total dipole moment (μ) value has been high because the component μy, so the designed molecule may be optically active in the y-direction.

The average polarizability and its components are another important molecule feature in the electronic characteristics. Therefore, calculating these properties is very important to assess the nonlinear optical potential of molecules. Values of average polarizability <α>, anisotropic polarizability Δα and major contributing tensors (αxx, αyy and αzz) of the five phases are shown in Figure 5.

Figure 4
Figure 4. Shows the electric dipole moment (µtot) and dipole components (µx, µy and µz) with the dihedral angles for the design NLO structure (isoindoline-1,3-dione-fullerene- isoindoline-1,3-dione) with five phases s (−. .−. .− IDF1 , − − −− IDF² , −.−.−.− IDF³ , − − − − − IDF4 and _______ IDF⁵ , respectively) using the B3LYB/6−31+G(d,p) method. Source: Elaborated by the authors

Elaborated by the authors

Order of average polarizability of designed structures was: IDF5> IDF4> IDF3> IDF2> IDF1, but its values are still high compared with the results by Resan et al. (2020). Likewise, this will also increase the molecule’s refractive index Resan et al. (2020), and increases the π-electron deficiency and polarizability (Migalska-Zalas et al., 2018). Generally, sizeable linear polarizability is required to obtain large hyperpolarizabilities(Ye et al., 2022). Anisotropic polarizability Δα, which is dependent on the direction of the electric field, is one of the other NLO qualities for this design NLO structures showed inverse behaviour with average polarizability results. Where the five phases show values were lower than the average polarizability. When compared to perpendicular polarizability, it can be seen that torsional angles have no effect on polarizability parallel to the symmetry axes of the molecule. The values of the two compounds αzz and αyy are higher than the αxx and showed more sensitivity to rotation dihedral angles.

Figure 5
Figure 5: Shows the average polarizability <α>, the anisotropic polarizability ∆α, and its major components (αxx, αyy and αzz) with the dihedral angles for the design NLO structure (isoindoline-1,3-dione-fullerene- isoindoline-1,3-dione) with five phases (−. .−. .− IDF¹ , − − −− IDF² , −.−.−.− IDF³ , − − − − − IDF⁴ and _______ IDF⁵ , respectively) using the B3LYB/6−31+G(d,p) method. Source: Elaborated by the authors.

Elaborated by the authors

Figure 6
Figure 6: Shows the anisotropy (κ) with the dihedral angles for the design NLO structure (isoindoline-1,3-dione-fullerene-isoindoline-1,3-dione) with five phases (−. .−. .− IDF¹ , − − −− IDF² , −.−.−.− IDF³ , − − − − − IDF⁴ and _______ IDF⁵ , respectively) using the B3LYB/6−31+G(d,p) method. Source: Elaborated by the authors

Elaborated by the authors

Figure 6 shows the results of the anisotropy (κ), which gives a measure of deviations from spherical symmetry, with rotational angles for five phases, and order due to high spherical symmetry as IDF⁴> IDF⁵> IDF¹> IDF²> IDF³ respectively. In other words, the spherically symmetric charge distribution hasn’t been an interesting response to the rotation angles as the phases change from one to five. Generally, IDF⁴ and IDF⁵ have the highest spherical symmetric charge distribution(Labidi 2016; Resan et al. 2020).

A molecule's hyperpolarizability and intramolecular charge transfer are the key factors for determining its nonlinear optical characteristics. Hyperpolarizability is a measurement of the nonlinear optical response of materials. The theoretical resolution of molecular hyperpolarizability (βtot) is useful in explaining the link between molecular manufacturing and nonlinear optical properties. Figure 7 shows the order βtot of five phases: IDF⁴> IDF³> IDF²> IDF⁵> IDF¹, respectively. The dihedral angle rotation has little influence on IDF⁴> IDF⁵ only. Figure 7 shows the results of hyperpolarizability along the dipole moment (βµ), where the relationship between hyperpolarizability and dipole moment does not change with dihedral angle rotation. Moreover, the order βµ of five phases is similar to βtot. The variation of XY-plane hyperpolarizability (βxy-plane) illustrates the amount of βtot in the XY-plane of the isoindoline-1,3-dione-fullerene-isoindoline-1,3-dione, as shown in Figure 7. Where IDF⁴ and IDF⁵ have been the highest βxy-plane value with little fluctuation along the rotation angle.

According to the results of the hyperpolarizability components (βx, βy, βz) were examined the variation with the dihedral angles for the five phases of the NLO structure (isoindoline-1,3-dione-fullerene-isoindoline-1,3-dione), as illustrated in Figure 8, where only βx for two phases, IDF⁴ and IDF⁵ respectively, have high value with significant switching behaviour for IDF⁴ by the rotation angle, which is the same as for μx component in Figure 4.

Figure 7
Figure 7: Shows the hyperpolarizability (βtot), the hyperpolarizability along the direction of the dipole moment (βµ) and the XY-plane hyperpolarizability (βxy−plane) with the dihedral angles for the design NLO structure (isoindoline-1,3-dione-fullerene-isoindoline-1,3-dione(C20)) with five phases (−..−..− IDF¹ , − − −− IDF² , −.−.−.− IDF³ , − − − − − IDF⁴ and _______ IDF⁵ , respectively) using the B3LYB/6−31+G(d,p) method. Source: Elaborated by the authors.

Elaborated by the authors

In addition, IDF³ and IDF⁴ have the highest values along the component βy and more than the values of βx. From the results of βy, it is expected that the third phase will improve the charge transfer along the plane of the central molecule axis that is between the x and y-axis. Generally, the hyperpolarizability maybe depends on the component βy, so that βtot ~ βy. Moreover, the biggest value of hyperpolarizability is seen in the βy direction, which suggests that the significant delocalization of the electron cloud is greater in that direction. In addition, this is attributed to the strong electron-donating ability of the fullerene molecule, which creates an efficient push-pull system, owing to a greater degree of charge transfer from the donor to the acceptor. Higher NLO response may be due to enhancing charge transfer from donor to acceptor, but not along the x-axis, as reported (Khalid et al. 2023).

Figure 8
Figure 8: Shows the hyperpolarizability components (βx, βy, βz) with the dihedral angles for the design NLO structure (isoindoline-1,3-dione-fullerene- isoindoline-1,3-dione) with five phases (−. .−. .− IDF¹ , − − −− IDF² , −.−.−.− IDF³ , − − − − − IDF⁴ and _______ IDF⁵ , respectively) using the B3LYB/6−31+G(d,p) method. Source: Elaborated by the authors.

Elaborated by the authors

Figure 9 shows that the energy levels of designed molecules are susceptible to the change from IDF⁵ to IDF¹. Where all these phases have little sensitivity to rotation angles except IDF⁵ and IDF². NLO design IDF¹, which is an A-D-A structure at the same straight line (φ1=180o), is an excellent group for building low band gap dyes because it significantly the lower ELUMO level, see Figure 9 (Kadhim et al. 2022). While the EHOMO level of IDF¹ is highest than the other phases. The variation of the HOMO–LUMO energy gap with the dihedral angles for the NLO structures (isoindoline-1,3-dione-fullerene-isoindoline-1,3-dione) is shown in Figure 9. As a result, maybe a red-shifted absorption spectrum can be obtained. The HOMO–LUMO energy gap value is found to be 1.148 eV is the lowest one for IDF¹, while the energy gap of 2.887eV corresponds to the IDF⁴, which has the highest hyperpolarizability. Maybe the inverse relationship between the energy gaps and the hyperpolarizability depends on the structure kind of the donor-acceptor group than their linking position effects. Generally, the change in energy gap with the dihedral angles was slightly compared to the value reported by Alyar et al. (2006), and Resan et al. (2020).

Figure 9
Figure 9: Shows the ELUMO, EHOMO and energy gap (Eg) with the dihedral angles for the design NLO structure (isoindoline-1,3-dione-fullerene- isoindoline-1,3-dione) with five phases (−. .−. .− IDF¹ , − − −− IDF² , −.−.−.− IDF³ , − − − − − IDF⁴ and _______ IDF⁵ , respectively) using the B3LYB/6−31+G(d,p) method. Source: Elaborated by the authors.

Elaborated by the authors

Frontier molecular orbital diagrams have been shown in Figure 10 for phases IDF¹, IDF², IDF³, IDF⁴ and IDF⁵ respectively. Where the highest-occupancy molecular orbitals (HOMOs) and the lowest-occupancy molecular orbitals (LUMOs) are distributed over the whole fullerene and are extended to the isoindoline-1,3-dione group for all phases. Hyperpolarizability is associated with molecular electronic distribution under the influence of the electrical field depending on loosely or tightly bound electrons. Therefore, there is no clear approximate inverse relationship between the hyperpolarizability and HOMO–LUMO energy gaps.

Figure 10
Figure 10: Frontier Molecular Orbitals of the design NLO structure (isoindoline-1,3-dione-fullerene-isoindoline-1,3-dione) with five phases (IDF¹ , IDF² , IDF³ , IDF⁴ and IDF⁵ , respectively) using the B3LYB/6-31+G(d,p) method. Source: Elaborated by the authors.

Elaborated by the authors

4 CONCLUSIONs

The study of isoindoline-1,3-dione-fullerene-isoindoline-1,3-dione, with five phases, were carried out. The equilibrium geometry, of the compound, has been investigated with the help of B3LYP density functional theory (DFT) using 6-31+G(d,p) as the basis set. The calculated HOMO and LUMO energies show that charge transfer occur within the molecule. From the results obtained, we observed that the rotation phase does have any effect on the optoelectronic properties of the molecule (title compound). At the B3LYP/6-31+G(d,p) level, there is some effect on the dipole moment, average polarizability, find first molecular hyperpolarizability for((isoindoline-1,3-dione-fullerene- isoindoline-1,3-dione).t can be concluded that A-D-A materials through investigation have high nonlinearity properties. The large value of βtot of the donor-acceptor molecules is associated with intermolecular charge transfer resulting from an electron cloud movement from an electron donor to electron acceptor groups, which is a measure of the NLO activity of the molecules. These properties of the organic compounds have made them propitious materials for electronics and optoelectronic devices with a high nonlinear optical (NLO) effect. Due to the fact that there are no corresponding experimental data in the existing literature, neither are there any other ab initio and DFT calculations on the optoelectronics properties of these molecules, we are optimistic that these results will serve as a benchmark for other studies.

Authors' Contribution

S. Resan and M. Al-Anber: conceived and designed the study. Calculations, data collection and analysis, were performed by S. Resan. The second author (M. Al-Anber) wrote the initial draft of the manuscript, and R. Resan commented and added some discussions on the initial version. Both authors read and approved the final manuscript.

Authors' declaration

Conflicts of Interest: None. We hereby confirm that all the Figures in the manuscript are ours. Ethical Clearance: The authors declare that the theoretical calculations were done according to ethical committee permission from the University of Basrah.

REFERENCES

Alyar, H., M. Bahat, ERGÜN Kasap & Z. Kantarci (2006). Torsional Barriers and Nonlinear Optical Properties of 2-, 3-, 4-Phenylpyridine Molecules. Czechoslovak Journal of Physics 56 (4), 349–358. https://doi.org/10.1007/s10582-006-0096-4.

Bahat, Mehmet & Emre Yörük (2009). A Computational Study on Structural, Electronic and Nonlinear Optical Properties of Furylpyridine Molecules. Proceedings of the 9th WSEAS International Conference on Applied Computer Science, 40–44.

Bahrani, F., Hameed, R., Resan S. & Al-anber, M. (2022). Impact of Torsion Angles to Tune Efficient Dye-Sensitized Solar Cell/Donor-π-Acceptor Model Containing Triphenylamine: DFT/TD-DFT Study. Acta Physica Polonica A 141 (6), 561–568. https://doi.org/10.12693/APhysPolA.141.561.

Dennington, Roy, Todd Keith & John Millam (2009). Semichem Inc. Shawnee Mission KS, GaussView, Version 5 (8).

Dikomang, Kago (2019). Elucidation of Structure− Property Correlations Using Experimental Nuclear Magnetic Resonance (NMR) and Fourier Transform Infrared (FTIR) Specroscopy Analysis, and Density Functional Theory (DFT) Calculations of the Difurylmethane (DFM) Molecule, and difurylmethane-methanol clusters.

Frisch, M J, G W Trucks, H B Schlegel, G E Scuseria, M A Robb, J R Cheeseman, G Scalmani, V Barone, B Mennucci & GAea Petersson. (2009). Gaussian 09, Revision A.02, Gaussian, Inc., Wallingford CT. See Also: URL: http://www. Gaussian. com.

Günay, N, Ö Tamer, D Avci, E Tarcan & Y Atalay (2020). Molecular Modelling, Spectroscopic Characterization and Nonlinear Optical Analysis on N-Acetyl-DL-Methionine. Revista Mexicana de Física, 66 (6), 749–760.

Hashemi, Daniel, Xiao Ma, Ramin Ansari, Jinsang Kim & John Kieffer (2019). Design Principles for the Energy Level Tuning in Donor/Acceptor Conjugated Polymers. Physical Chemistry Chemical Physics, 21 (2), 789–799. https://doi.org/10.1039/c8cp03341b.

Janjua, Muhammad Ramzan Saeed Ashraf (2017). Nonlinear Optical Response of a Series of Small Molecules: Quantum Modification of π-Spacer and Acceptor. Journal of the Iranian Chemical Society, 14 (9), 2041–2054. https://doi.org/10.1007/s13738-017-1141-x.

Kadadevarmath, J S, G H Malimath, N R Patil, H S Geethanjali & R M Melavanki (2013). Solvent Effect on the Dipole Moments and Photo Physical Behaviour of 2, 5-Di-(5-Tert-Butyl-2-Benzoxazolyl) Thiophene Dye. Canadian Journal of Physics, 91 (12), 1107–1113.

Kadhim, Mustafa M, Taleeb Zedan Taban, Sallal A H Abdullaha, Zainab Sabri Abbas, Safa K Hachim, Ahmed Mahdi Rheima, Haider Abdulkareem Almashhadani & Abdol Ghaffar Ebadi (2022). Design New D-π-A Materials for Sensitizers for Dye-Sensitized Solar Cells: Quantum Chemical Study. Optik, 270,169863.

Khalid, Muhammad, Riaz Hussain, Ajaz Hussain, Bakhat Ali, Farrukh Jaleel, Muhammad Imran, Mohammed Ali Assiri, et al. (2019). Electron Donor and Acceptor Influence on the Nonlinear Optical Response of Diacetylene-Functionalized Organic Materials (DFOMs): Density Functional Theory Calculations. Molecules 24 (11), https://doi.org/10.3390/molecules24112096.

Khalid, Muhammad, Mashal Khan, Iqra Shafiq, Khalid Mahmood, Muhammad Nadeem Akhtar, Javed Iqbal, Mohammed A Assiri, Muhammad Imran & Ataualpa Albert Carmo Braga (2023). Role of Donors in Triggering Second Order Non-Linear Optical Properties of Non-Fullerene FCO-2FR1 Based Derivatives: A Theoretical Perspective. Heliyon, 9 (2), e13033.

Khalid, M., Zafar, M., Hussain, S., Asghar, M. A., Khera, R. A., Imran, M., Abookleesh, F. L., Akram, M. Y., & Ullah, A. (2022). Influence of End-Capped Modifications in the Nonlinear Optical Amplitude of Nonfullerene-Based Chromophores with a D-π-A Architecture: A DFT/TDDFT Study. ACS omega, 7(27), 23532–23548. https://doi.org/10.1021/acsomega.2c02052

Kiven, Dinyuy Emmanuel, Nyiang Kennet Nkungli, Stanley Numbonui Tasheh & Julius Numbonui Ghogomu (2023). In Silico Screening of Ethyl 4-[(E)-(2-Hydroxy-4-Methoxyphenyl)Methyleneamino]Benzoate and Some of Its Derivatives for Their NLO Activities Using DFT. Royal Society Open Science, 10 (1). https://doi.org/10.1098/rsos.220430.

Labidi, Nouar Sofiane (2016). Semi Empirical and Ab Initio Methods for Calculation of Polarizability (α) and the Hyperpolarizability (β) of Substituted Polyacetylene Chain. Arabian Journal of Chemistry 9, S1252–S1259. https://doi.org/10.1016/j.arabjc.2012.01.007.

Mbala, G. F.Olinga, C. D.D. Mveme, Z. Ntieche, G. W. Ejuh, J. M.B. Ndjaka & M. T.Ottou Abe (2021). Effect of Chlorine and Bromine on the Nonlinear Optical, Electronic, Optoelectronic and Thermodynamic Properties on the BEDT-TTF Molecule: Ab-Initio and DFT Calculations. Optical and Quantum Electronics 53 (10), 1–18. https://doi.org/10.1007/s11082-021-03211-4.

Migalska-Zalas, A, · K El Korchi & T Chtouki (2018). Enhanced Nonlinear Optical Properties Due to Electronic Delocalization in Conjugated Benzodifuran Derivatives. Optical and Quantum Electronics 50, 389. https://doi.org/10.1007/s11082-018-1659-x.

Parol, Vinay, A. N. Prabhu, Md Abu Taher, Sri Ram G. Naraharisetty, N. K. Lokanath & V. Upadhyaya (2020). A Third-Order Nonlinear Optical Single Crystal of 3,4-Dimethoxy-Substituted Chalcone Derivative with High Laser Damage Threshold Value: A Potential Material for Optical Power Limiting. Journal of Materials Science: Materials in Electronics 31 (12). https://doi.org/10.1007/s10854-020-03443-2.

Resan, S., R. Hameed, A. Al-Hilo & M. Al-Anber (2020). The Impact of Torsional Angles to Tune the Nonlinear Optical Response of Chalcone Molecule: Quantum Computational Study. Revista Cubana de Fisica 37 (2), 95–100.

Shettigar, Venkataraya, P S Patil, S Naveen, S M Dharmaprakash, M A Sridhar & J Shashidhara Prasad (2006). Crystal Growth and Characterization of New Nonlinear Optical Chalcone Derivative: 1-(4-Methoxyphenyl)-3-(3, 4-Dimethoxyphenyl)-2-Propen-1-One. Journal of Crystal Growth, 295 (1), 44–49.

Shkir, Mohd., S Omar, M Aslam Manthremmel, M. A. Khan & Ahmad Umar (2022). DFT Studies on a Novel 1-Ethyl-4-Phenyl-1,5-Benzodiazepin-2-Thione Molecule and Its Derivatives for Opto-Nonlinear Applications. Engineered Science. https://doi.org/10.30919/es8d714.

Soscún, Humberto, Olga Castellano, Yaneth Bermúdez, Carlos Toro-Mendoza, Arı́stides Marcano & Ysaı́as Alvarado (2002). Linear and Nonlinear Optical Properties of Pyridine N-Oxide Molecule. Journal of Molecular Structure: THEOCHEM, 592 (1), 19–28. https://doi.org/https://doi.org/10.1016/S0166-1280(02)00222-1.

Wagnière, G. (1986). Theoretical Investigation of Kleinman Symmetry in Molecules. Applied Physics B Photophysics and Laser Chemistry, 41 (3), 169–172. https://doi.org/10.1007/BF00693918.

Wazzan, Nuha A, Ohoud S Al-Qurashi & Hassan M Faidallah (2016). DFT/and TD-DFT/PCM Calculations of Molecular Structure, Spectroscopic Characterization, NLO and NBO Analyses of 4-(4-Chlorophenyl) and 4-[4-(Dimethylamino) Phenyl]-2-Oxo-1, 2, 5, 6-Tetrahydrobenzo [h] Quinoline-3-Carbonitrile Dyes. Journal of Molecular Liquids, 223, 29–47.

Yang, Zhenqing, Chunmeng Liu, Changjin Shao, Xiaofei Zeng, Dapeng Cao, and D. Screening Yang, Z.; Liu, C.; Shao, C.; Zeng, X.; Cao (2016). Screening π-Conjugated Bridges of Organic Dyes for Dye-Sensitized Solar Cells with Panchromatic Visible Light Harvesting. Nanotechnology, 27 (26), 265701. https://doi.org/10.1088/0957-4484/27/26/265701.

Ye, Jin Ting, Xing Yi Chen & Yong Qing Qiu (2022). First Hyperpolarizabilities of Intramolecular Charge-Transfer Architectures Based on Acenaphthene Derivatives in Gas, Solution, and Solid States. Journal of Physical Chemistry A 126 (41), 7432–7441. https://doi.org/10.1021/acs.jpca.2c04380.