The standard thermodynamic functions of ethyl-, formyl-, and benzoylferrocene in the ideal gas state

The standard enthalpies of formation of gaseous ethylferrocene, formylferrocene, and benzoylferrocene were calculated at the B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d,p) level of density functional theory using the method of isodesmic reactions. The absolute entropy, heat capacity, and enthalpy were calculated for ethylferrocene and formylferrocene by statistical thermodynamics methods over the temperature range 100–1000 K. The molecular constants necessary for calculations (structural parameters, vibrational frequencies, and internal rotation potential) were estimated at the B3LYP/6-31G(d,p) level. Empirical estimates were used to check the reliability and mutual consistency of literature experimental data on the enthalpies of formation and entropies of several ferrocene derivatives. Unreliable data were revealed and more reliable estimates of their errors were obtained.     

A study of factors affecting the accurac of combined calculation methods

By the example of combined calculation methods (CCMs) corresponding to calculations in MP4/6–311+G(fd,p)//MP2/6-31G(d,p) and MP4/aug-cc-pvTZ//MP2/cc-pvDZ approximations the factors affecting their accuracy are considered. By total energy calculations made for a number of compounds it is shown that the accuracy of CCMs reduces when valence double- and triple-split basis sets are applied together, and also when different methods are used to allow for electron correlation. The use of empirical corrections taking into account the character of the electron distribution in a molecule allows an enhancement of the accuracy of CCMs. The mentioned factors taken into account enable us to obtain CCM, for which the mean absolute deviation of calculation data is 1.0 kJ/mol and the largest maximum deviation is 8.0 kJ/mol in the total energy calculation in the MP4/6-311+G(f,d,p) approximation.     

A quantum chemical study of 1,2,4-triazine reactivity in reactions with electrophilic and nucleophilic reagents

B3LYP/6-311+G(_d,p) and MP2/6-31G(_d,p) methods are used to study the electronic structure of a 1,2,4-triazine molecule. Quantitative characteristics of the reactivity are obtained in the form of atomic charges and Fukui coefficients. Thermodynamic characteristics are determined for the protonation reaction and addition of the hydride anion to a 1,2,4-triazine molecule.     

Protolytic properties of cyanic and thiocyanic acids and their isoforms

The electronic structure of HOCN, HSCN, HNCO, and HNCS molecules and [OCN]^? and [SCN]^? anions has been studied by ab initio calculations at HF/6-31G(d), HF/6-31G(d, p), MP2/6-31G(d)//HF/6-31G(d), and MP2/6-31G(d, p)//HF/6-31G(d, p) levels of theory. The HNCO and HNCS molecules are shown to have higher thermodynamic stability than HOCN and HSCN, respectively. The protolyte strength series are substantiated: HSCN > HOCN, HNCS > HNCO, HOCN > HNCO, HSCN > HNCS. Computations including electron correlation [MP2/6-31G(d)//HF/6-31G(d) and MP2/6-31G(d, p)//HF/6-31G (d, p)] reproduce the general sequence of proton-donor properties: HSCN > HOCN > HNCS > HNCO, which coincides with the hydrophobicity series for the compounds. The relative proton-donor capacity of these acids in water solutions is generally governed by the electronic structure and by the size of their molecules and [OCN]^? and [SCN]^? anions, but not by medium effects.     

DFT investigation of the mechanism of CH2CO + O(3P) reaction

A detailed theoretical survey of the potential energy surface (PES) for the CH₂CO + O(³P) reaction is carried out at the QCISD(T)/6‐311+G(3df,2p)//B3LYP/6‐311+G(d,p) level. The geometries, vibrational frequencies, and energies of all stationary points involved in the reaction are calculated at the B3LYP/6‐311+G(d,p) level. More accurate energy information is provided by single‐point calculations at the QCISD(T)/6‐311+G(3df,2p) level. Relationships of the reactants, transition states, intermediates, and products are confirmed by the intrinsic reaction coordinate (IRC) calculations. The results suggest that P1(CH₂+CO₂) is the most important product. This study presents highlights of the mechanism of the title reaction. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Quantum-Chemical Study on Positional Selectivity in the Trimethylsilylation and Sulfonation of Pyrrole and <Emphasis Type="Italic">N</Emphasis>-Alkylpyrroles

Positional selectivity (α:β ratio) of electrophilic substitution in pyrrole, N-methylpyrrole, and N-tert-butylpyrrole was analyzed by ab initio [RHF/6-31G(d), MP2/6-31G(d)//RHF/6-31G(d)] and DFT [B3LYP/6-31G(d)] calculations of some σ-complexes derived from the substrates. The results of calculations with the use as model electrophilic species of trimethylsilyl cation [MP2/6-31G(d)//RHF/6-31G(d) and B3LYP/6-31G(d)] and SO₃ molecule [B3LYP/6-31G(d)] instead of proton are fairly consistent with the experimental data, according to which trimethylsilylation of pyrrole and its N-substituted derivatives with trimethylsilyl trifluoromethanesulfonate, as well as sulfonation with pyridine-sulfur trioxide complex, gives the corresponding β-substituted products.     

Blue‐shifted and red‐shifted hydrogen bonds: Theoretical study of the CH3CHO· · ·HNO complexes

The blue‐shifted and red‐shifted H‐bonds have been studied in complexes CH₃CHO…HNO. At the MP2/6‐31G(d), MP2/6‐31+G(d,p) MP2/6‐311++G(d,p), B3LYP/6‐31G(d), B3LYP/6‐31+G(d,p) and B3LYP/6‐311++G(d,p) levels, the geometric structures and vibrational frequencies of complexes CH₃CHO…HNO are calculated by both standard and CP‐corrected methods, respectively. Complex A exhibits simultaneously red‐shifted C? H…O and blue‐shifted N? H…O H‐bonds. Complex B possesses simultaneously two blue‐shifted H‐bonds: C? H…O and N? H…O. From NBO analysis, it becomes evident that the red‐shifted C? H…O H‐bond can be explained on the basis of the two opposite effects: hyperconjugation and rehybridization. The blue‐shifted C? H…O H‐bond is a result of conjunct C? H bond strengthening effects of the hyperconjugation and the rehybridization due to existence of the significant electron density redistribution effect. For the blue‐shifted N? H…O H‐bonds, the hyperconjugation is inhibited due to existence of the electron density redistribution effect. The large blue shift of the N? H stretching frequency is observed because the rehybridization dominates the hyperconjugation. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Theoretical study on the insertion reactions of the germylenoid H2GeClMgCl with RH (R = F, OH, NH2)

The insertion reactions of the germylenoid H₂GeClMgCl with RH (R = F, OH, NH₂) have been studied by using the DFT B3LYP and QCISD methods. The geometries of the stationary points on the potential energy surfaces of the reactions were optimized at the B3LYP/6-311+G(d, p) level of theory. The calculated results indicate that all the mechanisms of the three insertion reactions are identical to each other. The QCISD/6-311++G(d, p)//B3LYP/6-311+G(d, p) calculated potential energy barriers for the three insertion reactions of R = F, OH, and NH₂ are 164.62, 193.30, and 200.73 kJ mol^?1, and the reaction energies for the three reactions are ?57.46, ?35.65, and ?22.22 kJ mol^?1, respectively. Under the same situation, the insertion reactions should occur easily in the following order H-F > H-OH > H-NH₂. In THF solvent the insertion reactions get more difficult than in gas phase.     

Proton affinities of ketones, vicinal diketones and α-keto esters: a computational study

The proton affinities of seven different ketones, vicinal diketones, and α-keto esters (acetophenone, 2,2,2-trifluoroacetophenone, 2,3-butanedione, 1-phenyl-1,2-propanedione, methyl pyruvate, ethyl benzoylformate, and ketopantolactone) have been evaluated theoretically using the conventional ab initio HF and several post-HF methods (MP2, MP4, CCSD), density functional methods with the B3LYP hybrid functional, as well as some ab initio model chemistries [CBS-4M, G2(MP2), and G3(MP2)//B3LYP]. The chemical compounds studied are frequently used substrates in the asymmetric hydrogenation over chirally modified platinum catalysts where the protonation properties of the chiral modifier and the substrates are of great interest. In most cases, the proton affinities (PAs) evaluated with the CCSD/6-311+G(d,p)//B3LYP/TZVP and G2(MP2) methods are in good agreement with the existing experimental ones. However, the previously reported PA of 2,3-butanedione seems to be too high by 10-15 kJ mol⁻¹. The B3LYP/TZVP//B3LYP/TZVP and MP2/6-311+G(d,p)//B3LYP/TZVP model chemistries predict proton affinities that are systematically higher and lower than the experimental PAs, respectively. If proton affinities are evaluated as the average of the PAs calculated with these two theoretical methods a very good agreement with the experimental results is obtained. The mean absolute deviation (MAD) from experiment of this combination method for the PAs of 13 test molecules is 4.0 kJ mol⁻¹. For 9 molecules composed only of first-row atoms the MAD is 2.5 kJ mol⁻¹. The B3LYP/TZVP//B3LYP/TZVP and MP2/6-311+G(d,p)//B3LYP/TZVP methods provide significant savings in computational time and disk space compared to the CCSD/6-311+G(d,p)//B3LYP/TZVP and G2(MP2) models. Therefore, it is suggested that if no experimental or highly accurate theoretical data is available (due to computational cost), the proton affinities of similar compounds as investigated in this paper, can be evaluated with the combination method. For the studied molecules, this method gives the following PAs (in kJ mol⁻¹): 788 (2,3-butanedione, exptl 802); 798 (2,2,2-trifluoroacetophenone, exptl 799); 811 (ketopantolactone); 813 (methyl pyruvate); 825 (1-phenyl-1,2-propanedione); 862 (acetophenone, exptl 861); 865 (ethyl benzoylformate).     

Why betaine crystallizes in high local Cs symmetry. An ab initio MO and DFT study of anhydrous betaine and betaine monohydrate

A theoretical study of the structure, charge distribution, rotational barrier and fundamental vibrations of anhydrous betaine (CH₃)₃NCH₂COO (trimethylglycine) was carried out and compared with available experimental data. Calculations were carried out at HF, MP2 and B3LYP levels using a 6-31+G(d,p) basis set. The calculated rotational barrier of the betaine carboxylic group is 40.5 kJ/mol at the MP4(SDQ)/6-311G(d,p)//HF/6-31+G(d,p) level of theory. The rotation of the carboxylic group changes the molecule from a highly symmetric (C _s ) conformation into a twisted conformation resulting in shortening of the molecule by about 50 pm. Natural population analysis (NPA) indicates intramolecular interaction between the carboxylic oxygen and the nearest methyl hydrogens resulting in internal hydrogen bonding. MP4(SDQ)/6-311G(d,p) single-point NPA calculations on a betaine monohydrate model taken from the X-ray geometry show an expected weakening in the internal hydrogen bond. Calculations explain why betaine preferentially crystallizes in high local C _s symmetry. Received: 24 March 1998 / Accepted: 3 September 1998 / Published online: 7 December 1998     

Theoretical study in gas phase of linear and nonlinear optical properties of the ortho-, meta- and para-nitrophenol isomers

The linear (α), and nonlinear (β, γ) optical NLO properties of ortho-, meta- and para-nitrophenol (ONP, MNP and PNP) isomers have been calculated in gas phase by using ab initio (HF, MP2 and MP4) and density functional theory (DFT) (B3LYP, CAM-B3LYP) methods, with the 6-31+G(d,p) and 6-311+G(3d,3p) standard and the Sadlej specialized basis sets. These properties were evaluated both at static and at dynamic regime within the finite field FF numerical techniques and the time-dependent-Hartree–Fock approach at 1,910 nm, respectively. Additional calculations were performed for the β static hyperpolarizability of these isomers in presence of p-dioxane solvent with the Onsager Model and the SCRF-PCM approach, using B3LYP/6-31+G(d,p) and MP2/6-31+G(d,p) levels of theory. Additionally, CCSD/6-31+G(d,p) calculations were performed for the α, β and γ properties of PNP isomer. The B3LYP and MP2 α _ave results of the nitrophenol isomers are comparable to the experimental α _ave reports; while the tendency for the β _v calculated values (β _v PNP > β _v MNP > β _v ONP), that can be explained in terms of the O _x atomic charge of the –NO₂ group, does not follow exactly the experimental ones. The B3LYP γ _ave results are in correspondence to the experimental measurements, the correlation of which is r ² = 0.99. The use of FF methodology in conjunction with the B3LYP and MP2 methods and the 6-31+G(d,p) basis set show to be appropriate approaches to predict qualitative optical properties of Push–Pull like organic molecules, provided are considered the solvent effects or frequency dependence. However, to have a clear picture of the NLO properties of an isolated molecule, higher order correlation effects combined with specialized basis sets, frequency and solvent effects should be employed. We have demonstrated that MP4/Sadlej level of theory is able to reproduce NLO properties that can be considered equivalent to those from more sophisticated approaches, such as CCSD together with extended basis sets.     

Proton affinities of molecules containing nitrogen and oxygen: comparing ab initio and semiempirical results to experiments

We report a comparison of theoretical and experimental proton affinities at nitrogen and oxygen sites within a series of small molecules. The calculated proton affinities are determined using the semiempirical methods AM 1, MNDO , and PM 3; the ab initio Hartree–Fock method at the following basis levels: 3-21G //3-21G , 3-21+G //3-21G , 6-31G *//6-31G *, and 6-31+G (d, p)//6-31G *; and Møller–Plesset perturbation calculations: MP 2/6-31G *//6-31G *, MP 3/6-31G *//6-31G *, MP 2/6-31G +(d, p)//6-31G *, MP 3/6-31G +(d, p)//6-31G *, and MP 4(SDTQ )/6-31G +G (d, p)//6-31G *. The semiempirical methods have more nonsystematic scatter from the experimental values, compared to even the minimal 3-21G level ab initio calculations. The thermodynamically corrected 6-31G *//6-31G * proton affinities provide acceptable results compared to experiment, and we see no significant improvement over 6-31G *//6-31G * in the proton affinities with any of the higher-level calculations. © 1992 John Wiley & Sons, Inc.     

Structure and conformational stability of CH2CHCH2X (X=F, Cl and Br) molecules—post Hartree–Fock and density functional theory methods

The molecular structure and conformational stability of CH₂CHCH₂X (X=F, Cl and Br) molecules were studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 3-fluoropropene were optimized employing BLYP and B3LYP levels of theory of DFT method implementing 6-311+G(d,p) basis set. The MP2/6-31G^*, BLYP and B3LYP levels of theory of ab initio and DFT methods were used to optimize the 3-chloropropene and 3-bromopropene molecules. The structural and physical parameters of the molecules are discussed with the available experimental values. The rotational potential energy surface of the above molecules were obtained at MP2/6-31G^* and B3LYP/6-311+G(d,p) levels of theory. The Fourier decomposition of the rotational potentials were analyzed. The HF/6-31G^* and MP2/6-31G^* levels of theory have predicted the cis conformer as the minimum energy structure for 3-fluoropropene, which is in agreement with the experimental values, whereas the BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory reverses the order of conformation. The ΔE values calculated for 3-chloropropene at MP2/6-31G^*, BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory show that the gauche form is more stable than the cis form, which is in agreement with the experimental value. The same levels of theory have also predicted that the gauche form is stable than cis for 3-bromopropene molecule. The maximum hardness principle has been able to predict the stable conformer of 3-fluoropropene at HF/6-31G^* level of theory, but the same level of theory reverses the conformational stability of 3-chloropropene and 3-bromopropene molecules and MP2/6-31G^* level of theory predicted the stable conformer correctly.     

Theoretical investigation on the GaCl3-catalyzed ring-closing metathesis reaction of N-2,3-butadienyl-2-propynyl-1-amine: three-membered ring versus four-membered ring mechanism

The gallium chloride (GaCl(3))-catalyzed ring-closing metathesis reaction mechanism of N-2,3-butadienyl-2-propynyl-1-amine has been studied at the Becke three-parameter hybrid functional combined with Lee-Yang-Parr correlation functional (B3LYP)/6-31G(d), B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p)//B3LYP/ 6-31G(d) and the second-order M?ller-Plesset perturbation (MP2)/6-311++G(d,p)//B3LYP/6-31+G(d,p) levels. It was found that the final metathesis product can be yielded via a three-membered or four-membered ring mechanism. The three-membered ring pathway is favorable due to its low energy barrier at the rate determining step. The whole reaction is stepwise and strongly exothermic.     

Theoretical study of the conformational energy hypersurface of cyclotrisarcosyl

The multidimensional Conformational Potential Energy Hypersurface (PEHS) of cyclotrisarcosyl was comprehensively investigated at the DFT (B3LYP/6-31G(d), B3LYP/6-31G(d,p) and B3LYP/6-311++G(d,p)), levels of theory. The equilibrium structures, their relative stability, and the Transition State (TS) structures involved in the conformational interconversion pathways were analyzed. Aug-cc-pVTZ//B3LYP/6-311++G(d,p) and MP2/6-31G(d)//B3LYP/6-311++G(d,p) single point calculations predict a symmetric cis-cis-cis crown conformation as the energetically preferred form for this compound, which is in agreement with the experimental data. The conformational interconversion between the global minimum and the twist form requires 20.88 kcal mol-1 at the MP2/6-31G(d)//B3LYP/6-311++G(d,p) level of theory. Our results allow us to form a concise idea about the internal intricacies of the PEHSs of this cyclic tripeptide, describing the conformations as well as the conformational interconversion processes in this hypersurface. In addition, a comparative analysis between the conformational behaviors of cyclotrisarcosyl with that previously reported for cyclotriglycine was carried out     

A quantum-chemical study of nitroguanidine and tetramethyl nitroguanidine: Stable structures,force fields,and the interpretation of the vibrational spectrum

Stable structures of nitroguanidine and tetramethyl nitroguanidine were searched by the MP2 and B3LYP methods with the use of the 6–31G(d,p) and 6–311G(d,p) basis sets. For nitroguanidine, the force fields calculated by these methods were scaled using the experimental spectra for seven isotopomers. The obtained scaling factors made it possible to theoretically evaluate and assign the frequencies in the spectrum of tetramethyl nitroguanidine.     

丁烯和丁烷自由基阳离子的分子结构和超精细结构的 理论研究

用密度函数理论B3LYP方法和6-31G(d,p),6-311G(d,p)及6-311+G(d,p)基组，分别对1-C4H^+~8,2-C4H^+~8和C4H^+~10进行了构型优化和频率分析计算，预言1-C4H^+~8具有非平面构型，与以往报道的从头算和密度函数理论计算结果不同。在各自由基阳离子的B3LYP构型上，进行了B3LYP、MP2及MRSDCI方法的超精细偶合常数计算，得到了比以往更好的结果，特别是MP2/B3LYP计算值是至今与实验值符合得最好的理论计算结果。     

A new scheme for determining the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of glycine and alanine peptides

In this paper a new scheme was proposed to calculate the intramolecular hydrogen-bonding energies in peptides and was applied to calculate the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of the glycine and alanine peptides. The density-functional theory B3LYP6-31G(d) and B3LYP6-311G(d,p) methods and the second-order Moller-Plesset perturbation theory MP26-31G(d) method were used to calculate the optimal geometries and frequencies of glycine and alanine peptides and related structures. MP26-311++G(d,p), MP26-311++G(3df,2p), and MP2/aug-cc-pVTZ methods were then used to evaluate the single-point energies. It was found that the B3LYP6-31G(d), MP26-31G(d), and B3LYP6-311G(d,p) methods yield almost similar structural parameters for the conformers of the glycine and alanine dipeptides. MP2/aug-cc-pVTZ predicts that the intramolecular seven-membered ring N-H...O=C hydrogen-bonding strength has a value of 5.54 kcal/mol in glycine dipeptide and 5.73 and 5.19 kcal/mol in alanine dipeptides, while the steric repulsive interactions of the seven-membered ring conformers are 4.13 kcal/mol in glycine dipeptide and 6.62 and 3.71 kcal/mol in alanine dipeptides. It was also found that MP26-311++G(3df,2p) gives as accurate intramolecular N-H...O=C hydrogen-bonding energies and steric repulsive interactions as the much more costly MP2/aug-cc-pVTZ does.     

已烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p),6-31+G(d,p),6-311G(d,p)及6- 311+G(d,p)基组,分别对1-C_6H_(12)~+,2-C_6H_(12)~+和3-C_6H_(12)~+的各种构 象进行了几何构型优化,并在B3LYP/6-311G(d,p)水平上进行了频率分析计算,在 各优化构型上,使用B3LYP和MP2(full)方法进行了超精细结构的计算。计算的3- C_6H_(12)~+的超精细偶合常数比以往的计算结果更好;1-C_6H_(12)~+和2-C_6H_ (12)~+的超精细偶合常数目前尚无实验数据报道,本计算预言了它们的超精细偶合 常数和最稳定构型。     

Density functional theory insight into Eu(III) and Am(III) complexes with two 2,6-dicarboxypyridine diamide-type ligands

Extraction complexes of Eu(III) and Am(III) with two 2,6-dicarboxypyridine diamide-type ligands L–A and L–B (Fig. 1) are studied by density functional theory (DFT). At both B3LYP/6-31G(d)/RECP and MP2/6-31G(d)/RECP levels of theory, the geometrical optimizations of the structures of the complexes can achieve the same accuracy and obtain the same geometrical configuration. At the B3LYP/6-311G(d,p)/RECP level of theory Eu³⁺ and Am³⁺ prefer to form [ML]³⁺ complexes under the solvation conditions, and the Am(III) complexes with L–A are more stable than the corresponding Eu(III) complexes. In the system with the ligand L–B, both [ML]³⁺ and [ML(NO₃)₃] species are very unstable.