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.     

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).     

Sodium microsolvation in ethanol: common features of Na(HO-R)n (R=H, CH3, C2H5) clusters

Ethanol clusters are generated in a continuous He seeded supersonic expansion and doped with sodium atoms in a pick-up cell. By this method clusters of the type Na(C(2)H(5)OH)(n) are formed and characterized by determining size selectively their ionization potentials (IPs) for n = 2-40 in photoionization experiments. A continuous decrease to 3.1 eV is found from n = 2 to 6 and a constant value of 3.07 ± 0.06 eV for n = 10-40. This IP evolution is similar to the sodium-water and the sodium-methanol system. Quantum chemical calculations (B3LYP and MP2) of the IPs indicate adiabatic contributions to the photoionization process for the cluster sizes n = 4 and 5, which is similar to the sodium-methanol case. The results of the extrapolated IPs and the vertical binding energies (VEBs) of cluster anions are compared with the recently reported VEBs of solvated electrons in liquid water, methanol, and ethanol solutions in the range of 3.1-3.4 eV. The new results imply that the extrapolated VBEs of solvated electrons in anionic clusters match the VBE in liquid water, while they are about 0.5 eV too low for methanol. The influence of the presence of counterions on these findings is discussed.     

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.     

Electronic structure and thermochemical properties of silicon-doped lithium clusters Li(n)Si0/+, n = 1-8: new insights on their stability

A theoretical investigation on small silicon-doped lithium clusters Li(n)Si with n = 1-8, in both neutral and cationic states is performed using the high accuracy CCSD(T)/complete basis set (CBS) method. Location of the global minima is carried out using a stochastic search method and the growth pattern of the clusters emerges as follows: (i) the species Li(n)Si with n ≤ 6 are formed by directly binding one Li to a Si of the smaller cluster Li(n-1)Si, (ii) the structures tend to have an as high as possible symmetry and to maximize the coordination number of silicon. The first three-dimensional global minimum is found for Li(4)Si, and (iii) for Li(7)Si and Li(8)Si, the global minima are formed by capping Li atoms on triangular faces of Li(6)Si (O(h)). A maximum coordination number of silicon is found to be 6 for the global minima, and structures with higher coordination of silicon exist but are less stable. Heats of formation at 0 K (Δ(f)H(0)) and 298 K (Δ(f)H(298)), average binding energies (E(b)), adiabatic (AIE) and vertical (VIE) ionization energies, dissociation energies (D(e)), and second-order difference in total energy (Δ(2)E) of the clusters in both neutral and cationic states are calculated from the CCSD(T)/CBS energies and used to evaluate the relative stability of clusters. The species Li(4)Si, Li(6)Si, and Li(5)Si(+) are the more stable systems with large HOMO-LUMO gaps, E(b), and Δ(2)E. Their enhanced stability can be rationalized using a modified phenomenological shell model, which includes the effects of additional factors such as geometrical symmetry and coordination number of the dopant. The new model is subsequently applied with consistency to other impure clusters Li(n)X with X = B, Al, C, Si, Ge, and Sn.     

Stability and structure of Li(n)H molecules (n=3-6): experimental and density functional study

The molecules Li(3)H and Li(4)H have been identified in mass-spectrometric measurements over solutions of hydrogen in liquid Li, and the gaseous equilibria of the reactions: Li(3)H+Li=Li(2)H+Li(2), Li(3)H+Li(2)=Li(2)H+Li(3), Li(3)H+Li=LiH+Li(3), Li(3)H+LiH=2Li(2)H, and Li(4)H+Li(2)=Li(3)H+Li(3) have been measured. Density functional calculations of Li(n)H molecules (n=3-6) provide structures, vibrational frequencies, ionization energies, and free energy functions of these molecules, and these are used to estimate the enthalpies of these reactions and the atomization energies of Li(3)H (119.4 kcal/mol) and Li(4)H (151.8 kcal/mol).     

Size resolved infrared spectroscopy of Na(CH3OH)n (n = 4-7) clusters in the OH stretching region: unravelling the interaction of methanol clusters with a sodium atom and the emergence of the solvated electron

Size resolved IR action spectra of neutral sodium doped methanol clusters have been measured using IR excitation modulated photoionisation mass spectroscopy. The Na(CH(3)OH)(n) clusters were generated in a supersonic He seeded expansion of methanol by subsequent Na doping in a pick-up cell. A combined analysis of IR action spectra, IP evolutions and harmonic predictions of IR spectra (using density functional theory) of the most stable structures revealed that for n = 4, 5 structures with an exterior Na atom showing high ionisation potentials (IPs) of ~4 eV dominate, while for n = 6, 7 clusters with lower IPs (~3.2 eV) featuring fully solvated Na atoms and solvated electrons emerge and dominate the IR action spectra. For n = 4 simulations of photoionisation spectra using an ab initio MD approach confirm the dominance of exterior structures and explain the previously reported appearance IP of 3.48 eV by small fractions of clusters with partly solvated Na atoms. Only for this cluster size a shift in the isomer composition with cluster temperature has been observed, which may be related to kinetic stabilisation of less Na solvated clusters at low temperatures. Features of slow fragmentation dynamics of cationic Na(+)(CH(3)OH)(6) clusters have been observed for the photoionisation near the adiabatic limit. This finding points to the relevance of previously proposed non-vertical photoionisation dynamics of this system.     

Infrared spectra and theoretical calculations of lithium hydride clusters in solid hydrogen, neon, and argon

A matrix isolation IR study of laser-ablated lithium atom reactions with H2 has been performed in solid para-hydrogen, normal hydrogen, neon, and argon. The LiH molecule and (LiH)(2,3,4) clusters were identified by IR spectra with isotopic substitution (HD, D(2), and H(2) + D(2)) and comparison to frequencies calculated by density functional theory and the MP2 method. The LiH diatomic molecule is highly polarized and associates additional H(2) to form primary (H(2))(2)LiH chemical complexes surrounded by a physical cage of solid hydrogen where the ortho and para spin states form three different primary complexes and play a role in the identification of the bis-dihydrogen complex and in characterization of the matrix cage. The highly ionic rhombic (LiH)(2) dimer, which is trapped in solid matrices, is calculated to be 22 kcal/mol more stable than the inverse hydrogen bonded linear LiH-LiH dimer, which is not observed here. The cyclic lithium hydride trimer and tetramer clusters were also observed. Although the spontaneous reaction of 2 Li and H(2) to form (LiH)(2) occurs on annealing in solid H(2), the formation of higher clusters requires visible irradiation. We observed the simplest possible chemical reduction of dihydrogen using two lithium valence electrons to form the rhombic (LiH)(2) dimer.     

Molecular structures and energetics of the (TiO2)n (n = 1-4) clusters and their anions

The (TiO2)n clusters and their anions for n = 1-4 have been studied with coupled cluster theory [CCSD(T)] and density functional theory (DFT). For n > 1, numerous conformations are located for both the neutral and anionic clusters, and their relative energies are calculated at both the DFT and CCSD(T) levels. The CCSD(T) energies are extrapolated to the complete basis set limit for the monomer and dimer and calculated up to the triple-zeta level for the trimer and tetramer. The adiabatic and vertical electron detachment energies of the anionic clusters to the ground and first excited states of the neutral clusters are calculated at both levels and compared with the experimental results. The comparison allows for the definitive assignment of the ground-state structures of the anionic clusters. Anions of the dimer and tetramer are found to have very closely lying conformations within 2 kcal/mol at the CCSD(T) level, whereas that of the trimer does not. In addition, accurate clustering energies and heats of formation are calculated for the neutral clusters and compared with the available experimental data. Estimates of the titanium-oxygen bond energies show that they are stronger than the group VIB transition metal-oxygen bonds except for tungsten. The atomization energies of these clusters display much stronger basis set dependence than the clustering energies. This allows the calculation of more accurate heats of formation for larger clusters on the basis of calculated clustering energies.     

Photoionization and ab initio study of Ba(H2O)n (n = 1-4) clusters

An experimental and theoretical study of the photoionization energies (IE's) of Ba(H(2)O)(n) clusters containing up to n = 4 water molecules has been performed. The clusters were generated by a pick-up source combining laser vaporization with pulsed supersonic expansion, and then photoionized by radiation of 272.5-340 nm. The experimentally determined IE(e)'s for n = 1 to 4 are 4.56 ± 0.05, 4.26 ± 0.05, 3.90 ± 0.05 and 3.71 ± 0.05 eV. This cluster size dependence of IE is reproduced within ±0.06 eV employing the mPW1PW91 density-functional and CCSD(T, Full) quantum-chemical methods combined with the 6-311++G(d,p) basis set for the H and O atoms and three different relativistic effective core potentials for Ba atoms. The calculations indicate that the lowest energy hydration structures represent the most relevant contributions to both the vertical and adiabatic ionization energies. Experimental and theoretical evidence correlates with the progressive surface-delocalization of the electron from the hydration cavity around the Ba atom and suggests that the intra-cluster electron transfer is possible even for small aggregates.     

Vacuum ultraviolet (VUV) photoionization of small water clusters

Tunable vacuum ultraviolet (VUV) photoionization studies of water clusters are performed using 10-14 eV synchrotron radiation and analyzed by reflectron time-of-flight (TOF) mass spectrometry. Photoionization efficiency (PIE) curves for protonated water clusters (H2O)(n)H+ are measured with 50 meV energy resolution. The appearance energies of a series of protonated water clusters are determined from the photoionization threshold for clusters composed of up to 79 molecules. These appearance energies represent an upper limit of the adiabatic ionization energy of the corresponding parent neutral water cluster in the supersonic molecular beam. The experimental results show a sharp drop in the appearance energy for the small neutral water clusters (from 12.62 +/- 0.05 to 10.94 +/- 0.06 eV, for H2O and (H2O)4, respectively), followed by a gradual decrease for clusters up to (H2O)23 converging to a value of 10.6 eV (+/-0.2 eV). The dissociation energy to remove a water molecule from the corresponding neutral water cluster is derived through thermodynamic cycles utilizing the dissociation energies of protonated water clusters reported previously in the literature. The experimental results show a gradual decrease of the dissociation energy for removal of one water molecule for small neutral water clusters (3 相似文献   

Stochastic search of the quantum conformational space of small lithium and bimetallic lithium-sodium clusters

In this paper we report the results obtained by an implementation and application of the simulated annealing optimization procedure to the exploration of the conformational space of small neutral and charged lithium clusters (Li(n)(q), n = 5, 6, 7; q = 0, +/-1) and of the bimetallic lithium/sodium clusters (Li5Na) in their lowest spin states. Our methodology eliminates the structure guessing procedure in the process of generating cluster configurations. We evaluate the quantum energy, typically with the Hartree-Fock Hamiltonian, of randomly generated points in the conformational space and use a modified Metropolis test in the annealing algorithm to generate candidate structures for atomic clusters. The structures are further optimized by analytical methods (gradient following) at the M?ller-Plesset second order perturbation theory level (MP2), in conjunction with basis sets including polarization functions with and without diffuse functions. High accuracy ab initio energies at the coupled clusters level, with single, double, and triple substitutions from the Hartree-Fock determinant (CCSD(T)), on the MP2 geometries were calculated and used to establish the relative stability of the isomers within each potential energy surface. Various cluster properties were computed and compared to existing values in order to validate our methods. Our results show excellent agreement with previous experimental and theoretical reports. Even at these small sizes, evidence for 10 new structures never reported before for the lithium clusters and four new structures for the bimetallic clusters is presented.     

A study of the reactions of molecular hydrogen with small gold clusters

This work presents a study of reactions between neutral and negatively charged Au(n) clusters (n=2,3) and molecular hydrogen. The binding energies of the first and second hydrogen molecule to the gold clusters were determined using density functional theory (DFT), second order perturbation theory (MP2) and coupled cluster (CCSD(T)) methods. It is found that molecular hydrogen easily binds to neutral Au(2) and Au(3) clusters with binding energies of 0.55 eV and 0.71 eV, respectively. The barriers to H(2) dissociation on these clusters with respect to Au(n)H(2) complexes are 1.10 eV and 0.59 eV for n=2 and 3. Although negatively charged Au(n) (-) clusters do not bind molecular hydrogen, H(2) dissociation can occur with energy barriers of 0.93 eV for Au(2) (-) and 1.39 eV for Au(3) (-). The energies of the Au(2)H(2) (-) and Au(3)H(2) (-) complexes with dissociated hydrogen molecules are lower than the energies of Au(2) (-)+H(2) and Au(3) (-)+H(2) by 0.49 eV and 0.96 eV, respectively. There is satisfactory agreement between the DFT and CCSD(T) results for binding energies, but the agreement is not as good for barrier heights.     

Structural evolution and stability of hydrogenated Li(n) (n = 1-30) clusters: a density functional study

The structure and electronic and optical properties of hydrogenated lithium clusters Li(n)H(m) (n = 1-30, m ≤ n) have been investigated by density functional theory (DFT). The structural optimizations are performed with the Becke 3 Lee-Yang-Parr (B3LYP) exchange-correlation functional with 6-311G++(d, p) basis set. The reliability of the method employed has been established by excellent agreement with computational and experimental data, wherever available. The turn over from two- to three-dimensional geometry in Li(n)H(m) clusters is found to occur at size n = 4 and m = 3. Interestingly, a rock-salt-like face-centered cubic structure is seen in Li(13)H(14). The sequential addition of hydrogen to small-sized Li clusters predicted regions of regular lattice in saturated hydrogenated clusters. This led us to focus on large-sized saturated clusters rather than to increase the number of hydrogen atoms monotonically. The lattice constants of Li(9)H(9), Li(18)H(18), Li(20)H(20), and Li(30)H(30) calculated at their optimized geometry are found to gradually approach the corresponding bulk values of 4.083. The sequential addition of hydrogen stabilizes the cluster, irrespective of the cluster size. A significant increase in stability is seen in the case of completely hydrogenated clusters, i.e., when the number of hydrogen atoms equals Li atoms. The enhanced stability has been interpreted in terms of various electronic and optical properties like adiabatic and vertical ionization potential, HOMO-LUMO gap, and polarizability.     

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.     

Small clusters of aluminum and tin: highly correlated calculations and validation of density functional procedures

We present results of molecular electronic structure treatments of multireference configuration interaction (MRCI) type for clusters Al(n) and Sn(n) in the range up to n = 4, and of coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)) type in the range up to n = 10. Basis sets of quadruple zeta size are employed, computed energy differences, such as cohesive energies, E(coh), or dissociation energies for the removal of a single atom, D(e), differ from the complete basis set limit by only a few 0.01 eV. MRCI and CCSD(T) results are then compared to those obtained from density functional theory (DFT) treatments, which show that all computational procedures agree with the general features of D(e) and E(coh). The best agreement of DFT with CCSD(T) is found for the meta-GGA (generalized gradient approximation) TPSS (Tao, Perdew, Staroverov, Scuseria) for which D(e) differs from CCSD(T) by at most 0.15 eV for Al(n) and 0.21 eV for Sn(n). The GGA PBE (Perdew, Burke, Ernzerhof) is slightly poorer with maximum deviations of 0.23 and 0.24 eV, whereas hybrid functionals are not competitive with GGA and meta-GGA functionals. A general conclusion is that errors of D(e) and/or energy differences of isomers computed with DFT procedures may easily reach 0.2 eV and errors for cohesive energies E(coh) 0.1 eV.     

Theoretical study of the spectroscopically relevant parameters for the detection of HNP(q) and HPN(q) (q = 0, +1, -1) in the gas phase

High level ab initio electronic structure calculations at different levels of theory have been performed on HNP and HPN neutrals, anions, and cations. This includes standard coupled cluster CCSD(T) level with augmented correlation-consistent basis sets, internally contacted multi-reference configuration interaction, and the newly developed CCSD(T)-F12 methods in connection with the explicitly correlated basis sets. Core-valence correction and scalar relativistic effects were examined. We present optimized equilibrium geometries, harmonic vibrational frequencies, rotational constants, adiabatic ionization energies, electron affinities, vertical detachment energies, and relative energies. In addition, the three-dimensional potential energy surfaces of HNP(-1,0,+1) and of HPN(-1,0,+1) were generated at the (R)CCSD(T)-F12b∕cc-pVTZ-F12 level. The anharmonic terms and fundamentals were derived using second order perturbation theory. For HNP, our best estimate for the adiabatic ionization energy is 7.31 eV, for the adiabatic electron affinity is 0.47 eV. The higher energy isomer, HPN, is 23.23 kcal∕mol above HNP. HPN possesses a rather large adiabatic electron affinity of 1.62 eV. The intramolecular isomerization pathways were computed. Our calculations show that HNP(-) to HPN(-) reaction is subject to electron detachment.     

Structures of small Li(NH3)n and Li(NH3)n+ clusters (n = 1-5): evidence from combined photoionization efficiency measurements and ab initio calculations

Photoionization threshold measurements have been carried out for small Li(NH3)n clusters (n = 1-5) and have been combined with ab initio calculations to determine structural information. The calculated adiabatic ionization energy for the lowest-energy isomer of each cluster is found to be in excellent agreement with the corresponding experimental photoionization threshold, providing evidence that the calculated structures are correct. The combination of the photoionization efficiency curve and the calculated adiabatic ionization energies also confirms the tentative assignment of the infrared spectrum of Li(NH3)4 reported by Salter and co-workers (J. Chem. Phys. 2006, 125, 34302); i.e., the 3 + 1 isomer does not contribute and the spectrum is due solely to the 4 + 0 isomer. The findings are consistent with an inner solvation shell that can hold a maximum of four ammonia molecules around the central lithium atom.     

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.     

Structural and electronic properties of boron-doped lithium clusters: ab initio and DFT studies

The lowest-energy structures and electronic properties of the BLi(n) (n = 1-7) clusters are reported using the B3LYP, MP2, and CCSD(T) methods with the aug-cc-pVDZ basis set. Though the results at the B3LYP level agree well with those at the CCSD(T) level, the MP2 method is rather unsatisfactory. The first three-dimensional ground state in the BLi(n) clusters occurs for BLi(4), and the impurity B atom is seen to be trapped in a Li cage from the BLi(6) cluster onwards. The evolution of the binding energies, vertical ionization potentials, and polarizability with size of cluster shows the BLi(5) cluster to be most stable among the BLi(n) clusters. Besides, the BLi(5) cluster is also found to have the largest reaction enthalpy (49.8 kcal/mol) upon losing a Li atom, which is different from the previous prediction. The unique stability of the 8-valence electron BLi(5) can be understood from the cluster electronic shell model (CSM). However, in contradiction to the prediction of the CSM, the 2s level is filled prior to the 1d level in the BLi(n) clusters.