Inhibitive action of Chamomile extract on the

17 downloads 0 Views 660KB Size Report
Herein, Fe will act as a Lewis acid while the inhibitors act as a Lewis base. .... According to this principle, soft bases inhibitor is the most effective ones for metals ...
Moroccan Journal of Chemistry ISSN: 2351-812X http://revues.imist.ma/?journal=morjchem&page=login

Ghassab & al. / Mor. J. Chem. 6 N°1 (2018) 195-202

Inhibitive action of Chamomile extract on the corrosion of Iron: Density Functional Theory Ghassab M. Al-Mazaideh*1, Saleh A. Al-Quran2 1

Department of Chemistry and Chemical Technology, Faculty of Science, Tafila Technical University, P.O. Box 179, Tafila 66110, Jordan 2 Department of Biological Sciences, Faculty of Science, Mutah University, P.O. Box 7, Mutah 61710, AlKarak, Jordan

Abstract DFT calculations have been performed on major compounds of Chamomile extract as a green source of environmentally friendly corrosion inhibitors for Fe using B3LYP/631G* level (d, p). Several global Quantum parameters were calculated on inhibitors and thermodynamic Gibbs function ∆Gads of adsorption of Fe has been calculated and used * Corresponding author: to evaluate the inhibitive action of the each inhibitor. Quercetin has been found the [email protected]

highest anti-corrosion efficiency as compared to other compounds. Consequently,

Received 14 May 2017, calculated global quantum parameters and thermodynamic function ∆Gads show Revised

01

Sept

2017, spontaneous physical adsorption of Apigenin-7-glucoside, luteolin-7-glucoside and

Accepted 11 Dec 2017.

Quercetin on Fe, wherase other inhibitors show non- spontaneous chemical adsorption.

Keywords: Fe, DFT, Inhibitors, Chamomile

Mor. J. Chem. 6 N°1 (2018) 195-202

195

1. Introduction There are lots of reports of a retun to natural products such as herb plants in pharmaceutical chemistry in order to find drugs and develop potential compounds to find new drugs, antimicrobial, cytotoxic activities and anticancer [1]. Furthermore, the natural products is proved as new cleaning inhibitors for green environment. Several studies have been showed the use of natural products as friendly corrosion inhibitors for importanr metals like Fe, Al and Cu alloys in aggressive media [2-5]. These cleaning inhibitors are biodegradable, nontoxic and available in different plants. Recent research showed that the extracts of every part of plant such as seeds, fruits, flowers, leaves, stems and roots have used as corrosion inhibitors [6-9]. The extract solution of Chamomile (Chamaemelum mixtum L.) showed an excellent effective inhibitionfor the corrosion of steel in H2SO4. The main chemical compositions of Chamomile extracts are shown in Fig 1[10].

Figure 1: Structures of main chemical constituents of Chamomile To the best of our knowledge [11-17], DFT model has been used to compute a spread array of properties like optimized structure, thermochemistry and kinetic of different types of inhibitors. Such information will help in explain the effect of these compounds as inhibitors for Fe corrosion. The target of this paper is test the theoretical chemistry and molecular dynamics simulations of the inhibitors. In addition to discuss their mechanism of inhibitive action. These methods were based on quantum chemical calculations and have been proved to be a very suitable tools and cost effective for finding the corrosion inhibition mechanism. Mor. J. Chem. 6 N°1 (2018) 195-202

196

Density functional theory (DFT) model has classified as an important tool to investigate the efficiency of each compounds as a corrosion inhibitor of Fe. Predict the systematic way for the analysis of the interaction of inhibitor with Fe surface is very interesting and very aim of this work. Theoretical calculations have been used to predict the inhibitive activity of inhibitors and the efficiency of each one on the corrosion of Fe. In the present work, the chemical reactivity of the inhibitor has been reprted as a function of different quantum parameters such as EHOMO (highest occupied molecular orbital), ELUMO (the lowest unoccupied molecular orbital), ΔEgap (energy gap), η (global hardness), σ (softness), X (electronegativity), ΔN (the fraction of electron transferred), and ω (electrophilicity index).

2. Computational Methodology DFT model were used with B3LYP as hybrid functional with 6-31G* (d, p) basis set. Gaussian 09 (G09) [18] with DFT-B3LYP was used to calculate the quantum chemical parameters by optimizing the inhibitor structure under no constraint. Moreover, the results of these parameters were used to calculate thermodynamic function of adsorption ∆Gads (Gibbs free energy of adsorption) of the Fe and the studied inhibitors. In this work, the theoretical values of XFe = 4.026 eV and ηFe = 3.875 eV were used for calculating the number of electrons transferred [19]. All calculations have been performed in gas phase.

3.

Results and Discussion

The optimized structure of Inhibitors under investigation and together with their HOMO and LUMO orbitals are shown in Figure .2 which obtained by B3LYP/6-31G* level. The values of EHOMO and ELUMO for Fe metal [20] have been compared to the calculated values for inhibitors to determine the type of the interaction. The calculated results of EHOMO and ELUMO for inhibitors are shown in Table .1. The adsorption ability of the inhibitor over Fe surface can be calculated by frontier molecular orbital theory [21]. The capability of an inhibitor to donate electrons to a suitable acceptor of inhibitor with unoccupied d-orbital of Fe can be determined by EHOMO. While ELUMO shows its ability to accept electrons. Consequently, Better inhibitor efficiency means lower ELUMO value which leads to more ability of the molecule to accept electrons [22]. In this case, increasing EHONO and decreasing ELUMO level lead to increase in the binding capability of the inhibitor to the Fe surface. Herein, Fe will act as a Lewis acid while the inhibitors act as a Lewis base. So, inhibitor utilizes the HOMO to initiate the reaction with LUMO of Fe. The interaction has a certain amount of ionic character because the values of ELUMO, Fe EHOMO, inh gap approximately fall between 5 to 5.9 eV (Table 2). Strong covalent bond can be obtained only if the ELUMO, inh - EHOMO, Fe gap is approximately zero [23]. In this case, all inhibitors act as an anodic inhibitor. ∆Egap is a good parameter to predict the adsorption of inhibitor on metallic surface. i.e the reactivity of the inhibitors incrases with decrasing ∆Egap which leads to increase the inhibitor efficiency [24]. Consequently, Quercetin (lowest ∆Egap) has the strongest ability for co-ordinate bonds the vacant d-orbitals of Fe by donating and accepting electrons. The ∆Egap values of the other inhibitors (Table 3) increases in the order Quercetin < trans-spiroeth < luteolin-7glucoside < Apigenin 7-glucoside < cis-spiroeth < Bisabolol oxides = A Bisabolol oxides B. The effectiveness of Chamomile compounds under investigation as inhibitors has been further studied by evaluating the global reactivity parameters such as: η (the global chemical hardness), X

(the electronegativity), σ (the global

Inhibitor

softness), ΔN (the fraction of electrons transferred) and ω (the electrophililcity) are shown in Table 3 After Optimization

HOMO

LUMO

Mor. J. Chem. 6 N°1 (2018) 195-202

197

Apigenin 7-

glucoside

luteolin-7-

glucoside

cis-spiroeth trans-spiroeth Quercetin Bisabololoxides B Bisabololoxides A Figure 2: Optimized structure with HOMO and LUMO of studied inhibitors

Table 1. Calculated EHOMO and ELUMO of the inhibitors by the DFT Substance

EHOMO (eV)

ELUMO (eV)

Mor. J. Chem. 6 N°1 (2018) 195-202

198

Fe

-7.902

Apigenin 7-glucoside

-0.151

-5.920

-1.780

-5.870

-1.770

cis-spiroeth

-5.240

-0.640

trans-spiroeth

-5.180

-1.360

Quercetin

-5.560

-1.910

Bisabolol oxides A

-6.060

1.000

Bisabolol oxides B

-6.070

0.990

luteolin-7-glucoside

Table 2. Calculated energy gap interaction of Fe- inhibitor by the DFT method Inhibitors

ELUMO, inh – EHOMO, Fe Cu(eV)

ELUMO, Fe – EHOMO, inh (eV)

Apigenin 7-glucoside

6.122

5.769

6.132

5.719

cis-spiroeth

7.262

5.089

trans-spiroeth

6.542

5.029

Quercetin

5.992

5.409

Bisabolol oxides A

8.902

5.909

Bisabolol oxides B

8.892

5.919

luteolin-7-glucoside

Table 3. Calculated quantum chemical parameters for the inhibitors Apigenin Inhibitors

luteolin

7-

cis-

trans-

Quercetin

Bisabolol

Bisabolol

oxides A

oxides B

7-glucoside

glucoside

spiroeth

spiroeth

EHOMO

-5.920

-5.870

-5.240

-5.180

-5.560

-6.060

-6.070

ELUMO

-1.780

-1.770

-0.640

-1.360

-1.910

1.000

0.990

∆Egap I= - EHOMO

4.140

4.100

4.600

3.820

3.650

7.060

7.060

5.920

5.870

5.240

5.180

5.560

6.060

6.070

A= - ELUMO

1.780

1.770

0.640

1.360

1.910

-1.000

-0.990

χ = (I+A) / 2

3.850

3.820

2.940

3.270

3.735

2.530

2.540

η =(I-A) / 2

2.070

2.050

2.300

1.910

1.825

3.530

2.540

σ = 1/η

0.483

0.488

0.435

0.524

0.548

0.283

0.394

0.015

0.017

0.088

0.065

0.026

0.101

0.116

ω = μ2/2η

1.913

1.883

1.115

1.380

1.800

0.826

0.832

µ (Debye)

2.630

4.030

3.770

2.800

1.070

1.340

1.830

ΔN = (χ Fe - χ inh) / 2 (η Fe + η inh)

I , A , χ, η, σ, ω [25] / ΔN, μ (-X) [26,27]

Hard- Soft- Acid- Base principle says that Hard acids tend to co-ordinate hard bases while soft acids tend to coordinate soft bases. So that, hard molecules have high ∆Egap wherase small molecules have small ∆Egap [24]. According to this principle, soft bases inhibitor is the most effective ones for metals [28]. For this reason, Quercetin Mor. J. Chem. 6 N°1 (2018) 195-202

199

which has the small ∆Egap and the highest softness (σ) is expected to have the good inhibitor efficiency as compared to other compounds. This result is confirmed by calculating softness (σ) of the inhibitor and so its reactivity. Table 3 shows that Quercetin has the highest value compared to other inhibitors. In addition, it can be observed (Table 3) that the hardness (η) of all inhibitors are the largest values than Quercetin. This tendency is the revers of what has been obtained for softness (σ). Consequently, Quercetin with highest values of softness (σ) and smallest value of hardness (η) is the best inhibitor. Soft molecule is more reactive than a hard molecule [29]. The fraction of transferred electrons (ΔN) is also calculated and shows that most electrons transferred to the Fe surface comes from the Quercetin molecules when compared to other inhibitors. The electrophililcity index (ω) indicates the ability of the inhibitor to accept electrons from Fe. Quercetin exhibits the highest electrophililcity index (ω) value as compared to other inhibitor because of low ELUMO. As a result, Quercetin has high capacity to accept electrons from Fe. That is, Fe acts as Lewis base while all inhibitors act as Lewis acids (cathodic inhibitor). Fe will act as anodic inhibitor, because it can accept electrons from all inhibitor to form a co-ordinated bond. In addition, the inhibitor can accept electrons from Fe atom to form back-donating bond depending on the orientation of inhibitors’optimized structure on the spatial. In this case, the inhibition efficiency will be increase, because of donation and back-donation processes strengthen the adsorption of inhibitors onto the Fe surface. By assuming that Fe act as anode and inhibitors act as cathode, the ∆Gads values of the inhibitors on Fe surface are calculated by the following relations: ∆Gads (cathodic inhibitor) = μ (inhibitor) - μ (Fe) ∆Gads (anodic inhibitor) = μ (Fe) - μ (inhibitor) Where μ indicates the chemical potential and is equal to the negative value of the absolute electronegativity (X) as in ref. [30]. The calculated values of ∆Gads for inhibitors are given in Table 4. Table 4. Calculated ∆Gads (kJ /mol) and Dipole moment of the investigated inhibitors Inhibitors

Anodic Inhibitors

Apigenin-7-glucoside luteolin-7-glucoside cis-spiroeth trans-spiroeth Quercetin Bisabolol oxides A Bisabolol oxides B

-16.981 -19.876 -104.783 -72.943 -28.077 -144.342 -143.377

µ (Debye) 2.630 4.030 3.770 2.800 1.070 1.340 1.830

The adsorption of the inhibitors on Fe, shows that all inhibitors have negative ∆G ads values. The value of ∆Gads changes depends on the differences between chemical and physical adsorption [31,32]. So that, physical adsorption when ∆Gads value is in the range of 0 to - 40 kJ.mol-1, whereas chemical adsorption when ∆Gads value is in the range of -80 to -400 kJ.mol-1 [33]. The suggested mechanism for Apigenin-7-glucoside, luteolin-7-glucoside and Quercetin is spontaneous physical adsorption because ∆Gads value is in the range of 0 to -40 kJ.mol-1 wherase other inhibitors show non- spontaneous chemical adsorption.

4. Conclusion The theoretical calculations gave a good picture about the effect of main chemical constituents of Chamomile extract by suggesting that Quercetin shows the most inhibition efficiency as compared to the other compounds, because they have Mor. J. Chem. 6 N°1 (2018) 195-202

200

low ELUMO. Therefore, it is expected that inhibitors can form a strong interaction with Fe to act as cathodic inhibitors. All inhibitors show spontaneous energy and suggest that Chamomile extract have inhibitive action which is mainly from Quercetin.

References [1] A. T. Bull, J. E. M. Stach, Trends Microbiol, (2007); 15: 491. [2] A. M. Abdel-Gaber, E. Khamis, H. Abo-ElDahab, Sh. Adeel, Mater. Chem. Phys. ,(2008); 109: 297. [3] Y. Li, P. Zhao, Q. Liang, B. Hou, Appl. Surf. Sci., (2005); 252: 1245. [4] R. Kanojia, G. Singh, Surf. Eng., (2005); 21:180. [5] E. Khamis, N.M. Al Andis, Materialwiss. Werkst., (2002); 9: 550. [6] F. Zucchi, I.H. Omar, Surf. Technol., (1985); 24: 391. [7] I.H. Farooqi, M.A. Quraishi, P.A. Saini, Corros. Prev. Control, (1999); 46: 93. [8] M. Kliškic´ , J. Radoševic´ , S. Gudic´ , V. Katalinic´ , J. Appl. Electro. Chem., (2000); 30: 823. [9] A. Minhaj, P.A. Saini, M.A. Quraishi, I.H. Farooqi, Corros. Prev. Control, (1999): 32. [10] A.M. Abdel-Gaber, B.A. Abd-El-Nabey, I.M. Sidahmed, A.M. El-Zayady, M. Saadawy, Corrosion Science, (2006); 48: 2765. [11] G. M. Al-Mazaideh, K. A. Abu-Sbeih, S. M. Khalil, J. Chem., Biolo., Phy. Sci., (2017); 7(2): 398. [12] G. M. Al-Mazaideh, W. A. Al-Zereini, A. H. Al-Mustafa, S. M. Khalil, Adv. Env. Biol., (2016); 10 (8): 159. [13] A. M. Al-Msiedeen , G. M. Al-Mazaideh, S. M. Khalil, Amer. Chem. Sci. J., (2016); 13(4): 1. [14] G. M. Al-Mazaideh, R. A. Enwisry, S. M. Khalil, Inter. Res. J. Pure & Appl. Chem., (2016); 12(3): 1. [15] G. M. Al-Mazaideh, W. B. Ejlidi, S. M. Khalil, Inter. Res. J. Pure & Appl. Chem., (2016); 12(4): 1. [16] G. M. Al-Mazaideh, T. S. Ababneh, K. H. Abu-Shandi, R. M. A. Q. Jamhour, H. J. Ayaal Salman , A. M. AlMsiedeen, S. M. Khalil. Phy. Sci. Inter. J., (2016); 12(1): 1. [17] S. M. Khalil, G. M. Al-Mazaideh, N. M. Ali, Inter. J. Biochem. Res. & Rev., (2016), 14(2): 1. [18] M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. arone, B. Mennucci, G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P. Hratchian,A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, J. A. Montgomery, Jr., J. E. Peralta, F. Ogliaro, M. Bearpark, J. J. Heyd, E. Brothers, K. N. Kudin, V. N. Staroverov, R. Kobayashi, J. Normand, K. Raghavachari, A. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi,M. Cossi, N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, R. L. Martin, K. Morokuma, V. G. Zakrzewski, G. A. Voth, P. Salvador, J. J. Dannenberg, S. Dapprich, A. D. Daniels, O. Farkas, J. B. Foresman, J. V. Ortiz, J. Cioslowski, and D. J. Fox, Gaussian, Inc., Wallingford CT, 2009. [19] R.G. Pearson. Inorg. Chem., (1988); 27:734. [20] D. R. Lide. CRC Handbook of Chemistry and Physics, 88th ed., 2007-2008. [21] D. Humpola, H. S. Odetti, A. E. Fertitta, J. L. Vicente, J. Chil. Chem. Soc., (2013); 58(1):1541; E. E. Ali-Shattle, M. H. Mami, M. M. Alnaili, Asian J. Chem., (2009); 21: 5431. [22] G. Gece, S. Bilgic. Corros. Sci., (2009); 51 (8):1876. [23] S. M. Khalil. Z. Naturforsch., (2008); 63a: 42. [24] M. K. Awad, M. R. Mustafa, M. M. Abo Elnga. J. Mol. Struct. (Theorem)., (2010); 959: 66. Mor. J. Chem. 6 N°1 (2018) 195-202

201

[25] W. Kohn and L. J. Sham, Physical Re- view, (1965); 140 (4): A1133. [26] R. G. Parr, L. V. Szentpaly, S. Liu, J. Am. Chem. Soc., (1999); 121:1922; Lesar, I. Miosev, Chem. Phys., (2009); 483: 198. [27] N. I. Levine, “Quantum Chemistry,” Prentice Hall, En- glewood Cliffs, 2000. [28] X. Li, S. Deng, H. Fu, T. Li. Electrochem. Acta., (2009); 54: 4089. [29] T. Arslan, F. Kandemirli, E. E. Love. Alemu H. Corros. Sci., (2009); 51(1): 35. [30] M. Mihit, K. Laarej, H. Abou El Makarim, L. Bazzi, R. Salghi, Arab. J. Chem., (2010); 3: 55. [31] Singh, K. P., Malik, A., Sinha, S., Ojha, P., J. Hazard Mat., (2008); 150: 626. [32] Srihan, V., Das, A., Desalination, (2008); 225: 220. [33] Humpola, P. D., Odetti, H. S., Fertitta, A. F., Vicente, J. L., J. Chil. Chem. Soc., (2013);58:1541.

Mor. J. Chem. 6 N°1 (2018) 195-202

202