A first-principle DFT study of solvent effects on metiamide tautomers and imaginary interactions with H<Subscript>2</Subscript>-receptors

Metiamide is a class of medications called H₂-receptor antagonist or H₂ blockers. It decreases the amount of acid made in the stomach. It is commonly used in the treatment for peptic ulcer disease and gastroesophageal reflux disease. In this study, the metiamide tautomer stability, structural data, HOMO and LUMO (energies and shapes), ΔΕ _HOMO–LUMO gaps, UV–visible data and graphs, dipole moments, Mulliken charges, and thermodynamic and kinetic stabilities in aqueous media as a biological solvent and some of the different media (vacuum, H₂O, Et-OH and DMSO) were investigated for the tautomers of metiamide by the density functional theory DFT-B3LYP/6-31G* method. The results presented that the probability of the adaptability and compatibility of which tautomer (T1–T6) are better than the other tautomers for interactions with the pattern and structural map of the H₂-receptor. The diversities of the interaction points and mosaic patterns of the T3 and T4 tautomers in H₂O media with the imaginary H₂-receptor were investigated.     

Pair-excitation MCSCF treatment of small molecules in an optimized Slater–Transform–Preuss basis set

Pair-excitation multiconfigurational self-consistent field (PEMCSCF ) treatment of 11 small molecules (LiH, BeH₂, BH₃, BF, CH₄, C₂H₄, C₂H₂, CH₂O, NH₃, H₂O, and HF) has been carried out in a minimum basis set of Slater Transform Preuss functions as fitted by six cartesian gaussians (STP -6G). The advantages of accuracy without using a split basis are shown by comparison to familiar 4-31G and 6-31G calculations using molecular geometries optimized with STO -6G basis sets. A benefit is shown for the use of minimum basis fitted to STP functions: they overemphasize long-tail radial dependence to achieve long range basis sensitivity without increasing the basis size at the AO -to-MO transformation step in the configuration interaction portion of the MCSCF algorithm. Fully optimized STP -6G parameters are given and appear to be transferable as shown for acrolein. A FORTRAN listing of the full least squares fitting algorithm is available* for in situ generation of STP -6G orbitals energetically superior to 4-31G, or a less accurate STP -6G 1S, 2S, and 2P basis may be scaled directly as if they were STO -6G functions, but with considerably lower energy that with an STO -6G basis.     

Theoretical Investigation of the Additivity of Structural Substituent Effects in Benzene Derivatives

The additivity of substituent effects in 1,3- and 1,4-disubstituted C₆H₄X₂, and 1,3,5-trisubstituted C₆H₃X₃ (X=F, Cl, CN, NO₂, CH₃, CF₃, NH₂, OH) benzene derivatives on the ring geometry has been investigated. The analysis is based on ab initio calculations at the MP2/6-31G** level of theory. The substituent impacts on the benzene ring are generally in good agreement with the results reported in earlier experimental and lower level theoretical studies. The impacts determined in the monosubstituted benzenes were used to estimate the ring distortions in the di- and trisubstituted derivatives. The estimated ring CC bond distances agree generally within 0.001 Å and the estimated CCC bond angles within 0.3 degree, with the optimized ones. The best agreement (deviations up to only 0.0003 Å and 0.03 deg.) between the estimated and optimized geometrical parameters was obtained for the CH₃ derivatives. Generally, the para-disubstituted derivatives showed the best compliance with additivity, somewhat poorer agreement characterized the meta derivatives while the trisubstituted derivatives showed angular distortions of up to about 0.4°.     

Isomeric tropane analogues of histamine H2-receptor antagonists

In an attempt to investigate stereochemical requirements of antagonists of the histamine H₂-receptor, tropane analogues of cimetidine and metiamide have been synthesized possessing axial and equatorial N-methylthiourea and N-methyl-N-cyanoguanidine moieties at the 3-position of the tropane system. The tropane analogues of this study have been fully characterized with regard to configuration at C-3 and to conformation of the piperidine ring of the tropane nucleus.     

HNNH3, a new possible isomer of N2H4: An ab initio study

Ab initio studies applying the 3-21G, 6-31G, and 6-31G** basis sets and also including the MP2 correction were carried out on H₂NNH₂, HNNH₃, and the transition state of the reaction H₂NNH₂(DOUBLE BOND)HNNH₃. First, the geometries of molecules were optimized using the theoretical methods mentioned in the restricted Hartree–Fock (RHF) scheme. The energies of the molecules corresponding to RHF/6-31G** geometries were subsequently calculated including electron-correlation effects at the level of the second-order Møller–Plesset (MP2) perturbation theory. The vibrational frequencies, net charges, and dipole moments were obtained from the theoretical calculations. The results of our calculations indicate unambiguously that H₂NNH₂ is thermodynamically more stable than is HNNH₃. On the other hand, an isolated HNNH₃ molecule once created would be stable since barriers for its unimolecular isomerization and decomposition are relatively high. But HNNH₃ is unlikely to be isolated in measurable amounts because of bimolecular tautomerization. Nevertheless, HNNH₃ can be considered as an intermediate in chemical processes involving N₂H₄. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 64 : 447–452, 1997     

The 6-31G++ basis set: An economical basis set for correlated wavefunctions

The 6-31G ⁺⁺ basis set is described. This basis set is very similar to the existing 6-31G ^** set but is somewhat smaller through the use of five (rather than six) second-order Gaussians (d functions) and has polarization function exponents optimized for correlated rather than Hartree–Fock wavefunctions. The performance of 6-31G ⁺⁺ is compared with that of the 6-31G ^** and 6-31G ^** basis sets through calculation of the geometries and atomization energies for the set of molecules LiH, FH, H₂O, NH₃, CH₄, N₂, CO, HCN, and HCCH.     

Ab initio GB study of chemical intermediates in solution; Ethylenesulfonium ion in hydrolysis of 2-chloroethyl methyl sulfide

Ab initio MO theory including solvent effects has been applied to the structure and reactivity of methyl ethylenesulfonium ion, 1 , in aqueous solution as a model of the three-membered cyclic sulfonium intermediate expected in the toxic action of sulfur mustard. The 6-31 + G* geometry optimization of the cyclic sulfonium ion 1 suggested that the ring size of 1 is expanded slightly by solvation. The contour lines map of the interaction energy between 1 and Cl⁻ has a very shallow and wide well at 5–6 Å distance from 1 . This is the solvent-separated ion pair, and the contact ion pair was not found between 1 and Cl⁻. The calculated energy diagrams for the S_N2-type reactions of 1 with Cl⁻, H₂O, and OH⁻ that give ring-opened compounds indicated the following: (1) The energy of the 1 + Cl⁻ system is similar to that of chloroethyl methyl sulfide (CEMS, 2 ), and the interconversion between 1 + Cl⁻ and 2 occurs easily in aqueous solution. The 3-21 + G(*) and 6-31 + G* activation energies for the 2 → 1 + Cl⁻ reaction, 20–22 kcal/mol, agree well with the experimental enthalpy of activation for the hydrolysis of 2 . (2) The reaction of 1 with OH⁻ gives a very stable hydroxyl compound, 4 , and no transition state was found. (3) The reaction of 1 with H₂O gives an unstable addition product that is expected to be converted to 4 with the assistance of another H₂O molecule. This mechanism is consistent with that proposed by Bartlett and Swain in their pioneering work on the hydrolysis of sulfur mustard. © 1998 John Wiley & Sons, Inc. Heteroatom Chem 9:503–510, 1998     

苯分子与Si6O18H12和Al6O24H30团簇模型相互作用的理论研究

粘土矿已经被广泛用来去除有机物,修复和净化被石油碳氢化合物污染的土壤和地下水. 我们选择高岭石作为研究对象,构造了Si₆O₁₈H₁₂和Al₆O₂₄H₃₀两个团簇模型分别代表高岭石的硅氧层表面和铝氧层表面,在MP2/6-31G(d,p)//B3LYP/6-31G(d,p)的理论水平上系统地研究了气态下苯分子和高岭石团簇模型的相互作用. 并进一步分析了苯分子和高岭石表面相互作用的各种气态性质,比如：优化的几何构型、结构参数、吸附能、自然键轨道电荷分布、振动频率变化、静电势、电子密度性质(次级氢键的电子密度和拉普拉斯算符值)和电子密度差分等. 优化的几何构型表明苯分子吸附在高岭石表面的本质可能是次级氢键的形成. 其他性质的结果进一步验证了次级氢键的存在,并指出苯更倾向于吸附在高岭石的铝氧层表面,且苯环和铝氧层表面形成近似90°的夹角.     

A non-empirical molecular orbital study of valence tautomers of C2H3N

The valence tautomers of C₂H₃N have been examined by non-empirical molecular orbital calculations using two split-valence shell basis sets. All geometries were fully optimized using the 4–31G basis set and these structures were then used in 6–31G basis set calculations. The order of stability of the three possible cyclic isomers is 1-azirine > cyclic carbene > 2-azirine. The profiles for conversion of vinylmethylene into cyclopropene, vinylnitrene into 1-arizine, and iminomethylene into 2-azirine have all been shown to have barriers.     

Eclipsed and staggered conformations of (SIH3)2F+: An ab initio study

Ab initio calculations at 6–31G**, 6–31++G**, and MP2/6–31G** levels were performed on disilyl–fluoronium, (SiH₃)₂F⁺, with the SiH₃ group eclipsed or staggered. Optimized geometries, total energies, dipole moments, atomic charges, electronic density, and vibrational frequencies were computed. The results were compared with calculated structural parameters and vibrational frequencies of H₃SiF, H₂SiF⁺, H₂SiF^?, and H₄SiF⁺ ions. The basis-set effects were studied. Several thermochemistry parameters—ZPE, thermal energy, rotational constants, and entropies—were also calculated. © 1994 John Wiley & Sons, Inc.     

DFT study on the hydrogen bonds of phenol–cyclohexanone and phenol–H2O2 in the Baeyer–Villiger oxidation

Hydrogen bonds of phenol–cyclohexanone and phenol–H₂O₂ in the studied Baeyer–Villiger (B–V) oxidation have been investigated by HF, B3LYP, and MP2 methods with various basis sets. The accurate single‐point energies were performed using CCSD(T)/6‐31+G(d,p) and CCSD(T)/aug‐cc‐pVDZ on the optimized geometries of MP2/6‐31+G(d,p). It has been confirmed that B3LYP/6‐31+G(d,p) could be used to study such hydrogen bonds. Energetic analysis of complexes was carried out using the Xantheas method with BSSE corrected by CP method. Orbital energy order (ε) illuminated that phenol with good hydrogen donor‐acceptor property can interact with cyclohexanone or H₂O₂ to form hydrogen bound complexes, and the binding energies (BE) range from ?4.38 to ?14.06 kcal mol^?1. NBO analysis indicated that the redistribution of atomic charges in the complexes facilitated nucleophilic attack of H₂O₂ on cyclohexanone. The calculated results match remarkably well with the experimental phenomena. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Proton and sodium ion affinities of glycine and its sodium salt in the gas phase. Ab initio calculations

The energies of protonation and Na⁺ cationization of glycine (GLY) and its (GLY ? H + Na) salt in the gas phase were calculated using ab inltio calculations. The proton affinity of GLY, valued at the MP2/6–31G*//3-21G level, is 937 kJ mol^?1. The amino function is confirmed to be the most favourable site of protonation: ‘proton affinities’ of the carbonyl and hydroxyl functions are calculated to be 75 and 180 kJ mol^?1, respectively, lower than that of NH₂ at the MP2/6-31G*//3–21G level. Calculations performed up to the MP2/6–31G*//3–21G level give the Na⁺ affinity of GLY as 189 kJ mol^?1 and the H⁺ and Na⁺ affinities of (GLY – H + Na) as 1079 and 298 kJ mol^?1, respectively. The geometries of all neutral and protonated species optimized with the 3–21G basis set are described. Both H* and Na⁺ cations complex preferably between the nitrogen atom and the carbonyl oxygen atom, leading to pseudo-five-membered ring structures in which Na? O and Na? N bonds lengths are greater than 2 Å.     

Rearrangement of azirine intermediates to nitriles: Theoretical study of cleavage of 3,4-dihydro-1aH-azirine[2,3-c]pyrrol-2-one to cyanoketene–formaldimine complex

A detailed investigation of the reaction path for the thermal rearrangement of 3,4-dihydro-1aH-azirine[2,3-c]pyrrol-2-one to yield a cyanoketene–formaldimine complex is carried out at the MP2/6-31G* and B3LYP/6-31G* levels of theory. The ring opening of the five-membered pyrrolinone ring and the formation of the nitrile group takes place in a concerted manner, presenting a significant strain energy release and allowing for an electronic stabilization by coarctate conjugation of the transition structure (TS). These two factors make possible a moderate energy barrier. Although the structural features B3LYP/6-31G* theoretical levels, it is found that the MP2 energy barrier (28.8) CCSD(T)/6-31G*//MP2/6-31G* value (17.1 kcal/mol). The complex electronic rearrangement can be rationalized using the theory of coarctate transition structures developed by Herges as the evolution of an azirine structure without referring to a hypothetical vinyl nitrene intermediate. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 912–922, 1998     

Bowl-shaped hydrocarbons related to C60

Ab initio studies at the HF/6-31G* and B3LYP/6-31G* levels are reported for two bowl-shaped hydrocarbons related to C₆₀: C₃₀H₁₂ and C₃₆H₁₂, of C₃ and C_3v symmetry, respectively. The former has an approximate heat of formation of 211 kcal/mol. Bowl-to-bowl interconversion may occur through a planar (C_3h) form of ca. 64 kcal/mol greater energy having one imaginary vibrational frequency. The larger C₃₆H₁₂ bowl has a calculated ΔH°_f of 265 kcal/mol. Its HF/6-31G*, B3LYP/6-31G*, and MM3 bond lengths are in good agreement with a recent X-ray structure. Chemical shifts for both compounds calculated by the GIAO method are in good agreement with the measured NMR spectra. The observed ¹³C chemical shifts increase with the extent of pyramidalization. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 189–194, 1998     

A computational mechanistic study of the deamination reaction of melamine

A detailed computational study of the deamination reaction of melamine by OH^–, n H₂O/OH^–, n H₂O (where n = 1, 2, 3), and protonated melamine with H₂O, has been carried out using density functional theory and ab initio calculations. All structures were optimized at M06/6‐31G(d) level of theory, as well as with the B3LYP functional with each of the basis sets: 6‐31G(d), 6‐31 + G(d), 6‐31G(2df,p), and 6‐311++G(3df,3pd). B3LYP, M06, and ω B97XD calculations with 6‐31 + G(d,p) have also been performed. All structures were optimized at B3LYP/6‐31 + G(d,p) level of theory for deamination simulations in an aqueous medium, using both the polarizable continuum solvation model and the solvation model based on solute electron density. Composite method calculations have been conducted at G4MP2 and CBS‐QB3. Fifteen different mechanistic pathways were explored. Most pathways consisted of two key steps: formation of a tetrahedral intermediate and in the final step, an intermediate that dissociates to products via a 1,3‐proton shift. The lowest overall activation energy, 111 kJ mol^?1 at G4MP2, was obtained for the deamination of melamine with 3H₂O/OH^?.     

Theoretical Investigations on Fluorine-Substituted Ethylene Dications C2HnF4-n 2+(n = 0–4)

The optimized geometries and energies of fluorine-substituted ethylene dications C₂H_nF_4-n ²⁺ (n = 0–4) have been investigated by means of ab initio methods. At the MP3/6-31G**//6-31G* + zero-point energy level of theory, the results predict that C₂F₄²⁺ and C₂HF₃²⁺ are planar, while C₂H₄²⁺, C₂H₃F²⁺ and 1,1—C₂H₂F₂²⁺ prefer a perpendicular geometry. For 1,2—C₂H₂F₂²⁺ an energy difference of only 0.3 kcal/mol is found between the (trans) planar and perpendicular structure. The stabilizations attributed to hyperconjugation, fluorine lone-pair donation, and (C? F) double-bond conjugation are discussed. A comparison is made for the C? C and C? F stretching frequencies determined at 6-31G*//6-31G* between the neutral and dicationic species. The theoretically determined ionization energies for the vertical process N⁺ → N²⁺ at the MP3/6-31G*//3-21G level are compared with experimental Q_min values.     

The molecular,electronic structures and IR and Raman spectra of metal-free,N,N′-dideuterio,and magnesium tetra-2,3-pyrazino-porphyrazines: Density functional calculations

A theoretical investigation of the fully optimized geometries and electronic structures of the metal-free (TPyPzH₂), N,N′-dideuterio (TPyPzD₂), and magnesium (TPyPzMg) tetra-2,3-pyrazino-porphyrazine has been conducted based on density functional theory. The optimized geometries at density functional theory level for these compounds are reported here for the first time. A comparison among the different molecules, including tetra-2,3-pyridino-porphyrazine (TPdPzH₂), phthalocyanine (H₂Pc) compounds, for the geometry, molecular orbital, and atomic charge is made. The substituent effect of the N atoms on these properties of these compounds is discussed.The IR and Raman frequencies and intensities for the three compounds have also been calculated at density functional B3LYP level using the 6-31G(d) basis set. Detailed assignments of the N–H, N–M, and pyrazine ring vibrational bands in the IR and Raman spectra have been made based on assistance of animated pictures. The isotope effect of D atoms are also discussed.     

Theoretical investigation of unimolecular decomposition channels of furan4

Density functional theory and high-level ab initio calculations were carried out to investigate three unimolecular decomposition channels of furan. All equilibrium and transition state structures along the proposed decomposition channels are fully optimized by B3LYP/6-31G** and characterized at the same level of theory by vibrational and intrinsic reaction coordinate analyses. Relative energies of the optimized structures were evaluated at theoretical levels up to QCISD(T)/6-311++G**. The theoretical results suggest that the unimolecular decomposition channel of isoxazole, proposed in an experimental study and implied to be the main decomposition channel of furan, is responsible only for the formation of HC(TRIPLE BOND)CH and H₂O(DOUBLE BOND)C(DOUBLE BOND)O, minor products of furan thermal decomposition. A new decomposition mechanism, proposed in the present study, is shown to be more likely responsible for the formation of CH₃C(TRIPLE BOND)CH and CO, major products of furan thermal decomposition. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 240–249, 1998     

The effect of including polarization functions on the geometrical parameters calculated for benzene,fluorobenzene and cyanobenzene

The geometrical parameters for benzene, fluorobenzene, and cyanobenzene have been calculated using the 6–31G*(5D) and 6–31G** basis sets, and, in addition, the 6–31 + G*(5D) basis set in the case of fluorobenzene. Compared to previous results obtained using the 6-31G basis set there are minor changes in the magnitude of the bond lengths and angles in the ring, but the relative values remain unaltered. The values for the ipso angles in fluorobenzene and cyanobenzene are again somewhat less than those reported from microwave and/or electron diffraction studies. The distortion of the ring is characterized as either an elongation or a flattening with respect to the F–C₁ ⃛C₄ and NC C₁C₄ axes, and the shape is characterized as either a broadening or a narrowing across the ring just below the F and Cn group, i.e., an increase or a decrease in the C₂ ⃛C₆ internuclear distance, relative of the C₃ ⃛C₅ distance.     

FT-IR and FT-Raman vibrational spectra and molecular structure investigation of nicotinamide: A combined experimental and theoretical study

In this work, the experimental and theoretical spectra of nicotinamide (C₆H₆N₂O) are studied. FT-IR and FT-Raman spectra of title molecule in the liquid phase have been recorded in the region 4000–100 cm^?1. The structural and spectroscopic data of the molecule in the ground state have been calculated by using Hartree–Fock and density functional method (B3LYP) with the 6-31+G*(d, p) and 6-31++G* (d, p)basis set. The vibrational frequencies have been calculated and scaled values have been compared with the experimental FT-IR and FT-Raman spectra. The observed and calculated frequencies are found in good agreement. The DFT-B3LYP/6-31++G (d, p) calculations have been found are more reliable than the ab initio HF/6-31+G (d, p) calculations for the vibrational study of nicotinamide. The optimized geometric parameters (bond lengths and bond angles) are compared with experimental values of the molecule. The alteration of vibrational bands due to the substitutions in the base molecule is also investigated from their characteristic region of linked spectrum.