List of density functional .
The use of computational chemistry as a tool in the design and development of organic corrosion inhibitors has been greatly enhanced by the development of density functional theory (DFT) and Molecular dynamic simulation (MD). Experimentally corrosion inhibitor development requires lots of money and time. Thus, in the era of hardware and software development, corrosion scientist can select a potential inhibitor on the basis of theoretical analysis of molecular properties of inhibitor molecules, which have reduced the cost. DFT and MD are capable to accurately predict the inhibition characteristics of inhibitor molecules using molecular/electronic properties and reactivity indices. The purpose of this book chapter is to summarize some important features related to DFT and MD, giving a brief background to the selected DFT/MD-based chemical reactivity concepts, calculations and their applications to organic corrosion inhibitor design. The impact of this book chapter is to illustrate the enormous power of DFT and MD.
Corrosion of active metals such as iron (Fe), copper (Cu) and Aluminum (Al) can be reduced by modifying their surface using organic/polymeric corrosion inhibitors. The inhibitor molecules adsorbed onto the metal surface and form a protective barrier against corrosion . The corrosion protection ability of inhibitors depends upon the extent of interaction between inhibitor molecules and metal surface . The basics mechanism through which adsorbed inhibitor molecules alter the corrosion process consists of either changing the cathodic and anodic corrosion reactions or physically blocking the active sites present on metal surface or by .
Perusal of the literature clearly reveals that huge number of publication exits where experimental works have been carried out to understand the corrosion inhibition process, however use of theoretical studies including density functional theory (DFT) and molecular dynamic simulation (MD) appeared in recent years. In traditional experimental work scientists has to test huge number of organic compounds for selecting them as a potential corrosion inhibitor. This kind of search takes more time and money. However, in the era of development in the field of hardware and software technology make us capable to select organic compounds as potential corrosion inhibitor without expensing huge amount of money on chemicals. The corrosion inhibition efficiency of organic compounds based on molecular and electronic properties can be predicted accurately using MD and DFT calculations. In the recent years, scientist have devoted their positive efforts in the field of corrosion using theoretical calculation [4, 5, 6]. To our knowledge, the review exists in the literature only describing separately DFT and MD and not both in relation to corrosion inhibition study [7, 8]. Here we are trying to describe both DFT and MD on the same platform simultaneously.
In this chapter we will try to explore some recent studies and concept which can provide some new concept in the field of computational chemistry related to corrosion inhibition phenomenon. Here, we summarize the adsorption of organic molecules onto the metal surface especially iron, aluminum and copper in order to study their corrosion inhibition properties using DFT and MD simulations.
2. Basic concepts related to molecular modeling
In literature, large number of publication/papers exits where corrosion inhibition property of inhibitors has been studied using theoretical study. In general organic molecules have the ability to donate electrons to empty d-orbital of metal surface and form covalent bond which is known as coordinate bond. However, they also undergo back donation, where they accept the free electrons presenting in filled metal orbitals and this is called as retrodonating bonds. Thus, these kinds of bonds make the organic compounds a potential candidate for corrosion inhibition. The theoretical parameters like highest occupied molecular orbital (HOMO), lowest unoccupied molecular orbital (LUMO), dipole moment (
DFT study separately represents the organic compound and metal surface [12, 13, 14, 15, 16]. In DFT, the interaction between the inhibitor molecule and metal surface was not directly modeled. However, in MD the most reactive fragment of the molecule that has greater affinity to interact with the metal surface is study using ab initio method. This method provide a real picture which is actually going in the corrosion process like interaction between inhibitor molecule and metal surface, orientation of molecule towards metal and organization of molecules.
3. Inhibitor/surface interaction study using DFT
3.1. The basics of DFT: the Hohenberg-Kohn theorem
Quantum electrochemistry is the field which includes quantum mechanics, electrochemistry and electrodynamics. In general, quantum electrochemistry is an application of density functional theory (DFT) for studying the electrochemical processes like transfer of electron towards electrode surfaces .
There are so many books and articles that exist in literature which describe about the basis set used in DFT [18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31]. So, here we provide only very few pints of DFT without using mathematics.
The functional is a function of a function, i.e., electron. Therefore, on the basis of fundamental quantum mechanics and from the experimental results different functional have been developed and are listed in Table 1 .
|Xα||X alpha||Exchange only|
|HFS||Hartree-Fock Slater||HF with LDA exchange|
|VWN||Vosko, Wilks, and Nusair||LDA|
|BLYP||Becke correlation functional with Lee, Yang, Parr exchange||Gradient-corrected|
|B3LYP||Becke 3 term with Lee, Yang and Parr||Hybrid exchange|
|PW91||Perdew and Wang 1991||Gradient-corrected|
|B96||Becke 1996||Gradient corrected|
|B3P86||Becke exchange, Perdew correlation||Hybrid|
|B3PW91||Becke exchange, Perdew and Wang||Hybrid correlation|
There are several generations of density functional. First generation is known as Xα method and it is the simplest one. This functional was developed by J.C. Slater, who was working on the approximation of Hartree-Fock but unwillingly discovered the simplest for DFT. In
The functional of second generation uses both the density and its gradients. In 1986 first gradient-corrected energy functional was proposed in which Becke, Perdew and Wang proposed exchange and Perdew correlation [32, 33]. In the present time the most popular functionalis of Becke and Perdew for exchange and correlation respectively , also Lee et al.  for correlation and Perdew and Wang for exchange and correlation both . These functional are collectively known as “generalized gradient approximations (GGAs).”
The third generation is hybrid functional and it is more advanced than GGA. In this, the functional consists of both Hartree-Fock type exchange and DFT exchange calculated from the orbitals. In 1998, a hybrid functional called as B3LYP functional was introduced. In 1994, Gaussian introduced this functional in in computational package for first time and it is written by below equation:
3.3. Basis sets
In DFT calculation various basis sets are used . The simplest one is STO-3G. In this basis set, 1s is given by three Gaussians and 2s, 2px, 2py, and 2pz each by another three. Improvement in the basis set can obtain by incorporating two 1s functions for hydrogen and for 2nd row atoms like carbon two 2s and two 2p functions has incorporated. These functional are called as split-valence basis sets. Here valence and core orbitals are represented by two sets of functions and single set of functions respectively. An example for this:
Carbon 3-21G: Three Gaussians for 1s combination, same two Gaussians for 2s and 2p combination; plus 2s’ and 2p’ the same one Gaussian.
Carbon 6-31G: Six Gaussians for 1s combination, same three Gaussians for 2s and 2p combination; plus 2s’ and 2p’ the same one Gaussian.
In general it can be represented as
Further advancement generates triple-split-valence basis sets, for example 6-311G. In this basis set core functional consists of six Gaussians and three sets of valence functions containing three, one, and one Gaussians, respectively.
Multiple zeta basis sets such as double zeta (DZ) and triple zeta (TZ) basis sets has also introduced. The difference between multiple and split-valence basis sets is that in multiple basis sets all of the orbitals are separated into either two or three sets of functions, and so forth. Also, multiple basis sets used different
So, improvement in the electronic structure calculation can be achieved by adding functions corresponding to orbitals with a higher angular momentum. This concept can explained as follows: for hydrogen atom, p functions are added; for carbon, nitrogen, oxygen, etc., d functions are added; and for transition metals, f functions are introduced. These are represented by adding an asterisk to the basis set and also specifying the p, d, f functions. This can be understood as per below example:
6-31G∗ or 6-31G(d): Adds d functions to 2nd row elements (C, N, O, etc.).
6.31G∗∗ or 6-31G(d,p): Adds d function to 2nd row elements (C, N, O, etc.) and p functions to H.
The lone pair of electrons on the heteroatoms can be taken into confederation by adding diffuse functions and these are represented in the basis set using + and ++ signs. For hydrogen atom polarization and diffuse functions are not very import but for other atoms like carbon, nitrogen, oxygen it is important, for example,
6-31+G: represents that diffuse functions has been added to 2nd row elements (C, N, O, etc.).
6-31++G: represents that diffuse functions has added to 2nd row elements and H.
In the case of transition elements which contain large number of electrons, requires more time for calculation. So, the LanL2DZ basis set is used for transition metal calculations. This basis set is derived by combining the valence electrons using double zeta functions and Los Alamos ECP. The software used in corrosion inhibition study with a wide range of basis sets are Spartan, Material Studio, Gaussian 03, Gaussian 09 and most recently Gaussian 14 has been introduced. Table 2 represents commonly used basis sets .
|6-31G||+ ++ ∗ ∗∗||H–Cl|
|6-311G||+ ++ ∗ ∗∗||H–Ar|
|D95||+ ++ ∗ ∗∗||H–Cl|
|D95V||+ ++ ∗ ∗∗||H–Ne|
|SEC||+ ++ ∗ ∗∗||H–Cl|
|CEP-4G||+ ++ ∗ ∗∗||H–Cl|
|CEP-31G||+ ++ ∗ ∗∗||H–Cl|
|CEP-121G||+ ++ ∗ ∗∗||H–Cl|
4. DFT parameters and their application in corrosion inhibition study
4.1. Frontier molecular orbitals
The optical and electronic property of organic compounds is explained by analyzing their frontier molecular orbitals (FMO) like highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) . The interaction of organic compounds as corrosion inhibitor with the metal surface depends upon the electron donation ability of the organic compounds, i.e., ionization potential (
In literature Singh et al.  has correlated some important quantum chemical parameters of organic corrosion inhibitors using/B3LYP methods using 6-311G (d, p) basis set with their corrosion inhibition efficiency. They studied three triazine-based hydrazone derivatives and their molecular, orbital pictures and quantum chemical parameters are represented in Figures 1
From the figure, it could be said that HOMO regions are distributed on the triazine ring, phenyl ring containing electron donating group, nitrogen and oxygen atoms. These regions correspond to the electron donating areas and inhibitor molecules interact with the metal surface through these regions and undergo adsorption over it. However, the LUMO are distributed over the triazine and phenyl ring respectively and these are corresponding to the electron accepting regions. From these LUMO regions inhibitor accept electron from metal. Order of inhibition efficiency is HT-1 > HT-2 > HT-3.
4.2. Electronegativity and the electronic chemical potential
One of the important global molecular property calculated by density functional theory (DFT) which provide an idea for chemical reactivity and selectivity of organic compound as corrosion inhibitor are electronegativity (
4.3. Global hardness and softness
According to Parr et al. global hardness (
The reciprocal of global hardness is global softness (
4.4. Electric dipole polarizability
The electric dipole polarizability (
The polarizability (
4.5. Electrophilicity index
The global electrophilicity index (
The electrophilicity index accounts the ability of the molecule to accept electrons. So, higher the
4.6. The fraction of electrons transferred
When the inhibitor and metal interaction occurs, flow of electron takes place from the lower electronegativity molecule to the higher electronegativity metal. This transfer of electron continues until the chemical potential becomes equal. The fraction of electron transfer is given as follows :
where metal and inhibitor molecule are represented by
The work function (
4.7. Fukui function and local softness
In the field of corrosion research inhibition of corrosion is governed through donation and acceptance of electron, which basically involves the nucleophilic and electrophilic reactions. This can be achieved through the evaluation of the Fukui indices [61, 62, 63]. Thus, by invoking the HSAB principle in a local sense, one may establish the behavior of the different sites with respect to hard or soft reagents.
The Fukui function generated using the finite difference approximation is as follows:
Castro et al.  studied the series of isomers of the Schiff base (E)-2-(2-hydroxybenzylideneamino) phenylarsonic acid as corrosion inhibitor. The author has pictorially presented the negative Fukui, i.e.,
In Figure 4, the black region represents the area where inhibitor contains the major portion of electrons for donating to the empty d orbitals of Fe metal. The Fukui function represents the change of the electronic density in a given point with respect to the change in the number of electrons
Local softness is calculated is using the equation :
The integration of equation gives global softness value:
Local softness and Fukui function are related as per below equation:
The local softness contains information similar to those obtained from Fukui function plus additional in-formation about the molecular softness, which is related to the global reactivity with respect to a reaction partner.
4.8. Charge transfer: donation and back donation
The energy change in the molecule is attributed because of two processes either electron donation to the metal or acceptance from the metal, i.e., back donation. Gomez et al.  proposed a simple charge transfer model for donation and back donation. In this model when the molecule receives certain amount of charge (∆
However, when the molecule back donate certain amount of charge (∆
If the total change in energy (∆
In a situation when the total energy change becomes a minimum with respect to ∆
and change in total energy becomes as follows:
According to Guo et al.  the change in total energy are negative, i.e., ∆
5. Inhibitor/surface interaction study using molecular dynamic simulation (MD)
Molecular dynamics is a computer-based modeling technique by which the evolution as a function of time or trajectory of a molecule is described by the principles of classical Newtonian mechanics . In molecular dynamics, intramolecular movements of inhibitor molecules are simulated in order to visualize the real picture of corrosion inhibition process. Molecular dynamics simulation calculation has two serve two purposes:
Intramolecular movements’ simulation and calculation of thermodynamic properties such as entropy, free energy, etc.
Optimization of the molecular structures efficiently by avoiding multiple minima by annealing the simulation process. The range of simulation time is approximately between 10−14 and 10−10s per simulation.
Molecular dynamics simulation involves the implementation of initial steps: In the first very step build molecular structure are introduced or immersed in a box having solvent molecules, after that simulation process starts. The total energy of the molecular system should be minimizing before the start of the simulation, in order to avoid the generation of aberrant trajectories (when the initial forces are too great). For achieving this, the steepest slope/steepest descent method and that of the conjugated gradients have been frequently used.
6. Corrosion inhibitors studied using DFT
6.1. Adsorption behavior of inhibitors on mild steel surface
In literature, several DFT calculations are observed where inhibition efficiency and molecular structure/electronic properties of organic corrosion inhibitors have been correlated. Xu et al.  has studied corrosion inhibition property of Schiff’s base derivatives via DFT calculation using GGA/BLYP method. The result indicates that 2-PCT molecules adsorb strongly onto iron surface and provide higher inhibition efficiency.
Zhang et al.  investigated the inhibition of methionine (Met) and proline (Pro) in PCMs solution using DFT and B3LYP/6-31G (d) method. The inhibition efficiency of Met and Pro follow the order of Met > Pro.
The corrosion inhibition effects of thiourea (TU), methylthiourea (MTU) and phenylthiourea (PTU) on mild steel in 0.1MH2SO4 was investigated by Ozcan et al. . The electronic properties were calculated using ab initio RHF/6-31G(d) method. The energy of
The corrosion inhibition performance of three naphthyridines derivatives were explored for mild steel in 1 M HCl by Singh et al.  using DFT methods. The inhibition efficiencies of N-3 and N-2 are higher due to the presence of electron-donating –OCH3 and –CH3 substituents respectively as compared to N-1, which has no substituents.
Feng et al. , studied inhibition activity of an imidazoline derivative, namely 1-[N,N′-bis(hydroxylethylether)- aminoethyl]-2-stearicimidazoline (HASI) for carbon steel in 5% NaCl saturated with Ca(OH)2 solution using DMol3 module. DFT calculations reveals that N=C-N region in imidazoline ring is the most active reaction site.
6.2. Adsorption behavior of inhibitors on aluminum surface
The inhibitive action of poly ethylene glycol (PEG) against the corrosion of aluminum surface in 1 M HCl was investigated by Awad et al.  using B3LYP/6–31 + G(d,p) basis set. DFT was performed on repeating units of 1, 2, 3, 4 and 5. The inhibition performance increases with increase in number of repeating unit.
The author Li et al.  studied the corrosion inhibition effect of three oxime compounds on the corrosion of aluminum in HCl solution using BLYP in conjunction with double numerical plus d-functions (DND) basis set. In these inhibitors, HOMO and LUMO electron densities are localized principally on C=N–OH, which suggest that C=N–OH functional group could both accept and donate the electron.
Khaled et al. , investigated the corrosion effect of and adsorption characteristics of three imidazole derivatives on aluminum in 1.0 M HCl using local density functional (LDF) method with a double numeric polarization (DNP) basis set, and a Becke–Perdew (BP) functional. A good correlation between the rate of corrosion and EHOMO, as well as with energy gap (∆
Kaya et al.  applied HF and DFT/B3LYP methods with SDD, 6-31++G (d, p) and 6-31 G basis sets both in gaseous and aqueous phase for studying the adsorption and corrosion inhibition properties of some benzotriazole and phospono derivatives on aluminum. The inhibition efficiency ranking of these molecules as: PBA > PBTA>PAA > TBTA.
6.3. Adsorption behavior of inhibitors on copper surface
Corrosion inhibition of indazole (IA) and 5-aminoindazole (AIA) for copper in NaCl solution was studied by Qiang et al. . The molecular structures of IA and AIA were geometrically optimized by density functional theory (DFT) using B3LYP functional with 6-311++G(d,p) basis set in aqueous phase. The higher value of HOMO and lower value of ∆
Khalid and Madkour et al. [13, 80, 81] probe the interactions between benzotriazole, methionine, phenol derivatives and copper via quantum calculations respectively. The obtained results suggests that quantum parameters are very useful in characterizing organic compounds as an adsorbate.
7. Corrosion inhibitors studied using MD
Feng et al.  performed the MD simulation of 1-[N,N′-bis(hydroxylethylether)- aminoethyl]-2-stearicimidazoline (HASI). The MD calculation suggests that adsorbed imidazoline molecule is parallel to the iron surface in order to maximize its contact with the surface (Figure 5). The interaction energy for Fe atom is −284 kJ/mol, for Fe3O4 is −226 kJ/mol and for Fe2O3 is −157 kJ/mol. So, inhibitor adsorbed on Fe surface more strongly than Fe3O4 and Fe2O3.
MD calculation of adsorption of indazole (IA) and 5-aminoindazole (AIA) on copper surface was studied Qiang et al. . The optimized equilibrium configuration of inhibitors molecules are shown in Figure 6 and inhibitors adsorption occurs through parallel orientation. The adsorption energy between Cu (1 1 1) surface and inhibitor molecule are −250.86 kJ/mol for IA and −307.86 kJ/mol for AIA.
Molecular level MD calculation of Benzotriazole and Phospono derivatives was studied by Kaya et al.  on aluminum surface. The best adsorption configuration of inhibitor molecules on Al (111) surface are represented in Figure 7. The parallel configuration of inhibitor molecules suggests the stronger adsorption.
8. Conclusion and future work
Computational program provide an atomistic view for understanding the corrosion inhibition process and significantly contribute for the selection of adequate inhibitor molecules. DFT calculations has made possible to characterize a suitable inhibitor molecules using their intrinsic properties such as energies of HOMO and LUMO orbitals, energy gap (∆E), electronegativity, hardness, softness, etc. It helps corrosion scientist to understand the inhibitor-metal interaction through the locations of HOMO and LUMO orbitals. The examination of local reactivity reveals the nucleophilic and electrophilic centers which react with the metal surfaces.
The literature survey suggests that with help of quantum chemistry a good understanding of corrosion system can be achieved when various metal surfaces and various inhibitor molecules are compared. DFT has ability to perform calculations for complex or large organic molecules that can be applied as a potential corrosion inhibitor. In addition, DFT provide a platform to scientist for designing a unique inhibitor molecule and understand their molecular structure at atomistic level.
In recent years, computational modeling like DFT and MD provide real inhibitor-metal interaction picture by performing calculation in aqueous phase. Such an approach act as a benchmark for experimental study in the field of corrosion research.
Molecular dynamic simulation (MD) gives a representative image of inhibitor molecule orientation on the metal surface. This helps to analyze the inhibition activity of inhibitor molecules, i.e., inhibitor with parallel or flat orientation with respect to metal surface will cover more surface area and provide a better corrosion inhibition protection as compared to the inhibitor which have vertical orientation.
Future research includes modeling the electrochemical potential, generation of alloy surface layers, addition of thin oxide film covering onto the metal surface, calculation of more complex molecule, comprehension of solvent molecule role in adsorption, understanding the defect role in adsorption process, etc.
The authors are thankful to the Sichuan 1000 Talent Fund and financial assistance provided by the National Natural Science Foundation of China (No. 51274170).
Conflict of interest
There are no conflicts to declare.
Mao ZJ, Jun L. Corrosion inhibition performance of carbon steel in brine solution containing H2S and CO2 by novel Gemini surfactants. Acta Physico-Chimica Sinica. 2012; 28:623-629
Obot IB, Ebenso EE, Kabanda MM. Metronidazole as environmentally safe corrosion inhibitor for mild steel in 0.5 M HCl: Experimental and theoretical investigation. Journal of Environmental Chemical Engineering. 2013; 1:431-439
Trabanelli G. Chemical Industries Corrosion Mechanism. Vol. 3. New York: Marcel Dekker; 1987
Gutiérreza E, Rodrígueza JA, Borbollaa JC, Rodrígueza JGA, Thangarasu P. Development of a predictive model for corrosion inhibition of carbonsteel by imidazole and benzimidazole derivatives. Corrosion Science. 2016; 108:23-35
Shaban SM. N-(3-(dimethyl benzyl ammonio)propyl) alkanamide chloride derivatives as corrosion inhibitors for mild steel in 1 M HCl solution: Experimental and theoretical investigation. RSC Advances. 2016; 6:39784-39800
Ramaganthan B, Gopiraman M, Olasunkanmi LO, Kabanda MM, Yesudass S, Bahadur I, et al. Synthesized photo-cross-linking chalcones as novel corrosion inhibitors for mild steel in acidic medium: Experimental, quantum chemical and Monte Carlo simulation studies. RSC Advances. 2015; 5:76675-76688
Obot IB, Macdonald DD, Gasem ZM. Density functional theory (DFT) as a powerful tool for designing new organic corrosion inhibitors. Part 1: An overview. Corrosion Science. 2015; 99:1-30
Verma C, Lgaz H, Verma DK, Ebenso EE, Bahadur I, Quraishi MA. Molecular dynamics and Monte Carlo simulations as powerful tools for study of interfacial adsorption behavior of corrosion inhibitors in aqueous phase: A review. Journal of Molecular Liquids. 2018; 260:99-120
Yang W, Parr RG, Pucci R. Electron density, Kohn–Sham frontier orbitals, and Fukui functions. The Journal of Chemical Physics. 1954; 81:2862-2863
Parr RG, Yang W. Density functional approach to the frontier-electron theory of chemical reactivity. Journal of the American Chemical Society. 1984; 106:4049-4050
Armaković S, Armaković SJ, Vraneš M, Tot A, Gadžurić S. DFT study of 1-butyl-3-methylimidazolium salicylate: A third-generation ionic liquid. Journal of Molecular Modeling. 2015; 21:246-256
Guo L, Zhang ST, Lv TM, Feng WJ. Comparative theoretical study on the corrosion inhibition properties of benzoxazole and benzothiazole. Research on Chemical Intermediates. 2015; 41:3729-3742
Madkour LH, Elshamy IH. Experimental and computational studies on the inhibition performances of benzimidazole and its derivatives for the corrosion of copper in nitric acid. International Journal of Industrial Chemistry. 2016; 7:195-221
Mousavi M, Mohammadalizadeh M, Khosravan A. Theoretical investigation of corrosion inhibition effect of imidazole and its derivatives on mild steel using cluster model. Corrosion Science. 2011; 53:3086-3091
Camacho RL, Montiel E, Jayanthi N, Pandiyan T, Cruz J. DFT studies of alpha-diimines adsorption over Fe-n surface (n = 1, 4, 9 and 14) as a model for metal surface coating. Chemical Physics Letters. 2010; 485(1–3):142-151
Wang Y, Cheng X, Yang X, Yang X. DFT study of solvent effects for some organic molecules using a polarizable continuum model. Journal of Solution Chemistry. 2006; 35:869-878
Dogonadze RR. Theory of molecular electrode kinetics. In: Hush NS, editor. Reactions of Molecules at Electrodes. London: Interscience Pub.; 1971. pp. 135-227
Sholl D, Steckel JA. Density Functional Theory: A Practical Introduction. New York: Wiley; 2009
Enge E l, Dreizler RM. Density Functional Theory: An Advanced Course. Berlin, Heidelberg: Springer-Verlag; 2011
Parr RG, Weitao Y. Density-Functional Theory of Atoms and Molecules. USA: Oxford University Press; 1989
Marques MAL, Maitra NT, Nogueira FMS, Gross EKU, Rubio A. Fundamentals of Time-Dependent Density Functional Theory. Berlin, Heidelberg: Springer-Verlag; 2012
Fiolhais C, Marques MAL, Nogueira F. A Primer in Density Functional Theory. Berlin, Heidelberg: Springer-Verlag; 2003
Sahni V. Quantal Density Functional Theory. Berlin, Heidelberg: Springer Science & Business Media; 2004
Tsuneda T. Density Functional Theory in Quantum Chemistry. Japan: Springer; 2014
Gross E, Dreizler RM. Density Functional Theory: An Approach to the Quantum Many-Body Problem. Berlin, Heidelberg: Springer; 1990
Gross EKU, Dreizler RM. Density Functional Theory. New York: Springer Science & Business Media; 1995
Amador C. Density Functional Theory. Berlin, Heidelberg: Springer-Verlag; 1983
Hohenberg P, Kohn W. Inhomogeneous electron gas. Physics Review B. 1964; 136:864-871
Koch W, Holthausen MC. A Chemist’s Guide to Density Functional Theory. Weinheim: Wiley-VCH; 2000
Bel SL, Dines TJ, Chowdhry BZ, Withnall R. Computational chemistry using modern electronic structure methods. Journal of Chemical Education. 2007; 84:1364-1368
Kohn W, Sham L. Self-consistent equations including exchange and correlation effects. Physics Review. 1965; 140:A1133-A1138
Perdew JP. Density functional approximation for the correlation energy of an inhomogeneous electron gas. Physical Review B. 1986; 33:8822-8824
Becke AD. Density functional calculations of molecular bond energies. The Journal of Chemical Physics. 1986; 84:4524-4529
Lee C, Yang W, Parr RG. Development of the Colle-Salvetti correlation energy formula into a functional of the electron density. Physical Review B. 1988; 37:785-789
Perdew J, Yang W. Accurate and simple density functional for the electronic exchange energy: Generalized gradient approximation. Physical Review B. 1986; 33:8800-8802
Becke AD. Density–functional exchange-energy approximation with correct asymptotic behavior. Physical Review A. 1988; 38:3098-3100
Ramachandran KI, Deepa G, Namboori K. Computational chemistry and molecular modeling, Principles and Applications. Berlin, Heidelberg: Springer; 2008;137
Arshad MN, Bibi A, Mahmood T, Asiri AM, Ayub K. Crystal structures and spectroscopic properties of triazine based hydrazone derivatives; a comparative experimental-theoretical study. Molecules. 2015; 20:5851-5874. DOI: 10.3390/molecules20045851
Zhan CG, Nichols JA, Dixon DA. Ionization potential, electron affinity, electronegativity, hardness, and electron excitation energy: Molecular properties from density functional theory orbital energies. The Journal of Physical Chemistry. A. 2003; 107:4184-4195
Singh A, Ansari KR, Kumar A, Liu W, Songsong C, Lin Y. Electrochemical, surface and quantum chemical studies of novel imidazole derivatives as corrosion inhibitors for J55 steel in sweet corrosive environment. Journal of Alloys and Compounds. 2017; 712:121-133
Tian H, Li W, Cao K, Hou B. Potent inhibition of copper corrosion in neutral chloride media by novel non-toxic thiadiazole derivatives. Corrosion Science. 2013; 73:281-291
Hegazy MA, Badawi AM, Abd El Rehim SS, Kamel WM. Corrosion inhibition of carbon steel using novel N-(2-(2-mercaptoacetoxy)ethyl)-N,N-dimethyl dodecan-1-aminium bromide during acid pickling. Corrosion Science. 2013; 69:110-122
Kovacevic N, Kokalj A. The relation between adsorption bonding and corrosion inhibition of azole molecules on copper. Corrosion Science. 2013; 73:7-17
Deng S, Li X, Xie X. Hydroxymethyl urea and 1,3-bis(hydroxymethyl) urea as corrosion inhibitors for steel in HCl solution. Corrosion Science. 2014; 80:276-289
Singh A, Ansari KR, Haque J, Dohare P, Lgaz H, Salghi R, et al. Effect of electron donating functional groups on corrosion inhibition of mild steel in hydrochloric acid: Experimental and quantum chemical study. Journal of the Taiwan Institute of Chemical Engineers. 2018; 82:233-251
Mendoza-Huizar LH, Rios-Reyes CH. Chemical reactivity of atrazine employing the Fukui function. Journal of the Mexican Chemical Society. 2011; 55:142-147
Chermette H. Chemical reactivity indexes in density functional theory. Journal of Computational Chemistry. 1999; 20:129-154
Pearson RG. Absolute electronegativity and hardness correlated with molecular orbital theory. Proceedings of the National Academy of Sciences. 1986; 83:8440-8441
Parr RG, Gadre SR, Bartolotti LJ. Local density functional theory of atoms and molecules. Proceedings of the National Academy of Sciences of the United States of America. 1979; 76:2522-2526
Haque J, Ansari KR, Srivastava V, Quraishi MA, Obot IB. Pyrimidine derivatives as novel acidizing corrosion inhibitors for N80 steel useful for petroleum industry: A combined experimental and theoretical approach. Journal of Industrial and Engineering Chemistry. 2017; 49:176-188
Zevallos J, Labbé AT. A theoretical analysis of the Kohn-Sham and Hartree-Fock orbitals and their use in the determination of electronic properties. Journal of the Chilean Chemical Society. 2003; 48:39-47
Ghanty TK, Ghosh SK. Correlation between hardness, polarizability, and size of atoms, molecules, and clusters. The Journal of Physical Chemistry. 1993; 97:4951-4953
Chattaraj PK, Sengupta S. Popular electronic structure principles in a dynamical context. The Journal of Physical Chemistry. 1996; 100:16126-16130
Hohm U. Is there a minimum polarizability principle in chemical reactions? The Journal of Physical Chemistry. A. 2000; 104:8418-8423
Parr RG, Sventpaly L, Liu S. Electrophilicity index. Journal of the American Chemical Society. 1999; 121:1922-1924
Geerlings P, De Proft F, Langenaeker W. Conceptual density functional theory. Chemical Reviews. 2003; 103:1793-1873
Pearson RG. Absolute electronegativity and hardness: Application to inorganic chemistry. Inorganic Chemistry. 1988; 27:734-740
Kovacˇević N, Kokalj A. DFT study of interaction of azoles with Cu(111) and Al(111) surfaces: Role of azole nitrogen atoms and dipole-dipole interactions. Journal of Physical Chemistry C. 2011; 115:24189-24197
Kokalj A. On the HSAB based estimate of charge transfer between adsorbates and metal surfaces. Chemical Physics. 2012; 393:1-12
Cao Z, Tang Y, Cang H, Xu J, Lu G, Jing W. Novel benzimidazole derivatives as corrosion inhibitors of mild steel in the acidic media. Part II: Theoretical studies. Corrosion Science. 2014; 83:292-298
Kokalj A. Is the analysis of molecular electronic structure of corrosion inhibitors sufficient to predict the trend of their inhibition performance. Electrochimica Acta. 2010; 56:745-755
Lewars EG. Computational Chemistry Introduction to the Theory and Applications of Molecular and Quantum Mechanics. London: Springer; 2011
Oláh J, Alsenoy CV. Condensed Fukui functions derived from stockholder charges: Assessment of their performance as local reactivity descriptors. The Journal of Physical Chemistry. A. 2002; 106:3885-3890
Castro ME, Percino MJ, Cerón M, Soriano G, Chapela VM. Theoretical inhibition efficiency study of Schiff Base (E)-2-(2-hydroxybenzylideneamino) phenylarsonic acid and its isomers. International Journal of Electrochemical Science. 2014; 9:7890-7903
López P, Méndez F. Fukui function as a descriptor of the imidazolium protonated cation resonance hybrid structure. Organic Letters. 2004; 6:1781-1783
Yang W, Parr RG. Hardness, softness and the Fukui function in the electronic theory of metals and catalysis. Proceedings of the National Academy of Sciences of the United States of America. 1985; 82:6723-6726
Gomez B, Likhanova NV, Dominguez-Aguilar MA, Martinez- Palou R, Vela A, Gazquez JL. Quantum chemical study of the inhibitive properties of 2-pyridyl-azoles. The Journal of Physical Chemistry. B. 2006; 110:8928-8934
Guo L, Ren X, Zhou Y, Xu S, Gong Y, Zhang S. Theoretical evaluation of the corrosion inhibition performance of 1,3-thiazole and its amino derivatives. Arabian Journal of Chemistry. 2017; 10:121-130
Rives JT, Jorgensen WL. Molecular dynamics of proteins with the OPLS potential functions. Simulation of the third domain of silver pheasant ovomucoid in water. Journal of the American Chemical Society. 1990; 112:2773-2781
Xu B, Yang W, Liu Y, Yin X, Gong W, Chen Y. Experimental and theoretical evaluation of two pyridinecarboxaldehyde thiosemicarbazone compounds as corrosion inhibitors for mild steel in hydrochloric acid solution. Corrosion Science. 2014; 78:260-268
Zhang Z, Tian N, Zhang W, Huang X, Ruan L, Wu L. Inhibition of carbon steel corrosion in phase-change-materials solution by methionine and proline. Corrosion Science. 2016; 111:675-689
Ozcan M, Dehri I, Erbil M. Organic sulphur-containing compounds as corrosion inhibitors for mild steel in acidic media: Correlation between inhibition efficiency and chemical structure. Applied Surface Science. 2004; 236:155-164
Singh P, Ebenso EE, Olasunkanmi LO, Obot I, Quraishi M. Electrochemical, theoretical, and surface morphological studies of corrosion inhibition effect of green naphthyridine derivatives on mild steel in hydrochloric acid. Journal of Physical Chemistry C. 2016; 120:3408-3419
Feng L, Yang H, Wang F. Experimental and theoretical studies for corrosion inhibition of carbon steel by imidazoline derivative in 5% NaCl saturated Ca(OH)2 solution. Electrochimica Acta. 2011; 58:427-436
Awad M, Metwally M, Soliman S, El-Zomrawy A. Experimental and quantum chemical studies of the effect of poly ethylene glycol as corrosion inhibitors of aluminum surface. Journal of Industrial and Engineering Chemistry. 2014; 20:796-808
Li X, Deng S, Xie X. Experimental and theoretical study on corrosion inhibition of oxime compounds for aluminium in HCl solution. Corrosion Science. 2014; 81:162-175
Khaled K, Amin MA. Electrochemical and molecular dynamics simulation studies on the corrosion inhibition of aluminum in molar hydrochloric acid using some imidazole derivatives. Journal of Applied Electrochemistry. 2009; 39:2553-2568
Kaya S, Banerjee P, Saha SK, Tüzün B, Kaya C. Theoretical evaluation of some benzotriazole and phospono derivatives as aluminum corrosion inhibitors: DFT and molecular dynamics simulation approaches. RSC Advances. 2016; 6:74550-74559
Qiang Y, Zhang S, Xu S, Li W. Experimental and theoretical studies on the corrosion inhibition of copper by two indazole derivatives in 3.0% NaCl solution. Journal of Colloid and Interface Science. 2016; 472:52-59
Khaled K, Amin MA, Al-Mobarak N. On the corrosion inhibition and adsorption behaviour of some benzotriazole derivatives during copper corrosion in nitric acid solutions: A combined experimental and theoretical study. Journal of Applied Electrochemistry. 2010; 40:601-613
Khaled K. Corrosion control of copper in nitric acid solutions using some amino acids—A combined experimental and theoretical study. Corrosion Science. 2010; 52:3225-3234