Quantum Theory Of Atoms In Molecules Calculations Research Articles

4,4-Bis(isopropylthio)-1,1-diphenyl-2-azabuta-1,3-diene Adducts with Cadmium(II), Mercury(II) and Copper(I) Iodides: Crystal, Molecular and Electronic Structures of d10 Transition Metal Chelate Complexes

The thioether-functionalized 2-azabutadiene (iPrS)2C=C(H)-N=CPh2 L ligates to CdI2 and HgI2 to form the chelate compounds [CdI2{(iPrS)2C=C(H)-N=CPh2] (1) and [HgI2(iPrS)2C=C(H)-N=CPh2] (2). Their crystal structures were solved via X-ray diffraction. Both crystallize in the non-centrosymmetric space groups: monoclinic P21 (1) and orthorhombic P212121 (2), respectively. The closed-shell d10 metal centers are four-coordinated (two iodides and S and N coordinating atoms from the ligand L) in both complexes. The geometrical indexes τ indicate that a highly distorted trigonal pyramidal is adopted for 1 and a seesaw geometry for 2. The comparative nature of metal–ligand bonds is discussed on the basis of metric parameters and of QT-AIM (quantum theory of atoms in molecules) calculations. L was also treated with CuI to obtain the dinuclear species [LCu(μ2-I2)CuL] (3), in which the two Cu(I) centers are linked by a short metal–metal bond. The geometric and electronic properties of 3 are compared with those of 1 and 2.

Evaluating the conformational space of the active site of D2 dopamine receptor. Scope and limitations of the standard docking methods.

We report here for the first time the potential energy surfaces (PES) of phenyletilamine (PEA) and meta-tyramine (m-OH-PEA) at the D2 dopamine receptor (D2DR) binding site. PESs not only allow us to observe all the critical points of the surface (minimums, maximums, and transition states), but also to note the ease or difficulty that each local minima have for their conformational inter-conversions and therefore know the conformational flexibility that these ligands have in their active sites. Taking advantage of possessing this valuable information, we analyze how accurate a standard docking study is in these cases. Our results indicate that although we have to be careful in how to carry out this type of study and to consider performing some extra-simulations, docking calculations can be satisfactory. In order to analyze in detail the different molecular interactions that are stabilizing the different ligand-receptor (L-R) complexes, we carried out quantum theory of atoms in molecules (QTAIM) computations and NMR shielding calculations. Although some of these techniques are a bit tedious and require more computational time, our results demonstrate the importance of performing computational simulations using different types of combined techniques (docking/MD/hybrid QM-MM/QTAIM and NMR shielding calculations) in order to obtain more accurate results. Our results allow us to understand in details the molecular interactions stabilizing and destabilizing the different L-R complexes reported here. Thus, the different activities observed for dopamine (DA), m-OH-PEA, and PEA can be clearly explained at molecular level.

Halogen Interactions in Halogenated Oxindoles: Crystallographic and Computational Investigations of Intermolecular Interactions.

X-ray structural determinations and computational studies were used to investigate halogen interactions in two halogenated oxindoles. Comparative analyses of the interaction energy and the interaction properties were carried out for Br···Br, C-H···Br, C-H···O and N-H···O interactions. Employing Møller–Plesset second-order perturbation theory (MP2) and density functional theory (DFT), the basis set superposition error (BSSE) corrected interaction energy (Eint(BSSE)) was determined using a supramolecular approach. The Eint(BSSE) results were compared with interaction energies obtained by Quantum Theory of Atoms in Molecules (QTAIM)-based methods. Reduced Density Gradient (RDG), QTAIM and Natural bond orbital (NBO) calculations provided insight into possible pathways for the intermolecular interactions examined. Comparative analysis employing the electron density at the bond critical points (BCP) and molecular electrostatic potential (MEP) showed that the interaction energies and the relative orientations of the monomers in the dimers may in part be understood in light of charge redistribution in these two compounds.

Probing the intermolecular interactions, binding affinity, charge density distribution and dynamics of silibinin in dual targets AChE and BACE1: QTAIM and molecular dynamics perspective

Alzheimer’s disease (AD) is the grievous neurodegenerative disorder. Reportedly, many enzymes are responsible for this disease, in which notably, acetylcholinesterase (AChE) and β-secretase (BACE1) are largely involved for AD. An experimental study reports that silibinin molecule inhibits both AChE and BACE1 enzymes. Present study aims to understand the dual binding mechanism of silibinin in the active site of AChE and BACE1 from the intermolecular interactions, conformational flexibility, charge density distribution, binding energy and the stability of molecule. To obtain the above information, the molecular docking, molecular dynamics (MD) and QTAIM (quantum theory of atoms in molecules) calculations have been performed. The molecular docking reveals that silibinin molecule is forming strong and weak intermolecular interactions with the catalytic site of both enzymes. The QTAIM analysis for the binding pockets of both complexes shows the charge density distribution of intermolecular interactions. The electrostatic potential map displays the electronegative/positive regions at the interaction zone of silibinin with AChE and BACE1 complexes. The MD simulation confirms that the silibinin molecule is stable in the active site of AChE and BACE1 enzymes. The binding free energies of silibinin with both enzymes are more favorable to have the interactions. Communicated by Ramaswamy H. Sarma

Novel Sulfonamide-Based Carbamates as Selective Inhibitors of BChE.

A series of 14 target benzyl [2-(arylsulfamoyl)-1-substituted-ethyl]carbamates was prepared by multi-step synthesis and characterized. All the final compounds were tested for their ability to inhibit acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) in vitro, and the selectivity index (SI) was determined. Except for three compounds, all compounds showed strong preferential inhibition of BChE, and nine compounds were even more active than the clinically used rivastigmine. Benzyl {(2S)-1-[(2-methoxybenzyl)sulfamoyl]-4-methylpentan-2-yl}carbamate (5k), benzyl {(2S)-1-[(4-chlorobenzyl)sulfamoyl]-4-methylpentan-2-yl}carbamate (5j), and benzyl [(2S)-1-(benzylsulfamoyl)-4-methylpentan-2-yl]carbamate (5c) showed the highest BChE inhibition (IC50 = 4.33, 6.57, and 8.52 µM, respectively), indicating that derivatives 5c and 5j had approximately 5-fold higher inhibitory activity against BChE than rivastigmine, and 5k was even 9-fold more effective than rivastigmine. In addition, the selectivity index of 5c and 5j was approx. 10 and that of 5k was even 34. The process of carbamylation and reactivation of BChE was studied for the most active derivatives 5k, 5j. The detailed information about the mode of binding of these compounds to the active site of both BChE and AChE was obtained in a molecular modeling study. In this study, combined techniques (docking, molecular dynamic simulations, and QTAIM (quantum theory of atoms in molecules) calculations) were employed.

Combined MD/QTAIM techniques to evaluate ligand-receptor interactions. Scope and limitations.

In medicinal chemistry, it is extremely important to evaluate, as accurately as possible, the molecular interactions involved in the formation of different ligand-receptor (L-R) complexes. Evaluating the different molecular interactions by quantum mechanics calculations is not a simple task, since formation of an L-R complex is a dynamic process. In this case, the use of combined techniques of molecular dynamics (MD) and quantum calculations is one the best possible approaches. In this work we report a comparative study using combined MD and QTAIM (Quantum Theory of Atoms In Molecules) calculations for five biological systems with different levels of structural complexity. We have studied Acetylcholinesterase (AChE), D2 Dopamine Receptor (D2DR), beta Secretase (BACE1), Dihydrofolate Reductase (DHFR) and Sphingosine Kinase 1 (SphK1). In these molecular targets, we have analyzed different ligands with diverse structural characteristics. The inhibitory activities of most of them have been previously measured in our laboratory. Our results indicate that QTAIM calculations can be extremely useful for in silico studies. It is possible to obtain very accurate information about the strength of the molecular interactions that stabilize the formation of the different L-R complexes. Better correlations can be obtained between theoretical and experimental data by using QTAIM calculations, allowing us to discriminate among ligands with similar affinities. QTAIM analysis gives fairly accurate information for weak interactions which are not well described by MD simulations. QTAIM study also allowed us to evaluate and determine which parts of the ligand need to be modified in order to increase its interactions with the molecular target. In this study we have discussed the importance of combined MD/QTAIM calculations for this type of simulations, showing their scopes and limitations.

Stability Investigation of Se- and Te-Substituted Tetrahedrite

Density functional theory (DFT) and quantum theory of atoms in molecules (QTAIM) calculations have been performed to investigate the influence of Se and Te substitutions on the stability of tetrahedrite. The formation energies, the Laplacian of the electron density, the total and local volume changes and the site charge changes have been determined. The DFT results show that the single-atom substitution is the most preferred, and that selenium is preferred for sulphur replacement, whereas antimony can be equally replaced by selenium and tellurium. However, as selenium substitutes for either sulphur or antimony, QTAIM calculations evidenced large structural and chemical changes in the structures. Additional experimental results, in agreement with QTAIM results, allow us to conclude that the investigation of alloy stability cannot be limited to DFT calculations of formation energies and total volume changes. QTAIM appears to be a valuable complementary method for obtaining information on phase stability.

X-ray diffraction and QTAIM calculations of the non-covalent intermolecular fluorine-fluorine interactions in tris(trifluoroacetylacetonato)-manganese(III)

An X-ray study was conducted on tris(trifluoroacetylacetonato)-manganese(III). The two different molecules in the unit cell exhibit a distorted octahedral geometry and moderate compression Jahn–Teller distortion respectively. A number of F⋯H and F⋯F intermolecular interactions were observed in the packing of the crystals in the solid state. The crystal packing effects could lead to the observed compressed Jahn–Teller distortion, instead of the expected elongation Jahn–Teller distortion. A Quantum Theory of Atoms in Molecules (QTAIM) study on selected molecular pairs identified intermolecular bond paths (BP) between fluorine atoms with bond critical points (BCP) which have positive electron density (ρ) and Laplacian of electron density (∇2ρ) values, which are indicative of weak F···F non-covalent bonding interactions.

Li and Na Adsorption on Graphene and Graphene Oxide Examined by Density Functional Theory, Quantum Theory of Atoms in Molecules, and Electron Localization Function.

Adsorption of Li and Na on pristine and defective graphene and graphene oxide (GO) is studied using density functional theory (DFT) structural and electronic calculations, quantum theory of atoms in molecules (QTAIM), and electron localization function (ELF) analyses. DFT calculations show that Li and Na adsorptions on pristine graphene are not stable at all metal coverages examined here. However, the presence of defects on graphene support stabilizes both Li and Na adsorptions. Increased Li and Na coverages cause metal nucleation and weaken adsorption. Defective graphene is associated with the presence of band gaps and, thus, Li and Na adsorptions can be used to tune these gaps. Electronic calculations show that Li– and Na–graphene interactions are Coulombic: as Li and Na coverages increase, the metal valences partially hybridize with the graphene bands and weaken metal–graphene support interactions. However, for Li adsorption on single vacancy graphene, QTAIM, ELF, and overlap populations calculations show that the Li-C bond has some covalent character. The Li and Na adsorptions on GO are significantly stronger than on graphene and strengthen upon increased coverages. This is due to Li and Na forming bonds with both carbon and oxygen GO atoms. QTAIM and ELF are used to analyze the metal–C and metal–metal bonds (when metal nucleation is present). The Li and Na clusters may contain both covalent and metallic intra metal–metal bonds: This effect is related to the adsorption support selection. ELF bifurcation diagrams show individual metal–C and metal–metal interactions, as Li and Na are adsorbed on graphene and GO, at the metal coverages examined here.

253 - Inhibition of toll-like receptor (TLR-2) activity by the nitrone spin trap 5,5-dimethyl-1-pyrroline N-oxide (DMPO) through direct interaction with the toll/interleukin-1 receptor domain (TIR)

Inflammatory activation of macrophages throughout toll-like Receptors (TLRs) is a fundamental step in the development of immune response triggered by bacterial and fungal infections. TLR deregulated signaling has been implicated in a great number of inflammatory diseases, thus their inhibition has become a therapeutic target. Our group has studied the effects of the nitrone spin trap 5,5-dimethyl-1-pirroline N-oxide (DMPO) on lipopolysaccharide (LPS)-primed RAW264.7 cells at transcriptomic and functional levels. Using microarray technology we found that the dampening effects of DMPO could be explained by impaired signaling of TLR-2 and -4. Because TLRs have structural similarities in their intracellular domain responsible for signal transduction called TIR, we hypothesize that the effects of DMPO in these receptors are caused by a direct binding of the spin trap to their TIR domain. Currently, TLR2-TIR domain is the only one crystallized, therefore we use combined techniques of docking, molecular dynamics simulations and QTAIM (Quantum Theory of Atoms In Molecules) calculations to determine the interaction between DMPO and this receptor. Our results show that DMPO could bind to four specific residues in a key region pointed as responsible for signal transduction known as BB-loop. To functionally corroborate these results we used an experimental model based on hTLR2.6-expressing HEKs cells and observed that DMPO can block zymozan-triggered-TLR2-mediated NF-kb activation. Because TLRs bind to adaptor protein MyD88 we used co-immunoprecipitation to test whether DMPO can prevent TLR2-MyD88 binding and found no effect. Taking together our results show that DMPO blocks TLR2 signaling without preventing completely the coupling of MyD88 to its TLR2-TIR domain. This phenomena could be due to DMPO disrupting proper coupling of TLRs with MyD88 by direct binding to the BB-loop region responsible for signal transduction. These data encourage DMPO derivatives as mechanism-based TLR inhibitors. Supported by PICT3369, PIP916 and PROICO-2018 to DCR and SEGM

Theoretical models to predict the inhibitory effect of ligands of sphingosine kinase 1 using QTAIM calculations and hydrogen bond dynamic propensity analysis.

We report here the results of two theoretical models to predict the inhibitory effect of inhibitors of sphingosine kinase 1 that stand on different computational basis. The active site of SphK1 is a complex system and the ligands under the study possess a significant conformational flexibility; therefore for our study we performed extended simulations and proper clusterization process. The two theoretical approaches used here, hydrogen bond dynamics propensity analysis and Quantum Theory of Atoms in Molecules (QTAIM) calculations, exhibit excellent correlations with the experimental data. In the case of the hydrogen bond dynamics propensity analysis, it is remarkable that a rather simple methodology with low computational requirements yields results in excellent accord with experimental data. In turn QTAIM calculations are much more computational demanding and are also more complex and tedious for data analysis than the hydrogen bond dynamic propensity analysis. However, this greater computational effort is justified because the QTAIM study, in addition to giving an excellent correlation with the experimental data, also gives us valuable information about which parts or functional groups of the different ligands are those that should be replaced in order to improve the interactions and thereby to increase the affinity for SphK1. Our results indicate that both approaches can be very useful in order to predict the inhibiting effect of new compounds before they are synthesized.

Estimating the energy of intramolecular bifurcated (three-centered) hydrogen bond by X-ray, IR and 1H NMR spectroscopy, and QTAIM calculations

The structure of a series of the 2-substituted pyrroles containing the intramolecular NH⋯OC hydrogen bond and the 2,5-disubstituted pyrroles having the three-centered intramolecular CO⋯H(N)⋯OC hydrogen bond was studied by X-ray, IR, NMR spectroscopy, and quantum chemical calculations, including calculations in terms of Quantum Theory of Atoms in Molecules (QTAIM). The X-ray diffraction data show that the intramolecular rH … O and RN … O distances increase upon conversion from a two-centered hydrogen bond to a three-centered hydrogen bond, indicating a weakening of the last interaction. The topological parameters of the hydrogen bond obtained from the QTAIM calculations (electron density ρBCP and potential energy density V at the critical point of the hydrogen bond) made it possible to quantify the hydrogen bond energy in the studied compounds. The energy of intramolecular NH⋯OC hydrogen bond in the 2-substituted pyrroles is in the range from 5.5 to 6.0 kcal mol−1 (moderate hydrogen bond), whereas the energy of the components of the three-centered intramolecular CO⋯H(N)⋯OC hydrogen bond in the 2,5-disubstituted pyrroles is reduced to 2.5–3.5 kcal mol−1 (weak hydrogen bond). For weak components of the three-centered hydrogen bond, the principle of the additivity of effects on the spectral parameters influenced by the hydrogen bond (change in the vibration frequency of the NH bond, ΔνN‒H, and the chemical shift of bridge proton, ΔδNH) was proposed. For this reason, the ½ΔνN-H and ½ΔδNH values also were used for estimating the energy of the two-centered components of the three-centered CO⋯H(N)⋯OC hydrogen bond. These estimates are in good agreement with those obtained from the ρBCP and V parameters. Consequently, the ½ΔνN-H and ½ΔδNH parameters can be used as an energy descriptor for two-centered components in the case of a symmetric intramolecular three-centered hydrogen bond.

A heuristic approach to evaluate peri interactions versus intermolecular interactions in an overcrowded naphthalene.

Octachloronaphthalene (OCN), a serious environmental pollutant, has been investigated by charge density analysis to unravel several unexplored factors responsible for steric overcrowding. The topological features of the enigmatic peri interactions contributing to steric overcrowding are qualified and quantified from experimental and theoretical charge-density studies. A new facet in the fundamental understanding of peri interactions is revealed by NCI (non-covalent interaction) analysis. The potential role of these interactions in deforming the molecular geometry and subsequent effect on aromaticity are substantiated from NICS (Nuclear Independent Chemical Shift) and QTAIM (Quantum Theory of Atoms in Molecules) calculations. The eye-catching dissimilarity in the out-of-plane twisting of OCN renders the molecule in an asymmetric geometry in the crystalline phase compared with symmetric geometry in the optimized solvated phase. This is uniquely characterized by their molecular electrostatic potential (MESP), respectively, and is explained in terms of conflict between two opposing forces - peri interactions, and symbiotic intermolecular Cl⋯Cl and Cl⋯π contacts.

Quantum chemical investigation on complexation of palladium with iminopyridyl ligands: Structural, thermochemical, and electronic aspects

ABSTRACTThe main purpose of this research is to investigate computationally the structural, thermochemical and electronic properties in complexation process of dichloride {N-[(5-methylthiophen-2-yl)methylidene]−2-(pyridine-2-yl)ethanamine-κ2N,N′}palladium(II) complex. In the first step, we have concentrated on comparative survey of ability of density functional theory (DFT) and also semi-empirical approaches to reproduce the crystal structure of palladium(II) complex. Comparison of our calculated structural parameters of aforementioned complex with the available crystallographical data reveals that both functionals (B3LYP and M06) can well-reproduce x-ray structure of the complex with a near accuracy while PM6-D2 semi-empirical calculated values are not in a reliable agreement with the crystallographical data.In the next step, we have shown the thermodynamical superiority in using THF as a polar solventin complexation reaction via polarized continuum model (PCM) computations which is in confirmation with experimental observations. Additionally, the bond orders of some selected key bonds in C13H14N2Sligand andPdCl2(C13H14N2S) complex have been evaluated comparatively to analyze the electronic behavior of coordination.Finally, we focused on topological analysis of electron density function via quantum theory of atoms in molecules (QTAIM) approach to explore the strength and nature of metal-ligand interactions on bond and ring critical points (BCPs).Strictly speaking, QTAIM calculations have been performed to determine the electronic density, its Laplacian and other electronic energy density indicators on some key BCPs to interpret the electronic features of complexation.

E versus Z isomers of Fischer aminocarbene complex [Mo(CO)4(PPh3){C(NHCy)(2-furyl)}]: N[sbnd]H⋯O versus C[sbnd]H⋯O intramolecular hydrogen bonds

A density functional theory (DFT) calculated conformation analysis of the twelve possible conformations of [Mo(CO)4(PPh3){C(NHCy)(2-furyl)}], 1, utilizing different DFT methods, showed that both the cis-syn,Z (2), 1-Z, and the cis-syn,Z (2), 1-E, conformations of this aminocarbene complex have similar electronic and Gibbs free energies. The solid-state crystal structure of the 1-Z is presented in this study and compared to the structure of the 1-E conformation. The stability of 1-Z, versus the previously published solid-state structure of conformer 1-E, is ascribed to stabilization of the NH⋯OFu intramolecular bond. This is evidenced by a shorter H⋯O intramolecular bond in the experimental crystal structure of 1-Z. The stability of the 1-Z and 1-E conformations were analyzed and compared by density functional theory, quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) methods. The NBO analysis reveals similar but stronger donor–acceptor interactions in 1-Z than in 1-E. QTAIM calculations indicate 6 and 5 bond critical points related to intramolecular bonds stabilizing the orientation of 2-furyl and NHCy in 1-Z and 1-E respectively.

The electronic donation and frequency shifts on the YCCH⋯BH4− boron-bonded complexes (Y = H, CH3, CF3 and CCl3)

In this theoretical work, the tetrahydroborate ion (BH4−) was used as proton acceptor in the formation of the YCC–H⋯BH4− complexes (Y=H, CH3, CCl3 and CF3). Using B3LYP/6-311++G(d,p) level of theory, the results of structure corroborate with the analyses of infrared spectra showing that the changes in the bond lengths are in good agreement with the frequency shifts of the HCC–H, H3CCC–H, Cl3CCC–H and F3CCC–H proton donors. Based on the calculations carried out by the Quantum Theory of Atoms in Molecules (QTAIM), the reductions of electronic density corroborate with the red shifts in the frequencies of the C–H bonds. In addition to that, the C–H bonds are polarized because the contributions of s orbital diminish whereas of p increase. In line with this, the variations on the atomic radii computed via QTAIM calculations show that carbon outweigh hydrogen as follows (ΔrC>ΔrH). This scenario is indirectly supported by the Bent’s rule of the chemical bonding. Although the interaction energies (corrected with BSSE and ZPE) vary between −19 and −67kJmol−1, these complexes interact without covalent character.

Single crystal architecture and absorption spectra of octathio[8]circulene and sym-tetraselenatetrathio[8]circulene: QTAIM and TD-DFT approach

The single crystal architecture of the high-symmetry octathio[8]circulene and sym-tetraselenatetrathio[8]circulene is studied at the density functional theory (DFT) level with the quantum theory of atoms in molecules (QTAIMs) approach to the electron density distribution analysis. The presence of stabilizing intermolecular C---C, C---S and C---Se contacts in the longitudinal and transversal projections of the single crystals is postulated on the grounds of the previous high-resolution X-ray data for octathio[8]circulene; it is supported by the present QTAIM calculations and also predicted in some new details for both circulenes. We suggest that the appearance of the observed red color for the monocrystalline octathio[8]circulene is caused by strong intermolecular interactions between the molecules in the single crystal. However, the intermolecular interactions for the sym-tetraselenatetrathio[8]circulene crystal fragment are weaker and molecular layers are more friable in comparison to octathio[8]circulene crystal structure. These lead to the absence of visible absorption for the sym-tetraselenatetrathio[8]circulene crystal.

Structure, energy, vibrational spectrum, and Bader's analysis of π⋯H hydrogen bonds and H−δ⋯H+δdihydrogen bonds

In this paper, the intermolecular structural study asserted by the vibrational analysis in the stretch frequencies of hydrogen bonds (π···H) and dihydrogen bonds (H(-δ)···H(+δ)) have definitively been revisited by means of calculations carried out by Density Functional Theory (DFT) and topological parameters derived from the classic treatise of the Quantum Theory of Atoms in Molecules (QTAIM). As a matter of fact the π···H hydrogen bond is formed between the hydrofluoric acid and the C≡C bond of the acetylene, but the QTAIM calculations revealed a distortion in this interaction due to the formation of the ternary complex C(2)H(2)···2(HF). Although the π bonds of ethylene (C(2)H(4)), propylene (C(2)H(3)(CH(3))), and t-butylene (C(2)H(2)(CH(3))(2)) are considered proton acceptors, two hydrogen-bond types--π···H and C···H--can be observed. Over and above the analysis of the π hydrogen bonds, theoretical arguments also were used to discuss the red-shifts in the stretch frequencies of the binary dihydrogen complexes formed by BeH(2)···HX with X = F, Cl, CN, and CCH. Although a vibrational blue-shift in the stretch frequency of the H-C bond of HCF(3) due to the formation of the BeH(2)···HCF(3) dihydrogen complex was obtained, unmistakable red-shifts were detected in LiH···HCF(3), MgH(2)···HCF(3), and NaH···HCF(3). Moreover, the alkali-halogen bonds were identified in relation to the formation of the trimolecular systems NaH···2(HCF(3)) and NaH···2(HCCl(3)). At last, theoretical calculations and QTAIM molecular integrations were used to study a novel class of dihydrogen-bonded complexes (mC(2)H(5)(+)···nMgH(2) with m = 1 or 2 and n = 1 or 2) based in the insight that MgH(2) can bind with the non-localized hydrogen H(+δ) of the ethyl cation (C(2)H(5)(+)). In an overview, QTAIM calculations were applied to evaluate the molecular topography, charge density, as well as to interpret the shifted frequencies either to red or blue caused by the formation of the hydrogen bonds and dihydrogen bonds.