Searching for stable hept-C62X2 (X = F, Cl, and Br): structures and stabilities of heptagon-containing C62 halogenated derivatives

To determine the geometries of the most stable hept-C(62)X(2) (X = F, Cl, and Br) isomers, all 967 possible hept-C(62)F(2) isomers have been orderly optimized using AM1, HF/STO-3G, B3LYP/3-21G, and B3LYP/6-31G* methods, and chlorofullerenes and bromofullerenes, which are isostructural with five most stable hept-C(62)F(2) isomers, were regarded as candidates of the most stable isomer, and optimized at the B3LYP/6-31G* level. The results reveal that 2,9- and 9,62-hept-C(62)X(2) (X = F, Cl, and Br) are the two most stable isomers with slight energy difference. The halogenation releases strain energy of hept-C(62), and all halogenated fullerenes are more chemically stable than hept-C(62) with lower E(HOMO) and higher E(LUMO). All five most stable hept-C(62)X(2) (X = F, Cl, and Br) isomers are energetically favorable, and their thermodynamic stability decreases along with the increase of sizes of addends. Only hept-C(62)F(2) isomers show high thermodynamic stability, and they are potentially synthesized in experiments. 59,62-squ-C(62)X(2) (X = F, Cl, and Br) were computed for comparison, and they are found to be more stable than their heptagon-containing isomers.     

Theoretical investigation of C56 fullerene isomers and related compounds

All the 924 classical isomers of fullerene C(56) have been investigated by PM3, and some most stable isomers are refined with HCTH/3-21G and B3LYP6-31G(d) methods. D(2):003 with the least number of adjacent pentagons is predicted to be the most stable isomer at B3LYP/6-31G(d) level, while C(s):022 and C(2):049 possess nearly degenerate energies with relative energies of 0.03 and 3.90 kcal/mol, respectively. However, as to dianionic C(56)(2-) fullerene, C(2v):011 is predicted to be the most stable isomer. Investigations also show that the encapsulation of Ca atom in C(56) fullerene is exothermic and the metallofullerenes Ca@C(56) can be described as Ca(2+)@C(56)(2-). The computed relative stabilities show that the D(2):003 behaves more thermodynamically stable than other isomers in a wide temperature interval, and C(2v):011 should also be an important component. The electronic isomerization of C(56) (C(2v):011) and C(50) (D(5h):002) indicates that this phenomenon might be rather general in fullerenes and causes different properties, thus bringing about new possible applications of fullerenes. The static second-order hyperpolarizabilities of the three most stable isomers are slightly larger than that of C(60).     

Stability of Fullerene Cages： C30 to C70 and Their Dependence on Shared Pentagon Bonds

1 INTRODUCTION Since the discovery of buckministerfullerene (C60)[1], experimental and theoretical studies of the structure and properties of fullerenes have spread worldwide. Many experimental studies appeared on their syntheses[2], isolation[3], and characterization[4]. Theoretically, except studies on their chemical pro- perties[5], most of the work was concentrated on their interesting geometric[6～8] and electronic structures[9, 10]. For example, the ground state of C48 has C2 symme…     

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…     

Structures, stabilities, electronic, and optical properties of C64 fullerene isomers, anions (C64(2-) and C64(4-)), metallofullerene Sc2@C64, and Sc2C2@C64

The 3465 classical isomers of C(64) fullerene have been investigated by quantum chemical methods PM3, and the most stable isomers have been refined with HCTH/3-21G//SVWN/STO-3G, B3LYP/6-31G(d)//HCTH/3-21G, and B3LYP/6-311G(d)//B3LYP/6-31G(d) level. C(64)(D(2):0003) with the lowest e(55) (e(55) = 2), the number of pentagon-pentagon fusions, is predicted to be the most stable isomer and it is followed by the C(64)(C(s):0077) and C(64)(C(2):0103) isomers within relative energy of 20.0 kcal/mol. C(64)(D(2):0003) prevails in a wide temperature range according to energy analysis with entropy contribution at B3LYP/6-31G(d) level. The simulated IR spectra and electronic spectra help to identify different fullerene isomers. All the hexagons in the isomers with e(55) = 2 display local aromaticity. The relative stabilities of C(64) isomers change with charging in ionic states. Doping also affects the relative stabilities of fullerene isomers as demonstrated by Sc(2)@C(64)(D(2):0003) and Sc(2)@C(64)(C(s):0077). The bonding of Sc atoms with C(64) elongates the C-C bond of two adjacent pentagons and enhances the local aromaticity of the fullerene cages. Charging, doping, and derativization can be utilized to isolate C(64) isomers through differentiating the electronic and steric effects.     

Prediction of low-energy isomers of large fullerenes from C132 to C160

To predict energetically favored isomers, we used a topological scheme as a prescreening tool to select candidate isomers for each fullerene from C(106) to C(160). Comparison with the PM3 and tight-binding (TB) potential calculated results and few published data for the low-energy isomers of C(106) to C(130) indicates that the prescreening approach is feasible. For each fullerene from C(132) up to C(160), the selected 1000 candidate isomers were further optimized by PM3 and TB potential. The analysis of the semiempirical PM3 and TB results of C(106) to C(160) provides some qualitative features of the large fullerenes. Furthermore, calculations at the B3LYP/6-31G*//B3LYP/3-21G level of theory were carried out on the top ten PM3 and TB low-energy isomers of C(132) to C(160) to accurately predict the stable isomers, and the HOMO-LUMO gap, the ionization energy, and electron affinity of the lowest-energy isomers were also investigated at the same level.     

N5H5异构体的结构与稳定性的理论研究

采用密度泛函理论的B3LYP方法在6-311++G**基组水平上对N5H5氮氢化合物异构体可能存在的构型进行了几何优化, 得到23种稳定异构体, 并研究了这些异构体间可能的互变异构情况. 为了讨论N5H5异构体作为含能材料候选物质的可能性, 还采用了G3B3方法计算了能量, 并且计算了异构体的生成热(⊿Hf,298).结果表明, 在23种异构体中链状异构体最稳定, 四元环四氮烷异构体最不稳定, 存在一个N=N双键的异构体较同类异构体能量低, 较为稳定; N5H5氮氢化合物的生成热均为正, 其中异构体E1生成热最高. 估算了N5H5的摩尔体积, 由密度公式ρ=MT/Vmol,得到E1 的密度最大.     

C50(D5h)衍生物-异质富勒烯C48P2的分子行为理论研究

用半经验的AM1和MNDO方法以及密度泛函B3LYP/3-21G方法对C₅₀(D_5h)的衍生物C₄₈P₂的所有可能的异构体进行了系统的理论研究. 优化了稳定构型, 计算了生成热、HOMO-LUMO能级差、NICS值、红外光谱及电子光谱, 并与C₄₈X₂(X=B, N)的分子行为进行了比较. 计算结果表明: (1) C₄₈P₂的最稳定异构体是异构体C₄₈P₂-78, 该异构体对应于赤道位置六元环内的1,4-取代产物; (2) 决定C₄₈P₂异构体稳定性的主要因素是碳笼的张力, 而稳定性和芳香性之间没有明显的相关性; (3) 相对较稳定的C₄₈P₂异构体的HOMO-LUMO能级差比C₅₀的HOMO-LUMO能级差大; (4) 计算出的红外光谱以及电子光谱可以供实验参考; 计算出的NICS 值也可以用来表征C₄₈P₂各异构体. (5) C₄₈P₂与C₄₈X₂(X=B, N)具有相同的取代选择性规律及稳定性决定因素, 并且相对较稳定的异构体均具有比C₅₀本体大的HOMO-LUMO能级差.     

The C-H and alpha(C-X) bond dissociation enthalpies of toluene, C6H5-CH2X (X = F, Cl), and their substituted derivatives: a DFT study

The homolytic C-H bond dissociation enthalpies (BDEs) of toluene and its para- and meta-substituted derivatives have been estimated by using the (RO)B3LYP/6-311++G(2df,2p)//(U)B3LYP/6-311G(d,p) procedure. The performance of two other hybrid functionals of DFT, namely, B3PWP91 and O3LYP, has also been evaluated using the same basis sets and molecules. Our computed results are compared with the available experimental values and are found to be in good agreement. The (RO)B3LYP and (RO)O3LYP procedures are found to produce reliable BDEs for the C-H bonds in toluene and the C-X (X = F, Cl) bond in alpha-substituted toluene (C6H5-CH2X) and their substituted derivatives. The substituent effect on the BDE values has been analyzed in terms of the ground-state effect and the radical effect. The effect of polarization of the C-H bond on the substituent effect is also analyzed. The BDE(C-H) and BDE(C-X) values for alpha-substituted (X = F and Cl) toluenes with a set of para substituents are presented for the first time.     

Fused five-membered rings determine the stability of C60F60

The structure and stability of a set of (CF)60 isomers have been computed at the B3LYP/6-31G(d) density functional theory level. The most stable isomer (6, F4@C60F56) has tube-like structure with four endo C-F bonds and fused five-membered rings at the end of the tube, while the reported most stable cage structure (2, F8@C60F52) with eight endo C-F bonds is higher in energy by 22.6 kcal/mol. This is in contrast to the isolated pentagon rule for the stability of fullerenes. The mean bond dissociation energy of 6 is larger than those of the experimental known C60F36, C60F48, and graphite fluoride. The relative energy per CF unit of 6 to graphite fluoride (CF)n is 3.7 kcal/mol, which is smaller than that of C60 fullerene per carbon to graphite (about 9-10 kcal/mol).     

Structures, stabilities, and electronic and optical properties of c(58) fullerene isomers, ions, and metallofullerenes.

The 1205 classical isomers of fullerene C58, as well as one quasi-fullerene C58 isomer with a heptagonal ring (labeled as Cs:hept) have been investigated by the quantum chemical methods PM3, HCTH/3-21G, and B3LYP/6-31G(d). Isomer C3v:0001, which has the lowest number of adjacent pentagons, is predicted to be the most stable isomer, but the quasi-fullerene isomer Cs:hept is only 2.50 kcal mol-1 higher in energy. Systematic investigations of the electronic properties of C3v:0001 and Cs:hept find that the C3v:0001 isomer has high vertical electron affinity (3.19 eV). The nucleus-independent chemical shifts (NICS) value at the center of Cs:hept (-5.1 ppm) is more negative than that of C60 (-2.8 ppm). The NICS value at the center of the heptagonal ring in Cs:hept (-2.5 ppm) indicates weakly aromatic character. In contrast, the C58(6-) and C58(8-) ions of the C3v:0001 and Cs:hept geometries possess large aromatic character, with NICS values between -14.0 and -26.2 ppm. To clarify the thermodynamic stabilities of C58 isomers at different temperatures, the entropy contributions are taken into account on the basis of the Gibbs energy at the B3LYP/6-31G(d) level. The C3v:0001 isomer prevails in a wide range of temperatures, and the Cs:hept isomer is also an important component around 2800 K. The IR spectra of C58 isomers are simulated to facilitate experimental identification of different isomers. In addition, the electronic spectra and the second-order hyperpolarizabilities are predicted by ZINDO and the sum-over-states model. The static second-order hyperpolarizability of the C3v:0001 isomer is 96.5 % larger than that of C60, and its second-order hyperpolarizabilities at external field frequencies are at least nine times larger than those of C60.     

Ab initio and DFT energetics of silylenic X–CNSi (X=H, F, Cl, and Br)

Singlet–triplet energy splitting for 24 silylenic reactive intermediates, X–CNSi (where X=H, F, Cl and Br), are compared and contrasted at 11 levels of theory: B1LYP/6-31++G**, B3LYP/6-31++G**, B1LYP/6-311++G**, B3LYP/6-311++G**, MP3/6-31G*, MP3/6-311++G**, MP2/6-31+G**, MP2/6-311++G**, MP4 (SDTQ)/6-311++G**, QCISD(T)/6-311++G** and CCSD(T)/6-311++G**. Each X-substituted silylenic species may either be singlet (s) or triplet (t), with one of the following three structures: 3-X-2-aza-1-silacyclopropenylidene (1_s-X, 1_t-X); [(X-imino)methylene]silylene (2_s-X, 2_t-X); and X-cyanosilylene (3_s-X, 3_t-X). For all X–CNSi species studied, orders of singlet–triplet energy separations (ΔE_s-t,X), appear as a function of electro-negativity (F>Cl>Br>H). For the six H–CNSi isomers (X=H), stability order is: 3_s-H>1_s-H>2_t-H>3_t-H>2_s-H>1_t-H. Likewise, stability order for the six isomers with X=F, is: 3_s-F>3_t-F>1_s-F>1_t-F>2_s-F>2_t-F. For X=Cl, the order of stability is: 3_s-Cl>1_s-Cl>3_t-Cl>2_t-Cl>1_t-Cl>2_t-Cl. Finally, the order of stability for six isomers of Br–CNSi is: 3_s-Br>3_t-Br>1_s-Br>2_s-Br>2_t-Br>1_t-Br. The lowest energy minimum, among all 24 species scrutinized, appears to be the singlet acyclic 3_s-X. Triplet silylene 2_t-H is suggested to be more stable than its corresponding 2_s-H at MP3, MP2 and DFT levels of theory. Comparisons between relative stabilities; multiplicities and geometrical parameters of 1–3 are discussed.     

Structures, stabilities, and electronic and optical properties of C62 fullerene isomers

The 2385 classical isomers and four nonclassical isomers of fullerene C62 have been studied by PM3, HCTH/3-21G//SVWN/STO-3G, B3LYP/6-31G(d)//HCTH/3-21G, and B3LYP/6-31G(d)//B3LYP/6-31G(d). The Cs:7mbr isomer, with a chain of four adjacent pentagons surrounding a heptagon, is predicted to be the most stable isomer, followed by C2v:4mbr which is 3.15 kcal/mol higher in energy. C2:0032 with three pairs of adjacent pentagons is the most stable isomer in the classical framework. To clarify the relative stabilities of C62 isomers at high temperatures, the entropy contributions are taken into account on the basis of the Gibbs energy at the B3LYP/6-31G(d) level. Analyses reveal that Cs:7mbr prevails in a wide temperature range. The vibrational frequencies of the five most stable C62 fullerene isomers are also predicted at the B3LYP/6-31G(d) level, and the simulated IR spectra show important differences in positions and intensities of the vibrational modes for different isomers. The nucleus-independent chemical shift and the density of states of the three most stable isomers show that the square in C2v:4mbr and the adjacent pentagons in Cs:7mbr and C2:0032 possess high chemical reactivity. In addition, the electronic spectra and second-order hyperpolarizabilities are determined by means of ZINDO and the sum-over-states mode. The intensity-dependent refractive index gamma(-omega; omega, omega, -omega) at omega = 2.3305 eV of Cs:7mbr is very large because of resonance with the external field. The second-order hyperpolarizabilities of the five most stable isomers of C62 are predicted to be larger than those of C60.     

Theoretical studies of the substitution patterns of boron-nitrogen (BN) fullerenes: from C50 up to C20B15N15 CBN ball

Study on the patterns of successive BN pair substitution in C50 fullerene and the chemical and electronic properties of these substitutionally doped heterofullerenes has been carried out with semiempirical (AM1 and MNDO) and density functional (B3LYP/3-21G) theories. The BN units prefer to stay together following "single bond", "hexagon filling", and "continuity and equatorial belt" rules. The driving force governing the stabilities of these BN-doped fullerenes is the strain of the cage. Compared with C50, the redox activity of C50-2x(BN)x (x = 1-15) isomers decreases and becomes weaker by increasing the number of BN units, while the aromaticity of the C50-2x(BN)x derivatives decreases and is independent of the number of BN units but related to the substitution positions. The main infrared absorptions are similar for all the C50-2x(BN)x isomers and the infrared spectrum becomes simpler and stronger with increasing the number of BN groups.     

Theoretical studies of structures and stabilities of a new odd-numbered fullerene dimer: C141

The possible isomers of a newly synthesized C(141) molecule are calculated using MNDO, AM1, PM3, B3LYP/3-21G, and B3LYP/6-31G(d) methods. The geometry optimizations showed that the isomer 8-8 has the lowest total energy in all 64 possible structures of C(141). Unlike those of C(130), C(140), etc., the C(141) 8-8 shows a new structure: two C(70) side cages open [6.6] ring junctions located at the equator (instead of cap) area to create new chemical bonds for the bridge atom. Theoretical measurements of the average length of the long and short axes of C(70) side cages in the C(141) molecule reveal that when two C(70) cages are connected with each other at the equators, their geometric shapes become more spherical compared with the pristine C(70); this leads to a reduction of the molecular polarizability. Analysis of the local and global strain indicates that the global strain of C(70) monomer in the C(141) 8-8 is greatly reduced compared to the pristine C(70). The stable C(70) derivatives that are formed with reacted C-C bonds in the equator area may put new insights into fullerene chemistry, in particular, for C(70) to react with a large molecule. The results are discussed together with the experimental data.     

Violating the isolated pentagon rule (IPR): endohedral non-IPR C98 cages of Gd2@C98

The geometric, electronic structure, and thermodynamic stability of large gadolinium-containing endohedral metallofullerenes, Gd(2)@C(98), have been systematically investigated by comprehensive density functional theory calculations combined with statistical mechanics treatments. The Gd(2)@C(2)(230924)-C(98) structure, which satisfies the isolated-pentagon rule (IPR), is determined to possess the lowest energy followed with some stable non-IPR isomers. In order to clarify the relative stabilities at elevated temperatures, entropy contributions are taken into account on the basis of the Gibbs energy at the B3LYP level for the first time. Interestingly, a novel non-IPR Gd(2)@C(1)(168785)-C(98) isomer which has one pair of pentagon adjacency is more thermodynamically stable than the lowest energy IPR species within a wide temperature interval related to fullerene formation. Therefore, the Gd(2)@C(1)(168785)-C(98) is predicted to be the most proper isomer obtained experimentally, which is the largest non-IPR carbon cage found so far. Our findings demonstrate that interaction between metals and carbon cages could stabilize the fused pentagons effectively, and thus, the non-IPR isomers should not be ignored in some cases of endohedral metallofullerenes. The IR features of Gd(2)@C(98) are simulated to assist its future experimental characterization.     

Theoretical study on the isomeric structures and the stability of silylenoid (Tsi)Cl2SiLi (Tsi = C(SiMe3)3)

The structures and isomerization of silylenoid (Tsi)Cl(2)SiLi (Tsi = C(SiMe(3))(3)) were studied by density functional theory (DFT) at the B3LYP/6-31G(d) level. Four equilibrium structures and three isomeric transition states were located. The three-membered ring and p-complex structures, 1 and 2, are the two most stable forms. Two other local minima, the sigma-complex 3 and tetrahedron structure 4, should rearrange to 1 with very low barriers, and then to the most stable isomer 2. To exploit further the stability of silylenoid (Tsi)Cl(2)SiLi, the insertion reactions of 2 and silylene (Tsi)ClSi into the HF molecule have been investigated at the B3LYP/6-31G(d) level, respectively. The results show that the insertion of 2 into HF is very similar to that of (Tsi)ClSi into HF, but the latter is more favorable. To probe the influence of the substituent Tsi on the stability of silylenoid (Tsi)Cl(2)SiLi, the isomers and insertion reaction of silylenoid CH(3)Cl(2)SiLi were investigated in a similar way of those with (Tsi)Cl(2)SiLi. The results indicate that silylenoid containing very bulky group Tsi exhibits unusual stability because of the severe steric hindrance produced by Tsi at the center to which it is attached.     

Structure of 1-(arylselanyl)naphthalenes. 2. G dependence in 8-G-1-(p-YC(6)H(4)Se)C(10)H(6).

The structures of 8-G-1-(p-YC(6)H(4)Se)C(10)H(6) (1 (G = Cl) and 2 (G = Br): Y = H (a), OMe (b), Me (c), Cl (d), Br (e), COOEt (f), and NO(2) (g)) were investigated by X-ray crystallographic analysis, NMR spectroscopy, and ab initio MO calculations. The structures of all members in 1 and 2 are concluded to be type B, which is in striking contrast to the type A structure for 4d-g (4 (g(n)), where G = H). The Se-C(i) bond of the p-YC(6)H(4)Se group in 8-G-1-(p-YC(6)H(4)Se)C(10)H(6) is almost perpendicular to the naphthyl plane in type A, and it is located on the plane in type B. The chlorine and bromine substitution at the 8-position in 1 and 2 dramatically changes the type A structure of 4 (g(n)) to type B. The nonbonded G- - -Se-C 3c-4e type interaction must contribute to stabilize the type B structure. The type B structure in 1 and 2 should also be more stabilized than the same structure in 4 by the 3c-4e type interaction: The structure of 4b is type B in the crystals and type B would be more stable for 4c and might be for 4a in solutions. Ab initio MO calculations are performed on 8-G-1-(p-YC(6)H(4)Se)C(10)H(6), 8-G-C(10)H(6)SeH-1, and models HG- - -SeH(2), where G = Cl, Br, and F, to clarify the reason for the dramatic change in the structures. The type B structure is optimized to be more stable than the type A for all species examined, which supports the observations. The energy differences between type B and type A are larger for the models than for the naphthalenes. While the superiority of the type B for the former is Br > Cl > F, that of the latter is Br approximately Cl >/= F. These results show that the main factor of the structural change from type A to type B is the nonbonded G- - -Se-C 3c-4e interaction. The electronic effect of halogens through the naphthalene pi-framework would also contribute to some extent, although the direct comparison of the evaluated values between the naphthalene systems and the models is not so easy. Factors to stabilize the two structures of 1, 2, 4, and 8-(MeSe)-1-(p-YC(6)H(4)Se)C(10)H(6) are reexamined from a viewpoint of the nonbonded G- - -Se-C 3c-4e interaction (G dependence), together with the electronic effect of Y (Y dependence).     

Low-valent group-13 chemistry. Theoretical investigation of the structures and relative stabilities of donor-acceptor complexes R(3)E[bond]E'R' and their isomers R(2)E[bond]E'RR'

The results of quantum chemical calculations at the gradient-corrected density functional theory (DFT) level with the B3LYP functional of the donor-acceptor complexes R(3)E[bond]E'R' and their isomers R(2)E[bond]E'RR', where E, E' = B[bond]Tl and R, R' = H, Cl, or CH(3), are reported. The theoretically predicted energy differences between the donor-acceptor form R(3)E[bond]E'R' and the classical isomer R(2)E[bond]E'RR' and the bond dissociation energies of the E[bond]E' bonds are given. The results are discussed in order to show which factors stabilize the isomers R(3)E[bond]E'R'. There is no simple correlation of the nature of the group-13 elements E, E' and the substituents R, R' with the stability of the complexes. The isomers R(3)E[bond]'R' come stabilized by pi donor groups R', while the substituents R may either be sigma- or pi-bonded groups. Calculations of Cl(3)B[bond]BR' [R' = Cl, cyclopentadienyl (Cp), or Cp*] indicate that the Cp* group has a particularly strong effect on the complex form. The calculations show that the experimentally known complex Cl(3)B[bond]BCp* is the strongest bonded donor-acceptor complex of main-group elements that has been synthesized until now. The theoretically predicted B[bond]B bond energy is D(o) = 50.6 kcal/mol. However, the calculations indicate that it should also be possible to isolate donor-acceptor complexes R(3)E[bond]E'R' where R' is a sigma-bonded bulky substituent. Possible candidates that are suggested for synthetic work are the borane complexes (C(6)F(5))(3)B[bond]E'R' and (t)Bu(3)B[bond]E'R' (E' = Al[bond]Tl) and the alane complexes Cl(3)Al[bond]E'R' (E' = Ga[bond]Tl).     

C40异构体的结构和稳定性的理论研究

利用Gaussian98程序,采用密度泛函(DFT)方法中的B3LYP,选用6-31G基组对富勒烯(Fullerene)C40的6种异构体[D5d,Td,D2h,C3v,D2(Ⅰ),D2(Ⅱ)]进行了几何构型优化,其中,对于Td对称性的C40由于易发生Jahn-Teller畸变,则降低其对称性为D2d,再进行优化.对它们的平衡几何和电子结构进行了比较具体的分析,同时,根据计算得到的总能量推断出这6种异构体的稳定性顺序是D2(Ⅰ)>D5d>Td>C3v>D2h>D2(Ⅱ).