Structures,vibrational frequencies,topologies, and energies of hydrogen bonds in cysteine‐formaldehyde complexes

The hydrogen bonding interactions between cysteine (Cys) and formaldehyde (FA) were studied with density functional theory regarding their geometries, energies, vibrational frequencies, and topological features of the electron density. The quantum theory of atoms in molecules and natural bond orbital analyses were employed to elucidate the interaction characteristics in the Cys‐FA complexes. The intramolecular hydrogen bonds (H‐bonds) formed between the hydroxyl and the N atom of cysteine moiety in some Cys‐FA complexes were strengthened because of the cooperativity. Most of intermolecular H‐bonds involve the O atom of cysteine/FA moiety as proton acceptors, while the strongest H‐bond involves the O atom of FA moiety as proton acceptor, which indicates that FA would rather accept proton than providing one. The H‐bonds formed between the CH group of FA and the S atom of cysteine in some complexes are so weak that no hydrogen bonding interactions exist among them. In most of complexes, the orbital interaction of H‐bond is predominant during the formation of complex. The electron density (ρ_b) and its Laplace (?²ρ_b) at the bond critical point significantly correlate with the H‐bond parameter δR, while a linearly relationship between the second‐perturbation energy E(2) and ρ_b has been found as well. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Density functional theory and topological analysis on the hydrogen bonding interactions in cysteine‐thymine complexes

Hydrogen bonding interactions between amino acids and nucleic acid bases constitute the most important interactions responsible for the specificity of protein binding. In this study, complexes formed by hydrogen bonding interactions between cysteine and thymine have been studied by density functional theory. The relevant geometries, energies, and IR characteristics of hydrogen bonds (H‐bonds) have been systematically investigated. The quantum theory of atoms in molecule and natural bond orbital analysis have also been applied to understand the nature of the hydrogen bonding interactions in complexes. More than 10 kinds of H‐bonds including intra‐ and intermolecular H‐bonds have been found in complexes. Most of intermolecular H‐bonds involve O (or N) atom as H‐acceptor, whereas the H‐bonds involving C or S atom usually are weaker than other ones. Both the strength of H‐bonds and the structural deformation are responsible for the stability of complexes. Because of the serious deformation, the complex involving the strongest H‐bond is not the most stable structures. Relationships between H‐bond length (ΔR_X‐H), frequency shifts (Δv), and the electron density (ρ_b) and its Laplace (?²ρ_b) at bond critical points have also been investigated. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

β‐Aminoacrolein: An ab initio,AIM and NBO study

The molecular structure and the intramolecular hydrogen bonding of β‐aminoacrolein and its simple derivatives were investigated at the MP2 and B3LYP levels of theory using the standard 6‐311++G(d, p) basis set. The “atoms in molecules” or AIM theory of Bader which is based on topological properties of the electron density (ρ), was used. Additionally, an analysis of the critical points was performed to study the nature hydrogen bonding in these systems. Natural bond orbital (NBO) analysis was also carried out for to better comprehend the nature of the intramolecular interactions in β‐aminoacrolein and its derivatives. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Halogen Bonding： An AIM Analysis of the Weak Interactions

A series of complexes formed between halogen-containing molecules and ammonia have been investigated by means of the atoms in molecules （AIM） approach to gain a deeper insight into halogen bonding. The existence of the halogen bond critical points （XBCP） and the values of the electron density （Pb） and Laplacian of electron density （V2pb） at the XBCP reveal the closed-shell interactions in these complexes. Integrated atomic properties such as charge, energy, polarization moment, volume of the halogen bond donor atoms, and the corresponding changes （△） upon complexation have been calculated. The present calculations have demonstrated that the halogen bond represents different AIM properties as compared to the well-documented hydrogen bond. Both the electron density and the Laplacian of electron density at the XBCP have been shown to correlate well with the interaction energy, which indicates that the topological parameters at the XBCP can be treated as a good measure of the halogen bond strength In addition, an excellent linear relationship between the interatomic distance d（X…N） and the logarithm of Pb has been established.     

Inequivalence of substitution pairs in hydroxynaphthaldehyde: A theoretical measurement by intramolecular hydrogen bond strength,aromaticity, and excited‐state intramolecular proton transfer reaction

The inequivalence of substitution pair positions of naphthalene ring has been investigated by a theoretical measurement of hydrogen bond strength, aromaticity, and excited state intramolecular proton transfer (ESIPT) reaction as the tools in three substituted naphthalene compounds viz 1‐hydroxy‐2‐naphthaldehyde (HN12), 2‐hydroxy‐1‐naphthaldehyde (HN21), and 2‐hydroxy‐3‐naphthaldehyde (HN23). The difference in intramolecular hydrogen bond (IMHB) strength clearly reflects the inequivalence of substitution pairs where the calculated IMHB strength is found to be greater for HN12 and HN21 than HN23. The H‐bonding interactions have been explored by calculation of electron density ρ(r) and Laplacian ?²ρ(r) at the bond critical point using atoms in molecule method and by calculation of interaction between σ* of OH with lone pair of carbonyl oxygen atom using NBO analysis. The ground and excited state potential energy surfaces (PESs) for the proton transfer reaction at HF (6‐31G**) and DFT (B3LYP/6‐31G**) levels are similar for HN12, HN21 and different for HN23. The computed aromaticity of the two rings of naphthalene moiety at B3LYP/6‐31G** method also predicts similarity between HN12 and HN21, but different for HN23. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Synthesis and Theoretical Study of Intramolecular Hydrogen Bond at Two Possible Positions in Pyrazolo[1,2‐b]phthalazine

Properties of dimethyl 3‐(alkylamino)‐5,10‐dioxo‐5,10‐dihydro‐1H‐pyrazolo[1,2‐b]phthalazine‐1,2‐dicarboxylate and its derivatives were studied by means of ab initio method. NO₂ derivative of title compound was synthesized and the nature of its intramolecular hydrogen bond (HB) was investigated. Furthermore, the topological properties of the electron density distributions for N? H···O intramolecular bridges were analyzed in terms of the Bader theory of atoms in molecules (AIM). The electron density (ρ) and Laplacian (?²ρ) properties, estimated by AIM calculations, indicated that O···H bond possesses low ρ and positive ?²ρ values which are in agreement with electrostatic character of the HBs, whereas N? H bonds have covalent character (?²ρ<0). Moreover, steric effect of the t‐Bu group on structure and topological parameters of pyrazolo[1,2‐b]phthalazine conformers was studied. Finally, the powerful method of Espinosa was used to obtain the H‐bond energy.     

Hydrogen Bonding Interaction of Formic Acid-, Formaldehyde-, Formylfluoride-Nitrosyl Hydride： Theoretical Study on the Geometries, Interaction Energies and Blue- or Red-Shifted Hydrogen Bonds

The hydrogen bonding interaction of formic acid-, formaldehyde-, formylfluoride-nitrosyl hydride complexes was investigated by the density functional theory （DFT） and ab inito method in conjunction with 6-311＋＋G（2d,2p） basis set. The geometries, vibrational frequencies and interaction energies of the complexes were calculated by both standard and CP-corrected methods respectively. Moreover, G3B3 method was employed to estimate the interaction energies. There are C--H…O, N--H…O, N--H…F blue-shifted H-bonds and red-shifted O----H…O H-bond in the complexes. Electron density redistribution and rehybridization contribute to the N--H and C--H blue shifts. All geometric reorganizations contribute to the N--H blue shifts and partial geometric reorganizations contribute to the C--H blue shifts. The geometric reorganizations of the complex C except ZH（5）-O（4）-C（1） contribute to the O----H red shift. For the N--H blue shifts, the effect of r（N--O） variation on the N--H blue shifts is larger than that of ZH-N-O variation. Rehybridization plays a dominant role in the degree of N--H blue shifts, whereas the electron density redistribution contributes more to the degree of C--H blue shifts than the other effects do.     

Assessment of intermolecular interactions at three sites of the arylalkyne in phenylacetylene‐containing lithium‐bonded complexes: Ab initio and QTAIM studies

The intermolecular interactions existing at three different sites between phenylacetylene and LiX (X = OH, NH₂, F, Cl, Br, CN, NC) have been investigated by means of second‐order Møller?Plesset perturbation theory (MP2) calculations and quantum theory of “atoms in molecules” (QTAIM) studies. At each site, the lithium‐bonding interactions with electron‐withdrawing groups (? F, ? Cl, ? Br, ? CN, ? NC) were found to be stronger than those with electron‐donating groups (? OH and ? NH₂). Molecular graphs of C₆H₅C?CH···LiF and πC₆H₅C?CH···LiF show the same connectional positions, and the electron densities at the lithium bond critical points (BCPs) of the πC₆H₅C?CH···LiF complexes are distinctly higher than those of the σC₆H₅C?CH···LiF complexes, indicating that the intermolecular interactions in the C₆H₅C?CH···LiX complexes can be mainly attributed to the π‐type interaction. QTAIM studies have shown that these lithium‐bond interactions display the characteristics of “closed‐shell” noncovalent interactions, and the molecular formation density difference indicates that electron transfer plays an important role in the formation of the lithium bond. For each site, linear relationships have been found between the topological properties at the BCP (the electron density ρ_b, its Laplacian ?²ρ_b, and the eigenvalue λ₃ of the Hessian matrix) and the lithium bond length d(Li‐bond). The shorter the lithium bond length d(Li‐bond), the larger ρ_b, and the stronger the π···Li bond. The shorter d(Li‐bond), the larger ?²ρ_b, and the greater the electrostatic character of the π···Li bond. © 2012 Wiley Periodicals, Inc.     

DFT study of nitroxide radicals: explicit modeling of solvent effects on the structural and electronic characteristics of 4‐amino‐2,2,6,6‐tetramethyl‐piperidine‐N‐oxyl

An explicit DFT modeling of water surroundings on the electron paramagnetic resonance properties of 4‐amino‐2,2,6,6‐tetramethyl‐piperidine‐N‐oxyl (TA) has been performed. A stepwise hydration of TA is accompanied with certain changes in geometrical parameters (bond lengths and angles) and redistribution of partial electric charges in TA. An aqueous cluster of 45 water molecules can be considered as an appropriate model for a complete aqueous shell around TA, although most of the structural and electronic characteristics of TA already converge at about 10 water molecules. Water surroundings induce an increase in electron spin density on the nitrogen atom of the nitroxide fragment due to stabilization of the polar resonance structure > N^+?? O^? at the expense of less polar structure > N? O^?. The water‐induced rise of the isotropic splitting constant a_iso, calculated from the contact term of the hyperfine interaction, comprises Δa_iso(ρ_N2) = 2.2–2.5 G, which is typical of experimental value for TA. There are two contributions to the solvent effect on the a_iso(ρ_N2) value: the redistribution of spin density in the nitroxide fragment (polarity effect) and water‐induced distortions of TA geometry. Microscopic variations in a hydrogen‐bonded water network cause noticeable fluctuations of the splitting constant a_iso(ρ_N2). Calculations of the atomic spin density (σ_N2) allowed us to compute the splitting constant from the relationship a_iso(σ_N2) = Qσ_N2, where Q = 36.2 G. A practical advantage of using this relationship is that it gives ‘smoothed’ values of the splitting constant, which are sensitive to the environment polarity but remain tolerant to microscopic fluctuations of the hydrogen‐bonded water network around a spin‐label molecule. Copyright © 2010 John Wiley & Sons, Ltd.     

Nature of Noncovalent Carbon‐Bonding Interactions Derived from Experimental Charge‐Density Analysis

In an effort to better understand the nature of noncovalent carbon‐bonding interactions, we undertook accurate high‐resolution X‐ray diffraction analysis of single crystals of 1,1,2,2‐tetracyanocyclopropane. We selected this compound to study the fundamental characteristics of carbon‐bonding interactions, because it provides accessible σ holes. The study required extremely accurate experimental diffraction data, because the interaction of interest is weak. The electron‐density distribution around the carbon nuclei, as shown by the experimental maps of the electrophilic bowl defined by a (CN)₂C?C(CN)₂ unit, was assigned as the origin of the interaction. This fact was also evidenced by plotting the Δ²ρ(r) distribution. Taken together, the obtained results clearly indicate that noncovalent carbon bonding can be explained as an interaction between confronted oppositely polarized regions. The interaction is, thus electrophilic–nucleophilic (electrostatic) in nature and unambiguously considered as attractive.     

Single‐electron halogen bond: Ab initio study

Ab initio calculations have been performed on single‐electron halogen bonds between methyl radical and bromine‐containing molecules to gain a deeper insight into the nature of such noncovalent interactions. Bader's atoms in molecules (AIM) theory have also been applied to the analysis of the linking of the single‐electron halogen bond. Various characteristics of the R? Br…CH₃ interaction, i.e., binding energies, geometrical parameters and topological properties of the electron density have been determined. The presence of the bond critical points (BCPs) between the bromine atom and methyl radical and the values of electron density and Laplacian of electron density at these BCPs indicate the closed‐shell interactions in the complexes. The single‐electron halogen bonds, which are significantly weaker than the normal halogen bonds, exhibit equally bond strength as compared to the single‐electron hydrogen bond. It has been also found that plotting of the binding energies versus topological properties of the electron density at the BCPs gives two straight lines. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Substituent effect on structure,electron density,and intramolecular hydrogen bonding in nitroso‐oxime methane

Density functional calculations with Beck's three‐parameter hybrid method using the correlation functional of Lee, Yang, and Parr (B3LYP) were carried out for investigation of the intramolecular hydrogen bond strength in Nitroso‐oxime methane and its derivatives. Also, vibrational frequencies for them were calculated at the same level of theory. The π‐electron delocalization parameter (Q) and as a geometrical indicator of a local aromaticity, the geometry‐based harmonic oscillator measure of aromaticity index has been applied. Additionally, the linear correlation coefficients between substituent constants and selected parameters in R position have calculated. The obtained results show that the hydrogen bond strength is mainly governed by the resonance variations inside the chelate ring induced by the substituent groups. The topological properties of the electron density distributions for O? H ··· O intramolecular bridges have been analyzed in terms of the Bader theory of atoms in molecules (AIM). Correlations between the H‐bond strength and topological parameters have been also studied. The electron density (ρ) and Laplacian (?²ρ) properties, estimated by AIM calculations, show that O ··· H bond have low ρ and negative (?²ρ) values (consistent with covalent character of the HBs), whereas O? H bond have positive (?²ρ) Furthermore, the analysis of hydrogen bond in this molecule and its derivatives by quantum theory of natural bond orbital (NBO) methods fairly support the ab initio results. Natural population analysis data, the electron density, and Laplacian properties as well as υ(O? H) and γ(O? H) were further used for estimation of the hydrogen bonding interactions and the forces driving their formation. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Surface Charge‐Transfer Doping of Graphene Nanoflakes Containing Double‐Vacancy (5‐8‐5) and Stone–Wales (55‐77) Defects through Molecular Adsorption

The adsorption of six electron donor–acceptor (D/A) organic molecules on various sizes of graphene nanoflakes (GNFs) containing two common defects, double‐vacancy (5‐8‐5) and Stone–Wales (55‐77), are investigated by means of ab initio DFT [M06‐2X(‐D3)/cc‐pVDZ]. Different D/A molecules adsorb on a defect graphene (DG) surface with binding energies (ΔE_b) of about ?12 to ?28 kcal mol^?1. The ΔE_b values for adsorption of molecules on the Stone–Wales GNF surface are higher than those on the double vacancy GNF surface. Moreover, binding energies increase by about 10 % with an increase in surface size. The nature of cooperative weak interactions is analyzed based on quantum theory of atoms in molecules, noncovalent interactions plot, and natural bond order analyses, and the dominant interaction is compared for different molecules. Electron density population analysis is used to explain the n‐ and p‐type character of defect graphene nanoflakes (DGNFs) and also the change in electronic properties and reactivity parameters of DGNFs upon adsorption of different molecules and with increasing DGNF size. Results indicate that the HOMO–LUMO energy gap (E_g) of DGNFs decreases upon adsorption of molecules. However, by increasing the size of DGNFs, the E_g and chemical hardness of all complexes decrease and the electrophilicity index increases. Furthermore, the values of the chemical potential of acceptor–DGNF complexes decrease with increasing size, whereas those of donor–DGNF complexes increase.     

Comparative analysis of blue‐shifted hydrogen bond versus conventional hydrogen bond in methyl radical complexes

Methyl radical complexes H₃C…HCN and H₃C…HNC have been investigated at the UMP2(full)/aug‐cc‐pVTZ level to elucidate the nature of hydrogen bonds. To better understand the intermolecular H‐bond interactions, topological analysis of electron density at bond critical points (BCP) is executed using Bader's atoms‐in‐molecules (AIM) theory. Natural bond orbital (NBO) analysis has also been performed to study the orbital interactions and change of hybridization. Theoretical calculations show that there is no essential difference between the blue‐shift H‐bond and the conventional one. In H₃C…HNC complex, rehybridization is responsible for shortening of the N? H bond. The hyperconjugative interaction between the single electron of the methyl radical and N? H antibonding orbital is up to 7.0 kcal/mol, exceeding 3.0 kcal/mol, the upper limit of hyperconjugative n(Y)→σ*(X–H) interaction to form the blue‐shifted H‐bond according to Alabugin's theory. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Structure‐Based Design of Platinum(II) Complexes as c‐myc Oncogene Down‐Regulators and Luminescent Probes for G‐Quadruplex DNA

A series of platinum(II) complexes with tridentate ligands was synthesized and their interactions with G‐quadruplex DNA within the c‐myc gene promoter were evaluated. Complex 1 , which has a flat planar 2,6‐bis(benzimidazol‐2‐yl)pyridine (bzimpy) scaffold, was found to stabilize the c‐myc G‐quadruplex structure in a cell‐free system. An in silico G‐quadruplex DNA model has been constructed for structure‐based virtual screening to develop new Pt^II‐based complexes with superior inhibitory activities. By using complex 1 as the initial structure for hit‐to‐lead optimization, bzimpy and related 2,6‐bis(pyrazol‐3‐yl)pyridine (dPzPy) scaffolds containing amine side‐chains emerge as the top candidates. Six of the top‐scoring complexes were synthesized and their interactions with c‐myc G‐quadruplex DNA have been investigated. The results revealed that all of the complexes have the ability to stabilize the c‐myc G‐quadruplex. Complex 3 a ([Pt^II L2R ] ⁺ ; L2 =2,6‐bis[1‐(3‐piperidinepropyl)‐1H‐enzo[d]imidazol‐2‐yl]pyridine, R =Cl) displayed the strongest inhibition in a cell‐free system (IC₅₀=2.2 μM ) and was 3.3‐fold more potent than that of 1 . Complexes 3 a and 4 a ([Pt^II L3R ]⁺; L3 =2,6‐bis[1‐(3‐morpholinopropyl)‐1H‐pyrazol‐3‐yl]pyridine, R =Cl) were found to effectively inhibit c‐myc gene expression in human hepatocarcinoma cells with IC₅₀ values of ≈17 μM , whereas initial hit 1 displayed no significant effect on gene expression at concentrations up to 50 μM . Complexes 3 a and 4 a have a strong preference for G‐quadruplex DNA over duplex DNA, as revealed by competition dialysis experiments and absorption titration; 3 a and 4 a bind G‐quadruplex DNA with binding constants (K) of approximately 10⁶–10⁷ dm³ mol^?1, which are at least an order of magnitude higher than the K values for duplex DNA. NMR spectroscopic titration experiments and molecular modeling showed that 4 a binds c‐myc G‐quadruplex DNA through an external end‐stacking mode at the 3′‐terminal face of the G‐quadruplex. Intriguingly, binding of c‐myc G‐quadruplex DNA by 3 b is accompanied by an increase of up to 38‐fold in photoluminescence intensity at λ_max=622 nm.     

Theoretical investigation on charge‐assisted halogen bonding interactions in the complexes of bromocarbons with some anions

A series of dimeric complexes formed between bromocarbon molecules and two anions (Br^? and CN^?) have been investigated by using MP2 method. The quantum theory of atoms in molecules (QTAIM) and the second‐order perturbation natural bond orbital (NBO) approaches were applied to analyze the electron density distributions of these complexes and to explore the nature of charge‐assisted halogen bonding interactions. As anticipated, these interactions are significantly stronger relative to the corresponding neutral ones. The results derived from ab initio calculations described herein reveal a major contribution from the electrostatic interaction on the stability of the systems considered. Beside the electrostatic interaction, the charge‐transfer force and the second‐order orbital interaction also play an important role in the formation of the complexes, as a NBO analysis suggested. The presence of halogen bonds in the complexes has been identified in terms of the QTAIM methodology, and several linear relationships have been established to provide more insight into charge‐assisted halogen bonding interactions. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Hydrogen‐bonded complexes of nicotine with simple alcohols

Ab initio and density functional theory studies have been performed on the hydrogen‐bonded complexes of neutral and protonated nicotine with ethanol, methanol, and trifluromethanol to explore their relative stability in a systematic way. Among all the hydrogen‐bonded nicotine complexes considered here, protonated forms in nicotine–ethanol and nicotine–methanol, and neutral form in nicotine–trifluromethanol complexes have been found to be the most stable. In the former two complexes, the proton attached to the pyrrolidine nitrogen acts as a strong hydrogen bond donor, whereas the pyrrolidine nitrogen atom acts as a hydrogen bond acceptor in the latter case. Neutral complex of nicotine with trifluromethanol has been found to possess a very short hydrogen bond (1.57 Å) and basis set superposition error corrected hydrogen bond energy value of 19 kcal/mol. The nature of the various hydrogen bonds formed has been investigated through topological aspects using Bader's atoms in molecules theory. From the calculated topological results, excellent linear correlation is shown to exist among the hydrogen bond length, electron density, and its Laplacian at the bond critical points for all the complexes considered. The natural bond orbital analysis has been carried out to investigate the charge transfer in the nicotine alcohol complexes. In contrast to the blue shifting behavior that is generally exhibited by other C? H···O hydrogen bonds involving sp³ carbon atom, the C? H···O hydrogen bond in the protonated nicotine–ethanol and methanol complexes has been found to be proper with red shifting in nature. © 2011 Wiley Periodicals, Inc.     

Comparative Analysis of Electron‐Density and Electron‐Localization Function for Dinuclear Manganese Complexes with Bridging Boron‐ and Carbon‐Centered Ligands

Bonding in borylene‐, carbene‐, and vinylidene‐bridged dinuclear manganese complexes [MnCp(CO)₂]₂X (X=B‐tBu, B=NMe₂, CH₂, C?CH₂) has been compared by analyses based on quantum theory of atoms in molecules (QTAIM), on the electron‐localization function (ELF), and by natural‐population analyses. All of the density functional theory based analyses agree on the absence of a significant direct Mn? Mn bond in these complexes and confirm a dominance of delocalized bonding via the bridging ligand. Interestingly, however, the topology of both charge density and ELF related to the Mn‐bridge‐Mn bonding depend qualitatively on the chosen density functional (except for the methylene‐bridged complex, which exhibits only one three‐center‐bonding attractor both in ??²ρ and in ELF). While gradient‐corrected functionals provide a picture with localized two‐center X? Mn bonding, increasing exact‐exchange admixture in hybrid functionals concentrates charge below the bridging atom and suggests a three‐center bonding situation. For example, the bridging boron ligands may be described either as substituted boranes (e.g., at BLYP or BP86 levels) or as true bridging borylenes (e.g., at BHLYP level). This dependence on the theoretical level appears to derive from a bifurcation between two different bonding situations and is discussed in terms of charge transfer between X and Mn, and in the context of self‐interaction errors exhibited by popular functionals.     

A comparative study of some lithium and hydrogen‐bonded complexes: Ab initio and QTAIM studies

The nature of the interactions of cyanide with lithium and hydrogen halides was investigated using ab initio calculations and topological analysis of electron density. The computed properties of the lithium‐bonded complexes RCN···LiX (R = H, F, Cl, Br, C?CH, CH?CH₂, CH₃, C₂H₅; X = Cl, Br) were compared with those of corresponding hydrogen‐bonded complexes RCN···HX. The results show that both types of intermolecular interactions are “closed‐shell” noncovalent interactions. The effect of substitution on the interaction energy and electron density at the bond critical points of the lithium and hydrogen bonding interactions is similar. In comparison, the interaction energies of lithium‐bonded complexes are more negative than those of hydrogen‐bonded counterparts. The electrostatic interaction plays a more important role in the lithium bond than in the hydrogen bond. On complex formation, the net charge and energy of the Li atom decrease and the atomic volume increases, while the net charge and energy of the H atom increase and the atomic volume decreases. © 2013 Wiley Periodicals, Inc.     

Study on the Cation-π Interactions between Ammonium Ion and Aromatic π Systems

The nature and strength of the cation-π interactions between NH4^＋ and toluene, p-cresol, or Me-indole were studied in terms of the topological properties of molecular charge density and binding energy decomposition. The results display that the diversity in the distribution pattern of bond and cage critical points reflects the profound influence of the number and nature of substituent on the electron density of the aromatic rings. On the other hand, the energy decomposition shows that dispersion and repulsive exchange forces play an important role in the organic cation （NH4^＋）-π interaction, although the electrostatic and induction forces dominate the interaction. In addition, it is intriguing that there is an excellent correlation between the electrostatic energy and ellipticity at the bond critical point of the aromatic π systems, which would be helpful to further understand the electrostatic interaction in the cation-π complexes.