QTAIM charge-charge flux-dipole flux models for the infrared fundamental intensities of the fluorochloromethanes

The molecular dipole moments, their derivatives, and the fundamental IR intensities of the fluoro-, chloro-, and fluorochloromethanes are determined from QTAIM atomic charges and dipoles and their fluxes at the MP2/6-311++G(3d,3p) level. Root-mean-square (rms) errors of 0.01 D and 5.6 km mol(-1) are found for the dipole moments and fundamental IR intensities calculated using QTAIM parameters when compared with those obtained directly from the MP2/6-311++(3d,3p) calculations and 0.04 D and 23.1 km mol(-1) when compared with the experimental values. Charge, charge flux, and dipole flux contributions are calculated for all the normal vibrations of these molecules. A large negative correlation coefficient of -0.92 is calculated between the charge flux and dipole flux contributions and indicates that electron transfer from one side of the molecule to the other during vibrations is accompanied by relaxation with electron density polarization in the opposite direction. The CF, CCl, and CH stretching normal modes of these molecules are shown to have characteristic sets of charge, charge flux, and dipole flux contributions. Although the FCF and ClCCl deformation normal modes can also be discriminated from one another based on the sizes and signs of these contributions, some HCH deformations have contributions that are similar to those for some of the ClCCl deformations.     

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.     

Ab initio studies of H2PXYH molecules (X,Y = O,S)

Ab initio molecular orbital theory is applied to the study of P? O and P? S bonding in the hypervalent phosphinic (H₂POOH), phosphinothioic (H₂POSH), and phosphinodithioic (H₂PSSH) acid molecules. Intramolecular proton exchange reactions are followed using the intrinsic reaction coordinate and Self-Consistent-Field energy localized orbitals. The P? O and PS bonds are characterized via force constants, phosphorus d orbital populations, and localized orbitals and are best described as either normal single bonds or dative bonds augmented by π back donation from the oxygen or sulfur lone pairs. The anions of these acids are also investigated, and they are found to contain only dative bonds to sulfur and oxygen.     

CH2X(X=H,F,Cl)与臭氧反应机理的量子化学研究

用量子化学UMP2方法,在6-311++G**基组水平上研究了CH2X(X=H,F,Cl)与臭氧反应机理,全参数优化了反应过程中反应物、中间体、过渡态和产物的几何构型,在UQCISD(T)/6-311++G**水平上计算了它们的能量,并对它们进行了振动分析,以确定中间体和过渡态的真实性.从CH2X(X=H,F,Cl)与O3的反应机理的研究结果看,它们与O3反应的活性都比较强,相对而言,活性大小顺序为CH2F＞CH3＞CH2Cl,也就是说,CH2F自由基与臭氧间的反应活性最强,对大气臭氧的损耗将是最大的.同时研究还发现CH2X(X=H,F,Cl)系列自由基与O3的反应都是强放热反应.     

C_4H_4Y(Y=O,S,Se)与BX_3(X=H,F,Cl)作用的计算化学研究

采用量子化学的密度泛函B3LYP和二阶微扰MP2(full)方法对C4H4Y(Y=O,S,Se)与BX3(X=H,F,Cl)形成的电子授受型复合物进行了研究,所得18个复合物的构型包括BX3位于C=C双键上方的π-p作用型和B与O,S,Se直接作用的n-p作用型.体系C4H4Y-BH3以n-p作用型较为稳定,体系C4H4Y-BF3,C4H4Y-BCl3的π-p和n-p作用型复合物稳定性相当.对各复合物的几何构型、振动频率和自然键轨道分析表明,复合物的形成过程中均存在几何构型的改变、电荷的转移和振动频率的变化,它们的变化规律与复合物稳定性的变化规律基本一致,即按H,F,Cl的顺序依次降低.     

Comparative Computational Study of Hydrogen Abstraction Reactions of CY3H + XO− (X,Y = F,Cl, and Br)

The effects of halogen substituents on the reactivity are characterized by the hybrid B3LYP and BHandHLYP functionals of density functional theory using the aug‐cc‐pVDZ basis set. The species XO⁻ and CY₃H, where X, Y = F, Cl, and Br, have been chosen as model reactants in this work. Also, the mechanism of the hydrogen abstraction (HAT) reaction has been used to study the chemical reactivity of these anionic reactions. Our theoretical findings suggest that the relative reactivity of the CY₃H + XO⁻ reactions increases as Y goes from F to Br and decreases as X goes from F to Br. Moreover, among all reactions investigated in this study, the special role of the Y has very dominant effect on activation of the C–H bond in CY₃H when XO⁻ attacks the CY₃H. Again, through the transition state theory the rate constants at 298–1000 K are also evaluated for the HAT reactions, indicating the lower the temperature the faster is the chemical reaction.     

Halogen-hydride interaction between Z-X (Z = CN, NC; X = F, Cl, Br) and H-Mg-Y (Y = H, F, Cl, Br, CH3)

Halogen-hydride interactions between Z-X (Z = CN, NC and X = F, Cl, Br) as halogen donor and H-Mg-Y (Y = H, F, Cl, Br, CH(3)) as electron donor have been investigated through the use of Becke three-parameter hybrid exchange with Lee-Yang-Parr correlation (B3LYP), second-order M?ller-Plesset perturbation theory (MP2), and coupled-cluster single and double excitation (with triple excitations) [CCSD(T)] approaches. Geometry changes during the halogen-hydride interaction are accompanied by a mutual polarization of both partners with some charge transfer occurring from the electron donor subunit. Interaction energies computed at MP2 level vary from -1.23 to -2.99 kJ/mol for Z-F···H-Mg-Y complexes, indicating that the fluorine interactions are relatively very weak but not negligible. Instead, for chlorine- and bromine-containing complexes the interaction energies span from -5.78 to a maximum of -26.42 kJ/mol, which intimate that the interactions are comparable to conventional hydrogen bonding. Moreover, the calculated interaction energy was found to increase in magnitude with increasing positive electrostatic potential on the extension of Z-X bond. Analysis of geometric, vibrational frequency shift and the interaction energies indicates that, depending on the halogen, CN-X···H interactions are about 1.3-2.0 times stronger than NC-X···H interactions in which the halogen bonds to carbon. We also identified a clear dependence of the halogen-hydride bond strength on the electron-donating or -withdrawing effect of the substituent in the H-Mg-Y subunits. Furthermore, the electronic and structural properties of the resulting complexes have been unveiled by means of the atoms in molecules (AIM) and natural bond orbital (NBO) analyses. Finally, several correlative relationships between interaction energies and various properties such as binding distance, frequency shift, molecular electrostatic potential, and intermolecular density at bond critical point have been checked for all studied systems.     

EPR of Mn2+-doped single crystals of [Co(H2O)6]XY6 (X = Si,Sn, Pt; Y = F,Cl)

EPR studies have been carried out in Mn²⁺-doped paramagnetic single crystals of [Co(H₂O)₆]SiF₆, [Co(H₂O)₆]SnF₆ and [Co(H₂O)₆]PtCl₆ at different temperatures using X-band microwave frequency. Reasonably sharp resonance lines of Mn²⁺ have been obtained at room temperature and analysed using the spin Hamiltonian for trigonal symmetry. From the observed linewidths of Mn^{2 +} in single-crystal [Co(H₂O)₆]SiF₆ the T₁ of Co^{2 +} has been estimated to be = 3.3 × 10⁻¹¹ s. A st phase transition in [Co(H₂O)₆]SiF₆ was observed and the low-temperature phase studied.     

Application of the FAKE molecular-orbital method to diatomic molecules XY (X,Y = H,F, Cl,Br, I)

Theoretical Chemistry Accounts - The FAKE (fast, accurate kinetic energy) method of semiempirical molecular orbital calculation is applied to diatomic molecules XY (X, Y= H, F, Cl, Br, I). The...     

C70X2(X=H,F, Cl)的稳定性和电子光谱

用INDO方法研究C70H2四种异构体的稳定性, 表明其最稳定异构体为1, 9-C70H2和7, 8-C70H2, 两者能量差为16.3KJ.mol^-^1, 与实验值及ab initio计算值接近; 光谱计算表明, 其特征吸收峰与实验值一致。在此基础上预测C70F2和C70Cl2的稳定性和电子光谱, 表明C70F2四种异构体的稳定性顺序与C70H2一致, 而C70Cl2则以21, 42-异构体最为稳定。二者的电子光谱与C70H2极其相似只是在500nm以上有细微差别。     

钠卤类硅烯的结构与稳定性

本文用RHF/STO-3G解析梯度方法研究了钠卤类硅烯H2SiNaX(X=F,Cl)势能面的主要特征, 得到了它们的几种平衡构型及其异构化的过渡态构型。与锂卤类硅烯相似, 三元环构型和P-配合物构型是最稳定的构型; 其它的两种构型-σ-配合物和经典"四面体", 也是势能面上的极值点, 但能量相对较高, 不稳定。文中对各构型的特点进行了分析。     

OXO(X=Cl,Br)与H(2S)反应机理的量子化学研究

卤素二氧化物自由基(OXO)是重要的同温层物种，在极地的夜间含量尤其丰富，由于其光解产生的卤素原子会造成臭氧的严重损耗，因此愈来愈受到人们的广泛关注嘲．OXO自由基可以通过自由基之间的反应或微波放电等方式获得，有关它与其它大气物种的反应已经进行了许多研究，Bemand等在实验上研究了Cl，H，O与OClO的反应动力学，对于OBrO与H之间的反应目前尚未见实验研究的报道，只有Guhua等在理论上对部分XBrO2(X=H，Cl，Br)同分异构体进行了研究．     

A challenge for density functional theory: the XONO and XNO2 (X=F, Cl, and Br) molecules

The equilibrium geometries, harmonic frequencies, dipole moments, infrared intensities, and relative energies of the cis-XONO, trans-XONO, and XNO₂ (X=F, Cl, and Br) have been investigated using four functionals in common use in Kohn-Sham density functional theory (DFT) calculations. Two of the functionals include non-local or gradient correction terms, while the other two also incorporate some exact Hartree-Fock exchange and are labeled hybrid functionals. The quality of the results obtained from the functionals is determined by comparison to previously published high-level coupled-cluster calculations. The hybrid functionals perform better for prediction of the equilibrium geometries, where the two gradient corrected functionals yield qualitatively incorrect molecular structures for cis-FONO and cis-ClONO. None of the functionals perform well in predicting all six harmonic frequencies, showing that the correlation between equilibrium geometries and harmonic frequencies is not as strong for these DFT methods as it is for conventional wavefunction ab initio methods, such as coupled-cluster theory. Results from the various functionals generally come into better agreement with each other and also with the coupled-cluster results moving down the periodic table. Received: 12 February 1997 / Accepted: 25 March 1997     

Use of oxidation-state differences and molecular orbitals to interpret bonding in the series ONXYZ (X, Y, Z=H, F, Cl), HNNX3, HNNX2Y, and HNNXY2 (X, Y=H, F) and OCX3-, OCX2Y-, and OCXY2- (X, Y=H, F)

Bonding in the series ONXYZ (X, Y, Z=H, F, Cl), HNNX3, HNNX2Y, HNNXY2 (X, Y=H, F), and OCX3-, OCX2Y-, OCXY2- (X, Y=H, F) shows evidence of a significant ionic contribution modifying the underlying covalent bonding. Increased ionic character can be correlated with oxidation-state differences between the bound atoms and is expressed in terms of shorter bond lengths. All members of the series, with the exception of ONH3, HNNH3, and OCH3-, possess a multiple O-N, N-N, or C-O bond modified by the ionic character of that bond. The O-N, N-N, and O-C single bonds in ONH3, HNNH3, and OCH3-, respectively, show some variation in length relative to "typical" single bonds of these types due to differences in ionic character. The two highest-occupied molecular orbitals in the ONXYZ or OCXYZ- (X, Y, Z=H, F) series which are piNO or piCO (when X=Y=Z=H) exhibit a distinct shift in their nodal plane as hydrogen is replaced by fluorine. The nodal plane moves from a location between the oxygen and the nitrogen or carbon to between the nitrogen or carbon and the fluorines impacting on the nature and length of the bonds joining these atoms. The pattern of N-F and C-F bond lengths in the series, ONH3-ONF3 and OCH3--OCF3-, respectively, lends support to the idea of resonance structures of the form ONXY+ F- or OCXY F- (where X, Y=H, F).     

Trends of the bonding effect on the performance of DFT methods in electric properties calculations: A pattern recognition and metric space approach on some XY2 (X = O,S and Y = H,O, F,S, Cl) molecules

A test set of 10 molecules (open and ring forms of ozone and sulfur dioxide as well as water and hydrogen sulfide and their respective fluoro‐ and chloro‐substituted analogs) of specific atmospheric interest has been formed as to assess the performance of various density functional theory methods in (hyper)polarizability calculations against well‐established ab initio methods. The choice of these molecules was further based on (i) the profound change in the physics between isomeric systems, e.g., open (C_2v) and ring (D_3h) forms of ozone, (ii) the relation between isomeric forms, e.g., open and ring form of sulfur dioxide (both of C_2v symmetry), and (iii) the effect of the substitution, e.g., in fluoro‐ and chloro‐substituted water analogs. The analysis is aided by arguments chosen from the information theory, graph theory, and pattern recognition fields of Mathematics: In brief, a multidimensional space is formed by the methods which are playing the role of vectors with the independent components of the electric properties to act as the coordinates of these vectors, hence the relation between different vectors (e.g., methods) can be quantified by a proximity measure. Results are in agreement with previous studies revealing the acceptable and consistent behavior of the mPW1PW91, B3P86, and PBE0 methods. It is worth noting the remarkable good performance of the double hybrid functionals (namely: B2PLYP and mPW2PLYP) which are for the first time used in calculations of electric response properties. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Theoretical Studies of Inorganic Compounds. 19 1) Quantum Chemical Investigations of the Phosphane Complexes X3B‐PY3 and X3Al‐PY3 (X = H,F, Cl; Y = F,Cl, Me,CN)

We report about quantum chemical ab initio calculations at the MP2/6‐311+G(2d)//MP2/6‐31G(d) level and DFT calculations at BP86/TZP of the geometries and bond dissociation energies of the borane‐phosphane complexes X₃B‐PY₃ and the alane‐phosphane complexes X₃Al‐PY₃ (X = H, F, Cl; Y = F, Cl, Me, CN). The nature of the B‐P and Al‐P bonds is analyzed with a bond energy partitioning method. The calculated bond dissociation energies D_e of the borane adducts X₃B‐PY₃ show for the phosphane ligands the trend PMe₃ > PCl₃ ∼ PF₃ > P(CN)₃. A similar trend PMe₃ > PCl₃ > PF₃ > P(CN)₃ is predicted for the alane complexes X₃Al‐PY₃. The order of the Lewis acid strength of the boranes depends on the phosphane Lewis base. The boranes show with PMe₃ and PCl₃ the trend BH₃ > BCl₃ > BF₃ but with PF₃ and P(CN)₃ the order is BH₃ > BF₃ > BCl₃. The bond energies of the alane complexes show always the trend AlCl₃ ≥ AlF₃ > AlH₃. The bonding analysis shows that it is generally not possible to correlate the trend of the bond energies with one single factor which determines the bond strength. The preparation energy which is necessary to deform the Lewis acid and Lewis base from the equilibrium form to the geometry in the complex may have a strong influence on the bond energies. The intrinsic interaction energies may have a different order than the bond dissociation energies. The trend of the interaction energies are sometimes determined by a single factor (Pauli repulsion, electrostatic attraction or covalent bonding) but sometimes all components are important. The higher Lewis acid strength of BCl₃ compared with BF₃ in strongly bonded complexes is not caused by the deformation energy of the fragments but it is rather caused by the intrinsic interaction energy. P(CN)₃ is a weaker Lewis base than PF₃, PCl₃ and PMe₃ mainly because of its weaker electrostatic attraction. The complex H₃B‐P(CN)₃ is predicted to have a bond dissociation energy D_o = 14.8 kcal/mol which should be sufficient to synthesize the compound as the first adduct with the ligand P(CN)₃. The calculated bond energies at the BP86 level are in most cases very similar to the MP2 results. In a few cases significantly different absolute values have been found which are caused by the method and not by the quality of the basis set.     

Structures and electron attachment properties of halomethanes (CX(n)Y(m), X=H, F; Y=Cl, Br, I; n=0,4; m=4-n)

Theoretical studies of structures of neutral molecules and their anions as well as dissociative electron attachment properties are presented for the halomethanes of general formula CX(n) Y(m); X=H, F; Y=Cl, Br, I; n=0,4; m=4-n. The dissociative electron attachment seems to be the primary process resulting in toxicity of these species. The halomethane anions containing hydrogens are formed as radical-anion adducts. When H is replaced by F, these species become true sigma* radicals. The electron affinities are computed using a variety of computational techniques including the four-order M?ller-Plesset (MP4) technique that included 250 basis functions. It is challenging to compare the computed results with experiment due to dearth of experimental data and uncertainties in the existing experimental data. Thus in certain cases larger differences are found between the computed and experimental values.     

An MP2 study on pre-reactive complexes (CH_2)_2O … XY (X,Y = H,F,Cl,Br,I)

The structures, energies, atomic chaiges and IR spectra of complexes (CH2)2O…XY (X, Y = H, F, Cl, Br, and I) have been examined by means of ab initio molecular orbital theory at the second-order level of Moller-Plesset perturbation correction. It is found that the hydrogen bond O…H-X is non-linear. The attraction between X and the H atoms in oxirane ring causes O…H-X bond bending. The hydrogen bond slighdy weakens the bond strength of C-O, and leads the bending and stretching mode of IR to shift to the red. The calculation results show that there is no evidence of a significant extent of proton transfer to give (CH2)2OH …X- in the isolated complexes.     

Diboron Bonds Between BX3 (X=H,F, CH3) and BYZ2 (Y=H,F; Z=CO,N2, CNH)

The ability of B atoms on two different molecules to engage with one another in a noncovalent diboron bond is studied by ab initio calculations. Due to electron donation from its substituents, the trivalent B atom of BYZ₂ (Z=CO, N₂, and CNH; Y=H and F) has the ability to in turn donate charge to the B of a BX₃ molecule (X=H, F, and CH₃), thus forming a B⋅⋅⋅B diboron bond. These bonds are of two different strengths and character. BH(CO)₂ and BH(CNH)₂, and their fluorosubstituted analogues BF(CO)₂ and BF(CNH)₂, engage in a typical noncovalent bond with B(CH₃)₃ and BF₃, with interaction energies in the 3–8 kcal/mol range. Certain other combinations result in a much stronger diboron bond, in the 26–44 kcal/mol range, and with a high degree of covalent character. Bonds of this type occur when BH₃ is added to BH(CO)₂, BH(CNH)₂, BH(N₂)₂, and BF(CO)₂, or in the complexes of BH(N₂)₂ with B(CH₃)₃ and BF₃. The weaker noncovalent bonds are held together by roughly equal electrostatic and dispersion components, complemented by smaller polarization energy, while polarization is primarily responsible for the stronger ones.     

Structures and Electron Attachment Properties of Halomethanes (CX n Y m,X = H,F; Y = Cl,Br, I; n = 0,4; m = 4n)

Abstract

Theoretical studies of structures of neutral molecules and their anions as well as dissociative electron attachment properties are presented for the halomethanes of general formula CX _n Y _m ; X = H, F; Y = Cl, Br, I; n = 0,4; m = 4 – n. The dissociative electron attachment seems to be the primary process resulting in toxicity of these species. The halomethane anions containing hydrogens are formed as radical-anion adducts. When H is replaced by F, these species become true σ? radicals. The electron affinities are computed using a variety of computational techniques including the four-order M?ller-Plesset (MP4) technique that included 250 basis functions. It is challenging to compare the computed results with experiment due to dearth of experimental data and uncertainties in the existing experimental data. Thus in certain cases larger differences are found between the computed and experimental values.