Ab initio calculation of bowl, cage, and ring isomers of C20 and C20-

High-level ab initio calculations have been carried out to reexamine relative stability of bowl, cage, and ring isomers of C(20) and C(20)(-). The total electronic energies of the three isomers show different energy orderings, strongly depending on the hybrid functionals selected. It is found that among three popular hybrid density-functional (DF) methods B3LYP, B3PW91, PBE1PBE, and a new hybrid-meta-DF method TPSSKCIS, only the PBE1PBE method (with cc-pVTZ basis set) gives qualitatively correct energy ordering as that predicted from ab initio CCSD(T)/cc-pVDZ [CCSD(T)-coupled-cluster method including singles, doubles, and noniterative perturbative triples; cc-pVDZ-correlation consistent polarized valence double zeta] as well as from MP4(SDQ)/cc-pVTZ [MP4-fourth-order Moller-Plesset; cc-pVTZ-correlation consistent polarized valence triple zeta] calculations. Both CCSD(T) and MP4 calculations indicate that the bowl is most likely the global minimum of neutral C(20) isomers, followed by the fullerene cage and ring. For the anionic counterparts, the PBE1PBE calculation also agrees with MP4/cc-pVTZ calculation, both predicting that the bowl is still the lowest-energy structure of C(20)(-) at T=0 K, followed by the ring and the cage. In contrast, both B3LYP/cc-pVTZ and B3PW91/cc-pVTZ calculations predict that the ring is the lowest-energy structure of C(20)(-). Apparently, this good reliability in predicting the energy ordering renders the hybrid PBE method a leading choice for predicting relative stability among large-sized carbon clusters and other carbon nanostructures (e.g., finite-size carbon nanotubes, nano-onions, or nanohorns). The relative stabilities derived from total energy with Gibbs free-energy corrections demonstrate a changing ordering in which ring becomes more favorable for both C(20) and C(20)(-) at high temperatures. Finally, photoelectron spectra (PES) for the anionic C(20)(-) isomers have been computed. With binding energies up to 7 eV, the simulated PES show ample spectral features to distinguish the three competitive C(20)(-) isomers.     

B14: an all-boron fullerene

Experiments revealed that small boron cluster anions and cations are (quasi-)planar. For neutral boron cluster, (quasi-)planar motifs are also suggested to be global minimum by many theoretical studies, and a structural transformation from quasi-planar to double-ring tubular structures occurs at B(20). However, a missing opportunity is found for neutral B(14), which is a flat cage and more stable than the previous quasi-planar one by high level ab initio calculations. The B(14) cage has a large HOMO-LUMO gap (2.69 eV), and NICS values reveal that it is even more aromatic than the known most aromatic quasi-planar B(12) and double-ring B(20), which indicates a close-shell electronic structure. Chemical bonding analysis given by AdNDP reveals that the B(14) cage is an all-boron fullerene with 18 delocalized σ-electrons following the 2(n+1)(2) rule of spherical aromaticity. The geometry and bonding features of the B(14) cage are unique denying conversional thinking.     

An Ab Initio Molecular Orbital Study of Valence Bond Isomers of Silabenzene and Phosphabenzene

Ab initio molecular orbital calculation at HF/6-31G*, HF/6-31G**, HF/6-311G**, HF/6-311++G**, RMP2-FC/6-31G*, and B3LYP/6-31G* levels of theory for geometry optimization and MP4(SDQ)/6-31G* for a single point total energy calculation are reported for silabenzene ( 7 ), phosphabenzene ( 8 ) and 16 valence bond isomers of silabenzene and phosphabenzene ( 9-24 ). The calculated energy difference (19.78 kcal mol m 1 ) between silabenzene and the most stable valence bond isomer of silabenzene (1-silabenzvalene, 9 ) is much smaller than the difference (73.60 kcal mol m 1 ) between benzene and benzvalene ( 2 ). The energy difference between phosphabenzene and the most stable valence bond isomer of phosphabenzene (1-phosphabenzvalene, 17 ) is calculated to be 43.29 kcal mol m 1 .     

Quantum-chemical Study on the Structures and Properties of Uracil-BX3 （X = F, Cl） Complexes

1 INTRODUCTION The intermolecular interaction of bases in DNA or RNA is of immense interest and significance to che- mists and biologists alike. The interactions of these bases with metal cations, solvent molecules and other small molecules or ions would affect the struc- ture and biological properties or recognition process,which has been investigated widely[1～8]. Boron contained compounds are electron deficient com- pounds and have been extensively used as catalysts in chemical react…     

Molecular structure of the octatetranyl anion, C8H-: a computational study

The equilibrium molecular structure of the octatetranyl anion, C8H(-), which has been recently detected in two astronomical environments, is investigated with the aid of both ab initio post-Hartree-Fock and density functional theory (DFT) calculations. The model chemistry adopted in this study was selected after a series of benchmark calculations performed on molecular acetylene for which accurate gas-phase structural data are available. Geometry optimizations performed at the CCSD/6-311+G(2d,p), QCISD/6-311+G(2d,p), and MP4(SDQ)/6-311+G(2d,p) levels of theory yield for C8H(-) an interesting polyyne-type structure that defies the chemical formula displaying a simple alternation of triple and single carbon-carbon bonds, [:C[triple bond]C-C[triple bond]C-C[triple bond]C-C[triple bond]CH](1-). In the optimized geometry of C8H(-), as one proceeds from the naked carbon atom on one side of the chain to the CH unit on the opposite side of the chain, the short (formally triple) carbon-carbon bonds decrease in length from 1.255 to 1.213 A whereas the long (formally single) carbon-carbon bonds increase (albeit only slightly) in length from 1.362 to 1.378 A (CCSD results). In striking contrast, both MP2 and DFT (B3LYP and PBE0) calculations fail in reproducing the pattern of the carbon-carbon bond lengths obtained with the CCSD, QCISD, and MP4 methods. The structures of three shorter n-even chains, C(n)H(-) (n = 2, 4, and 6), along with those of four n-odd compounds (n = 3, 5, 7, and 9) are also investigated at the CCSD/6-311+G(2d,p) level of theory.     

NC2S+离子的结构和稳定性的理论研究

在密度泛函和从头算理论水平下计算了单重态的NC2S+离子的结构、能量、光谱以及稳定性. 在B3LYP/6-311G(d)水平下, 得到8个异构体, 它们由15个过渡态相连接. 在CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE水平下, 得到能量最低的异构体是直线型的具有1Σ电子态的NCCS+(1)(0.0 kJ/mol), 其次是直线型的异构体CNCS+(2)(54.8 kJ/mol). 两个低能量的异构体1和2及另外一个高能量的直线型异构体CCNS+(3)(323.8 kJ/mol)都具有相当大的动力学稳定性, 这三个异构体在具备一定条件的实验室和星际条件下是可以进行观测的. 分析了这3个异构体的成键性质.     

Theoretical study on stability and properties of NC2O isomers

The structures, energetics, spectroscopies, and stabilities of the doublet NC(2)O radical are explored at density functional theory and ab initio levels. Nine minimum isomers are located connected by 22 interconversion transition states. At the CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE level, the lowest-lying isomer is bent NCCO 1 (0.0 kcal/mol) with (2)A' state followed by bent isomer CNCO 2 (16.7). Two isomers (1 and 2) and another high-lying species CCNO 4 (99.4) with bent structure are considerably stabilized by a barrier of at least 20 kcal/mol. All of the three isomers should be experimentally or astrophysically observable. This result is consistent with their indication of neutralization-reionization mass spectrometry experiments. Also, the calculated spectroscopic properties and bond distances of known NCCO 1 are consistent with recent experimental observations and theoretical studies. The bonding natures of the isomers 1, 2, and 4 are analyzed. Their molecular properties including the heats of formation, adiabatic ionization potentials, and adiabatic electronic affinities are calculated at the higher levels G3//B3LYP, G3(MP2)//B3LYP, QCISD, and CCSD(T) (single-point). Possible formation strategies of the isomers 1, 2, and 4 in laboratory and space are also discussed in detail.     

Ring inversion, structural stability and vibrational assignments of sulfolane c-C4H8SO2 and 3-sulfolene c-C4H6SO2

The structural stability of sulfolane (tetrahydrothiophene1,1-dioxide) and 3-sulfolene (dihydrothiophene1,1-dioxide) was investigated by DFT-B3LYP and ab initio MP2 calculations with 6-311+G**) basis set. The calculated symmetric ring-puckering potential of 3-sulfolene at the B3LYP level is consistent with a flat minimum that corresponds to a planar ring but at the MP2 level with a double minimum with a low barrier of about 193calmol(-1) to ring planarity in reasonable agreement with experimental results. From the calculations at the two levels of theory sulfolane was predicted to exist predominantly in the twist conformation. The vibrational wavenumbers were calculated at the MP2/6-31G** level of theory and the potential energy distributions PED among the symmetry coordinates of the normal modes were computed for the low-energy structure of the molecules. Complete vibrational assignments were provided on the basis of the calculated PED values. The experimental infrared and Raman spectra of the two molecules were compared to the calculated ones.     

HBO2异构体的结构和相对稳定性

用ab initio方法在MP2/6 311＋＋G(d,p)水平下优化得到了HBO2体系的若干异构体和过渡态，并在QCISD(t)/6 311＋＋G(3df,2p)//MP2/6 311＋＋G(d,p)水平下进行了单点能量校正.对计算结果的分析表明，无论是在热力学还是在动力学上，具有链状结构的HOBO异构体(E1)是势能面上最稳定的结构，并对E1的电子结构进行了分析；另一具有C2v对称性的HBO(O)结构的异构体(E2)的能量比E1高381.72 kJ•mol－1，由于E2处于一个较深的势垒中，因此是比较稳定的，可以推断，在适合的实验中应该可以观察到异构体E2.     

Gold apes hydrogen. The structure and bonding in the planar B7Au2- and B7Au2 clusters

We produced the B7Au2- mixed cluster and studied its electronic structure and chemical bonding using photoelectron spectroscopy and ab initio calculations. The photoelectron spectra of B7Au2- were observed to be relatively simple with vibrational resolution, in contrast to the complicated spectra observed for pure B7-, which had contributions from three isomers (Alexandrova et al. J. Phys. Chem. A 2004, 108, 3509). Theoretical calculations show that B7Au2- possesses an extremely stable planar structure, identical to that of B7H2-, demonstrating that Au mimics H in its bonding to boron, analogous to the Au-Si bonding. The ground-state structure of B7Au2- (B7H2-) can be viewed as adding two Au (H) atoms to the terminal B atoms of a higher-lying planar isomer of B7-. The bonding and stability in the planar B7Au2- (B7H2-) clusters are elucidated on the basis of the strong covalent B-Au (H) bonding and the concepts of aromaticity/antiaromaticity in these systems.     

A photoelectron spectroscopy and ab initio study of B21-: negatively charged boron clusters continue to be planar at 21

The structures and chemical bonding of the B(21)(-) cluster have been investigated by a combined photoelectron spectroscopy and ab initio study. The photoelectron spectrum at 193 nm revealed a very high adiabatic electron binding energy of 4.38 eV for B(21)(-) and a congested spectral pattern. Extensive global minimum searches were conducted using two different methods, followed by high-level calculations of the low-lying isomers. The global minimum of B(21)(-) was found to be a quasiplanar structure with the next low-lying planar isomer only 1.9 kcal/mol higher in energy at the CCSD(T)/6-311-G* level of theory. The calculated vertical detachment energies for the two isomers were found to be in good agreement with the experimental spectrum, suggesting that they were both present experimentally and contributed to the observed spectrum. Chemical bonding analyses showed that both isomers consist of a 14-atom periphery, which is bonded by classical two-center two-electron bonds, and seven interior atoms in the planar structures. A localized two-center two-electron bond is found in the interior of the two planar isomers, in addition to delocalized multi-center σ and π bonds. The structures and the delocalized bonding of the two lowest lying isomers of B(21)(-) were found to be similar to those in the two lowest energy isomers in B(19)(-).     

Investigation on Structure and Bonding of As_4S_4 Isomers

Mineral medicine, especially those containing heavy metals, is one of the characteristics of traditional Chinese medicine. A famous mineral medicine, realgar, containing heavy metal arsenic with a chemical formula of As_4S_4, has the function of detoxification, killing bacteria and viruses, and eliminating dampness and phlegm. Different As_4S_4 isomers are likely to have different drug effects and pharmacological actions. Therefore, it is of great scientific significance to find more stable As_4S_4 isomers. In view of this, ab initio molecular orbital theory and density functional theory(DFT) have been used to study ten isomers of As_4S_4 at the B3LYP/6-31 G*, B3LYP/6-311+G*, B3LYP/6-311+G(3 df, 2 p) and MP2/(6-311+G*, LanL2 MB) levels of theory. In addition to the two isomers having been studied previously, eight new isomers were investigated in the present paper. All the ten As_4S_4 isomers were proved to be true local minima on their potential energy surfaces. The calculated NICS values and molecular orbital analyses showed that, the D_(2d) symmetric As_4S_4, isomer 1, may be s-aromatic. The study proves that ten As_4S_4 isomers are stable thermodynamically, and are highly desirable for the future theoretical study of realgar.     

Peroxide linkage N2O4 molecule: Prediction of new ONOONO isomers

Ab initio and density functional theory (DFT) have been applied to predict the existence of the ONOONO molecules. Two new isomers, cis‐perp‐cis and cis‐perp‐trans ONOONO, are found to be minima on the potential energy hypersurface calculated at MP2 level with 6‐311G(d) basis set and Becke3P86 and Becke3LYP levels with the cc‐PVDZ basis set. Vibrational frequencies are calculated at the optimized structures. Calculations at both ab initio and DFT methods show that cis‐perp‐cis is energetically favored in all three O O bonded ONOONO isomers. © 1999 John Wiley & Sons, Inc. Int J Quant Chem 76: 77–82, 2000     

An ab initio correlated study of the potential energy surface for the HOBr.H2O complex

The potential energy surface (PES) for the HOBr.H(2)O complex has been investigated using second- and fourth-order M?ller-Plesset perturbation theory (MP2, MP4) and coupled cluster theory with single and doubles excitations (CCSD), and a perturbative approximation of triple excitations (CCSD-T), correlated ab initio levels of theory employing basis sets of triple zeta quality with polarization and diffuse functions up to the 6-311++G(3dp,3df ) standard Pople's basis set. Six stationary points being three minima, two first-order transition state (TS) structures and one second-order TS were located on the PES. The global minimum syn and the anti equilibrium structure are virtually degenerated [DeltaE(ele-nuc) approximately 0.3 kcal mol(-1), CCSD-T/6-311++G(3df,3pd) value], with the third minima being approximately 4 kcal mol(-1) away. IRC analysis was performed to confirm the correct connectivity of the two first-order TS structures. The CCSD-T/6-311++G(3df,3pd)//MP2/6-311G(d,p) barrier for the syn<-->anti interconversion is 0.3 kcal mol(-1), indicating that a mixture of the syn and anti forms of the HOBr.H(2)O complex is likely to exist.     

Ab initio molecular-orbital study of structures and energetics of Si3H3 neutral and anion

The geometric structures and isomeric stabilities of various stationary points in Si(3)H(3) neutral and its anion are investigated at the coupled-cluster singles, doubles (triples) [CCSD(T)] level of theory. For geometrical surveys, the basis sets used are of the Dunning's correlation consistent basis sets of triple-zeta quality for the neutral. To the anions, the Dunning's correlation consistent basis sets of double-zeta quality with diffuse functions are applied. For the three lower-lying anion isomers, the Dunning's correlation consistent basis sets of triple-zeta quality with diffuse functions (aug-cc-pVTZ) are also used. The final energies for the optimized stationary points are calculated at the CCSD(T) level of theory with the aug-cc-pVTZ basis sets. The basis sets of 6-311++G(3df,2pd) were also used for the lower-lying anion isomers. The Gaussian-2 method was performed only for the lower-lying anion isomers to clarify the relative stabilities. The global minimum neutral 1 (C(1):(2)A) has an unsymmetrical hydrogen-bridged bond; the conformer 2 in C(s) symmetry is a saddle point connecting the two equivalent isomers 1. Two lower-lying isomers (3 and 4) are also predicted within the energy range of 20 kJmol. In the anion, however, the conformer 4 (C(s):(1)A(')) with five formal valence electrons is a global minimum. Two more isomers (2 and 3) lie within 20 kJmol as in the neutral; the conformer 1 converts to the isomer 2. The quartets for the neutrals and diradical triplets for the anions were further studied; lower-lying quartets and triplets, competing with the corresponding doublet and singlet, respectively, were not found in the present systems. The vertical and adiabatic electron affinities of the global minimum neutral 1, producing the second lowest-lying anion isomer 2, amount to 2.18 and 2.35 eV, respectively, at the CCSD(T)/aug-cc-pVTZ level of theory. The electron addition to the third lowest-lying neutral isomer 4 produces the largest vertical electron affinities of 2.48 eV. The D(3h) structure, being the global minimum in the corresponding Si(3)H(3) (+) cation (trisilacyclopropenyl cation), converts to the isomer 8 (C(s)) or 11 (C(2)) due to the Jahn-Teller effect in the Si(3)H(3) neutral.     

Density functional theory and post-Hartree-Fock studies on molecular structure and harmonic vibrational spectrum of formaldehyde

Density functional theory (DFT) with the Becke's three-parameter exchange correlation functional and the functional of Lee, Yang and Parr, gradient-corrected functionals of Perdew, and Perdew and Wang [the DFT(B3LYP), DFT(B3P86) and DFT(B3PW91) methods, respectively], and several levels of conventional ab initio post-Hartree-Fock theory (second- and fourth-order perturbation theory M?ller-Plesset MP2 and MP4(SDTQ), coupled cluster with the single and double excitations (CCSD), and CCSD with perturbative triple excitation [CCSD(T)], configuration interaction with the single and double excitations [CISD], and quadratic configuration interaction method [QCISD(T)], using several basis sets [ranging from a simple 6-31G(d,p) basis set to a 6-311+ +G(3df, 2pd) one], were applied to study of the molecular structure (geometrical parameters, rotational constants, dipole moment) and harmonized infrared (IR) spectrum of formaldehyde (CH₂O). High-level ab initio methods CCSD(T) and QCISD(T) with the 6-311+ +G(3df, 2pd) predict correctly molecular parameters, vibrational harmonic wavenumbers and the shifts of the harmonic IR spectrum of ¹²CH₂ ¹⁶O upon isotopic substitution. Received: 30 January 1997 / Accepted: 7 May 1997     

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.     

On the vibrational spectra and structural parameters of methyl, silyl, and germyl azide from theoretical predictions and experimental data

The infrared and Raman spectra of methyl, silyl, and germyl azide (XN3 where X=CH3, SiH3 and GeH3) have been predicted from ab initio calculations with full electron correlation by second order perturbation theory (MP2) and hybrid density function theory (DFT) by the B3LYP method with a variety of basis sets. These predicted data are compared to previously reported experimental data and complete vibrational assignments are provided for all three molecules. It is shown that several of the assignments recently proposed [J. Mol. Struct. (Theochem.) 434 (1998) 1] for methyl azide are not correct. Structural parameters for CH3N3 and GeH3N3 have been obtained by combining the previously reported microwave rotational constants with the ab initio MP2/6-311+G(d,p) predicted values. These "adjusted r0" parameters have very small uncertainties of +/-0.003 A for the XH distances and a maximum of +/-0.005 A for the heavy atom distances and +/-0.5 degrees for the angles. The predicted distance for the terminal NN bond which is nearly a triple bond is much better predicted by the B3LYP calculations, whereas the fundamental frequencies are better predicted by the scaled ab initio calculations. The results are discussed and compared to those obtained for some similar molecules.     

顺、反二氟乙烯分子结构和稳定性的ab initio研究

The geometries and their harmonic frequencies of the difluoroethylene isomers were gradient optimized using ab initio methods at the different levels. High-level methods give very good structural parameters in comparison to experiment. The single energy of stable structure and the frontier orbital energy gap were calculated using MP4SDTQ/6-311++G(2df,2pd) method, and it is found that the single energy of stable structure of cis isomer is lower than that of trans isomer. This is consistent with the experimental result. Mulliken population analyses were also performed to interpret the internal nature of cis preference of difluoroethylene isomers.