Theoretical Study on the Stable Structures and Aromaticities for Pnictogen Dianionic Sb4^2-, Bi4^2-, and （SbBi）2^2- Clusters

Through the theoretical calculation of structural optimization, vibrational frequencies and atomization energies with one method of density functional theory （B3LYP） and two post- Hartree-Fock approaches （MP2, CCSD（T））, several stable isomers for new three pnictogen dianionic Sb4^2-, Bi4^2-, and （SbBi）2^2- species were determined. For two homoatomic Sb4^2- and Bi4^2- species, there are three stable isomers： square （D4h）, roof-shaped （C2v-1）, and C2v-2 structure with the square isomer being the ground state. For the heteroatomic dian- ionic （SbBi）2^2- species, there are also three stable isomers： rhombus （D2h）, roof-shaped （C1）, and C2v structures with the rhombic isomer being the ground state. The calculated NICS values show that nucleus-independent chemical shifts （NICS） values of roof-shaped isomers for Sb4^2-, Bi4^2-, and （SbBi）2^2- species are all negative, consequently indicating that these roof-shaped isomers possess aromaticities. NICS values for the planar ring isomers are all positive, suggesting that these three planar ring isomers have antiaromatic characters. The aromaticity for the two stable roof-shaped and square isomers are preliminarily explained and discussed with MO analysis.     

Theoretical Study on the Electronic Structures of Small TinNm （n ＋ m = 5, 6） Clusters

45 isomers of TinNm （n ＋ m = 5, 6） clusters, including linear, some planar and some stero configurations, have been predicted by density functional theory method. For five-atom clusters Ti3N2 and Ti2N3, the most stable structures are trigonal bipyramid in D3h symmetry, and for TiaN cluster, the isomer with one nitrogen atom occupying the center of quasi-tetrahedron is the most stable. In the isomers of Ti4N2 and Ti3N3, the planar networks are more stable, but for Ti2N4, the six-membered ring configuration is the most favorable. Most linear structures can form weak-strong bonds alternately with higher energy. As regards to planar structures, the more Ti-N bonds are formed, the more stable they will be; for stero closed polyhedral isomers, their energies are lower.     

B20N20团簇的结构与稳定性

B3LYP/6-31G* density functional theory calculations have been carried out on the structure and stability of ten B20N20 clusters. It was found that two new proposed isomers with two octagons, twelve hexagons, eight squares in Cab and C2 symmetry were more stable than the isomer with sixteen hexagons and six squares in C2 symmetry which was previously deemed to the most stable by 79.5 and 13.8 kJ/mol respectively. The isomer with two decagons in S10 symmetry is much higher in energy than the most stable structure in C4h symmetry by 637.2 kJ/mol.     

Density Functional Theory Study on the Adsorption of HCNH and CNH2 on Cu（100） Surface

The HCNH and CNH2 adsorption on different coordination sites of Cu（100） was theoretically studied considering the cluster approach. The present calculations show that the bridge site is the most favorite for CNH2 perpendicularly adsorbed on the Cu（100） surface via the C atom. For HCNH absorbed on the Cu（100） surface, the parallel adsorption mode with the C and N atoms nearly directly above the adjacent top sites of Cu（100） surface is the most favored. Both CNH2 and HCNH are strongly bound to the Cu（100） surface with CNH2 which is lightly stable （2.51 kJ·mol^-1）, indicating that both species may be co-adsorbed on the Cu（100） surface.     

Structures and stabilities of HPS_2 isomers

The potential energy surface of HPS2 system containing nine isomers and fifteen transition states is obtained at MP2/6-311++G(d, p) and QCISD(t)/6-311++G(3df, 2p)(single-point) levels. On the potential energy surface, the lowest-lying frans-HSPS(EI) is found to be thermodynami-cally the most stable isomer followed by cis-HSPS(E2) and HP(S)S(C2v, E3) at 3.43 and 14.17 kJ/mol higher, respectively. The computed results show that species E1, E2, E3, stereo HP(S)S(Cs, E4) with PSS three-membered ring, isomers trans-HPSS(E5) and cis-HPSS(E6) which coexist with E4 are kinetically stable isomers. The products E6 and E5 in the reaction of HP with S2 can be isomerized into higher kinetic stable isomer E4 with 65.75 and 71.73 kJ/mol reaction barrier height, respectively. The predicated results may correct the possible inaccurate conclusion in that the product was experimentally assigned as isomer cis-HPSS(E6).     

Infrared Photodissociation Spectroscopy of Ti＋（CO2）2Ar and Ti＋（CO2）n （n=3-7） Complexes

Ti＋（CO2）2Ar and Ti＋（CO2）n （n=3-7） complexes are produced by laser vaporization in a pulsed supersonic expansion. The ion complexes of interest are each mass-selected in a time- of-flight spectrometer, and studied with infrared photodissociation spectroscopy. For each complex, a sharp band in the CO stretching frequency region is observed, which confirms the formation of the OTi＋CO（CO2）~_l oxide-carbonyl species. Small OTi＋CO（CO2）~_1 complexes （n_〈5） exhibit CO stretching and antisymmetric CO2 stretching vibrational bands that are blue-shifted from those of free CO and CO2. The experimental observations indicate that the coordination number of CO and CO2 molecules around TiO＋ is five. Evidence is also observed for the presence of another electrostatic bonding Ti＋（CO2）2 structural isomer for the Ti＋（CO2）2Ar complex, which is characterized to have a bent OCO-Ti＋-OCO structure stabilized by argon coordination.     

Semi-Empirical and DFT Studies on Structures and Spectra for C78（CH2）2

Eighteen possible isomers of C78（CH2）2 weTe investigated by the INDO method. It was indicated that the most stable isomer was 42,43,62,63-C78（CH2）2, where the -CH2 groups were added to the 6/6 bonds located at the same hexagon passed by the longest axis of C78 （C2v）, to form cyclopropane structures. Based on the most stable four geometries of C78（CH2）2 optimized at B3LYP/3-21G level, the first absorptions in the electronic spectra calculated with the INDO/CIS method and the IR frequencies of the C-C bonds on the carbon cage computed using the AM1 method were blue-shifted compared with those of C78 （C2v） because of the bigger LUMO-HOMO energy gap and the less conjugated carbon cage after the addition. The chemical shifts of ^13C NMR for the carbon atoms on the added bonds calculated at B3LYP/3-21G level were moved upfield thanks to the conversion from sp^2-C to sp^3-C.     

Theoretical Prediction on the Germylenoid HB=GeLiF and Its Insertion Reaction with R-H（R = F,OH,NH_2）

In the present work we investigated a novel triplet ground-state germylenoid HB=GeLiF as well as its insertion reactions with RH(R = F,OH and NH2) using the DFT B3LYP and QCISD methods for the first time.Geometry optimization calculations show that the triplet HB=GeLiF has three equilibrium structures,in which the four-membered ring structure is the most stable with the lowest energy.All mechanisms of the three insertion reactions of germylenoid HB=GeLiF with RH(R = F,OH,and NH2) are identical to each other.Based on the calculated results,it is concluded that under the same conditions the insertion reactions should occur easily in the order of H-F > H-OH > H-NH2.In THF solvent the insertion reactions get more difficult than in the gas phase.     

Structures and stabilities of HPO_2 isomers

The potential energy surface of HPO2 system including eight isomers and twelve transition states is predicated at MP2/6-311++G(d, p) and QCISD(t)/6-311++G(3df,2p)(single-point) levels of theory. On the potential energy surface, cis-HOPO(E1) is found to be thermodynamically and kinetically most stable isomer followed by trans-HOPO(E2) and HPO(O)(C2v, E3) at 10.99 and 48.36 kJ/mol higher, respectively. Based on the potential energy surface, only E1 and E3 are thermodynamically stable isomers, and should be experimentally observable. The products cis-HPOO(E5) and frans-HPOO(E6) in the first-step reaction of HP with O2 can isomerize into isomer E1 that has higher stability. The reaction of OH with PO will directly lead to the formation of isomer E1. The computed results are well consistent with the previous experimental studies.     

Structure and stability for (AlN)_n~+ and (AlN)_n~- (n=1-15) clusters

Using density functional theory (DFT) method with 6-31G* basis set, we have carried out the optimizing calculation of geometry, vibrational frequency and thermodynamical stability for(AIN)n+ and (AIN)n- (n =1-15) clusters. Moreover, their ionic potential (IP) and electron affinity(EA) were discussed. The results show that the electrical charge condition of the cluster has a relatively great impact on the structure of the cluster and with the increase of n, this kind of impactis reduced gradually. There are no AI-AI and N-N bonds in the stable structure of (AIN)n+ or (AIN)n-, and the AI-N bond is the sole bond type. The magic number regularity of (AIN)n+, and (AIN)n- is consistent with that for (AIN)n, indicating that the structure with even n such as 2, 4,6, … is more stable. In addition, (AIN)10 has the maximal ionization power (9.14 eV) and the minimal electron affinity energy (0.19 eV), which manifests that (AIN)10 is more stable than other clusters.     

Structure and stability of B5C and C5B clusters

Geometry optimizations and vibrational frequencies of B5C and C5B clusters were calculated with the Becke-3LYP method using the 6-311+G(d) basis set and some stable configurations of B5C and C5B clusters have been found. The most stable structure of B5C is a planar six-membered ring. However, for C5B clusters, the most stable structure is linear with a boron atom in position 3. Various configurations of B5C clusters containing three-membered boron rings have predominance in energy, whereas various configurations of C5B clusters containing three-membered carbon rings are disadvantageous in energy. In B5C clusters, isomer2 can be converted into isomer1 by surmounting an energy barrier of 43.83 kJ.mol(-1). In C5B clusters, the conversions of isomer5 into isomer2 and isomer7 into isomer2 have energy barriers of 19.66 and 20.57 kJ.mol(-1), respectively.     

Predicting Potential Stable Isomers on the Singlet Surface of the [H,P,C,S] System by the MP2 and QCISD(T) Methods

A detailed singlet potential energy surface of [H,P, C,S] system is investigated by means of the MP2 and QCISD(T) methods. Eight isomers are located on the potential energy surface, and at the final QCISD(T)/6-311++G (3df,2p)//MP2/6-311++G(d,p) level with zero-point energy correction, the chainlike isomer HPCS is found to be kinetically and thermodynamically the most stable species followed by the chainlike HSCP, planar three-membered ring HC(S)P, chainlike HCPS, and stereo three-membered ring HP(C)S, which are predicted to be also kinetically stable isomers and should be experimentally observable provided that accurate experimental conditions are available. The dissociation processes from the kinetically and thermodynamically most stable species HPCS to the low-lying molecular dissociation fragments are not more favorable in energy than the isomerization process from HPCS to HSCP. Therefore, the experimental observation for potential isomer HSCP with C ≡ P triple bond is possible by means of photoisomerization technology using HPCS as precursor.     

磷炔R-C≡P(R=－BH2,－CH3,－NH2,－OH)及其异构体的稳定性

采用B3LYP和CCSD(T)方法对R－C≡P(R=－BH2, －CH3,－NH2, －OH)体系进行了理论研究.结果表明，含C≡P三键的异构体BH2－C≡P和CH3－C≡P在各自的体系中分别是热力学最稳定的结构.而在HO－C≡P和NH2－C≡P体系中，热力学最稳定的结构却是H－P=C=O和含C≡N三键的N≡C－PH2.动力学理论研究表明，没有相关实验研究的R－C≡P(R=－BH2, NH2)体系中共有5种异构体是动力学稳定的. 在HO－C≡P体系的2种动力学稳定的异构体中， H－P=C=O连接方式的异构体已被实验所证实，而另外一种HO－C≡P连接方式的异构体的动力学稳定性较高，实验中可以观察到.对于CH3C≡P体系，研究所预示的2种动力学稳定的异构体中CH3－C≡P已被实验证实，从理论上推测另一种动力学稳定性较高的异构体HC≡C－PH2在实验中也可以检测到.     

(MN)nHm(M=Ga,In; n=1-4; m=1, 2)团簇的结构与稳定性

采用密度泛函理论(DFT)的B3LYP方法, 在6-31G**和Lanl2dz水平上分别对(MN)nHm(M=Ga, In; n=1-4; m=1, 2)进行了优化和振动频率计算. 得到了上述团簇的最稳定构型、H原子的结合能以及它们的能隙. 结果表明, (MN)nH(M=Ga, In; n=1-4)的基态构型均为双重态, (MN)nH2(M=Ga, In; n=1-4)的基态构型均为单重态; 当氢的个数为1时, 加在N原子上比加在M(M=Ga, In)原子上稳定, 如有N3单元, 那么加在N3单元两侧的构型是相同的, 且它是最稳定的; 当氢的个数为2时, 除n=1外, 分别加在两个N原子上的构型是最稳定的, 如有N3单元, 那么分别加在N3单元分离最远的两个N原子的构型是最稳定的. GaNH、(GaN)3H 和InNH的结合能和能隙都很大, 说明这些团簇都有很高的稳定性.     

Structure and Electronic Properties of Fullerenes C52q+: Is C522+ an Exception to the Pentagon Adjacency Penalty Rule?

The structure, vibrational spectra and electronic properties of the neutral, singly and doubly charged C52 fullerenes were studied by means of the Hartree-Fock method and density functional theory. Different isomers were considered, in particular those with the lowest possible number (five or six) of adjacent pentagons, and an isomer with a four-atom ring. For neutral and singly charged species, the most stable isomer is that with the lowest number of adjacent pentagons, namely five. However, for C(52)2+, the most stable structure has six adjacent pentagons. This finding, which contradicts the pentagon adjacency penalty rule, is a consequence of complete filling of the HOMO pi shell and the near-perfect sphericity of the most stable isomer. The simulated vibrational spectra show important differences in the positions and intensities of the vibrations for the different isomers.     

A photoelectron spectroscopic and computational study of sodium auride clusters, NanAun- (n = 1-3)

Negatively charged sodium auride clusters, NanAun- (n = 1-3), have been investigated experimentally using photoelectron spectroscopy and ab initio calculations. Well-resolved electronic transitions were observed in the photoelectron spectra of NanAun- (n = 1-3) at several photon energies. Very large band gaps were observed in the photoelectron spectra of the anion clusters, indicating that the corresponding neutral clusters are stable closed-shell species. Calculations show that the global minimum of Na2Au2- is a quasi-linear species with Cs symmetry. A planar isomer of D2h symmetry is found to be 0.137 eV higher in energy. The two lowest energy isomers of Na3Au3- consist of three-dimensional structures of Cs symmetry. The global minimum of Na3Au3- has a bent-flake structure lying 0.077 eV below a more compact structure. The global minima of the sodium auride clusters are confirmed by the good agreement between the calculated electron detachment energies of the anions and the measured photoelectron spectra. The global minima of neutral Na2Au2 and Na3Au3 are found to possess higher symmetries with a planar four-membered ring (D2h) and a six-membered ring (D3h) structure, respectively. The chemical bonding in the sodium auride clusters is found to be highly ionic with Au acting as the electron acceptor.     

The conversion among various B4C clusters: a density functional theoretical study

Geometry optimizations and vibration frequencies of B4C clusters were performed with Becke-3LYP method using 6-31G(d) basis set. We have found 14 stable isomers, and the most stable structure among them is the five-member ring containing two three-member boron rings. We also analyzed these stable isomers in detail, and the results show that the structures containing three-member boron rings are predominant in energy for B4C clusters. In terms of MO and NBO analysis, the three-centered bond and the pi-electron delocalization play an important role in stabilizing the planar five-member rings of these B4C clusters. Our calculations suggest that isomer4 can be converted into isomer7 with only an energy barrier of 0.31 kJ mol(-1) at the B3LYP/6-311G+(3df) level. Although the planar structures of the five-member rings (isomers12-14) can be converted with each other, the conversions of isomer14 to isomer13 and isomer13 to isomer12 have high-energy barriers of 70.99 and 68.51 kJ mol(-1) at the B3LYP/6-31G(d) level, respectively.     

Low-symmetry structures of Au32Z (Z = +1, 0, -1) clusters

In this work, we have explored new stable structures of the Au32Z (Z = +1, 0, -1) clusters. Theoretical calculations using density functional theory within the generalized-gradient approximation were performed. Our results show that, in the anion state (Au32-), low-symmetry (disordered) structures are preferred over the caged fullerene-like isomer. In addition, the cationic cluster (Au32+) also exhibits a disordered low-symmetry structure as its lowest energy configuration, but it is much closer in energy to the fullerene-like isomer. These results, obtained at T = 0 K, indicate that disordered structures for the Au32- and Au32+ clusters may be detected not only at room temperature, as was experimentally verified for the Au32- one, but also at much lower temperatures.     

Theoretical study on the aromaticity of the bimetallic clusters X2M2 (X=Si, Ge, M=Al, Ga)

Four neutral bimetallic clusters X₂M₂ (X=Si, Ge, M=Al, Ga) are investigated using density functional theory (DFT) and post-HF methods. The calculated results show that each of four X₂M₂ species has two energetically close stable isomers with rhombic structure (D_2h symmetry) and trapezoidal structure (C_2v symmetry) respectively. For the Ge₂Al₂ species the rhombic (D_2h) isomer is the ground state, whereas for other three species Ge₂Ga₂, Si₂Al₂, and Si₂Ga₂, the trapezoidal (C_2v) isomers are the ground states. The calculated magnetic susceptibility anisotropy (χ_anis) and nucleus-independent chemical shift (NICS) indicate that a strong diatropic ring current exists in the two heterocyclic planar isomers, suggesting they are highly aromatic. A detailed molecular orbital analysis further reveals that both heterocyclic isomers possess multiple aromaticity derived from one delocalized π MOs and two delocalized σ MOs.