Accurate molecular electrostatic potentials based on modified PRDDO/M wave functions. I. Electrostatic potential derived atomic charges

A new approach for the calculation of electrostatic potential derived atomic charges is presented. Based on molecular orbital calculations in the PRDDO/M approximation, the new parametrized electrostatic potential (PESP) method is parametrized against ab initio MP2/6-31G** calculations. For a data set of 820 atoms in 145 molecules containing H, C, N. O, F, P, S, Cl, and Br (including hypervalent species), the PESP method achieves a mean absolute error of 0.037 e⁻ with a correlation coefficient of 0.990. Unlike other approximate approaches, no scaling factor is required to improve the agreement between PESP charges and the underlying ab initio results. PESP calculations are an order of magnitude faster than the simplest ab initio calculation (STO-3G) on large molecules while achieving a level of accuracy that rivals much more elaborate ab initio methods. © 1997 by John Wiley & Sons, Inc. J Comput Chem 18: 955–969, 1997     

Accurate molecular electrostatic potentials based on modified PRDDO/M wave functions: II. Electrostatic potentials inside the molecular van der Waals envelope

In part I of this series, the PESP (parameterized electrostatic potential) method was described and applied to the calculation of electrostatic-potential-derived charges for a wide variety of organic and inorganic systems. Based on PRDDO/M wave functions and parameterized against ab initio MP2/6-31G** calculations, PESP is an order of magnitude faster than ab initio STO-3G calculations, while achieving a level of accuracy that rivals that of far more sophisticated ab initio methods. In this study, the application of the PESP method to the high potential regions of molecules containing H, C, N, O, F, P, S, Cl, and Br is described. For a collection of 48 molecules and 55 distinct lone pair minima, PESP yields the location and depth of lone pair minima to an average accuracy (relative to MP2/6-31G**) of 0.03 Å and 2.5 kcal/mol, respectively. Similarly, the location and well depths of minima in the π regions of organic molecules are calculated to an accuracy of 0.08 Å and 1.5 kcal/mol. PESP electrostatic potential maps are, in some cases, virtually indistinguishable from those obtained at the MP2/6-31G** level. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1682–1693, 1997     

金属离子(Na~+、K~+、Ca~(2+)、Mg~(2+)、Zn~(2+))与鸟嘌呤异构体配合物的稳定性

在B3LYP/6-311++G**水平上用极化连续介质模型(PCM)系统研究了金属离子(M+/2+=Na+,K+,Ca2+,Mg2+,Zn2+)和十三种鸟嘌呤异构体形成的配合物GnxM+/2+(n为鸟嘌呤异构体的编号,x表示M+/2+与鸟嘌呤异构体的结合位点)在气(g)液(a)两相中的稳定性顺序.着重探讨了液相中配合物的稳定性差异,并且从溶质-溶剂效应、结合能、形变能及异构体的相对能量等几个方面分析了造成稳定顺序发生变化的原因.报道了溶液中这五种金属离子与鸟嘌呤异构体结合形成的六种基态配合物:aG1N2,N3Na+,aG1N2,N3K+,aG1O6,N7Ca2+,aG1N2,N3Mg2+(aG1O6,N7Mg2+),aG2N3,N9Zn2+.可以看出,除了在Zn2+配合物中鸟嘌呤异构体为G2外,构成其余四种金属离子配合物的鸟嘌呤异构体都是G1,但结合位点不同.同时对气相中各类配合物稳定性也进行了系统的排序,并报道了几种较稳定的配合物,如:gG3N1,O6K+,gG5N1,O6K+,gG3N1,O6Ca2+/Mg2+,gG4O6,N7Ca2+/Mg2+.     

Effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water coordination on the structure of glycine and zwitterionic glycine

Interactions between metal ions and amino acids are common both in solution and in the gas phase. Here, the effect of metal ions and water on the structure of glycine is examined. The effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water on structures of Gly.Mn+(H2O)m and GlyZwitt.Mn+(H2O)m (m = 0, 2, 5) complexes have been determined theoretically by employing the hybrid B3LYP exchange-correlation functional and using extended basis sets. Selected calculations were carried out also by means of CBS-QB3 model chemistry. The interaction enthalpies, entropies, and Gibbs energies of eight complexes Gly.Mn+ (Mn+ = Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) were determined at the B3LYP density functional level of theory. The computed Gibbs energies DeltaG degrees are negative and span a rather broad energy interval (from -90 to -1100 kJ mol(-1)), meaning that the ions studied form strong complexes. The largest interaction Gibbs energy (-1076 kJ mol(-1)) was computed for the NiGly2+ complex. Calculations of the molecular structure and relative stability of the Gly.Mn+(H2O)m and GlyZwitt.Mn+(H2O)m (Mn+ = Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+; m = 0, 2, and 5) systems indicate that in the complexes with monovalent metal cations the most stable species are the NO coordinated metal cations in non-zwitterionic glycine. Divalent cations Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+ prefer coordination via the OO bifurcated bonds of the zwitterionic glycine. Stepwise addition of two and five water molecules leads to considerable changes in the relative stability of the hydrated species. Addition of two water molecules at the metal ion in both Gly.Mn+ and GlyZwitt.Mn+ complexes reduces the relative stability of metallic complexes of glycine. For Mn+ = Li+ or Na+, the addition of five water molecules does not change the relative order of stability. In the Gly.K+ complex, the solvation shell of water molecules around K+ ion has, because of the larger size of the potassium cation, a different structure with a reduced number of hydrogen-bonded contacts. This results in a net preference (by 10.3 kJ mol(-1)) of the GlyZwitt.K+H2O5 system. Addition of five water molecules to the glycine complexes containing divalent cations Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+ results in a net preference for non-zwitterionic glycine species. The computed relative Gibbs energies are quite high (-10 to -38 kJ mol(-1)), and the NO coordination is preferred in the Gly.Mn+(H2O)5 (Mn+ = Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) complexes over the OO coordination.     

Effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water coordination on the structure and properties of L-arginine and zwitterionic L-arginine

Interactions between metal ions and amino acids are common both in solution and in the gas phase. The effect of metal ions and water on the structure of L-arginine is examined. The effects of metal ions (Li(+), Na(+), K(+), Mg(2+), Ca(2+), Ni(2+), Cu(2+), and Zn(2+)) and water on structures of Arg x M(H2O)m , m = 0, 1 complexes have been determined theoretically by employing the density functional theories (DFT) and using extended basis sets. Of the three stable complexes investigated, the relative stability of the gas-phase complexes computed with DFT methods (with the exception of K(+) systems) suggests metallic complexes of the neutral L-arginine to be the most stable species. The calculations of monohydrated systems show that even one water molecule has a profound effect on the relative stability of individual complexes. Proton dissociation enthalpies and Gibbs energies of arginine in the presence of the metal cations Li(+), Na(+), K(+), Mg(2+), Ca(2+), Ni(2+), Cu(2+), and Zn(2+) were also computed. Its gas-phase acidity considerably increases upon chelation. Of the Lewis acids investigated, the strongest affinity to arginine is exhibited by the Cu(2+) cation. The computed Gibbs energies DeltaG(o) are negative, span a rather broad energy interval (from -150 to -1500 kJ/mol), and are appreciably lowered upon hydration.     

Effect of ring annelation on cations [M = H+, Li+, Na+, K+, Be2+, Mg2+, and Ca2+]…benzene interaction: A density functional theory investigation

Density functional theory calculation was carried out on cation‐π complexes formed by cations [M = H⁺, Li⁺, Na⁺, K⁺, Be²⁺, Mg²⁺, and Ca²⁺] and π systems of annelated benzene. The cation‐π bonding energy of Be²⁺ or Mg²⁺ with annelated benzene is very strong in comparison with the common cation‐π intermolecular interaction, and the bonding energies follow the order Be²⁺ > Mg²⁺ > Ca²⁺ > Li⁺ > Na⁺ > K⁺. Similarly, the interaction energies follow the trend 1‐M < 2‐M < 3‐M for all the metal cations considered. These outcomes may be due to the weak interactions of the metal cations with C? H and the interactions of metal cations with π in addition to the nature of a metal cation. We have also investigated on all the possible substituted sites, and find that the metal ion tends to interact with all ring atoms while proton prefers to bind covalently to one of the ring carbons. The binding of metal cations with annelated benzenes has striking effect on nuclear magnetic resonance chemical shifts using the gauge independent atomic orbital method. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Sandwich-like compounds based on the all-metal aromatic unit Al(4)2- and the main-group metals M (M=Li, Na, K, Be, Mg, Ca)

Inspired by the pioneering experimental characterisation of the all-metal aromatic unit Al(4)2- in the bimetallic molecules MAl4- (M=Li, Na, Cu) and by the very recent theoretical design of sandwich-type transition-metal complexes [Al4MAl4]q- (q=0-2; M=Ti, V, Cr, Zr, Nb, Mo, Hf, Ta, W), we used density functional theory (DFT) calculations (B3LYP/6-311+G(d) to design a series of novel non-transition-metal sandwich complexes based on the all-metal aromatic unit Al4(2-) and the main-group metals M (M=Li, Na, K, Be, Mg, Ca). The traditional homo-decked sandwich compounds [Al4MAl4]q- (without counterions) and (nM)q+[Al4MAl4]q- (with counterions M) (q=2-3, M=Li, Na, K, Be, Mg, Ca), although some of them are truly energy minima, have a much higher energy than many fused isomers. We thus concluded that it seems unlikely for Al4(2-) to sandwich the main-group metal atoms in the homo-decked sandwich form. Alternatively, we proposed a new type of sandwich complex, namely hetero-decked sandwich compounds [CpMAl4]q-, that are the ground-state structures for each M both with and without counterions. It was shown that with the rigid Cp- partner, the all-metal aromatic unit Al(4)2- might indeed act as a "superatom". These new types of all-metal aromatic unit-based sandwich complexes await future experimental verification.     

Geometries, stability and aromaticity of , [M(Al2P2)] (M = Li, Na, K, Cu) and N(Al2P2) (N = Be, Mg, Ca, Zn) clusters

Geometrical structure, aromaticity and other properties of , [M(Al₂P₂)]⁻ (M = Li, Na, K, Cu) and N(Al₂P₂) (N = Be, Mg, Ca, Zn) species are theoretically investigated with density functional theory (DFT) methods. Calculation results show that for species, the planar structure, with D_2h symmetry at the ¹A_g state, is the global minimum at the B3LYP/6-311+G_* level. Natural bond orbital (NBO) analysis indicates the existence of delocalization in the most stable species and its pyramidal complexes. Nucleus-independent chemical shift (NICS) and molecular orbital (MO) analysis further reveal that that pyramidal [M(Al₂P₂)]⁻ and N(Al₂P₂) species preserve the aromatic nature of the most stable unit.