A density functional theory study of the molecular structure, reactivity, and spectroscopic properties of 2-(2-mercaptophenyl)-1-azaazulene tautomers and rotamers

Five stable tautomer and rotamers of the 2-(2-Mercaptophenyl)-1-azaazulene (thiol, thione, R1, R2, and R3) molecules were studied using density functional theory (DFT). The geometries of the studied tautomer and rotamers were fully optimized at the B3LYP/6-31G(d,p) level. Thermodynamic calculations were performed at M06-2X/6-311G++(2d,2p) and ωB97XD/6-311G++(2d,2p) in the gas phase and ethanol solution conditions modeled by the solvation model based on density (SMD). The kinetic constant of tautomer and rotamers conversion was calculated in the temperature range 270–320 K using variational transition state theory (VTST) accompanied by one-dimensional wigner tunneling correction. Energy refinement at CCSD(T)/6–311++G(2d,2p) in the gas phase has been calculated. All the studied DFT methods qualitatively give similar tautomer stability orders in the gas phase. The ethanol solvent causes some reordering of the relative stability of 2-(2-Mercaptophenyl)-1-azaazulene conformers. The transition states for the 2-(2-Mercaptophenyl)-1-azaazulene tautomerization and rotamerization processes were also determined. The reactivity, electric dipole moment, and spectroscopic properties of the studied tautomer and rotamers were computed. The hyper-Rayleigh scattering (βHRS), and depolarization ratio (DR) exhibited promising optical properties when nonlinear optical properties were calculated.

Tautomerism provides knowledge of the variety of forms in which organic molecules can exist.Many biologically active compounds, including pyrazolone derivatives, take part in tautomeric transformation 1 .In numerous fields of chemistry, intramolecular proton transfer (IPT) is crucial 2 .Its mechanism takes place in both ground and excited states 3 .Compounds that exhibit intramolecular proton transfer can be used as laser dyes, fluorescent probes, high-energy radiation detectors, memory storage devices, and polymer protectors 4 .Hence, molecules that have intramolecular proton transfer have received great interest from an experimental and theoretical perspective 5,6 .Numerous topics in organic chemistry and biochemistry can be studied more effectively by understanding the nature of tautomeric equilibrium 7 .
Thione-thiol tautomerization is a remarkable subject that has attracted a lot of study.The thione tautomer is mostly detected in the gas phase under standard experimental conditions.Thiol tautomers of thione compounds, which can be stabilized in inert matrices at lower temperatures, can be subjected to UV light.These compounds include (1H)-pyridinethione 8 , 4(3H)-pyrimidinethione and 3(2H)-pyridazinethione 9 , 2(1H)-quinolinethione 10 , methimazole 11 , 2-thiobenzimidazole 12 , thiourea 13 , and thioacetamide 14 .Thiophenol is a notable exception that only exists as a thiol tautomer 15 .Organosulfur compounds containing thiocarbonyl groups have attracted interest because of their diverse chemistry 16 , various biological applications 17,18 , and the ability of the C = S bond to oxidize quickly to form the corresponding thiols and disulfides 19 .

Computational details
The commonly used B3LYP method has provided the correct molecular geometries, activation energies, and energy differences between pairs of tautomers, similar to MP2-level predictions 30 .It has been concluded that the small absolute errors and the cost-effectiveness of the computations make the B3LYP method a good choice for these systems 30 .Our computations were carried out using the gradient-corrected hybrid density functional B3LYP DFT method 31 .For each structure, a full geometry optimization without symmetry constraints was performed using this functional along with the 6-31G(d,p) basis set 32 as implemented by the Gaussian 09 package 33 .All geometries were visualized either using GaussView 5.0.9 34 or ChemCraft 1.6 35 programs.
The solvation effect has been modeled in ethanol using the solvation model based on density (SMD) 41 at the B3LYP/6-31+G(d,p) level for all the tautomers and rotamers.We estimated solvation energies in ethanol as single point energy calculations on geometries optimized in ethanol at B3LYP using M06-2X, B3LYP, and ωB97XD functionals with 6-311++G(2d,2p) basis sets.
On the basis of the statistical thermodynamics principle 42 , the entropy (S) and enthalpy (H) of the tautomers in all environments at standard conditions of T = 298.15K and P = 1 atm have been obtained.Gibbs free energy values, which were calculated at the same levels of theory, were used for the evaluation of the tautomeric equilibrium constants.For the tautomerization and rotamerization transformation processes, rate coefficients were calculated in the temperature range 270-320 K (k uni , in s −1 ) using the variational transition state theory (CVT) 43 that includes generalized transition state rates k GT with wigner tunneling correction ( χ(T)) 44 according to (1-3) equations.
where σ, k b , T, h, R, and s refer to the reaction path degeneracy, Boltzmann constant, temperature in Kelvin, Planck constant, ideal gas constant, and the distance along the minimum energy reaction path (MEP) in isoinertial coordinates, respectively.
Global chemical reactivity descriptors 45 have been calculated using the energies of the highest occupied and lowest unoccupied molecular orbitals (HOMO and LUMO), respectively, to better comprehend the reactivity and stability of the structures under study.Accordingly, for B3LYP/6-31G(d,p) in the gas phase, the ionization potential (IP), electron affinity (EA), absolute hardness (η), softness (S), and electronegativity (χ) were calculated.The following formulas have been used to calculate the global chemical reactivity descriptors: Nucleus-independent chemical shift (NICS) computations have been performed for each compound in the gas phase at the same level of optimization (B3LYP/6-31G(d,p)).Particularly, the computed NICS values are NICS(0)iso, which refers to the absolute isotropic shielding at the center of the ring; NICS(0)zz, which is the total contribution to the orthogonal to the plane component of the NICS tensor, assuming the molecule lies in the xyplane; as well as NICS(1)iso and NICS(1)zz, which are analogue indices computed at 1 A above the center of the ring 46,47 .The geometry-based index harmonic oscillator measure of aromaticity (HOMA) 48,49 , which makes use of bond lengths in accordance with Krygowski's 48 instructions, is another criterion of aromaticity.Electrostatic potential (ESP) surfaces have been derived from the B3LYP/6-31G(d,p) calculation.
UV-vis spectra for the investigated structures were computed at TD-DFT PBE 50,51 (TD-PBE/6-311+G(d,p)) using the SMD method at the B3LYP/6-31G(d,p) optimized gas-phase geometry in acetonitrile.The electron excitations of 2OHPhAZ were precisely determined by the Perdew-Burke-Ernzerhof technique (PBE) 22 .Gauss-Sum program 52 was used to simulate the ultraviolet-visible (UV-vis) spectra.Instead of focusing on the canonical orbitals, natural transition orbitals (NTOs) 53 were estimated for each electron excitation.We create the frontier orbitals and NTOs using Chemcraft.Additionally, we use the optimized gas-phase geometry to calculate the NMR spectra of the various tautomers and rotamers relative to the 13C and 1H isotropic chemical shielding of tetramethylsilane (TMS) at the B3LYP/6-31G(d,p) in chloroform using the gauge-independent atomic orbital (GIAO) method 54,55 .
We concentrated on the hyper-Rayleigh scattering (β HRS ) and depolarization ratio (DR) among second-order NLO characteristics 56 , and the complete equations for calculating the magnitude of the total dipole moment μ tot , the average polarizability α tot and the first hyperpolarizability β tot , using the x, y and z components are as follows: 57,58 (1)
2-(2-Mercaptophenyl)-1-azaazulene could exist as thione-thiol tautomerism based on the movement of hydrogen atom among nitrogen and sulfur, and three rotamers form R1, R2, and R3.The optimized geometries of the studied tautomeric and rotameric forms and their corresponding inter-conversion transition states are shown in Fig. 2, and their geometrical parameters are listed in Table 1 at the B3LYP/6-31G(d,p) level.Also, Table S1 in the supplementary information (SI) reported their coordinates.The rotamer R1 is formed by rotating the OH group of thiol, while R2 and R3 are formed by rotating the phenyl ring of thiol and R1, respectively.To the best of our knowledge, there are no theoretical and experimental reports on the geometry of the titled compound conformers in the literature for comparison.The calculated geometrical parameters are in good agreement with the corresponding values reported for a similar compound 22 .
The thione form in the gas phase shows that it is almost planar with a dihedral angle (N1C2C1'C2') of 0.0 o , but in the case of the thiol form, the terminal phenyl ring is out of plane with a dihedral angle (N1C2C1'C2') of-13.9°.The thione-thiol tautomerization occurs via an IPT reaction.The intramolecular hydrogen bond generates a six-membered ring in both the thiol and thione tautomers.The intramolecular H-bonding is between the S1-H1' donor and N1 acceptor (1.808 Å) of the thiol tautomer, but in the case of the thione tautomer, the interaction is between the N-H donor and the thiophenol S1 acceptor, where the hydrogen-bond lengths are about 2.015 Å.In the IPT, the C2′-S1 bond length decreases from 1.779 in the thiol form to 1.735 Å in the thione form.Also, the calculated S-H1' …N hydrogen-bond angles decrease from 143.3° in the thiol form to 142.2° in the thione form.
The thiol-thione tautomer passes through the transition state (TST) (see Fig. 2).The extracted observations of the proton shift process include: (1) The S1′-H1′ bond length of 0.968 Å in thiol elongated by 0.203 Å in the TST and eventually separated to release proton; (2) The proton then migrated and approached N1 with a distance of 1.329 Å in the TST before the formation of the N1-H1' bond (1.056 Å) in thione; (3) In concerted processes, the thiol S1C2'C1' and N1 C2 C1' angles closed up by 6.2° and 1.7°, respectively, in the TST and finally settling at 126.5• and 122.7•, respectively, in thione.
In the process of rotamerization, the observed slightly increased C2'-S1 bond distance for R2 and R3 (1.794 Ǻ) compared to R1 (1.789 Ǻ) is due to the decrease in double bond character in R2 and R3 because there is no conjugation of the sigma lone pair in R2 and R3.Additionally, there is no hydrogen bond in the R1 rotamer, and the donor-acceptor distance of 2.773 Ǻ is smaller than the corresponding distance in thiol of 3.021 Ǻ.The S1 − H1' bond length in TSR1 and R1 is similar but shorter than that in thiol by 0.02 Å due to the presence of HB in the latter structure.The H1' − S1 − C2' angle is rotated by 3.4° from R1.The donor-acceptor distance in TSR1 is elongated to 3.249 Å.The calculated dihedral angles show that the planarity of rotamers is distorted (Table 1).The non-rigid structure of 2-(2-Mercaptophenyl)-1-azaazulene tautomers and rotamers gives them the ability to fit with receptors in biological systems.
Tautomerism and thermodynamic properties.This section deals with the study of equilibrium between the studied tautomers and rotamers using the DFT/M06-2X and ωB97XD in comparison to the more accurate CCSD (T) calculations.Table 2 presents the theoretically computed corrected relative energies (ΔE), Gibbs free energies (ΔG), and enthalpies (ΔH) of the five tautomers and rotamers at the ωB97XD/6-311++G(2d,2p), M06-2X/6-311++G(2d,2p), and CCSD(T)/6-311++G(2d,2p) levels in the gas and ethanol phases.For tautomerization, it is found that the thiol form is more stable than the thione form in the gas phase.The relative energy difference between the thiol and thione in the gas phase is 3.30 and 5.30 kcal/mol at ωB97XD and CCSD(T), respectively.The thiol stability in the gas phase is expected by having a lower and stronger IHB than thione and by having the higher aromaticity of the phenyl ring that has been broken by the formation of the thione form (as will be discussed in Section "Aromaticity").
In the rotamerization process, the order of stability of the rotameric structures is R2 > R3 > R1.Rotamer R2 is slightly more stable than R3 and R1 by about 0.42 and 0.44 kcal/mol at ωB97XD/6-311 + +G(2d,2p), respectively.R1 is the least stable rotamer at all calculated levels.The stability of R2 and R3 can be attributed to the possible Vol.:(0123456789) www.nature.com/scientificreports/intramolecular H bonding interaction of the C6′-H6′…N1 as shown in Fig. 1.The slightly higher stability of R3 than R2 can be due to the anti-position present between the S1-H1' group and the nitrogen of the azaazulene ring, given the possibility of H3 interaction with the sulfur atom.Mezey et al. 63,64 proposed that any species with a higher energy than 10 kcal mol −1 above the most stable form would not exist in any appreciable concentration.So, we can conclude that all the studied tautomers and rotamers of 2-(2-Mercaptophenyl)-1-azaazulene can be observed and exist experimentally by having a low difference in   2).A chemical equilibrium between tautomers and rotamers is characterized by the equilibrium constant (K) that can be calculated using: where the gas constant (R) is 1.987*10 -3 kcal/mol and the temperature (T) is 298.15K and the ΔG is the difference between the Gibbs free energies of a given tautomer or rotamer relative to thiol.Table 3 collects the K values for the four possible equilibrium reactions in the gas phase and ethanol at ωB97XD/6-311 ++G(2d,2p) and CCSD(T)/6-311++G(2d,2p).
In accordance with the highest stability of thiol among the other tautomers and rotamers in the gas phase, K 1 has the lowest values.This implies the presence of thiol alone in the gas phase.The values of K 2 , K 3 and K 4 (see Table 3) indicate that the equilibrium is effective and the rotamers with the thiol are present in comparable quantities in both the gas phase and ethanol.The calculated rate constant for the studied transformation with VTST and wigner tunneling is given in Table 4.The results show that the tautomerization transformation has a higher rate and a higher tunneling correction effect than rotamerization in the gas phase during the applied temperature range.The rate of rotamerization seems to slightly increase as the temperature is raised.
The heterocyclic tautomerism is very sensitive to the solvent nature.Therefore, more efforts have been made to study the effect of solvents such as ethanol on the stability of the tautomers and rotamers using the same levels of theory.Energy results obtained from these calculations are collected in Table 2.In contrast to the gas phase, the higher stability of the thione form than thiol in ethanol has been observed at all levels of calculations.The stability order is reversed in ethanol, and the thione becomes more stable relative to thiol by 4.50 kcal/mol at the ωB97XD/6-311++G(2d,2p) level (Table 2).Also, the higher value of the equilibrium constant (1.208*10 4 ) emphasized the higher stability of thione in the ethanol solvent (Table 3).
The calculated dipole moment values for all tautomers and rotamers at the ωB97XD/6-311++G(2d,2p) level of theory in the gas and ethanol phases are given in Table 2.It can be noted that thione has the highest polarity.As a result, it has a significant impact on the stability order in ethanol as a polar solvent.Obviously, a larger dipole moment leads to greater stabilization.As can be seen in Table 2, the thione has 16.31 D, while the thiol, R1, R2 and R3 have 7.97, 4.49, 5.76 and 3.76 D, respectively, at ωB97XD/6-311++G(2d,2p) level of theory in the ethanol phase.Comparing the previous findings with those of the ethanol phase, it is observed that the relative energy difference is positive for all the tautomers and rotamers.This result suggests that the presence of ethanol helps stabilize them, with a clear stabilization of five tautomers and rotamers.For the given tautomers and rotamers, the results show that the thiol form has the lowest formation enthalpy values (Table 2) in the gas phase.This means that it is more stable than the thione form.The lower positive reaction enthalpies mean that the equilibrium reaction is an endothermic one.
Chemical reactivity.Theory of frontier molecular orbitals.HOMO and LUMO orbitals play a fundamental role in the qualitative interpretation of chemical reactivity 71 .The energy gap between HOMO and LUMO frontier orbitals is one of the significant characteristics of molecules and plays an important role in electric properties, electronic spectra, and photochemical reactions.The HOMO-LUMO energy gap helps characterize the chemical reactivity and kinetic stability of the molecule 72 .A molecule with a high energy gap (ΔE) is less polarizable and is generally associated with low chemical reactivity and high kinetic stability 73 .Table 5 contains the energies of the HOMO LUMO boundary orbitals and the HOMO-LUMO energy gaps in gas and ethanol phases.The energy gap values decrease in the ethanol (Table 5) compared with those obtained in the gas phase.Therefore, the tautomers and rotamers are more reactive in ethanol.The obtained value of the thione tautomer (ΔE gap = 1.87 and 0.08 eV) is smaller than those obtained with thiol (ΔE gap = 2.98 and 0.12 eV) in gas and ethanol, respectively.So, thione is the most polarizable and has the highest interactions for intramolecular charge transfer.The large energy gap of the thiol confirms the high stability of the thiol tautomer.The rotamers (R1-R3) are hard molecules and have the least ability to polarize because they have a higher energy gap than tautomers.R3 has the largest value of the energy gap, 3.44 and 0.13 eV in the gas phase and ethanol, respectively.
Global indices of reactivity.The study of the global reactivity of molecules is based on the calculation of global indices deduced from electronic properties.The global indices of the reactivity in the gas phase and ethanol of the studied conformers are recorded in Table 5.Chemical hardness (η) and global softness (S) express the resistance of a system to a change in its number of electrons.In a given series of molecules, when η is weak, the molecule is called soft and when it is high, the molecule is called hard.This is quite the opposite of softness, which evolves in the opposite direction of hardness 74 .The value of the chemical hardness of thione is the lowest, where η = 0.93 and 0.04 eV in the gas phase and ethanol, respectively.Also, we note that the thione form has a lower electronegativity value (χ = 3.77 eV) than other tautomers and rotamers in the gas phase, which supports the stability of thiol in the gas phase.

Electrostatic potential (ESP).
The electrostatic potential (ESP) surface is another useful tool to understand the chemical reactivity of a molecule.The ESP study shows that H-donor and H-acceptor properties of molecules are revealed by positive and negative regions, respectively 75 .
By using ESP maps, the charged areas of the molecule can be shown.In ESP maps, the regions with the highest negative electrostatic potential are shown in red, while those with the most positive electrostatic potential are shown in blue.Figure 3 shows the ESP surfaces of the examined structures derived using the B3LYP/6-31G (d,p).The ESP of thiol, thione, and R1 shows that a substantial amount of negative charge is localized on the S atom, with just a tiny amount of negative charge localized on the sulfur atom for R1.
As a result, the S atom is believed to have the highest electron donation ability towards the formation of a H-bond or metal ions.The blue colour around the H atoms indicates their positive charge.For R2 and R3 the negativity charge becomes localized on the N atom.Therefore, R2 and R3 have a high tendency to bind with metal ions on the N side.www.nature.com/scientificreports/Aromaticity.The aromaticity of organic compounds is one of the most important characteristics related to their specific chemical reactive structure 76 .In this study, we used some aromaticity indicators that depend on geometry.HOMA is one of the best ways to describe the change in aromaticity 48,49 .When HOMA is unity 66 , then the compound is completely aromatic, while if HOMA equals zero (0), the compound or ring is completely non-aromatic.More aromatic rings have more π-electrons delocalized, and this is indicated by rings with higher HOMA values.HOMA can be calculated according to the following equation: where α i is a normalization constant (α CC = 257.7 and α CN = 93.52)and n is the total number of bonds in the molecule.The optimized bond length, or R opt , is 1.388 Å for a C-C bond and 1.334 Å for a C-N bond.The results of HOMA calculations at the B3LYP/6-31G(d,p) level in the gas and ethanol phases are shown in Table 6.Evidently, Table 6's HOMA results offer a good match with the relative energies and the structure of the thiol-thione tautomer in both gas and ethanol solvent.The strength of HB of thiol in the gas phase can be attributed to the highest value of HOMA (0.89) of the thiophenyl ring.Through the formation of thione tautomer, the five-and six-membered rings exhibit reasonably higher and lower delocalization, with HOMA values of 0.61 and 0.67, respectively.The five-membered ring's lower HOMA values in the thiol form (0.48) compared to the thione (0.61) tautomer indicate that its electron delocalization to the IHB ring in the thiol form is involved.With a slightly higher HOMA value than the five-membered ring due to tautomerization, the seven-membered ring indicates a greater aromatic character.From the HOMA values of all rings, thiol seems to have strong π-stacking ability, hence having more aromaticity and stability in gas phase.In ethanol, the aromaticity of the five-and sevenmembered rings of thione is increased (0.66 and 0.76, respectively), which raises the summation of HOMA.Thus, the aromaticity and stability of thione are noticed in ethanol solvent.The comparison of HOMA values for the thiophenyl ring and five-membered ring for R1 in the gas phase (0.92 and 0.44), R2 (0.94 and 0.46) and R3 (0.93 and 0.47) seems to suggest a stronger delocalization and the higher stability of the R2 structure.A slight change observed in ethanol preserves the higher stability of R2 and R3 than R1.Another aromaticity investigation that can be assessed through magnetic criteria is the NICS 46,47 .That can be defined as the negative value of the absolute magnetic shielding calculated at some selective points.The more negative the NICS values, the more aromatic the ring.The NICS(0) and NICS(1) values are computed at the center and 1A˚ above the aromatic ring, respectively.Aromatic systems possess negative NICS values, as that indicates the presence of an induced diatropic ring current.Positive NICS values refer to paratropicity and are present in antiaromatic systems.Very small NICS values refer to nonaromatic systems.The NICS indices calculated as a single point at the B3LYP/6-31 G(d,p) level in the gas phase are given in Table 6.According to previous studies of NICS calculation 46,47,[77][78][79] , the five-membered ring is more aromatic than the six-membered ring followed by the seven-membered ring.As can be seen from Table 6, the pyrrole ring is more aromatic than benzene, and the aromaticity is less in the seven-membered ring.In the gas phase, the thiophenyl ring of thiol exhibits a larger negative NICS(1)zz value (-18.96ppm) than thione (NICS(1)zz (-14.57)).In the presence of a polar solvent like ethanol, very minute variations in magnetic indices of aromaticity are seen.Due to the absence of IHB production and an increase in π-electron delocalization, the investigated rotamers exhibit somewhat larger NICS(1)zz values of the thiophenyl ring than thiol.www.nature.com/scientificreports/NMR analysis.The high correlation between the experimental 21 and theoretical 22 NMR chemical shifts of 2-(2-hydroxyphenyl)-1-azaazulene (2OHPhAZ) at the B3LYP/6-31G(d,p) level that has been previously recorded 22 imparts trust in the computational techniques utilized.The calculated 1 H and 13 C NMR chemical shifts δ(ppm) values for the 2-(2-Mercaptophenyl)-1-azaazulene tautomers and rotamers calculated in CHCl 3 at the B3LYP/6-31G (d,p) level of theory are listed in Table 7.For the 1 H NMR chemical shift, the hydrogen atoms bonded to carbon atoms in aromatic rings, ranging from 6 to 8 ppm [80][81][82] .The thione NH1 proton and the thiol SH1 proton have high chemical shifts (high downfield shift/low-field proton signal) at δ 18.4 and 13.6 ppm, respectively.While for rotamers, the SH 1 proton was assigned at lower chemical shifts in the range of δ 3.9 to 4.6 ppm.This is due to the presence of HB in the thiol and thione structures and the absence of HB in the structures of rotamers.The azaazulenyl ring and phenyl ring protons in the studied compounds showed a low chemical shift in the range of 7.2 to 8.9 ppm.The computed 13 C NMR chemical shifts for 2-(2-Mercaptophenyl)-1-azaazulene tautomers and rotamers are in the region 104.0-173.8ppm (in CHCl 3 ).The highest chemical shift values given in thiol are C2 (163.4 ppm) and C2' in thione (173.8 ppm).This low field shift is due to the fact that C2 and C2' are bonded to the high electronegativity N and S atoms and the high polarity of the NH and C = S groups in thiol and thione, respectively.Similarly, C8a has a high chemical shift (greater than 140 ppm) for thiol and all the studied rotamers by connecting to the N atom.The highest chemical shifts of C2' and H1 in the studied thione correlated well with the highest chemical shift of C2' and H1 in the keto of 2OHPhAZ ( δ 171.3 and 19.4 ppm, respectivel) 22 .The calculated chemical shift of the thiocarbonyl group (C = S) (173.8 ppm) is in good agreement with the reported Table 7. 13 C and 1 H NMR chemical shifts (in ppm) calculated for C2' and H1-S/N at the B3LYP/6-31G(d,p) level of theory of 2-(2-Mercaptophenyl)-1-azaazulene tautomer and rotamers using the GIAO method in CHCl 3 .www.nature.com/scientificreports/chemical shifts of 169.9 and 177.7 ppm for 3-Hydroxy-2-methyl-4-p-pyridinethione 83 and 5-(3-pyridyl)-4H-1,2,4-triazole-3-thione 84 , respectively.
UV-Vis spectral analysis.We have found that our previous calculation of 2OHPhAZ 22 using TDDFT-PBE/SMD with a 6-311+G(d,p) basis set level exhibits better quantitative agreement regarding the first and second maximum excitation peaks with the available experimental data 21 than other levels 22 .The calculated UV-vis parameters of the intense peaks in acetonitrile for the studied tautomers and rotamers with the TD-PBE method are presented in Table 8.The calculated UV-vis spectra for all the studied compounds are given in Fig. 4. The thione form is accompanied by a large red shift (extended to 1000 nm), followed by the thiol form (extended to 700 nm), and the lower shift has been found with the rotamers (extended to 600 nm), as depicted in Fig. 4. As can be seen in Table 5, the energy differences between HOMOs and LUMOs are 2.98 eV (for thiol), 1.87 eV (for    www.nature.com/scientificreports/thione), 3.24 eV (for R1), 2.33 eV (for R2) and 2.23 eV (for R3) in the gas phase.As expected from the lower E g value of thione, the thione maxima in the electronic absorption spectra are shifted bathochromically by 16, 20, 22 and 47 nm in comparison with thiol, R1 , R2, and R3, respectively.The strong electronic absorption of thione is attributed to the HOMO-3 to LUMO and HOMO-2 to LUMO + 1 transitions.The transitions from to HOMO-2 to LUMO and LUMO + 1 and HOMO-1 to LUMO + 1 cause maximum absorption peak for thiol, which appears at a lower wavelength than thione.However, the strong electronic absorption of the studied rotamers has a different contribution, as shown in Table 8.The NTOs for the highintensity excited states of the investigated systems are shown in Fig. 5 to analyze the nature of absorption.The occupied and unoccupied NTOs are referred to as "hole" and "particles" transition orbitals, respectively.The NTOs generally give a simpler description of the excited state than the canonical orbitals.As displayed in Fig. 6, where the canonical orbitals were used, the dominant transitions are π-π* for the excitations, with some contribution from n-π* excitation.This makes the analysis of excitations cumbersome.However, as depicted in Fig. 6, the hole NTOs contributing to the illustrated band in Fig. 6 and Table 8 of all studied structures are delocalized over the whole molecular skeleton, while the particle NTOs are mainly delocalized over either azaazulene or benzene rings.This suggests π-π* excitation.

Nonlinear optical (NLO) analysis.
It is known that heterocyclic compounds have enhanced charge transfer, which could lead to large nonlinear optical properties.In this section, we will test this information and provide the scientific community any basic information about the hyperpolarizability of these tautomer's and rotamers.Non-linear optical properties (NLO) are the ability of any compound to convert light [with an intense electric field (LASER)] of a longer wavelength into light of a shorter wavelength.One of the non-linear optical phenomena is the second harmonic generation (SHG), where intense light of a longer wavelength is converted to half of the incident value upon absorption by the non-linear optical material.
Most applications of single crystals of any nonlinear material are evident in the fields of semiconductors, infrared detectors, solid state lasers, photosensitive materials, and crystalline thin films for microelectronics [85][86][87] .The investigation of the relationship between the electronic structure and NLO parameters of the studied compounds is calculated theoretically using ωB97XD/6-311++G(2d,2p) and B3LYP/6-311++g(2d,2p) (Table 9 and S2, respectively).Total static dipole moment (µ), the mean polarizability, the anisotropy of the polarizability ∆, α the mean first-order hyperpolarizability (β), the hyper-Rayleigh scattering (β HRS ) and the depolarization ratio (DR) of the studied compounds are listed in Table 9.In this study, P-nitro aniline (PNA), a standard prototype molecule used in NLO studies, was chosen as a reference as there were no experimental values for the NLO properties of the studied compounds.The values of, α in Table 9 show that the order of increasing α with respect to PNA is: R2, Thiol, and R1 are ~ 3.5 and 3 times higher than PNA, respectively, whereas rotamer R3, and Thione is ~ 4 times higher than the standard PNA at ωB97XD/6-311 + +G(2d,2p).The calculated first order hyperpolarizability β of p-nitroacetanilide (PNA) is 15.5 × 10 -30 esu as reported by T. Gnanasambandan et al. [88][89][90] The analysis of the β parameter shows that Thiol, R2, R3, and R1 are ~ 2 times higher than PNA, while rotamer Thione is ~ 2.5 times higher than the reference at ωB97XD/6-311 + +G(2d,2p).Furthermore, for the investigated compound, the lowest value of DR, and the highest value of β HRS confirm short bond length, indicating increased selectivity.Therefore, the studies tautomer's and rotamers show promising optical properties.

Conclusions
Five tautomers and rotamers (thiol, thione, R1, R2, and R3) of 2-(2-Mercaptophenyl)-1-azaazulene were studied using density functional theory (DFT).At the B3LYP/6-31G(d,p), the geometrical structural parameters and vibrational frequencies have been discussed.According to the results, B97XD/6-311++G(2d,2p) performs better than CCSD(T)/6-311++G(2d,2p).A very good agreement is seen when compared to the similar molecule 2-(2-hydroxyphenyl)-1-azaazulene.All levels of calculations show that the thiol form is the most stable tautomer in the gas phase.The order of stability is reversed in ethanol as a polar solvent, whereas thione becomes the most stable tautomer.It can be seen that the thione NH1 proton and the thiol SH1 proton have high chemical shifts at δ 18.4 and 13.6 ppm, respectively.HOMA and NICS results for aromaticity provide a good match with the relative energies and the structure of the thiol-thione tautomer in both gas and ethanol solvent.The UV-vis absorption spectra of the studied compound have been calculated using the TD-PBE method with a 6-311+G(d,p) basis set for 2-(2-Mercaptophenyl)-1-azaazulene.UV-vis computations indicate that the thione form has a significant red-shift (extended to 1000 nm), followed by the thiol form (extended to 700 nm), and that the rotamers have the lowest red-shift (extended to 600 nm).NTOs are used to indicate the π − π* nature of the transitions.In comparison to the PNA molecule, the studied molecules offer good benefits in technology-related applications.

Figure 3 .
Figure 3. Molecular ESP surfaces of the investigated structures at the B3LYP/6-31G(d,p) in gas phase.

Figure 6 .
Figure 6.Natural transition orbitals (NTOs) for the excitation with significant and small but non-negligible oscillation strengths for the studied structures at the PBE/6-311+G(d,p) level with solvent effects of acetonitrile through SMD.The displayed occupied (holes) and unoccupied (electrons) NTO pairs are the only that have contributed more than 50% to each excited state (λ is eigenvalues of the pairs). 2

Table 5 .
Global chemical descriptor (eV) of the studied structures at B3LYP/6-31G (d,p) in the gas phase and ethanol.

Table 6 .
The NICS(1) zz (in ppm) and HOMA index of the studied tautomers and rotamers calculated in gas phase and ethanol at the B3LYP/6-31G(d,p).Five-and seven-membered rings of azaazulene and phenyl rings are presented in Fig.2.