An evaluation of various computational methods for the treatment of organoselenium compounds

A reliable computational method for the prediction of organoselenium geometries and bond dissociation energies (BDEs) has been determined on the basis of the performance of density functional theory (DFT: B3LYP and B3PW91) and ab initio molecular orbital procedures (Hartree-Fock (HF)) in conjunction with various Pople basis sets including (but not limited to) the 6-31G(d), 6-31G(d,p), 6-311G(d), 6-311G(d,p), 6-311G(2df,p), and 6-311G(3df,3pd) sets. Predicted geometries and BDEs are compared with available experimental data and quadratic configuration interaction including single and double substitutions (QCISD) results. The B3PW91/6-311G(2df,p) level of theory is recommended for the prediction of the geometries and energetics of organoselenium compounds.     

A new 3D-aromatic organoboron species on the basis of CB6 unit: Two states of carbon hypercoordination and structural isomerism of non-classical forms

Calculations using the DFT B3LYP/6-311+G(df,2p) and wb97XD/6-311+G(df,2p) methods predict the stability of a new three-dimensional aromatic organoboron species based on the CB₆ unit. Derived systems containing carbon atoms in two hypercoordination states (planar hexacoordination and inverted-umbrella heptacoordination) can be used as building blocks for constructing various stable condensed systems with many hypercoordinated carbon centers.     

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.     

DFT-GIAO calculations of 19F NMR chemical shifts for perfluoro compounds

The (19)F NMR shieldings for 53 kinds of perfluoro compounds were calculated by the B3LYP-GIAO method using the 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-31++G(d,p), 6-311G(d,p), 6-311++G(d,p), 6-311G(2d,2p), 6-311++G(2d,2p), 6-311++G(2df,2p), 6-311++G(3d,2p), and 6-311++G(3df,2p) basis sets. The diffuse functions markedly reduce the difference between the calculated and experimental chemical shifts. The calculations using the 6-31++G(d,p) basis set give the chemical shifts within 10 ppm deviations from experimental values except for the fluorine nuclei attached to an oxygen atom, a four- and a six-coordinated sulfur atom, and FC(CF(3))(2) attached to a sulfur atom.     

Theoretical studies on dynamics and thermochemistry of the reactions CF(3)CHCl(2)+Cl-->CF(3)CCl(2)+HCl and CF(3)CHFCl+Cl-->CF(3)CFCl+HCl

A dual-level direct dynamics study has been carried out for the two hydrogen abstraction reactions CF(3)CHCl(2)+Cl and CF(3)CHFCl+Cl. The geometries and frequencies of the stationary points are optimized at the BHLYP/6-311G(d,p), B3LYP/6-311G(d,p), and MP2/6-31G(d) levels, respectively, with single-point calculations for energy at the BHLYP/6-311++G(3df,2p), G3(MP2), and QCISD(T)/6-311G(d,p) levels. The enthalpies of formation for the species CF(3)CHCl(2), CF(3)CHFCl, CF(3)CCl(2), and CF(3)CFCl are evaluated at higher levels. With the information of the potential energy surface at BHLYP/6-311++G(3df,2p)//6-311G(d,p) level, we employ canonical variational transition-state theory with small-curvature tunneling correction to calculate the rate constants. The agreement between theoretical and experimental rate constants is good in the measured temperature range 276-382 K. The effect of fluorine substitution on reactivity of the C-H bond is discussed.     

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.     

Basis set dependence of the vibrational wavenumbers of the out-of-plane modes of conjugated π-electron ring systems

Basis set dependence of the vibrational wavenumbers of out-of-plane modes calculated at the MP2 level of ab initio molecular orbital theory is examined for benzene, p-benzoquinodimethane, p-benzoquinone, furan, and thiophene. Various polarization functions up to (3df,p) are used in combination with the 6-31G and 6-311G basis sets. It is shown that, especially in the case of normal modes with alternate out-of-plane motions of the carbon atoms (such as the ν₄ (b_2g) mode of benzene), the calculated wavenumbers depend strongly on the exponents (_d) of the d functions on the carbon atoms. It is therefore necessary to include d functions with an optimum exponent (_d0.4) on the carbon atoms to obtain reasonable out-of-plane vibrational force fields. In a few cases (such as the ν₁₆ (a₂) mode of furan), inclusion of a set of f functions on the carbon atoms has some effects on the calculated wavenumbers of out-of-plane modes. However, unless the basis set contains an optimum set of d functions, inclusion of a set of f functions does not improve the agreement between the observed and calculated vibrational wavenumbers. As a case with an exaggerated effect of basis set, it is shown that the wavenumber of the ν₄₁ (b_2g) mode of the planar optimized structure of p-benzoquinodimethane is calculated to be imaginary by using the 6-311G(d,p) or 6-311G(df,p) basis set at the MP2 level. For all the molecules treated in the present study, reasonable out-of-plane vibrational force fields are obtained by using the 6-31G(2df,p) and 6-311G(2df,p) basis sets.     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

Theoretical Studies on Structures and Stabilities of C4H2+ Isomers

The structures, energies, stabilities and spectroscopies of doublet C₄H₂⁺ cations were explored at the DFT/B3LYP/6-311G(d,p), CCSD(T)/6-311+G(2df,2pd)(single-point), and G3B3 levels. Ten minimum isomers including the chainlike, three-member-ring, and four-member-ring structures are interconverted by means of 15 interconversion transition states. The potential energy surface was investigated. At the CCSD(T)/6-311+G(2df,2pd) and G3B3 levels, the global minimum isomer was found to be a linear HCCCCH. The structures of the stable isomer and its relevant transition state are further optimized at the QCISD/6-311G(d,p) level. The bonding nature and structure of isomer HCCCCH were analyzed.     

Conformational studies in the cyclohexane series. 2. Phenylcyclohexane and 1-methyl-1-phenylcyclohexane

The structures and relative energies of the conformers of phenylcyclohexane, and 1-methyl-1-phenylcyclohexane have been calculated at theoretical levels including HF/6-31G, B3LYP/6-311G, MP2/6-311G, MP2/6-311(2df,p), QCISD/6-311G, and QCISD/6-311G(2df,p). The latter gives conformational enthalpy (DeltaH degrees ), entropy (DeltaS degrees ), and free energy (DeltaG degrees ) values for phenylcyclohexane that are in excellent agreement with the experimental data. The calculations for 1-methyl-1-phenylcyclohexane find a free energy difference of 1.0 kcal/mol at -100 degrees C, favoring the conformation having an axial phenyl group, that is in only modest agreement with the experimental value of 0.32 +/- 0.04 kcal/mol. The origin of the phenyl rotational profiles for the conformers of phenylcyclohexane and 1-methyl-1-phenylcyclohexane is discussed.     

Scaling factors for vibrational frequencies and zero-point vibrational energies of some recently developed exchange-correlation functionals

The scaling factors for the vibrational frequencies and zero-point vibrational energies evaluated at various combinations of recently developed exchange-correlation functionals and various basis sets are reported. The exchange-correlation functionals considered are B972, B98, HCTH, OLYP, O3LYP, G96LYP, PBE0 and VSXC functionals; the basis sets employed are 3-21G, 6-31G*, 6-31G**, 6-31+G, 6-311G*, 6-311G**, 6-311G(df,p), 6-311+G(df,p), cc-pVDZ and aug-cc-pVDZ. The experimental harmonic frequencies of 122 small molecules and the zero-point vibrational energies of 39 small molecules are used to determine the scaling factors through the least-square fitting procedure. It was found that the scaling factors do not depend significantly on the basis sets considered. The vibrational frequency scaling factors evaluated by using the B98 and PBE0 functionals are recommended on the basis of smallest root mean square error. The zero-point vibrational energy scaling factors evaluated from the B972 functional with Pople's double-zeta basis set and the HCTH functional with Pople's triple-zeta basis set are recommended on the basis of smallest root mean square error.     

Inductive effects in radicals calculated from DFT energies; substituted bicyclo[2.2.2]octan-1-yloxy radicals

Energies of a series of 4-substituted 1-oxybicyclo[2.2.2]octan-1-yloxy radicals with 18 various substituents were calculated within the framework of the DFT theory at the levels UB3LYP/6-311+G(d,p)//UB3LYP/6-311+G(d,p) and UB3LYP/6-311++G(2df,p)//UB3LYP/6-311+G(d,p) and compared with similar series of the parent alcohols, their deprotonated and protonated forms calculated at the levels B3LYP/6-311+G(d,p)//B3LYP/6-311+G(d,p) and B3LYP/6-311++G(2df,p)//B3LYP/6-311+G(d,p). The two levels are of the same performance and both are sufficient for molecules of this type according to comparison with scarce experimental gas-phase acidities and basicities. The substituent effects were analyzed in terms of isodesmic equations. In addition to strong dependence on the substituent inductive effect, a slight dependence on the electronegativity of the first atom of the substituent was proven in certain cases. In all aspects, there is no qualitative difference between the effects on radicals and on similar closed shell species. Radicals behave as slightly electron deficient; the substituent effect is weaker than that on the ions but stronger than on neutral molecules.     

Molecular structure and vibrational spectra of 3-amino-5-hydroxypyrazole by density functional method

In this work, we report a combined experimental and theoretical study on molecular and vibrational structure of 3-amino-5-hydroxypyrazole (3A5HP). The Fourier transform infrared and Fourier transform Raman spectra of 3A5HP were recorded in the solid phase. The molecular geometry and vibrational frequencies of 3A5HP in the ground state have been calculated by using the density functional method B3LYP with basis sets, 6-311++G(d,p), 6-311+G(3df,2p), 6-311+G(3df,2pd), CC-pVDZ, aug-CC-pVDZ and CC-pVTZ. The optimized geometrical parameters obtained by B3LYP show best agreement with the experimental values. The theoretical spectrograms for FT-IR and FT-Raman spectra of the title molecule have been constructed.     

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.     

Li+粘贴CHnCl3-nO2和CFnCl3-nO2( n = 1~3 )自由基：Li+–CHnCl3-nO2和 Li+–CFnCl3-nO2的结构与性质的研究

The structures and stabilities of these still experimentally unknown CH_nCl_(3-n)O_(2-)Li~ and CF_nCl_(3-n)O_(2-)Li~ ionshave been theoretically investigated by ab initio molecular orbital theory and density functional theory(DFT)inconjunction with the 6-311G(d,p),6-311 G(d,p),6-311 G(2d,p)and 6-311 G(2df,2p)basis sets.The optimizedgeometries,chemical bonding and NBO analysis indicate that these complexes of CH_nCl_(3-n)O_(2-)Li~ and CF_nCl_(3-n)O_(2-)Li~ exist as ion-dipole molecules.The calculated affinity energies of these species exceed 41.9 kJ/mol,which arelarge enough to suggest the possibility that these title complexes could be detected as stable species in gas phase byLi~ ion attachment mass spectrometry.     

TDDFT investigation on nucleic acid bases: comparison with experiments and standard approach

A comprehensive theoretical study of electronic transitions of canonical nucleic acid bases, namely guanine, adenine, cytosine, uracil, and thymine, was performed. Ground state geometries were optimized at the MP2/6-311G(d,p) level. The nature of respective potential energy surfaces was determined using the harmonic vibrational frequency analysis. The MP2 optimized geometries were used to compute electronic vertical singlet transition energies at the time-dependent density functional theory (TDDFT) level using the B3LYP functional. The 6-311++G(d,p), 6-311(2+,2+)G(d,p), 6-311(3+,3+)G(df,pd), and 6-311(5+,5+)G(df,pd) basis sets were used for the transition energy calculations. Computed transition energies were found in good agreement with the corresponding experimental data. However, in higher transitions, the Rydberg contaminations were also obtained. The existence of pisigma* type Rydberg transition was found near the lowest singlet pipi* state of all bases, which may be responsible for the ultrafast deactivation process in nucleic acid bases.     

优化几何构型对高级别从头算能量的影响

对《ＣＲＣ物理与化学手册》（第７７版）中第三周期以前的无机双原子分子,当其理论优化构型的相对误差大于２％时,分别在实验构型和最大偏差的理论构型下,计算了ＱＣＩＳＤ（Ｔ）／６－３１１＋Ｇ（３ｄｆ,２ｐ）能量并作了比较。结果表明,大多数能量的差别在４．２ｋＪ．ｍｏｌ＾－１以内。由此说明,目前一般采用的构型优化理论方法,多数情况下不至于明显影响单点高级别从头算的计算精度。同时还发现,Ｇ２（ＱＣＩ）的高级     

Boron,carbon, and aluminum supertetrahedral graphane analogues

The electronic and spatial structures of carbon, boron, and aluminum supertetrahedrane models of graphane have been studied by means of density functional theory methods in the supermolecular approximation (B3LYP/6-311G(df,2p)) and with imposing periodic boundary conditions (PBEPBE/6-311G (d,p), HSEH1PBE/6-311G (d,p)). Calculations predict that pure boron and aluminum structures are narrow-gap semiconductors. For supertetrahedral carbon graphane, calculations predict properties intermediate between the semiconductor and insulator properties. All bonds in the carbon system are two-center two-electron (2с–2е), while for the boron system, intratetrahedrane bonds are three-center two-electron (3с–2е), and intertetrahedrane bonds are common two-center two-electron bonds (2с–2е).     

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.     

Polynuclear sandwich derivatives of [10]annulene: A quantum chemical study

Density-functional theory (DFT) calculations at B3LYP/6-311+G(df,p) and M06/6-311+G(df,p) levels of theory predict the stability of new polynuclear sandwich complexes [Cu₅(C₁₀H₁₀)₂]⁺ and Ni₅(C₁₀H₁₀)₂ based on flat [10]annulene cycles.