E. Armon,
A. Bekkerman,
B. Tsipinyuk and
E. Kolodney*
Schulich Faculty of Chemistry, Technion–Israel Institute of Technology, Haifa 3200003, Israel. E-mail: eliko@technion.ac.il
First published on 31st July 2025
Earlier measurements of the velocity correlated cluster emission (VCCE) effect in keV C60− impact with a silver target (Agn+ (n = 1–9), E. Armon, A. Bekkerman, Y. Cohen, J. Bernstein, B. Tsipinyuk, E. Kolodney, Phys. Rev. Lett., 2014, 113, 027604) are extended to much larger clusters up to Ag21+. This is the largest range of cluster size demonstrated so far for the impact induced emission of velocity correlated clusters from any target. Measurements of the kinetic energy distributions (KEDs) of all emitted Agn+ (n = 1–21) cluster ions showed that the VCCE effect as exhibited by a linear increase of most probable energies of the KEDs as a function of cluster size, Emp(n) is still valid up to Ag21+. The extension of the cluster size range was achieved mainly by applying a higher extraction/acceleration voltage for the surface emitted cluster ions than done before for the smaller range (Uext. = 100 V as compared with 25 V before). The preferential collection and detection of some off-normal lower energy cluster ions resulted in a rather moderate and nearly constant 1–2 eV decrease of Emp(n) as compared with the earlier (Uext. = 25 V) results. Namely, the slope of the linear Emp(n) relation (representing the strength of the VCCE effect) is nearly unaffected by the extraction voltages. It shows that the distortion of the VCCE effect by the extraction field is rather small and specifically, that the effect is getting more pronounced with decrease of the extraction voltage. The observation reported here is important for extending the range of validity of the VCCE effect and for gaining a deeper understanding of the nature and size of the (hypothesized) superhot moving precursor serving as the source for all the emitted clusters.
Recently, based on KEDs measurements of complete families of emitted clusters in collisions of keV C60− with different targets, we have reported a novel velocity correlated cluster emission (VCCE) effect.19,24–27 The velocities of all outgoing cluster ions from each given target are correlated in the sense that they reflect a thermal distribution superimposed on a common center-of-mass velocity VCM. This implies gradual shift of the KEDs to higher energies (as given by their most probable energy values Emp(n)) and broadening with increase of cluster size n. The VCCE effect was measured for the metallic targets Au, Ag, Cu and Al (measuring KEDS of Aun+ (n = 1–15),19 Agn+ (n = 1–9),24,26 Cun+ (n = 1–9)19,26 and Aln+ (n = 1–14)25 cluster ions, respectively) and pre-grown thin carbidic films TaC and NbC (measuring KEDS of TanCn+ (n = 1–10)24,26 and NbnCn+ (n = 1–9)26 cluster ions, respectively). The behavior of the KEDs of the emitted clusters following the C60− impact is very different (actually the opposite one) from that observed for the KEDs of both cationic and neutral clusters (sizes typically up to n = 7) emitted following surface impact of a heavy monoatomic ion (at the same energy), where the KEDs are getting narrower and shifted to lower energies with increase in cluster's size.28–35
Each complete family of measured KEDs of clusters emitted from a given target could be reproduced rather well within a model in which the source of the clusters is a superhot precursor moving with a center-of-mass velocity VCM (in the lab coordinate system). This precursor state can be partially or fully detached from the surface.24 The energy distributions of each complete family of emitted clusters of size n are therefore expressed by a shifted Maxwell function f(E, n):
![]() | (1) |
A crucial aspect of the VCCE effect is the time window for the emission following a single C60− impact. In the absence of usable keV ion pulses (sufficiently intense, jitter free and synchronized) on the ps timescale, direct time resolved (pump–probe) measurements are not yet possible. This unfortunate situation implies that the intriguing time resolved dynamics of the thermal spike over the 0.1–20 ps timescale (evolution and decay) is experimentally inaccessible. This includes the dynamics of clusters emission by different mechanisms, energy transfer from the (collapsed) lattice to the electronic sub-system, etc. The VCCE effect is driven by the collisionally induced, exceptional high energy density produced in a relatively shallow subsurface region, at an early stage of the spike dynamics. We believe that the characteristics of the measured KEDS provide sufficient evidence (although indirect) showing that the correlated emission occurs within the first ps from impact, at the very early phase of the evolution of the thermal spike and well before formation of the crater. This is based on the velocity correlation effect itself, the high energies for emission (up to Emp(15) = 10.5 eV for the Ag15+ cluster19 and the very high temperatures consistently derived from the Maxwellians fit (in the range of kT ∼ 1.1–1.3 eV). Also, detailed molecular dynamics (MD) calculations focusing on the sub-ps decay of the atomic number density within the spike nano volume support this conclusion.19,26,27 In case that this decay is too fast to sustain the accumulation of a sufficiently high subsurface pressure as required for triggering the VCCE effect (efficiently propelling the precursor), we indeed measure a weakening of the effect (Al25) down to its complete disappearance (Be27). We have therefore concluded that the VCCE effect is at work within an ultra-short time window of only a few hundred fs, roughly extending between 200–500 fs from impact. The lower bound is dictated by the need for sufficient time for building up an intense pressure spike which is then decaying within the next few hundred fs. The KEDS of the emitted clusters therefore carry unique and valuable experimental information, not available so far, about the extreme conditions and state of matter existing within the first ps of keV ion-solid impact. Note that the ultra-short time scale for the clusters emission also implies that the translational temperature may not be fully equilibrated with the microcanonical internal (vibrational–rotational) temperature.
Using a simple thermodynamic model (the so-called piston mechanism)19 based on the C60− impact induced isochoric heating of the subsurface nanovolume, followed by its adiabatic expansion, we have arrived at a simple relation expressing the accumulated sub-surface pressure P triggering the VCCE effect (movement of the precursor at VCM) in terms of the experimental ε value. The model does not assume any specific equation of state (EOS) for the spike zone and resulted in the relation P = nprec.ε with nprec. as the initial number density of the precursor. We assume nprec. ≈ ncrystal. for the initial non-perturbed target based on the nature of the large clusters. Using this relation combined with the experimental ε value and MD based density calculations we arrive at an average pressure of 5 GPa for gold (over the 200–400 fs time window, from impact).19
The VCCE effect as measured for atomically heavy metallic target (Au, Ag, Cu) is fully developed for the emission of large clusters (n > 5), where the probability of other contributing mechanisms seems negligible. While for a gold target the VCCE effect is nicely obeyed for a relatively large cluster size range (n = 6–15),19 the situation is less convincing for a silver target where measurements are limited to Ag9+. Here we extend the range of Agn+ cluster sizes for which good quality KEDs are measured up to Ag21+ and show that the main characteristic of the correlated emission, a linear Emp(n) relation, remains valid even for the largest size range measured so far. This is achieved by using a higher extraction voltage than before, also demonstrating that the effect of the extraction voltage on the shape of the KEDs is rather minor.
The C60− ion beam was produced using a homemade capillary source.19,26 The electron pickup process inside the capillary is very gentle resulting in fragmentation free C60− beam with mass purity of 99.9% or better. Kinetic energy spread is 0.55 eV or below (measurement limited). The beam is stable (below 2% drift per hour), intense (up to 80 nA at 10 keV) and quiet (current fluctuations below 0.1%). A complete description of the source along with its performance analysis will be published.37 Performance of the C60− gun is superior to existing C60+ sources which are electron ionization based, and it is better suited for the KEDs measurements. During measurement, the beam was rastered over a 0.5 × 0.5 mm2 area with a beam spot diameter of 30 μm and the target was kept at room temperature. The VCCE effect was also studied at elevated surface temperatures (up to 950 K for Nb target) but no change was found.26 All applied potentials are accurate to within 0.1 volt or better as was also verified by measuring directly along the ion path by a precision digital multimeter. About 5–6 repeating KEDs measurements were averaged and smoothed for each cluster. A routinely used control and calibration measurement is that of the energy distribution of thermally emitted (surface ionized) K+ which serves as a near-zero energy marker.
The ionization process of the metallic clusters emitted in the VCCE mode as induced by bombardment with keV C60− ions, is not well understood. Based on a simplified model of the formation of the VCCE effect (the “piston” model),19 we believe that due to the extremely high temperature (kT ∼ 1 eV) of the moving precursor (from which the detected cluster ions are emitted) some fraction of its constituent atoms and clusters are already in an ionized state. During the very early phase of the evolution of the thermal spike (within the 0.2–0.4 ps time-frame after C60− impact, well before formation of the crater), the nanovolume between these clusters (moving in the detector direction) and the undamaged part of the metal target contains atoms with a very high temperature and a greatly reduced averaged density (down to ∼0.5–0.4 of the density of the undamaged metal).19,26 At such a low density, the energy structure of the electron system in this nanovolume should greatly differ from that of the undamaged metal and therefore hinders charge exchange between the metal and cluster ion (e.g., its neutralization). A competing ionization channel via delayed thermionic emission (during flight) from a neutral Agn parent is highly unlikely, since the ionization energies of the neutral Agn clusters are considerably higher (by about 3–4 eV38) than the dissociation energies for the most favorable fragmentation channels: emission of a single neutral atom or a dimer. Note that calculated binding energy per atom for n = 2–13, 20 are about the same (1–3 eV) for both neutral and cationic clusters39 approaching the bulk cohesive energy of 2.95 eV40 for large clusters. Pronounced odd–even ion intensity alternations are observed here (14 keV C60− bombardment) favoring the odd-size cluster ions, similarly to that observed before for sputtered Agn+ cluster ions using a heavy monoatomic projectile (10 keV Xe+ (ref. 41)). The same behavior was also observed and studied in detail in collision induced dissociation (CID)42,43 and photofragmentation measurements.44 It was concluded that these odd–even intensity alternations reflect corresponding alternations in dissociation energies (measured to be in the range of 1–3 eV for n = 2–2543) resulting in preferential stability of the odd-size cluster ions (relative to the even-size ones) against delayed unimolecular (statistical) decay of the thermally excited ions, mainly via monomer evaporation.42–44 We have therefore focused here on the higher intensity odd-numbered clusters. However, in order to get a more complete picture, given a sufficient signal we have also measured KEDs of even numbered clusters (n = 2, 4, 8, 14, 16).
Fig. 2 shows KEDs of odd numbered Agn+ (n = 1, 3, 5, …, 21) clusters emitted from an Ag target following 14 keV impact of C60− projectile ion incident at 45°. Extraction of the cluster ions is along the normal at 100 V. For the sake of clarity, measured KEDs of even numbered Agn+ clusters were omitted and the KEDS in the figure are shown with alternating solid and empty circles. The main features of the KEDs are in line with the earlier measurements of Agn+ (n = 1–9),24,26 namely, gradual shift of the KEDs to higher energies and broadening with increase of cluster size. Extending the range of measured KEDs from the former (n = 1–9) up to n = 21 here was made possible by several improvements in signal to noise ratio but mainly by the higher extraction voltage of 100 V as compared with 25 V before,24,26 thus enabling better collection of ions emitted at off-normal angles. We did not attempt to fit the KEDs with a shifted Maxwellian energy distribution as done for all the former cluster correlated emission measurements (for the purpose of extracting ε and kT). The mixing/superimposing of mainly normally emitted clusters with a preferentially large fraction of off-normal low energy (field bent) trajectories can potentially bias the fitting process and reduce its reliability. Already when visually comparing the present KEDs (100 V) with the former ones (25 V) one can see that the low energy tail is indeed more pronounced now, in line with the field induced enrichment in low energy trajectories. In early studies the weak field 25 V extraction was chosen as a compromise between minimal field induced distortion (mainly on the very low energy side) and good signal to noise ratio enabling the measurement of high quality KEDS of large clusters. The present results support this choice.
Fig. 3 shows the most probable energy values of all measured KEDS from both the earlier measurements (Uext. = 25 V, empty circles)24,26 and the present measurements (Uext. = 100 V, solid triangles), denoted Emp(n, 25 V) and Emp(n, 100 V) correspondingly. The error bars are ± one standard deviation over 5–6 independent repeating measurements. In the earlier Uext. = 25 V measurements, the complete family of KEDs was fitted by shifted Maxwellians (eqn (1)) according with the precursor model resulting in ε = 0.40 eV and kT = 1.13 eV as a single pair of fitting parameters reproducing the experiment reasonably well. The solid line is given by (connecting) the most probable energies Emp(n) for the shifted Maxwellian
![]() | (2) |
Emp(n) = n·ε + 2kT. | (3) |
![]() | ||
Fig. 3 Most probable kinetic energies of emitted Agn+ cluster ions (Emp(n)) following 14 keV impact of a C60− projectile ion with an Ag target at a 45° incidence. Shown are measured values for 25 V and 100 V extraction voltages. The Emp(n, 25 V) values (empty circles) are taken from ref. 24 along with calculated values (thick solid line and dashed line) based on a shifted Maxwellian for the KEDs (moving thermal precursor model, see text). The Emp(n, 100 V) values (solid triangles) are taken from Fig. 1 along with additional even numbered KEDs up to Ag21+. The dense dashed line is a best fit to the measured values. |
Both relations are plotted using the experimentally derived (fitted) ε, kT values as given above. Although the cluster size n is a discrete variable, for the sake of clarity we treat it here as a continuous one. Especially for the asymptotic expression eqn (3), it is convenient for extracting ε value as the slope and the kT value from the y-intercept. Note that Emp(n) values derived from both expressions are quite similar already from n = 7 (and for other targets such as NbC and TaC even for smaller clusters24,26). Eqn (3) is interesting since it presents Emp(n) (in the limit of large n) as a sum of two independent kinetic energy terms. The n·ε term is associated with the CM movement of the precursor (VCM) and the 2kT term is associated with the thermal motion in the CM frame. This separability is useful as it provides a good approximation for large values of n.
The observations obtained from the Emp(n, 100 V) measurements can be summarized as follows: (1) a nearly linear dependence is observed, in line with the precursor model predictions (eqn (2) and (3)). The dense dashed line is a best fit to the measured Emp(n, 100 V) values. (2) The slope of this nearly linear dependence is about the same as that of the Emp(n, 25 V) dependence (especially above n = 4), this implies that the ε value (see eqn (3)) is nearly the same for both dependences (at least as determined by Emp(n) values only). (3) The increased collection and detection of the off-normal lower energy cluster ions by the 100 V extraction resulted in a rather moderate and nearly constant 1–2 eV decrease of Emp(n, 100 V) values as compared with the earlier Emp(n, 25 V) results. (4) The Emp value of the Ag2+ dimer is exceptionally high in both the 25 V and 100 V measurements. The ratio of dimer to monomer Emp values is about 3 for the earlier 25 V measurement and about 2.6 for the present 100 V measurements. We tentatively attribute this behavior to the contribution of an additional mechanism, which is at work mainly for the dimers. A plausible formation mechanism for the Ag2+ dimer is a recombination between two independently emitted atoms moving with a mutual velocity mismatch corresponding to CM kinetic energy below the dimer dissociation energy.45 Such a mechanism usually requires the two-body association to occur very close to the surface, still within the range of the surface potential, in order to stabilize the dimer by draining excess vibrational energy. Interestingly, this irregularly high Emp value of the emitted dimer ion is observed for all three coin metals26 but is most pronounced for Ag. This behavior is different than that observed for emitted carbidic clusters from ultrathin NbC and TaC layers24,26 where the dimers Emp values obey the linear dependence as in eqn (3).
The conclusions derived from these observations are as follows: (1) the VCCE effect as an outcome of the precursor mechanism and as described by shifted Maxwellian KEDs is shown to be valid till Ag21+. This is a meaningful extension of the validity range beyond the former rather limited range of up to Ag9+. This is also the largest range of correlated cluster emission ever observed for any target beyond Au15+ and Al14+ as measured for 14 keV C60− impact with gold19 and aluminum25 targets respectively. We note that large cluster emissions are especially important for demonstrating the VCCE effect. (2) Increasing the extraction voltage leads to increased signal by collecting a larger fraction of clusters emitted at off-normal angles (by bending their initial outgoing trajectories). These clusters are preferentially of lower energy (with respect to the normally emitted ones) resulting in Emp(n, 100 V) values somewhat lower than the Emp(n, 25 V) values. (3) The effect caused by the fraction of off-normal low energy clusters collected by the 100 V extraction is rather similar for all cluster sizes (within error bars) and therefore the down shift of all Emp(n, 100 V) values to lower energies is about the same and the slope of the linear dependence (ε value, reflecting the strength of the VCCE effect) is similar to that exhibited by the Emp(n, 25 V) dependence. This fraction is relatively modest as an increase of the extraction voltage by a factor of 4 results in a down shift of Emp(n) values by only 1–2 eV.
We have already emphasized the importance of extending the clusters size range of KEDs measurements since the VCCE effect is best manifested for the large clusters where possible contributions by other statistical mechanisms are negligible. Also, since the clusters originate from the superhot moving precursor, their maximal size provide information about the nature and size of the precursor. An additional reason to focus on the large clusters is related with the fact that due to their high translational temperatures (kT ∼ 1 eV), the most probable velocities Vmp of the relatively small (up to n = 4–7 per target) emitted clusters are actually strongly mass dependent, approaching asymptotically the precursor velocity VCM for large clusters.26 The shifted Maxwell velocity distribution is given by with
and ε = mVCM2/2. Recall (see Experimental) that within this presentation both V and VCM are nearly parallel and axially directed along the surface normal and the RFA/QMS axis. The extended range as measured here describes the gradual approach of the Vmp value of the large Agn+ clusters to the common precursor velocity. The behavior of large clusters according to the precursor model also emphasizes the fact that an alternative model of gas dynamic acceleration in adiabatic free jet like expansion (seeded beam) cannot properly describe the results. Here we measure an asymptotical convergence of the clusters Vmp values to a final velocity (VCM) with increasing cluster size, while seeded beams are characterized by a velocity mismatch (cluster velocity vs. terminal flow velocity) which is increasing with the mass of the seed species (velocity slip effect).46–49 We also note that the axial velocity distribution usually assumed in free jet expansion49–51 (neglecting perpendicular velocity components) is similar to the shifted Maxwell distribution as given above (for the precursor model) with similar Vmp expression, but with terminal flow velocity instead of the precursor CM velocity (drifting Maxwellian distribution). This implies that similarly to the case of the precursor mechanism and due to the high temperatures involved, also within a free jet model the small clusters range will not show any velocity (Vmp) uniformity but a strong mass dependence.
This journal is © the Owner Societies 2025 |