Computational and Experimental Study on Molecular Structure of Benzo [ g ] pyrimido [ 4 , 5b ] quinoline Derivatives : Preference of Linear over the Angular Isomer

A series of 5-aryl-2-methylthio-5,12-dihydrobenzo[g]pyrimido[4,5-b]quinoline-4,6,11(3H) -trione was synthesized through an environmental friendly multicomponent methodology and characterized with FT-IR (Fourier Transform infrared spectroscopy), 1H NMR (Nuclear Magnetic Resonance ), 13C NMR and GC-MS (gas chromatography-mass spectrometry). The 5-(4-methoxyphenyl)2-methylthio-5,12-dihydrobenzo[g]pyrimido[4,5-b]quinoline-4,6,11(3H)-trione 4c compound was characterized by X-ray single crystal diffraction. The geometry of 4c has been fully optimized using DFT (Density functional theory), B3LYP functional and 6-31G(d,p) basis set, thus establishing the ground state energy and thermodynamic features for the mentioned compound, which are in accordance with the experimental data and the crystal structure. The experimental results reveal a strong preference for the regioselective formation of 4c linear four fused rings over the angular four fused and suggest a possible kinetic control in product formation.


Introduction
The use of environmental friendly procedures is nowadays a standard tool in synthetic applications because of their benefits both in chemistry and environmental point of views.In this sense, reactions assisted by microwave (MW) radiation are well known because of their advantages over conventional methods, becoming an important and almost regular methodology for the synthesis of organic compounds [1][2][3][4][5][6][7][8].On the other hand, multicomponent reaction (MCR), most running with atom economy, offer convenient procedures for the introduction of structural diversity of heterocyclic compounds which are prepared by a straightforward manner, in a single synthetic step [9][10][11][12][13][14][15][16][17].Combining the advantages of MCR with those of microwave assisted organic synthesis under solvent-free conditions provides fast and efficient methods for the synthesis of heterocyclic systems, which are among the most common scaffolds in compounds with diverse applications.
In this sense, we have recently reported about simple and environmental friendly MCR one-pot procedures for the synthesis of heterofused-quinolines such as: pyrazolo [3,4-b]quinolindiones [23], prepared by MW assisted synthesis under solvent-free/catalyst-free conditions, or pyrimido [4,5-b]quinolindiones [36] by heating in ethanol.In these reactions, also reported by other researches, the dihydroderivatives are obtained as final products (Scheme 1) [37,38].
Our research group has made great efforts to develop synthetic strategies to obtain highly functionalized heterocycles [13,14,23,[27][28][29]35,36].In this sense, we have recently reported about simple and environmental friendly MCR one-pot procedures for the synthesis of heterofusedquinolines such as: pyrazolo [3,4-b]quinolindiones [23], prepared by MW assisted synthesis under solvent-free/catalyst-free conditions, or pyrimido [4,5-b]quinolindiones [36] by heating in ethanol.In these reactions, also reported by other researches, the dihydroderivatives are obtained as final products (Scheme 1) [37,38].Continuing with our studies for the synthesis of aza-heterocycles, we have decided to investigate the MCR of 6-aminopyrimidin-4-ones 1, naphthalene-1,2,4(3H)-trione 2 and aromatic aldehydes 3 under MW irradiation and solvent-free/catalyst free conditions in order to obtain benzopyrimidoquinolines, and so providing results for a comparative analysis with previous reported experiments [39] and the better knowledge of the reaction type involved in the formation of carbon-carbon bonds [40,41].

Materials and Methods
All reagents used in this work were purchased commercially without further purification.Identifications of compounds and measurements of properties were carried out by general procedures employing the following equipment: Microwave irradiation was carried out with microwave oven CEM Discover (CEM Corporation, Matthews, NC, USA) with controlled power and temperature.Melting points were determined in a Büchi Melting Point Apparatus (BUCHI Latinoamérica, Valinhos, SP, Brasil) and are reported without any corrections.The 1 H and 13 C NMR (Nuclear Magnetic Resonance) spectra were measured at RT (Room Temperature) on a Bruker Avance 400 spectrometer (Bruker, Rheinstetten, Germany) operating at 400 and 100 MHz, respectively, and using Dimethyl Sulfoxide deuterated (DMSO-d6) as solvent and tetramethylsilane (TMS) as internal standard.The mass-spectra were scanned on a Hewlett Packard HP Engine-5989 spectrometer (Hewlett Packard, Palo Alto, CA, USA), equipped with a direct inlet probe which was operating at 70 eV.Elemental analyses were obtained using a LECO CHNS-900 elemental analyzer (LECO Corporation, Saint Joseph, MI, USA).
Continuing with our studies for the synthesis of aza-heterocycles, we have decided to investigate the MCR of 6-aminopyrimidin-4-ones 1, naphthalene-1,2,4(3H)-trione 2 and aromatic aldehydes 3 under MW irradiation and solvent-free/catalyst free conditions in order to obtain benzopyrimidoquinolines, and so providing results for a comparative analysis with previous reported experiments [39] and the better knowledge of the reaction type involved in the formation of carbon-carbon bonds [40,41].

Materials and Methods
All reagents used in this work were purchased commercially without further purification.Identifications of compounds and measurements of properties were carried out by general procedures employing the following equipment: Microwave irradiation was carried out with microwave oven CEM Discover (CEM Corporation, Matthews, NC, USA) with controlled power and temperature.Melting points were determined in a Büchi Melting Point Apparatus (BUCHI Latinoamérica, Valinhos, SP, Brasil) and are reported without any corrections.The 1 H and 13 C NMR (Nuclear Magnetic Resonance) spectra were measured at RT (Room Temperature) on a Bruker Avance 400 spectrometer (Bruker, Rheinstetten, Germany) operating at 400 and 100 MHz, respectively, and using.
Dimethyl Sulfoxide deuterated (DMSO-d 6 ) as solvent and tetramethylsilane (TMS) as internal standard.The mass-spectra were scanned on a Hewlett Packard HP Engine-5989 spectrometer (Hewlett Packard, Palo Alto, CA, USA), equipped with a direct inlet probe which was operating at 70 eV.Elemental analyses were obtained using a LECO CHNS-900 elemental analyzer (LECO Corporation, Saint Joseph, MI, USA).

Chemistry
For the preparation of benzo[g]pyrimido [4,5-b]quinolines, three-component reactions assisted by MW irradiation were carried out using as reagents 6-amino-2-(methylthio)pyri,idin-4(3H)-one 1, naphthalene-1,2,4(3H)-trione 2 and benzaldehyde 3a, in equimolar amounts.Taking as a model reaction the synthesis of 4a, we have tested several reaction conditions, combining diverse solvents (ethanol, acetic acid and ethylenglycol), temperatures and power of the microwave source in order to find the optimal conditions (Table 1).It was found that the product 4a is formed with higher yield under solvent/catalyst-free conditions, (Table 1, entry 2).The synthesis of 4a was carried out at reflux (using ethanol, acetic acid, ethylenglycol or a mixture of both) obtaining low yields and longer reaction times in contrast to MW method.When ethanol was used as the solvent by conventional heating, the desired product 4a was obtained after 2 h in low yield (30%, Table 1, entry 7).In the case of using AcOH or AcOH/Ethylenglycol mixtures moderate similar yields were obtained (45-50%, entry 8-10).
The extension of reaction times, keeping the same reaction conditions, did not lead to the formation of aromatized derivative.On the contrary, a decrease was observed in the reaction yield and double recrystallization of the crude reaction is required, first using a N,N-Dimethylformamide (DMF)-Ethanol mixture (ratio 1-9 v/v) and the isolated solid is then recrystallized from ethanol.
Considering the optimal conditions, that is, solvent/catalysis-free conditions, and in order to show the generality and scope of this protocol, we tried various aromatic aldehydes, 6-amino-2-(methylthio)pyrimidin-4(3H)-one and 6-amino-3-methyl-2-(methylthio)pyrimidin-4(3H)-one, to prepare a series of compounds 4a-q (Scheme 2).The extension of reaction times, keeping the same reaction conditions, did not lead to the formation of aromatized derivative.On the contrary, a decrease was observed in the reaction yield and double recrystallization of the crude reaction is required, first using a N,N-Dimethylformamide (DMF)-Ethanol mixture (ratio 1-9 v/v) and the isolated solid is then recrystallized from ethanol.
Considering the optimal conditions, that is, solvent/catalysis-free conditions, and in order to show the generality and scope of this protocol, we tried various aromatic aldehydes, 6-amino-2-(methylthio)pyrimidin-4(3H)-one and 6-amino-3-methyl-2-(methylthio)pyrimidin-4(3H)-one, to prepare a series of compounds 4a-q (Scheme 2).In all cases, the starting materials were completely consumed, in times ranged between 5 to 9 min, to afford the 5-aryl-2-methyltio-5,12-dihydrobenzo[g]pyrimido [4,5-b]quinoline-4,6,11 (3H)-trione derivatives 4a-q with high yields (70-81%) after easy work-up (Table 2).This reaction permits the use of aromatic aldehydes with electron withdrawing group (EWG) or electron releasing group (ERG), but also to 2-or/and 3-substitutedpyrimidin-4-ones with good yields.As shown in Table 2, the difference between the yields of the compounds obtained are not significant, In all cases, the starting materials were completely consumed, in times ranged between 5 to 9 min, to afford the 5-aryl-2-methyltio-5,12-dihydrobenzo[g]pyrimido [4,5-b]quinoline-4,6,11(3H)-trione derivatives 4a-q with high yields (70-81%) after easy work-up (Table 2).This reaction permits the use of aromatic aldehydes with electron withdrawing group (EWG) or electron releasing group (ERG), but also to 2-or/and 3-substitutedpyrimidin-4-ones with good yields.As shown in Table 2, the difference between the yields of the compounds obtained are not significant, results suggest that the electronic nature of the substituent on aromatic aldehydes and the pyrimidin-4-ones have no significant effect on rate, yields and the course to lead to aromatization of pyridinic ring.
This protocol is a like a Mannich-type reaction, which compared with conventional methods, does not require long reaction times, and the use of flammable solvents, or expensive and toxic catalysts (ionic liquids, superacid, Brønsted acids), to obtain regioselectivelly the reaction products [40,41,43,44].
The structures of 5-aryl-2-methylthio-5,12-dihydrobenzo[g]pyrimido[4,5-b]quinoline-4,6,11(3H) -trione 4a-q derivatives were elucidated from the standart spectroscopic and analytical methods ( 1 H, 13 C NMR, IR, mass spectra and Elemental Analysis (EA).The naphtoquinone system formed by the two non-equivalent nucleophilic center of 6-aminopyrimidinone 1 (-NH 2 group and C-5), is evidenced from 1 H NMR spectra.All derivatives 4a-q showed a singlet between 5.17-5.53ppm assigned to methine proton (stereogenic centre) and other single between 9.21-9.89ppm corresponding to NH of the pyridine ring.All protons exhibit chemical shift signals in their respective expected region.In this reaction, two possible regioisomer could be formed depending on the carbonyl group of naphthalene-1,2,4(3H)-trione 2 involved in the cyclocondensation with the amino group (Figure 1).results suggest that the electronic nature of the substituent on aromatic aldehydes and the pyrimidin-4-ones have no significant effect on rate, yields and the course to lead to aromatization of pyridinic ring.This protocol is a like a Mannich-type reaction, which compared with conventional methods, does not require long reaction times, and the use of flammable solvents, or expensive and toxic catalysts (ionic liquids, superacid, Brønsted acids), to obtain regioselectivelly the reaction products [40,41,43,44].
The structures of 5-aryl-2-methylthio-5,12-dihydrobenzo[g]pyrimido[4,5-b]quinoline-4,6,11 (3H)-trione 4a-q derivatives were elucidated from the standart spectroscopic and analytical methods ( 1 H, 13 C NMR, IR, mass spectra and Elemental Analysis (EA).The naphtoquinone system formed by the two non-equivalent nucleophilic center of 6-aminopyrimidinone 1 (-NH2 group and C-5), is evidenced from 1 H NMR spectra.All derivatives 4a-q showed a singlet between 5.17-5.53ppm assigned to methine proton (stereogenic centre) and other single between 9.21-9.89ppm corresponding to NH of the pyridine ring.All protons exhibit chemical shift signals in their respective expected region.In this reaction, two possible regioisomer could be formed depending on the carbonyl group of naphthalene-1,2,4(3H)-trione 2 involved in the cyclocondensation with the amino group (Figure 1).NOESY (Nuclear Overhauser Effect Spectroscopy) experiments show no correlation between the NH from dihydropyridine (DHP) ring and aromatic protons of the naphtoquinone system as expected for the non-linear derivatives 5.It is clear that the lack of correlation is not sufficient to rule out a structure, but this is confirmed with the rest of spectroscopic analysis.In similar work, we have changed the amino-heterocyclic component (5-amino-1-NH-pyrazole by 6-amino-pyrimidinone), both bearing two non-equivalent nucleophilic center (-NH2 group and C-4 in pyrazole or C-5 in pyrimidinone).In order to obtain pyrazolo [3,4-b]quinolindiones derivatives, the reaction with (1-NH)-5-aminopyrazole, which has an additional nucleophilic center lead regioselectivelly to the formation of non-linear tetracyclic dihydrocompounds.The difference between the works is the change of process regioselectivity, ortho-quinone derivatives (angular) are obtained with aminopyrazoles [23] and para-quinone derivatives (lineal) are obtained with aminopyrimidinones [37].Now, as compared to similar work previously reported [39], the absolute configuration of the products 4 was determined by both the NMR and the X-ray analysis, showing that the linear isomer is isolated [45] including a computational study of the molecular structure of 4 compounds.

Theoretical Calculations
To unambiguously define the structures of 4, we attempted to obtain a crystal of suitable for Xray analysis (Figure 2).Crystal of 5-(4-methoxyphenyl)-2-methylthio-5,12-dihydrobenzo[g] NOESY (Nuclear Overhauser Effect Spectroscopy) experiments show no correlation between the NH from dihydropyridine (DHP) ring and aromatic protons of the naphtoquinone system as expected for the non-linear derivatives 5.It is clear that the lack of correlation is not sufficient to rule out a structure, but this is confirmed with the rest of spectroscopic analysis.In similar work, we have changed the amino-heterocyclic component (5-amino-1-NH-pyrazole by 6-amino-pyrimidinone), both bearing two non-equivalent nucleophilic center (-NH 2 group and C-4 in pyrazole or C-5 in pyrimidinone).In order to obtain pyrazolo [3,4-b]quinolindiones derivatives, the reaction with (1-NH)-5-aminopyrazole, which has an additional nucleophilic center lead regioselectivelly to the formation of non-linear tetracyclic dihydrocompounds.The difference between the works is the change of process regioselectivity, ortho-quinone derivatives (angular) are obtained with aminopyrazoles [23] and para-quinone derivatives (lineal) are obtained with aminopyrimidinones [37].Now, as compared to similar work previously reported [39], the absolute configuration of the products 4 was determined by both the NMR and the X-ray analysis, showing that the linear isomer is isolated [45] including a computational study of the molecular structure of 4 compounds.

Theoretical Calculations
To unambiguously define the structures of 4, we attempted to obtain a crystal of suitable for X-ray analysis (Figure 2).Crystal of 5-(4-methoxyphenyl)-2-methylthio-5,12-dihydrobenzo[g]pyrimido [4,5-b]quinoline-4,6,11(3H)-trione 4c was grown in ethanol (95%) using the slow evaporation method at room temperature (RT) [45].In accordance with the X-ray structural analysis of 5-(4-methoxyphenyl)-2-methylthio-5, 12-dihydrobenzo[g]pyrimido[4,5-b]quinoline-4,6,11(3H)-trione 4c, the computational studies indicate that the dihydropyridine ring in these compounds has a boat conformation due to the presence of a sp 3 hybridized C5-atom in which the aryl ring occupies the pseudo-axial position.The orientation of the additional substituent on the C5-aryl ring (CH3O) located on the para-positions with respect to the heterocyclic ring may be syn-orientated to the C5-H bond (synperiplanar, C5-H sp), or may lie above the heterocyclic ring and anti-orientated to the C5-H bond (antiperiplanar, C5-H ap) (Figure 3).The optimized structures of the compound considered in the present study are illustrated in Figure 3, in which the CH3O group has the more stable planar conformation.Structural and bonding analysis of these compounds is started by comparison of the selected bond lengths, and bond dihedral angles.Orientation of the CH3S group with respect to the heterocyclic ring pyrimidine is denoted by a dihedral angleτ1 (N1-C2-S19-C20), inter-ring dihedral angles defining orientations of the aryl ring towards the heterocyclic ring are indicated as τ2 (C14-C5-C24-C29), and the inner dihedral angles of the dihydropyridine ring are denoted by τ3 (N12-C13-C14-C5) and τ4 (N12-C18-C15-C5), and dihedral angle In accordance with the X-ray structural analysis of 5-(4-methoxyphenyl)-2-methylthio-5,12dihydrobenzo[g]pyrimido[4,5-b]quinoline-4,6,11(3H)-trione 4c, the computational studies indicate that the dihydropyridine ring in these compounds has a boat conformation due to the presence of a sp 3 hybridized C5-atom in which the aryl ring occupies the pseudo-axial position.The orientation of the additional substituent on the C5-aryl ring (CH 3 O) located on the para-positions with respect to the heterocyclic ring may be syn-orientated to the C5-H bond (synperiplanar, C5-H sp), or may lie above the heterocyclic ring and anti-orientated to the C5-H bond (antiperiplanar, C5-H ap) (Figure 3).The optimized structures of the compound considered in the present study are illustrated in Figure 3, in which the CH 3 O group has the more stable planar conformation.Structural and bonding analysis of these compounds is started by comparison of the selected bond lengths, and bond dihedral angles.Orientation of the CH 3 S group with respect to the heterocyclic ring pyrimidine is denoted by a dihedral angleτ1 (N 1 -C 2 -S 19 -C 20 ), inter-ring dihedral angles defining orientations of the aryl ring towards the heterocyclic ring are indicated as τ2 (C 14 -C 5 -C 24 -C 29 ), and the inner dihedral angles of the dihydropyridine ring are denoted by τ3 (N 12 -C 13 -C 14 -C 5 ) and τ4 (N 12 -C 18 -C 15 -C 5 ), and dihedral angle τ5 (C 28 -C 27 -O 44 -C 45 ) the 4-metoxi group the 5-aryl, which are extracted from the optimized structures of compound 4c (Figure 3) and listed in Table 3.These orientations may influence the molecular energy content and the changes of some characteristic bond lengths, bond angles and bond dihedral angles.To study the vibrational characteristics and structural parameters of 5-(4-methoxyphenyl)-2methylthio-5,12-dihydrobenzo[g]pyrimido [4,5-b]quinoline-4,6,11(3H)-trione 4c, we performed a theoretical study at DFT level, which allows us to set the geometric and energetic parameters of the title compound.The optimized structure with B3LYP/6-31G(d,p) level is shown in Figure 2 [46][47][48][49][50].The calculated geometric parameter (bond lengths, bond angles and dihedral angles) at same levels of calculation for title compound was compared with the experimental parameters, see Table 3, showing very good correlation.
The assignment of experimental vibrational bands to normal vibration modes is based on the comparison with related molecules and with results of the calculations obtained.We consider the B3LYP/6-31G(d,p) calculations because the scale factors used are well defined for this base set and also because the compared molecule was optimized at the same calculation level (0.961 ± 0.045) [51].Vibrational spectral assignments (Table 4) were performed on the recorded FT-IR spectra, based on the theoretically predicted wavenumbers by DFT methods, using 6-31G(d,p) basis set [50].The observed spectra are in good agreement with the simulated spectra.
The assignment of experimental vibrational bands to normal vibration modes is based on the comparison with related molecules and with results of the calculations obtained.We consider the B3LYP/6-31G(d,p) calculations because the scale factors used are well defined for this base set and also because the compared molecule was optimized at the same calculation level (0.961 ± 0.045) [51].Vibrational spectral assignments (Table 4) were performed on the recorded FT-IR spectra, based on the theoretically predicted wavenumbers by DFT methods, using 6-31G(d,p) basis set [50].The observed spectra are in good agreement with the simulated spectra.The assignment of experimental vibrational bands to normal vibration modes is based on the comparison with related molecules and with results of the calculations obtained.We consider the B3LYP/6-31G(d,p) calculations because the scale factors used are well defined for this base set and also because the compared molecule was optimized at the same calculation level (0.961 ± 0.045) [51].

N-H vibrations
The bands at 3459 cm −1 and 3450 cm −1 are assigned to the symmetrical stretches of the N-H groups, N 12 and N 3 , respectively.These bands appear overlapping in the experimental spectrum at 3272 cm −1 .

C=O groups vibrations
The C=O stretches are observed as medium intensity bands at 1674 cm −1 , 1650 cm −1 and 1620 cm −1 .The DFT calculation gives the stretch wave number at 1784 cm −1 for stretching the C 4 =O 21 group, 1756 cm −1 for C 11 =O 23 and 1745 cm −1 for C 6 =O 22 after scaling.The difference between the calculated and experimental wave numbers can be attributed to the conjugation of C=O bonds to the phenyl ring, which is expected to decrease the wave numbers of the stretches.

C=N vibrations
Stretch bands C=N expected in the range of 1672-1566 cm −1 .For compound 4c, the stretching mode C=N is assigned to 1496 cm −1 and 1450 cm −1 in FT-IR spectrum and theoretically at 1495 (1436) and 1467 (1409) cm −1 .
Two factors may be responsible for the discrepancies between the experimental and computed wavenumbers of the compound.The first is caused by the environment (gas and solid phase) and the second is because the experimental values are inharmonic wavenumbers while the calculated values are harmonic ones.Therefore, the calculated values are very close to experimental measures.
Natural bond orbital (NBO) analysis allows a detailed description of the electronic structure of compound 4c, in terms of occupancy and composition of the group of NBOs Lewis (see Table 5) and not Lewis (see   In Lewis type orbitals, as seen in Table 5, the σ bond (C 4 -O 21 ) is formed from the sp 2.01 hybrid on carbon, mixture of s (33.23%) p (66.68%) d (0.09%); on the other hand, the π bond (C 4 -O 21 ) is formed from the sp 1.00 hybrid on carbon and oxygen, mixture of s (0.01%) p (99.81%) d (0.18%); the electronic densities or maximum occupations calculated for σ (C 4 -O 21 ), π (C 4 -C 21 ), σ (C 6 -C 22 ) and π (C 6 -N 22 ) are 1.99441, 1.98431, 1.99483 and 1.94776, respectively.Therefore, these results allow to infer that these bonds are controlled by the character p of the hybrid orbitals.
The energy values E 2 for the interaction between the filled orbital i (donors) and the vacant orbital j (acceptors) or other molecular subsystem, calculated according to the theory of the second order perturbation, Equation (1), predicts the occurrence of delocalization or hyperconjugation [55,56].The higher the value of E 2 , the more intense the interaction between electron donors and the greater the degree of conjugation of the whole system [57].
q i = Occupation or electronic density of the donor orbital.E i and E j are the diagonal elements and F(i, j) is the element of the Fock NBO diagonal matrix.Therefore, we calculated the hyperconjugative interaction and density transfer of lone pair electron (LP) of the N 3 and N 12 atoms to neighboring π antibonding orbitals.Similarly, the conjugation of carbonyl groups to neighboring atoms via π* → π*, see Table 7.In Table 7, the important contribution of LP N → π* interactions to the overall stability of the system is observed.Stability explained by the intramolecular electron transfer that occurs with these interactions.As evidence of this orbital phenomenon, it is the lowest electronic density on N3 and N12 compared to N1 (Table 6).The energy involved in the hyperconjugative interactions of the bonds between carbon and oxygen atoms of the carbonyl groups of molecule 4c, give the most intense interactions and correspond to the highest degree of conjugation occurring in the α,β-unsaturated intramolecular system of aromatic ketones.In the analysis of the nonlinear optical properties (NLO) for 4c, the polarization of the molecule, induced by an external radiation field, it often resembles the generation of a dipole moment induced by an external electric field.Under the weak polarization condition, a Taylor serial development in the electric field components can be used to demonstrate the dipole interaction with the electric field of external radiation.
The first static hyperpolarizability (β 0 ) and its related properties (β, α 0 and ∆α) were calculated using the B3LYP/6-31G(d,p) level based on the finite field approach.The total static dipole moment µ, average polarizability α 0 , anisotropy polarizability ∆α and hyperpolarizability of first average β 0 , using the components x, y and z is defined as: (2) The polarizability and hyperpolarizability data are reported in atomic units (a.u), by Gaussian 09, reason why were converted into electrostatic units (esu) (α: 1 au = 0.1482 × 10 −24 esu; β: 1 au = 8.639 × 10 −33 esu).The mean polarization α 0 and total polarizability ∆α for 4c are 45.46 × 10 −24 esu and 30.66 × 10 −24 esu, respectively.The total molecular dipole momentum and first order hyperpolarizability are 1.88 Debye and 3768.07 × 10 −33 esu, respectively, see Table 8.These analysis frequently, are compared with the family of urea, as a reference for characterizing organic NLO materials [58,59].The total dipole moment of the molecule is about 2.22 times lower than urea, and hyperpolarizability first order is 5.89 times higher than urea (µ and β for urea are 4.1775 Debye and 638.906 × 10 −33 esu at the same calculation level).This result indicates the good non-linearity of the molecule.
The intramolecular charge transfer process, is determined by the separation E between the HOMO (highest energy occupied molecular orbital) and LUMO (lowest energy unoccupied molecular orbital) levels.Energy susceptible to adjust to λ in the range of visible light or tunable laser, by controllable variations due to the nature donor-acceptor of the groups.As an indirect method for determining potential NLO properties, the basal and excited state of a molecule is calculated and these results are used in the qualitative estimation of the efficiency of the intramolecular charge transfer process [60,61].Which is related to the intramolecular electronic transfer in 5-arylaryl-2-methylthio-5,12-dihydrobenzo[g]pyrimido [4,5-b]quinoline-4,6,11(3H)-trione derivatives, the carbonyl groups being the acceptors of this charge transfer.

Figure 3 .
Figure 3. Conformation of the dihydropyridine ring in the structure of compound 4c.

Figure 3 .
Figure 3. Conformation of the dihydropyridine ring in the structure of compound 4c.

Table 2 .
Results of microwave (MW)-assisted synthesis of 4a-q compounds, in solvent/catalyst-free conditions.
a Compounds which showed decomposition.* Isolated yield.

Table 2 .
Results of microwave (MW)-assisted synthesis of 4a-q compounds, in solvent/catalyst-free conditions.
a Compounds which showed decomposition.* Isolated yield.

Table 3 .
X-ray experimental analysis data and lengths, angles and dihedral angles bonds of 4c calculated by Density functional theory (DFT) method.

Table 4 .
Fourier Transform infrared spectroscopy (FT-IR) Experimental and computed vibrational

Table 4 .
Fourier Transform infrared spectroscopy (FT-IR) Experimental and computed vibrational

Table 6 )
. 224 NBOs were calculated for this compound, representative NBOs are highlighted by finding three carbonyl-like bonds (C 4 =O 21 , C 6 =O 22 , C 11 =O 23 ), a non-aromatic double bond (C 15 =C 18 ), four heteroatom bonds (C 2 -S 19 , C 27 -O 44 , N 12 -H 31 and N 3 -H 30 ) and solitary pairs of N and O (N 1 , N 3 , N 12 , S 19 , O 21 , O 22 , O 23 and O 44 ) of the expected Lewis structure.For each of the molecular orbitals, the corresponding percentage of s, p, d character is included.

Table 5 .
Occupation of natural orbitals and hybrids for 4c calculated by the DFT/B-3LYP/6-31G(d,p) method for the representative atoms.
a Lone pair on natural Lewis structure.

Table 6 .
Occupation of natural bond orbitals (NBO) no Lewis and hybrids for 4c calculated by the DFT/B-3LYP/6-31G(d,p) method for the representative atoms.

Table 7 .
Analysis perturbation theory second order in Fock matrix in NBO by level calculation B3LYP/6-31G (d,p) for compound 4c.

Table 9 .
Ground state energies (enthalpy, Gibbs free energy and entropy) calculated for 4c and 5c compounds.