Optical absorption in boron clusters B6_{6} and B6+_{6}^{+} : A first principles configuration interaction approach

The linear optical absorption spectra in neutral boron cluster B$_{6}$ and cationic B$_{6}^{+}$ are calculated using a first principles correlated electron approach. The geometries of several low-lying isomers of these clusters were optimized at the coupled-cluster singles doubles (CCSD) level of theory. With these optimized ground-state geometries, excited states of different isomers were computed using the singles configuration-interaction (SCI) approach. The many body wavefunctions of various excited states have been analysed and the nature of optical excitation involved are found to be of collective, plasmonic type.Comment: 22 pages, 38 figures. An invited article submitted to European Physical Journal D. This work was presented in the International Symposium on Small Particles and Inorganic Clusters - XVI, held in Leuven, Belgiu

Ionic and electronic structure of sodium clusters up to N=59

We determined the ionic and electronic structure of sodium clusters with even electron numbers and 2 to 59 atoms in axially averaged and three-dimensional density functional calculations. A local, phenomenological pseudopotential that reproduces important bulk and atomic properties and facilitates structure calculations has been developed. Photoabsorption spectra have been calculated for $\mathrm{Na}_2$, $\mathrm{Na}_8$, and $\mathrm{Na}_9^+$ to $\mathrm{Na}_{59}^+$. The consistent inclusion of ionic structure considerably improves agreement with experiment. An icosahedral growth pattern is observed for $\mathrm{Na}_{19}^+$ to $\mathrm{Na}_{59}^+$. This finding is supported by photoabsorption data.Comment: To appear in Phys. Rev. B 62. Version with figures in better quality can be requested from the author

Ionic structure and photoabsorption in medium sized sodium clusters

We present ground-state configurations and photoabsorption spectra of Na-7+, Na-27+ and Na-41+. Both the ionic structure and the photoabsorption spectra of medium-size sodium clusters beyond Na-20 have been calculated self-consistently with a nonspherical treatment of the valence electrons in density functional theory. We use a local pseudopotential that has been adjusted to experimental bulk properties and the atomic 3s level of sodium. Our studies have shown that both the ionic structure of the ground state and the positions of the plasmon resonances depend sensitively on the pseudopotential used in the calculation, which stresses the importance of its consistent use in both steps.Comment: 4 pages, 3 figures. Accepted for publication in PRB, tentatively July 15th, 1998 some typos corrected, brought to nicer forma

Electronic-structure-induced deformations of liquid metal clusters

Ab initio molecular dynamics is used to study deformations of sodium clusters at temperatures $500\cdots 1100$ K. Open-shell Na$_{14}$ cluster has two shape isomers, prolate and oblate, in the liquid state. The deformation is stabilized by opening a gap at the Fermi level. The closed-shell Na$_8$ remains magic also at the liquid state.Comment: REVTex, 11 pages, no figures, figures (2) available upon request (e-mail to hakkinen at jyfl.jyu.fi), submitted to Phys. Rev.

Thermal expansion in small metal clusters and its impact on the electric polarizability

The thermal expansion coefficients of $\mathrm{Na}_{N}$ clusters with $8 \le N \le 40$ and $\mathrm{Al}_{7}$, $\mathrm{Al}_{13}^-$ and $\mathrm{Al}_{14}^-$ are obtained from {\it ab initio} Born-Oppenheimer LDA molecular dynamics. Thermal expansion of small metal clusters is considerably larger than that in the bulk and size-dependent. We demonstrate that the average static electric dipole polarizability of Na clusters depends linearly on the mean interatomic distance and only to a minor extent on the detailed ionic configuration when the overall shape of the electron density is enforced by electronic shell effects. The polarizability is thus a sensitive indicator for thermal expansion. We show that taking this effect into account brings theoretical and experimental polarizabilities into quantitative agreement.Comment: 4 pages, 2 figures, one table. Accepted for publication in Physical Review Letters. References 10 and 23 update

Optical response of small silver clusters

The time-dependent local density approximation is applied to the optical response of the silver clusters, Ag_2, Ag_3, Ag_8 and Ag_9^+. The calculation includes all the electrons beyond the closed-shell Ag^{+11} ionic core, thus including for the first time explicitly the filled d-shell in the response. The excitation energy of the strong surface plasmon near 4 eV agrees well with experiment. The theoretical transition strength is quenched by a factor of 4 with respect to the pure s-electron sum rule in Ag_8 due to the d-electrons. A comparable amount of strength lies in complex states below 6 eV excitation. The total below 6 eV, about 50% of the s sum rule, is consistent with published experiments.Comment: 13 pages RevTex and 9 Postscript figure

Isomer specific spectroscopy of metal clusters trapped in a matrix: Ag9

Clusters of metal atoms at a fixed size can assume different structural arrangements, known as isomers, which may have nearly the same energy. Therefore at given experimental conditions distribution of these isomers can be present. While the size selection is a relatively common technique, the isomer selection is not; it is therefore more difficult to obtain information about a single isomer. We report here on isomer specific spectroscopy of Ag9 clusters together with ab initio calculations allowing to identify the isomer responsible for the measured excitation pattern and fluorescence

Optical absorption spectra of finite systems from a conserving Bethe-Salpeter equation approach

We present a method for computing optical absorption spectra by means of a Bethe-Salpeter equation approach, which is based on a conserving linear response calculation for electron-hole coherences in the presence of an external electromagnetic field. This procedure allows, in principle, for the determination of the electron-hole correlation function self-consistently with the corresponding single-particle Green function. We analyze the general approach for a "one-shot" calculation of the photoabsorption cross section of finite systems, and discuss the importance of scattering and dephasing contributions in this approach. We apply the method to the closed-shell clusters Na_4, Na^+_9 and Na^+_(21), treating one active electron per Na atom.Comment: 9 pages, 3 figure

An orbital-free molecular dynamics study of melting in K_20, K_55, K_92, K_142, Rb_55 and Cs_55 clusters

The melting-like transition in potasium clusters K_N, with N=20, 55, 92 and 142, is studied by using an orbital-free density-functional constant-energy molecular dynamics simulation method, and compared to previous theoretical results on the melting-like transition in sodium clusters of the same sizes. Melting in potasium and sodium clusters proceeds in a similar way: a surface melting stage develops upon heating before the homogeneous melting temperature is reached. Premelting effects are nevertheless more important and more easily established in potasium clusters, and the transition regions spread over temperature intervals which are wider than in the case of sodium. For all the sizes considered, the percentage melting temperature reduction when passing from Na to K clusters is substantially larger than in the bulk. Once those two materials have been compared for a number of different cluster sizes, we study the melting-like transition in Rb_55 and Cs_55 clusters and make a comparison with the melting behavior of Na_55 and K_55. As the atomic number increases, the height of the specific heat peaks decreases, their width increases, and the melting temperature decreases as in bulk melting, but in a more pronounced way.Comment: LaTeX file. 6 pages with 17 pictures. Final version with minor change

Evolution of electronic and ionic structure of Mg-clusters with the growth cluster size

The optimized structure and electronic properties of neutral and singly charged magnesium clusters have been investigated using ab initio theoretical methods based on density-functional theory and systematic post-Hartree-Fock many-body perturbation theory accounting for all electrons in the system. We have systematically calculated the optimized geometries of neutral and singly charged magnesium clusters consisting of up to 21 atoms, electronic shell closures, binding energies per atom, ionization potentials and the gap between the highest occupied and the lowest unoccupied molecular orbitals. We have investigated the transition to the hcp structure and metallic evolution of the magnesium clusters, as well as the stability of linear chains and rings of magnesium atoms. The results obtained are compared with the available experimental data and the results of other theoretical works.Comment: 30 pages, 10 figures, 3 table