# Prediction of the Linear and Nonlinear Optical Properties of a Schiff Base Derivatives via DFT

Authors: Clodoaldo Valverde; Ítalo Nuta Ribeiro; João Victor B. Soares; Basílio Baseia; Francisco A. P. Osório
Journal: Advances in Condensed Matter Physics
 (2019)

---

## Abstract

In this work, the density functional theory (DFT) calculation combined with a polarizable continuum model (PCM) was used to study the solvent media effects on the electrical and geometrical behaviors of the Schiff-base derivative, (E)-4-[({4-[(pyridin-2-ilmetilideno)amino]phenyl}amino)-metil]fenol (EPAF). The linear and nonlinear optics parameters, as the dipole moment, linear polarizability, and first and second hyperpolarizabilities, were calculated at DFT/B3LYP/6-311+G(d) level, for the EPAF molecule in several solvent media. The dynamic behavior of the Hyper-Rayleigh Scattering (HRS) first hyperpolarizability was studied as function of the electric field frequency. The results presented for HRS first hyperpolarizability suggest that the studied crystal has good nonlinear optical properties. In addition, the gap energy was calculated from the HOMO-LUMO energies difference in several solvent media. The EPAF crystal intermolecular interactions were studied by the Hirshfeld surface analysis. The third-order electric susceptibility χ(3) of the crystal EPAF was also calculated, indicating the EPAF crystal as a promising candidate for NLO applications in photonic and optoelectronic devices.


---

## Body

## 1. Introduction

In recent years, the use of organic crystals as nonlinear optical (NLO) materials has been growing motivated by the easy manipulation of these crystals, which allow controlling the material NLO properties [1]. Studies of nonlinear optical processes contribute significantly to the development of photonics [2, 3], spectroscopy [4, 5], fiber optic lines [4], optical switches [6], frequency converters [7], electrooptic modulators [8], and data transmission network, among others [9] and still in numerous applications in the medical and pharmacological sectors [10–12]. The compounds that have a high nonlinearity are of great interest for the field of nonlinear optics, since they make up the manufacture of devices that operate with high speed [9, 13].Schiff bases are aldehyde or ketone-like compounds in which the carbonyl group is replaced by animine or azomethine group. They are widely used for industrial purposes and also exhibit a broad range of biological activities, as antibacterial, anticancer, anti-inflammatory, and antitoxic properties as stated by Lozier et al. [14]. The compounds of this group have a potential for use in optical memory device, because they form coordinating or grid polymers. Recently, a Schiff base derivative (E)-4-[({4-[(pyridin-2-ilmetilideno)amino]phenyl}amino)-metil]fenol (EPAF) with molecular formula C19H17N3O has been synthesized, crystallized, and structurally characterized by Faizi et al. [15]. Also, these authors used the density functional theory (DFT) at B3LYP/6-311G(d,p) level to calculate geometrical parameters of the EPAF single molecule and compare them with the obtained X-ray results and verified a good agreement with the experimental data.In this work, we study the effect of several solvent media on geometry and the electrical parameters of the EPAF molecule. The EPAF geometry optimization was performed in the gas phase and in various solvent media using DFT at level CAM-B3LYP/6-311+G(d). The overlap between the X-ray data for EPAF molecules with the DFT results in several solvent media has shown a big variation of the torsion angles between the rings. The solvent media effects on EPAF molecule static electrical parameters as dipole moment, linear polarizability, first hyperpolarizability, and the second hyperpolarizability are studied. The behaviors of the EPAF Hyper-Rayleigh Scattering (HRS) first hyperpolarizability in several solvent media as function of the electric field frequencies and of the static dielectric constant value were analyzed. Also the energies of the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO), as the gap energies in several solvent media, have been calculated. Also a Hirshfeld surface analysis is presented to complement the EPAF crystal study of [15].


## 2. Methodology

The EPAF compound was crystallized in a monoclinic space group P21/c [15] (see Figure 1) with the following crystallographic data: a = 10.55652 (7), b = 7.9136 (6), c = 20.8153 (13), α = 90°, γ = 90°, and β = 118.408°, unit cell volume 1530.77 Å, with four molecules in the unit cell (Z = 4). The EPAF molecule is nonplanar, with the phenolic ((a)) and pyridine rings ((c)) being inclined to the central benzene ring ((b)).Figure 1
A view of the EPAF molecule showing the atom labeling scheme.


### 2.1. Hirshfeld Surface

Hirshfeld surface (HS) analysis serves as a powerful tool for obtaining crucial information about the intermolecular interaction of molecular crystals. The size and shape of the Hirshfeld surface allow the visualization and investigation of both qualitative and quantitative intermolecular crystal bonds. Thus, the HS surface is obtained and two distances are defined, the distance from the point to the nearest atom off the surface (de) and the distance to the nearest atom within the surface (di). The identification of the regions of particular importance to intermolecular interactions is obtained by mapping normalized contact distance (dnorm), defined by (1)dnorm=di-rivdWrivdW+de-revdWrevdW,where rivdW and revdW are the van der Waals radii of the atoms. Studies of the π…π interactions were performed analyzing the surface shape index, where this interaction type can be identified through the red and blue triangle, one facing the other [18].


### 2.2. Computational Details

The solvent media effects on the electric parameters of the EPAF molecule were studied employing the polarizable continuum model (PCM). First the geometry optimization calculations were performed in gas phase and after in nineteen solvent media, both method using the DFT at B3LYP/6-311+G(d) level.In the present study, the total dipole moment, average linear polarizabilityα, and anisotropy of the linear polarizability ∆α of the title compounds have been calculated using the following expression:(2)μ=μx2+μy2+μz21/2,(3)α=αxx+αyy+αzz3,(4)∆α=2-1/2αxx-αyy2+αyy-αzz2+αzz-αxx2+6αxz2+αxy2+αyz21/2.The total and HRS (Hyper Rayleigh Scattering) molecular first hyperpolarizabilities are given by(5)βtotal=βxxx+βxyy+βxzz2+βyyy+βyxx+βyzz2+βzxx+βzyy+βzzz21/2,and(6)βHRS=βZZZ2+βXZZ2,whereX-direction is assumed as the fundamental light beam propagation and polarized in theZ-direction and βZZZ2 and βXZZ2 are macroscopic averages calculated from the first hyperpolarizability components (βijk)[19, 20] through the expressions. Hyper-Rayleigh Scattering (HRS) components (βijk) [19, 20] were calculated through the expressions(7)βZZZ2=121030δ1+12δ2+δ3+δ5+6δ4+δ6+2δ7+δ8+δ11+4δ9+δ10,βXZZ2=12106δ1-δ3-δ5+δ7+8δ2+18δ4+4δ6-δ8-2δ9+3δ10-δ11,in which the coefficients δn are defined in Table 1. In this case, we adopted the laboratory system of reference by theX,Y, andZ coordinates, and the molecular system of reference by thex,y, andz coordinates.Table 1
The HRS first hyperpolarizabilityδn coefficients.
δ 1 =  ∑ i β i i i 2 , δ 7 =  ∑ i , j , k β i j j β i k k ,  δ 2 =  ∑ i , j β i i i β i j j, δ 8 =  ∑ i , j , k β j i j + β j j i β k i k + β k k i ,  δ 3 =  ∑ i , j β i i i β j i j + β j j i , δ 9 =  ∑ i , j , k β i j j β k i k + β k k i ,  δ 4 =  ∑ i , j β i j j 2, δ 10 =  ∑ i , j , k β i j k + β i k j 2 ,  δ 5 =  ∑ i , j β i j j β j i j + β j j i , δ 11 =  ∑ i , j , k β i j k + β i k j β j j k + β j k i .  δ 6 =  ∑ i , j β j i j + β j j i 2The average second molecular hyperpolarizability is given by(8)γ=115∑i,j=x,y,zγiijj+γijij+γijji.Using the Kleymann symmetry, theγ-value can be calculated through the following expression:(9)γ=15γxxxx+γyyyy+γzzzz+2γxxyy+γxxzz+γyyzz.All computational calculations related to the linear and nonlinear electric parameters of the compounds were performed in theGAUSSIAN09 program [21].


### 2.3. Frontiers Molecular Orbital

In order to verify the EPAF molecule stability in the several solvent media, the highest occupied orbital energy (HOMO), which has an electron-donor character, and the lowest unoccupied molecular orbital energy (LUMO), which has an electron-acceptor character, were calculated. The molecule ability to donate or receive electrons is greater as the HOMO-value or the LUMO-value, respectively [22]. The energy difference between HOMO and LUMO gives the GAP energy, a parameter which is directly related to the stability of the compound [23].


### 2.4. Solvent Media

The solvent medium polarity is related to it solvability. Defining the polarity concept quantitatively is very difficult, and controversial; in this work, we will use the concept of polarity defined by the scaleETN (normalized transition energy) of Dimroth and Reichardt [16]. The ETN-value is based on the transition energy for the longest-wavelength solvatochromic absorption band of the pyridinium N-phenolate betaine dye (see Table 2). Here, solvent media with static dielectric constant (ε) value smaller than 5 will be considered as nonpolar.Table 2
Solvent media polarity parameter and static dielectric constant, as in ref. [16].
Solvent medium E T N ε Water 1.00 78.355 protic Formamide 0.775 108.940 protic Methanol 0.762 32.613 protic FormicAcid 0.728 51.100 protic n-MethylFormamide-mixture 0.722 181.560 protic Ethanol 0.654 24.852 protic 1-Butanol 0.586 17.332 protic Acetonitrile 0.460 35.688 aprotic DMSO 0.444 46.826 aprotic 2-Methyl-2-Propanol 0.389 12.470 protic Acetone 0.355 20.493 aprotic DiChloroEthane 0.327 10.125 aprotic DichloroMethane 0.309 8.930 aprotic Chloroform 0.259 4.711 nonpolar Tetrahydrofuran 0.207 7.426 aprotic ChloroBezene 0.188 5.697 aprotic Toluene 0.099 2.374 nonpolar Heptane 0.012 1.911 nonpolarAll computational calculations related to the linear and nonlinear electric parameters of the compounds were performed in theGAUSSIAN09 program [21].


## 2.1. Hirshfeld Surface

Hirshfeld surface (HS) analysis serves as a powerful tool for obtaining crucial information about the intermolecular interaction of molecular crystals. The size and shape of the Hirshfeld surface allow the visualization and investigation of both qualitative and quantitative intermolecular crystal bonds. Thus, the HS surface is obtained and two distances are defined, the distance from the point to the nearest atom off the surface (de) and the distance to the nearest atom within the surface (di). The identification of the regions of particular importance to intermolecular interactions is obtained by mapping normalized contact distance (dnorm), defined by (1)dnorm=di-rivdWrivdW+de-revdWrevdW,where rivdW and revdW are the van der Waals radii of the atoms. Studies of the π…π interactions were performed analyzing the surface shape index, where this interaction type can be identified through the red and blue triangle, one facing the other [18].


## 2.2. Computational Details

The solvent media effects on the electric parameters of the EPAF molecule were studied employing the polarizable continuum model (PCM). First the geometry optimization calculations were performed in gas phase and after in nineteen solvent media, both method using the DFT at B3LYP/6-311+G(d) level.In the present study, the total dipole moment, average linear polarizabilityα, and anisotropy of the linear polarizability ∆α of the title compounds have been calculated using the following expression:(2)μ=μx2+μy2+μz21/2,(3)α=αxx+αyy+αzz3,(4)∆α=2-1/2αxx-αyy2+αyy-αzz2+αzz-αxx2+6αxz2+αxy2+αyz21/2.The total and HRS (Hyper Rayleigh Scattering) molecular first hyperpolarizabilities are given by(5)βtotal=βxxx+βxyy+βxzz2+βyyy+βyxx+βyzz2+βzxx+βzyy+βzzz21/2,and(6)βHRS=βZZZ2+βXZZ2,whereX-direction is assumed as the fundamental light beam propagation and polarized in theZ-direction and βZZZ2 and βXZZ2 are macroscopic averages calculated from the first hyperpolarizability components (βijk)[19, 20] through the expressions. Hyper-Rayleigh Scattering (HRS) components (βijk) [19, 20] were calculated through the expressions(7)βZZZ2=121030δ1+12δ2+δ3+δ5+6δ4+δ6+2δ7+δ8+δ11+4δ9+δ10,βXZZ2=12106δ1-δ3-δ5+δ7+8δ2+18δ4+4δ6-δ8-2δ9+3δ10-δ11,in which the coefficients δn are defined in Table 1. In this case, we adopted the laboratory system of reference by theX,Y, andZ coordinates, and the molecular system of reference by thex,y, andz coordinates.Table 1
The HRS first hyperpolarizabilityδn coefficients.
δ 1 =  ∑ i β i i i 2 , δ 7 =  ∑ i , j , k β i j j β i k k ,  δ 2 =  ∑ i , j β i i i β i j j, δ 8 =  ∑ i , j , k β j i j + β j j i β k i k + β k k i ,  δ 3 =  ∑ i , j β i i i β j i j + β j j i , δ 9 =  ∑ i , j , k β i j j β k i k + β k k i ,  δ 4 =  ∑ i , j β i j j 2, δ 10 =  ∑ i , j , k β i j k + β i k j 2 ,  δ 5 =  ∑ i , j β i j j β j i j + β j j i , δ 11 =  ∑ i , j , k β i j k + β i k j β j j k + β j k i .  δ 6 =  ∑ i , j β j i j + β j j i 2The average second molecular hyperpolarizability is given by(8)γ=115∑i,j=x,y,zγiijj+γijij+γijji.Using the Kleymann symmetry, theγ-value can be calculated through the following expression:(9)γ=15γxxxx+γyyyy+γzzzz+2γxxyy+γxxzz+γyyzz.All computational calculations related to the linear and nonlinear electric parameters of the compounds were performed in theGAUSSIAN09 program [21].


## 2.3. Frontiers Molecular Orbital

In order to verify the EPAF molecule stability in the several solvent media, the highest occupied orbital energy (HOMO), which has an electron-donor character, and the lowest unoccupied molecular orbital energy (LUMO), which has an electron-acceptor character, were calculated. The molecule ability to donate or receive electrons is greater as the HOMO-value or the LUMO-value, respectively [22]. The energy difference between HOMO and LUMO gives the GAP energy, a parameter which is directly related to the stability of the compound [23].


## 2.4. Solvent Media

The solvent medium polarity is related to it solvability. Defining the polarity concept quantitatively is very difficult, and controversial; in this work, we will use the concept of polarity defined by the scaleETN (normalized transition energy) of Dimroth and Reichardt [16]. The ETN-value is based on the transition energy for the longest-wavelength solvatochromic absorption band of the pyridinium N-phenolate betaine dye (see Table 2). Here, solvent media with static dielectric constant (ε) value smaller than 5 will be considered as nonpolar.Table 2
Solvent media polarity parameter and static dielectric constant, as in ref. [16].
Solvent medium E T N ε Water 1.00 78.355 protic Formamide 0.775 108.940 protic Methanol 0.762 32.613 protic FormicAcid 0.728 51.100 protic n-MethylFormamide-mixture 0.722 181.560 protic Ethanol 0.654 24.852 protic 1-Butanol 0.586 17.332 protic Acetonitrile 0.460 35.688 aprotic DMSO 0.444 46.826 aprotic 2-Methyl-2-Propanol 0.389 12.470 protic Acetone 0.355 20.493 aprotic DiChloroEthane 0.327 10.125 aprotic DichloroMethane 0.309 8.930 aprotic Chloroform 0.259 4.711 nonpolar Tetrahydrofuran 0.207 7.426 aprotic ChloroBezene 0.188 5.697 aprotic Toluene 0.099 2.374 nonpolar Heptane 0.012 1.911 nonpolarAll computational calculations related to the linear and nonlinear electric parameters of the compounds were performed in theGAUSSIAN09 program [21].


## 3. Results and Discussions

### 3.1. Structural Commentary

The optimized geometry in gas phase and in several solvent media of the EPAF structure was analyzed through the root mean square deviation (RMSD) of the overlap between the molecular geometry determined by X-ray and the theoretical results obtained in the presence of solvent media; the H-atoms were disregarded in view of their uncertainties in X-ray position refinement. All the optimized geometries of the EPAF molecule calculated in various solvent media are in the Supplementary Materials (see TablesS1 – S20). The optimized geometry results in chloroform using the ring (a) as anchorage are presented in Figure 2; the RMSD is 0.5460 a.u. Table 3 shows the RMSD for the optimized geometry in several solvent media; as can be observed, the RMSD-values for nonpolar solvent media (argon, heptane, and toluene) are approximately 0.50 a.u. and for the others it is 0.54 a.u.Table 3
Comparison of selected geometric data for EPAF (°) from calculated DFT-Solvent and X-ray data.
ε RMSD C6-C1-O1 C8-N1-C7 N1-C7-C4 C14-N2-C11 N2-C14-C15 N3-C15-C14 EPAF-X-Ray   120.3 123.4 112.3 121.5 122.2 115.9 Gas-phase 1.000 0.498 122.9 122.7 110.8 121.5 122.1 115.4 Argon 1.430 0.502 122.8 122.7 110.8 121.6 122.1 115.4 Heptane 1.911 0.502 122.8 122.8 110.8 121.6 122.2 115.4 Toluene 2.374 0.501 122.8 122.8 110.8 121.6 122.2 115.4 Chloroform 4.711 0.546 122.8 122.7 110.8 121.6 122.3 115.3 ChloroBezene 5.697 0.542 122.8 122.7 110.8 121.6 122.3 115.3 Tetrahydrofuran 7.426 0.544 122.8 122.7 110.7 121.6 122.4 115.3 DichloroMethane 8.930 0.544 122.7 122.8 110.7 121.6 122.4 115.3 DiChloroEthane 10.125 0.544 122.7 122.8 110.7 121.6 122.4 115.3 2-Methyl-2-Propanol 12.470 0.543 122.7 122.7 110.7 121.6 122.4 115.3 1-Butanol 17.332 0.543 122.7 122.7 110.7 121.6 122.4 115.3 Acetone 20.493 0.543 122.7 122.7 110.7 121.6 122.4 115.3 Ethanol 24.852 0.543 122.7 122.7 110.7 121.6 122.4 115.3 Methanol 32.613 0.543 122.7 122.7 110.7 121.6 122.4 115.3 Acetonitrile 35.688 0.543 122.7 122.7 110.7 121.7 122.4 115.3 DMSO 46.826 0.543 122.7 122.7 110.7 121.7 122.4 115.2 FormicAcid 51.100 0.543 122.7 122.7 110.7 121.7 122.4 115.2 Water 78.355 0.543 122.7 122.7 110.7 121.7 122.4 115.2 Formamide 108.940 0.543 122.7 122.7 110.7 121.7 122.4 115.2 n-MethylFormamide-mixture 181.560 0.542 122.7 122.7 110.7 121.7 122.4 115.2Figure 2
Overlap between the molecular structure determined by X-ray (in red) and in Chloroform (in blue). The ring (a) was used as anchorage.The solvent media presence causes a significant deviation of the RMSD parameter between the X-ray geometry data and the theoretical results, as can be seen in Figure2 for the chloroform (RMSD=0.544). The C4-C7-N1-C8 and C11-N2-C14-C15 torsion angles change due to the solvent medium effects (chloroform) from -166.3° (X-ray data) to +176.24 and +176.4 (X-ray data) to -177.72°, respectively. This effect occurs due to the negative charge transfer to the hydroxyl bond to the terminal phenolic ring (C1–C6). Table S21 (Supplementary Materials) shows the charges of the atoms of the molecule under the effect of several solvent media using the PCM method calculated via DFT (CAM-B3LYP/6-311+G(d)) using the electrostatic model ChelpG. The ChelpG fit of the pyridine ring and hydroxyl charges are 0.04e and -0.238e both in chloroform medium (see the ChelpG charges for all atoms in Table S21 of the Supplementary Materials).DFT results for the EPAF optimized geometry in gas phase and in several solvent media are shown in Table3. The C7—N1—C8 angle presents a reduction of 0.59% for both polar and nonpolar solvent. In the solvent media, the C6—C1—O1 angle increased around 2.0% and the N1—C7—C4 angle remains invariable in gas phase, toluene, heptane, and argon. This angle has increased 0.08% and 0.16% and in chloroform and in the polar solvent medium, respectively. Also the C14—N2—C11 angle increased 0.08% and 0.16% for solvent media with ε < 37 and ε > 37, respectively. The X-ray and theoretical results in the solvent media for the angles, N2—C14—C15 and N3—C15—C14, are practically the same.


### 3.2. HOMO-LUMO

In this section, the HOMO and LUMO results for the EPAF molecule in several solvent media are presented. The HOMO and LUMO orbital and the gap energy values in the solvent medium n-methyl formamide mixture and in gas phase are shown in Figure3. All HOMO-LUMO figures are in the Supplementary Materials (see Figures S1 – S18).Figure 3
HOMO and LUMO in n- methyl formamide mixture.As can be seen in Figure4, the gap energy decreases with the increasing of the dielectric constant value of the solvent medium. The gap energy EG goes from 3.466 eV (ε=1.0) to 3.349 eV (ε=182.4) indicating that the maximum absorption wavelengths (λmax=hc/EG) fall into the ultra violet (UV) region. The gap energy is an important parameter that determines the chemical properties of a molecule as the kinetic stability, chemical reactivity, optical polarizability, etc. Larger values of the gap energy are directly related to greater stability of the molecule.Figure 4
Gap energy as function of the dielectric constant values.All values of HOMO-LUMO and gap energy of the EPAF molecule calculated in various solvent media are in the Supplementary Materials (see TablesS22).


### 3.3. Nonlinear Optical Properties in Solvent Medium

Table4 shows the DFT/B3LYP/6-311+G(d) static results for the total dipole moment (μtot), the isotropic linear polarizability α(0;0), the first hyperpolarizability parallel to the dipole moment (β||z(0;0,0)), and the second hyperpolarizability γ(0;0,0,0) both in gas-phase and in several solvent media. As can be noted, the values of μ, α(0;0) and γ(0;0,0,0) increase with the increasing of the ε-values (Figures 5(a), 5(b), and 5(d)). In all solvent media, the β||z(0;0,0)-values are negative, and for ε≥5 they oscillate around 0.6×10-30esu; however, the β||z(0;0,0) increases with the increasing of the ε-values, as can be seen in Figure 5(c).Table 4
DFT static results for the EPAF electrical parameters.
ε γ ( 0 ; 0,0 , 0 ) β | | z ( 0 ; 0,0 ) α ( 0 ; 0 ) μ (D) 10−36esu 10−30esu 10−24esu 1.000 170 -0.230 42 2.86 1.430 220 -0.200 44 3.01 1.911 254 -0.157 46.4 3.12 2.374 280 -0.127 47.7 3.19 4.711 351 -0.611 51 3.44 5.697 367 -0.581 51.7 3.48 7.426 387 -0.602 52.6 3.53 8.930 399 -0.602 53.1 3.57 10.125 406 -0.603 53.4 3.59 12.470 416 -0.595 53.9 3.61 17.332 435 -0.605 54.7 3.67 20.493 440 -0.607 54.9 3.68 24.852 446 -0.609 55.2 3.70 32.613 448 -0.609 55.3 3.70 35.688 452 -0.609 55.5 3.72 46.826 454 -0.608 55.6 3.72 51.100 458 -0.612 55.6 3.73 78.355 461 -0.604 55.9 3.74 108.940 464 -0.604 56.0 3.75Figure 5
DFT static results for the EPAF electric parameter as function of theε-value.
(a) (b) (c) (d)From Table4, it can be verified that the first and second hyperpolarizabilities present the greater variation due to the solvent medium presence, from heptane (ε=1.911) to n-methyl formamide mixture (ε=181.56); the values of β||z(0;0,0) and γ(0;0,0,0) present a percentage increasing of 285% and 83% (see Figures 5(c) and 5(d) ). However, the character, protic or aprotic, seems to not influence the results, but the nonpolar solvent media present a smaller value of the total dipole moment, and this can be understood taking into account that the polar solvent media the OH present greater increasing of the negative charge, and similar effect can be verified for other electric parameters (Figure 5).The DFT static results for theβHRS-values obtained from (6) in a solvent medium are several orders of magnitude higher than the values of β||z(0;0,0). While as shown in Table 3 for the solvent media the absolute values of β||z(0;0,0) go from 0.2 to 0.6 (in units of 10-30esu), βHRS(0,0,0) go from 24.4 to 53.0 (in units of 10-30esu). Table S23 of the Supplementary Material shows the static and dynamic βHRS-values for the gas phase and all the solvent media. The behaviors of the static and dynamic HRS first hyperpolarizability in the solvent media are similar; therefore, in Figure 6, the results for βHRS as function of the field frequency (ω) are shown for only four solvents, namely, heptane (ε=1.911, nonpolar), chloroform (ε=4.71; larger RMSD parameter), DMSO (ε=46.8; aprotic), and water ( ε=78.4; ETN=1). From Figure 6 and from Table S23, two resonant regions can be identified in the region 0.06a.u.<ω<0.09a.u. So we will work away from the transition region, because in this region it would be complicated for the experiment to work due to deleterious effects. We will study βHRS on the frequency of ω=0.0428 a.u.(1064 nm).Figure 6
Dynamic HR Scattering first hyperpolarizability as function of the field frequency (ω) in a.u..
(a) (b) (c) (d)We can also highlight that the protic or aprotic character and also the polarity of the solvent did not alter the HRS first hyperpolarizability behavior as a function of the frequency. TheβHRS-values increase with the increasing of the dielectric constant value (or with the decreasing of the gap energy) as can be seen in Figure 7. This effect can be explained by the increase of the hydroxyl group charge, which acts as an electrons receptor of the benzene ring, decreasing the electron density and the resonant stability, making it difficult to the electrophilic attack. In this case, the OH group acts as a disabling or metamanager. The values of βHRS in chloroform are 3 to 5 times higher than the 2,4,6-tris(benzylamino)-1,3,5-triazine [24] derivatives also measured in chloroform.Figure 7
Static and dynamic (1064nm) HRS first hyperpolarizability as function of theε-value.


### 3.4. Hirshfeld Surface Analyses

The molecular structure of (E)-4-[({4-[(piridin-2-ilmetilideno) amino]fenil}amino)-metil]fenol (EPAF), C19H17N3O, shows H bonds of the type π- π, C-H... π, C-H..C, O-H...N e N-H...O which contribute to the structure stability. Figure 8 shows the normalized Hirshfeld surface. Figure 8(a) shows the interactions O1-H1...N3 (D...A = 2.796 Å, D-H...A = 178.75°) and C14-H14...C1 (D...A = 3.640 Å, D-H...A = 146.17°). Figure 8(b) identifies the interactions N1-H1A...O1 (D...A = 2.982 Å, D-H...A = 170.07°), O1-H1...N3 (D...A = 2.796 Å, D-H...A = 178.75°), and C14-H14...C1 (D...A = 3.640 Å, D-H...A = 146.17°). Figure 8(c) shows the interactions C17-H17...C8 (D...A = 3,655 Å, D-H...A = 144.90°).Figure 8
View of EPAFdnorm surface showing the interactions O1-H1...N3 e C14-H14...C1 (a), N1-H1A...O, O1-H1...N3, and C14-H14...C1 (b), and C17-H17...C8 (c).Figure9(a) shows the interaction C17-H17...C8 (D...A = 3.655 Å, D-H...A = 144.90°), the same shown in Figure 9(c). The interaction N1-H1A...O1 (D...A = 2.982 Å, D-H...A = 170.07°) is shown in Figure 9(b) and the interactions C7-H7A...C9 (D...A = 3.723 Å, D-H...A = 151.80°) and C7-H7A...C10 (D...A = 3.681Å, D-H...A = 156.93°) are shown in Figure 9(c).Figure 9
View of EPAFdnorm surface showing the interactions (a) O1-H1...N3 and C14-H14...C1, (b) N1-H1A...O, O1-H1...N3, and C14-H14...C1, and (c) C17-H17...C8.Figure10 shows the shape index surface for the EPAF, where it is possible to see two triangles, one blue and the other red, one facing the other, representing the π- π interaction (π...π = 3,779 Å).Figure 10
View of the shape index surface for EPAF showing the highlightedπ...π interaction.The fingerprints plots were used to analyze the amount of intermolecular interactions according to the nature and percentage of the interaction. Figures11(a)–11(g) show the fingerprints for the EPAF; as we can see, the interaction with greater percentage in the whole crystal is the H..H with 45.8% of the interactions followed by the C…H with 34.2% of the interactions.Figure 11
Fingerprint showing the 2D interactions of EPAF crystal.
(a)
C...C 2.1% (b)
C...H 34.2% (c)
H...H 45.8% (d)
N...C 0.9% (e)
N...H 10.6% (f)
O...C 0.1% (g)
O...H 6.4% (h)


### 3.5. Nonlinear Optical Properties of the EPAF Crystal

The supermolecular (SM) method was used to simulate the EPAF crystal; details of the SM approach are given in the work of C. Valverde et al. [25]. To employ the SM approach, we have used the (x-ray) experimental geometry of the asymmetric unit of the EPAF. The packaging effects of the EPAF were modeled by constructing a bulk with the unit cells in a 9×9×9 configuration. Each unit cell contains four asymmetric units, totaling 2916 molecules in the monoclinic system, each atom surrounding the EPAF molecule (blue) being treated as a point charge (see scheme in Figure 12).Figure 12
Scheme of the bulk representing the embedded molecule.For the calculation, the average linear polarizabilityα can be related to the linear refractive index (n) of the crystal via the Clausius-Mossotti relation, given by [26]: (10)n2-1n2+2=4πN3α,whereN stands for the number of molecules per unit cell volume. The experimental quantity, the third-order electric susceptibility χ(3), is related to the second hyperpolarizability by (11)χ3=f4NγϵoV,where N is the number of molecules per unit cell volume (V) andf is the Lorentz local field correction factor given by(12)f=n2+23.The EPAF crystal presents a rapid convergence in the dipole moment through the SM approach [25] (see Figure 13).Figure 13
Evolution of values of the dipole moment of the EPAF crystal with the respective iteration numbers.The quantum molecular calculations were performed with the Gaussian 09 program package [21]. For the calculation of χ3-ω;ω,ω,-ω, we used a small frequency [27, 28], estimate of the second frequency-dependent hyperpolarizability (γ-ω;ω,ω,-ω≅2γ-ω;ω,0,0-γ0;0,0,0) associated to a nonlinear optical process [29] of the intensity dependent refractive index (IDRI) from dc-K results. Table 5 shows the values of the EPAF when we consider the environment via the SM.Table 5
The values of the linear polarizability (in10-24 esu) were calculated CAM-B3LYP/6-311++G(d,p) and for values of the second hyperpolarizability in (in 10-36 esu) we used CAM-B3LYP/6-311++G(d,p) for EPAF embedded molecule (λ=532nm).
α ( - ω , ω ) γ ( 0 ; 0,0 , 0 ) γ ( - ω ; ω , 0,0 ) γ - ω ; ω , ω , - ω 46.50 174.00 628.195 1082.39Table6 presents the DFT prediction and the experimental results of the macroscopic quantities studied. The values of EPAF crystal are therefore 24.51, 2953, 3413, and 2866 times higher, respectively, than the values found experimentally by Prabhu et al. [17] (see Table 6).Table 6
DFT/CAM-B3LYP/ 6-311++G(d,p) results for the linear refractive index and third-order nonlinear susceptibility (10-22m2/V2) for the case dynamic (λ=532nm) of the EPAF crystal.
Sample n ( ω ) χ 3 - ω ; ω , ω , - ω EPAF (this work) 2.02 6793.08 (2E)-3-(3-methylphenyl)-1-(4-nitrophenyl)prop-2-en-1-one (3MPNP) [17] 1.418 277.1 (2E)-1-(4-bromophenyl)-3-[4-methylsulfanyl) phenyl]prop-2-en-1-one (4Br4MSP) [17] 1.363 2.30 (2E)-1-(3-bromophenyl)-3-[4(methylsulfanyl) phenyl]prop-2-en-1-one (3Br4MSP) [17] 1.365 1.99 (2E)-3[4(methylsulfanyl)phenyl]-1-(4-nitrophenyl)prop-2-en-1-one (4N4MSP) [17] 1.360 2.37The high value of the third-order electric susceptibilityχ(3) indicates the crystal as a promising candidate for NLO applications in photonic and optoelectronic devices.


## 3.1. Structural Commentary

The optimized geometry in gas phase and in several solvent media of the EPAF structure was analyzed through the root mean square deviation (RMSD) of the overlap between the molecular geometry determined by X-ray and the theoretical results obtained in the presence of solvent media; the H-atoms were disregarded in view of their uncertainties in X-ray position refinement. All the optimized geometries of the EPAF molecule calculated in various solvent media are in the Supplementary Materials (see TablesS1 – S20). The optimized geometry results in chloroform using the ring (a) as anchorage are presented in Figure 2; the RMSD is 0.5460 a.u. Table 3 shows the RMSD for the optimized geometry in several solvent media; as can be observed, the RMSD-values for nonpolar solvent media (argon, heptane, and toluene) are approximately 0.50 a.u. and for the others it is 0.54 a.u.Table 3
Comparison of selected geometric data for EPAF (°) from calculated DFT-Solvent and X-ray data.
ε RMSD C6-C1-O1 C8-N1-C7 N1-C7-C4 C14-N2-C11 N2-C14-C15 N3-C15-C14 EPAF-X-Ray   120.3 123.4 112.3 121.5 122.2 115.9 Gas-phase 1.000 0.498 122.9 122.7 110.8 121.5 122.1 115.4 Argon 1.430 0.502 122.8 122.7 110.8 121.6 122.1 115.4 Heptane 1.911 0.502 122.8 122.8 110.8 121.6 122.2 115.4 Toluene 2.374 0.501 122.8 122.8 110.8 121.6 122.2 115.4 Chloroform 4.711 0.546 122.8 122.7 110.8 121.6 122.3 115.3 ChloroBezene 5.697 0.542 122.8 122.7 110.8 121.6 122.3 115.3 Tetrahydrofuran 7.426 0.544 122.8 122.7 110.7 121.6 122.4 115.3 DichloroMethane 8.930 0.544 122.7 122.8 110.7 121.6 122.4 115.3 DiChloroEthane 10.125 0.544 122.7 122.8 110.7 121.6 122.4 115.3 2-Methyl-2-Propanol 12.470 0.543 122.7 122.7 110.7 121.6 122.4 115.3 1-Butanol 17.332 0.543 122.7 122.7 110.7 121.6 122.4 115.3 Acetone 20.493 0.543 122.7 122.7 110.7 121.6 122.4 115.3 Ethanol 24.852 0.543 122.7 122.7 110.7 121.6 122.4 115.3 Methanol 32.613 0.543 122.7 122.7 110.7 121.6 122.4 115.3 Acetonitrile 35.688 0.543 122.7 122.7 110.7 121.7 122.4 115.3 DMSO 46.826 0.543 122.7 122.7 110.7 121.7 122.4 115.2 FormicAcid 51.100 0.543 122.7 122.7 110.7 121.7 122.4 115.2 Water 78.355 0.543 122.7 122.7 110.7 121.7 122.4 115.2 Formamide 108.940 0.543 122.7 122.7 110.7 121.7 122.4 115.2 n-MethylFormamide-mixture 181.560 0.542 122.7 122.7 110.7 121.7 122.4 115.2Figure 2
Overlap between the molecular structure determined by X-ray (in red) and in Chloroform (in blue). The ring (a) was used as anchorage.The solvent media presence causes a significant deviation of the RMSD parameter between the X-ray geometry data and the theoretical results, as can be seen in Figure2 for the chloroform (RMSD=0.544). The C4-C7-N1-C8 and C11-N2-C14-C15 torsion angles change due to the solvent medium effects (chloroform) from -166.3° (X-ray data) to +176.24 and +176.4 (X-ray data) to -177.72°, respectively. This effect occurs due to the negative charge transfer to the hydroxyl bond to the terminal phenolic ring (C1–C6). Table S21 (Supplementary Materials) shows the charges of the atoms of the molecule under the effect of several solvent media using the PCM method calculated via DFT (CAM-B3LYP/6-311+G(d)) using the electrostatic model ChelpG. The ChelpG fit of the pyridine ring and hydroxyl charges are 0.04e and -0.238e both in chloroform medium (see the ChelpG charges for all atoms in Table S21 of the Supplementary Materials).DFT results for the EPAF optimized geometry in gas phase and in several solvent media are shown in Table3. The C7—N1—C8 angle presents a reduction of 0.59% for both polar and nonpolar solvent. In the solvent media, the C6—C1—O1 angle increased around 2.0% and the N1—C7—C4 angle remains invariable in gas phase, toluene, heptane, and argon. This angle has increased 0.08% and 0.16% and in chloroform and in the polar solvent medium, respectively. Also the C14—N2—C11 angle increased 0.08% and 0.16% for solvent media with ε < 37 and ε > 37, respectively. The X-ray and theoretical results in the solvent media for the angles, N2—C14—C15 and N3—C15—C14, are practically the same.


## 3.2. HOMO-LUMO

In this section, the HOMO and LUMO results for the EPAF molecule in several solvent media are presented. The HOMO and LUMO orbital and the gap energy values in the solvent medium n-methyl formamide mixture and in gas phase are shown in Figure3. All HOMO-LUMO figures are in the Supplementary Materials (see Figures S1 – S18).Figure 3
HOMO and LUMO in n- methyl formamide mixture.As can be seen in Figure4, the gap energy decreases with the increasing of the dielectric constant value of the solvent medium. The gap energy EG goes from 3.466 eV (ε=1.0) to 3.349 eV (ε=182.4) indicating that the maximum absorption wavelengths (λmax=hc/EG) fall into the ultra violet (UV) region. The gap energy is an important parameter that determines the chemical properties of a molecule as the kinetic stability, chemical reactivity, optical polarizability, etc. Larger values of the gap energy are directly related to greater stability of the molecule.Figure 4
Gap energy as function of the dielectric constant values.All values of HOMO-LUMO and gap energy of the EPAF molecule calculated in various solvent media are in the Supplementary Materials (see TablesS22).


## 3.3. Nonlinear Optical Properties in Solvent Medium

Table4 shows the DFT/B3LYP/6-311+G(d) static results for the total dipole moment (μtot), the isotropic linear polarizability α(0;0), the first hyperpolarizability parallel to the dipole moment (β||z(0;0,0)), and the second hyperpolarizability γ(0;0,0,0) both in gas-phase and in several solvent media. As can be noted, the values of μ, α(0;0) and γ(0;0,0,0) increase with the increasing of the ε-values (Figures 5(a), 5(b), and 5(d)). In all solvent media, the β||z(0;0,0)-values are negative, and for ε≥5 they oscillate around 0.6×10-30esu; however, the β||z(0;0,0) increases with the increasing of the ε-values, as can be seen in Figure 5(c).Table 4
DFT static results for the EPAF electrical parameters.
ε γ ( 0 ; 0,0 , 0 ) β | | z ( 0 ; 0,0 ) α ( 0 ; 0 ) μ (D) 10−36esu 10−30esu 10−24esu 1.000 170 -0.230 42 2.86 1.430 220 -0.200 44 3.01 1.911 254 -0.157 46.4 3.12 2.374 280 -0.127 47.7 3.19 4.711 351 -0.611 51 3.44 5.697 367 -0.581 51.7 3.48 7.426 387 -0.602 52.6 3.53 8.930 399 -0.602 53.1 3.57 10.125 406 -0.603 53.4 3.59 12.470 416 -0.595 53.9 3.61 17.332 435 -0.605 54.7 3.67 20.493 440 -0.607 54.9 3.68 24.852 446 -0.609 55.2 3.70 32.613 448 -0.609 55.3 3.70 35.688 452 -0.609 55.5 3.72 46.826 454 -0.608 55.6 3.72 51.100 458 -0.612 55.6 3.73 78.355 461 -0.604 55.9 3.74 108.940 464 -0.604 56.0 3.75Figure 5
DFT static results for the EPAF electric parameter as function of theε-value.
(a) (b) (c) (d)From Table4, it can be verified that the first and second hyperpolarizabilities present the greater variation due to the solvent medium presence, from heptane (ε=1.911) to n-methyl formamide mixture (ε=181.56); the values of β||z(0;0,0) and γ(0;0,0,0) present a percentage increasing of 285% and 83% (see Figures 5(c) and 5(d) ). However, the character, protic or aprotic, seems to not influence the results, but the nonpolar solvent media present a smaller value of the total dipole moment, and this can be understood taking into account that the polar solvent media the OH present greater increasing of the negative charge, and similar effect can be verified for other electric parameters (Figure 5).The DFT static results for theβHRS-values obtained from (6) in a solvent medium are several orders of magnitude higher than the values of β||z(0;0,0). While as shown in Table 3 for the solvent media the absolute values of β||z(0;0,0) go from 0.2 to 0.6 (in units of 10-30esu), βHRS(0,0,0) go from 24.4 to 53.0 (in units of 10-30esu). Table S23 of the Supplementary Material shows the static and dynamic βHRS-values for the gas phase and all the solvent media. The behaviors of the static and dynamic HRS first hyperpolarizability in the solvent media are similar; therefore, in Figure 6, the results for βHRS as function of the field frequency (ω) are shown for only four solvents, namely, heptane (ε=1.911, nonpolar), chloroform (ε=4.71; larger RMSD parameter), DMSO (ε=46.8; aprotic), and water ( ε=78.4; ETN=1). From Figure 6 and from Table S23, two resonant regions can be identified in the region 0.06a.u.<ω<0.09a.u. So we will work away from the transition region, because in this region it would be complicated for the experiment to work due to deleterious effects. We will study βHRS on the frequency of ω=0.0428 a.u.(1064 nm).Figure 6
Dynamic HR Scattering first hyperpolarizability as function of the field frequency (ω) in a.u..
(a) (b) (c) (d)We can also highlight that the protic or aprotic character and also the polarity of the solvent did not alter the HRS first hyperpolarizability behavior as a function of the frequency. TheβHRS-values increase with the increasing of the dielectric constant value (or with the decreasing of the gap energy) as can be seen in Figure 7. This effect can be explained by the increase of the hydroxyl group charge, which acts as an electrons receptor of the benzene ring, decreasing the electron density and the resonant stability, making it difficult to the electrophilic attack. In this case, the OH group acts as a disabling or metamanager. The values of βHRS in chloroform are 3 to 5 times higher than the 2,4,6-tris(benzylamino)-1,3,5-triazine [24] derivatives also measured in chloroform.Figure 7
Static and dynamic (1064nm) HRS first hyperpolarizability as function of theε-value.


## 3.4. Hirshfeld Surface Analyses

The molecular structure of (E)-4-[({4-[(piridin-2-ilmetilideno) amino]fenil}amino)-metil]fenol (EPAF), C19H17N3O, shows H bonds of the type π- π, C-H... π, C-H..C, O-H...N e N-H...O which contribute to the structure stability. Figure 8 shows the normalized Hirshfeld surface. Figure 8(a) shows the interactions O1-H1...N3 (D...A = 2.796 Å, D-H...A = 178.75°) and C14-H14...C1 (D...A = 3.640 Å, D-H...A = 146.17°). Figure 8(b) identifies the interactions N1-H1A...O1 (D...A = 2.982 Å, D-H...A = 170.07°), O1-H1...N3 (D...A = 2.796 Å, D-H...A = 178.75°), and C14-H14...C1 (D...A = 3.640 Å, D-H...A = 146.17°). Figure 8(c) shows the interactions C17-H17...C8 (D...A = 3,655 Å, D-H...A = 144.90°).Figure 8
View of EPAFdnorm surface showing the interactions O1-H1...N3 e C14-H14...C1 (a), N1-H1A...O, O1-H1...N3, and C14-H14...C1 (b), and C17-H17...C8 (c).Figure9(a) shows the interaction C17-H17...C8 (D...A = 3.655 Å, D-H...A = 144.90°), the same shown in Figure 9(c). The interaction N1-H1A...O1 (D...A = 2.982 Å, D-H...A = 170.07°) is shown in Figure 9(b) and the interactions C7-H7A...C9 (D...A = 3.723 Å, D-H...A = 151.80°) and C7-H7A...C10 (D...A = 3.681Å, D-H...A = 156.93°) are shown in Figure 9(c).Figure 9
View of EPAFdnorm surface showing the interactions (a) O1-H1...N3 and C14-H14...C1, (b) N1-H1A...O, O1-H1...N3, and C14-H14...C1, and (c) C17-H17...C8.Figure10 shows the shape index surface for the EPAF, where it is possible to see two triangles, one blue and the other red, one facing the other, representing the π- π interaction (π...π = 3,779 Å).Figure 10
View of the shape index surface for EPAF showing the highlightedπ...π interaction.The fingerprints plots were used to analyze the amount of intermolecular interactions according to the nature and percentage of the interaction. Figures11(a)–11(g) show the fingerprints for the EPAF; as we can see, the interaction with greater percentage in the whole crystal is the H..H with 45.8% of the interactions followed by the C…H with 34.2% of the interactions.Figure 11
Fingerprint showing the 2D interactions of EPAF crystal.
(a)
C...C 2.1% (b)
C...H 34.2% (c)
H...H 45.8% (d)
N...C 0.9% (e)
N...H 10.6% (f)
O...C 0.1% (g)
O...H 6.4% (h)


## 3.5. Nonlinear Optical Properties of the EPAF Crystal

The supermolecular (SM) method was used to simulate the EPAF crystal; details of the SM approach are given in the work of C. Valverde et al. [25]. To employ the SM approach, we have used the (x-ray) experimental geometry of the asymmetric unit of the EPAF. The packaging effects of the EPAF were modeled by constructing a bulk with the unit cells in a 9×9×9 configuration. Each unit cell contains four asymmetric units, totaling 2916 molecules in the monoclinic system, each atom surrounding the EPAF molecule (blue) being treated as a point charge (see scheme in Figure 12).Figure 12
Scheme of the bulk representing the embedded molecule.For the calculation, the average linear polarizabilityα can be related to the linear refractive index (n) of the crystal via the Clausius-Mossotti relation, given by [26]: (10)n2-1n2+2=4πN3α,whereN stands for the number of molecules per unit cell volume. The experimental quantity, the third-order electric susceptibility χ(3), is related to the second hyperpolarizability by (11)χ3=f4NγϵoV,where N is the number of molecules per unit cell volume (V) andf is the Lorentz local field correction factor given by(12)f=n2+23.The EPAF crystal presents a rapid convergence in the dipole moment through the SM approach [25] (see Figure 13).Figure 13
Evolution of values of the dipole moment of the EPAF crystal with the respective iteration numbers.The quantum molecular calculations were performed with the Gaussian 09 program package [21]. For the calculation of χ3-ω;ω,ω,-ω, we used a small frequency [27, 28], estimate of the second frequency-dependent hyperpolarizability (γ-ω;ω,ω,-ω≅2γ-ω;ω,0,0-γ0;0,0,0) associated to a nonlinear optical process [29] of the intensity dependent refractive index (IDRI) from dc-K results. Table 5 shows the values of the EPAF when we consider the environment via the SM.Table 5
The values of the linear polarizability (in10-24 esu) were calculated CAM-B3LYP/6-311++G(d,p) and for values of the second hyperpolarizability in (in 10-36 esu) we used CAM-B3LYP/6-311++G(d,p) for EPAF embedded molecule (λ=532nm).
α ( - ω , ω ) γ ( 0 ; 0,0 , 0 ) γ ( - ω ; ω , 0,0 ) γ - ω ; ω , ω , - ω 46.50 174.00 628.195 1082.39Table6 presents the DFT prediction and the experimental results of the macroscopic quantities studied. The values of EPAF crystal are therefore 24.51, 2953, 3413, and 2866 times higher, respectively, than the values found experimentally by Prabhu et al. [17] (see Table 6).Table 6
DFT/CAM-B3LYP/ 6-311++G(d,p) results for the linear refractive index and third-order nonlinear susceptibility (10-22m2/V2) for the case dynamic (λ=532nm) of the EPAF crystal.
Sample n ( ω ) χ 3 - ω ; ω , ω , - ω EPAF (this work) 2.02 6793.08 (2E)-3-(3-methylphenyl)-1-(4-nitrophenyl)prop-2-en-1-one (3MPNP) [17] 1.418 277.1 (2E)-1-(4-bromophenyl)-3-[4-methylsulfanyl) phenyl]prop-2-en-1-one (4Br4MSP) [17] 1.363 2.30 (2E)-1-(3-bromophenyl)-3-[4(methylsulfanyl) phenyl]prop-2-en-1-one (3Br4MSP) [17] 1.365 1.99 (2E)-3[4(methylsulfanyl)phenyl]-1-(4-nitrophenyl)prop-2-en-1-one (4N4MSP) [17] 1.360 2.37The high value of the third-order electric susceptibilityχ(3) indicates the crystal as a promising candidate for NLO applications in photonic and optoelectronic devices.


## 4. Conclusions

In this work, using the DFT/B3LYP/6-311+G(d) calculation level, the effects of several solvent media on the geometric and the electric parameters of a Schiff base derivative (E)-4-[({4-[(piridin-2-ilmetilideno)amino]fenil}amino)-metil]fenol (EPAF) [15] were studied. A Hirshfeld surface analysis is presented to complement the EPAF crystal study of [15]. The fingerprints plots show that the interaction with greater percentage in the whole crystal is the H..H with 45.8% of the interactions followed by the C…H with 34.2% of the interactions. The geometry optimization was performed in the gas phase and in various solvent media. The overlap between the X-ray data for the EPAF and the DFT results in several solvent media showed a significant deviation of the RMSD parameter. Also in chloroform, the torsion angles C4-C7-N1-C8 and C11-N2-C14-C15 change from -166.3 (X-ray data) to +176.24 and of +176.4 (X-ray data) to -177.72°, respectively. This effect occurs for all solvents studied here and was due to the negative charge transfer to the hydroxyl bond to the terminal phenolic ring (C1–C6).The NLO properties of the EPAF molecule in several solvent media were studied, and the static electric parameters values of the total dipole moment, linear polarizability, first hyperpolarizability, and the second hyperpolarizability increase with the increasing of the dielectric constant value of the solvent medium. The dispersion relation of the HRS first hyperpolarizability showed two resonant regions forω>0.06 a.u.. The gap energies were calculated from the HOMO-LUMO energy difference in several solvent media calculated and the values go from 3.466 eV (ε=1.0) to 3.349 eV (ε=182.4).The p-nitroaniline (pNA) is used as a critical parameter for comparative studies because it has good NLO properties. The first-order dynamic hyperpolarizability of the molecule p-nitroaniline dissolved in chloroform isβHRSPNA=17.5×10-30cm3/esu at 1064 nm [30, 31]; the EPAF in chloroform is βHRSEPAF=80.75×10-30cm3/esu, about 5 times greater than PNA.The value of the third-order electric susceptibilityχ(3) of the EPAF crystal is (6793.08×10-22m2/V2) for the dynamic case (λ=532nm), thus 24.51 times higher than the values found experimentally by Prabhu et al. [17].Thus, based on the magnitude of the dynamic HRS first hyperpolarizability in chloroform and the third-order electric susceptibility of the EPAF crystal, we can conclude that the EPAF offers potential applications to the development of materials with NLO properties.



---

*Source: 8148392-2019-04-01.xml*