Pseudo-Stilbene-and Azobenzene-Type Systems for Optical Frequency Conversion: Estimating the First-Order Molecular Hyperpolarizability

: This study investigates the potential of a set of pseudo-stilbene and azobenzene molecular structures to become optical frequency converters for optical communications based on a detailed exploration of the first-order molecular hyperpolarizability ( β HRS ), which is the microscopic counterpart of second harmonic generation (SHG). β HRS values were obtained via quantum chemical calculations using the Gaussian 16 software package in solvent and gas-phase media at different wavelengths, i.e., 1064 nm, 1310 nm, and 1510 nm. The latter two wavelengths are of particular interest for optical communications. Our study focused on discerning how the molecular structure influences the β HRS response, explicitly highlighting the influence of the azomethine group (CH=N). The results revealed that the molecular planarity, affected by this group, plays a crucial role in modulating the optical properties. The highest β HRS value in a solvent medium using the CAM-B3LYP/6-311+G(2d,p) level of theory achieved in this work was around 1400 × 10 − 30 cm 4 startvolt − 1 , four orders of magnitude higher than KDP (0.2 × 10 − 30 cm 4 startvolt − 1 ), which is a reference in SHG experiments at 1064 nm. The highest calculated β HRS value at the same level of theory and solvent at 1310 nm and 1550 nm was 631 × 10 − 30 cm 4 startvolt − 1 and 456 × 10 − 30 cm 4 startvolt − 1 , respectively. All these values belong to molecular structures with azo-coupling with donor (4-NMe 2 ) and acceptor (4 ′ -NO 2 ) peripheral groups, designated as AB-3.


Introduction
Photonics is a well-established field, transversing various science and technology domains, such as optical communications [1,2].Since 2020, internet traffic has rocketed, mainly due to the increase in many digital activities because of the lockdowns and restrictions provoked by the coronavirus pandemic [3].In telecommunications, the broadband concept means wide bandwidth data transmission that transports multiple signals at a wide range of frequencies.Suppose such a signal is an optical rather than an electric signal.In this case, a photonic device that converts the light frequency, i.e., an optical frequency converter (OFC), is essential to increase the optical communications broadband.To choose a material to be used in developing an OFC, firstly, we need to characterize it in terms of nonlinear optics (NLO) [4], particularly by estimating the second harmonic generation (SHG) phenomenon [5].In other words, we need to assess the material's ability to double the frequency of incident light.Several materials can generate SHG efficiently [6][7][8][9][10][11], but organic compounds are best suited for optical communications due to their intrinsic ultra-fast polarization response.In fact, optical nonlinearities, such as optical switching and regeneration, can provide a massive speed advantage over electronic techniques.For example, landmark experiments have reported all-optical switching at well over Terabaud rates [12].Therefore, fundamental research to understand the SHG response of organic materials is crucial to developing novel OFC.
Numerous organic materials have already investigated [13][14][15][16][17][18], but the search for novel or modified organic molecular structures with significant NLO responses continues.However, attempting to synthesize a novel molecular structure may be futile and expensive without prior fundamental research and understanding of its inherent optical properties.Moreover, NLO characterization using ultra-short pulses is costly due to the high-end photonics infrastructure needed.For instance, utilizing the hyper-Rayleigh scattering (HRS) method with a femtosecond laser system stands out as a trusted experimental approach to ascertain the first-order molecular hyperpolarizability in solution, denoted as β HRS .The emanated light from β HRS results in incoherent second harmonic generation (ISHG), a crucial spectroscopic measure when assessing the prospects of a proposed OFC.Conversely, adopting a theoretical strategy, specifically quantum chemical calculations (QCCs), offers a more cost-effective means to estimate β HRS .To illustrate, evaluating the β HRS value for a molecular structure comprising roughly 40 atoms can be completed in under a day with a 16-core setup using the Gaussian software package.This timeframe encompasses the entire spectrum of tasks, from molecular design and optimization to the post-processing of the resulting QCC data.Naturally, the duration can vary based on the chosen theoretical level, as noted in reference [19].
Still, not all organic systems are suitable for exploring NLO properties.Usually, a π-conjugated molecular structure containing an electron donor (D) and an electron acceptor (A) is widely known as a push-pull system (D-π-A).These molecular systems have become intensely explored [20][21][22], and the interaction between groups D and A through a bridge of π-electrons is called intramolecular charge transfer (ICT) [23][24][25].For example, when charge transfer occurs between electron-donating and electron-accepting groups connected by a conjugated π-bridge, the chromophore exhibits a significant increase in dipole moment [25].Thus, the impact of ICT on the optical properties of π-conjugated molecules can be observed by the geometry change, the absorption spectrum, and their dependence on the first-order hyperpolarizabilities of molecules [18,[26][27][28][29].
Chemically, stilbenes are a common and integral feature of a variety of natural products and synthetic structures.They represent a significant "privileged scaffold" and are also of great interest in medicinal chemistry [30].These scaffolds are also present in a wide range of molecules which have critical applications in material sciences [31,32].Considering the biological relevance and importance of stilbenes in material science, many related synthetic structures can be found in the literature, known as pseudo-stilbenes.These structures have 4 and 4 ′ positions of the two azo rings with electron-donating and electron-withdrawing groups at opposite ends of the aromatic system [33], which include azomethines (-CH=N-) [34][35][36][37], azobenzenes (-N=N-) [20,[38][39][40][41][42], and azo-coupling (-N=N-X-Y-) [26,37,[43][44][45].Such structures are a well-known class of materials with reported NLO responses due to their potential for applications in photonic devices [46].Furthermore, these systems have been widely reported experimentally and theoretically due to their interesting physicochemical properties, particularly with excellent values of β HRS [34,[45][46][47][48]. Thus, new organic materials containing the core of these essential groups with highly delocalized π-electrons and their theoretical calculations play a crucial role in NLO applications [26,29,49].
Encouraged by the preceding works, we have decided to estimate the magnitude of the β HRS of a series of pseudo-stilbenes (PS) and azobenzene-based (AB) compounds.In fact, to the best of our knowledge, only the static first-order molecular hyperpolarizability (β 0 ) has been reported [49], but not using the Gaussian software package.Also, we have used the most common wavelengths in the field of optical communications (λ = 1064 nm, 1310 nm, and 1550 nm) to establish an empirical correlation for future OFC development.Therefore, in this work, QCCs were carried out at the time-dependent density functional theory (TD-DFT) level to estimate the influence of functional groups on the magnitude of the β HRS using Gaussian 16 and Hyper-QCC [19], which is post-processing software.

Investigated Compounds
From a synthetic point of view, stilbenes (-CH=CH-) are accessible via organic compound carbon-carbon coupling reactions [50].On the other hand, pseudo-stilbenes can be prepared from the coupling of two different anilines in the case of azobenzenes (-N=N-) [51] and azo-coupling (-N=N-X-Y-), or from the condensation reaction of aldehyde and aniline in the case of azomethines (-CH=N-) [52].A series of synthetically feasible stilbenes and pseudo-stilbenes was planned for this study.More details regarding the synthesis of the investigated compounds can be found in the Supplementary Information (SI).
The investigated molecules are shown in Figure 1, which are a series of pseudostilbenes (System I) and azobenzene-based compounds (System II) containing a donor dimethylamine group (4-NMe 2 ) and an acceptor nitro group (4 ′ -NO 2 ).

Computational Methodology Details
Dynamic and static first-order hyperpolarizabilities were estimated by using QCCs performed using TD-DFT via the Gaussian 16 (rev.B.01) [53] software package.Three methods were used, i.e., the exchange-correlation functional Becke; the three-parameter, Lee-Yang-Parr (B3LYP) method [54]; and its long-range corrected version including the Coulomb-attenuating method (CAM-B3LYP) [55].Also, a variation of a hybrid functional by Truhlar and Zhao (M06-2X) [56] optimized to predict non-covalent interactions complemented the level of theory.All functionals were combined with a triple-zeta split-valence Pople basis set [57] with two additional polarization functions for non-hydrogen atoms, one additional polarization function on hydrogen atoms, and diffuse functions on hydrogen and non-hydrogen atoms: 6-311++G(2d,p).
The computational methodology consisted of three steps: molecular structure geometry optimization, calculating the one-photon vertical excitation energies, and computing the tensor components of the first-order molecular hyperpolarizability.Then, Hyper-QCC [19], a post-processing software dedicated to computing third-order ranking tensor components in a fast and reliable way, was used for the final calculation of β 0 and β HRS values.
The first computational step was the geometry optimization calculation in gas-phase and solvent medium (dimethyl sulfoxide-DMSO).The solvation effect was obtained by employing a polarizable continuum model (PCM) using the integral equation formalism variant (IEF-PCM) [58].The level of theory CAM-B3LYP was the only method used to calculate the energy minima of all compounds.The second step was to obtain the vertical excitation energies in DMSO using the same level of theory of geometry optimization.The vertical excitation energies allow us to verify the spectral transparency window of all compounds.The third step was to calculate the static (ω = 0 a.u.) and dynamic (ω ̸ = 0 a.u.) tensor components, i.e., β ijk .As mentioned before, the frequencies used to achieve the dynamic values are equivalent to the wavelengths of 1064 nm (0.0428 a.u.), 1310 nm (0.0348 a.u.), and 1550 nm (0.0294 a.u.).All three hybrid functionals (B3LYP, CAM-B3LYP, and M06-2X) were used to obtain the hyperpolarizability tensor components in the solvent and gas phase.The Gaussian output files provided in the third step were post-processed in Hyper-QCC software v.1.3.2 to gather the final values of β 0 and β HRS .
The β HRS tensor is decomposed as the sum of dipolar (J = 1) and octupolar (J = 3) tensorial components.The relationships between the dipolar and octupolar components and the cartesian components of β HRS are described in Equations ( 2) and (3).

Results and Discussion
Both investigated systems (type-I and -II) are push-pull systems in which the linear and nonlinear optical properties depend on the functional and peripheral groups present in each molecular structure, as well as on the increase in pi-conjugate bonds [26,61].Thus, am understanding of the strength of the donor and acceptor substituents and the nature of the π-bridge is essential to discuss how ICT will affect the β HRS values.As mentioned before, the β HRS values of all compounds were calculated in the gas phase and medium solvent at their minimum energies.Section S1 of the SI presents all optimized cartesian coordinates.The 3D illustration in Figure 2 shows the visualization in different planes of the molecules: (a) front view, (b) planar view, and (c) side view.This visualization was generated using GaussView 6, a graphical interface used with Gaussian.
With the minimum energy obtained, we investigated planarity (geometric) parameters, such as the bond length, bond angle, and twist (torsion) angle for all compounds in solvent medium (DMSO).Planarity implies that all atoms within the molecule are arranged in the same plane, meaning that the angles between atoms are at 180 • , with no significant deviations from the ideal planar structure.The planarity of the molecule can be experimentally evaluated through measurements like the average deviation of atoms from the plane defined by a set of reference atoms or through theoretical calculations.In the compounds investigated, central torsion angles (-X = Y-) and torsion angles near the NMe2 substituent were selected to evaluate each molecule's planar structure.Planar conformations can affect the interactions of donor and acceptor groups throughout the molecule, resulting in varying electronic mobilities and optical nonlinearities.Third-party software was used to obtain parameters related to bond length ( Ǻ), bond angle ( Ǻ), and torsion angle ( • ).The methodology employed for these determinations is outlined as follows: (i) import the input file containing the XYZ coordinates of the optimized compound; (ii) examine the molecular structure in three dimensions to measure the distance, angles, and torsion with tools available in the software; (iii) calculate distances by selecting pairs of atoms; (iv) calculate angles by selecting three atoms that make up the angle of interest; and (v) calculate molecular torsion by selecting four adjacent atoms along the main chain and near the substituents at the end of the molecular chain.All the methodology results are presented in Table 1.The planarity parameters were compared from the central core of the molecules to the connection between the donor and acceptor groups, aiming to analyze the main structural changes, which can be a reasonable indication of charge transfer.In all investigated systems, it is possible to observe favorable planarity; however, the molecular structures containing the azomethine group (CH=N), i.e., PS-2 and AB-2, deviate slightly from molecular planarity.This variation suggests that it is due to three main factors: increased conjugation and the interaction between the functional group and both the electron-donating (-NMe 2 ) and -accepting (-NO 2 ) groups.On the other hand, for compounds AB-1 and AB-3, we noticed that, despite the increase in the π-conjugated chain, the planarity of the molecules was preserved.Consequently, this can directly influence the optical properties of these systems.
In this context, when examining the spectroscopic parameters listed in Table 2, it is noteworthy that compounds with elevated dipole moment values exhibited minimal deviations from their planar configuration.We analyzed each structure's main torsions and bond angles to verify this result.Initially, we observed some deviations in the torsions between the functional groups that bind the phenyls, which presented the following values: PS-1 (179.99 • ), PS-2 (178.07 • ), PS-3 (179.96• ), AB-1 (179.99 • ), AB-2 (179.60 • ), and AB-3 (179.96• ).From the analysis, we noticed that the PS-2 compound shows a twisting of 1.93 • in relation to other planar structures (PS-1, PS-3, AB-1, and AB-3), while compound AB-2 exhibits a variation of 0.4 • .Therefore, when analyzing the torsions within the plane of the functional groups -NO 2 and -NMe2, we observed that the greatest torsions occur towards the donor group, with the following values: PS-1 (179.98 Ǻ), PS-2 (169.59 • ), PS-3 (169.67 • ), AB-1 (179.92 • ), AB-2 (−166.35• ) and AB-3 (169.43 • ).The lowest values were for compounds AB-2, PS-2, AB-3, and PS-3, respectively.This observation is relevant, as it can influence the NLO response, given that these structural changes may lead to modifications in the dipole moment values and the first molecular hyperpolarizability [62][63][64][65].In general, larger dipole moments tend to result in larger NLO responses.Therefore, organic compounds with larger dipole moments exhibit more robust nonlinear optical properties.However, the relationship between dipole moments and NLO responses can be affected by a variety of factors, including the molecular structure and the chemical environment [18].When analyzing the bond angles, we observed similar behavior in most compounds, except for PS-2.In the specific case of PS-2, we noted that the C3-C1-N2 angle, toward the phenyl that binds the NO 2 group, is 122.04Ǻ, while the C1-N3-C9 angle towards the phenyl that binds the NMe 2 group is 121.07 degrees, resulting in a deviation of 0.97 degrees between the angles.The other central bond angles present slightly higher values towards the plane of the NMe 2 group, as highlighted in the Table 1 (in bold).Furthermore, the systems exhibit deviations of 1.65 Ǻ, 1.09 Ǻ, 0.89 Ǻ, 0.95 Ǻ, and 1.23 Ǻ, for PS-1, PS-3, AB-1, AB-2, and AB-3, respectively.
The molecular structures PS-1 and PS-2 were the only ones found deposited in the Cambridge Crystallographic Data Center (CCDC) database, which reports the single crystal X-ray analysis, recorded with crystallographic information numbers (CIF) 234915 (PS-1) and 1217889 (PS-2) (see Supplementary Materials Table S1).CCDC data allowed us to compare our theoretical results for PS-1 and PS-2 with their experimental analogs.Such a comparison revealed a slight difference between the gas phase and the solid phase, with an emphasis on the central torsion and the torsion caused by the acceptor and donor groups, which can be justified due to the crystalline packing of these compounds, clearly evidencing a significant difference between the molecular and crystalline structure [66][67][68].
Table 2. Spectroscopic parameters regarding the lowest energy transitions retrieved from calculated excitation energies of all compounds in DMSO medium using CAM-B3LYP/6-311++G(2d,p) level of theory.λ max , is the maximum absorption wavelength; f osc is the oscillator strength; µ (D) is the permanent dipole moment.Understanding linear optical effects, like one-photon absorption (OPA), is crucial for determining the transparency range of the samples under investigation.Specifically, we aimed to ascertain if the sample possesses a clear transparency window at approximately 532 nm and beyond.Therefore, before evaluating β HRS , we aimed to explore the characteristics of the vertical transitions linked to the compounds in question.Also, by using TD-DFT calculations, we determined the vertical excitation energies in the DMSO medium, presented in Table 2.It can thus be inferred that all compounds offer a distinct transparency window commencing at ca. 500 nm.This suggests that any ISHG emissions prompted by incident light at 1064 nm, or greater wavelengths, would not be reabsorbed within the medium.Figure 3 illustrates all the vertical transitions associated with OPA simulations using the CAM-B3LYP/6-311++G(2d,p) level of theory.Moreover, a zero or very low probability of vertical transitions around 532 nm will also reduce the probability of two-photon absorption transitions, which could influence the β HRS results, as reported by L. Sciuti et al. [69].The permanent dipole moment (µ) value, also calculated via TD-DFT, can be used, for example, to compare the values of the investigated compounds with reference SHG compounds, such as paranitroaniline (PNA), urea, and potassium dideuterium phosphate (KDP).All investigated compounds presented dipole moments ranging from 10.52 D to 14.03 D, which are greater than the ones calculated for PNA (7.86 D), urea (5.91 D), and KDP (10.28 D) in solvent medium.Among the investigated molecular structures, the lowest µ values belong to PS-2 (~10.5 D) and AB-2 (~10.6 D), while the highest ones belong to PS-1 (~13.9D) and AB-2 (~14.0D).Since the peripheral groups are the same on both molecules, the azomethine functional group (CH=N) present in PS-2 and AB-2 reduces the µ value.On the other hand, PS-1 and AB-2, with an azo functional group (N=N), have higher µ values.PS-3 and AB-3, with a vinyl functional group (CH=CH), have intermediate µ values in this work, with ~12.6 D and ~13.2 D, respectively.The calculated µ values are consistent with the ones reported in the literature [34,49].

Compounds
Finally, all calculated β 0 and β HRS values are presented in Supplementary Materials Tables S3-A-S3-I in the SI for the investigated and reference compounds, PNA, urea, and KDP, under the same computational conditions.The highest β HRS values were obtained using the B3LYP level of theory; this level of theory is well known to overestimate first-order hyperpolarizability calculations [15,70].HRS is not a resonant process, but the contribution of the permanent and transition dipole moments could strongly influence β HRS , according to the semi-classical time-dependent perturbation theory proposed by Orr and Ward [71].According to Orr and Ward, if the transition dipole moment between the ground and the first excited state (µ 01 ) exhibits a higher value and aligns more closely with the incident frequency λ µ 01 → λ HRS , this proximity indicates a scenario nearing resonance absorption.In such cases, it is expected that both β 0 and β HRS will exhibit greater values compared to situations with lower µ 01 and λ µ 01 .Therefore, this can justify β 0 and β HRS overestimations because B3LYP usually enhances oscillator forces and results in a bathochromic shift in their spectral position compared with the experimental values of one-photon absorption of organic compounds [47,72].
Choosing the most suitable level of theory to predict the true value of β HRS is not easy, especially when no experimental results are reported.However, CAM-B3LYP and M06-2X are probably the most suitable levels of theory to estimate the linear and nonlinear optical responses of these types of organic compounds, as previously reported and compared with experimental values [34,40,48,[72][73][74].Moreover, as mentioned in the introduction section, Krawczyk, P. [49] has calculated β 0 values for compounds PS-1, -2, and -3, which corroborate the observation of the tendency for overestimated values of the first-molecular hyperpolarizability when using B3LYP, as shown in Table 3.The β 0 values calculated by Krawczyk, P. show some variation compared to our results, which could be explained by the different approaches, i.e., Møller-Plesset perturbation (MP2), HF methods, BLYP, B3LYP, and the CAM-B3LYP functionals, used to calculate β 0 , as well as the different software QCC package (GAMESS 2023 R2 and Dalton 2.0).For all investigated compounds in our work, both β 0 and β HRS values are indeed overestimated when calculated with B3LYP.
In the photonics field context, particularly in optical communications, the most important results are associated with β HRS .Therefore, we direct the reader to our estimated β HRS results exclusively obtained with CAM-B3LYP and M06-2X, which returned similar values.The M06-2X level of theory, as far as we know, was first used to compute first-order molecular hyperpolarizabilities by Castet et al. [13].Figure 4 depicts β HRS values with CAM-B3LYP/6-311+G(2d,p) used at three distinct incident wavelengths, including the the two most used wavelengths, so far, in optical communications (1310 nm and 1550 nm).
Table 3. Calculated values of static first-order molecular hyperpolarizability of all compounds using the B3LYP and CAM-B3LYP level of theory using the 6-311++G(2d, 2p) basis set with an incident wavelength of 1064 nm.All values are presented in units of 10 −30 cm 4 statvolt −1 both in the gas phase and solvent medium.One can see that the compounds belonging to system type-2 (AB-1, -2, and -3) presented higher values than their analogs from system type-1 (PS-1, -2, and -3), regardless of the incident wavelength.For example, the highest β HRS magnitude difference is between AB-3 and PS-3, in which β HRS has tripled (β HRS ∼ 456).However, it is essential to note that the reported values are still significantly higher than PNA, urea, and KDP.As mentioned before, these are compounds used as references in many experimental works [10,[75][76][77][78] at least at an incident light of 1064 nm.Thus, we used PNA, urea, and KDP to establish a reliable benchmark, allowing us to draw significant conclusions about the NLO activity of the materials under investigation.For example, we obtained values of 0.2 × 10 −30 esu, 0.5 × 10 −30 esu, and 21.5 × 10 −30 esu for KDP, urea, and PNA, which are in good agreement with previous reports found in the literature.Nevertheless, the values obtained for the pseudo-stilbenes and azobenzenes investigated here are up to sixty-five, thirty-seven, and thirty times higher than PNA and hundreds of times greater than urea and KDP at the wavelengths of 1064 nm, 1330 nm, and 1550 nm, respectively.The results were compared under the same conditions of solvent medium (DMSO) and level of theory.The difference in the β HRS magnitude between the investigated compounds and the reference ones reflects the increased π-conjugation and molecular size in pseudo-stilbene and azobenzene compared to the reference compounds.The observed differences in β HRS values reflect the underlying fundamental principles of molecular structure and electronic distribution, affirming the validity and relevance of the investigated compounds as OFCs.
Additionally, it was found that the β HRS of the compounds studied is strongly influenced by the central torsion angle and the torsion angle close to the -NMe 2 substituent group, and is potentially improved by the presence of CH=CH and N=N functional groups and increased conjugation, as indicated in Table 4. Furthermore, it is observed that the largest torsion angles were noted for the compounds PS-2 and AB-2, corroborating the lowest values of dipole moment and first-order molecular hyperpolarizability, respectively.However, compound AB-3, despite having good planarity and the highest first hyperpolarizability value, exhibits a decrease in the dipole moment value.This suggests that the twist near the -NMe 2 group may have had a significant influence.The properties of NLO depend on several structural factors, such as molecular planarity, charge transfer between side groups, conjugation length, chemical environment, and the effect of intermolecular interactions [27].Compound AB-2 (system type-I) shows an increased conjugation length with intermolecular interactions C-H• • • π.However, it presents a non-planar conformation with a twisted structure, significantly reducing µ and β HRS values compared to the other molecules of system II.In addition, structural modifications to system type-II seem to improve the optical activity compared to system type-I, revealing the higher potential of this system for OFC applications.

Conclusions
This comprehensive investigation into the linear and nonlinear optical properties of two novel organic molecular systems (type-I and type-II) provides valuable insight into the intricacies of push-pull systems, mainly how they affect the β HRS magnitude.The correlation between molecular structure and β HRS response was evident, particularly the significance of the azomethine group (CH=N) in affecting the planarity and optical properties.The presence of this group led to reduced molecular planarity in some compounds, which in turn influenced the overall β HRS response.Predicted β HRS values, which are particularly important for photonics and optical communications applications, were significantly more significant for the type-II system (AB-1, -2, and -3) than their type-I counterparts (PS-1, -2, and -3).Especially when considering optical communication wavelengths of 1310 nm and 1550 nm, the type-II system consistently outperformed the type-I system.The compounds studied in this work exhibited β HRS values that substantially surpassed those of reference compounds PNA, urea, and KDP.This comparative analysis emphasized the potential of these molecules for applications requiring robust NLO properties.Moreover, it shows that the novel molecular design (type-II) presented in this work has the potential to be used for photonic crystals with the ability for nonlinear optical conversion.
However, there are significant challenges to predicting β HRS , such as the overestimation and underestimation of β HRS when utilizing different levels of theory.For example, using the functional B3LYP level of theory led to an overestimation of the magnitude compared with the other two functionals (CAM-B3LYP and M06-2X).Choosing an accurate level of theory to predict the true value of β HRS proved complex, especially when no experimental data were available for comparison.As far as we know, there is no reported work with experimental β HRS values measured via the HRS or EFISH techniques.Therefore, the comparison between experimental and theoretical β HRS was not possible.In this case, the absence of experimental results will not be thought of as a drawback of our work.On the contrary, it confirms the novelty of this current work.Moreover, such novelty could encourage chemical research groups to synthesize compounds of system type-II for future experiments, which will verify the true potential of these organic compounds in OFC devices.Organic OFCs with high responses have the potential to revolutionize communication and data transmission due to their potential to offer incredibly high speeds, greater bandwidths, and lower signal losses.

Figure 1 .
Figure 1.The molecular structures investigated in this work are pseudo-stilbenes (type I) and azobenzene-based compounds (type II).

Figure 2 .
Figure 2. Graphical representations regarding the views of all investigated molecular structures in (a) front view, (b) planar view, and (c) side view.All results were obtained after geometry optimization using the CAM-B3LYP/6-311++G(2d, p) level of theory.The colors in the illustration represent the atoms within the compounds, where carbon is depicted in gray, hydrogen in white, nitrogen in blue, and oxygen in red.

Figure 3 .
Figure 3. TD-DFT calculations were used to obtain the vertical transitions for system type I and II molecules using the CAM-B3LYP/6-311++G(2d,p) level of theory in solvent medium (DMSO).
The values were converted from atomic units to electrostatic unities using 1 a.u.= 8.641 × 10 −33 esu. a this work, using the Gaussian 16 (rev.B.01) software package; b values in parentheses were taken from ref.[49] using GAMESS and Dalton 2.0 with the 6-311++G(d,p) basis set.

Table 1 .
Bond length, bond angle, and torsion angle of all compounds investigated in this work after geometry optimization in solvent medium (DMSO).

Table 4 .
Relationship between geometric and optical parameters for all molecules using CAM-B3LYP/6-311+G(2d,p) level of theory.