CCComputational Chemistry2332-5968Scientific Research Publishing10.4236/cc.2021.93010CC-110514ArticlesChemistry&Materials Science Investigation of the Antioxidant and UV Absorption Properties of 2-(2’-hydroxy-5’-methylphenyl)-benzotriazole and Its Ortho-Substituted Derivatives via DFT/TD-DFT NumbonuiStanley Tasheh1AymardDidier Tamafo Fouegue2JuliusNumbonui Ghogomu3*Department of Chemistry, Higher Teacher Training College Bertoua, University of Ngaoundéré, Bertoua, CameroonResearch Unit of Noxious Chemistry and Environmental Engineering, Department of Chemistry, Faculty of Science, University of Dschang, Dschang, CameroonDepartment of Chemistry, Faculty of Science, The University of Bamenda, Bamenda, Cameroon0306202109031611967, June 202111, July 2021 14, July 2021© Copyright 2014 by authors and Scientific Research Publishing Inc. 2014This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

The demand and pursuit of chemical entities with UV filtration and antioxidant properties for enhanced photoprotection have been driven in recent times by acute exposure of humans to solar ultraviolet radiations. The structural, electronic, antioxidant and UV absorption properties of drometrizole (PBT) and designed ortho-substituted derivatives are reported via DFT and TD-DFT in the gas and aqueous phases. DFT and TD-DFT computations were performed at the M062x-D3Zero/6-311++G(d,p)//B97-3c and PBE0-D3(BJ)/def2-TZVP levels of theory respectively. Reaction enthalpies related to hydrogen atom transfer (HAT), single-electron transfer followed by proton transfer (SET-PT), and sequential proton loss electron transfer (SPLET) mechanisms were computed and compared with those of phenol. Results show that the presence of -NH2 substituent reduces the O-H bond dissociation enthalpy and ionization potential, while that of -CN increases the proton affinity. The HAT and SPLET mechanisms are the most plausible in the gas and aqueous phases respectively. The molecule with the -NH2 substituent (PBT1) was identified to be the compound with the highest antioxidant activity. The UV spectra of the studied compounds are characterized by two bands in the 280 - 400 nm regions. Results from this study provide a better comprehension antioxidant mechanism of drometrizole and present a new perspective for the design of electron-donor antioxidant molecules with enhanced antioxidant-photoprotective efficiencies for applications in commercial sunscreens.

Antioxidant Drometrizole Density Functional Theory Ultraviolet Absorption Hydrogen Atom Transfer Single Electron-Proton Transfer and Sequential Proton Loss Electron Transfer Mechanisms1. Introduction

Biological systems depend on solar radiation for energy to sustain life. Unfortunately, chronic exposure to solar ultraviolet radiation (UVR) is responsible for health problems such as premature aging, erythema, skin cancers and sunburn [1] [2] [3] [4] [5] . Melanin (naturally occurring UV-absorbing pigment in the human skin) provides the first line of defense to the harmful effects of UVR but does not provide enough shield against high levels of UVR. To this effect, organic sunscreens are frequently used as photoprotectors [4] . Recently, attention has been shifted towards the use of sunscreens which in addition to the screening of detrimental UVR equally have photoprotection constituents [6] . Antioxidants efficiently serve this purpose, since they scavenge UV-induced free radicals [2] [7] . Flavonoids, polyphenols, resveratrol, etc., have been explored for such purposes. To qualify as a potential organic molecule for sunscreen applications, the molecule must be able to block UVR and scavenge resultant free radicals from UVR exposure [4] [8] . Knowledge of sunscreens comprising UV filters and antioxidant ingredients serve as a novel photoprotective strategy and research in this field is ongoing [3] [4] [6] [9] .

2-(2’-hydroxy-5’-methylphenyl)-benzotriazole (PBT) commonly known as drometrizole (see Figure 1) is considered as a sunscreen ingredient and is utilized as a UV absorber and stabilizer [1] [10] [11] . Furthermore, it is known to

have antioxidant and stabilizing effects and has been approved as an indirect food additive [12] . Reference [1] reported that PBT is a potent Exited state intramolecular proton transfer (ESIPT)-mediated Non-linear optical (NLO) molecular switch which could be used as a chromophore. Despite all the studies on drometrizole from literature, there is no comprehensive theoretical study on the structure-antioxidant relationship as well as the antioxidant mechanism of drometrizole. Additionally, the presence of some substituents on PBT may enhance its antioxidant properties. Interestingly, electron-donating groups (EDGs) and electron-withdrawing groups (EWGs) are known to significantly affect the antioxidant activity of a phenolic antioxidant [13] . This insight has been a prime aid for the design and synthesis of new and more effective antioxidants [14] . Three major mechanisms to scavenge free radical by antioxidants have been recognized [15] [16] . These include hydrogen atom transfer (HAT), single-electron transfer followed by proton transfer (SET-PT), and sequential proton loss electron transfer (SPLET). The antioxidant activity of the phenolic antioxidants can be assessed thermodynamically by several physicochemical parameters such as O-H bond dissociation enthalpy (O-HBDE), proton dissociation enthalpy (PDE), proton affinity (PA), ionization potential (IP) and electron transfer enthalpy (ETE).

As can be seen in Figure 1, the structural features of drometrizole (pi-conjugation, presence of hydroxyl group) and its above-mentioned properties warrant it to be a potent antioxidant-UV filter. Such a study is conveniently achieved via the use of density functional theory (DFT) and time-dependent DFT (TD-DFT) methods. DFT has been utilized by several authors to predict the relationship between the structure-activity of a compound and also for the design of novel potential antioxidants [17] .

Within the framework of the design of more UV filter-antioxidant-based molecules, this work had as the main objective, to evaluate the thermodynamics of the antioxidant mechanism and UV-absorption properties of PBT using the DFT and TD-DFT methods. In addition, the effects of an electron-donating substituent (EDS) and an electron-accepting substituent (EAS) (see Figure 1) on the antioxidant activity and UV spectra of PBT were assessed. The effects of these substituents were explored with the intention of improving their anti-oxidant efficiencies. A comprehension of why a mechanism is favored over another requires in due course some kinetic studies to confirm speculations drawn from thermodynamic studies. However, this is not the main aim of this study. Here, we set out to specifically: 1) Find drometrizole derivatives with high antioxidant activity from a theoretical point of view; 2) Identify the substituent that are the most effective in reducing the BDE, IP, and PA; 3) Find the preferred mechanisms in the gas and solvent phases from a thermodynamic perspective; 4) Compute the UV spectra of these molecules. To achieve these, the physiochemical descriptors including O-HBDE, IP, PDE, PA and ETE were computed to examine the HAT, SET-PT and SPLET mechanisms. Also, the frontier molecular orbitals (FMO) and energies of the studied molecules were computed. Furthermore, the transition wavelengths λ (nm), oscillator strengths (f), the major components of the molecular orbital (MO) transitions and configuration interaction (CI) were computed to elucidate the UV filtering ability of the investigated compounds. This study will serve as a base to create an awareness of the photoprotective properties of the compounds studied possible design and applications as sunscreen products.

2. Computational Details

All calculations in this study were done with the ORCA 4.2.1 program suite [18] [19] . Initial guess geometries were conceived by means of the Avogadro 1.2.1 visualization program [20] . For all reported calculations, a tight self-consistent field (SCF) criterion was utilized with a grid6 sized numerical quadrature. The restricted Kohn-Sham (RKS) formalism and its unrestricted counterpart (UKS) were adopted for closed and open-shell molecules (i.e. for radicals and radical cation species) respectively. The geometry of each drometrizole molecule (ArOH), their corresponding radicals (ArO•), cationic radicals (ArOH•+) and anions (ArO−) were optimized in the gas phase with no constraints, using the low-cost B97-3c method. This method was used thanks to its efficient performance in geometry optimization and frequency calculations of organic, inorganic, and organo-metallic systems [21] . Only real frequencies were obtained from vibrational frequency calculations at the same level of theory as that used for geometry optimization, confirming the optimized geometries are minima on their potential energy surfaces. In all cases, no spin contamination was found.

The electronic energies from the B97-3c method were enhanced via single-point (SP) computations at a higher level of theory using the Minnesota exchange-correlation (XC) functional, M06-2X [22] , together with the 6-311++G(d,p) basis set. The M06-2X functional was used owing to its reliability in calculating thermodynamic parameters [22] . The D3 dispersion correction with zero damping (D3Zero) was employed to justify long-range dispersion interactions in all SP calculations with the M06-2X functional. The resolution-of-identity (RI-J) approximation [23] and the Chain-of-Spheres (COSX) approximation [24] , establishing the RIJCOSX approximation, were utilized with the suitable auxiliary basis sets to fasten the SP computations. Remarkably, the loss of accuracy from applying the RIJCOSX approximation is insignificant [25] . Reaction enthalpies were computed M06-2X/6-311++G(d,p) level of theory, by making use of the electronic energies obtained at this level with vibrational frequencies obtained from the B97-3c method, with the aid of the Kisthelp software.

To include solvent effects on the computed properties, a continuum solvation method termed Solvation Model Density (SMD) was employed [26] via SP calculations on the gas-phase optimized B97-3c geometries. SMD was used owing to its universal and applicable nature for any charged or uncharged solutes, with relatively low errors [27] . Water (dielectric constants, ε = 78.35) was used to mimic the aqueous environment.

Frontier molecular orbital energies of the studied molecules were obtained from the SP calculations. The HOMO-LUMO energy gaps (ΔE) of the studied compound were calculated based on the energies of the highest and lowest occupied molecular orbitals (EHOMO and ELUMO, respectively) as follows:

Δ E = E LUMO − E HOMO (1)

To compute the absorption spectra of the studied compounds, vertical electronic excitations from the minima of the ground-state structures were computed using the PBE0 functional [28] with the def2-TZVP Ahlrichs triple-ζ basis set [29] , employing the time-dependent density functional theory (TD-DFT) as implemented in the ORCA software. PBE0 has been utilized because it produces computed values which corroborate their experimental counterparts [1] . The def2-TZVP basis set was employed since it is appropriate to compute electronic properties [29] . To account for long-range dispersion interactions, the Grimme’s atom-pairwise dispersion correction with Becke-Johnson damping (D3BJ) [30] was employed. Using the TD-DFT method, electronic absorption properties (like the transition wavelengths λ (nm), oscillator strengths (f), the major components of the molecular orbital (MO) transitions and configuration interaction (CI)) of the studied compounds on the basis of the B97-3c optimized ground state geometries were obtained.

2.1. Antioxidant Potential

As stated before now the antioxidant capacity of phenolic compounds has been elucidated through several mechanisms. The commonly used ones are hydrogen atom transfer (HAT), single electron transfer followed by proton transfer (SET-PT), and sequential proton loss electron transfer (SPLET) [31] . These mechanisms individually offer an approach on how these molecules respond in the presence of a radical [15] [32] . It is worth noting that the previously mentioned mechanisms are liable to happen in parallel, however, at different rates. The above-mentioned mechanisms are presented in Figure 2 and have been briefly discussed below.

2.1.1. Hydrogen Atom Transfer (HAT) Mechanism

This mechanism is characterized by the capacity of a free radical (R•) to remove a hydrogen atom from the antioxidant (ArOH), yielding the radical form of phenolic antioxidant (ArO•).

ArOH → ArO • + H • (2)

2.1.2. Single Electron Transfer-Proton Transfer (SET-PT) Mechanism

The SET-PT mechanism is a two-step reaction. The first step is ruled by the ability of an antioxidant to transfer an electron, giving rise to a cationic radical (ArOH•+).

ArOH → ArOH • + + e − (3)

In the second step, the ArOH•+ previously formed is deprotonated into a phenolic radical (ArO•).

ArOH • + → ArO • + H + (4)

The numerical descriptors resulting from this mechanism are IP and PDE for the first and second reactions respectively.

2.2. Sequential Proton Loss Electron Transfer (SPLET) Mechanism

Like the SET-PT mechanism, the SPLET mechanism also consists of a two-step reaction. The first step corresponds to the deprotonation of the phenolic antioxidant, yielding the antioxidant anion (ArO−).

ArOH → ArO − + H + (5)

The second step corresponds to the transfer of an electron form ArO−, resulting in the formation of the phenoxyl radical (ArO•) of the antioxidant.

ArO − → ArO • + e − (6)

The numerical descriptors resulting from this mechanism are PA and ETE for the first and second steps respectively.

To investigate the likelihood of these mechanisms to occur, related thermodynamic parameters comprising: O-H bond dissociation enthalpy (O-HBDE), ionization potential (IP), proton dissociation enthalpy (PDE), proton affinity (PA), and the electron transfer enthalpy (ETE) were assessed according to Equations (7-8) as:

BDE = H ( BTO • ) + H ( H • ) − H ( BTOH ) (7)

IP = H ( BTOH • + ) + H ( e − ) − H ( BTOH ) (8)

PDE = H ( BTO • ) + H ( H + ) − H ( BTOH • + ) (9)

PA = H ( BTO − ) + H ( H + ) − H ( BTOH ) (10)

ETE = H ( BTO • ) + H ( e − ) − H ( BTO − ) (11)

where H ( BTOH ) , H ( BTO • ) , H ( H • ) , H ( BTOH • + ) , H ( BTO − ) , H ( e − ) , and H ( H + ) are enthalpies of the neutral molecule (antioxidant), neutral radical, hydrogen atom, radical cation, anion, electron, and proton, respectively. The gas phase H ( e − ) , and H ( H + ) values for this study were obtained from [33] [34] [35] , while those of H ( e − ) , and H ( H + ) in the aqueous phase from [32] . The H ( H • ) in the gaseous and aqueous phases for this study were computed at the M06-2x/6-311++G(d,P) level of theory.

3. Results and Discussion3.1. Investigating the Suitability of the Optimized Structures

In order to assess the suitability of the B97-3c method and the accuracy of the optimized structures, a selection of gas-phase calculated bond lengths and angles of PBT have been likened to the corresponding experimental values obtained elsewhere [36] . The computed bond lengths and angles of PBT are enumerated in Table S1 of the electronic Supplementary Information (ESI), alongside their corresponding experimental values. It can be perceived from Table S1 that the theoretical bond lengths and angles agree with experimental data, indicating the suitability of the method used. For better comparison of these values, Pearson’s linear regression was performed. Correlation coefficients (R2) of 0.99 and 0.93 were found for the bond angles and bond lengths respectively. From the high R2 values (very close to unity), the calculated bond lengths and bond angles of PBT agree perfectly with the experimental values, indicating the reliability of the optimized structures.

3.2. Geometrical Analyses

The fully optimized B97-3c geometries of the neutral forms of PBT, PBT1 and PBT2 in the gaseous phase are shown in Figure 3. Those of the corresponding radical, cationic radical and anionic forms are included in the accompanying ESI as visualized with the ChemCraft v1.7 software [37] . The Cartesian coordinates for all these structures are provided in the ESI. Selected geometric parameters of these molecules are listed in Table 1 (see Figure 1 for atom numbering). The selected parameters presented in Table 1 originate from the 4-methoxy phenolic ring, as changes caused by the ED and EASs at the ortho position (C17 atom) of this ring, may yield some loss of planarity of PBT. As can be observed from the values presented in Table 1 it is obvious that the type of substituents (i.e. EDS and EASs) on C17 of the 4-methoxy phenolic ring does not have a significant effect on the selected geometrical parameters. Also, it is inferred from the bond

Selected geometric parameters of PBT and ortho-substituted optimized with the B97-3C method gaseous phase
ParametersPBTPBT1PBT2
Bond lengths (Å)
N11-N121.3411.3421.341
N12-N131.3261.3251.325
C14-C151.4081.3991.408
C15-O161.3431.3511.333
C15-C171.3961.4101.409
C17-C181.3771.3861.393
C18-C191.3981.3961.389
C19-C201.3821.3801.386
Bond angles (˚)
N11-N12-N13116.3116.2116.4
C14-C15-C17117.4118.5117.1
C15-C17-C18121.3119.0120.9
C17-C18-C19121.5122.0121.5
C18-C19-C20117.8118.9117.9
O16-C15-C17118.0116.4118.2
Dihedral angles (˚)
O16-C15-C14-C20−179.8179.5−179.8
O16-C15-C17-C18179.8−179.4179.8
C14-C15-C17-C180.1−0.20.1
C14-C20-C19-C180.20.40.3

length values listed in Table 1 that for each molecule investigated, the numerical values are close to 120˚, suggesting that the atoms in the 4-methoxy phenol ring are nearly sp2 hybridized and are therefore very likely to adopt a highly planar arrangement. That notwithstanding, the planarity of these molecules is likely to be significantly reduced owing to free-rotation around the N12-C14 bond. Nevertheless, the rotation is highly hindered by the intramolecular hydrogen bonds O16-H22⋯N11 which also stabilize the molecules. Also, from the torsional angles in the 4-methoxy phenolic ring listed in Table 1 are close to 0˚ or 180˚. On the basis of the ongoing, the optimized molecular geometries of the studied molecules are highly planar and there is no significant effect with the presence of the EDS and EASs.

3.3. Antioxidant Properties of the Studied Molecules

In this study, the relative values are referred to those of phenol, generally designated as the zero compound [38] .

3.3.1. Evaluation of the HAT Mechanism

As earlier mentioned, O-HBDE is the numerical parameter related to the HAT mechanism. It characterizes the stability of the corresponding hydroxyl group. Lower O-HBDE values designate lower stability and the ease to which the equivalent O-H bond is broken, giving rise to a higher antioxidant activity [39] . The O-HBDE values of the studied molecules were calculated from Equation (2) and the results are presented in Table 2 in the gas and aqueous phases. As can be seen from Table 2, the O-HBDE values of the molecules under study are larger in the gaseous phase (by about 21 kJ/mol), implying that H-atom abstraction is more difficult in the gaseous phase.

It is also found in Table 2 that the difference between the O-HBDE values of PBT and that PBT1 is approximately 42 and 46 kJ/mol, while that with PBT2 is about −6 and −9 kJ/mol in the gas and aqueous phases respectively. This means that the presence of -NH2 lowers the O-HBDE values, while that of -CN increases this value as compared to that of PBT. This variation is attributed to the fact

Computed O-H bond dissociation energies (O-H<sub>BDE</sub>) (in kJ/mol) of the studied compounds obtained at the M06-2x-D3(BJ)/6-311++G(d,p) level of theory
HAT
Molecular speciesGasWater
Phenol372.2 (0.0)a376.7 (0.0)
PBT398.0 (25.9)376.9 (0.2)
PBT1356.1 (−16.0)331.1 (−45.6)
PBT2404.0 (31.9)385.6 (8.9)

aRefer to phenol for the relative O-HBDE values (in the parentheses).

that EASs in the ortho-position stabilize the parent molecule and destabilize the radical, causing a rise in the O-HBDE values. Contrariwise, EDSs in ortho-position have an opposite effect, leading to a reduction in the O-HBDE values. The O-HBDE values of the studied compounds were compared to that of phenol at the same theoretical level, and it is observed that PBT1 has a lower value (by about 16 and 45 kJ/mol in both study environments respectively). This clearly shows that PBT1 has a better antioxidant potential than phenol. In all investigated media, the lowest O-HBDE was recorded for PBT1, while the highest for PBT2. This implies that H-atom abstraction is easier in PBT1 and more difficult in PBT2. On the basis of the O-HBDE and the HAT mechanism, the antioxidant activity decreases as follows: PBT1 > PBT > PBT2 in all studied medium.

Spin density distributions

The spin density distribution is a well-grounded parameter that offers a better illustration of the reactivity of a radical system. Lower spin density distribution values yield a higher delocalization of the spin density in a given radical. The computed spin densities on the oxygen atom of the phenoxy radicals (PBT, PBT1 and PBT2) in the gas and aqueous phases are presented in Table 3. Results from this table indicate that the spin density on the oxygen atom of the PBT1 radical is lower than those of the PBT and PBT2 radicals in both studied environments, implying that radical PBT1 has the highest stability and consequently is the least reactive in radical reactions. This result is allied to the lowest O-HBDE values of PBT1 from Table 2, whereby PBT1 is the best antioxidant molecule.

3.3.2. Evaluation of the SET-PT Mechanism

As stated earlier, the ionization potential (IP) and proton dissociation enthalpy (PDE) are the descriptors associated with the two-step SET-PT mechanism. Molecules with low IP values are more susceptible to ionization and have potent antioxidant properties. The lower the IP value for a given molecule, the higher its electron donating ability and higher antioxidant activity [40] . For this study, the IP and PDE values of the molecules studied were computed from Equations (3) and (4) respectively and the values alongside that of phenol in the vacuum and aqueous phase are listed in Table 4. In this table, it is evident that the IP values of the compounds under study are less than that of phenol in both studied environments. Implying that the investigated molecules have a greater tendency to give out electrons than phenol and consequently have better antioxidant properties than the latter. From this table, the difference between the IP value of PBT

Spin density values on the O-radical of the studied compounds computed at the at the M06-2x-D3(BJ)/6-311++G(d,p) level of theory
Molecular speciesGasWater
PBT0.39360.3252
PBT10.32640.2609
PBT20.36990.3068
 Computed ionization potential (IP) and proton dissociation energies (PDE) (in kJ/mol) of the studied compounds obtained at the M06-2x-D3(BJ)/6-311++G(d,p) level of theory
Molecular speciesSET-PT
IPPDE
GasWaterGasWater
Phenol825.8 (0.0)b521.0 (0.0)857.5 (0.0)c18.7 (0.0)
PBT759.5 (−66.3)509.4 (−11.60)949.6 (92.1)30.5 (11.8)
PBT1693.7 (−132.1)443.3 (−77.5)973.6 (116.1)50.9 (32.2)
PBT2805.0 (−20.8)539.4 (18.35)910.2 (52.7)9.3 (−9.4)

bRefer to phenol for the relative IP and PDE values (in the parentheses); cRefer to phenol for the relative PDE values (in the parentheses).

and PBT1 is 66 kJ/mol in both media while that of PBT and PBT2 are −45 and −30 kJ/mol in the gas and aqueous respectively. Clearly, from Table 4, the lowest IP value was found to be for PBT1 while the highest for PBT2 in all studied environments. This indicates that PBT1 has the greatest ability to donate electrons while the reverse is true for PBT2. Additionally, Table 4 shows that the presence of the EDS reduces the IP values, while that of the EAS increases this value as compared to that of PBT in both media. The increase in the IP values is attributed to the fact the EAS stabilizes the parent molecule and destabilizes ArOH•+, while the decrease in IP values is due to the fact that the EDS tends to stabilize ArOH•+ and destabilize the parent molecule.

Also, it can be seen in Table 4 that the formation of ArOH•+ for the molecules under investigation is easier in water (by about 255 kJ/mol). This indicates that the formed radical cation, ArOH•+ is more stable in water due to the delocalization of π-electrons. This result agrees with that in [40] , which confirms that electron donation is easier in the solvent, particularly polar medium (water in this case). On the basis of the ongoing, the increasing IP trend is: PBT1 < PBT < PBT2.

PDE describes the potential of the cationic radical, ArOH•+ to transfer a proton to the free radical. Molecules with lower PDE values are more susceptible to proton loss [41] . As can be seen from Table 4, the highest PDE values are found for PBT1, whereas lowest for PBT2. This observation is as a result of the fact that the presence of the EDS in PBT1 stabilizes ArOH•+ but destabilizes the parent molecule, while the presence of the EAS in PBT2 have an opposite effect. These results are in accordance with a previous paper on substituted phenols [42] . From a thermodynamic view point, the first step (IP) is the most important. From the values reported in Table 4, the SET-PT mechanism is not preferred in the gas phase due to the fact that the IP and PDE values are higher than the corresponding aqueous phase values. The significant decrease in both parameters (predominantly the PDE values) in aqueous phase implies the ease in the proton transfer and electron abstraction processes respectively. This observation is similar to that obtained in [31] .

3.3.3. Evaluation of the SPLET Mechanism

As earlier mentioned, the proton affinity (PA) and electron transfer enthalpy (ETE) parameters are linked to the SPLET mechanism. The PA and ETE values of the compounds studied were calculated according to Equations (5) and (6) respectively in the gaseous and aqueous phases together with those of phenol are enumerated in Table 5. The proton loss (characterized by the PA values) is the underlying step in the SPLET mechanism. Lower PA values favor the deprotonation of the drometrizole antioxidants yielding the phenoxyl anion. It can be seen from Table 5 that the PA values for all studied compounds are higher in the gaseous phase than in the aqueous phase, suggesting that the proton loss is highly preferred in aqueous phase. In addition, it can be inferred from Table 5 that the PA values of the molecules under investigation are lower than those of phenol.

Electron transfer (characterized by the ETE values) from the phenoxide ion represents the second step in SPLET mechanism. ETE denotes the capacity of the antioxidant anion to transfer electrons to the free radical. The lower the ETE value for a given molecule, the more active is the resulting phenoxide anion. From Table 5, the ETE values in the gas and aqueous phases follow the same sequence: ETEPBT1 < ETEPBT < ETEPBT2. In aqueous media, the ETE values of the studied compounds are higher than the PA values but enough to account for the second step of SPLET mechanism. It is evident from Table 5 that the lowest ETE value corresponds to PBT1 while the highest to PBT2 in all studied environments. This shows that PBT1 has a greater tendency to transfer electrons to the free radical and the ease with which it forms the corresponding phenoxide anion, while the reverse is the case for PBT2. Also, it is obvious from Table 5 that the presence of the EDS in PBT1 decreases the ETE value, while that of the EAS in PBT2 increases this value as compared to PBT. It can also be seen in Table 5 that the formation of the phenoxide anion for the investigated molecules is easier in the aqueous phase than in the gaseous phase. On the basis of the ETE values, PBT1 is the best anti-oxidant molecule, while PBT2 is the least.

Computed proton affinities (PA) and electron transfer energies (ETE) (in kJ/mol) of the studied compounds obtained at the M06-2x-D3(BJ)/6-311++G(d,p) level of theory
Molecular speciesSPLET
PAETE
GasWaterGasWater
Phenol1453.1 (0.0)d152.3 (0.0)230.2 (0.0)e387.5 (0.0)
PBT1439.1 (−14)149.9 (−2.4)270.0 (39.8)390.1 (2.6)
PBT11432.6 (−20.5)152.9 (0.6)234.7 (4.5)341.3 (−46.2)
PBT21379.3 (−73.8)118.5 (−33.8)335.9 (105.7)430.2 (42.7)

dRefer to phenol for the relative PA values (in the parentheses); eRefer to phenol for the relative ETE values (in the parentheses)

To ease comparisons, the O-HBDE, IP and PA values for all the molecules investigated were plotted on the same axes and are presented on Figure 4(a) and Figure 4(b) for the gas and aqueous phases respectively. It is obvious from Figure 4(a) that the computed PA values of the studied molecules are higher than the equivalent IP and O-HBDE values. This illustrates that in the gas phase, the favored mechanism can be classified thermodynamically as: HAT < SET-PT < SPLET. In the same manner, it can be seen from Figure 4(b) that the PA values of the compounds investigated are remarkably lower than the corresponding O-HBDE and IP values. Therefore, in aqueous phase, the SPLET mechanism represents the most probable reaction pathway from thermodynamic point of view. In the aqueous phase, the trend is: SPLET < SET-PT < HAT. This trend ascertains the importance of this mechanism in a polar solvent as reported elsewhere [43] . Conclusively, among the three mechanisms investigated, the HAT mechanism is thermodynamically favored in the gas phase and SPLET in water.

3.4. Frontier Molecular Orbital (FMO) Analyses

FMOs are important parameters which portray electronic properties as well as reactivity of a molecule. The Highest Occupied Molecular Orbital (HOMO) and Lowest Unoccupied Molecular Orbital (LUMO) distribution can equally be utilized to describe the antioxidant capabilities of a molecule [44] . Knowledge of the HOMO distribution can be used to deduce the probable site for attack by free radicals in the molecule.

As can be seen from Table 6, the HOMO of PBT1 is mainly localized on the 4-methylphenolic moiety, while those PBT and PBT2, are appreciably distributed over their entire molecular structures (i.e. on the 4-methylphenolic and benzotriazole moieties). This suggests the probable sites for radical attack.

On the basis of FMO and the Fukui theories, FMO energies (i.e. EHOMO, and ELUMO are vital to predict chemical reactivity and stability of molecules. The EHOMO

HOMO distributions of the studied molecules obtained at the M062x-D3Zero/ 6-311++G(d,p) level of theory
SpeciesGas phaseAqueous phase
PBT
PBT1
PBT2

is another parameter to measure the electron-donating ability of antioxidants. The higher the EHOMO value, the greater the tendency of a molecule to donate electrons is and consequently a higher antioxidant is potential. The HOMO-LUMO energy gap, ΔE in turn provides information on the chemical activity of the molecule. The computed EHOMO, ELUMO and ΔE values of the studied molecules in the gas and aqueous phases are summarized in Table 7.

It can be seen in Table 7 that the presence of the EDS (-NH2) in PBT1 increases EHOMO (i.e. it becomes less negative), while the presence of EAS (-CN) in PBT2 reduces EHOMO (i.e. it becomes more negative) as compared to that of PBT. Consequently, drometrizoles with strong EDSs are better electron-donors, i.e., they easily partake in the SET-PT mechanism than those with EASs. Accordingly, the electron-donating ability of the studied molecules increases in the order: PBT2 < PBT < PBT1. From its relatively high EHOMO values, PBT1 is the most chemically reactive species among the molecules studied. This is attributable to the presence of the -NH2 group at the ortho-position of the -OH group, which acts as an EDS by means of its positive inductive (+I) effect which activates the 4-methoxy phenolic ring, making it more reactive. The IP trend is consistent with the data presented in Table 7, confirming that PBT1 is the best electron-donor among the investigated compounds. From this trend, we conclude that

Calculated E<sub>HOMO</sub>, E<sub>LUMO</sub> and ΔE (eV) of drometrizole and its ortho-substituted derivatives computed at the M06-2x-D3(BJ)/6-311++G(d,p) level of theory in gas and aqueous phases
Molecular speciesGasAqueous
SubstituentsNotationsEHOMOELUMOΔEEHOMOELUMOΔE
-HPBT−7.3891−1.42685.9623−7.4763−1.41616.0602
-NH2PBT1−6.9245−1.35105.5735−7.0732−1.38795.6853
-CNPBT2−7.9124−1.85546.0570−7.7558−1.61436.1415

EHOMO can be rapidly used to estimate the reaction enthalpies of the first step of the SET-PT mechanism. These dependences can be used to select appropriate substituents for the synthesis of novel drometrizole-based derivatives with improved antioxidant activity [45] .

ΔE provides valuable information on the reactivity and stability of molecules [46] [47] . The ΔE values of species investigated in this work were calculated according to Equation (1). It is obvious from the values reported in Table 7 that PBT1 has the lowest ΔE values while PBT2 has the highest. Accordingly, the transfer of electrons from HOMO to LUMO in PBT1 is relatively easier than in the other studied compounds. Based on the ΔE values of the studied compounds, the order of decreasing ΔE is: ΔEPBT2 > ΔEPBT > ΔEPBT1. This trend is consistent with those obtained from the EHOMO values. PBT1 was predicted to be the most reactive and consequently, has the best antioxidant activity amongst the studied molecules. This result is in perfect agreement with the predicted antioxidant activity as presented in Table 3.

3.5. Computed Electronic Absorption Spectra of the Studied Compounds

To assess the ability of PBT to screen UV radiations and the effect of ED and EASs, their electronic absorption spectra were simulated at the TD-DFT/PBE0-D3(BJ)/ def2-TZVP levels of theory in gas and aqueous (water solvent) phases. The spectral data obtained is presented in Table 8. Listed in Table 8 are the vertical absorption transition wavelengths (λabs in nm), the oscillator strengths (f), the configuration interaction (CI) or molecular orbital percentage contributions, and the prominent electronic transitions. As can be seen in Table 8, two absorption bands with the maximum band centered around 334, and the shoulder band around 282 nm are observed in the absorption spectrum of PBT in both the gas and water solvent. These bands correspond to S0→S1 and S0→S2 transitions respectively. Interestingly, these values are very close to the experimentally obtained absorption spectrum of PBT obtained in [1] , which shows two intense bands centered around 330 and 280 nm. This result confirms the suitability of the level of theory used. The relatively high oscillator strengths of the calculated electronic transitions indicate their high intensities. The observed transitions: S0→S1, and S0→S2 originate from the molecular orbitals: HOMO (H), HOMO-2 (H-2)

Transition wavelengths λ (nm), oscillator strengths (f), the major components of the molecular orbital (MO) transitions and configuration interaction (CI) for PBT and derivatives computed at PBE0-D3(BJ)/def2-TZVP level of theory
SpeciesMediaTransitionλabs (nm)fMO contributionCI (%)
PBTGasS0→S13400.3587H→L96.42
S0→S22830.4250H-2→L76.78
WaterS0→S13330.3713H→L95.88
S0→S22840.1431H-1→L88.87
PBT1GasS0→S13950.0411H→L98.54
S0→S23130.7123H-1→L97.71
WaterS0→S13870.0367H→L98.52
S0→S23100.7145H-1→L97.67
PBT2GasS0→S13360.4074H→L95.72
S0→S22880.0676H-1→L94.36
WaterS0→S13270.4735H→L94.92
S0→S22940.0331H-1→L96.28

and HOMO-1 (H-1), respectively, and terminate in the molecular orbital LUMO (L) in each case. In the gas phase, the maximum absorption is due to the transition from H→L and the shoulder peak due to transitions from H-2→L. In water, the maximum absorption is due to the transition from H→L and the shoulder peak due to the H-1→L transitions. It is evident from the results obtained that PBT has an absorption in the region 283 - 340 nm, and can effectively absorb the dangerous UV A and UV B radiations from sunlight. This result is in accordance with that obtained elsewhere [1] .

Further search for enhanced UV filter molecules based on the drometrizole moiety led to the absorption spectra of the designed drometrizole derivatives studied and the results are presented in Table 8.

From this table, two intense absorption bands with the maximum band centered around 361, and the shoulder band around 301 nm in both the gas and water solvent are visible. These bands like those for PBT correspond to the S0→S1 and S0→S2 transitions respectively. These transitions correspond to the H→L and H-1→L for the maximum absorption and the shoulder peak respectively in both environments studied. It is also observed from Table 8 that all the absorption bands are red shifted as compared to that of PBT by factors of 33 and 21 nm for the maximum absorption and the shoulder peaks respectively in both environments under study. This bathochromic shift results from the presence of the EDS (-NH2) at the ortho-position of the -OH group on the drometrizole moiety i.e. corresponding to PBT1. No significant change is observed for the molecules obtained from substitutions at the ortho-position with the presence of the EAS (-CN) i.e. the PBT2 molecule. It is important to point out that the derivatives of PBT have absorptions in the region 288 - 395 nm, implying that they can readily absorb the hazardous UV A and UV B radiations from sunlight.

On the basis of the ongoing, the studied compounds especially PBT1 are potential candidates for the manufacture of UV filters and sunscreen lotions to fight the devastating effects of UV A and UVB radiations.

4. Conclusion

The antioxidant and UV spectral characteristics of Drometrizole (PBT) and the effect of an electron-donating group (ED), -NH2 and an electron-attracting group (EA), -CN groups on these properties have been analyzed theoretically in the gas and aqueous phases. From a thermodynamic perspective, the HAT mechanism is preferred in the gas phase, while SPLET in the aqueous phase. Also, it is observed that the presence of -NH2 reduces the O-HBDE and IP values, while that of -CN increases the PA values. Consequently, the molecule with the ED substituent (PBT1) is identified to be the compound with the highest antioxidant activity among the studied molecules. The spin density values on the O-radical are the lowest for PBT1 confirming that it is the best antioxidant molecule. The UV spectra of these molecules exhibit two intense absorption bands within the 280 - 400 nm region. This corresponds to the UV A and UV B regions of the electromagnetic spectrum. The major transitions in the UV-Visible spectra of the studied compounds are between the HOMO, HOMO-1 and HOMO-2 with LUMO level. The presence of the ED substituent is found to cause a bathochromic shift in absorption wavelengths. The studied molecules especially PBT1 can thus be considered as a potent antioxidant and UV absorption agent. The antioxidant capability of these molecules is therefore an added advantage for photoprotection as they can scavenge the free radicals. Results from this study can further be used for the design and synthesis of compounds with more improved photoprotective properties for applications in sunscreen products. Also, this study gives way for a subsequent kinetic study to make perfect the results deduced from this thermodynamic study.

Acknowledgements

We gratefully acknowledge the annual research modernization grant to lecturers by the Ministry of Higher Education of Cameroon.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Tasheh, N.S., Fouegue, A.D.T. and Ghogomu, J.N. (2021) Investigation of the Antioxidant and UV Absorption Properties of 2-(2’-hydroxy-5’-methylphenyl)-benzotriazole and Its Ortho-Substituted Derivatives via DFT/TDDFT. Computational Chemistry, 9, 161-196. https://doi.org/10.4236/cc.2021.93010

Supplementary Information Selected gas-phase computed bond lengths and angles of PBT <xref ref-type="table" rid="table">Table </xref>S1
Bond typeBond lengthsa (Å)Bond typeBond anglesa (˚)
Exptb.Comptc.Exptb.Comptc.
C1-C21.4071.401C2-C3-C4121.1120.9
C2-C31.4021.421C6-C1-C2116.1116.9
C3-C41.4041.401C4-C5-C6121.8121.9
C4-C51.3611.374C5-C6-C1122.8122.0
C5-C61.4091.414C5-C4-C3116.8117.1
C6-C11.3581.374C2-C3-N13109.2108.7
C2-N111.3491.348N11-C2-C1130.3131.0
C3-N131.3521.350N13-C3-C2109.2108.7
N11-N121.3401.341N12-N13-C3102.8103.5
N12-N131.3301.326N12-N11-C2103.2103.8
N12-C141.4301.415N12-C14-C15120.5120.6
C14-C151.3981.408C14-C15-O16124.6124.6
C14-C201.3851.394C15-C17-C18121.0121.3
C15-O161.3611.343C17-C18-C19121.7121.5
C15-C171.3831.396C18-C19-C20116.9117.8
C17-C181.3791.377C19-C20-C14122.3121.4
C18-C191.3921.398C20-C14-C15120.1120.6
C19-C201.3801.382C18-C19-C25121.2121.1
C19-C251.4991.498

aRefer to Figure 1 for atom numbering. bExperimental values of selected bond lengths and angles. cComputed values of selected bond lengths and angles from this work using the B97-3c method.

Enthalpies (in kJ/mol) of the electron H ( e − ) , proton H ( H + ) and hydrogen atom H ( H • ) used in this study at the M06-2x/6-311++G(d,p) level of theory
MediumElectronProtonH-atom
Gas3.1466.197−1301.81
Water−77.5−1055.7−1296.26
 Enthalpies (in kJ/mol) of the neutral molecules H ( BTOH ) used in this studied as computed at the M06-2x/6-311++G(d,p) level of theory
SpeciesGasWater
Phenol−806,850.3−806,877.0
PBT−1,945,828.3−1,945,852.8
PBT1−2,091,120.0−2,091,155.6
PBT2−2,187,994.2−2,188,027.5
 Enthalpies (in kJ/mol) of the radical species H ( BTO • ) of the molecules under investigation used in this studied as computed at the M06-2x/6-311++G(d,p) level of theory
SpeciesGasWater
Phenol−805,176.3−805,204.1
PBT−1,944,128.5−1,944,179.6
PBT1−2,089,462.1−2,089,528.2
PBT2−2,186,288.4−2,186,345.7
 Enthalpies (in kJ/mol) of the cationic radical species H ( BTOH • + ) of the molecules under investigation used in this studied as computed at the M06-2x/6-311++G(d,p) level of theory
SpeciesGasWater
Phenol−806,027.6−806,278.5
PBT−1,945,071.9−1,945,265.8
PBT1−2,090,429.5−2,090,634.8
PBT2−2,187,192.3−2,187,410.6
 Enthalpies (in kJ/mol) of the anionic species H ( BTO − ) of the molecules under investigation used in this studied as computed at the M06-2x/6-311++G(d,p) level of theory
SpeciesGasWater
Phenol−805,403.3−805,669.0
PBT−1,944,395.3−1,944,647.2
PBT1−2,089,693.6−2,089,947.0
PBT2−2,186,621.1−2,186,853.3
 Cartesian coordinates of the gas phase optimised geometry of neutral form of PBT
28
C−3.055915419345231.956478806747360.18111160950319
C−1.845013032340371.313047284038050.46665934216414
C−1.78913919457020−0.089635882480160.68527513644604
C−2.94828738478569−0.874721606474010.62179918183079
C−4.12731004505138−0.228187677460710.34070160062671
C−4.180167463159081.168413947893890.12345326059736
H−3.098074486637613.023320166715710.01482466629639
H−2.90635820495959−1.941484797128580.78796218738119
H−5.04477250661527−0.797526487066950.28258553664493
H−5.136310546966811.62412877846508−0.09426844268632
N−0.592788295702191.797687119441280.58974806961426
N0.144275380929350.712523304404370.86739229309917
N−0.50941711403224−0.438869207765290.93735002645170
C1.541146166962570.806314641378561.07083382375431
C2.195934166830322.050705792378990.99212061649174
O1.585661472875283.217593496105300.72690596178347
C3.574434985834232.066823630462051.21082270436012
C4.267644961023370.910339648674311.49167890947933
C3.62160988441622−0.327669380862381.56717894898357
C2.25638645674016−0.355847914561771.35536644338310
H1.71147576017589−1.285942590375371.40761844181401
H0.622944405023943.026876241616990.60262577312658
H4.076378799696053.022102171457851.15650855493710
H5.336167503130680.964770407132311.65874099073661
C4.38847573547735−1.584264774609631.84203787124110
H3.73611858790326−2.382094035425082.18957891152585
H4.89383007695634−1.944524989022620.94426343652391
H5.15548934019067−1.425366093679572.59864431388954
 Cartesian coordinates of the gas phase optimised geometry of neutral form of PBT1
30
C−2.947662820287052.026480184584780.19831926100417
C−1.755704047826771.351265289656630.48957371737747
C−1.73582664485798−0.054189319763000.69348551118245
C−2.91255193731574−0.810442275778290.60919763034913
C−4.07306415983764−0.133091982526650.32257593449206
C−4.090059035314591.266263399736780.12003823591484
H−2.963122727598413.095879988004200.04376904841125
H−2.89791082508299−1.879582585804230.76450072174093
H−5.00321637420174−0.679710322079860.24871560727256
H−5.032685032136241.74735267889062−0.10248941910324
N−0.493161284815961.803476861929320.63144907665784
N0.214745608077690.697489812259580.90593014101485
N−0.46769574556743−0.437306004305700.95597642536766
C1.611743225295500.752238556675621.12415431770585
C2.279871481206551.980059438121331.07676639455490
O1.697542284855643.172375102740620.82273247224273
C3.673783853647031.997012761352961.29071270103694
C4.334951865228860.805424358910411.54405078082654
C3.66281435627999−0.417567697633191.59076028273350
C2.29852773503055−0.436351275217151.38411958796984
H1.74252933028501−1.359406488251391.41386340336209
H0.733516123521662.999072950027900.67423209685579
H5.405277511239710.831534431640431.71482578157341
C4.41876543030049−1.682544373122431.85984566803855
H3.75073325899644−2.538166111507331.91910935004764
H5.14745346845345−1.885588841030781.07387254242035
H4.97002026723605−1.625256543719522.79886465453390
N4.315989566154263.224843433946641.30680257391063
H3.817010214244003.952695977303200.82231835432913
H5.288835054789663.194698594958501.05324714617623
 Cartesian coordinates of the gas phase optimised geometry of neutral form of PBT2
29
C−3.019105053819411.975274377120610.18610123592083
C−1.814437332021071.318890469284750.47122741528945
C−1.77097728944808−0.086128397155770.68292320302841
C−2.93685273866720−0.861613579753850.61285039689658
C−4.10809811935556−0.202579995590470.33261309408181
C−4.148394768721091.196525987131080.12226453232015
H−3.051704688824043.043103137114620.02526976512964
H−2.90476228299184−1.929385579023070.77358999205479
H−5.03074192786794−0.762768587474010.26935496706957
H−5.100385286048721.66081466756971−0.09508920874691
N−0.558984742536121.791331959237260.59976089518870
N0.167151918790580.697935793025290.87342415183469
N−0.49560828593910−0.448111824108810.93613535957497
C1.564072766138510.777725416723411.08083416376983
C2.221855154912962.020500476152001.00648120095197
O1.621483731644743.181269308962120.74524940927613
C3.612802934155422.025479345379981.23099868028155
C4.295169877005220.844461552605701.51200127190174
C3.63574033631295−0.376660773426041.57864909836813
C2.26676219533413−0.387917778090491.36249715267377
H1.71286058390979−1.312976735502481.41017708417538
H0.654721945907642.994813781065870.61832987217556
H5.361806241732230.896034828846501.68078291231911
C4.38458955031283−1.645164203486891.84882014509548
H3.73465507926669−2.409453810826702.26821837004751
H4.81464566041381−2.050325556001700.93168653846849
H5.20501756884774−1.482080704379082.54472423137782
C4.320158096868973.251277124786961.17561288125592
N4.933298874685934.235129299813421.13963118821889
 Cartesian coordinates of the gas phase optimised geometry of PBT radical
27
C−2.953332691739091.74898350608147−0.34206885994238
C−1.790145736541811.223074040008580.24553283467541
C−1.826210565257760.001657090546880.98412343728246
C−3.02950912413525−0.703451077519261.15397086215762
C−4.15057721561081−0.169964388895390.57293491226653
C−4.113564160514721.04175141563141−0.16630868032711
H−2.923261047975852.67379654314569−0.89976873538339
H−3.06095674196065−1.624581940890721.71810182724924
H−5.09763108551453−0.681859529891140.67765798178955
H−5.033425138221221.41177036813396−0.59822383517486
N−0.531071853544481.687993269920280.24140166235750
N0.124659831747350.769937285982800.95710521998348
N−0.58666809997206−0.269732753154481.42315824197017
C1.508198430797020.839845762022001.16292253195775
C2.154925856391302.140662085333341.35426434139206
O1.531202609738493.209791517763041.47181925731569
C3.596552593289842.091618175657191.47612331457752
C4.280551719250140.916437558707211.52420976003788
C3.61175432312103−0.324253160130191.41593444357443
C2.22916723708456−0.330823450834921.23382807301429
H1.70262465385125−1.268390053098811.13066946435931
H4.094312960971953.045476606317431.57598267025652
H5.356898835669080.920425095359601.64735017009994
C4.37961371139903−1.598428603096021.49185054700550
H3.76819986657590−2.458006924103451.23219849297199
H5.24560870724612−1.579057429480250.83002624832809
H4.76356212385515−1.756561009516272.50209381620474
 Cartesian coordinates of the gas phase optimised geometry of PBT1 radical
29
C−3.054584701601531.927658589769040.94835684820104
C−1.796504978180801.313789446338860.84041154782556
C−1.66850871307665−0.016516349818290.34677826760822
C−2.79744988018331−0.74936454273257−0.05297213979047
C−4.01640970355543−0.129933705462300.06068839499912
C−4.143928395841261.192325859865560.55574569923467
H−3.149207439974022.937109052492931.32224768465664
H−2.70198421282288−1.75666930734280−0.43297524163904
H−4.91150355858208−0.66008442844466−0.23587749050070
H−5.131887147463711.627855652786730.62093165104266
N−0.571070405278421.778882602521471.14574955263696
N0.219593826725330.755756594675380.82794037902012
N−0.36466616834294−0.350887394545160.35134400509916
C1.612055312449670.796836731286291.04633228363515
C2.346028878700611.995700428374290.70283628167634
O1.870706269077193.006620096276590.15481295940671
C3.789259898290921.942944011178441.00102312203817
C4.387479672134310.788824585670691.49727066400933
C3.63499418901651−0.343813306802751.74989852254125
C2.24244490920756−0.320154423972331.52101736279456
H1.65186257641958−1.201379348839481.72622944430792
H5.454312251657210.779006790217301.68795364141538
C4.26414321623756−1.590353560479912.27709927415883
H4.10913984634618−2.426026571182131.59323165358796
H5.33412663956379−1.471173405596782.42571975945432
H3.82167253567494−1.880844458613333.23084115423379
N4.451606838136093.081816076833390.74209202712655
H3.923275601573593.806636012794520.27929050478964
H5.452932843691923.125318272750960.74308218643011
 Cartesian coordinates of the gas phase optimised geometry of PBT2 radical
28
C−3.098622446731291.926519029819700.78683406252558
C−1.841051421370411.299614418186550.74061444031888
C−1.71801447648128−0.082715848124480.39706835028972
C−2.85155152159610−0.851986004601880.08328273883643
C−4.06480645206162−0.217840115246520.13376633272925
C−4.187552134214501.154980424770000.48111070287808
H−3.188047296494832.971316266454341.04612421917260
H−2.76133527823090−1.89507616706838−0.18349680492655
H−4.96364921461538−0.77266227569966−0.09826898951279
H−5.174672143290471.595996723848400.49937853876246
N−0.619383509982401.794021329527600.98575151175934
N0.172950796064210.735933801967600.78290381250769
N−0.42064118325891−0.420555122890170.43886479105798
C1.554331687181150.793099983862730.98768100984780
C2.280014973385542.027768481546250.67235207197156
O1.773140521146302.998895663447580.10428867075031
C3.708252713737711.994705542236121.01363809321166
C4.315474843536350.847257147270631.46463858176893
C3.58531765522023−0.336228583063831.66341743566641
C2.21011457301288−0.335364165498781.41957308403111
H1.63408315901964−1.234095602195831.58408441844891
H5.378046940525200.855428226775481.66739057018141
C4.26519183102217−1.557351065538562.17953419901750
H3.67469876969890−2.452337640433632.00450824087381
H5.24331856588061−1.691121139827221.72058914486499
H4.43096389593128−1.478416538664833.25641654263802
C4.462746963375923.170048341848930.80986168182638
N5.112709189589944.120744887291820.66486254850241
 Cartesian coordinates of the gas phase optimised geometry of PBT cationic radical
28
C−3.033623405912111.986809775767120.18928645882103
C−1.833250484418121.323798385977490.47217802224838
C−1.79731419063745−0.099317517701720.69356307271014
C−2.97183772918257−0.873060923075120.63259635521362
C−4.12940856369712−0.204476716464770.35527016385007
C−4.159435930218181.206896734873650.13656142531903
H−3.066350891819393.052781425192100.02257816158186
H−2.94452809442709−1.939167408816960.79850338160791
H−5.05976137969944−0.750349792680660.29674201766132
H−5.111928955014891.66963947722362−0.07836038522363
N−0.580191204277971.782575740198100.59282127629760
N0.143154460311440.677753592365560.87022563793939
N−0.53767071453062−0.471497999143300.94129632838533
C1.514554041438970.745866703286821.06801267965097
C2.185294955832672.023492721886780.97618821951162
O1.561732531992313.139777392900710.71209847819793
C3.580546184556962.067676245867831.18152765490840
C4.267225552244100.926711200400441.45880231685109
C3.61180014808693−0.336987581424961.55069435920932
C2.24014478161936−0.398207005521461.35156158050373
H1.71768330370075−1.339593907015941.41486496061962
H0.574362951492242.931480366041820.59342675455312
H4.066334618809443.029016621979341.10963216147672
H5.336288560160990.970669277586241.61459240970442
C4.40277768254067−1.551371167950201.85478915013890
H3.79100351532538−2.445819077756021.89890663332839
H5.18162703386560−1.699800584912001.10395203814963
H4.92519121185717−1.440305979084582.80720885678397
 Cartesian coordinates of the gas phase optimised geometry of PBT1 cationic radical
30
C−2.926316796579652.049548172603040.19567262237572
C−1.742369324933321.352728021768670.47690998385092
C−1.73904452018544−0.065047247027170.65948359625533
C−2.92477504210929−0.808721711653180.56293881340222
C−4.06994034721000−0.109171986901150.28765376508465
C−4.069613097635921.299149200312100.10686613983241
H−2.932842208715003.120323077339440.05714346746693
H−2.92311554732839−1.879412578041940.70070867081158
H−5.00815062139716−0.638422206382380.20368484077336
H−5.007863003558971.79020979479073−0.10813500798932
N−0.476711181727471.780058743782010.62716185149165
N0.216310749357430.651529175856980.88422743717103
N−0.47999133838461−0.471056077434740.91564593445130
C1.601993707514390.702120984299691.10231600510446
C2.255431784808961.955809679663381.04719776534548
O1.653126567917663.093208430799710.80347152035538
C3.682231268574621.997659904873521.27303790690566
C4.380774889979420.811462067284631.53885127325867
C3.72325974938998−0.394483681306991.58726664486861
C2.32079048743401−0.437552618286631.36538391018286
H1.79880058233715−1.381678693762411.40198137058725
H0.649665846219002.874393641676940.66908245626702
H5.447960140696990.857264292415161.70595010885163
C4.44055923823745−1.666549080315821.86729168294625
H4.03239125494105−2.151492179811242.75462589743208
H4.31176448523439−2.370282928631381.04428605650968
H5.50332432571517−1.511426674949742.02110926234979
N4.273905376227993.191970074716641.21541128089554
H3.736492132800334.021408573100861.02174248794421
H5.263400442379163.295413829221251.36235225521764
 Cartesian coordinates of the gas phase optimised geometry of PBT2 cationic radical
29
C−2.991217260494642.007350035660680.19741842545736
C−1.800657820108951.333944342842650.47866495210380
C−1.77766009587501−0.090498465032510.69179467330546
C−2.96093438617254−0.856111054999940.62339746645479
C−4.11156952895674−0.177241082567230.34753615594132
C−4.126672822675621.234702892213060.13745377208851
H−3.017013559972793.074382286560850.03645257388037
H−2.94100737124012−1.923269371266450.78346233802948
H−5.04682328110815−0.713764606346060.28342311127520
H−5.073666591249981.70879638582146−0.07716588508976
N−0.540157189912121.779282420902690.60607616417246
N0.169237129002110.670010150664650.87788952848084
N−0.52466146854005−0.476757902737090.94040722255974
C1.542088759141570.718427962322841.08101499499121
C2.216093729235721.989501181474670.99729791592094
O1.604239458645923.106454073099110.73927744195500
C3.627770636811062.019482225738291.20913930337583
C4.298956649872740.852381106714821.48359154059384
C3.62775972261292−0.392349340027771.56383421094396
C2.25360820688867−0.434337743791891.35952703499048
H1.72033293824629−1.370178358660191.41549571780753
H0.611098801091822.904976748504250.61676812292719
H5.367176490654070.889255442421651.64180689137461
C4.40045509193420−1.621291579798131.86289683204565
H3.77360394959684−2.505572956477481.89739171754940
H5.17931495768801−1.777124652386231.11400617253019
H4.91886644854781−1.525992677048762.81888283874525
C4.301423668601043.260437117763371.12979884645179
N4.866754737735874.270505418434611.06757991913748
 Cartesian coordinates of the gas phase optimised geometry of the anionic form of PBT
27
C−2.964575528095451.71813826549758−0.36244382439448
C−1.793941563407371.221449161629050.22349266513345
C−1.814285208183390.024092258847340.99236657570215
C−3.00803303809804−0.683961329417951.18261212666825
C−4.15191988312790−0.178264769954640.59902543640672
C−4.130389250220391.00622390541892−0.16673447307182
H−2.945671490890402.62991426811387−0.94458448271000
H−3.02627687477200−1.592655047197251.77043023002555
H−5.09209015744976−0.700598971454020.72971441352790
H−5.053840568618331.36300195799587−0.60654810936375
N−0.529133168402481.689639310063640.18869244145431
N0.145190775423530.805987062041380.92477996500788
N−0.56592624388758−0.216764543111311.43957393617201
C1.533983146174220.899555497360991.13382953361176
C2.162407749298262.197434153478631.34142734338873
O1.566292774087943.298252388701011.37990597019757
C3.589332585833452.086916345905061.55341539586775
C4.263490099804390.896637443121431.59679146715818
C3.60772313254958−0.337862760309291.42624487731429
C2.24563151251619−0.297706995655601.20518260016767
H1.69386782821873−1.219335517027121.06872192170718
H4.110333633564403.025104013523431.70217917851387
H5.338504453427670.901059006496361.76562146999595
C4.36415889793349−1.630190565059191.46077774066048
H3.69026723793457−2.481145234584751.36304604689337
H5.09954972268367−1.706346144735850.65295620842336
H4.91682942570294−1.760463159687622.39641334554160
 Cartesian coordinates of the gas phase optimised geometry of the anionic form of PBT1
29
C−3.060996523196681.946974562267060.79520863954266
C−1.799992415836151.339069079319830.77610890995638
C−1.65238034823242−0.022317832581940.38913728545066
C−2.76644682097183−0.784827808864150.01597078527932
C−4.00166004612240−0.168124751138480.03919144739223
C−4.147956040207031.179701482473730.42683960600691
H−3.169138874663132.982823952977261.08846734917532
H−2.65574344573791−1.81948477307294−0.28218184684678
H−4.88278891956323−0.73011260036414−0.24661508633144
H−5.137845231505271.619807116648060.43134010530328
N−0.584894615980301.830911753913811.09482040881934
N0.231392796672670.797519707919420.87815231217537
N−0.34526677377062−0.344525552987300.45251959578446
C1.616088010314700.858182518525621.11755218027196
C2.339168941558562.074992753254200.87880578510921
O1.878049121109773.159120955707270.42709890195503
C3.772563560276991.968688493886451.14796691427479
C4.359671599631490.796062620195621.54796890068630
C3.60690598854439−0.382989065380451.73321724439683
C2.25160784750287−0.330623018439791.51217729961417
H1.64062646886493−1.212133299100891.64822978399330
H5.431235819828880.776251098700781.73713319466993
C4.27731986896745−1.642216771246382.18874694259525
H3.56377448054820−2.462870166834732.25686831580545
H5.07186701863013−1.956283641953141.50429655737327
H4.74236307929998−1.531897199784423.17378838771974
N4.471855258408613.167594268234350.97481165767601
H3.854710346232033.799368987887290.46937712540930
H5.353839849395313.051217129837980.50010129674171
 Cartesian coordinates of the gas phase optimised geometry of the anionic form of PBT2
28
C−2.906780257078201.63003537126730−0.50431023721673
C−1.754692965444861.168763942185120.14619351197600
C−1.814835003982410.070022046106071.04779845303852
C−3.03034810545064−0.577303975307621.30955168052773
C−4.15182679841722−0.109182972612470.66144115461421
C−4.090862923729150.97999538138691−0.23668092349078
H−2.857562192025972.46792523954964−1.18642902737294
H−3.07901090611893−1.409982103635211.99876368528368
H−5.10836628928912−0.583742413834480.84253255860593
H−5.001725799094231.30969671800949−0.72066904638046
N−0.476590193508371.592863874088900.07107855186471
N0.161327918133300.779910235506330.90739478004089
N−0.57548400745566−0.155086938580061.52836438963665
C1.557638089506470.860432516141331.11379820532550
C2.168669541913492.167678149759161.33748720926482
O1.544665350688753.239740893896041.39430028816647
C3.613247669372822.064025569044081.53709708380138
C4.286806307606690.844269840394561.55843935772025
C3.63567446232966−0.368035695139491.38306407819448
C2.25913781441391−0.325231388249201.16593711835287
H1.70772819104077−1.247030835794861.02407189713651
H5.359165803332310.853834808579151.72186576638486
C4.37703117374256−1.670671940237031.38331974465566
H3.89297492236622−2.414550207709432.01959833626948
H4.44882057673692−2.109489373047000.38361578530193
H5.39518446339706−1.541104228732881.74944017300732
C4.349584271870743.247271000652071.74644898278875
N5.002458885143064.195526486313511.92325644250224
 Cartesian coordinates of the gas phase optimised geometry of neutral form of Phenol
13
C−2.452374515935951.93997528076273−0.02582868805402
C−1.371435018210141.066182848070960.03876897963578
C−1.59124536634413−0.302724712907420.07701304743363
C−2.88169589721732−0.810262426266010.05146735356627
C−3.955309195327940.06904195554530−0.01305603553135
C−3.749494374406241.43757943509676−0.05172546405231
O−2.302310361484163.29787797496739−0.06610047861539
H−1.362075866057823.50951261975985−0.04184979507995
H−0.359612962802971.455326242027750.05904339577356
H−0.74404804430274−0.973912719958980.12715709031594
H−3.04953988813392−1.877831376018830.08140296612123
H−4.96736389853680−0.31327430682404−0.03357889547072
H−4.578904611239842.12888918574453−0.10192347604266
 Cartesian coordinates of the gas phase optimised geometry of Phenolic radical
12
C−2.735018208789671.95284375076834−0.05123045197369
C−1.540306299051451.140556081690420.01552096708406
C−1.62110259626482−0.224832181889520.07551306878256
C−2.87019080512006−0.860137411097890.07304761718924
C−4.05284862911457−0.111012242504750.00941817441481
C−4.002257065484461.25574381809561−0.05136676246012
O−2.675002169552963.20504937560642−0.10648800097934
H−0.591610755386561.658624231692070.01579403940128
H−0.72056359422903−0.823043350334780.12530581315976
H−2.92189493024017−1.939294788347370.12075953867259
H−5.00663486777962−0.622470122977210.00873272848750
H−4.896670078986621.86070283929867−0.10123673177866
 Cartesian coordinates of the gas phase optimised geometry of Phenolic cationic radical
13
C−2.443151369066741.93137544502077−0.02517593041211
C−1.327918606592191.049756535532560.04079854115914
C−1.55768458128042−0.293794580968290.07780587073701
C−2.88268984066563−0.787892279666230.05067371434882
C−3.984691686534520.08677297030023−0.01462523792434
C−3.780938382094091.43524177822578−0.05265117099545
O−2.316851486309233.23777666054159−0.06446838214934
H−1.388050820036163.52763593674504−0.04355576590399
H−0.322704457690151.449895534390180.06054004704033
H−0.73088387265492−0.986808579288780.12797810929560
H−3.04961371092455−1.855161237668350.08064186902499
H−4.98721557005567−0.31397890481796−0.03425601034227
H−4.593015616095692.14556072165345−0.10291565387841
 Cartesian coordinates of the gas phase optimised geometry of the Phenolic anion
12
C−2.445269091426472.02923402858075−0.02855701739006
C−1.374554115323441.065956953025610.03863837711111
C−1.59651996766629−0.296930826753640.07662704168464
C−2.88671236088660−0.831064025782790.05203232711885
C−3.956501474079490.06466993803139−0.01295241904460
C−3.757129373174941.43083086429520−0.05177867455199
O−2.252361256464543.28281722307097−0.06386231270575
H−0.362298996402741.457916733135000.05884008078921
H−0.74280522296212−0.970806938709630.12708522321310
H−3.05204749824726−1.902020680075190.08215199417783
H−4.97379653113973−0.32284738129657−0.03346020369565
H−4.603334112226382.10913411247891−0.10213441670671
ReferencesTasheh, N.S., Nkungli, N.K. and Ghogomu, J.N. (2019) A DFT and TD-DFT Study of ESIPT-Mediated NLO Switching and UV Absorption by 2-(2’-hydroxy-5’-me- thylphenyl) Benzotriazole. Theoretical Chemistry Accounts, 138, Article No. 100. https://doi.org/10.1007/s00214-019-2488-0Dunaway, S., Odin, R., Zhou, L., Ji, L., Zhang, Y. and Kadekaro, A.L. (2018) Natural Antioxidants: Multiple Mechanisms to Protect Skin from Solar Radiation. Frontiers in pharmacology, 9, Article No. 392. https://doi.org/10.3389/fphar.2018.00392Chinh, N.T., Le Anh, N.T., Thao, P.T. and Quang, D.D. (2020) Photoprotective Properties of Natural Antioxidant Flavonoids: A DFT and TD-DFT Study on Acridone Derivatives. Vietnam Journal of Chemistry, 58, 157-161. https://doi.org/10.1002/vjch.201900083Dao, D.Q., Phan, T. T. T., Nguyen, T.L.A., Trinh, P.T. H., Tran, T.T.V., Lee, J.S., et al. (2020) Insight into Antioxidant and Photoprotective Properties of Natural Compounds from Marine Fungus. Journal of Chemical Information and Modeling, 60, 1329-1351. https://doi.org/10.1021/acs.jcim.9b00964Abiola, T.T., Rodrigues, N.D.N., Ho, C., Coxon, D.J., Horbury, M.D., Toldo, J.M., et al. (2020) New Generation UV-A Filters: Understanding their Photodynamics on a Human Skin Mimic. The Journal of Physical Chemistry Letters, 12, 337-344. https://doi.org/10.1021/acs.jpclett.0c03004Shireen, P.A., Mujeeb, V.A. and Muraleedharan, K. (2017) Theoretical Insights on Flavanones as Antioxidants and UV Filters: A TDDFT and NLMO Study. Journal of Photochemistry and Photobiology B: Biology, 170, 286-294. https://doi.org/10.1016/j.jphotobiol.2017.04.021Galanakis, C.M., Tsatalas, P. and Galanakis, I.M. (2018) Phenols from Olive Mill Wastewater and Other Natural Antioxidants as UV Filters in Sunscreens. Environmental Technology & Innovation, 9, 160-168. https://doi.org/10.1016/j.eti.2017.12.002Lorigo, M. and Cairrao, E. (2019) Antioxidants as Stabilizers of UV Filters: An Example for the UV-B Filter Octylmethoxycinnamate. Biomedical Dermatology, 3, Article No. 11. https://doi.org/10.1186/s41702-019-0048-9Rajan, V.K., Ahamed, T.S. and Muraleedharan, K. (2018) Studies on the UV Filtering and Radical Scavenging Capacity of the Bitter Masking Flavanone Eriodictyol. Journal of Photochemistry and Photobiology B: Biology, 185, 254-261. https://doi.org/10.1016/j.jphotobiol.2018.06.017Liu, T. and Wu, D. (2011) Simultaneous Determination of Some Ultraviolet-Ab- sorbing Chemicals in Sunscreen Cosmetics Using a High-Performance Liquid Chromatography Method. International Journal of Cosmetic Science, 33, 408-415. https://doi.org/10.1111/j.1468-2494.2011.00643.xEl-Asmy, H.A., Butler, I.S., Mouhri, Z.S., Jean-Claude, B.J., Emmam, M.S. and Mostafa, S.I. (2014) Zinc (II), Ruthenium (II), Rhodium (III), Palladium (II), Silver (I), Platinum (II) and MoO2+ 2 Complexes of 2-(2’-hydroxy-5’-methylphenyl)-ben- zotriazole as Simple or Primary Ligand and 2, 2’-bipyridyl, 9, 10-phenanthroline or Triphenylphosphine as Secondary Ligands: Structure and Anticancer Activity. Journal of Molecular Structure, 1059, 193-201. https://doi.org/10.1016/j.molstruc.2013.11.039Lee, J.K., Kim, K.B., Lee, J.D., Shin, C.Y., Kwack, S.J., Lee, B.M. and Lee, J.Y. (2019) Risk Assessment of Drometrizole, a Cosmetic Ingredient Used as an Ultraviolet Light Absorber. Toxicological Research, 35, 119-129. https://doi.org/10.5487/TR.2019.35.2.119Lee, C.Y., Sharma, A., Semenya, J., Anamoah, C., Chapman, K.N. and Barone, V. (2020) Computational Study of Ortho-Substituent Effects on Antioxidant Activities of Phenolic Dendritic Antioxidants. Antioxidants, 9, Article No. 189. https://doi.org/10.3390/antiox9030189Javan, A.J., Javan, M.J. and Tehrani, Z.A. (2013) Theoretical Investigation on Antioxidant Activity of Bromophenols from the Marine Red Alga Rhodomela confervoides: H-Atom vs Electron Transfer Mechanism. Journal of agricultural and food chemistry, 61, 1534-1541. https://doi.org/10.1021/jf304926mGalano, A., Mazzone, G., Alvarez-Diduk, R., Marino, T., Alvarez-Idaboy, J.R. and Russo, N. (2016) Food Antioxidants: Chemical Insights at the Molecular Level. Annual Review of Food Science and Technology, 7, 335-352. https://doi.org/10.1146/annurev-food-041715-033206Vo, Q.V., Nam, P.C., Van Bay, M., Thong, N.M., Cuong, N.D. and Mechler, A. (2018) Density Functional Theory Study of the Role of Benzylic Hydrogen Atoms in the Antioxidant Properties of Lignans. Scientific Reports, 8, Article No. 12361. https://doi.org/10.1038/s41598-018-30860-5Wang, G., Liu, Y., Zhang, L., An, L., Chen, R., Liu, Y., et al. (2020) Computational Study on the Antioxidant Property of Coumarin-Fused Coumarins. Food Chemistry, 304, Article ID: 125446. https://doi.org/10.1016/j.foodchem.2019.125446Neese, F. (2012) The ORCA Program System. WIREs Computational Molecular Science, 2, 73-78. https://doi.org/10.1002/wcms.81Neese, F. (2017) Software Update: The ORCA Program System, Version 4.0. WIREs Computational Molecular Science, 8, e1327. https://doi.org/10.1002/wcms.1327Hanwell, M.D., Curtis, D.E., Lonie, D.C., Vandermeersch, T., Zurek, E. and Hutchison, G.R. (2012) Avogadro: An Advanced Semantic Chemical Editor, Visualization, and Analysis Platform. Journal of Cheminformatics, 4, Article No. 17. https://doi.org/10.1186/1758-2946-4-17Brandenburg, J.G., Bannwarth, C., Hansen, A. and Grimme, S. (2018) B97-3c: A Revised Low-Cost Variant of the B97-D Density Functional Method. The Journal of Chemical Physics, 148, Article ID: 064104. https://doi.org/10.1063/1.5012601Zhao, Y. and Truhlar, D.G. (2008) The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. Theoretical Chemistry accounts, 120, 215-241. https://doi.org/10.1007/s00214-007-0310-xNeese, F., (2003) An Improvement of the Resolution of the Identity Approximation for the Formation of the Coulomb Matrix. Journal of Computational Chemistry, 24, 1740-1747. https://doi.org/10.1002/jcc.10318Neese, F., Wennmohs, F., Hansen, A. and Becker, U. (2009) Efficient, Approximate and Parallel Hartree-Fock and hybrid DFT Calculations. A ‘Chain-of-Spheres’ Algorithm for the Hartree-Fock Exchange. Chemical Physics, 356, 98-109. https://doi.org/10.1016/j.chemphys.2008.10.036Goerigk, L. and Grimme, S. (2011) A Thorough Benchmark of Density Functional Methods for General Main Group Thermochemistry, Kinetics, and Noncovalent Interactions. Physical Chemistry Chemical Physics, 13, 6670-6688. https://doi.org/10.1039/C0CP02984JMarenich, A.V., Cramer, C.J. and Truhlar, D.G. (2009) Universal Solvation Model Based Onsolute Electron Density and on a Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface Tensions. The Journal of Physical Chemistry B, 113, 6378-6396. https://doi.org/10.1021/jp810292nBine, F.K., Tasheh, S.N. and Nkungli, N.K. (2020) Corrosion Inhibition of Aluminium in Gas and Acid Media by Some Chalcone-Based N-(3-Aminopropyl) Imidazoles: TD-DFT-Based FMO, Conceptual DFT, QTAIM and EDA Studies. Computational Chemistry, 9, 37-63. https://doi.org/10.4236/cc.2021.91003Grimme, S., Brandenburg, J.G., Bannwarth, C. and Hansen, A. (2015) Consistent Structures and Interactions by Density Functional Theory with Small Atomic Orbital Basis Sets. Journal of Chemical Physics, 143, Article ID: 054107. https://doi.org/10.1063/1.4927476Weigend, F. and Ahlrichs, R. (2005) Balanced Basis Sets of Split Valence, Triple Zeta Valence and Quadruple Zeta Valence Quality for H to Rn: Design and Assessment of Accuracy. Physical Chemistry Chemical Physics, 7, 3297-3305. https://doi.org/10.1039/B508541ACaldeweyher, E., Bannwarth, C. and Grimme, S. (2017) Extension of the D3 Dispersion Coefficient Model. The Journal of Chemical Physics, 147, Article ID: 034112. https://doi.org/10.1063/1.4993215Ahamed, T.S., Rajan, V.K., Sabira, K. and Muraleedharan, K. (2019) DFT and QTAIM Based Investigation on the Structure and Antioxidant Behavior of Lichen Substances Atranorin, Evernic acid and Diffractaic Acid. Computational biology and chemistry, 80, 66-78. https://doi.org/10.1016/j.compbiolchem.2019.03.009Markovi&#263, Z., Markovi&#263, Z., To&#353ovi&#263, J., Milenkovi&#263, D. and Markovi&#263, S. (2016) Revisiting the Solvation Enthalpies and Free Energies of the Proton and Electron in Various Solvents. Computational and Theoretical Chemistry, 1077, 11-17. https://doi.org/10.1016/j.comptc.2015.09.007Bartmess, J.E. (1994) Thermodynamics of the Electron and the Proton. The Journal of Physical Chemistry, 98, 6420-6424. https://doi.org/10.1021/j100076a029Bizarro, Cabral, B.C., Dos Santos, R.B. and Sim&#245es, J.M. (1999) Substituent Effects on the O-H Bond Dissociation Enthalpies in Phenolic Compounds: Agreements and Controversies. Pure and Applied Chemistry, 71, 1249-1256. https://doi.org/10.1351/pac199971081609Boulebd, H. (2020) The Role of Benzylic-Allylic Hydrogen Atoms on the Antiradical Activity of Prenylated Natural Chalcones: A Thermodynamic and Kinetic Study. Journal of Biomolecular Structure and Dynamics, 39, 1955-1964. https://doi.org/10.1080/07391102.2020.1740791Catalán, J., Torres, M.R. and Tornero, J.D. (1997) Molecular Structure of a Unique UV Stabilizer: Tinuvin P. Journal of the Chemical Society, Faraday Transactions, 93, 1691-1696. https://doi.org/10.1039/A607652AZhurko, G. and Zhurko, D. (2013) Chemcraft (Version 1.7, Build 375, 2013).Jin, R. (2010) A DFT Study on the Radical Scavenging Activity of Juglone and Its Derivatives. Journal of Molecular Structure: THEOCHEM, 939, 9-13. https://doi.org/10.1016/j.theochem.2009.09.024Alisi, I.O., Uzairu, A. and Abechi, S.E. (2019) In Silico Design of Hydrazone Antioxidants and Analysis of Their Free Radical-Scavenging Mechanism by Thermodynamic Studies. Beni-Suef University Journal of Basic and Applied Sciences, 8, Article No. 11. https://doi.org/10.1186/s43088-019-0011-2&#214zbak&#305r I&#351&#305n, D. (2016) Theoretical Study on the Investigation of Antioxidant Properties of Some Hydroxyanthraquinones. Molecular Physics, 114, 3578-3588. https://doi.org/10.1080/00268976.2016.1248514Mikulski, D., Eder, K. and Molski, M. (2014) Quantum-Chemical Study on Relationship between Structure and Antioxidant Properties of Hepatoprotective Compounds Occurring in Cynara scolymus and Silybum marianum. Journal of Theoretical and Computational Chemistry, 13, Article ID: 1450004. https://doi.org/10.1142/S0219633614500047Farmanzadeh, D. and Najafi, M. (2013) On the Antioxidant Activity of the Tryptophan Derivatives. Bulletin of the Chemical Society of Japan, 86, 1041-1050. https://doi.org/10.1246/bcsj.20130035Fouegue, A.D.T., Mama, D.B., Ghogomu, J.N., Elie, Y. and Etoh, M.A. (2018) The Substitution Effect on Reaction Enthalpies of Antioxidant Mechanisms of Juglone and Its Derivatives in Gas and Solution Phase: DFT Study. Journal of Chemistry, 2018, Article ID: 1958047. https://doi.org/10.1155/2018/1958047Akhtari, K., Hassanzadeh, K., Fakhraei, B., Fakhraei, N., Hassanzadeh, H. and Zarei, S.A. (2013) A Density Functional Theory Study of the Reactivity Descriptors and Antioxidant Behavior of Crocin. Computational and Theoretical Chemistry, 1013, 123-129. https://doi.org/10.1016/j.comptc.2013.03.015Najafi, M., Najafi, M. and Najafi, H. (2012) DFT/B3LYP Study of the Substituent Effects on the Reaction Enthalpies of the Antioxidant Mechanisms of Magnolol Derivatives in the Gas-Phase and Water. Bulletin of the Korean Chemical Society, 33, 3607-3617. https://doi.org/10.5012/bkcs.2012.33.11.3607Mishra, A., Verma, C., Srivastava, V., Lgaz, H., Quraishi, M., Ebenso, E.E. and Chung, I.-M. (2018) Chemical, Electrochemical and Computational Studies of Newly Synthesized Novel and Environmental Friendly Heterocyclic Compounds as Corrosion Inhibitors for Mild Steel in Acidic Medium. Journal of Bio- and Tribo- Corrosion, 4, Article No. 32. https://doi.org/10.1007/s40735-018-0147-yObot, I., Kaya, S., Kaya, C. and B. Tüzün (2016) Density Functional Theory (DFT) Modeling and Monte Carlo Simulation Assessment of Inhibition Performance of Some Carbohydrazide Schiff Bases for Steel Corrosion. Physica E: Low-Dimensional Systems and Nanostructures, 80, 82-90. https://doi.org/10.1016/j.physe.2016.01.024