Electronic structures and thermochemical properties of the small silicon-doped boron clusters B(n)Si (n=1-7) and their anions

We perform a systematic investigation on small silicon-doped boron clusters B(n)Si (n=1-7) in both neutral and anionic states using density functional (DFT) and coupled-cluster (CCSD(T)) theories. The global minima of these B(n)Si(0/-) clusters are characterized together with their growth mechanisms. The planar structures are dominant for small B(n)Si clusters with n≤5. The B(6)Si molecule represents a geometrical transition with a quasi-planar geometry, and the first 3D global minimum is found for the B(7)Si cluster. The small neutral B(n)Si clusters can be formed by substituting the single boron atom of B(n+1) by silicon. The Si atom prefers the external position of the skeleton and tends to form bonds with its two neighboring B atoms. The larger B(7)Si cluster is constructed by doping Si-atoms on the symmetry axis of the B(n) host, which leads to the bonding of the silicon to the ring boron atoms through a number of hyper-coordination. Calculations of the thermochemical properties of B(n)Si(0/-) clusters, such as binding energies (BE), heats of formation at 0 K (ΔH(f)(0)) and 298 K (ΔH(f)([298])), adiabatic (ADE) and vertical (VDE) detachment energies, and dissociation energies (D(e)), are performed using the high accuracy G4 and complete basis-set extrapolation (CCSD(T)/CBS) approaches. The differences of heats of formation (at 0 K) between the G4 and CBS approaches for the B(n)Si clusters vary in the range of 0.0-4.6 kcal mol(-1). The largest difference between two approaches for ADE values is 0.15 eV. Our theoretical predictions also indicate that the species B(2)Si, B(4)Si, B(3)Si(-) and B(7)Si(-) are systems with enhanced stability, exhibiting each a double (σ and π) aromaticity. B(5)Si(-) and B(6)Si are doubly antiaromatic (σ and π) with lower stability.     

Electronic structure and thermochemical properties of small neutral and cationic lithium clusters and boron-doped lithium clusters: Li(n)(0/+) and Li(n)B(0/+) (n = 1-8)

The stability, electronic structure, and thermochemical properties of the pure Li(n) and boron-doped Li(n)B (n = 1-8) clusters in both neutral and cationic states are studied using electronic structure methods. The global equilibrium structures are established, and their heats of formation are evaluated using the G3B3 and CCSD(T)/CBS methods based on the density functional theory geometries. Theoretical adiabatic ionization energies (IE(a)) for the Li(n) clusters are in good agreement with experiment: Li(2) (G3B3, 5.21 eV; CCSD(T), 5.14 eV; expt, 5.1127 ± 0.0003 eV), Li(3) (4.16, 4.11, 4.08 ± 0.10), Li(4) (4.76, 4.68, 4.70 ± 0.05), Li(5) (4.11, 4.06, 4.02 ± 0.10), Li(6) (4.46, 4.32, 4.20 ± 0.10), Li(7) (4.07, 3.99, 3.94 ± 0.10), and Li(8) (4.49, 4.31, 4.16 ± 0.10). The Li(4) experimental IE(a) has been revised on the basis of the Franck-Condon simulations. Species Li(5)B, Li(6)B(+), Li(7)B, and Li(8)B(+) exhibit high stability as compared to their neighbors, which can be understood by considering the magic numbers of the phenomenological shell model (PSM).     

Singly and doubly lithium doped silicon clusters: geometrical and electronic structures and ionization energies

The geometric structures of neutral and cationic Si(n)Li(m)(0/+) clusters with n = 2-11 and m = 1, 2 are investigated using combined experimental and computational methods. The adiabatic ionization energy and vertical ionization energy (VIE) of Si(n)Li(m) clusters are determined using quantum chemical methods (B3LYP/6-311+G(d), G3B3, and CCSD(T)/aug-cc-pVxZ with x = D,T), whereas experimental values are derived from threshold photoionization experiments in the 4.68-6.24 eV range. Among the investigated cluster sizes, only Si(6)Li(2), Si(7)Li, Si(10)Li, and Si(11)Li have ionization thresholds below 6.24 eV and could be measured accurately. The ionization threshold and VIE obtained from the experimental photoionization efficiency curves agree well with the computed values. The growth mechanism of the lithium doped silicon clusters follows some simple rules: (1) neutral singly doped Si(n)Li clusters favor the Li atom addition on an edge or a face of the structure of the corresponding Si(n)(-) anion, while the cationic Si(n)Li(+) binds with one Si atom of the bare Si(n) cluster or adds on one of its edges, and (2) for doubly doped Si(n)Li(2)(0/+) clusters, the neutrals have the shape of the Si(n+1) counterparts with an additional Li atom added on an edge or a face of it, while the cations have both Li atoms added on edges or faces of the Si(n)(-) clusters.     

Structural and static electric response properties of highly symmetric lithiated silicon cages: Theoretical predictions

It is shown by density functional theory calculations that high symmetry silicon cages can be designed by coating with Li atoms. The resulting highly symmetric lithiated silicon cages (up to D_5d symmetry) are low‐lying true minima of the energy hypersurface with binding energies of the order of 4.6 eV per Si atom and moderate highest occupied molecular orbital–lowest unoccupied molecular orbital gaps. Moreover, relying on a systematic study of the electric response properties obtained by ab initio (Hartree–Fock, MP2, and configuration interaction singles (CIS)) and density functional (B3LYP, B2PLYP, and CAM‐B3LYP) methods, it is shown that lithium coating has a large impact on the magnitude of their second hyperpolarizabilities resulting to highly hyperpolarizable species. Such hyperpolarizable character is directly connected to the increase in the density of the low‐lying excited states triggered by the interaction between the Si cage and the surrounding Li atoms. © 2012 Wiley Periodicals, Inc.     

Theoretical study of mixed silicon-lithium clusters Si(n)Li(p)(+) (n=1-6, p=1-2)

Theoretical study on the structures of neutral and singly charged Si(n)Li(p)((+)) (n=1-6, p=1-2) clusters have been carried out in the framework of the density functional theory (DFT) with the B3LYP functional. The structures of the neutral Si(n)Li(p) and cationic Si(n)Li(p)(+) clusters are found to keep the frame of the corresponding Si(n), Li species being adsorbed at the surface. The localization of the lithium cation is not the same one as that of the neutral atom. The Li(+) ion is preferentially located on a Si atom, while the Li atom is preferentially attached at a bridge site. A clear parallelism between the structures of Si(n)Na(p) and those of Si(n)Li(p) appears. The population analysis show that the electronic structure of Si(n)Li(p) can be described as Si(n)(p)(-)+pLi(+) for the small sizes considered. Vertical and adiabatic ionization potentials, adsorption energies, as well as electric dipole moments and static dipolar polarizabilities, are calculated for each considered isomer of neutral species.     

Efficiency of random search procedures along the silicon cluster series: Si(n) (n=5-10, 15, and 20)

The efficiency of the simplest isomeric search procedure consisting in random generation of sets of atomic coordinates followed by density functional theory geometry optimization is tested on the silicon cluster series (Si(5-10, 15, 20)). Criteria such as yield, isomer distributions and recurrences are used to clearly establish the performance of the approach with respect to increasing cluster size. The elimination of unphysical candidate structures and the use of distinct box shapes and theoretical levels are also investigated. For the smaller Si(n) (n=5-10) clusters, the generation of random coordinates within a spherical box is found to offer a reasonable alternative to more complex algorithms by allowing straightforward identification of every known low-lying local minima. The simple stochastic search of larger clusters (i.e. Si(15) and Si(20)) is however complicated by the exponentially increasing number of both low- and high-lying minima leading to rather arbitrary and non-comprehensive results.     

Structures and stability of lithium monosilicide clusters SiLin (n = 4–16): What is the maximum number,magic number,and core number for lithium coordination to silicon?

In the coordination, hypervalent and cluster chemistry, three important characteristic properties are the maximum coordination number, magic number, and core coordination number. Yet, few studies have considered these three numbers at the same time for an ML(n) cluster with n larger than 8. In this article, we systematically studied the three properties of SiLi(n) (n = 4-16) clusters at the B3LYP/6-31G(d), B3LYP/6-311++G(2d), and CCSD(T)/6-311++G(3df)//B3LYP/6-311++G(2d) (for energy only) levels. Various isomeric forms with different symmetries were calculated. For each SiLi(n) (n = 4-9), silicon cohesive energy (cE) from SiLi(n) --> Si + Li(n) reaction, vertical ionization potential (vIP), and vertical electron affinity (vEA) were obtained for the lowest-energy isomer. We found that the maximum Li-coordination number of Si is 9, which is the largest number among the known MLi(n) clusters. All cE, vIP, and vEA values predicted that 6 is the magic Li-coordination number of Si. For small SiLi(n) (n < or = 6) clusters, Li atoms favor direct coordination to Si, whereas for larger SiLi(n) (n > or = 7) clusters, there is a core cluster that is surrounded by excessive Li atoms. The core Li-coordination number is 6 for SiLi(n) (n = 7,8), 7 for SiLi(n) (n = 9,10), 8 for SiLi(n) (n = 11-15) and 9 for SiLi(n) (n > or = 16). Through the calculations, we verified the relationship between the structure and stability of SiLi(n) with the maximum coordination number, magic number, and core coordination number.     

Ionized water clusters,n = 2 to 6: A high-accuracy study of structures and energetics

Ionized water clusters, , have been of remarkable interest owing to their crucial roles in many chemical and biological processes. Small cationic water clusters , n = 2 to 6 serve as reasonable models for understanding the nature of the ionized water. In this study, employing high-level ab initio quantum chemical methods, such as the density-fitted orbital-optimized linearized coupled-cluster doubles (DF-OLCCD), coupled-cluster singles and doubles (CCSD), and coupled-cluster singles and doubles with perturbative triples [CCSD(T)], a high-accuracy study of structures and energetics for cationic water clusters [, n = 2-6] is presented. In this study, 2 dimer, 8 trimer, 18 tetramer, 23 pentamer, and 25 hexamer clusters are reported. Most of the structures considered are reported for the first time. Relative, binding, and vertical attachment energies (VAEs), for the first time, are presented at the complete basis set (CBS) limit, extrapolating energies of the aug-cc-pVTZ and aug-cc-pVQZ basis sets, to provide the most accurate energetics to date. Our results demonstrate that as cluster size increases, the VAE value decreases, which indicates that large-size clusters better compensate for the electron deficiency compared with small-size clusters. The VAE values for pentamer and hexamer clusters are 118.5 to 165.5 and 121.9 to 153.7 kcal mol⁻¹, respectively. Further, our binding energy results, at the CCSD(T)/CBS level, indicate strong bindings in cationic clusters due to hydrogen bond interactions. The average binding energy per water molecule varies from −16.6 to −21.8 kcal mol⁻¹ for the clusters considered. Hence, we present the most extensive and accurate study on ionized water clusters to date. Further, our results indicate that the DF-OLCCD method is very promising for ionic molecular clusters, and its accuracy approaches the CCSD(T) quality. The inexpensive analytic gradients of DF-OLCCD, compared with CCSD(T), make it very helpful for high-accuracy studies of molecular geometries.     

Stabilisation energy of C(6)H(6)...C(6)X(6) (X = F, Cl, Br, I, CN) complexes: complete basis set limit calculations at MP2 and CCSD(T) levels

Stabilisation energies of stacked structures of C(6)H(6)...C(6)X(6) (X = F, Cl, Br, CN) complexes were determined at the CCSD(T) complete basis set (CBS) limit level. These energies were constructed from MP2/CBS stabilisation energies and a CCSD(T) correction term determined with a medium basis set (6-31G**). The former energies were extrapolated using the two-point formula of Helgaker et al. from aug-cc-pVDZ and aug-cc-pVTZ Hartree-Fock energies and MP2 correlation energies. The CCSD(T) correction term is systematically repulsive. The final CCSD(T)/CBS stabilisation energies are large, considerably larger than previously calculated and increase in the series as follows: hexafluorobenzene (6.3 kcal mol(-1)), hexachlorobenzene (8.8 kcal mol(-1)), hexabromobenzene (8.1 kcal mol(-1)) and hexacyanobenzene (11.0 kcal mol(-1)). MP2/SDD** relativistic calculations performed for all complexes mentioned and also for benzene[dot dot dot]hexaiodobenzene have clearly shown that due to relativistic effects the stabilisation energy of the hexaiodobenzene complex is lower than that of hexabromobenzene complex. The decomposition of the total interaction energy to physically defined energy components was made by using the symmetry adapted perturbation treatment (SAPT). The main stabilisation contribution for all complexes investigated is due to London dispersion energy, with the induction term being smaller. Electrostatic and induction terms which are attractive are compensated by their exchange counterparts. The stacked motif in the complexes studied is very stable and might thus be valuable as a supramolecular synthon.     

High-level ab initio predictions for the ionization energies and heats of formation of five-membered-ring molecules: thiophene, furan, pyrrole, 1,3-cyclopentadiene, and borole, C4H4X/C4H4X+ (X = S, O, NH, CH2, and BH)

The ionization energies (IEs) and heats of formation (ΔH°(f0)/ΔH°(f298)) for thiophene (C(4)H(4)S), furan (C(4)H(4)O), pyrrole (C(4)H(4)NH), 1,3-cyclopentadiene (C(4)H(4)CH(2)), and borole (C(4)H(4)BH) have been calculated by the wave function-based ab initio CCSD(T)/CBS approach, which involves the approximation to the complete basis set (CBS) limit at the coupled-cluster level with single and double excitations plus a quasi-perturbative triple excitation [CCSD(T)]. Where appropriate, the zero-point vibrational energy correction (ZPVE), the core-valence electronic correction (CV), and the scalar relativistic effect (SR) are included in these calculations. The respective CCSD(T)/CBS predictions for C(4)H(4)S, C(4)H(4)O, C(4)H(4)NH, and C(4)H(4)CH(2), being 8.888, 8.897, 8.222, and 8.582 eV, are in excellent agreement with the experimental values obtained from previous photoelectron and photoion measurements. The ΔH°(f0)/ΔH°(f298) values for the aforementioned molecules and their corresponding cations have also been predicted by the CCSD(T)/CBS method, and the results are compared with the available experimental data. The comparisons between the CCSD(T)/CBS predictions and the experimental values for C(4)H(4)S, C(4)H(4)O, C(4)H(4)NH, and C(4)H(4)CH(2) suggest that the CCSD(T)/CBS procedure is capable of predicting reliable IE values for five-membered-ring molecules with an uncertainty of ±13 meV. In view of the excellent agreements between the CCSD(T)/CBS predictions and the experimental values for C(4)H(4)S, C(4)H(4)O, C(4)H(4)NH, and C(4)H(4)CH(2), the similar CCSD(T)/CBS IE and ΔH°(f0)/ΔH°(f298) predictions for C(4)H(4)BH, whose thermochemical data are not readily available due to its reactive nature, should constitute a reliable data set. The CCSD(T)/CBS IE(C(4)H(4)BH) value is 8.868 eV, and ΔH°(f0)/ΔH°(f298) values for C(4)H(4)BH and C(4)H(4)BH(+) are 269.5/258.6 and 1125.1/1114.6 kJ/mol, respectively. The highest occupied molecular orbitals (HOMO) of C(4)H(4)S, C(4)H(4)O, C(4)H(4)NH, C(4)H(4)CH(2), and C(4)H(4)BH have also been studied by the natural bond orbital (NBO) method, and the extent of π-electron delocalization in these five-membered rings are discussed in correlation with their molecular structures and orbitals.     

Electronic structure and stability of AlnPn (n = 2-4) clusters

The geometry, electronic configurations, harmonic vibrational frequencies, and stability of the structural isomers of aluminum phosphide clusters have been investigated using the density functional theory. For dimers and trimers, the lowest energy structures are cyclic (IIs, IIIs) with D(nh) symmetry. The caged structure with Td symmetry (Xs) lie lowest in energy among the tetramers. The Al--P bond dominates the structures for many isomers so that one preferred dissociation channel is loss of the AlP monomer. The hybridization and chemical bonding in the different structures are also discussed. Comparisons with silicon and boron nitride clusters, the ground state structures of Al(n)P(n) clusters are analogous to those of their corresponding Si(2n) counterparts. This similarity follows the isoelectronic principle.     

(SiO2)n团簇的基态能量和结构

应用遗传算法对二氧化硅团簇（SiO2）n（n≤20）的结构进行了优化计算。分析讨论了结构和结合能随团簇尺寸的变化规律,发现（SiO2）n团簇系列不存在明显的幻数,并在n≤20的范围内,不呈现出相应大块物质的结构特征,     

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.     

Binding energies of small lithium clusters (Li(n)) and hydrogenated lithium clusters (Li(n)H)

Large coupled cluster computations utilizing the Dunning weighted correlation-consistent polarized core-valence (cc-pwCVXZ) hierarchy of basis sets have been conducted, resulting in a panoply of internally consistent geometries and atomization energies for small Li(n) and Li(n)H (n=1-4) clusters. In contrast to previous ab initio results, we predict a monotonic increase in atomization energies per atom with increasing cluster size for lithium clusters, in accordance with the historical Knudsen-effusion measurements of Wu. For hydrogenated lithium clusters, our results support previous theoretical work concerning the relatively low atomization energy per atom for Li(2)H compared to LiH and Li(3)H. The CCSD(T)/cc-pwCVQZ atomization energies for LiH, Li(2)H, Li(3)H, and the most stable isomer of Li(4)H, including zero-point energy corrections, are 55.7, 79.6, 113.0, and 130.6 kcal/mol, respectively. The latter results are not consistent with the most recent experiments of Wu.     

Structural and electronic properties of Si(n), Si(n)-, and PSi(n-1) clusters (2 < or = n < or = 13): Theoretical investigation based on ab initio molecular orbital theory

The geometric and electronic structures of Si(n), Si(n)-, and PSi(n-1) clusters (2 < or = n < or = 13) have been investigated using the ab initio molecular orbital theory formalism. The hybrid exchange-correlation energy functional (B3LYP) and a standard split-valence basis set with polarization functions (6-31+G(d)) were employed to optimize geometrical configurations. The total energies of the lowest energy isomers thus obtained were recalculated at the MP2/aug-cc-pVTZ level of theory. Unlike positively charged clusters, which showed similar structural behavior as that of neutral clusters [Nigam et al., J. Chem. Phys. 121, 7756 (2004)], significant geometrical changes were observed between Si(n) and Si(n)- clusters for n = 6, 8, 11, and 13. However, the geometries of P substituted silicon clusters show similar growth as that of negatively charged Si(n) clusters with small local distortions. The relative stability as a function of cluster size has been verified based on their binding energies, second difference in energy (Delta2 E), and fragmentation behavior. In general, the average binding energy of Si(n)- clusters is found to be higher than that of Si(n) clusters. For isoelectronic PSi(n-1) clusters, it is found that although for small clusters (n < 4) substitution of P atom improves the binding energy of Si(n) clusters, for larger clusters (n > or = 4) the effect is opposite. The fragmentation behavior of these clusters reveals that while small clusters prefer to evaporate monomer, the larger ones dissociate into two stable clusters of smaller size. The adiabatic electron affinities of Si(n) clusters and vertical detachment energies of Si(n)- clusters were calculated and compared with available experimental results. Finally, a good agreement between experimental and our theoretical results suggests good prediction of the lowest energy isomeric structures for all clusters calculated in the present study.     

A Gaussian-3 theoretical study of small silicon-lithium clusters: electronic structures and electron affinities of SinLi(-) (n = 2-8)

The molecular structures of neutral Si n Li ( n = 2-8) species and their anions have been studied by means of the higher level of the Gaussian-3 (G3) techniques. The lowest energy structures of these clusters have been reported. The ground-state structures of neutral clusters are "attaching structures", in which the Li atom is bound to Si n clusters. The ground-state geometries of anions, however, are "substitutional structures", which is derived from Si n+1 by replacing a Si atom with a Li (-). The electron affinities of Si n Li and Si n have been presented. The theoretical electron affinities of Si n are in good agreement with the experiment data. The reliable electron affinities of Si n Li are predicted to be 1.87 eV for Si 2Li, 2.06 eV for Si 3Li, 2.01 eV for Si 4Li, 2.61 eV for Si 5Li, 2.36 eV for Si 6Li, 2.21 eV for Si 7Li, and 3.18 eV for Si 8Li. The dissociation energies of Li atom from the lowest energy structures of Si n Li and Si atom from Si n clusters have also been estimated respectively to examine relative stabilities.     

Structure and stability of oxygen adsorption on Si(n) (n = 5-10) clusters

The structures, binding energies, and electronic properties of one oxygen atom (O) and two oxygen atoms (2O) adsorption on silicon clusters Si(n) with n ranging from 5 to 10 are studied systematically by ab initio calculations. Twelve stable structures are obtained, two of which are in agreement with those reported in previous literature and the others are new structures that have not been proposed before. Further investigations on the fragmentations of Si(n)O and Si(n)O2 (n = 5-10) clusters indicate that the pathways Si(n)O --> Si(n-1) + SiO and Si(n)O2 --> Si(n-2) + Si2O2 are most favorable from thermodynamic viewpoint. Among the studied silicon oxide clusters, Si8O, Si9O, Si5O2 and Si8O2 correspond to large adsorption energies of silicon clusters with respect to O or 2O, while Si8O, with the smallest dissociation energy, has a tendency to separate into Si7 + SiO. Using the recently developed quasi-atomic minimal-basis-orbital method, we have also calculated the unsaturated valences of the neutral Si(n) clusters. Our calculation results show that the Si atoms which have the largest unsaturated valences are more attractive to O atom. Placing O atom right around the Si atoms with the largest unsaturated valences usually leads to stable structures of the silicon oxide clusters.     

Structures and properties of neutral gallium clusters: a theoretical investigation

A systematic and unbiased structure search based on a genetic algorithm in combination with density functional theory (DFT) procedures has been carried out to locate low-energy isomers of Ga(n) up to n = 25. For the smaller clusters up to n = 8 results are checked by coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)) employing a quadruple zeta type basis set. The CCSD(T) calculations confirm a (3)Π(u) ground state for the dimer. Ga(3) has a doublet ground state 0.2 eV below two quartet states, whereas two isoenergetic triplet states are predicted for Ga(4) with D(4h) and a rhombus structure (D(2h)). Three low-lying isomers with doublet electronic states are found for Ga(5): a W-structure (C(2v)), a planar envelope (C(s)) at 0.015 eV, and a non-planar envelope (C(1)) 0.086 eV above the ground state. A triplet state for a trigonal prism (D(3h)) and a singlet for an open prism (C(2v)) are computed with virtually identical energy for Ga(6). The global minimum for Ga(7) is a capped trigonal prism (C(s)) and that for Ga(8) a distorted cube in D(2h). DFT provides a fair agreement with CCSD(T), deviations in dissociation energies are up to 0.2 eV for n ≤ 8. The structures for Ga(n) are mostly irregular for n ≥ 9, those for Ga(12) to Ga(17) can be derived from the truncated decahedron with D(5h) symmetry though highly distorted by Jahn-Teller effects, for example. For Ga(18) to Ga(23) we find stacks of five- and six-membered rings as global minima, e.g., 5-1-5-1-6 for Ga(18). Ga(24) and Ga(25) consist of layers with packing sequence ABCBA similar to those found for clusters of aluminum. The most important feature of computed cohesive energies is a rapid increase with n: for Ga(25) it reaches 2.46 eV, the experimental bulk value is 2.84 eV. Particularly stable clusters for Ga(n) are seen for n = 7, 14, and 20.     

CAl2Be3(2-) and its salt complex LiCAl2Be3-: anionic global minima with planar pentacoordinate carbon

Following the isoelectronic relationship in global minima planar pentacoordinate carbon (ppC) species (cationic CAl(5)(+), neutral CAl(4)Be, and monoanionic CAl(3)Be(2)(-)), we designed a dianionic ppC species C(2v) CAl(2)Be(3)(2-) (1a) and its salt complex C(2v) LiCAl(2)Be(3)(-) (2a) in this work. In combination with DFT and high-level ab initio calculations (CCSD(T)), the extensive exploration on their potential energy surfaces indicates that they are the global minima. Their kinetic stability was proved by two sets of 100 ps ab initio Born-Oppenheimer molecular dynamic simulations at the B3LYP/6-31+G(d) level. The detailed analyses indicate that the introduction of Li(+) into 1a only influences the electrovalent bonding (through changing of the charge distribution) and the σ aromaticity (through changing of the in-plane ring current), while the structures, the bonding properties, the π aromaticity, and so forth are almost unchanged. Nevertheless, the MO energy levels, the HOMO-LUMO gaps, and the values of vertical detachment energies (VDEs) all verify that the lithiation significantly improves the stability. We think the ppC dianion 1a is possible to detect directly in the gas-phase experiments, but it can be detected as its salt complex 2a more easily.     

SinP+m(n+m=5)结构和稳定性的理论研究

在实验的基础上 ,利用量子化学方法对 Sin P+ m( n+m=5 )的各种可能构型进行几何构型优化 ,预测各团簇的稳定结构 ,从中得出各个团簇稳定构型之间的基本关系 ,当 n>m时 ,团簇的稳定构型与 Si+ n 相似 ,而当 n相似文献