Why is the antipodal effect in closo-1-SB9H9 so large? A possible explanation based on the geometry from the concerted use of gas electron diffraction and computational methods

The molecular structure of 1-thia-closo-decaborane(9), 1-SB(9)H(9), has been determined by the concerted use of gas electron diffraction and quantum-chemical calculations. Assuming C(4v) symmetry, the cage structure was distorted from a symmetrically bicapped square antiprism (D(4d) symmetry) mainly through substantial expansion of the tetragonal belt of boron atoms adjacent to sulfur. The S-B and (B-B)(mean) distances are well determined with r(h1) = 193.86(14) and 182.14(8) pm, respectively. Geometrical parameters calculated using the MP2(full)/6-311++G** method and at levels reported earlier [MP2(full)/6-311G**, B3LYP/6-311G** and B3LYP/cc-pVQZ], as well as calculated vibrational amplitudes and (11)B NMR chemical shifts, are in good agreement with the experimental findings. In particular, the so-called antipodal chemical shift of apical B(10) (71.8 ppm) is reproduced well by the GIAO-MP2 calculations and its large magnitude is schematically accounted for, as is the analogous antipodal chemical shift of B(12) in the twelve-vertex closo-1-SB(11)H(11).     

The C7-C10 cycloalkanes revisited

The conformers of cycloheptane through cyclodecane have been examined at the B3LYP/6-311+G* and MP2/6-311+G* theoretical levels, with some additional calculations at the CCD/6-311+G* and CCSD(T)/6-311++G** levels. With cyclooctane, B3LYP predicts that the boat-chair and crown conformers have similar energies, whereas MP2 and CCSD(T) predict that the crown conformer is 2 kcal/mol higher in energy. The latter is in agreement with the electron diffraction data. With cyclononane, B3LYP predicts that two of the higher-energy conformers found in molecular mechanics calculations should convert to one of the lower-energy conformers. However, MP2/6-311+G* optimizations find them to be true minima on the potential energy surface. B3LYP systematically predicts larger C-C-C bond angles for these compounds than either MP2 or CCD. The results of molecular mechanics MM4 calculations are generally in good agreement with those obtained using MP2.     

Three rotor potential energy scans, conformational equilibrium constants and vibrational analysis of 3-fluoro-1-propanol CH(2)FCH(2)CH(2)OH

The conformational stability and the three rotor internal rotations in 3-fluoro-1-propanol were investigated by the DFT-B3LYP/6-311+G** and the ab initio MP2/6-311+G** levels of theory. The calculated potential energy curves of the molecule at both levels of theory were consistent with complex conformational equilibria of about 12 minima, all of which were predicted to have real frequencies at both the B3LYP and the MP2 levels. The lowest energy minimum in the potential curves of 3-fluoro-1-propanol was predicted to correspond to the Gauche-gauche-trans (Ggt) conformer in excellent agreement with microwave and electron diffraction results. The equilibrium constants for the conformational interconversion of the molecule were calculated and found to correspond to an equilibrium mixture of about 33% Ggt, 14% Ggg1 and 13% Gg1g and about 43% Ggt, 12% Ggg1 and 10% Gg1g distribution by the B3LYP/6-311+G** and the MP2/6-311+G** calculations, respectively, at 298.15K. The vibrational frequencies of each molecule in its three stable forms were computed at B3LYP level and complete vibrational assignments were made based on normal coordinate calculations and comparison with experimental data of the molecule.     

Structural studies on ethyl isovalerate by microwave spectroscopy and quantum chemical calculations

We observed the microwave spectrum of ethyl isovalerate by molecular beam Fourier transform microwave spectroscopy. The rotational and centrifugal distortion constants of the most abundant conformer were determined. Its structure was investigated by comparison of the experimental rotational constants with those obtained by ab initio methods. In a first step, the rotational constants of various conformers were calculated at the MP2/6-311++G** level of theory. Surprisingly, no agreement with the experimental results was found. Therefore, we concluded that in the case of ethyl isovalerate more advanced quantum chemical methods are required to obtain a reliable molecular geometry. Ab initio calculations carried out at MP3/6-311++G**, MP4/6-311++G**, and CCSD/6-311++G** levels and also density functional theory calculations using the B3LYP/6-311++G** method gave similar results for the rotational constants, but they were clearly distinct from those obtained at the MP2/6-311++G** level. With use of these more advanced methods, the rotational constants of the lowest energy conformer were in good agreement with those obtained from the microwave spectrum.     

Theoretical study of light-metal tetrahydroaluminates

Ab initio calculations of the structural, energetic, vibrational, and magnetic characteristics of the lowest-lying structures for isolated molecules and ions of light-metal tetrahydroaluminates (Li, Na, Be, Mg, and Al) have been performed by the perturbation theory (MP2), quadratic configuration interaction (QCISD(T)), coupled cluster (CCSD(T)), and density functional theory (B3LYP) methods using the 6-31G*, 6-31G**, 6-311+G**, and 6-311++G** basis sets. The trends in the behavior of the molecular characteristics have been analyzed in various related series of these compounds. The results are compared with the data on analogous light-metal tetrahydroborates calculated at the same levels of theory. The differences in structure and stability between analogous hydroborate and alanate complexes are examined. The economical approximation B3LYP/6-311++G**//B3LYP/6-311+G** has been shown to adequately reproduce the results obtained at the higher level of theory CCSD(T)/6-311++G**//MP2/6-31G* even though it requires considerably shorter CPU times and smaller amounts of memory.     

Density functional and experimental studies on the FT-IR and FT-Raman spectra and structure of 2,6-diamino purine and 6-methoxy purine

FT-IR and FT-Raman spectra of 2,6-diamino purine (DAP) and 6-methoxy purine (MP) have been recorded in the regions of 4000-400cm(-1) and 3500-100cm(-1), respectively. The spectra were interpreted with the aid of normal coordinate analysis following full structure optimizations and force field calculations based on density functional theory (DFT) using standard B3LYP/6-31G* and B3LYP/6-311+G** method and basis set combinations. Normal coordinate calculations were performed with the DFT force field corrected by a recommended set of scaling factors yielding fairly good agreement between observed and calculated frequencies.     

D-苯丙氨酸与Cu＋(1S0, 3d10)气相反应理论研究

用量子化学密度泛函(DFT)方法研究D-苯丙氨酸与一价基态金属阳离子Cu^＋在气相中反应的机理. 在B3LYP/6-31G*水平上, 优化了反应包含的4个反应通道的反应物、中间体、过渡态和产物的几何构型, 并采用B3LYP/DZVP, B3LYP/[6-311＋G**(C,H,O)＋Lanl2dz(Cu)], B3LYP/6-311＋G**, MP2/6-311＋G** 等方法对各驻点进行了单点能计算. 通过对计算结果的分析, 获得了其单重态反应势能面的一般轮廓、各驻点几何构型优化参数, 明确了其反应机理.     

Synthesis, characterization, crystal structure and DFT studies on 3-(1H-benzo[d][1,2,3]triazol-1-yl)-1-oxo-1-m-tolylpropan-2-yl-nicotinate

3-(1H-Benzo[d][1,2,3]triazol-1-yl)-1-oxo-1-m-tolylpropan-2-yl-nicotinate (BOTN) has been synthesized and characterized by elemental analysis, IR, UV-vis and fluorescence spectroscopy. Its crystal structure has also been determined by X-ray single crystal diffraction. For BOTN, density functional theory (DFT) calculations of the structure and vibrational frequencies have been performed at B3LYP/6-311G** level. The comparisons between the experimental vibrational frequencies and the predicted data show that B3LYP/6-311G** method can simulate the IR of BOTN on the whole. Based on the vibration analysis, thermodynamic properties of BOTN have been calculated. The correlative equations between the thermodynamic properties and the temperatures have also been listed. The experimental UV-vis spectra present two peaks and theoretical UV-vis spectra obtained by TD-DFT method exhibit three peaks. The comparison between them suggests that the B3LYP/6-311G** method can only approximately simulate the UV-vis spectra of BOTN. The fluorescence determination reveals two emission bands at 423 and 489 nm, respectively.     

Correlated ab initio quantum chemical calculations of di- and trisaccharide conformations

High level correlated quantum chemical calculations, using MP2 and local MP2 theory, have been performed for conformations of the disaccharide, beta-maltose, and the trisaccharide, 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranose. For beta-maltose, MP2 and local MP2 calculations using the 6-311++G** basis set are in good agreement, predicting a global minimum gas-phase conformation with a counterclockwise hydrogen bond network and the experimentally-observed intersaccharide hydrogen bonding arrangement. For conformations of 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranose, MP2/6-311++G**, and local MP2/6-311++G** calculations do not provide a consensus prediction of relative energetics, with the MP2 method finding large differences in stability between extended and folded trisaccharide conformations. Local MP2 calculations, less susceptible to intramolecular basis set superposition errors, predict a narrower range of trisaccharide energetics, in line with estimates from Hartree-Fock theory and B3LYP and BP86 density functionals. All levels of theory predict compact, highly hydrogen-bonded conformations as lowest in energy on the in vacuo potential energy surface of the trisaccharide. These high level, correlated local MP2/6-311++G** calculations of di- and trisaccharide energetics constitute potential reference data in the development and testing of improved empirical and semiempirical potentials for modeling of carbohydrates in the condensed phase.     

4-乙酰基-1,3-二氟苯硝化反应机理研究

以4-氨基-1,3-二氟苯、乙酸酐和浓硝酸等为原料,合成了4-乙酰基-5-硝基-1,3-二氟苯,并且采用FTIR和1H NMR对产物进行了表征.采用量子化学方法,探讨了在酸催化条件下4-乙酰基-1,3-二氟苯硝化反应过程中的微观机理,在密度泛函理论的B3LYP/6-311G**水平下搜索得到反应路径,并对各中间体、过渡态和产物进行几何优化、振动频率分析和标准热力学函数计算;反应路径经过IRC验证.在MP2/6-311G**水平上对所有中间体、过渡态和产物进行更精确的能量计算.     

On the aromatic character of the heterocyclic bases of DNA and RNA

Studies based on ab initio optimized geometries (at B3LYP/6-311+G** and MP2/6-311+G** levels) and on experimental structures retrieved from the Cambridge Structural Database (CSD) reveal that the nucleobases constituting DNA and RNA differ significantly in their aromatic character, as shown by the geometry-based index of aromaticity HOMA that ranges from 0.466 for thymine to 0.917 for adenine, based on B3LYP/6-311+G** calculations, and 0.495-0.926, respectively, if based on the MP2/6-311+G** level. Aromaticity of the bases decreases markedly with an increase of the number of double-bond C=X (X = N, O) substituents at the rings. H-bonds involving C=O groups in Watson-Crick pairs cause an increase of the aromatic character of the rings.     

Analysis of vibrational spectra of 3-halo-1-propanols CH(2)XCH(2)CH(2)OH (X is Cl and Br)

The conformational stability and the three rotor internal rotations in 3-chloro- and 3-bromo-1-propanols were investigated by DFT-B3LYP/6-311+G and ab initio MP2/6-311+G, MP3/6-311+G and MP4(SDTQ)//MP3/6-311+G levels of theory. On the calculated potential energy surface twelve distinct minima were located all of which were not predicted to have imaginary frequencies at the B3LYP level of theory. The calculated lowest energy minimum in the potential curves of both molecules was predicted to correspond to the Gauche-gauche-trans (Ggt) conformer in excellent agreement with earlier microwave and electron diffraction results. The equilibrium constants for the conformational interconversion of the two 3-halo-1-propanols were calculated at the B3LYP/6-311+G level of calculation and found to correspond to an equilibrium mixture of about 32% Ggt, 18% Ggg1, 13% Tgt, 8% Tgg and 8% Gtt conformations for 3-chloro-1-propanol and 34% Ggt, 15% Tgt, 13% Ggg1, 9% Tgg and 7% Gtt conformations for 3-bromo-1-propanol at 298.15K. The nature of the high energy conformations was verified by carrying out solvent experiments using formamide ( epsilon=109.5) and MP3 and MP4//MP3 calculations. The vibrational frequencies of each molecule in its three most stable forms were computed at the B3LYP level and complete vibrational assignments were made based on normal coordinate calculations and comparison with experimental data of the molecules.     

An exploration of conformational search of leucine molecule and their vibrational spectra in gas phase using ab initio methods

An extensive exploration of the conformational space has been carried out to characterize all possible gas phase structures of leucine. A total of 324 unique trial structures for canonical leucine were generated by considering all possible combinations of single bond rotamers. All trial structures were optimized at the B3LYP/6-311G* level of the DFT method. A total of 77 unique and stationary canonical conformers were found. Further, 15 most stable conformers were reoptimized at B3LYP/6-311++G** level and their respective relative energies, vertical ionization energies, hydrogen bonding patterns, rotational constants and dipole moments were calculated. A single point energy calculations for leucine conformers have also been done at both B3LYP/6-311++G(2df, p) and MP2/6-311++G(2df, p) levels. The good agreement between our estimates of rotational constants for two most stable conformers and available experimental measurements supports the reliability of the B3LYP/6-311++G** level of theory for describing the conformational behavior of leucine molecule. The proton affinity and gas phase basicity were also determined. Using the statistical approach, conformational distributions at various temperatures have also been performed and analyzed. Vibrational spectra were also calculated. It is also observed that zwitterions of leucine are not stable in gas phase.     

A study of internal rotations and vibrational spectra of oxiranemethanol (glycidol)

The conformational stability and the C–O and O–H internal rotations in oxiranemethanol were investigated at the DFT-B3LYP/6-311G**, MP2/6-311G** and MP4(SDQ)/6-311G** levels of theory. Three minima were predicted in the CCOH potential energy scans of the molecule to have relative energies of about 2 kcal/mol or less and all were calculated to have real frequencies upon full optimization of structural parameters at the DFT and the MP2 levels of calculations. The Cg1 (H bond inner) conformation was predicted to be the lowest energy conformation for oxiranemethanol in excellent agreement with an earlier microwave study. The equilibrium mixture was calculated from Gibb's free-energy changes to be about 79% Cg1, 17% G1g and 3% Gg1 at the B3LYP/6-311G** level and about 87% Cg1, 11% G1g and 2% Gg1 at the MP2/6-311G** level for oxiranemethanol at 298.15 K. No conclusive evidence was obtained for the presence of high-energy form in the liquid phase of oxiranemethanol. The vibrational frequencies of oxiranemethanol in its three stable forms were computed at the B3LYP level and complete vibrational assignments were made for the lowest energy Cg1 form on basis of calculated and experimental data of the molecule.     

Modeling zinc in biomolecules with the self consistent charge-density functional tight binding (SCC-DFTB) method: applications to structural and energetic analysis

Parameters for the zinc ion have been developed in the self-consistent charge density functional tight-binding (SCC-DFTB) framework. The approach was tested against B3LYP calculations for a range of systems, including small molecules that contain the typical coordination environment of zinc in biological systems (cysteine, histidine, glutamic/aspartic acids, and water) and active site models for a number of enzymes such as alcohol dehydrogenase, carbonic anhydrase, and aminopeptidase. The SCC-DFTB approach reproduces structural and energetic properties rather reliably (e.g., total and relative ligand binding energies and deprotonation energies of ligands and barriers for zinc-assisted proton transfers), as compared with B3LYP/6-311+G** or MP2/6-311+G** calculations.     

Structures of water octamers (H2O)8: exploration on ab initio potential energy surfaces by the scaled hypersphere search method

The potential energy surface (PES) of water octamers has been explored by the scaled hypersphere search method. Among 164 minima on the PES (based on MP2/6-311++G(3df,2p)//B3LYP/6-311+G(d,p) calculations), the cubic structure with D2d symmetry has been confirmed to be the global minimum. In a thermodynamic simulation using these 164 structures, the cubic structure with S4 symmetry has the highest population at low temperature, though double rings can become dominant as temperature going up, in good accord with a recent Monte Carlo simulation using an empirical potential. A transition temperature from cubic to noncubic has significantly been underestimated when potential energy data of B3LYP/6-311+G(d,p) calculations are employed in the simulation. This serious discrepancy between the MP2 and the B3LYP results suggests an importance of dispersion interactions for discussions on thermodynamics of water octamers.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要.在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型,B3LYP方法得到了三种构型(复合物Ⅰ,Ⅱ和Ⅲ),而MP2方法只能得到一种构犁(复合物Ⅱ).在复合物Ⅰ和Ⅲ中,HSO单元中的1H原子作为质子供体.与O3分子中的端基O原子作为质子受体相互作用,形成红移氢键复合物;而在复合物Ⅱ中,虽与复合物Ⅰ和Ⅲ中具有相间的质子供体和质子受体,却形成了蓝移氢键复合物.B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重甍误差(BSSE)和零点振动能(ZPVE)校正,其值在-3.37到-4.55 kJ·mol-1之间.采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查,并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

Theoretical study of light-metal tetrahydroborates

Ab initio calculations of the structural, energetic, vibrational, and magnetic characteristics of the lowest-lying structures for isolated molecules and ions of light-metal tetrahydroborates (Li, Na, Be, Mg, and Al) have been performed by the perturbation theory (MP2), quadratic configuration interaction (QCISD(T)), coupled cluster (CCSD(T)), and density functional theory (B3LYP) methods using the 6-31G*, 6-31G**, 6-311+G**, and 6-311++G** basis sets. The trends in the behavior of the structural parameters, the energies of different decomposition pathways, barriers to internal rotation of BH₄ groups, normal mode frequencies, magnetic shielding constants, and spin density distribution (in radicals) have been analyzed in various related series of these compounds. The results obtained by ab initio methods and at the DFT level are compared. The economical approximation B3LYP/6-311++G**//B3LYP/6-311+G** adequately reproduces the results obtained at the higher level of theory CCSD(T)/6-311++G*s*//MP2/6-31G* even though it requires considerably shorter CPU times and smaller amounts of memory.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要. 在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型, B3LYP方法得到了三种构型(复合物I, II和III), 而MP2方法只能得到一种构型(复合物II). 在复合物I和III中, HSO单元中的1H原子作为质子供体, 与O3分子中的端基O原子作为质子受体相互作用, 形成红移氢键复合物; 而在复合物II中, 虽与复合物I和III中具有相同的质子供体和质子受体, 却形成了蓝移氢键复合物. B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重叠误差(BSSE)和零点振动能(ZPVE)校正, 其值在-3.37到-4.55 kJ·mol-1之间. 采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查, 并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

Interaction of anions with perfluoro aromatic compounds

The complexes formed by a variety of anions with perfluoro derivatives of benzene, naphthalene, pyridine, thiophene, and furan have been calculated using DFT (B3LYP/6-31++G**) and MP2 (MP2/6-31++G** and MP2/6-311++G**) ab initio methods. The minimum structures show the anion interacting with the pi-cloud of the aromatic compounds. The interaction energies obtained range between -8 and -19 kcal mol(-1). The results obtained at the MP2/6-31++G** and MP2/6-311++G** levels are similar. However, the B3LYP/6-31++G** results provide longer interaction distances and smaller interaction energies than do the MP2 results. The interaction energies have been partitioned using an electrostatic, polarization, and van der Waals scheme. The AIM analysis of the electron density shows a variety of topologies depending on the aromatic system considered.