1,3-D,A取代方酸极化率和一阶超极化率的从头算研究

本文在RHF／6-311 G(3d，2p)水平上，对一系列1，3取代方酸衍生物的几何构型进行了全优化，得到分子最稳定构型，并用CPHF方法计算了分子线性极化率和一阶超极化率，分析了取代基对分子线性和非线性光学性质的影响。     

Experimental and computational studies of the phenyl radical reaction with propyne.

The kinetics for the gas-phase reaction of phenyl radical with propyne has been measured by cavity ring-down spectrometry (CRDS), and the mechanism and initial product branching have been elucidated with the help of quantum chemical calculations. Absolute rate constants measured by the CRDS technique can be expressed by the following Arrhenius equation: (k/cm(3) mol(-1) s(-1)): k(propyne)(T=301-428 K)=(3.68+/-0.92) x 10(11)exp[-(1685+/-80)/T]. The experiment is unable to distinguish between the possible reactive channels, but theory indicates that phenyl radicals preferably add to the unsaturated terminal carbon atom in propyne under our experimental conditions. Theoretical kinetic calculations, employing high-level G2M(RCC, RMP2) and G3 energetic and IRCMax(RCCSD(T)//B3LYP-DFT) molecular parameters, reproduce the total experimental rate constants within a factor of three. Calculated total and branching rate constants are provided for high-T kinetic modeling. Addition reactions of phenyl to C3H4 are estimated to be less important molecular-growth pathways in high-T conditions (T>1000 K) in comparison to the C6H5 + C2H2 reaction.     

Theoretical predictions of the structure,gas‐phase acidity,and aromaticity of tetrathiosquaric acid

Results of ab initio self‐consistent‐field and density functional theory calculations of the gas‐phase structure, acidity (free energy of deprotonation, ΔG⁰), and aromaticity of tetrathiosquaric acid (3,4‐dithiohydroxy‐3‐cyclobutene‐1,2‐dithione, H₂C₄S₄) are reported. The global minimum found on the potential energy surface of tetrathiosquaric acid presents a planar conformation. The ZZ isomer was found to have the lowest energy among the three planar conformers and the ZZ and ZE isomers are very close in energy. The optimized geometric parameters exhibit a bond length equalization relative to reference compounds, cyclobutanedithione, and cyclobutenedithiol. The computed aromatic stabilization energy by homodesmotic reaction is −18.4 (MP2(fu)/6‐311+G**//RHF/6‐311+G**) and −15.1 kcal/mol (B3LYP//6‐311+G**// B3LYP/6‐311+G**). The aromaticity of tetrathiosquaric acid is indicated by the calculated diamagnetic susceptibility exaltation (Λ) −11.77 (CSGT(IGAIM)‐RHF/6‐311+G**// RHF/6‐311+G**) and −18.08 (CSGT(IGAIM)‐B3LYP/6‐311+G**// B3LYP/6‐311+G**). Thus, tetrathiosquaric acid fulfils the geometric, energetic and magnetic criteria of aromaticity. The most reliable theoretical gas‐phase acidities are $\Delta G^{0}_{1(298\mathrm{K})}=303.7$ and $\Delta G^{0}_{2(298\mathrm{K})}=394.1$ kcal/mol. Hence, tetrathiosquaric acid is a stronger acid than squaric acid (3,4‐dihydroxy‐3‐ cyclobutene‐1,2‐dione, H₂C₄O₄). Comparisons of the computed results of tetrathiosquaric acid with squaric acid have also been made. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 443–449, 2000     

Anion-pi interactions in cyanuric acids: a combined crystallographic and computational study

Several structures of pi complexes of isocyanuric acid and of several thio derivatives with anions have been computed by using high level ab initio calculations. The nature of the complexes has been studied by means of the method of molecular interaction potential with polarization (MIPp) and Bader's theory of atoms-in-molecules. These molecules form favorable complexes with anions and can be used as binding units for building receptors for the molecular recognition of anions. In several cases, the anion-pi interaction has been demonstrated experimentally by means of X-ray crystallography.     

甲醛与甲酸相互作用的从头算研究

在MP2／6—31G(d)和MP2(FC)／6—311 G(d，p)水平上，对H2CO和HCOOH以及设计的4种构型H2CO…HCOOH复合物等进行几何全优化计算，经振动频率分析，确认它们为势能面上的稳定驻点．然后在MP2／6—311 G(2df，2p)水平上进行单点能计算和基组重叠误差(BSSE)校正以获得相互作用能，并利用自然键轨道理论探讨H2CO和H(X)OH相互作用的本质。     

某些含氧酸中原子净电荷与pK_a的关系

在量子化学从头算 (分别采用HF/STO 3G ,HF/STO 3G ,HF/ 3 2 1G和MP2 /STO 3G)水平下 ,对 11种含氧酸分子 (H2 O ,H2 O2 ,HNO2 ,HNO3 ,HClO ,HClO2 ,HClO3 ,HCO2 H ,CH3 CO2 H ,HBrO和HIO)进行几何构型全优化和电子结构计算 ,将获得的电子结构数据与 pKa 实验值进行逐步回归分析 ,结果发现含氧酸中羟基氢原子净电荷、含氧酸中与羟基氧相连非氢原子净电荷两项与 pKa 呈良好的二元线性关系 ,表明含氧酸中羟基氧所连原子的净电荷在决定含氧酸的酸性上起重要作用     

Ab initio study of the positronation of the CaO and SrO molecules including calculation of annihilation rates

Ab initio multireference single‐ and double‐excitation configuration interaction calculations have been performed to compute potential curves for ground and excited states of the CaO and SrO molecules and their positronic complexes, e⁺CaO, and e⁺SrO. The adiabatic dissociation limit for the ²Σ⁺ lowest states of the latter systems consists of the positive metal ion ground state (M⁺) and the OPs complex (e⁺O^?), although the lowest energy limit is thought to be e⁺M + O. Good agreement is found between the calculated and experimental spectroscopic constants for the neutral diatomics wherever available. The positron affinity of the closed‐shell X ¹Σ⁺ ground states of both systems is found to lie in the 0.16–0.19 eV range, less than half the corresponding values for the lighter members of the alkaline earth monoxide series, BeO and MgO. Annihilation rates (ARs) have been calculated for all four positronated systems for the first time. The variation with bond distance is generally similar to what has been found earlier for the alkali monoxide series of positronic complexes, falling off gradually from the OPs AR value at their respective dissociation limits. The e⁺SrO system shows some exceptional behavior, however, with its AR value reaching a minimum at a relatively large bond distance and then rising to more than twice the OPs value close to its equilibrium distance. © 2012 Wiley Periodicals, Inc.     

Ab initio and density functional theory studies of the structure,gas-phase acidity and aromaticity of tetraselenosquaric acid

Results of ab initio self-consistent-field (SCF) and density functional theory (DFT) calculations of the gas-phase structure,acidity (free energy of deprotonation,G0) and aro-maticity of tetraselenosquaric acid (3,4-diselenyl-3-cy-dobutene-1,2-diselenone,H2C4Se4) are reported.The global minimum found on the potential energy surface of tetraselenosquaric acid presents a planar conformation.The ZZ iso-mer was found to have the lowest energy among the three planar conformers and the ZZ and ZE isomers are very dose in energy.The optimized geometric parameters exhibit a bond length equalization relative to reference compounds,cyclobu-tanediselenone,and cydobutenediselenol.The computed aromatic stabilization energy (ASE) by homodesmotic reaction is -77.4 (MP2(fu)/6 - 311 G //RHF/6 - 311 G) and - 54.8 kJ/mol (B3LYP/6 - 311 G //B3LYP/6 -311 G).The aromaticity of tetraselenosquaric add is indicated by the calculated diamagnetic susceptibility exaltation (A) - 19.13 (CSGT(IGAEM) - RHF/6 - 311 G// RHF/6-     

Toward the preparation of the HAuF6, HAu2F11, and HAu3F16 superacids: Theoretical study

We apply the B3LYP and QCISD theoretical methods with the 6-311++G(d,p) basis set and the LANL2DZ effective core potentials to investigate the reaction paths leading to the preparation of the HAuF₆, HAu₂F₁₁, and HAu₃F₁₆ superacids. The Gibbs free energies of deprotonation of these systems are calculated to estimate their acid strength. The thermodynamic stability of the corresponding anionic precursors ((AuF₆)^–, (Au₂F₁₁)^–, and (Au₃F₁₆)^–) and their vertical excess electron detachment energies are evaluated and discussed. The suggested route of the HAu₂F₁₁ preparation involves the F^– attachment to the Au₂F₁₀ reactant which results in the formation of the (Au₂F₁₁)^– anion whose protonation yields the HAu₂F₁₁ superacid. The suggested HAuF₆ superacid preparation route is based on a qualitatively similar scheme and involves the conversion of the AuF₅ reactant into the (AuF₆)^– anion followed by its protonation. The proposed path for the HAu₃F₁₆ superacid preparation involves the attachment of AuF₅ to HAu₂F₁₁.     

Direct ab initio dynamics study on the rate constants and kinetic isotope effect for the reactions of H atoms with GeDn(CH3)4-n (n = 1-4)

Direct ab initio dynamic calculations are performed on the reactions of atomic hydrogen with GeD(n)(CH(3))(4-n) (n = 1-4) over the temperature range 200-2000 K at the PMP4SDTQ/6-311 +G(3df,2p)//MP2/6-31 +G(d) (for n = 2-4) and G2//MP2/6-31 +G(d) (for n = 1) levels. The corresponding k(H)/k(D) ratios are then calculated in order to determine the kinetic isotope effect for the four reactions. For the simplest GeD(4) +H reaction, the only one that has available experimental data, the calculated canonical variational transition state theory incorporates small-curvature tunneling correction (CVT/SCT) thermal rate constants, and the k(H)/k(D) values are in good agreement with the experimental values within the experimental temperature range 293-550 K. For the four GeD(n)(CH(3))(4-4) (n = 1-4) reactions, the variational effect is small over the whole temperature range, whereas the small-curvature effect is important in the lower temperature range. Finally, the overall rate constants are fitted to the three-parameter expression over the whole temperature range 200-2000 K as 5.8 x 10(8)T(1.68)exp(-929/T), 1.7 x 10(8)T(1.80)exp(-691/T), 2.58 x 10(8)T(1.71)exp(-706/T), and 1.0 x 10(7)T(2.08)exp(-544/T) cm(3) mol(-1) s(-1) for the n = 4, 3, 2, and 1 reactions. Our work may represent the first theoretical study of the kinetic isotope effect for the H-attack on the G-H bonding.     

1, 2-二硒方酸气相酸性和芳香性的理论研究

在RHF/6-311G^*^*水平优化得到1,2-二硒方酸(3,4-二羟基-3-环丁烯-1,2-二硒酮)三种平面构象异构体的平衡几何构型。进一步用MP2(full)/6-311G^*//RHF/6-311G^*^*方法计算三种异构体的单点能量,发现ZZ型异构体是能量最低构象,且ZZ和ZE型能量非常接近。用优化的最稳定构象ZZ型异构体在RHF/6-311G^*^*//RHF/6-311G^*^*,RHF/6-311+G^*^*//RHF/6-311+G^*^*,MP2(full)/6-311+G^*^*//RHF/6-311+G^*^*和B3LYP/6-311+G^*^*//B3LYP/6-311+G^*^*水平计算其气相酸性[ΔGⅲ~(~2~9~8~K~)]和同键反应芳香性稳定化能(HASE)。用基团加和法(groupincrementapproach)在RHF/6-311+G^*^*//RHF/6-311+G^*^*和B3LYP/6-311+G^*^*//B3LYP/6-311+G^*^*水平计算其磁化率增量(Λ)。计算结果指出标题化合物的键长发生了平均化,同键反应芳香性稳定化能和磁化率增量均为负值,表明它具有芳香性,实现了标题化合物芳香性的几何、能量和磁性的判定。     

Gas-phase protonation and deprotonation of acrylonitrile derivatives N[triple chemical bond]C--CH==CH--X (X=CH(3), NH(2), PH(2), SiH(3))

A combined experimental and theoretical study on the gas-phase basicity and acidity of a series of cyanovinyl derivatives is presented. The gas-phase basicities and acidities of (N[triple chemical bond]C--CH==CH--X, X=CH(3), NH(2)) were obtained by means of Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry techniques. The corresponding calculated values were obtained at the G3B3 level of theory. The effects of exchanging CH(3) for SiH(3), and NH(2) for PH(2), were analyzed at the same level of theory. For the neutral molecules, the Z isomer is always the dominant species under standard gas-phase conditions at 298 K. The loss of the proton from the substituent X was found systematically to be much more favorable than deprotonation of the HC==CH linking group. The corresponding isomeric E ion is much more stable than the Z ion, so that only the former should be found in the gas phase. The most significant structural changes upon deprotonation occur for the methyl and amino derivatives because, in both cases, deprotonation of X leads to a significant charge delocalization in the corresponding anion. Protonation takes place systematically at the cyano group, whereby the isomeric E ion is again more stable than the Z ion. Push-pull effects explain the preference of aminoacrylonitrile to be protonated at the cyano group, which also explains the high basicity of this derivative relative to other members of the analyzed series that present rather similar gas-phase basicities, GB approximately 780 kJ mol(-1), indicating that the different nature of the substituents has only a weak effect on the intrinsic basicity of the cyano group. The cyanovinyl derivatives have a significantly stronger gas-phase acidity than that of the corresponding vinyl compounds CH(2)==CH--X. This acidity-strengthening effect of the cyano group is attributed to the greater stabilization of the anion with respect to the corresponding neutral compound.     

Encasing of Na+ ion in dimer‐formed acetic acid clusters

Peaks with anomalous abundance found in the mass spectra are associated with ions with enhanced stability. Among the scientific community focused on mass spectrometry, these peaks are called ‘magic peaks’ and their stability is often because of suggestive symmetric structures. Here, we report findings on ionised Na‐acetic acid clusters [Na⁺‐(AcA)_n] produced by Na‐doping of (AcA)_n and UV laser ionisation. Peaks labelled n = 2, 4, 8 are clearly distinguishable in the mass spectra from their anomalous intensity. Ab initio calculations helped elucidate cluster structures and energetic. A plausible interpretation of the magic peaks is given in terms of (AcA)_n formed by dimer aggregation. The encasing of Na⁺ by twisted dimers is proposed to be the origin of the enhanced cluster stability. A conceivable dimer‐formed tube‐like closed structure is found for the Na⁺‐(AcA)₈. Copyright © 2015 John Wiley & Sons, Ltd.     

Kinetics of the hydrogen abstraction reactions of 1,1‐ and 1,2‐difluoroethane with hydroxyl radical: an ab initio study

The hydrogen abstraction reactions of 1,1‐ and 1,2‐difluoroethane with the OH radical have been investigated by the ab initio molecular orbital theory. The geometries of the reactants, products, and transition states have been optimized at the (U)MP2=full level of theory in conjunction with 6‐311G(d,p) basis functions. Single‐point (U)MP2=full with larger basis set, such as 6‐311G(3d,2p), and QCISD(T)=full/6‐311G(d,p) calculations have also been carried out to observe the effects of basis sets utilized and higher order electron correlation. Three and four reaction channels have been identified for 1,1‐ and 1,2‐difluoroethane, respectively. In the case of 1,1‐difluoroethane, hydrogen abstraction from the α‐carbon has been found to be easier than that from the β‐carbon. The barriers of the four reaction channels for 1,2‐difluoroethane are close to each other. Weak hydrogen bonding interactions have been observed between hydroxyl hydrogen and a fluorine atom in the transition states. Rate constants for the reactions of 1,1‐ and 1,2‐difluoroethane with the OH radical have been calculated using the standard transition state theory and found to be in good agreement with the experimental results. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 1305–1318, 2000     

Ab initio direct dynamics studies on the reaction of H atom with CH3CH2Cl

The multiple channel reaction H + CH(3)CH(2)Cl --> products has been studied by the ab initio direct dynamics method. The potential energy surface information is calculated at the MP2/6-311G(d,p) level of theory. The energies along the minimum energy path are further improved by single-point energy calculations at the PMP4(SDTQ)/6-311+G(3df,2p) level of theory. For the reaction, four reaction channels (one chlorine abstraction, one alpha-hydrogen abstraction, and two beta-hydrogen abstractions) have been identified. The rate constants for each reaction channel are calculated by using canonical variational transition state theory incorporating the small-curvature tunneling correction in the temperature range 298-5000 K. The total rate constants, which are calculated from the sum of the individual rate constants, are in good agreement with the experimental data. The calculated temperature dependence of the branching fractions indicates that for the title reaction, H-abstraction reaction is the major reaction channel in the whole temperature range 298-5000 K.     

Strong dissimilarities between the gas-phase acidities of saturated and alpha,beta-unsaturated boranes and the corresponding alanes and gallanes

The effect that unsaturation has on the intrinsic acidity of boranes, alanes, and gallanes, was analyzed by B3 LYP and CCSD(T)/6-311+G(3df,2p) calculations on methyl-, ethyl-, vinyl-, and ethynylboranes, -alanes and -gallanes, and on the corresponding hydrides XH3. Quite unexpectedly, methylborane, which behaves as a carbon acid, is predicted to have an intrinsic acidity almost 200 kJ mol(-1) stronger than BH3, reflecting the large reinforcement of the C--B bond, which upon deprotonation becomes a double bond through the donation of the lone pair created on the carbon atom into the empty p orbital of the boron. Also unexpectedly, and for the same reason, the saturated and alpha,beta-unsaturated boranes are much stronger acids than the corresponding hydrocarbons, in spite of being carbon acids as well. The Al derivatives also behave as carbon acids, but in this case the most favorable deprotonation process occurs at C beta, leading to the formation of rather stable three-membered rings, again through the donation of the C beta lone pair into the empty p orbital of Al. For Ga-containing compounds the deprotonation of the GaH2 group is the most favorable process. Therefore only Ga derivatives behave similarly to the analogues of Groups 14, 15, and 16 of the periodic table, and the saturated derivatives exhibit a weaker acidity than the unsaturated ones. Within Group 13, boranes are stronger acids than alanes and gallanes. For ethyl and vinyl derivatives, alanes are stronger acids than gallanes. We have shown, for the first time, that acidity enhancement for primary heterocompounds is not only dictated by the position of the heteroatom in the periodic table and the nature of the substituent, but also by the bonding rearrangements triggered by the deprotonation of the neutral acid.