The infrared fundamental intensities and polar tensor of CH3NC

The molecular force field and polar tensor of methyl isocyanide have been determined from its gas phase vibrational frequencies and infrared intensities. Quantum chemical results from MP2(FC), B3LYP and quadratic configuration interaction calculation including single and double substitutions procedures using a 6-311 + +G(3d,3p) basis set have been used to determine the signs of the dipole moment derivatives with respect to the normal coordinates as well as estimate individual fundamental intensities of the overlapped v1-v5 and v3-v6 band systems. Principal component graphical representations of the A1 and E symmetry polar tensor elements were useful in determining preferred sets of tensor elements. The mean dipole moment derivative (GAPT charge) of the methyl carbon in CH3NC, 0.347 e, is between the corresponding values in CH3CN, 0.110 e, and CH3F, 0.541 e. The mean dipole moment derivatives obtained here indicate the correct 1s methyl carbon ionization energy as 293.35 eV which is 0.98 eV higher than the corresponding ionization energy of the terminal atom.     

The infrared vibrational intensities and polar tensors of HFCO and DFCO

The infrared vibrational intensities of HFCO and DFCO have been calculated at the B3LYP/cc-pVTZ, MP2(FC)/6-311++G(3d,3p) and QCISD(FU)/aug-cc-pVTZ levels. All calculations predict the isotopomers to have identical intensity sums, within about 1 km mol(-1). This is in contrast with experimental intensity sum results reported in the literature. Dipole moment derivative directions calculated by the three methods are in excellent agreement for all in-plane normal coordinates. All the theoretical polar tensor elements are also in good agreement with each other having standard deviations varying between 0.003 and 0.043 e. The oxygen and fluorine atoms have negative mean derivatives (approximately -0.6 e), whereas the carbon mean derivative is very positive (approximately +1.1e) and the hydrogen one is almost zero (approximately +0.03 e). The HFCO theoretical intensity sums calculated by all three methods as well as their carbon and oxygen mean dipole moment derivatives are in good agreement with those estimated from the experimental intensities and atomic mean derivatives of H2CO and F2CO.     

An atomic charge-charge flux-dipole flux atom-in-molecule decomposition for molecular dipole-moment derivatives and infrared fundamental intensities

The molecular dipole moment and its derivatives are determined from atomic charges, atomic dipoles, and their fluxes obtained from AIM formalism and calculated at the MP2(FC)/6-311++G(3d,3p) level for 16 molecules: 6 diatomic hydrides, CO, HCN, OCS, CO2, CS2, C2H2, C2N2, H2O, H2CO, and CH4. Root-mean-square (rms) errors of 0.052 D and 0.019 e are found for the dipole moments and their derivatives calculated using AIM parameters when compared with those obtained directly from the MP2(FC)/6-311++G(3d,3p) calculations and 0.097 D and 0.049 e when compared to the experimental values. The major deviations occur for the NaH, HF, and H2O molecules. Parallel polar tensor elements for the diatomic and linear polyatomic molecules, except H2, HF, LiH, and NaH, have values resulting from cancellations of substantial contributions from atomic charge fluxes and atomic dipole fluxes. These fluxes have a large negative correlation coefficient, -0.97. IR fundamental intensity sums for CO, HCN, OCS, CO2, CS2, C2H2, C2N2, H2CO, and CH4 calculated using AIM charges, charge fluxes, and atomic dipole fluxes have rms errors of 14.9 km mol(-1) when compared with sums calculated directly from the molecular wave function and 36.2 km mol(-1) relative to experimental values. The classical model proposed here to calculate dipole-moment derivatives is compared with the charge-charge flux-overlap model long used by spectroscopists for interpreting IR vibrational intensities. The utility of the AIM atomic charges and dipoles was illustrated by calculating the forces exerted on molecules by a charged particle. AIM quantities were able to reproduce forces due to a +0.1 e particle over a 3-8-A separation range for the CO and HF molecules in collinear and perpendicular arrangements. These results show that IR intensities do contain information relevant to the study of intermolecular interactions.     

Protonated nitro group: structure, energy and conjugation

Structure of protonated nitro compounds was investigated by calculations at the levels MP2(FC)/6-311++G(2d,2p)//MP2(FC)/6-311++G(2d,2p)(nitromethane and reference compounds) or B3LYP/6-311+G(d,p)//B3LYP/6-311+G(d,p)(nitrobenzene and its 18 meta- and para-substituted derivatives). The group NO2H+ reveals many similarities with the isoelectronic group CO2H as the preferred conformation, conformational equilibrium, and stabilization by interaction (resonance) within the group quantified by means of isodesmic reactions. However, there is a difference in the interaction with donor groups (for instance in 4-nitroaniline) that is much stronger with NO2H+ than with CO2H. This interaction may be called resonance and may be described by standard resonance formulas, but these formulas predict only partially the geometry and cannot explain the great interaction energy.     

The infrared fundamental intensities and polar tensor of CF4

Atomic polar tensors of carbon tetrafluoride are calculated from experimental fundamental infrared intensities measured by several research groups. Quantum chemical calculations using a 6-311 + + G(3d, 3p) basis set at the Hartree-Fock, M?ller-Plesset 2 and Density Functional Theory (B3LYP) levels are used to resolve the sign ambiguities of the dipole moment derivatives. The resulting carbon mean dipole moment derivative, pC = 2.051 e, is in excellent agreement with values estimated by a MP2/6-311 + + G(3d, 3p) theoretical calculation, 2.040 e, and by an empirical electronegativity model, 2.016 e. The pC value determined here is also in excellent agreement with the one obtained from the CF4 1s carbon ionization energy using a simple potential model, 2.059 e. Crawford's G intensity sum rule applied to the fundamental intensities of CH4, CH3F, CH2F2 and CHF3 results in a prediction of a 1249 km mol(-1) intensity sum for CF4 in good agreement with the experimental values of 1328 +/- 37.9, 1208.0 +/- 54.4 and 1194.8 +/- 7.4 km mol(-1) reported in the literature.     

QTAIM charge-charge flux-dipole flux models for the infrared fundamental intensities of difluoro- and dichloroethylenes

A quantum theory of atoms in molecules (QTAIM) charge-charge flux-dipole flux (CCFDF) decomposition of the MP2/6-311++G(3d,3p) level molecular dipole moment derivatives is reported for the cis-, trans-, and 1,1-difluoroethylenes and the cis- and trans-dichloroethylenes. Although the dipole moment derivatives and infrared fundamental intensities calculated at the MP2 level are overestimated for high-intensity bands corresponding to CF and CC stretching vibrations, the overall agreement is good with a root-mean-square (rms) error of 19.6 km mol-1 for intensities ranging from 0 to 217.7 km mol-1. The intensities calculated from the QTAIM/CCFDF model parameters are in excellent agreement with those calculated directly by the MP2/6-311++G(3d,3p) approach with only a 1.8 km mol-1 rms error. A high negative correlation (r=-0.91) is found between the charge flux and dipole flux contributions to the dipole moment derivatives. Characteristic values of charge, charge flux, and dipole flux contributions are found for CF, CCl, and CH stretching derivatives. The CH stretching derivatives provide especially interesting results with very high charge flux and dipole flux contributions with opposite signs. The charge, charge flux, and dipole flux contributions are found to be transferable from the cis to the trans isomers providing accurate predictions of the theoretical trans intensities with rms errors of 8.6 km mol-1 for trans-difluoroethylene and 5.9 km mol-1 for trans-dichloroethylene.     

Molecular Dipole Moment Computed with Ab Initio MKS Charges

1 INTRODUCTION Dipole moment (DM) can be measured by a variety of experimental methods[1] or computed with atomic charge distribution directly derived from molecular orbital calculations[2~10]. The average errors in DM computed with AM1, PM3 and MNDO methods for 125 organic compounds are 0.35 D, 0.45 D and 0.45 D, respectively[2]. Feller et al. reported that the average errors in the computed DM at the HF level are 0.65 D, 0.49 D and 0.30 D with the basis set STO-3G, 3-21G and…     

分子几何构型优化方法的系统性比较

对《ＣＲＣ物理化学手册》第７７版中收集的第三周期以前的所有已知实验构型的无机分子的构型,以ＭＰ２、Ｂ３ＬＹＰ、Ｂ３ＰＷ９１级别上进行了构型优化的系统性比较,优化采用基组为６－３１Ｇ（ｄ,ｐ）,６－３１１Ｇ（ｄ,ｐ）,６－３１１Ｇ（２ｄ,ｐ）。对大多数分子另比较了ＱＣＩＳＤ（Ｔ）方法,对多原子分子比较了ＢＰＷ９１方法,对含氢双原子分子也使用ＱＣＩＳＤ方法。结果表明,键长的平均绝对偏差（单位：ｐｐｍ）     

Theoretical study on the hydrogen abstraction reaction of CH3CH2F (HFC-161) with Cl atom

A direct dynamics method is employed to study the hydrogen abstraction reaction of CH₃CH₂F+Cl. Three distinct transition states are located, one for -H abstraction and two for β-H abstraction. The potential energy surface (PES) information is obtained at the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) and G2//MP2/6-311G(d,p) level. Based on the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) results, the rate constants of the three reaction channels are evaluated by using the canonical variational transition state theory (CVT) with small-curvature tunneling (SCT) contributions over the temperature range of 220–2800 K. The calculated results indicate that -H abstraction dominates the total reaction almost over the whole temperature range.     

Stacking interactions as the principal design element in acyl-transfer catalysts

The conformational properties and the stability of acylpyridinium intermediates formed in pyridine-catalyzed acylation reactions have been studied at the SCS-MP2(FC)/6-311 +G(d,p)//MP2(FC)/ 6-31G(d) level of theory. It has been shown that stacking interactions can play a decisive role in the stability as well as the conformational preferences of these transient intermediates.     

Modern valence-bond description of chemical reaction mechanisms: the 1,3-dipolar addition of diazomethane to ethene

The electronic mechanism for the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane (CH2N2) to ethene (C2H4) is described through spin-coupled (SC) calculations at a sequence of geometries along the intrinsic reaction coordinate obtained at the MP2/6-31G(d) level of theory. It is shown that the bonding rearrangements occurring during the course of this reaction follow a heterolytic pattern, characterized by the movement of three well-identifiable orbital pairs, which are initially responsible for the pi bond in ethene and the C-N pi bond and one of the N-N pi bonds in diazomethane and are retained throughout the entire reaction path from reactants to product. Taken together with our previous SC study of the electronic mechanism of the 1,3-dipolar cycloaddition of fulminic acid (HCNO) to ethyne (C2H2) (Theor. Chim. Acc. 1998, 100, 222), the results of the present work suggest strongly that most gas-phase concerted 1,3-dipolar cycloaddition reactions can be expected to follow a heterolytic mechanism of this type, which does not involve an aromatic transition state. The more conventional aspects of the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane to ethene, including optimized transition structure geometry, electronic activation energy, activation barrier corrected for zero-point energies, standard enthalpy, entropy and Gibbs free energy of activation, have been calculated at the HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), MP2/6-31G(d,p), QCISD/6-31G(d) and CCD/6-31G(d) levels of theory. We also report the CCD/6-311++G(2d, 2p)//CCD/6-31G(d), MP4(SDTQ)/6-311++G(2d,2p)//CCD/6-31G(d) and CCSD(T)/6-311++G(2d, 2p)//CCD/6-31G(d) electronic activation energies.     

Thermochemical properties of molecules and cations of MgNnHm (n=1,2 and m=1–6): enthalpies of formation, proton affinities and ionization energies

Ab initio molecular orbital calculations are reported for small neutral molecules and cations containing magnesium, nitrogen and hydrogen. Structures have been optimized using gradient techniques at B3LYP/6-31+G(d) and at MP2(full)/6-311++G(d,p). Single-point calculations are reported at QCISD(T)(full)/6-311++G(2df,p) and at CCSD(T)(full)/6-311++G(2df,p) levels using geometries optimized at MP2(full)/6-311++G(d,p). Standard enthalpies of formation at 298 K have been calculated at these two higher levels of theory. Other thermochemical properties calculated include ionization energies and proton affinities. The binding enthalpies of ammonia to Mg⁺, MgNH₂⁺ and MgNH₃⁺ are also reported.     

QTAIM charge-charge flux-dipole flux interpretation of electronegativity and potential models of the fluorochloromethane mean dipole moment derivatives

Infrared fundamental vibrational intensities and quantum theory atoms in molecules (QTAIM) charge-charge flux-dipole flux (CCFDF) contributions to the polar tensors of the fluorochloromethanes have been calculated at the QCISD/cc-pVTZ level. A root-mean-square error of 20.0 km mol(-1) has been found compared to an experimental error estimate of 14.4 and 21.1 km mol(-1) for MP2/6-311++G(3d,3p) results. The errors in the QCISD polar tensor elements and mean dipole moment derivatives are 0.059 e when compared with the experimental values. Both theoretical levels provide results showing that the dynamical charge and dipole fluxes provide significant contributions to the mean dipole moment derivatives and tend to be of opposite signs canceling one another. Although the experimental mean dipole moment derivative values suggest that all the fluorochloromethane molecules have electronic structures consistent with a simple electronegativity model with transferable atomic charges for their terminal atoms, the QTAIM/CCFDF models confirm this only for the fluoromethanes. Whereas the fluorine atom does not suffer a saturation effect in its capacity to drain electronic charge from carbon atoms that are attached to other fluorine and chlorine atoms, the zero flux electronic charge of the chlorine atom depends on the number and kind of the other substituent atoms. Both the QTAIM carbon charges (r = 0.990) and mean dipole moment derivatives (r = 0.996) are found to obey Siegbahn's potential model for carbon 1s electron ionization energies at the QCISD/cc-pVTZ level. The latter is a consequence of the carbon mean derivatives obeying the electronegativity model and not necessarily to their similarities with atomic charges. Atomic dipole contributions to the neighboring atom electrostatic potentials of the fluorochloromethanes are found to be of comparable size to the atomic charge contributions and increase the accuracy of Siegbahn's model for the QTAIM charge model results. Substitution effects of the hydrogen, fluorine, and chlorine atoms on the charge and dipole flux QTAIM contributions are found to be additive for the mean dipole derivatives of the fluorochloromethanes.     

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.     

Ab initio and density functional calculations on the pericyclic vs pseudopericyclic mode of conjugated nitrile ylide 1, 5-electrocyclizations

Ab initio (MP2/6-311+G and MP4(SDTQ)/6-311+G//MP2/6-311+G) and density functional (B3LYP/6-311+G) calculations on the ring closure reactions of conjugated nitrile ylides 1a-e, 3, and 6 to the corresponding oxazoles 2a, 5, 7, and 8; thiazoles 2b and 4; imidazole 2c; and pyrroles 2d and 2e, respectively, are reported. Vinyl nitrile ylides 1d and 1e cyclize with a substantially higher barrier than nitrile ylides containing a heteroatom. Geometric features as well as electronic structures as obtained by NBO analysis are indicative of a pericyclic, monorotatory 1, 5-electrocyclization of 1d and 1e. For nitrile ylides where X = heteroatom, a pseudopericyclic heteroelectrocyclization pathway, characterized by in-plane attack of the heteroatom's lone pair at the nitrile ylide group, is found. For 3 and 6, where two different cyclization products are possible, the calculated barriers and reaction energies are in line with the experimentally observed direction of reaction. Vinyl nitrile ylides 1d and 1e are characterized by an allene, acyl substituted derivatives 1a, 1b, 3, and 6 by a propargyl type structure. The nitrogen derivative 1c represents an intermediate case.     

Methyl and trifluoromethyl vinyl sulfoxides. Steric and electronic structure

Potential functions of internal rotation about the $C_{sp^2 } $ -S bond in H₂C=CHSOCY₃ (Y = H or F) were determined and stationary points were identified by vibrational analysis at the MP2(full)/6-31+G(d), B3PW91/6-31+G(d), and B3PW91/6-311+G(3df,p) levels. Energetically favorable conformations were established, and rotation barriers and molecular geometry parameters were evaluated. Wave functions [MP2(full)/6-31+G(d)] were analyzed by the NBO method. Energies and dipole moments of bond and lone-pair orbitals, principal types of donor-acceptor interactions, bond orders, and atomic charges were determined.     

A theoretical study of the H-abstraction reactions from HOI by moist air radiolytic products (H, OH, and O (3P)) and iodine atoms (2P(3/2))

The rate constants of the reactions of HOI molecules with H, OH, O ((3)P), and I ((2)P(3/2)) atoms have been estimated over the temperature range 300-2500 K using four different levels of theory. Geometry optimizations and vibrational frequency calculations are performed using MP2 methods combined with two basis sets (cc-pVTZ and 6-311G(d,p)). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the cc-pVTZ, cc-pVQZ, 6-311+G(3df,2p), and 6-311++G(3df,3pd) basis sets. Reaction enthalpies at 0 K were calculated at the CCSD(T)/cc-pVnZ//MP2/cc-pVTZ (n = T and Q), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), and CCSD(T)/6-311++G(3df,3pd)//MP2/6-311G(d,p) levels of theory and compared to the experimental values taken from the literature. Canonical transition-state theory with an Eckart tunneling correction is used to predict the rate constants as a function of temperature. The computational procedure has been used to predict rate constants for H-abstraction elementary reactions because there are actually no literature data to which the calculated rate constants can be directly compared. The final objective is to implement kinetics of gaseous reactions in the ASTEC (accident source term evaluation code) program to improve speciation of fission products, which can be transported along the reactor coolant system (RCS) of a pressurized water reactor (PWR) in the case of a severe accident.     

几种含芯电子相关能修正的G2和G2(ACI)方法

针对冻结芯电子近似,在MP2/6-311G(d,p)级别上对G2、MP2/6-311G(d)和MP2/6-311G(d,p)级别上对G2(QCI)方法进一步考虑了芯电子相关能修正,尝试建立了G2(fu2)、G2(QCI/ful)和G2(QCI/fu2)方法。G2-l test set 的反应能量计算结果表明,这些方法进一步减小了经验修正量;G2(QCI/ful)和G2(QCI/fu2)也比G2(QCI)的总体精度有所提高;但G2(fu2)在G2基础上,总体精度没有改善。G2(fu2)、G2(QCI/ful)和G2(QCI/fu2)计算G2-l test set反应能量的平均绝对偏差分别为5.11、4.74和4.81kJ mol-1,G2和G2(QCI)分别为5.09和4.97kJ mol-1.     

N(4S)+CH3X(X=Cl、Br)反应机理的理论研究

The reaction of N(4S)+CH3X(X=Cl、Br) was studied by the ab initio method. The geometries of the reactants, transition states and products were optimized at the MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single-point calculations for all the stationary points were carried out at the MP2/6-311++G(3df,2p) and the QCISD(T)/6-311+G(d,p) levels using the MP2/6-311+G(d,p) optimized geometries. The energies of all the stationary points were calculated by the G2MP2 method. The results of this theoretical study indicate that the reaction has three reaction channels: H abstraction reaction channel a, Cl or Br abstraction reaction channel b and substitution reaction channel c. For the N(4S)+CH3Cl reaction, reaction channel a is the main reaction channel. Reaction channels b and c may have a slight contribution in the reaction. For the N(4S)+CH3Br reaction, reaction channel a is the main reaction channel. Reaction channels b and c may have some contribution in the reaction.     

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.