Ab initio and analytic intermolecular potentials for Ar-CH(3)OH

Ab initio calculations at the CCSD(T)/aug-cc-pVTZ level of theory were used to characterize the Ar-CH(3)OH intermolecular potential energy surface (PES). Potential energy curves were calculated for four different Ar + CH(3)OH orientations and used to derive an analytic function for the intermolecular PES. A sum of Ar-C, Ar-O, Ar-H(C), and Ar-H(O) two-body potentials gives an excellent fit to these potential energy curves up to 100 kcal mol(-1), and adding an additional r(-n) term to the Buckingham two-body potential results in only a minor improvement in the fit. Three Ar-CH(3)OH van der Waals minima were found from the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVTZ calculations. The structure of the global minimum is in overall good agreement with experiment (X.-C. Tan, L. Sun and R. L. Kuczkowski, J. Mol. Spectrosc., 1995, 171, 248). It is T-shaped with the hydroxyl H-atom syn with respect to Ar. Extrapolated to the complete basis set (CBS) limit, the global minimum has a well depth of 0.72 kcal mol(-1) with basis set superposition error (BSSE) correction. The aug-cc-pVTZ basis set gives a well depth only 0.10 kcal mol(-1) smaller than this value. The well depths of the other two minima are within 0.16 kcal mol(-1) of the global minimum. The analytic Ar-CH(3)OH intermolecular potential also identifies these three minima as the only van der Waals minima and the structures predicted by the analytic potential are similar to the ab initio structures. The analytic potential identifies the same global minimum and the predicted well depths for the minima are within 0.05 kcal mol(-1) of the ab initio values. Combining this Ar-CH(3)OH intermolecular potential with a potential for a OH-terminated alkylthiolate self-assembled monolayer surface (i.e., HO-SAM) provides a potential to model Ar + HO-SAM collisions.     

The effect of electron correlation on the conformational space of melatonin

The pineal gland hormone melatonin regulates several physiological processes including circadian rhythm and also alleviates oxidative stress‐induced degenerative diseases. In spite of its important biological roles, no high level ab initio conformational study has been conducted to reveal its structural features. In this work, the conformational flexibility of melatonin was investigated using correlated ab initio calculations. Conformers, obtained previously at the Hartree‐Fock level (HF/6‐31G*), were fully optimized using second order Møller‐Plesset perturbation theory applying the frozen core approximation (MP2(FC)/6‐31G*). Furthermore, single‐point MP4(SDQ,FC)/6‐31G*//MP2(FC)/6‐31G* computations were performed to investigate the effect of higher order perturbation terms. The HF and MP2 conformational spaces are considerably different: the initial 128 structures converged into 102 different local minima as confirmed by frequency calculations; 28 new minima appeared and 26 previous HF local minima disappeared; no “all‐trans” C3 side chain conformations are seen at the MP2(FC) level. The MP2 global minimum conformation is stabilized by an aromatic‐side chain interaction. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2008     

Representation of molecules by atomic charges: A new population analysis

For the calculation of atomic charges in molecules, a new fast procedure based on extended Mulliken population analysis is presented. The reliability of the new population analysis is tested within the AM1 approximation and the results are compared in detail with those obtained by different methods reported in the literature.     

Semispectroscopic and quantitative structure-property relationship estimates of the equilibrium and vibrationally averaged structure and dipole moment of 1-buten-3-yne

Systematic quantum chemical calculations have been performed to obtain precise estimates of the equilibrium and vibrationally averaged molecular structure and electric dipole moment of vinylacetylene (VA, 1-buten-3-yne). Anharmonic (cubic and semi-diagonal quartic) MP2/cc-pVTZ force fields in normal coordinates were computed to account for anharmonic vibrational effects, including zero-point contributions to the rotational constants and the electric dipole moment. A simultaneous weighted least-squares structural refinement was performed, resulting in the best semispectroscopic estimate of the re structure of VA. The refinement was based on experimentally measured ground-state rotational constants of two isotopologs of VA corrected to equilibrium values using MP2/cc-pVTZ vibration-rotation interaction constants and all-electron CCSD(T)/aug-cc-pVTZ structural constraints. The semispectroscopic re structure of VA agrees excellently with the high-level CCSD(T)/aug-cc-pVTZ ab initio structure. The most dependable, CCSD(T)/cc-pVQZ//CCSD(T)/aug-cc-pVTZ equilibrium electric dipole moment of VA, in D, is mua= 0.4088, mub= 0.0004, and muc= 0. The vibrationally corrected a-component of 0.4214 D is in excellent agreement with one of the available experimental values. The present analysis shows that mub is negligible even after vibrational correction. A simple quantitative structure-property relationship (QSPR) model resulted in a highly similar estimate, 0.45 D, for the electric dipole moment of VA.     

Vector algebra and molecular symmetry: a tribute to Professor Josiah Willard Gibbs

Vector algebra, as developed by Josiah Willard Gibbs, is much simpler than matrix or tensor algebra, therefore, it is more suitable to introduce the students of chemistry into the wonderful world of molecular symmetry. A program based on elementary vector algebra has been written to determine all symmetry elements and symmetry operations of rigid molecular structures. The program also contains data for 57 point groups common in chemistry. Therefore, it automatically supplies the particular point group to which the structure belongs. Since the locations of the nuclei related to the symmetry elements are also revealed by the program, even the detailed notation of the framework group of the molecular structure can be deduced. The program can be a great help in determining the symmetries of the normal modes of vibration, too.     

优化固氮活性模型信息研究—双齿配体DPPE和DPPM对固氮酶促反应的影响

报道了某些双齿配体对固氮酶促反应影响的研究结果．发现１，２－双（二苯基磷）－乙烷（ＤＰＰＥ）对固氮酶酶促乙炔还原反应有促进作用，但同时又抑制了固氮酶的放氢反应。而比ＤＰＰＥ少一个－ＣＨ２－链的双（二苯基磷）－甲烷（ＤＰＰＭ）却不能表现出对固氮酶促乙炔还原活力的促进作用．对照固氮酶ＭｏＦｅ－蛋白Ｘ－光衍射结构分析结果和量子化学近似计算所导出的固氮酶活性中心结构模型提出了ＤＰＰＥ促进酶促乙炔还原反应的一种可能的解释．     

About the generalized LM-inverse and the weighted Moore-Penrose inverse

The recursive method for computing the generalized LM-inverse of a constant rectangular matrix augmented by a column vector is proposed in Udwadia and Phohomsiri (2007) [16] and [17]. The corresponding algorithm for the sequential determination of the generalized LM-inverse is established in the present paper. We prove that the introduced algorithm for computing the generalized LM-inverse and the algorithm for the computation of the weighted Moore-Penrose inverse developed by Wang and Chen (1986) in [23] are equivalent algorithms. Both of the algorithms are implemented in the present paper using the package MATHEMATICA. Several rational test matrices and randomly generated constant matrices are tested and the CPU time is compared and discussed.