Theoretical study of the inhibitive properties of a group of benzimidazoles

In this work, a group of benzimidazoles, namely benzimidazole (BIM), 2‐hexyl benzimidazole (2‐HBIM), and 2‐benzyl chloride benzimidazole (2‐ClBBIM), have been studied using density functional theory (DFT) at the level of B3LYP with the 6‐31G (d) and 6‐311G (d, p) base sets in order to elucidate the different inhibition efficiencies and reactive sites of these compounds as corrosion inhibitors. Based on the calculated results, the studied molecules interact with Fe‐atoms mainly in their stable pyridine‐N protonated forms. It is found that the inhibition efficiencies of the stable protonated molecules involved in this study have a certain relation with some parameters, such as the highest occupied molecular orbital energy (E_HOMO) and the global hardness (η). The results indicate that the C (7) atoms, as the reactive sites, receive electrons from Fe‐atoms, benzene ring denote electrons to vacant orbital of Fe‐atoms. Additionally, Fukui indices and the orbital analyses suggest that 2‐ClBBIM has the highest reaction activity among the three molecules, the efficiency order of three inhibitors is found to be 2‐ClBBIM > 2‐HBIM > BIM, which accords with experimental results. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Quantum chemical studies on some thiadiazolines as corrosion inhibitors for mild steel in acidic medium

Density functional theory at the B3LYP/6-31G(d,p) basis set level was performed on three thiadiazolines, namely 4-chloro-N-(5-phenyl-1,3,4-thiadiazol-2(3H)-ylidene)aniline (TD01), 4-chloro-N-(5-(4-methoxyphenyl)-1,3,4-thiadiazol-2(3H)-ylidene)aniline (TD02), and 2-(5-(4-chlorophenylimino)-4,5-dihydro-1,3,4-thiadiazol-2-yl) phenol (TD03), and the inhibitive effect of these thiadiazolines against the corrosion of mild steel in acidic medium is elucidated. The calculated quantum chemical parameters correlated to the inhibition efficiency are E_HOMO (highest occupied molecular orbital energy), E_LUMO (lowest unoccupied molecular orbital energy), the energy gap (ΔE) hardness (η), softness (S), dipole moment (μ), electron affinity (EA) ionization potential (IE), the absolute electro negativity (χ), and the fraction of electron transferred (ΔN). The decreasing order of %IE of the thiadiazolines studied was found to be in agreement with experimental corrosion inhibition efficiencies. The local reactivity has been analyzed through the condensed Fukui function and local softness indices using population analysis.     

Theoretical studies of some sulphonamides as corrosion inhibitors for mild steel in acidic medium

Density functional theory (DFT) at the B3LYP/6‐31G (d,p) and BP86/CEP‐31G* basis set levels and ab initio calculations using the RHF/6‐31G (d,p) methods were performed on four sulfonamides (namely sulfaacetamide (SAM), sulfapyridine (SPY), sulfamerazine (SMR), and sulfathiazole (STI)) used as corrosion inhibitors for mild steel in acidic medium to determine the relationship between molecular structure and their inhibition efficiencies (%IE). The order of inhibition efficiency obtained was SMR > SPY > STI > SAM which corresponded with the order of most of the calculated quantum chemical parameters namely E_HOMO (highest occupied molecular orbital energy), E_LUMO (lowest unoccupied molecular orbital energy), the energy gap (ΔE), the Mulliken charges on the C, O, N, S atoms, hardness (η), softness (S), polarizability (α), dipole moment (μ), total energy change (ΔE_T), electrophilicity (ω), electron affinity (A), ionization potential (I), the absolute electronegativity (χ), and the fraction of electrons transferred (ΔN). Quantitative structure activity relationship (QSAR) approach has been used and a correlation of the composite index of some of the quantum chemical parameters was performed to characterize the inhibition performance of the sulfonamides studied. The results showed that the %IE of the sulfonamides was closely related to some of the quantum chemical parameters but with varying degrees/order. The calculated %IE of the sulfonamides studied was found to be close to their experimental corrosion inhibition efficiencies. The experimental data obtained fits the Langmuir adsorption isotherm. The negative sign of the E_HOMO values and other thermodynamic parameters obtained indicates that the data obtained supports physical adsorption mechanism. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Theoretical approach to the corrosion inhibition efficiency of some quinoxaline derivatives of steel in acid media using the DFT method

Corrosion inhibition efficiencies of 1,4-dihydroquinoxaline-2,3-dione (Q1) and 2-phenylthieno[2,3-b]quinoxaline (Q2) as corrosion inhibitors against the corrosion of steel surface in hydrochloric acid is studied by means of density functional approach B3LYP/6-31G calculations. Quantum chemical parameters such as highest occupied molecular orbital energy (E _HOMO), lowest unoccupied molecular orbital energy (E _LUMO), energy gap (ΔE), dipole moment (μ), electronegativity (χ), electron affinity (A), global hardness (η), softness (σ), ionization potential (I), the fraction of electrons transferred (?N), the global electrophilicity ω, and the total energy were calculated. All calculations have been performed by considering density functional theory using the GAUSSIAN03W suite of programs.     

苯并咪唑类缓蚀剂缓蚀性能的理论评价

采用量子化学计算和分子动力学模拟相结合的方法, 对2-巯基苯并咪唑(A)、2-氨基苯并咪唑(B)、2-甲基苯并咪唑(C)和苯并咪唑(D)等四种缓蚀剂抑制HCl对碳钢腐蚀的性能进行理论评价, 并对其缓蚀机理进行分析. 全局活性指数的计算表明, 四种分子中, 2-巯基苯并咪唑分子具有最强的反应活性; 对于其他三种分子, Fukui指数和全电子密度分布指出, 2-氨基苯并咪唑具有两个亲电攻击中心, 可在金属表面形成双中心吸附, 其缓蚀性能应优于2-甲基苯并咪唑和苯并咪唑; 缓蚀剂分子与三层铁原子表面相互作用的分子动力学模拟进一步确认2-甲基苯并咪唑比苯并咪唑在金属表面吸附更稳定. 综合量子化学计算和分子动力学模拟的计算结果, 四种缓蚀剂分子缓蚀效率的顺序应为A>B>C>D, 缓蚀性能的理论评价结论与实验结果相吻合.     

Interfacial Behavior of Cysteine between Mild Steel and Sulfuric Acid as Corrosion Inhibitor

Interfacial behavior of cysteine (Cys) between mild steel and sulfuric acid solution as a corrosion inhibitor has been studied with electrochemical AC (alternating current) and DC (direct current) techniques at (25.0±0.1) °C. The AC impedance results were evaluated using equivalent circuits in which a constant phase element (CPE) has been replaced with double layer capacitance (C_dl) to represent the frequency distribution of experimental data. Changes in impedance parameters (charge transfer resistance and double layer capacitance) indicated that cysteine molecules acted by accumulating at the metal/solution interface. The fractional coverage of the metal surface () was determined using AC impedance results and it was found that the adsorption of cysteine on the mild steel surface followed a Langmuir isotherm model with a standard free energy of adsorption (ΔG⁰_ads) of −35.1 kJ·mol⁻¹.To clarify the type of interaction between mild steel surface and cysteine molecules with a molecular orbital approach, electronic properties, such as, the highest occupied molecular orbital (HOMO) energy, the lowest unoccupied molecular orbital (LUMO) energy, and the frontier molecular orbital coefficients have been calculated. Energy gaps for the interaction of mild steel surface and cysteine molecules (E_{LUMO Fe}−E_{HOMO Cys} and E_{LUMO Cys}−E_{HOMO Fe}) were used to determine whether cysteine molecules acted as electron donors or electron acceptors when they interacted with the mild steel surface. The local reactivity was evaluated through the condensed Fukui indices. Theoretical calculations were carried out using the density functional theory (DFT) at B3LYP level with the 6-311++G(d,p) basis set for all atoms by Gaussian 03W program.     

New monodentate and bidentate silylene ligands by DFT

Considering the importance of silylene ligands in transition metal-mediated catalytic reactions, we have scrutinized eight novel monodentate ( 1 – 4 ) and bidentate ( 1 ^′ - 4 ^′ ) derivatives of 2,5-diX-cyclopentasilylene-2,4-dienes (X = NH₂, OH, PH₂, and SH), at M06/6-311++G** level of theory. To probe the complexation ability of our scrutinized silylene ligands with Rh atom ( 1 _Rh - 4 _Rh and 1 ^′ _Rh - 4 ^′ _Rh , respectively), thermodynamic and structural parameters such as complexation energy (ΔE_Com), singlet-triplet energy gap (ΔE_s-t), bond length, along with NBO and atoms in molecules analyses are provided. In going from less electron donating groups (EDGs) to more EDGs (NH₂ > OH > PH₂ > SH), the σ -donor strength and ligand flexibility increase. Structures 1 and 1 ^′ turn out as the most nucleophilic species for showing the highest nucleophilicity (N = 5.47 and 5.40 eV, respectively). Furthermore, they exhibit the highest proton affinity values (PA = 271.46 and 271.23 kcal/mol, respectively). The results indicate that bidentate coordination mode of silylene leads to a stronger Si-Rh complex. The overall orders of σ -donation ability for monodentate and bidentate silylene ligands are 1 > 2 > 3 > 4 and 1 ^′ > 2 ^′ > 3 ^′ > 4 ^′ , respectively.     

Electronic structure and quantum chemical analysis of the corrosion inhibition efficiency of quinoxalines

A theoretical investigation of the corrosion inhibition effectiveness of 1-[4-acetyl-2-(4-chlorophenyl) quinoxalin-1(4H)-yl] acetone (Q1), 2-(4-(2-ethoxy-2-oxoethyl)-2-p-tolylquinoxalin-1(4H)-yl) acetate (Q2) and 2-(4-methylphenyl)-1,4-dihydroquinoxaline (Q3) was evaluated by using quantum chemical parameters from density functional theory (DFT) with 6–311++ G (d, p) basis set at B3LYP level. Several quantum chemical parameters were determined to evaluate the array of selected molecules such as (E_LUMO), energy gap (ΔE), (E_HOMO), hardness, ionization potential, electronegativity, dipole moment (μ), the fraction of the electrons transferred from the inhibitor to the metal surface (ΔN), the softness (σ) and the total energy (TE). Theoretical data were found to confirm experimental results. By using these different quantitative chemical parameters to determine corrosion inhibition efficiency, we compare the results of the recent corrosion investigation of these inhibitors.     

Mismatch base pairing of the mutagen 8‐oxoguanine and its derivatives with adenine: A theoretical search for possible antimutagenic agents

Molecular geometries of 8‐oxoguanine (8OG), those of its substituted derivatives with the substitutions CH₂, CF₂, CO, CNH, O, and S in place of the N7H7 group, adenine (A), and the base pairs of 8OG and its substituted derivatives with adenine were optimized using the RHF/6‐31+G* and B3LYP/6‐31+G* methods in gas phase. All the molecules and their hydrogen‐bonded complexes were solvated in aqueous media employing the polarized continuum model (PCM) of the self‐consistent reaction field (SCRF) theory using the RHF/6‐31+G* and B3LYP/6‐31+G* methods. The optimized geometrical parameters of the 8OG‐A base pair at the RHF/6‐31+G* and B3LYP/6‐31+G* levels of theory agree satisfactorily with those of an oligonucleotide containing the base pair found from X‐ray crystallography. The pattern of hydrogen bonding in the CF₂‐ and O‐substituted 8OG‐A base pair is of Watson–Crick type and that in the unsubstituted and CH₂‐, CNH‐, and S‐substituted base pairs is of Hoogsteen type. In the CO‐substituted base pair, the hydrogen bonding pattern is of neither Watson–Crick nor Hoogsteen type. The CF₂‐substitution appears to introduce steric hindrance for stacking of DNA bases. On the basis of these results, it appears that among all the substituted 8OG molecules considered here, the O‐substituted derivative may be useful as an antimutagenic drug. It is, however, subject to experimental verification. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Theoretical Study of Molecular Structure, Tautomerism, and Geometrical Isomerism of N-Methyl- and N-Phenyl-Substituted Cyclic Imidazolines, Oxazolines, and Thiazolines

The geometries of various tautomers and isomers of 2-methylamino-2-imidazoline, 2-methylamino-2-oxazoline, 2-methylamino-2-thiazoline, 2-phenylamino-2-imidazoline, 2-phenylamino-2-oxazoline, and 2-phenylamino-2-thiazoline have been studied using the Becke3LYP/6–31+G(d,p) DFT, ONIOM(Becke3LYP/6–31+G(d,p):HF/3–21G^*) and ONIOM(Becke3LYP/6–31+G(d,p):AM1) methods. The optimized geometries indicate that these molecules show a distinctly nonplanar configuration of the cyclic moieties. In the gas phase, the amino tautomers (with exception of 2-phenylamino-2-imidazoline) are computed to be more stable than the imino tautomers. Of the two possible (E and Z) isomers of methyl and phenyl derivatives of imino-oxazolidine and imino-thiazolidine species, the (Z) isomers have the lowest energy. The iminozation free energies in the gas phase were found to be 5 – 15 kJ/mol. Absolute values of K _T depend strongly on the accuracy of the method used for calculation of free energy. Solvation (using the MD simulations) causes, in most cases, a shift in tautomeric preference toward the imino species.     

An experimental and theoretical investigation of the regio- and stereoselectivity of the polar [3+2] cycloaddition of azomethine ylides to nitrostyrenes

The regio- and stereochemical polar [3+2] cycloaddition of the azomethine ylides, which were generated in situ by the reaction of isatin derivatives and proline, with trans-β-nitrostyrene and (E)-1-phenyl-2-nitropropene were studied using experimental and theoretical methods. In comparison with trans-β-nitrostyrene, when the reactions were performed with (E)-1-phenyl-2-nitropropene, a remarkable inversion in the regioselectively was observed. The regioselectivity of the reactions was investigated using global and local reactivity indices and frontier molecular orbital (FMO) analysis at the B3LYP/6-31G(d,p) level of theory. The effects of the electronic and steric factors on the regioselectivity of the reactions were discussed. The inspection of geometries and energetics of transition states revealed the importance of weak interactions in regioselectivity of the cycloaddition reactions.     

A theoretical study on electronic structure and structure–activity properties of novel drug precursor 6‐acylbenzothiazolon derivatives

In this work, electronic properties and structure–activity relationship (SAR) parameters of 20 novel drug precursor 6‐acylbenzothiazolon derivatives with analgesic activity have been investigated theoretically by performing Austin Model‐1 (AM1) and DFT‐B3LYP/6‐31G (d) calculations with the aim to correlate the properties of each substance—particularly electronic properties and SAR parameters—with the biological interactions that are linked to their pharmacological effects. Their molecular properties were related to the biological activity of these drug precursor molecules. The relationship between octanol–water partition coefficient (log P) and each of the SAR parameters [E_LUMO–HOMO, molecular volume (V_m), ionization potential (IP), electron affinity, electronegativity (χ), chemical hardness (η), chemical softness (S), electrophilic index (ω), and molar refractivity] present linear correlation except for IP and χ. This result suggests that there are future prospects for designing or developing new drugs based on the correlation between the theoretically calculated parameters. According to AM1 calculation, the values of heat of formation of 6‐acylbenzothiazolon derivatives are negative (exothermic), which shows that these molecules are thermodynamically stable. E_LUMO–HOMO energy levels of the studied molecules are 4–5 eV, which also indicate that they are kinetically unstable. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Steric effects on the dialkyl substituted X2C2Si silylenes: A theoretical study

With the aim of recognizing the steric effects on the silylenic H₂C₂Si structures, ab initio and DFT calculations are carried out on 24 structures of X₂C₂Si (where X is hydrogen (H), methyl (Me), isopropyl (i‐pro), and tert‐butyl (tert‐Bu)). These species are at either triplet (t) or singlet (s) states. They are confined to the following three sets of structures ( 1 _X, 2 _X and 3 _X). Structures 1 _X include silacyclopropenylidenes ( 1 _s‐H and 1 _t‐H) and their 2,3‐disubstituted derivatives ( 1 _t‐Me, 1 _s‐Me; 1 _t‐i‐pro, 1 _s‐i‐pro; 1 _{t‐tert‐Bu}, 1 _{s‐tert‐Bu}). Structures 2 _X include vinylidenesilylenes ( 2 _s‐H and 2 _t‐H) and their 3,3‐disubstituted derivatives ( 2 _t‐Me, 2 _s‐Me; 2 _t‐i‐pro, 2 _s‐i‐pro; 2 _{t‐tert‐Bu}, 2 _{s‐tert‐Bu}). Structures 3 _X include ethynylsilylenes ( 3 _s‐H and 3 _t‐H) and their 1,3‐disubstituted derivatives ( 3 _t‐Me, 3 _s‐Me; 3 _t‐i‐pro, 3 _s‐i‐pro; 3 _{t‐tert‐Bu}, 3 _{s‐tert‐Bu}). Singlet–triplet energy separations (Δ E_{s‐t, X}) and relative energies for the above structures are acquired at HF/6‐31G*, B₁LYP/6‐31G*, B3LYP/6‐31G*, MP2/6‐31G*, HF/6‐31G**, B1LYP/6‐31G**, B3LYP/6‐31G**, and MP2/6‐31G** levels of theory. The highest Δ E_{s‐t, X} is encountered for 1 _X. All singlet states of X₂C₂Si, are more stable than their corresponding triplet states. Linear correlations are found between the LUMO–HOMO energy gaps of the singlet 1 _s‐X and 2 _s‐X with their corresponding singlet–triplet energy separations calculated at B3LYP/6‐31G**. The seven structures 2 _s‐Me, 2 _t‐Me, 3 _s‐Me, 1 _t‐Me, 1 _s‐Me, 1 _{s‐tert‐Bu}, and 3 _{t‐tert‐Bu} do not appear to be real isomers. Different stability orders are obtained as a function of the substituents (X). The order of stability for six isomers of H₂C₂Si is 1 _s‐H > 2 _s‐H > 3 _s‐H > 2 _t‐H > 3 _t‐H > 1 _t‐H. Replacing hydrogen atoms by methyl group (X = Me) presents a new stability order: 1 _s‐Me > 3 _s‐Me > 2 _s‐Me > 3 _t‐Me > 2 _t‐Me > 1 _t‐Me; and for (i‐pro)₂C₂Si is 1 _s‐i‐pro > 2 _s‐i‐pro ≈ 3 _s‐i‐pro > 3 _t‐i‐pro ≈ 2 _t‐i‐pro > 1 _t‐i‐pro. Using the larger tert‐butyl group as a substituent (X), yet it offers a more different stability order for six structures of (tert‐Bu)₂C₂Si: 1 _{s‐tert‐Bu} > 3 _{s‐tert‐Bu} > 2 _{s‐tert‐Bu} > 3 _{t‐tert‐Bu} > 1 _{t‐tert‐Bu} > 2 _{t‐tert‐Bu}. Among eight levels employed, B3LYP/6‐31G** appears as the method of choice. © 2006 Wiley Periodicals, Inc. Heteroatom Chem 17:619–633, 2006; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/hc.20204     

Conformational equilibrium study of calix[4]tetrolarenes using Density Functional Theory (DFT) and Molecular Dynamics simulations

The present work deals with the theoretical study of conformational equilibrium of calix[4]tetrolarenes (1,2,3 trimethoxy substituted calix[4]arene) in gaseous and solvation phases. A total of 64 calculations (four conformations at eight level of theories) were performed using Density Functional Theory (DFT) functionals viz. wB97XD, B97D, B3LYP, CAM-B3LYP with diffused (6–31 + G(d)) and non-diffused basis sets (6-31G(d)). It has been found that the conformational profile of calix[4]tetrolarene changes from cone to 1,3-Alternate upon incorporating – OMe group. The B97D and wB97XD functionals gave the most accurate result having rmsd value ~0.73Å followed by B3LYP and CAM-B3LYP > 0.78Å. Furthermore, NBO calculations demonstrated that reduction in charges at lower rim oxygens reduces the chances of hydrogen bonding. Moreover, global reactivity parameters and molecular dynamics also complement the observed trend. Owing to the methoxy substitutions, anion binding study of these new molecules indicates towards the promising capability to bind Cl^? and F^? ions.     

Theoretical study of a new group of corrosion inhibitors

In the present work, the molecular interactions of four amino acid compounds were simulated through the density functional theory (DFT) indexes to study their inhibitive properties. The prototype inhibitors previously synthesized 2-amino-N-decylacetamide (G), 2-amino-N-decylpropionamide (A), 2-amino-N-decyl-3-methylbutyramide (V), and 2-amino-N-decyl-3-(4-hydroxyphenyl)propionamide (T) were used to test the accuracy of this calculation. The generalized gradient approximation (GGA) was the ab initio approach used to optimize the ground state of the molecules. The simulation of molecular dynamics with force field (AMBER) was calculated to obtain the interaction energy between the metallic surface and the inhibitor molecules. A strong correlation of the global and local indexes with the inhibiting capacity was observed. The inhibitive properties of compounds on mild steel in an aqueous hydrochloric acid solution agreed well with those derived from the reactivity and selectivity indexes in gaseous phase.     

Conformation and reactivity of α-oxo-ketenes: Ab initio and semiempirical (AM1, PM3) calculations

Ab initio MP2/6-31G*//MP2/6-31G* and semiempirical AM1 and PM3 calculations on a series of differently substituted α-oxo-ketenes are used to investigate E/Z-isomerism and rotational barriers in these molecules. Sterically crowded derivatives are found to exist solely as s-E conformers. The unusual stability of these derivatives thus can be attributed to their inability to adopt the s-Z conformation required for the normal α-oxo-ketene reactions. With respect to structures and energies, the PM3 method (especially in the case of highly crowded molecules) is found to be less reliable than AM1. Ab initio HF/3-21G and PM3 vibrational frequencies appear to be of sufficient accuracy for a distinction between s-Z and s-E conformers. In this respect, the AM1 method appears less reliable. © 1994 by John Wiley & Sons, Inc.