The solvation of Li and Na in acetonitrile from ab initio-derived many-body ion–solvent potentials

Several Li⁺- and Na⁺-acetonitrile models were derived from ab initio calculations at the counterpoise-corrected MP2/TZV++(d,p) level for distorted ion-(MeCN)_n clusters with n=1, 4 and 6. Two different many-body ion-acetonitrile models were constructed: an effective three-body potential for use with the six-site effective pair model of Böhm et al., and an effective polarizable many-body model. The polarizable acetonitrile model used in the latter model is a new empirical model which was also derived in the present paper. Mainly for comparative purposes, two ion-acetonitrile pair potentials were also constructed from the ab initio cluster calculations: one pure pair potential and one effective pair potential. Using all these potential models, MD simulations in the NPT ensemble were performed for the pure acetonitrile liquid and for Li⁺(MeCN) and Na⁺(MeCN) solutions with 1 ion in 512 solvent molecules and with a simulation time of at least 120 ps per system. Thermodynamic properties, solvation-shell structure and the self-diffusion coefficient of the ions and of the solvent molecules were calculated and compared between the different models and with experimental data, where available. The Li⁺ ion is found to be four-coordinated when the new many-body potentials are used, in contrast to the six-coordinated structure obtained for the pure pair and effective pair potentials. The coordination number of Na⁺ is close to six for all the models derived here, although the coordination number becomes slightly smaller with the many-body potentials. For both ions, the solvent molecules in the first shell point their nitrogen ends towards the cation, while in the second shell the opposite orientation is the most common.     

An ab initio potential energy surface for the υ4 vibrational mode of hydrogen peroxide

In this paper, the levels and the torsional microstates of hydrogen peroxide are determined from fully optimized ab initio calculations using a nuclear model in one dimension. Calculations have been performed at the MP2 level with the 6-311 G(2df,2pd), 6-31 1+G(2df,2pd), cc-pVTZ and AUG-cc-pVTZ basis sets including polarization orbitals and diffuse functions. The most stable conformation, calculated with the MP2/AUG-cc-pVTZ approach, is a trans–gauche conformer lying at 67.5° from the trans structure. By using the same level of calculations, the heights of the trans and cis barriers have been determined to be 386.5 and 2643.8 cm⁻¹ in a good agreement with the experimental data. The variational torsional levels split into four components by the tunnelling effect of the barriers. The splitting of the fundamental level caused by the trans barrier has been found to be 11.8683 cm⁻¹, whereas the splitting caused by the cis barrier is insignificant under n=2. Current ab initio energies confirm the experimental assignments and verify the separability of the torsion from the rest of the vibrations. However, the experimental relation of dependence on the torsion of the rotational constants cannot be reproduced in one-dimension and depends on several additional vibrational effects.     

Applicability of the H NMR chemical shifts computed by the ab initio/GIAO-HF methodology to the study of geometrical features of Zn-porphyrin dimers

Several porphyrin dimers have been newly designed and synthesized to construct assemblies with 1,4-diazabicyclo[2.2.2]octane (DABCO) as a bidentate binding ligand. Semi empirical (AM1) and ab initio calculations have been used to study the assemblies generated by the organization of dimers and DABCO, including the computation of ¹H NMR complexation-induced chemical shifts using the ab initio/GIAO methodology. The diagnostic capacity of the theoretical method has been applied to explain experimental results and geometrical features of the complexes.     

NMR chemical shifts in the low–pH form of a-chymotrypsin. A QM/MM and ONIOM–NMR study

 The relationship between hydrogen bonding and NMR chemical shifts in the catalytic triad of low-pH α-chymotrypsin is investigated by combined use of the effective fragment potential [(2001) J Phys Chem A 105:293] and ONIOM–NMR [(2000) Chem Phys Lett 317:589] methods. Our study shows that while the His57 N^δ1−H bond is stretched by a relatively modest amount (to about 1.060 ?) this lengthening, combined with the polarization due to the molecular environment, is sufficient to explain the experimentally observed chemical shifts of 18.2 ppm. Furthermore, the unusual down-field shift of H^ɛ1 (9.2 ppm) observed experimentally is reproduced and shown to be induced by interactions with the C=O group of Ser214 as previously postulated. The free-energy cost of moving H^δ1 from His57 to Asp102 is predicted to be 5.5 kcal/mol. Received: 26 September 2001 / Accepted: 6 September 2002 / Published online: 21 January 2003 Contribution to the Proceedings of the Symposium on Combined QM/MM Methods at the 222nd National Meeting of the American Chemical Society, 2001 Correspondence to: J. H. Jensen e-mail: jan-jensen@uiowa.edu Acknowledgements. This work was supported by a Research Innovation Award from the Research Corporation and a type G starter grant from the Petroleum Research Fund. The calculations were performed on IBM RS/6000 workstations obtained through a CRIF grant from the NSF (CHE-9974502) and on supercomputers at the National Center for Supercomputer Applications at Urbana-Champaign. The authors are indebted to Visvaldas Kairys for help with the CHARMM program, and to Daniel Quinn for many helpful discussions.     

Quantum chemistry study on B32 and B32H2–32

The structures of B32 and B₃₂H^2–₃₂ with I_h symmetry have been investigated by means of ab initio calculations at STO-3G level. The relationship between molecular orbitals of them has been analyzed and their bonding properties have been discussed. Then the possibility of their existence, as well as the similarity and difference between B₃₂ (B₃₂H^2–₃₂) and C₆₀ (C₆₀H₆₀) have been inferred.     

A Quantum Chemical Study of N14 Cluster

Ab initio molecular orbital theory and density functional theory have been employed to study N₁₄ cluster with low spin at the HF/6-31G*, B3LYP/6-31G*, B3PW91/6-31G*, BP86/6-31G*, and BHLYP/6-31G* levels of theory. Twelve isomers were studied, including one previously investigated cage molecule. The most stable isomer of N₁₄ is a C _2h -symmetric molecule that contains two separated five-membered nitrogen rings connected by a —N=N—N=N— bridge. The second, third, and fifth most stable isomers each have one five-membered nitrogen ring. The theoretical results suggest that the five-membered nitrogen ring gives rise to a particularly stable structural unit, and the more side chains that the five-membered nitrogen ring links with, the less stable the structure will become.     

The mixed valent tellurate SrTe3O8: Electronic lone pair effect of Te

A novel mixed valent tellurium oxide, SrTe₃O₈, has been synthesized and its crystal structure was determined ab initio from powder X-ray diffraction data. This oxide, which crystallizes in a tetragonal unit-cell, P4₂/m space group, with very close a and c cell parameters (6.8257(1) and 6.7603(1) Å, respectively), exhibits a very original structure built up of corner-sharing TeO₆ (Te⁶⁺) octahedra and Te₂O₈ (Te⁴⁺) twin-pyramidal units. The latter ones form [Te₃O₈]_∞ chains running along the [001] and the [110] directions. Besides the four sided tunnels where the Sr²⁺ cations are located, there are very large four sided tunnels running along the c-axis which are obstructed by the electronic lone pairs of the Te⁴⁺ cations.     

Vibrational spectrum of BiH3: Six-dimensional variational calculations on high-level ab initio potential energy surfaces

We report a theoretical study of the ground electronic state of BiH₃. The potential energy surface (PES) is obtained from coupled cluster CCSD(T) calculations with a large basis set (289 contracted Gaussian functions). The previously available quartic force field (P4) is extended by adding the dominant quintic and sextic stretching terms to yield improved potential functions in symmetry coordinates (P6) and Morse-type coordinates (M4). Second-order rovibrational perturbation calculations on the P4-PES and full variational calculations on the P6-PES and M4-PES yield almost identical vibrational term values which is rationalized by considering the local mode behavior of BiH₃ and the Morse-type character of the M4-PES. The remaining deviations between the computed and observed vibrational term values must thus be caused by imperfections in the CCSD(T) surface. A refinement of this ab initio surface by a restrained fit to experimental data allows an essentially perfect reproduction of the observed vibrational term values. Variational calculations on this refined surface provide predictions for several overtone and combination bands that have not yet been observed. Dedicated to Hermann Stoll on the occasion of his 60th birthday     

基于三维静电势参数研究C60溶解性的构效关系

在量子化学从头算基础上,对一系列溶剂化合物分子进行了结构优化和三维静电势参数计算,运用多元线性回归分析和神经网络方法对C60在121种不同溶剂的溶解性与计算的结构参数进行了构效关系研究.用建立起来的QSPR 关系式对15种不同结构类型溶剂进行了预测,并阐述了C60溶质与不同溶剂之间的相互作用,获得了满意的结果.     

C12簇及其碳管烷C12H12的理论研究

C60和C70等碳笼烯的发现及其在新物质、新材料等研究领域中的重要性，激励人们去探索合成更多新型多面体碳簇化合物与碳元素类似物[1，2]．最近，两类新型多面体碳管笼烯的设想提出来了[3，4]应用HMO和MNDO方法对其稳定性变化规律、结构和成键特征进行了讨论.关于多苯的vanderWaals簇实验上已有广泛的调查[5，6]，理论上对其二聚物（C6H6）。的不同几何构型与稳定性进行了深入的研究[7]．本文采用abinitio计算有效势（effectivecorePotential）方法，对C12碳管元素簇和C12H12，碳管烷的平衡结构、稳定性和价键特征进行了理论预测．1…