首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
The structural and electronic characteristics of the Si2C2H2 and Si2C2 clusters are studied by ab initio calculations based on coupled cluster and density functional theory using the hybrid B3LYP functional. In addition, similar species, such as SiC2H2 and Si3C2H2, are also studied for comparison. It is illustrated that the lowest energy structures of all three hydrogenated clusters, which have the general form Si(n)(CH)2, n = 1, 2, 3, are fully analogous to the structures of the corresponding organometallic isovalent carboranes. The most stable structure of Si2C2H2 is obtained by attaching two hydrogens onto the carbon atoms of a higher energy (+1.5 eV) planar trapezoidal structure of Si2C2, followed by geometry optimization which leads to puckering of the planar structure. Furthermore, it is demonstrated that Si2C2H2 and the other two "similar" hydrogenated clusters are much more stable than the corresponding bare nonhydrogenated clusters. Comparison of Si2C2H2 and C2B2H4 shows that their structural and bonding similarity includes also nuclear rearrangement similarity. The two species are isomerizable with an energy difference between their lowest energy puckered 1,2- and 1,3-isomers of about +/-0.3 eV. It is suggested that SiC2H2, Si2C2H2, and Si3C2H2 are special cases of a larger class of stable clusters. It is speculated on the basis of the calculated infrared spectrum that Si2C2H2 and perhaps other members of this class of clusters could be found in appreciable abundance in interstellar space.  相似文献   

2.
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.  相似文献   

3.
Zdetsis AD 《Inorganic chemistry》2008,47(19):8823-8829
The aromatic, bonding, and structural characteristics of the Si 4C 2H 2-C 2B 4H 6, Si 2C 4H 4-C 4B 2H 6, and other Si n C 2H 2-C 2B n H n+2 ( n = 1, 2, 3, 5) isovalent pairs are studied using density functional theory (DFT) and coupled cluster methods to fully illustrate the homology of the two species. This homology, which is based on the replacement of the carborane B-H units by isovalent Si atoms, is extended to all three characteristics (structural, electronic, and aromatic) and includes all three lowest-energy structures of the isovalent pairs. This type of "boron connection", which has been tested for silicon clusters recently, seems to be a valid and extremely useful concept. For the aromatic properties of the Si n C 2H 2-C 2B n H n+2 species, expressed through the nucleus independent chemical shifts (NICS), a strange odd-even effect with respect to the number of Si atoms is observed which seems rather difficult to explain. To help possible future identification and characterization of the Si n C 2H 2 clusters, their infrared, Raman, and optical excitation spectra are calculated within the framework of DFT, using the 6-311+G(2d, p) basis set. It is expected that the present results would facilitate the exploitation of the well-known carborane and metallacarborane chemical properties and applications for the design and development of novel silicon-carbon-based composite materials.  相似文献   

4.
The structural, electronic, vibrational, optical, magnetic, and aromatic characteristics of Si(n), Si(n) (1-), Si(n) (2-), and Si(n) (1+), clusters have been calculated very accurately with a variety of high level ab initio techniques. These calculations have been performed with the aim to clarify existing ambiguities in the literature and to bring up the fluxional and aromatic characteristics of these species. The fluxional behavior, according to earlier conjecture of the present author, could be connected to the magic property. In addition such behavior could also explain the existence of conflicting results. The ab initio techniques include quadratic configuration interaction, coupled cluster, and multireference second order perturbation theory, together with density functional theory ("static" and time dependent) with the hybrid B3LYP functional. Various high quality correlation-consistent basis sets, ranging from 2Z up to 5Z quality, were employed. It is demonstrated that Si(6) is fluxional, fluctuating around a symmetric D(4h) structure. Si(10) is also fluxional but to a lesser degree, in contrast to Si(10) (1-) anion which is highly fluxional. For both clusters, in full agreement with Wade's and Lipscomb's rules for deltahedral boranes, the corresponding dianions have higher symmetry (O(h) and D(4d), respectively) and lower energy than the neutral clusters. The aromatic behavior of Si(6) fits better to a mixed conflicting aromaticity picture. This type of aromatic and fluxional behavior has also been observed in stable "magic" carbon clusters as C(6) and carbon fullerenes such as C(20). The present results, which support possible connection of fluxional and magic properties, are in excellent agreement with experimental measurements of ionization energies, electron affinities, and vibrationally resolved photoelectron spectra.  相似文献   

5.
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.  相似文献   

6.
In a previous communication, we showed that a single Au atom behaves like H in its bonding to Si in a series of Si-Au clusters, SiAu(n) (n = 2-4) (Kiran et al. Angew. Chem., Int. Ed. 2004, 43, 2125). In this article, we show that the H analogy of Au is more general. We find that the chemical bonding and potential energy surfaces of two disilicon Au clusters, Si(2)Au(2) and Si(2)Au(4), are analogous to Si(2)H(2) and Si(2)H(4), respectively. Photoelectron spectroscopy and ab initio calculations are used to investigate the geometrical and electronic structures of Si(2)Au(2)(-), Si(2)Au(4)(-), and their neutral species. The most stable structures for both Si(2)Au(2) and Si(2)Au(2)(-) are found to be C(2)(v), in which each Au bridges the two Si atoms. For Si(2)Au(4)(-), two nearly degenerate dibridged structures in a cis (C(2)(h)) and a trans (C(2)(v)) configuration are found to be the most stable isomers. However, in the neural potential energy surface of Si(2)Au(4), a monobridged isomer is the global minimum. The ground-state structures of Si(2)Au(2)(-) and Si(2)Au(4)(-) are confirmed by comparing the computed vertical detachment energies with the experimental data. The various stable isomers found for Si(2)Au(2) and Si(2)Au(4) are similar to those known for Si(2)H(2) and Si(2)H(4), respectively. Geometrical and electronic structure comparisons with the corresponding silicon hydrides are made to further establish the isolobal analogy between a gold atom and a hydrogen atom.  相似文献   

7.
Ab initio calculations in the framework of the density functional theory, with B3LYP functional, are performed to study the lowest-energy isomers of silicon sodium clusters Si(n)Na(p)(+) (n 相似文献   

8.
Density functional theory (DFT) at the hybrid B3LYP level has been applied to the polyhedral boranes B(n)H(n)(z) (n = 8 and 11, z = -2, -4, and -6) for comparison with isoelectronic germanium clusters Ge(n)(z). The energy differences between the global minima and other higher energy borane structures are much larger relative to the case of the corresponding bare germanium clusters. Furthermore, for both B(8)H(8)(2-) and B(11)H(11)(2-), the lowest energy computed structures are the corresponding experimentally observed most spherical deltahedra predicted by the Wade-Mingos rules, namely the D(2)(d) bisdisphenoid and the C(2)(v) edge-coalesced icosahedron, respectively. Only in the case of B(8)H(8)(2-) is there a second structure close (+2.6 kcal/mol) to the D(2)(d) bisdisphenoid global minimum, namely the C(2)(v) bicapped trigonal prism corresponding to the "square" intermediate in a single diamond-square-diamond process that can lead to the experimentally observed room temperature fluxionality of B(8)H(8)(2-). Stable borane structures with 3-fold symmetry (e.g., D(3)(h), C(3)(v), etc.) are not found for boranes with 8- and 11-vertices, in contrast to the corresponding germanium clusters where stable structures derived from the D(3)(d) bicapped octahedron and D(3)(h) pentacapped trigonal prism are found for the 8- and 11-vertex systems, respectively. The lowest energy structures found for the electron-rich boranes B(8)H(8)(4-) and B(11)H(11)(4-) are nido polyhedra derived from a closo deltahedron by removal of a relatively high degree vertex, as predicted by the Wade-Mingos rules. They relate to isoelectronic species found experimentally, e.g., B(8)H(12) and R(4)C(4)B(4)H(4) for B(8)H(8)(4-) and C(2)B(9)H(11)(2-) for B(11)H(11)(4-). Three structures were found for B(11)H(11)(6-) with arachno type geometry having two open faces in accord with the Wade-Mingos rules.  相似文献   

9.
The authors present theoretical results describing the adsorption of H2 and H2S molecules on small neutral and cationic gold clusters (Au(n)((0/+1)), n=1-8) using density functional theory with the generalized gradient approximation. Lowest energy structures of the gold clusters along with their isomers are considered in the optimization process for molecular adsorption. The adsorption energies of H2S molecule on the cationic clusters are generally greater than those on the corresponding neutral clusters. These are also greater than the H2 adsorption energies on the corresponding cationic and neutral clusters. The adsorption energies for cationic clusters decrease with increasing cluster size. This fact is reflected in the elongations of the Au-S and Au-H bonds indicating weak adsorption as the cluster grows. In most cases, the geometry of the lowest energy gold cluster remains planar even after the adsorption. In addition, the adsorbed molecule gets adjusted such that its center of mass lies on the plane of the gold cluster. Study of the orbital charge density of the gold adsorbed H2S molecule reveals that conduction is possible through molecular orbitals other than the lowest unoccupied molecular orbital level. The dissociation of the cationic Au(n)SH2+ cluster into Au(n)S+ and H2 is preferred over the dissociation into Au(m)SH2+ and Au(n-m), where n=2-8 and m=1-(n-1). H2S adsorbed clusters with odd number of gold atoms are more stable than neighboring even n clusters.  相似文献   

10.
Hydroperoxide anion (HOO(-)), the conjugate base of hydrogen peroxide (HOOH), has been relatively little studied despite the importance of HOOH in commercial processes, atmospheric science, and biology. The anion has been shown to exist as a stable species in alkaline water. This project explored the structure of gas phase (HOO(-))(H(2)O)(n) clusters and identified the lowest energy configurations for n ≤ 8 at the B3LYP/6-311++G** level of theory and for n ≤ 6 at the MP2/aug-cc-pVTZ level of theory. As a start toward understanding equilibration between HOO(-) and HOOH in an alkaline environment, (HOOH)(OH(-))(H(2)O)(n-1) clusters were likewise examined, and the lowest energy configurations were determined for n ≤ 8 (B3LYP/6-311++G**) and n ≤ 6 (MP2/aug-cc-pVTZ). Some studies were also done for n = 20. The two species have very different solvation behaviors. In low energy (HOOH)(OH(-))(H(2)O)(n-1) clusters, HOOH sits on the surface of the cluster, is 4-coordinated (each O is donor once and acceptor once), and donates to the hydroxide ion. In contrast, in low energy (HOO(-))(H(2)O)(n) clusters, (HOO(-)) takes a position in the cluster center surrounded on all sides by water molecules, and its optimum coordination number appears to be 7 (one O is donor-acceptor-acceptor while the other is a 4-fold acceptor). For n ≤ 6 the lowest (HOOH)(OH(-))(H(2)O)(n-1) cluster lies 1.0-2.1 kcal/mol below the lowest (HOO(-))(H(2)O)(n) cluster, but the lowest clusters found for n = 20 favor (HOO(-))(H(2)O)(20). The results suggest that ambient water could act as a substantial kinetic brake that slows equilibration between (HOOH)(OH(-)) and (HOO(-))(H(2)O) because extensive rearrangement of solvation shells is necessary to restabilize either species after proton transfer.  相似文献   

11.
The geometries, stabilities, and electronic properties of TiSin (n=2-15) clusters with different spin configurations have been systematically investigated by using density-functional theory approach at B3LYP/LanL2DZ level. According to the optimum TiSin clusters, the equilibrium site of Ti atom gradually moves from convex to surface, and to a concave site as the number of Si atom increases from 2 to 15. When n=12, the Ti atom in TiSi12 completely falls into the center of the Si outer frame, forming metal-encapsulated Si cages, which can be explained by using 16-electron rule. On the basis of the optimized geometries, various energetic properties are calculated for the most stable isomers of TiSin clusters, including the average binding energy, the highest occupied molecular orbital and lowest unoccupied molecular orbital (HOMO-LUMO) gap, fragmentation energy, and the second-order difference of energy. It is found that at size n=6,8,12 the clusters are more stable than neighboring ones. According to the Mulliken charge population analysis, charges always transfer from Si atoms to Ti atom. Furthermore, the HOMO-LUMO gaps of the most stable TiSin clusters are usually smaller than those of Sin clusters.  相似文献   

12.
The growth mechanisms of small cationic silicon clusters containing up to 11 Si atoms, exohedrally doped by V and Cu atoms, are described. We find that as dopants, V and Cu follow two different paths: while V prefers substitution of a silicon atom in a highly coordinated position of the cationic bare silicon clusters, Cu favors adsorption to the neutral or cationic bare clusters in a lower coordination site. The different behavior of the two transition metals becomes evident in the structures of Si(n)M(+) (n = 4-11 for M = V, and n = 6-11 for M = Cu), which are investigated by density functional theory and, for several sizes, confirmed by comparison with their experimental vibrational spectra. The spectra are measured on the corresponding Si(n)M(+)·Ar complexes, which can be formed for the exohedrally doped silicon clusters. The comparison between experimental and calculated spectra indicates that the BP86 functional is suitable to predict far-infrared spectra of these clusters. In most cases, the calculated infrared spectrum of the lowest-lying isomer fits well with the experiment, even when various isomers and different electronic states are close in energy. However, in a few cases, namely Si(9)Cu(+), Si(11)Cu(+), and Si(10)V(+), the experimentally verified isomers are not the lowest in energy according to the density functional theory calculations, but their structures still follow the described growth mechanism. The different growth patterns of the two series of doped Si clusters reflect the role of the transition metal's 3d orbitals in the binding of the dopant atoms.  相似文献   

13.
We have used photoionization efficiency spectroscopy to determine ionization potentials (IP) of the niobium-carbide clusters, Nb3C(n) (n = 1-4) and Nb4C(n) (n = 1-6). The Nb3C2 and Nb4C4 clusters exhibit the lowest IPs for the two series, respectively. For clusters containing up to four carbon atoms, excellent agreement is found with relative IPs calculated using density functional theory. The lowest energy isomers are mostly consistent with the development of a 2 x 2 x 2 face-centered cubic structure of Nb4C4. However, for Nb3C4 a low-lying isomer containing a molecular C2 unit is assigned to the experimental IP rather than the depleted 2 x 2 x 2 nanocrystal isomer. For Nb4C5 and Nb4C6, interpretation is less straightforward, but results indicate isomers containing molecular C2 units are the lowest in energy, suggesting that carbon-carbon bonding is preferred when the number of carbon atoms exceeds the number of metal atoms. A double IP onset is observed for Nb4C3, which is attributed to ionization from the both the lowest energy singlet state and a meta-stable triplet state. This work further supports the notion that IPs can be used as a reliable validation for the geometries of metal-carbide clusters calculated by theory.  相似文献   

14.
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.  相似文献   

15.
C(2n)N7(-) cluster ions are produced by laser ablating on the K(3)[Fe(CN)6] sample. DFT calculations have been performed for these cluster anions. Various isomeric structures of these clusters are optimized and their energies are compared to find the most stable isomers. The most stable structure for C8N7(-) is similar to that of adenine by theoretical calculation, which is in agreement with the collision-induced dissociation (CID) experimental results. With the increasing even numbers of C atoms from 8 to 16, the N atoms in the double-ring structure are gradually substituted by C atoms from the six-membered ring to the five-membered ring. All these C(2n)N7(-) (n = 3-9) clusters exhibit planar aromatic characters. The energy difference and incremental binding energy analyses show that C(2n)N7(-) (n = 4-8) clusters are more stable than C6N7(-) and C18N7(-), which are consistent with the observed mass spectrum.  相似文献   

16.
Ab initio all-electron molecular-orbital calculations are carried out to study the structures and relative stability of low-energy silicon clusters (Si(n),n = 12-20). Selected geometric isomers include those predicted by Ho et al. [Nature (London) 392, 582 (1998)] based on an unbiased search with tight-binding/genetic algorithm, as well as those found by Rata et al. [Phys. Rev. Lett. 85, 546 (2000)] based on density-functional tight-binding/single-parent evolution algorithm. These geometric isomers are optimized at the M?ller-Plesset (MP2) MP2/6-31G(d) level. The single-point energy at the coupled-cluster single and double substitutions (including triple excitations) [CCSD(T)] CCSD(T)/6-31G(d) level for several low-lying isomers are further computed. Harmonic vibrational frequency analysis at the MP2/6-31G(d) level of theory is also undertaken to assure that the optimized geometries are stable. For Si12-Si17 and Si19 the isomer with the lowest-energy at the CCSD(T)/6-31G(d) level is the same as that predicted by Ho et al., whereas for Si18 and Si20, the same as predicted by Rata et al. However, for Si14 and Si15, the vibrational frequency analysis indicates that the isomer with the lowest CCSD(T)/6-31G(d) single-point energy gives rise to imaginary frequencies. Small structural perturbation onto the Si14 and Si15 isomers can remove the imaginary frequencies and results in new isomers with slightly lower MP2/6-31G(d) energy; however the new isomers have a higher single-point energy at the CCSD(T)/6-31G(d) level. For most Si(n) (n = 12-18,20) the low-lying isomers are prolate in shape, whereas for Si19 a spherical-like isomer is slightly lower in energy at the CCSD(T)/6-31G(d) level than low-lying prolate isomers.  相似文献   

17.
The molecular structures, electron affinities, and dissociation energies of the Si(n)H/Si(n)H- (n = 4-10) species have been examined via five hybrid and pure density functional theory (DFT) methods. The basis set used in this work is of double-zeta plus polarization quality with additional diffuse s- and p-type functions, denoted DZP++. The geometries are fully optimized with each DFT method independently. The three different types of neutral-anion energy separations presented in this work are the adiabatic electron affinity (EA(ad)), the vertical electron affinity (EA(vert)), and the vertical detachment energy (VDE). The first Si-H dissociation energies, D(e)(Si(n)H --> Si(n) + H) for neutral Si(n)H and D(e)(Si(n)H- --> Si(n)- + H) for anionic Si(n)H- species, have also been reported. The structures of the ground states of these clusters are traditional H-Si single-bond forms. The ground-state geometries of Si5H, Si6H, Si8H, and Si9H predicted by the DFT methods are different from previous calculations, such as those obtained by Car-Parrinello molecular dynamics and nonorthogonal tight-binding molecular dynamics schemes. The most reliable EA(ad) values obtained at the B3LYP level of theory are 2.59 (Si4H), 2.84 (Si5H), 2.86 (Si6H), 3.19 (Si7H), 3.14 (Si8H), 3.36 (Si9H), and 3.56 (Si10H) eV. The first dissociation energies (Si(n)H --> Si(n) + H) predicted by all of these methods are 2.20-2.29 (Si4H), 2.30-2.83 (Si5H), 2.12-2.41 (Si6H), 1.75-2.03 (Si7H), 2.41-2.72 (Si8H), 1.86-2.11 (Si9H), and 1.92-2.27 (Si10H) eV. For the negatively charged ion clusters (Si(n)H- --> Si(n)- + H), the dissociation energies predicted are 2.56-2.69 (Si4H-), 2.80-3.01 (Si5H-), 2.86-3.06 (Si6H-), 2.80-3.03 (Si7H-), 2.69-2.92 (Si8H-), 2.92-3.18 (Si9H-), and 2.89-3.25 (Si10H-) eV.  相似文献   

18.
Density functional theory (DFT) at the hybrid B3LYP level has been applied to the germanium clusters Ge(11)(z) (z = -6, -4, -2, 0, +2, +4, +6) starting from eight different initial configurations. The global minimum within the Ge(11)(2-) set is an elongated pentacapped trigonal prism distorted from D(3)(h) to C(2v) symmetry. However, the much more spherical edge-coalesced icosahedron, also of C(2v) symmetry, expected by the Wade-Mingos rules for a 2n + 2 skeletal electron system and found experimentally in B(11)H(11)(2-) and isoelectronic carboranes, is of only slightly higher energy (+5.2 kcal/mol). Even more elongated D(3)(h) pentacapped trigonal prisms are the global minima for the electron-rich structures Ge(11)(4-) and Ge(11)(6-). For Ge(11)(4-) the C(5v) 5-capped pentagonal antiprism analogous to the dicarbollide ligand C(2)B(9)H(11)(2-) is of significantly higher energy (approximately 28 kcal/mol) than the D(3h) global minimum. The C(2v) edge-coalesced icosahedron is also the global minimum for the electron-poor Ge(11) similar to its occurrence in experimentally known 11-vertex "isocloso" metallaboranes of the type (eta(6)-arene)RuB(10)H(10). The lowest energy polyhedral structures computed for the more hypoelectronic Ge(11)(4+) and Ge(11)(6+) clusters are very similar to those found experimentally for the isoelectronic ions E(11)(7-) (E = Ga, In, Tl) and Tl(9)Au(2)(9-) in intermetallics in the case of Ge(11)(4+) and Ge(11)(6+), respectively. These DFT studies predict an interesting D(5h) centered pentagonal prismatic structure for Ge(11)(2+) and isoelectronic metal clusters.  相似文献   

19.
Infrared photodissociation (IRPD) spectra of clusters composed of protonated phenol (C(6)H(7)O(+)) and several ligands L are recorded in the O-H and C-H stretch ranges using a tandem mass spectrometer coupled to a cluster ion source. The C(6)H(7)O(+)-L(n) complexes (L=Ar/N(2), n=1-6) are generated by chemical ionization of a supersonic expansion. The IRPD spectra of mass selected C(6)H(7)O(+)-L(n) clusters obtained in various C(6)H(7)O(+)-L(m) fragment channels (m相似文献   

20.
Extensive semiempirical calculations of the hexaanions of IPR (isolated pentagon rule) and non-IPR isomers of C(68)-C(88) and IPR isomers of C(90)-C(98) followed by DFT calculations of the lowest energy structures were performed to find the carbon cages that can provide the most stable isomers of M(3)N@C(2n) clusterfullerenes (M = Sc, Y) with Y as a model for rare earth ions. DFT calculations of isomers of M(3)N@C(2n) (M = Sc, Y; 2n = 68-98) based on the most stable C(2n)(6-) cages were also performed. The lowest energy isomers found by this methodology for Sc(3)N@C(68), Sc(3)N@C(78), Sc(3)N@C(80), Y(3)N@C(78), Y(3)N@C(80), Y(3)N@C(84), Y(3)N@C(86), and Y(3)N@C(88) are those that have been shown to exist by single-crystal X-ray studies as Sc(3)N@C(2n) (2n = 68, 78, 80), Dy(3)N@C(80), and Tb(3)N@C(2n) (2n = 80, 84, 86, 88) clusterfullerenes. Reassignment of the carbon cage of Sc(2)@C(76) to the non-IPR Cs: 17490 isomer is also proposed. The stability of nitride clusterfullerenes was found to correlate well with the stability of the empty 6-fold charged cages. However, the dimensions of the cage in terms of its ability to encapsulate M(3)N clusters were also found to be an important factor, especially for the medium size cages and the large Y(3)N cluster. In some cases the most stable structures are based on the different cage isomers for Sc(3)N and Y(3)N clusters. Up to the cage size of C(84), non-IPR isomers of C(2n)(6-) and M(3)N@C(2n) were found to compete with or to be even more stable than IPR isomers. However, the number of adjacent pentagon pairs in the most stable non-IPR isomers decreases as cage size increases: the most stable M(3)N@C(2n) isomers have three such pairs for 2n = 68-72, two pairs for n = 74-80, and only one pair for n = 82, 84. For C(86) and C(88) the lowest energy IPR isomers are much more stable than any non-IPR isomer. The trends in the stability of the fullerene isomers and the cluster-cage binding energies are discussed, and general rules for stability of clusterfullerenes are established. Finally, the high yield of M(3)N@C(80) (Ih) clusterfullerenes for any metal is explained by the exceptional stability of the C(80)(6-) (Ih: 31924) cage, rationalized by the optimum distribution of the pentagons leading to the minimization of the steric strain, and structural similarities of C(80) (Ih: 31924) with the lowest energy non-IPR isomers of C(760(6-), C(78)(6-), C(82)(6-), and C(84)(6-) pointed out.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号