Plasmon response in K, Na and Li clusters: systematics using the separable random-phase-approximation with pseudo-Hamiltonians

The systematics of the plasmon response in spherical K, Na and Li clusters in a wide size region is studied. Two simplifying approximations whose validity has been established previously are considered: (a) a separable approach to the random-phase-approximation, involving an expansion of the residual interaction into a sum of separable terms, (b) the electron-ion interaction is modeled within the pseudo-Hamiltonian jellium model (PHJM) including nonlocal effects by means of realistic atomic pseudoHamiltonians. In cases where nonlocal effects turn out to be negligible, the Structure Averaged Jellium Model (SAJM) has been used. The leading role of Landau damping in forming the plasmon width in medium and large clusters is demonstrated. Good agreement with available experimental data is achieved for K, Na (using the SAJM) and small Li clusters (invoking the PHJM). The trends for peak position and width are generally well reproduced, even up to details of the Landau fragmentation in several clusters. Less good agreement, however, is found for large Li clusters. The possible reasons of the discrepancy are discussed. Received: 22 April 1998 / Accepted: 24 July 1998     

Nanoplasma dynamics in Xe clusters driven by ultraintense laser fields

We present a theoretical and computational study of the properties and the response of the nanoplasma and of outer ionization in Xe_n clusters (n = 55–2171, initial cluster radius R₀ = 8.7–31.0 ?) driven by ultraintense near-infrared laser fields (peak intensity I_M = 10¹⁵–10²⁰ Wcm^-2, temporal pulse length τ= 10–100 fs, and frequency ν= 0.35 fs^-1). The positively charged high-energy nanoplasma produced by inner ionization nearly follows the oscillations of the fs laser pulse and can either be persistent (at lower intensities of I_M = 10¹⁵–10¹⁶ Wcm^-2 and/or for larger cluster sizes, where the electron energy distribution is nearly thermal) or transient (at higher intensities of I_M = 10¹⁸–10²⁰ Wcm^-2 and/or for smaller cluster sizes). The nanoplasma is depleted by outer ionization that was semiquantitatively described by the cluster barrier suppression electrostatic model, which accounts for the cluster size, laser intensity and pulse length dependence of the outer ionization yield. The electrostatic model was further utilized for estimates of the laser intensity and pulse width dependence of the border radius R₀ ^(I) for the attainment of complete outer ionization at , while at R₀ > R₀ ^(I) a persistent nanoplasma prevails. R₀ ^(I) establishes an interrelationship between electron dynamics and nuclear Coulomb explosion dynamics in ultraintense laser-cluster interactions.     

Towards spectral pattern of spin polarized sodium clusters: The example of Na

We investigate spin modes in the ground state and the polarized first isomer of the Na₁₂ cluster describing the valence electrons in time-dependent local-spin-density approximation (TDLSDA) and the detailed ionic background using local pseudopotentials. The spin modes show a collective redshift compared to the unperturbed particle-hole excitations. They are strongly fragmented and the average energy of the modes along the principal axes are related to the underlying geometry (triaxial or axially symmetric). For the polarized isomer, we find significant cross talk between the spin modes and the dipole plasmon, which hints at a possible spectroscopic identification. Received: 22 June 1998 / Accepted: 29 July 1998     

Electronic structure and magnetism in 4d-transition metal clusters

We analyze the stability of magnetic states obtained within the tight-binding model for cubooctahedral (O_h) and icosahedral (I_h) clusters of early 4d (Y, Zr, Nb, Mo, and Tc) transition metals. Several metastable magnetic clusters are identified which suggests the existence of multiple magnetic solutions in realistic systems. A bulk-like parabolic behavior is observed for the binding energy of O_h and I_h clusters as a function of the atomic number along the 4 d-series. The charge transfer on the central atom changes sign, while the average magnetic moments present an oscillatory behavior as a function of the number of d electrons in the cluster. Our results are in agreement with other theoretical calculations. Received: 20 November 1997 / Received in final form: 9 March 1998 / Accepted: 30 March 1998     

Trianionic gold clusters

Using Penning-trap experiments and a shell-correction method incorporating ellipsoidal shape deformations, we investigate the formation and stability patterns of trianionic gold clusters. Theory and experiment are in remarkable agreement concerning appearance sizes and electronic shell effects. In contrast to multiply cationic clusters, decay of the trianionic gold clusters occurs primarily via electron autodetachment and tunneling through a Coulomb barrier, rather than via fission. Received 9 January 2001     

Charge transfer between alkali cluster ions and atoms in the 1 to 10 keV collisional energy range

The cross-sections for collisional charge transfer between singly charged free clusters M _n ⁺ (M = Li, Na; n=1...50) and atomic targets A (cesium, potassium) have been measured as a function of collisional relative velocity in laboratory energy range 1–10 keV. For each cluster size, the experimental values of the charge transfer cross-section are fitted with an universal parametric curve with two independent parameters and v_m, the maximum cross-section and the corresponding velocity. For small size clusters (), the characteristic parameters show strong variations with the number of atoms in the cluster. Abrupt dips observed for n=10 and n=22 are attributed to electronic properties. Charge transfer patterns observed for various collisional systems present similarities, which appear more sensitive to cluster quantum size effects than to collision energy defects. In their whole, the and v_m parameters show differences in both their size evolution and their absolute values discussed in term of projectile and target electronic structures. Received 13 April 2000 and Received in final form 29 June 2000     

Stability and fragmentation of multiply ionized clusters

The stability of neutral, singly and multiply ionized silicon clusters, (N = 2-7, M = 0, , , ), has been investigated using an ab initio density functional method. We show that the fragmentation effect significantly affects the structure of mass-spectra of multiply ionized silicon clusters. For clusters, the clusters with a large fragmentation energy are found to correspond to the high peaks at N = 4 and 6 in mass-spectra. For clusters, a peak at N = 5 in mass-spectra has been predicted to be especially high. Received: 9 June 1997 / Revised: 8 January 1998 / Accepted: 25 February 1998     

Relative stability of Sin and SinSc- clusters in the range n = 14–18

We present a first-principles pseudopotential optimization of the lower energy equilibrium structure of Si_nSc^- anions for n=14-18. We find that Si₁₆Sc^- is more stable than its neighbors clusters, in agreement with recent experimental mass spectra. We also optimize the geometry of pure Si_n neutral clusters in the range n=14-18, and compare our results with those from previous first-principles calculations.     

The effect of volumetric weighting in the interaction of intense laser fields with clusters

Silver clusters embedded in helium nanodroplets are exposed to intense femtosecond laser pulses (10¹³ - 10¹⁶ W/cm²). The signal of highly charged (q≤11) atomic fragments is maximized by delayed plasmon enhanced ionization using stretched laser pulses. Further details with respect to the dynamics of the charging process can be obtained, when the intensity distribution within the laser focus is taken into account. For the first time, the z-scan method is applied to clusters which offers a route to investigate the explicit dependence of the ion signals with respect to the laser intensity. By taking advantage of the volumetric weighting effect ionization thresholds are determined, yielding values well below 10¹⁴ W/cm² for Ag^q+ ions with q≤11.     

Ionic recoil energies in the Coulomb explosion of metal clusters

The photoionization of metal clusters in intense femtosecond laser fields has been studied. In contrast to an experiment on atoms, the interaction in this case leads to a very efficient and high charging of the particle where tens of electrons per atom are ejected from the cluster. The recoil energy distribution of the atomic fragment ions was measured which in the case of lead clusters exceeds 180 keV. Enhanced charging efficiency which we observed earlier for specific pulse conditions is not reflected in the recoil energy spectra. Both the average and the maximum energies decrease with increasing laser pulse width. This is in good agreement with molecular dynamics calculations. Received 20 December 2000     

Photoionization of argon, krypton and xenon clusters in the inner valence shell region

Photoionization of rare gas clusters in the innervalence shell region has been investigated using threshold photoelectron and photoion spectrometers and synchrotron radiation. Two classes of states are found to play an important role: (A) valence states, correlated to dissociation limits involving an ion with a hole in its innervalence ns shell, (B) Rydberg states correlated to dissociation limits involving an ion with a hole in its outervalence np shell plus an excited neutral atom. In dimers, class A states are “bright”, that is, accessible by photoionization, and serve as an entrance step to form the class B “dark” states; this character fades as the size of the cluster increases. In the dimer, the “Mulliken” valence state is found to present a shallow potential well housing a few vibrational levels; it is predissociated by the class B Rydberg states. During the predissociation a remarkable energy transfer process is observed from the excited ion that loses its innershell electron to its neutral partner. Received: 10 February 1998 / Revised: 17 July 1998 / Accepted: 31 July 1998     

Simulation of small positively charged MgO clusters and study of their electronic densities

With the help of ab initio methods the clusters [(MgO)₁₃Mg] ^Q+ are simulated for Q = 0, 1, 2. Then, vacancy clusters [(MgO)₁₂Mg₂] ^Q+ obtained by removing one oxygen atom are computed for Q running from 0 to 4. These clusters exhibit a slight sphericity and generally shorter interatomic distances than in the crystal. The electronic densities variations are studied in function of Q. In particular, it is observed that the electronic density in the oxygen vacancy goes to a maximum when Q = 2. The ionisation potentials vary from approximately 4 to 14 eV when Q varies from 0 to 3, with a more rapid increase from Q = 1 to Q = 2. The stability study of vacancy clusters show that they experience a phase transition when their charge becomes equal to 2, in accordance with the features mentioned above. Received 14 September 1999 and Received in final form 2 December 1999     

Interpretive approach to collective effects of electrons in multicenter problems

Interpretive theoretical tools prove valuable in guiding the analysis of experiments in the realm of atomic clusters. Here, we review basic elements of an analytic approach that makes it possible to find and visualize the effective electrostatic potential and Coulomb correlations in multicenter problems. To illustrate the utility of these concepts we apply them to exploring molecular-doped metallic clusters. This study is aiming at a systematic, visual assessment of changes induced in screening, Coulomb correlation and effective potential by varying the charge of the electronegative impurity and its position in the cluster cage.     

Geometrical and quantal fragmentation of optical response in

By comparing quantal and semi-classical calculations of optical response, we work out the part of the splitting of the plasmon spectra which is exclusively due to geometrical effects. We apply the analysis to the test case which exhibits an interesting geometry with strong prolate quadrupole deformation and a pronounced asymmetry in addition. We find a new type of resonance splitting which is due to geometrical effects but goes beyond the simple and well known deformation splitting. Received: 6 April 1998 / Accepted: 24 April 1998     

Stability of charged metal clusters in the vicinity of an ionic charge

Stability of highly charged metal clusters in the electric field of an external ion is investigated with the classical liquid drop model. We study the optimum shape of the cluster which has a local minimum of the total energy, taking account of the effects of the surface charge polarization on the Coulomb energy and the cluster deformation on the surface energy. We find that the cluster deformation greatly affects the total energy of the system and that a cluster with a fissility larger than some critical value 0.7-0.8 can become unstable against deformation. We investigate the local competition between the Coulomb force and the surface tension at the cluster surface and show that the surface charge polarization which is induced by the external electric field significantly affects the shape of the cluster and its stability. Received 5 November 2002 / Received in final form 27 January 2003 Published online 11 March 2003 RID="a" ID="a"e-mail: hamada@konan-u.ac.jp     

Temperature-dependent ionization potential of sodium clusters

The ionization potential of sodium clusters () at a finite temperature is studied using density functional theory and ab initio molecular dynamics. The threshold regions of the photoionization efficiency curves are deduced from the integrated IP distributions, which are obtained from the energy eigenvalues of the highest occupied Kohn-Sham states during molecular dynamics by applying a theoretically well-defined shift. The calculated ionization potentials are directly compared to the experimental values. The energetically best geometry of Na₅₅ is found to be a slightly distorted icosahedron. Received 16 April 1999 and Received in final form 6 July 1999     

Effect of melting on ionization potential of sodium clusters

The effect of melting transition on the ionization potential has been studied for sodium clusters with 40, 55, 142, and 147 atoms, using ab initio and classical molecular dynamics. Classical and ab initio simulations were performed to determine the ionization potential of Na₁₄₂ and Na₁₄₇ for solid, partly melted, and liquid structures. The results reveal no correlation between the vertical ionization potential and the degree of surface disorder, melting, or the total energy of the cluster obtained with the ab initio method. However, in the case of 40 and 55 atom clusters, the ionization potential seems to decrease when the cluster melts. Received 1st November 2002 Published online 24 April 2003 RID="a" ID="a"e-mail: ar@phys.jyu.fi     

The influence of O and C doping on the ionization potentials of Li-clusters

The influence of doping of Li-clusters by electronegative O and C atoms on the ionization potentials was investigated. Experimentally, we report ionization potentials for bare Li_n clusters deduced from photoionization efficiency spectra. The values are compared with the results for Li_nO and Li_nC clusters. Observed differences are largely attributed to a quantum size effect caused by the segregated molecular part around the impurity, which changes the electron work function. Theoretically, the Fermi and exchange-correlation energies which enter the work function, are calculated in the frame of the augmented plane wave (APW) method by taking explicitly into account the presence of the molecular core. The other contribution to the work function, the moment of the double layer at the cluster surface, is computed by solving the corresponding Poisson's equation. Received 9 September 1999 and Received in final form 7 February 2000     

Landau fragmentation and deformation effects in dipoSe response of sodium clusters

Recent experimental data on the dipole plasmon in axial sodium clusters Na _N ⁺ with 11 ≤ N ≤ 57 are analyzed within a self-consistent separable random-phase approximation (SRPA) based on the deformed Konh-Sham functional. Good agreement with the data is achieved. The calculations show that, while in light clusters plasmon properties (gross structure and width) are determined mainly by deformation splitting, in medium clusters with N τ 50 the Landau fragmentation becomes decisive. Moreover, in medium clusters shape isomers come to play with contributions to the plasmon comparable with the ground state one. As a result, commonly used methods of the experimental analysis of cluster deformation become useless and correct treatment of cluster shape requires microscopic calculations.