Super-valence phenomenon of carbon atoms in C20 molecule

Fourteen structures of C₂₀ are studied at DFT/B3LYP/6-31G* theoretical level. Except ring, bowl, cage and isomer 1 which have been studied before, other isomers have not been reported so far. Calculated results show that the ring has the lowest energy at this level and isomers 1, 2, 3 and 4 have lower energies than that of cage. Analyses of optimized bond lengths, electronic structure indicate that some carbon atoms express super-valence property. In addition, NICS value is consistent with molecular orbital character in denoting aromaticity of C₂₀ molecule. Delocalization character averts influence of curvature strain, which can well explain the stability of the cage.     

Comparison studies of fullerenes C72, C74, C76 and C78 by tight-binding Monte Carlo and quantum chemical methods

From C₇₂ to C₇₈, the top 20 low-energy isomers screened out from all isomers of each fullerene are optimized and computed by tight-binding Monte Carlo (TBMC), semi-empirical PM3, and ab initio B3LYP/6-31G*//HF/3-21G methods. The comparison results show that the TBMC method can efficiently optimize the structures and correctly predicate the low energy isomers. The relative energies computed by TBMC are in good agreement with the high-lever B3LYP calculation results. Our TBMC and B3LYP results show that the most energetically favorable structure of C₇₂ is not an isomer satisfying the isolated pentagon rule (IPR), which is different with the result by PM3. The symmetry of the most stable IPR isomer tends to low as the fullerene becomes large and several non-isolated-pentagon structures are found to have low symmetries and low energies close to the most stable isomer.     

Theoretical studies on the BN substituted fullerenes C70−2x(BN)x (x=1–3)—isoelectronic equivalents of C70

The equilibrium structures and relative stabilities of BN-doped fullerenes C_70−2x(BN)_x (x=1–3) have been studied at the AM1 and MNDO level. The most stable isomers of C_70−2x(BN)_x have been found out and their electronic properties have been predicted. The calculation results show that the BN substituted fullerenes C_70−2x(BN)_x have considerable stabilities, though they are less stable than their all carbon analog. For C₆₈BN, the isomers whose BN is located in the most chemically active bonds of C₇₀ (namely B and A) are among the most stable species, of which B is predicted to be the ground state. The stabilities of C₆₈BN decrease and the dipole moments increase with increasing the distance between the heteroatoms. For C₆₆(BN)₂, the lowest energy species is the isomer in which the B–N–B–N bond is formed; For C₆₄(BN)₃, the most stable species should have three BN units located in the same hexagon to form B–N–B–N–B–N ring. The ionization potentials and the affinity energies of the most stable species of BN-doped C₇₀ are almost the same as those of C₇₀ because of the isoelectronic relationship. The ionization potentials and affinity energies depend on the relative position of the heteroatoms in C₆₈BN, the chemical reactivities of the isomers whose heteroatoms are well separated should differ significantly from their all carbon analog.     

Isomers of C24. Density functional studies including gradient corrections

Density functional techniques are used to investigate the relative energies of seven different structural isomers of C₂₄. The traditional local density approximation yields the fullerene-like isomer to be the most stable. As in the case of C₂₀, the inclusion of gradient corrections has a dramatic effect on the relative energies. The gradient-corrected B-LYP method yields the monocyclic ring and graphite-like isomers to be almost isoenergetic (and most stable) while the bicyclic ring, fullerene-like, and bowl-like isomers are progressively higher in energy. The Hartree—Fock results are quite similar to the B-LYP results. Implications to fullerene growth mechanisms are pointed out.     

C60O3的结构和电子光谱的理论研究

用INDO系列方法对C₆₀O₃的可能构型进行研究,结果表明:环氧结构邻近的6-6键易发生进一步的加成反应.其中3个氧原子加在同一个六元环的6-6边上,形成环氧结构最稳定的C₃v构型,第3个氧原子加在2个环氧结构相邻的六元环的6-6边上的C₂、C_s构型也相当稳定,C₂、C_s构型的部分¹³C NMR谱与实验吻合.C₆₀O₃可能有较好的反应活性,其电子光谱属于理论预测.     

Gas-phase synthesis and the structural and electronic properties of C60S studied by mass spectrometry and molecular orbital calculation

C₆₀S⁺ was synthesized through the gas-phase ion-molecule reaction of C₆₀ with the plasmas of carbon disulfide under self-chemical-ionization (self-CI) conditions in the ion source of a mass spectrometer. Semi-empirical PM3-UHF and density functional B3LYP levels of theory with 6-31G(d) basis set calculations were performed on all the possible structures and electronic properties of the product. The results showed that the most stable structure among the possible isomers was the 6/6 closed derivative. The reaction energies of C₆₀+S⁺→C₆₀S⁺ and C₆₀+S→C₆₀S were also calculated to suggest the possibility of C₆₀S synthesis in condensed phase.     

Basis set effects on the stability of the Cl2O3 isomers using B3P86 and B3LYP methods of density functional theory

Four Cl₂O₃ isomers have been calculated using density functional theory with B3P86 and B3LYP functionals and various basis sets. The energy hypersurfaces of Cl₂O₃ are very flat and the relative energies of the isomers which have hypervalent characters such as ClOCl(O)O and ClClO₃ are strongly dependent on the basis sets. The stability for Cl₂O₃ isomers are in the order of (1)ClOOOCl(C_s), (2)ClOOOCl(C₂), (3)ClClO₃ and (4)ClOCl(O)O with ClOCl(O)O being most stable. We suggest that at least the cc-pV6Z(-ghi) basis set for Cl and the cc-pVTZ basis set for O are required to obtain reliable relative stabilities of Cl₂O₃ isomers with hypervalent characters.     

The interaction between C28 (Td) and X (X=H and F)

Hydrogen and fluorine addition reactions with C₂₈(T_d) have been investigated by the density function theory method at B3LYP/6-31G level. The interaction potential between C₂₈(T_d) and atom X (X=H and F) shows that there are three possible stable isomers of C₂₈(T_d)X (X=H and F) and the average binding energy calculations suggest that C₂₈(T_d)H₄ is the most stable hydrogen adduct among C₂₈(T_d)H_n (n=1–28). Furthermore, by comparisons of the energy between C₂₈(T_d)H and C₂₈(C_s)H we found that the former are more stable than the later, and the structural and energy analysis further indicate that C₂₈(C_s)H is only with a small distortion of C₂₈(T_d)H symmetry. In addition, the transition states, as well as reaction pathways of X transfer reactions between different key points on C₂₈(T_d) representative patch are given to explore the possible reaction mechanism.     

CH3F分子间相互作用的ab initio研究

在MP2/6-311+ +G(3d,3p)电子相关校正水平上,对CH₃F二聚体可能存在的几何构型进行全自由度能量梯度优化和频率验证,发现3种势能面上有极小点的构型.进一步在高级电子相关校正的MP4S-DTQ、CCSD(T)/6-311+ +G(3df,3pd)方法水平上,对其中总能量最小的构型进行精确计算,得到二聚体的结合能为-9.707kJ/mol.研究结果支持了由光谱实验结果推测的构型,解释了CH₃F二聚体的谐振频率的多样性.     

一些富勒烯C60双加成物稳定性的理论研究

Ｃ6 0 具有 30个等同的可参与化学反应的活泼双键 ,制备并表征富勒烯多加成产物是富勒烯化学中最前沿的课题之一 .Ｈｉｒｓｃｈ等[15 ]通过研究Ｃ6 0 亚甲基加成反应 ,提出了双加成物立体选择性的一般规律 .我们在Ｃ6 0 氧加成方面做了类似的研究工作[6 8].本文选取几类富勒烯环双加成衍生物即富勒烯氧化产物Ｃ6 0 Ｏ2 [6 ]( 1 )、富勒烯含氮衍生物Ｃ6 0 (ＮＨ) 2 [5 ]( 2 )、富勒烯含吡咯环衍生物Ｃ6 0 (ＣＨ2 ＮＨＣＨ2 ) 2 [9]( 3)和富勒醇前体Ｃ6 0 (ＳＯ4) 2 [10 ]( 4)作为模型分子进行理论研究 ,以探寻富勒烯多加成反应的一般性规律 .…     

Structure of small TiCn clusters: A theoretical study

A theoretical study of the TiC_n (n = 1–8) clusters has been carried out at the B3LYP/6-311+G(d) level. Molecular properties for three different isomers, namely linear, cyclic, and fan species, have been determined. The fan isomers, where the titanium atom is essentially side-bonded to the entire C_n unit, are predicted to be more stable than both linear and cyclic isomers. Only for the largest studied species, TiC₈, the cyclic isomer is located lower in energy. An even–odd parity effect in the incremental binding energies is observed for the three isomers, n-even species being in general more stable for linear and fan isomers, whereas for the cyclic species n-odd clusters are favoured. A topological analysis of the electronic charge density shows that all cyclic isomers correspond to true monocyclic rings, whereas for the fan species a variety of different connectivities has been observed.     

富勒烯C36及其衍生物C36H2n的量子化学研究

应用密度泛函B₃LYP方法和从头算(abinitio)HF方法,在3-21G理论水平上,对具有C_2v,D_2d,D_3h和D_6h对称性的四种C₃₆异构体以及在保持D_2d,D_3h和D_6h对称性条件下形成的不同氢化物进行了量子化学计算,研究了它们在不同量子态时的分子几何构型和电子结构.结果表明,电子相关效应在C₃₆簇化合物的电子结构中起着重要作用.它们最稳定的结构是C₃₆的³A_2u量子态D_6h异构体,而C₃₆的氢化物是在D_3h异构体的C₂位置等同碳原子上形成的¹A_1'量子态氢化物C₃₆H₁₂.     

Isomers of C20: an energy profile III

Semiempirical calculations, at the PM3 level provided within the Winmopac v2.0 software package, are used to geometrically optimize and determine the absolute energies (heats of formation) of a variety of C(20) isomers that are predicted to exist in and around the ring and cage isomers. Using the optimized Cartesian coordinates for the ring and the cage isomers, a saddle-point calculation was performed. The resulting energy profile, consisting of a series of peaks and valleys, is used as a starting point for the identification and location of fifteen additional isomers of C(20) that are predicted to be energetically stable, both via geometry optimizations and force constant analysis. These additional isomers were subsequently determined to lie adjacent to one another on the potential surface and establish a step-wise transformation between the ring and the cage. Transition-state optimization of the Cartesian coordinates at the saddle point between adjacent isomers was performed to quantify the energy of the transition state. The step-wise process from one isomer to another, which extends out over the three-dimensional surface, is predicted to require approximately 15% less energy than that of the direct, two-dimensional transformation predicted in the bowl-cage profile. However, the net atomic rearrangement for the step-wise process is about four times greater than that of the direct process. Although less in energy, the amount of atomic rearrangement in the step-wise process would make the occurrence of such a route prohibitive. Utilizing the direct distance separating the three primary isomers (ring, bowl, cage), the method of triangulation is performed to quantitatively position other C(20) structures on the potential surface, relative to the ring, bowl, and cage isomers.     

Isomers of C(20): an energy profile

Semiempirical calucaltions, at the PM 3 level, are used to geometrically optimize and determine the absolute energies (heats of formation) of a variety of C(20) isomers. Based on the geometrically optimized Cartesian coordinates of the ring and the bowl isomers, and the subsequent saddle-point calculation, a two-dimensional energy profile between these two isomers is generated. Performing geometry optimization on the Cartesian coordinates that correspond to energy minima within the ring-bowl profile, we have been able to identify several more isomers of C(20) that are predicted to be energitically stable. With these additional stable structures, we have identified pairs of isomers that lie adjacent to one another on the potential energy surface, as is evidenced by the form of their respective energy profiles. These adjacent pairs of isomers establish a step-wise transformation between the ring and the bowl. This process, which extends out over the three-dimensional surface, is predicted to require less energy than that of the direct, two-dimensional transformation predicted in the ring-bowl profile.     

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.     

Ab initio and DFT computer studies of complexes of trimethylphosphine and products of substituted phosphonium cations with hydroxide, chloride and fluoride anions: Phosphorus analogues of choline and acetylcholine

Theoretical calculations (DFT, MP2) are reported for up to four sets of reaction products of trimethylphosphine, (CH₃)₃P, each with H₂O, HCl and HF together with DFT calculations on up to three sets of reaction products of substituted phosphonium cations, (CH₃)₃P–R⁺. These products comprise (a) P(III) normal complexes (CH₃)₃PHY, (b) P(IV) ‘reverse’ complexes Y(H–CH₂)₃P–R, (c) P(IV) ylidic complexes YHCH₂(CH₃)₂P–R and (d) P(V) covalent compounds Y–P(CH₃)₃–R for Y=HO, Cl and F and R=H, CH₃, C₂H₅, C₂H₄OH and C₂H₄OC:OCH₃. Calculations are carried out at the B3LYP/6-31+G(d,p) level in all cases and also at the MP2/6-31+G(d,p) level for systems in which R=H. Minimum energy structures are determined for predicted complexes or structures and geometrical properties, harmonic vibrations and BSSE corrected binding energies are reported and compared with the limited experimental information available. Potential energy scans predict equilibria between covalent trigonal bipyramidal P(V) forms and reverse complexes comprising hydrogen bonded or ion pair, tetrahedral P(IV) forms separated by low potential energy barriers. Similar scans are also reported for equilibria between reverse complexes and ylidic complexes for Y=OH and R=CH₃, C₂H₅, C₂H₄OH and C₂H₄OC:OCH₃. Corrected binding energies, structures and values of harmonic modes are discussed in relation to bonding The names ‘pholine’ and ‘acetylpholine’ are suggested for phosphorus analogues to choline and acetylcholine.     

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.     

Sugars in the gas phase: The conformational properties of erythrose, threose, and erythrulose characterized by quantum chemistry calculations

The structural properties of the three open chain C₄H₈O₄ sugars, i.e. two aldoses (erythrose and threose) and one ketose (erythrulose), have been investigated by DFT and ab initio calculations to get accurate structures and relative energies. The structure of all the conformers predicted within 10 kJ/mol has been optimized at the B3LYP/6-311++G(d,p) level of the theory. Two types of intramolecular hydrogen bonds have been clearly identified. They are related to the hydroxyl and to the carbonyl oxygen atoms and are of weak and middle strength, respectively. The most stable structures have been optimized at the B3LYP/6-311++G(2df,p) and at the MP2/6-311++G(2df,p) levels of the theory in order to calculate accurate rotational parameters and dipole moment for their future detection in the microwave range in the gas phase. Their corresponding harmonic IR spectra have also been calculated and their fingerprint signature is discussed in the region of the OH stretching vibrations, of the torsion of the C–O bonds and of the deformation of the C–C skeleton.     

C_(2~6)烯烃在H-ZSM-5分子筛周期性模型上吸附的理论研究

采用密度泛函理论中的PBE-D方法研究了C2~6直链烯烃以及丁烯异构体在H-ZSM-5分子筛周期性模型上的吸附行为。结果表明,对于C2~6直链烯烃,随着碳数的增加,烯烃的吸附能以-12 kJ/mol的常数线性增大,体系中色散校正能ED以-13 kJ/mol的常数线性增大,而不考虑范德华力的能量E*随碳数增加变化不大,只在C5和C6时由于空间阻力凸显而减小。ED显著大于E*,说明烯烃在分子筛孔道中吸附时范德华力起主要作用,且碳数增加其影响增大。丁烯异构体吸附能大小顺序为:反式-2-丁烯顺式-2-丁烯正丁烯异丁烯。正丁烯的3种异构体ED相近且都比正丁烯的大,3种异构体的吸附能差异由E*不同引起。差分电荷密度分析表明,烯烃双键与酸性中心之间有电子聚集,且聚集程度与π配位作用相一致,而且酸性中心H原子有质子化趋向,O原子及其周围电子增加。     

Complexes of a bio-molecule and a C60 cage

C₆₀ is the most important fullerene cage and glycine is the simplest representative of a backbone unit of a protein. In this paper, the structures and the energies of glycine–C₆₀ complexes were calculated at the B3LYP/6-31G(d) level DFT. It was found that the binding of glycine to C₆₀ generated a slightly unstable complex via its amino nitrogen, a moderately unstable complex via its hydroxyl oxygen, and a very unstable complex via its carbonyl oxygen. This indicates that fullerene cages might be unable to form stable bindings to proteins via their amino nitrogen, hydroxyl oxygen and carbonyl oxygen active sites.