Structural and electronic properties of TaSi(n) (n=1-13) clusters: a relativistic density functional investigation

The TaSi(n) (n=1-13) clusters with doublet, quartet, and sextet spin configurations have been systematically investigated by a relativistic density functional theory with the generalized gradient approximation available in Amsterdam density functional program. The total bonding energies, equilibrium geometries, Mulliken populations as well as Hirshfeld charges of TaSi(n) (n=1-13) clusters are calculated and presented. The emphasis on the stabilities and electronic properties is discussed. The most stable structures of the small TaSi(n) (n=1-6) clusters and the evolutional rule of low-lying geometries of the larger TaSi(n) (n=7-13) clusters are obtained. Theoretical results indicate that the most stable structure of TaSi(n) (n=1-6) clusters keeps the similar framework as the most stable structure of Si(n+1) clusters except for TaSi(3) cluster. The Ta atom in the lowest-energy TaSi(n) (n=1-13) isomers occupies a gradual sinking site, and the site moves from convex, to flatness, and to concave with the number of Si atom varying from 1 to 13. When n=12, the Ta atom in TaSi(12) cluster completely falls into the center of the Si frame, and a cagelike TaSi(12) geometry is formed. Meanwhile, the net Mulliken and Hirsheld populations of the Ta atom in the TaSi(n) (n=1-13) clusters vary from positive to negative, manifesting that the charges in TaSi(n) (n>/=12) clusters transfer from Si atoms to Ta atom. Additionally, the contribution of Si-Si and Si-Ta interactions to the stability of TaSi(n) clusters is briefly discussed. Furthermore, the investigations on atomic averaged binding energies and fragmentation energies show that the TaSi(n) (n=2,3,5,7,10,11,12) clusters have enhanced stabilities. Compared with pure silicon clusters, a universal narrowing of highest occupied molecular orbital-lowest unoccupied molecular orbital gap in TaSi(n) clusters is found.     

Geometries, stabilities, and electronic properties of different-sized ZrSi(n) (n=1-16) clusters: a density-functional investigation

The ZrSi(n) (n=1-16) clusters with different spin configurations have been systematically investigated by using the density-functional approach. The total energies, equilibrium geometries, growth-pattern mechanisms, natural population analysis, etc., are discussed. The equilibrium structures of different-sized ZrSi(n) clusters can be determined by two evolution patterns. Theoretical results indicate that the most stable ZrSi(n) (n=1-7) geometries, except ZrSi3, keep the analogous frameworks as the lowest-energy or the second lowest-energy Si(n+1) clusters. However, for large ZrSi(n) (n=8-16) clusters, Zr atom obviously disturbs the framework of silicon clusters, and the localized position of the transition-metal (TM) Zr atom gradually varies from the surface insertion site to the concave site of the open silicon cage and to the encapsulated site of the sealed silicon cage. It should be mentioned that the lowest-energy sandwich-like ZrSi12 geometry is not a sealed structure and appears irregular as compared with other TM@Si12 (TM = Re,Ni). The growth patterns of ZrSi(n) (n=1-16) clusters are concerned showing the Zr-encapsulated structures as the favorable geometries. In addition, the calculated fragmentation energies of the ZrSi(n) (n=1-16) clusters manifest that the magic numbers of stabilities are 6, 8, 10, 14, and 16, and that the fullerene-like ZrSi16 is the most stable structure, which is in good agreement with the calculated atomic binding energies of ZrSi(n) (n=8-16) and with available experimental and theoretical results. Natural population analysis shows that the natural charge population of Zr atom in the most stable ZrSi(n) (n=1-16) structures exactly varies from positive to negative at the critical-sized ZrSi8 cluster; furthermore, the charge distribution around the Zr atom appears clearly covalent in character for the small- or middle-sized clusters and metallic in character for the large-sized clusters. Finally, the properties of frontier orbitals and polarizabilities of ZrSi(n) are also discussed.     

The structural and electronic properties of In(n)N(n = 1-13) clusters

The structural and electronic properties of In(n)N(n=1-13) clusters have been investigated by density-functional theory with the generalized gradient approximation. The results indicate that the equilibrium structures of In(n)N are linear for n=1,2, planar for n=3-5, and three dimensional for n=6-13. Maximum peaks were observed for In(n)N clusters at n=3,7,9 on the size dependence for second-order energy difference. These imply that these clusters possess relatively higher stability, which is consistent with the case of binding energy per atom. Moreover, the results show that the bonding in small In(n)N clusters has a little ionic character by Mulliken population analysis. The energy gap between the highest occupied and lowest unoccupied molecular orbitals, the vertical ionization potential and electron vertical affinity (VIP and VEA) form an even-odd alternating pattern with increasing cluster size. In general, the VIP tends to lower as the cluster size increases, while the VEA tends to increase as the cluster size increases.     

Structures and electronic properties of small FeB n (n=1-10) clusters

The geometries, stabilities, electronic properties, and magnetism of FeB(n) clusters up to n=10 are systematically studied with density functional theory. We find that our optimized structures of FeB(2), FeB(3), FeB(4), and FeB(5) clusters are more stable than those proposed in previous literature. The results show that it is favorable for the Fe atom to locate at the surface, not at the center of the cluster, and that FeB(4) and FeB(9) clusters exhibit high stability. For all the FeB(n) clusters studied, we find the charge transfer from Fe to B site and the coexistence of ionic and covalent bonding characteristics. The computed total magnetic moments of the lowest-energy structures oscillate with the cluster size and are quenched at n=4, 6, 8, and 10.     

Equilibrium geometries, stabilities, and electronic properties of the bimetallic M2-doped Au(n) (M = Ag, Cu; n = 1-10) clusters: comparison with pure gold clusters

The density functional method with relativistic effective core potential has been employed to investigate systematically the geometrical structures, relative stabilities, growth-pattern behaviors, and electronic properties of small bimetallic M(2)Au(n) (M = Ag, Cu; n = 1-10) and pure gold Au(n) (n ≤ 12) clusters. The optimized geometries reveal that M(2) substituted Au(n+2) clusters and one Au atom capped M(2)Au(n-1) structures are dominant growth patterns of the stable alloyed M(2)Au(n) clusters. The calculated averaged atomic binding energies, fragmentation energies, and the second-order difference of energies as a function of the cluster size exhibit a pronounced even-odd alternation phenomenon. The analytic results exhibit that the planar structure Ag(2)Au(4) and Cu(2)Au(2) isomers are the most stable geometries of Ag(2)Au(n) and Cu(2)Au(n) clusters, respectively. In addition, the HOMO-LUMO gaps, charge transfers, chemical hardnesses and polarizabilities have been analyzed and compared further.     

Optical properties of (GaAs)n clusters (n = 2-16)

The electronic and geometrical structures of the lowest triplet states of (GaAs) n clusters ( n = 2-16) are studied using density functional theory with generalized gradient approximation (DFT-GGA). It is found that the triplet-state geometries are different from the corresponding singlet-state geometries; for n = 2-8, 10, and 11, the triplets and singlets have different topologies, while the (GaAs) 9, (GaAs) 12, (GaAs) 15, and (GaAs) 16 triplets possess a reduced symmetry, due to Jahn-Teller distortions. Except for GaAs, the singlet states are the ground states. Excitation energies and oscillator strengths are computed for excitations from the ground state to ten singlet states of all (GaAs) n clusters using time-dependent density functional theory. The adiabatic singlet-triplet gap is compared to the vertical gap, and the difference in the eigenvalues of the highest-occupied and lowest-unoccupied molecular orbitals (the HOMO-LUMO gap). While these three values show large oscillations for small n, they approach each other as the cluster size grows. Thus, the HOMO-LUMO gap computed using the DFT-GGA approach presents a rather reliable estimate of the adiabatic singlet-triplet gap.     

Structure and properties of Mn(n), Mn(n)-, and Mn(n)+ clusters (n = 3-10)

Electronic and geometrical structures of Mn(3)-Mn(10) together with their singly negatively and positively charged ions are computed using density functional theory with generalized gradient approximation. The ground-state spin multiplicities in the neutral series are 16, 21, 4, 9, 6, 5, 2, and 5, for Mn(3)-Mn(10), respectively. Thus, there is a transition from a ferromagnetic ground state to a ferrimagnetic ground state at Mn(5). The energy difference between ferrimagnetic and ferromagnetic states in Mn(n) grows rapidly with increasing n and exceeds 2 eV in Mn(10). The corresponding change from ferro- to ferrimagnetic ground state occurs at Mn(6)(-) and Mn(3)(+) in the anionic and cationic series, respectively. Beginning with Mn(6), the ion spin multiplicities differ from that of the neutral by +/-1 (i.e., they obey the empirical "+/-1 rule"). We found that the energy required to remove an Mn atom is nearly independent of the charge state of an Mn(n) cluster and the number of atoms in the cluster, except for Mn(3). The results of our calculations are in reasonable agreement with experiment, except for the experimental data on the magnetic moments per atom, where, in general, we predict smaller values than the experiment.     

Geometries and electronic properties of the tungsten-doped germanium clusters: WGen (n = 1-17)

Geometries associated with relative stabilities, energy gaps, and polarities of W-doped germanium clusters have been investigated systematically by using density functional theory. The threshold size for the endohedral coordination and the critical size of W-encapsulated Gen structures emerge as, respectively, n = 8 and n = 12, while the fullerene-like W@Ge(n) clusters appears at n = 14. The evaluated relative stabilities in term of the calculated fragmentation energies reveal that the fullerene-like W@Ge(14) and W@Ge(16) structures as well as the hexagonal prism WGe(12) have enhanced stabilities over their neighboring clusters. Furthermore, the calculated polarities of the W@Ge(n) reveal that the bicapped tetragonal antiprism WGe(10) is a polar molecule while the hexagonal prism WGe(12) is a nonpolar molecule. Moreover, the recorded natural populations show that the charges transfer from the germanium framework to the W atom. Additionally, the WGe(12) cluster with large highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) gap, large fragmentation energy, and large binding energy is supposed to be suitable as a building block of assembly cluster material. It should be pointed out that the remarkable features of W@Ge(n) clusters above are distinctly different from those of transition metal (TM) doped Ge(n) (TM = Cu and Ni) clusters, indicating that the growth pattern of the TMGe(n) depends on the kind of doped TM impurity.     

The structural and electronic properties of Ag-adsorbed (SiO2)n (n=1-7) clusters

Equilibrium geometries, charge distributions, stabilities, and electronic properties of the Ag-adsorbed (SiO(2))(n) (n=1-7) clusters have been investigated using density functional theory with generalized gradient approximation for exchange-correlation functional. The results show that the Ag atom preferably binds to silicon atom with dangling bond in nearly a fixed direction, and the incoming Ag atoms tend to cluster on the existing Ag cluster leading to the formation of Ag islands. The adsorbed Ag atom only causes charge redistributions of the atoms near itself. The effect of the adsorbed Ag atom on the bonding natures and structural features of the silica clusters is minor, attributing to the tendency of stability order of Ag(SiO(2))(n) (n=1-7) clusters in consistent with silica clusters. In addition, the energy gaps between the highest occupied and lowest unoccupied molecular orbitals remarkably decrease compared with the pure (SiO(2))(n) (n=1-7) clusters, eventually approaching the near infrared radiation region. This suggests that these small clusters may be an alternative material which has a similar functionality in treating cancer to the large gold-coated silica nanoshells and the small Au(3)(SiO(2))(3) cluster.     

Structural and electronic properties of Si n, Si n+, and AlSi n-1 (n=2-13) clusters: theoretical investigation based on ab initio molecular orbital theory

The geometric and electronic structures of Si(n), Si(n) (+), and AlSi(n-1) clusters (2< or =n< or =13) have been investigated using the ab initio molecular orbital theory under the density functional theory formalism. The hybrid exchange-correlation energy function (B3LYP) and a standard split-valence basis set with polarization functions [6-31G(d)] were employed for this purpose. Relative stabilities of these clusters have been analyzed based on their binding energies, second difference in energy (Delta (2)E) and fragmentation behavior. The equilibrium geometry of the neutral and charged Si(n) clusters show similar structural growth. However, significant differences have been observed in the electronic structure leading to their different stability pattern. While for neutral clusters, the Si(10) is magic, the extra stability of the Si(11) (+) cluster over the Si(10) (+) and Si(12) (+) bears evidence for the magic behavior of the Si(11) (+) cluster, which is in excellent agreement with the recent experimental observations. Similarly for AlSi(n-1) clusters, which is isoelectronic with Si(n) (+) clusters show extra stability of the AlSi(10) cluster suggesting the influence of the electronic structures for different stabilities between neutral and charged clusters. The ground state geometries of the AlSi(n-1) clusters show that the impurity Al atom prefers to substitute for the Si atom, that has the highest coordination number in the host Si(n) cluster. The fragmentation behavior of all these clusters show that while small clusters prefers to evaporate monomer, the larger ones dissociate into two stable clusters of smaller size.     

Density-functional study of structural and electronic properties of Si(n)C(n) (n=1-10) clusters

Density-functional theory with generalized gradient approximation for the exchange-correlation potential has been used to calculate the structural and electronic structure of Si(n)C(n) (n=1-10) clusters. The geometries are found to undergo a structural change from two dimensional to three dimensional when the cluster size n equals 4. Cagelike structures are favored as the cluster size increases. A distinct segregation between the silicon and carbon atoms is observed for these clusters. It is found that the C atoms favor to form five-membered rings as the cluster size n increases. However, the growth motif for Si atoms is not observed. The Si(n)C(n) clusters at n=2, 6, and 9 are found to possess relatively higher stability. On the basis of the lowest-energy geometries obtained, the size dependence of cluster properties such as binding energy, HOMO-LUMO gap, Mulliken charge, vibrational spectrum, and ionization potential has been computed and analyzed. The bonding characteristics of the clusters are discussed.     

Density functional study of the structures and electronic properties of nitrogen-doped Ni(n) clusters, n = 1-10

We present a first-principles study of the equilibrium geometries, electronic structure, and related properties (binding energies, ionization potentials, electron affinities, and magnetic moments) of free-standing Ni(n) (n = 1-10) clusters doped with one impurity of N. Calculations have been performed in the framework of the density functional theory, as implemented in the SIESTA code within the generalized gradient approximation to exchange and correlation. We show that, in contrast to the molecular adsorption of N(2), the adsorption of a single N atom can dramatically change the structure of the host Ni(n) cluster, examples of which are Ni(5)N, Ni(7)N, and Ni(10)N, and that noticeable structure relaxations take place otherwise. Doping with a nitrogen impurity increases the binding energy as well as the ionization potential (except for Ni(6)N), which proves that N-doping works in favor of stabilizing the Ni clusters. We also find that the magnetic moments decrease in most cases upon N-doping despite the fact that the average Ni-Ni distance increases. The HUMO-LUMO gap for one spin channel strongly changes as a function of size upon N-doping, in contrast with the HUMO-LUMO gap for the other spin channel. This might have important implication in electronic transport properties through these molecular contacts anchored to source and drain electrodes.     

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.     

Density functional calculation of the structure and electronic properties of Cu(n)O(n) (n = 1-8) clusters

Ab initio simulations and calculations were used to study the structures and stabilities of copper oxide clusters, Cu(n)O(n) (n = 1-8). The lowest energy structures of neutral and charged copper oxide clusters were determined using primarily the B3LYP/LANL2DZ model chemistry. For n ≥ 4, the clusters are nonplanar. Selected electronic properties including atomization energies, ionization energies, electron affinities, and Bader charges were calculated and examined as a function of n.     

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.     

Geometries and electronic properties of the neutral and charged rare earth Yb-doped Si(n) (n = 1-6) clusters: a relativistic density functional investigation

The neutral and charged YbSi(n) (n = 1-6) clusters considering different spin configurations have been systematically investigated by using the relativistic density functional theory with generalized gradient approximation. The total bonding energies, equilibrium geometries, Mulliken populations (MP), Hirshfeld charges (HC), fragmentation energies, and highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) gaps are calculated and discussed. The optimized geometries indicate that the most stable YbSi(n) (n = 1-6) clusters keep basically the analogous frameworks as the low-lying Si(n)(+1) clusters, while the charged species deviate from their neutral counterparts, and that the doped Yb tends to occupy the substitutional site of the neutral and charged YbSi(n) isomers. The relative stabilities are investigated in terms of the calculated fragmentation energies, exhibiting enhanced stabilities for the remarkably stable neutral and charged YbSi2 and YbSi5 clusters. Furthermore, the calculated MP and HC values show that the charges of the neutral and charged YbSi(n) clusters transfer from the Yb atom to Si(n) atoms and the Yb atom acts as an electron donor, and that the f orbitals of the Yb atom in the neutral and charged YbSi(n) clusters behave as core without involvement in chemical bonding. The calculated HOMO-LUMO gaps indicate that the YbSi2 and YbSi4+ clusters have stronger chemical stabilities. Comparisons of the Yb-doped Si(n) (n = 1-6) with available theoretical results of transition-metal-doped silicon clusters are made. The growth pattern is investigated also.     

Low-energy isomer identification, structural evolution, and magnetic properties in manganese-doped gold clusters MnAu(n) (n = 1-16)

The size-dependent electronic, structural, and magnetic properties of Mn-doped gold clusters have been systematically investigated by using relativistic all-electron density functional theory with generalized gradient approximation. A number of new isomers are obtained for neutral MnAu(n) (n = 1-16) clusters to probe the structural evolution. The two-dimensional (2D) to three-dimensional (3D) transition occurs in the size range n = 7-10 with manifest structure competitions. From size n = 13 to n = 16, the MnAu(n) prefers a gold cage structure with Mn atom locating at the center. The relative stabilities of the ground-state MnAu(n) clusters show a pronounced odd-even oscillation with the number of Au atoms. The magnetic moments of MnAu(n) clusters vary from 3 μ(B) to 6 μ(B) with the different cluster size, suggesting that nonmagnetic Au(n) clusters can serve as a flexible host to tailor the dopant's magnetism, which has potential applications in new nanomaterials with tunable magnetic properties.     

Structure and stability of Al13H(n) (n=1-13) clusters: exceptional stability of Al13H13

We have performed density functional calculations for the structure and stability of Al(13)H(n) (n=1-13) clusters. Population analysis has shown significant charge transfer occurring from the Al cluster to the H atoms. The population for Al(13) varies from 0.24 (Al(13)H) to 2.83 (Al(13)H(13)). The shape of Al(13) moieties in the Al(13)H(n) (n>/=8) clusters is significantly distorted from the icosahedral structure of Al(13) and is a "cagelike" form. Al(13)H(13) has a capped icosahedron as the ground-state structure, similar to B(13)H(13), while the shape of B(13) (planar) is different from Al(13) (icosahedral). The Al(13)H(13) cluster is predicted to be exceptionally stable on the basis of the high stabilization energy and the negative nucleus independent chemical shift value.     

Trends in structural, electronic and energetic properties of bimetallic vanadium-gold clusters Au(n)V with n = 1-14

A systematic quantum chemical investigation on the electronic, geometric and energetic properties of Au(n)V clusters with n = 1-14 in both neutral and anionic states is performed using BP86/cc-pVTZ-PP calculations. Most clusters having an even number of electrons prefer a high spin state. For odd-electron systems, a quartet state is consistently favoured as the ground state up to Au(8)V. The larger sized Au(10)V, Au(12)V and Au(14)V prefer a doublet state. The clusters prefer 2D geometries up to Au(8)V involving a weak charge transfer. The larger systems bear 3D conformations with a more effective electron transfer from Au to V. The lowest-energy structure of a size Au(n)V is built upon the most stable form of Au(n-1)V. During the growth, V is endohedrally doped in order to maximize its coordination numbers and augment the charge transfer. Energetic properties, including the binding energies, embedding energies and second-order energy differences, show that the presence of a V atom enhances considerably the thermodynamic stability of odd-numbered gold clusters but reduces that of even-numbered systems. The atomic shape has an apparently more important effect on the clusters stability than the electronic structure. Especially, if both atomic shape and electronic condition are satisfied, the resulting cluster becomes particularly stable such as the anion Au(12)V(-), which can thus combine with the cation Au(+) to form a superatomic molecule of the type [Au(12)V]Au. Numerous lower-lying electronic states of these clusters are very close in energy, in such a way that DFT computations cannot clearly establish their ground electronic states. Calculated results demonstrate the existence of structural isomers with comparable energy content for several species including Au(9)V, Au(10)V, Au(13)V and Au(14)V.     

Size- and charge-dependent geometric and electronic structures of Bi n (Bi(-)n) clusters (n=2-13) by first-principles simulations

Neutral and negatively charged bismuth clusters, Bi n and Bi(-)n (n=2-13), are investigated by first-principles simulations with the scalar-relativistic projector-augmented wave potential and the spin-polarized generalized gradient approximation. Both types of clusters show size-dependent odd-even oscillations in stability, density of states, and vertical and adiabatic electron affinities, in close agreement with experiment. The negative charge thoroughly reverses the oscillations and considerably influences the geometric structures, particularly of the odd-sized clusters. We note that most atoms in the ground states and the low-lying isomers are three coordinated with a quasilayerlike growth mode based on pentagon units, due to a weak s-p hybridization. The Bi12 cluster is found to prefer a small elongated tubelike structure with the surface consists of six curved-pentagon rings and two triangular facets, which may be the basis for the formation of bismuth nanotubes experimentally reported.