Theoretically Calculated Deformation Density of Small Molecules

Analysis of the theoretical electron deformation density based on EHMO and ab initio calculations has been applied to the simple molecules F₂, H₂O and SO₂ The effects from varied basis sets for such deformation density were sought. The accumulation of electron density between the bonded atoms calculated from EHMO and ab initio methods with STO-3G is generally under-estimated. Such phenomena are significantly improved by using split-valence basis sets e.g. 3–21G and 4–31G. The addition of d polarization functions is apparently important for the sulfur atom in sulfur-related bonding. 3–21G or 3–21G^* basis sets were found to provide not only valuable deformation density distributions of molecules but also comparable orbital energy states with respect to the experimental values.     

Conformational stability,barriers to internal rotation,structural parameters,ab initio calculations,and vibrational assignment of cyclopropanecarboxaldehyde

The asymmetric torsional potential function, conformational energy difference, vibrational frequencies, and structural parameters of Cyclopropane-carboxaldehyde have been obtained from ab initio calculations at the 3–21G and/or 6-31G^* baiss set levels. These results have allowed for a reinterpretation or clarification of some of the corresponding results obtained from experiment. The conformations that have the oxygen atom oriented cis and trans to the three-membered ring are observed and calculated to be the most stable and high energy forms in the gaseous phase, respectively. From the ab initio calculations using the 6–31 G^* basis set, the energy difference between the two conformers is 114 cm^–1. For the liquid, the trans conformer is more stable and is the only rotamer present in the annealed solid. Based on a combination of results obtained from ab initio calculations, microwave spectroscopy, and the electron diffraction technique,r _o structural parameters have been obtained for both conformations.     

Predicting shielding constants in solution using gauge invariant atomic orbital theory and the effective fragment potential method

A method to approximate ab initio shielding constants is presented, in which the ab initio density matrix is replaced in the gauge invariant atomic orbital formalism with the density matrix resulting from an effective fragment potential calculation. The resulting first-order density matrix is then iterated to self-consistency. The method is compared with fully ab initio gauge invariant atomic orbital restricted Hartree-Fock calculations on hydrogen chloride, water, and ammonia solutes with up to nine solvent water molecules using the 6-31G, 6-31G(d,p), and 6-31+G(d,p) basis sets. Using the 6-31G(d,p) basis sets, the average of the average absolute deviations for the three environments tested is 0.34 ppm. This is sufficiently accurate to allow for the identification of specific (1)H nuclei in a solvated molecule when the chemical shift between nuclei is not less than 1 ppm. The success of the method at this level of approximation is due to a cancellation of errors between the paramagnetic and diamagnetic terms of the shielding constant: the diamagnetic term is underestimated by roughly the same amount that the paramagnetic term is overestimated.     

Locally dense basis sets for chemical shift calculations

Calculations of chemical shifts have been carried out using “locally dense” basis sets for the resonant atom of interest, and smaller, attenuated sets on other atoms in the molecule. For carbon, calculations involving a 6-311G(d) triply split valence set with polarization on the resonant atom and 3-21G atomic bases on other heavy atoms result in good agreement with experiment, and are virtually identical to those found employing the larger basis on all atoms. For species such as nitrogen, oxygen, and fluorine where standard balanced basis sets do not agree well with experiment, use of attenuated sets fail as well. The use of locally dense basis sets permits calculations previously impractical, and the successful application to carbon suggests that the chemical shift is most dependent on the local basis set, and less so on whether or not a balanced or unbalanced calculation is being carried out.     

The gas-phase molecular structure of 1,1-diethynylsilacyclobutane as determined by means of electron diffraction and ab initio calculations

As a continuation of our systematic investigation of the effect of substituents on the ring geometry and dynamics in silacyclobutanes and in order to explore the role of the silicon atom as a mediator for electronic interactions between the attached fragments, we studied the molecular structure of 1,1-diethynylsilacyclobutane (DESCB) by means of gas-phase electron diffraction and ab initio calculations. The structural refinement of the electron diffraction data yielded the following bond lengths (r_a) and bond angles (uncertainties are 3σ): r(Si–C)=1.874(2) Å, r(Si–C)=1.817(1) Å, (C–Si–C)=79.2(6)°, (C–Si–C)=106.5(6)°. The geminal Si–CC moieties were found to be bent outwards by 3.1(15)° and the puckering angle was determined to be 30.0(15)°. The evidently short Si–C bond length, which was also reproduced by the ab initio calculations, could be rationalized as being the consequence of the electronic interaction between the outer π charges of the triple bond and the 3pπ orbitals at the silicon atom. It is also likely that the conjugation of the geminal ethynyl groups leads to an enhancement of this bond contraction. Electrostatic interactions and the subsequent reduction of the covalent radius of the silicon atom may also contribute to this bond shortening. It has been found that the endocyclic Si–C bond length fits nicely within a scheme describing a monotonous decrease of the Si–C bond length with the increase of the electronegativity of the substituent in various geminally substituted silacyclobutanes.A series of related silacyclobutanes and acyclic diethynylsilanes have been studied by applying various ab initio methods and their optimized structures were compared to the structure of DESCB. Among these compounds are 1,1-dicyanosilacyclobutane (DCYSCB), which is isoelectronic to DESCB, 1,1-diethynylcyclobutane (DECB) which is isovalent to DESCB, monoethynylsilacyclobutane (MESCB) and monocyanosilacyclobutane (MCYSCB). Searching for reasonable support for the explanation of the structural results of DESCB we performed detailed natural population analysis as well as Mulliken population analysis (MPA) on DESCB and other related molecules. In contrast to the Mulliken charges, the natural atomic charges provided helpful information concerning the bonding properties in DESCB and the corresponding compounds. By varying the size of some basis sets, we could demonstrate the validity of the repeatedly discussed dependency of the Mulliken MPA on the basis set.For the performance of the quantum mechanical calculations we employed the following methods and basis sets: HF/6-31G(d,p), DFT/B3PW91/6-31G(d), DFT/B3PW91/6-311++G(d,p), MP2/6-31G(d,p) and MP2/6-311++G(d,p).     

Synthesis and vibrational study of platinum(II) and palladium(II) complexes of glyoxilic acid oxime

New platinum(II) and palladium(II) complexes of glyoxilic acid oxime (gao) have been prepared and characterised by infrared (4000–150 cm⁻¹) and Raman (4000–200 cm⁻¹) spectra. The gao acts as bidentate ligand bonding through the oxime nitrogen and carboxyl oxygen atoms to form neutral bis-chelate square-planar complexes. The lowest energy conformer of the gao ligand (ectt) was selected among 16 theoretically possible conformers on the basis of ab initio calculations at HF/3-21G*, HF/6-31G* and HF/6-311** levels of the theory from which structural parameters and conformational stabilities have been obtained. A complete vibrational assignment of the gao was performed for the lowest energy ectt conformer on the basis of ab initio optimised parameters and normal coordinate analysis calculations (PED). NCA calculations of the complexes studied were also performed.     

Scaling the Coulomb interaction in calculations of electron spectra of transition metal complexes

A simple technique of scaling two-electron integrals in ab initio calculations of the electronically excited states of transition metal complexes is proposed. This technique uses the fact that one-center two-electron integrals depend linearly on the scaling factor when Slater type functions are subjected to scaling transformation. This leads to a linear dependence of the d—d transition energy on the “scale” of Coulomb interaction, which allows one to affect the calculation result by varying the Slater exponential. To test the technique, ab initio configuration interaction and full active space calculations of the low excited states of the CrF ₆ ^3- , MnF ₆ ^2- , and VF ₆ ^3- complexes are performed. For transition elements, a basis of Slater type effective functions chosen from the optical spectra of the atoms and ions of transition elements is used. It is shown that in the STO-6G basis with effective exponentials, experimental transitions are reproduced with an accuracy of about 2000 cm^-1 even with the use of small active space determined by the orbitals localized on the central atom of the complex.     

Molecular mechanics (MM3) calculations on lithium amide compounds

The MM3 force field has been extended to deal with the lithium amide molecules that are widely used as efficient catalysts for stereoselective asymmetric synthesis. The MM3 force field parameters have been determined on the basis of the ab initio MP2/6-31G* and/or DFT (B3LYP/6-31G*, B3-PW91/6-31G*) geometry optimization calculations. To evaluate the electronic interactions specific to the lithium amides derived from the diamine molecules properly, the Lewis bonding potential term for the interaction between the lithium atom and the nonbonded adjacent electronegative atom such as nitrogen was introduced into the MM3 force field. The bond dipoles were evaluated correctly from the electronic charges on the atoms calculated by fitting to the electrostatic potential at points selected. The MM3 results on the molecular structures, conformational energies, and vibrational spectra show good agreement with those from the quantum mechanical calculations.     

Use of locally dense basis sets for nuclear magnetic resonance shielding calculations

The locally dense basis set approach to the calculation of nuclear magnetic resonance shieldings is one in which a sufficiently large or dense set of basis functions is used for an atom or molecular fragment containing the resonant nucleus or nuclei of interest and fewer or attenuated sets of basis functions employed elsewhere. Provided the dense set is of sufficient size, this approach is capable of determining chemical shieldings nearly as well as a calculation with a balanced basis set of quality equal to the locally dense set, but with considerable savings of CPU time. Detailed comparisons are provided of locally dense and balanced calculations in the gauge including atomic orbital (GIAO) method for the individual principal values, the isotropic shieldings, and the tensor orientations for hydrogen, carbon, nitrogen, oxygen, fluorine, and phosphorus nuclei. It is seen that chemical functional groups can often define the appropriate molecular fragment to be taken locally dense. While the present test cases are for the most part small molecules, the value of the method is that it will allow calculations on systems that would otherwise presently be computationally expensive or inaccessible. © John Wiley & Sons, Inc.     

Theoretical Study of the Inner/sphere Reorganization Energy and Activation Energy of Self/exchange Electron Transfer Reaction for M(H2O)26+/3+ (M=V, Cr, Mn, Fe and Co) Systems

通过建立电子转移过程的活化模型和重组模型, 提出了用量子化学从头算方法研究电子转移过程内层重组能和活化能的新方法. 在UMP26/311G水平上获得了5对过渡金属水合离子体系M(H     

Matrix isolation FT-IR, FT-Raman spectroscopy, conformational ab initio calculations, and vibrational frequencies of meso and racemic-2,4-pentanediol

The infrared spectra of meso-2,4-pentanediol and racemic-2,4-pentanediol were measured in an argon matrix at 20 K. The Raman spectra of the pure liquids (meso and racemic) were measured at room temperature. The spectra were obtained using a Fourier transform spectrophotometer and a cryostat for the low temperature matrix. The meso and racemic forms of the diol were separated by means of a spinning band distillation column. The energies of nine possible conformers of the meso form and nine conformers of the racemic form were calculated. Extensive ab initio calculations using B3LYP, MP2 and HF methods with several basis sets consistently gave the lowest energy for the TT conformer of the meso form and the G⁻T (=TG⁻) conformer of the racemic form. Ab initio calculations at the B3LYP/6-31G** level were performed for the lowest energy conformer of meso and racemic pentanediol to obtain the equilibrium geometry, vibrational frequencies, and infrared and Raman intensities. Calculated and experimental frequencies were compared to make vibrational assignments.     

Determination of ab initio polarizabilities of polymers: Application to polyethylene and polysilane

An ab initio method for calculating the longitudinal linear polarizability of polymeric chains is described. This method is equivalent to an uncoupled Hartree–Fock scheme. It is applied to polyethylene and polysilane in minimal STO-3G and extended 4-31G basis sets. The study describes important techniques for solving the difficulties met in actual calculations: band reordering of the band structures, calculation of analytical derivatives of the energy bands ?_n(k) and LCAO coefficients c_np(k), and errors caused by the improper lattice sum truncations of the Hartree–Fock matrix.     

Proton and sodium ion affinities of glycine and its sodium salt in the gas phase. Ab initio calculations

The energies of protonation and Na⁺ cationization of glycine (GLY) and its (GLY ? H + Na) salt in the gas phase were calculated using ab inltio calculations. The proton affinity of GLY, valued at the MP2/6–31G*//3-21G level, is 937 kJ mol^?1. The amino function is confirmed to be the most favourable site of protonation: ‘proton affinities’ of the carbonyl and hydroxyl functions are calculated to be 75 and 180 kJ mol^?1, respectively, lower than that of NH₂ at the MP2/6-31G*//3–21G level. Calculations performed up to the MP2/6–31G*//3–21G level give the Na⁺ affinity of GLY as 189 kJ mol^?1 and the H⁺ and Na⁺ affinities of (GLY – H + Na) as 1079 and 298 kJ mol^?1, respectively. The geometries of all neutral and protonated species optimized with the 3–21G basis set are described. Both H* and Na⁺ cations complex preferably between the nitrogen atom and the carbonyl oxygen atom, leading to pseudo-five-membered ring structures in which Na? O and Na? N bonds lengths are greater than 2 Å.     

Infrared spectra,vibrational assignment and ab initio calculations for 2,2,2-trifluoroethanimidamide

The infrared spectra of gaseous and solid 2,2,2-trifluoroethanimidamide, CF₃(NH₂)C=NH, have been recorded from 4000 to 80 cm^–1. A vibrational assignment for the normal modes is proposed based on group frequencies and normal coordinate calculations utilizing C₁ symmetry. The structures for both the cis [hydrogen atom of the =NH group is cis to the NH₂ group] and trans geometric isomers have been determined from ab initio Hartree-Fock gradient calculations employing the GAUSSIAN-82 program with the 3–21G basis set. The most stable conformer at this level of calculation is found to be a C₁, structure with a partially rotated CF₃ group and the hydrogen atom of the imine group trans to the NH₂ group. The calculated structural parameters have only very small differences between the conformers. Barriers to internal rotation for the NH₂ and CF₃ groups and vibrational frequencies have been calculated for the C₁ form. The results of this investigation are compared with similar data on some corresponding molecules.Taken in part from the thesis of T. G. Sheehan which was submitted to the Department of Chemistry in partial fulfillment of the Ph.D. degree, May 1990.     

主成分分析方法在核酸碱基量子化学计算数据处理中的应用

由于碱基在核酸中的重要性 ,多年来一直有关于碱基的理论计算报道[1～ 7] 。本文将化学计量学中的主成分分析方法[8] 用来分析五种碱基 :腺嘌呤 (Ａ)、鸟嘌呤 (Ｇ)、胞嘧啶 (Ｃ)、尿嘧啶 (Ｕ)和胸腺嘧啶 (Ｔ)计算结果的几何参数 ,以期取得有用的结构信息。1　方法通过ＡＣＤ ＣｈｅｍＳｋｅｔｃｈ 3 .5 [9] 的三维优化 (分子力学方法ＣＨＡＲＭＭ力场 )获得碱基的起始几何结构 ,其原子编号见图 1。所有的计算均采用Ｇａｕｓｓｉａｎ 94程序[10 ] 在ＩＢＭ ＰＣ兼容机上完成。首先 ,对 5种碱基作了 6种半经验方法 (ＡＭ1、ＰＭ3、ＭＮＤＯ、…     

Ab initio quantum-mechanical calculations of magnetic shielding constants of 13C, 15N and 19F in CF 3X(X = H,F,CN,CI,NO)

The magnetic shielding constants for ¹³C, ¹⁵N and ¹⁹F nuclei are calculated in the title compounds with several GIAO basis sets. Better overall agreement with experimental data is obstained with the 6-311G basis set.     

Theoretical Study of the Interplay between Lithium Bond and Hydrogen Bond in Complexes Involved with HLi and HCN

The lithium‐ and hydrogen‐bonded complex of HLi? NCH? NCH is studied with ab initio calculations. The optimized structure, vibrational frequencies, and binding energy are calculated at the MP2 level with 6‐311++G(2d,2p) basis set. The interplay between lithium bonding and hydrogen bonding in the complex is investigated with these properties. The effect of lithium bonding on the properties of hydrogen bonding is larger than that of hydrogen bonding on the properties of lithium bonding. In the trimer, the binding energies are increased by about 19 % and 61 % for the lithium and hydrogen bonds, respectively. A big cooperative energy (?5.50 kcal mol^?1) is observed in the complex. Both the charge transfer and induction effect due to the electrostatic interaction are responsible for the cooperativity in the trimer. The effect of HCN chain length on the lithium bonding has been considered. The natural bond orbital and atoms in molecules analyses indicate that the electrostatic force plays a main role in the lithium bonding. A many‐body interaction analysis has also been performed for HLi? (NCH)_N (N=2–5) systems.     

Ab initio molecular orbital study on the gas phase SN2 reaction F + CH3Cl → CH3F + Cl

Ab-initio molecular orbital (MO) and direct ab initio dynamics calculations have been applied to the gas phase S_N2 reaction F⁻ + CH₃Cl → CH₃F + Cl⁻. Several basis sets were examined in order to select the most convenient and best fitted basis set to that of high-quality calculations. The Hartree–Fock (HF) 3−21+G(d) calculation reasonably represents a potential energy surface calculated at the MP2/6−311++G(2df,2pd) level. A direct ab initio dynamics calculation at the HF/3−21+G(d) level was carried out for the S_N2 reaction. A full dimensional ab initio potential energy surface including all degrees of freedom was used in the dynamics calculation. Total energies and gradients were calculated at each time step. Two initial configurations at time zero were examined in the direct dynamics calculations: one is a near collinear collision, and the other is a side-attack collision. It was found that in the near collinear collision almost all total available energy is partitioned into two modes: the relative translational mode between the products (40%) and the C − F stretching mode (60%). The other internal modes of CH₃F were still in the ground state. The lifetimes of the early- and late-complexes F⁻ … CH₃Cl and FCH₃ … Cl⁻ are significantly short enough to dissociate directly to the products. On the other hand, in the side-attack collision, the relative translation energy was about 20% of total available energy.     

Study of hydrogen bonds in N-methylacetamide by DFT calculations of oxygen, nitrogen, and hydrogen solid-state NMR parameters

Solid-state nuclear magnetic resonance (NMR) parameters of ¹⁷O, ¹⁴N/¹⁵N, and ²H/¹H nuclei were evaluated in two available neutron crystalline structures of N-methylacetamide (NMA) at 250 and 276 K, NMA-I and NMA-II, respectively. Density functional theory calculations were performed by B3LYP method and 6-311++G** and IGLO-II type basis sets to calculate the electric field gradient (EFG) and chemical shielding (CS) tensors at the sites of mentioned nuclei. In order to investigate hydrogen bonds (HBs) effects on NMR tensors, calculations were performed on four-model systems of NMA: an optimized isolated gas-phase, crystalline monomers, crystalline dimers, and crystalline trimers. Comparing the calculated results reveal the influence of N–H···O=C and C–H···O=C HB types on the NMR tensors which are observable by the evaluated parameters including quadrupole coupling constant, C _Q, and isotropic CS, σ _iso. Furthermore, the results demonstrate more influence of HB on the EFG and CS tensors of NMA at 276 K rather than that of 250 K.