A Monte Carlo study of the structure and thermodynamic behaviour of small Lennard-Jones clusters

Melting and Freezing of Metal Clusters

Recent developments allow heat capacities to be measured for size-selected clusters isolated in the gas phase. For clusters with tens to hundreds of atoms, the heat capacities determined as a function of temperature usually have a single peak attributed to a melting transition. The melting temperatures and latent heats show large size-dependent fluctuations. In some cases, the melting temperatures change by hundreds of degrees with the addition of a single atom. Theory has played a critical role in understanding the origin of the size-dependent fluctuations, and in understanding the properties of the liquid-like and solid-like states. In some cases, the heat capacities have extra features (an additional peak or a dip) that reveal a more complex behavior than simple melting. In this article we provide a description of the methods used to measure the heat capacities and provide an overview of the experimental and theoretical results obtained for sodium and aluminum clusters.

Surface Melting and Surface Diffusion on Clusters

Surface melting on clusters is investigated by a combination of analytic modeling and computer simulation. Homogeneous, argon-like clusters bound by Lennard-Jones forces and Cu-like clusters bound by ‘embedded atom’ potentials are the systems considered. Molecular dynamics (MD) calculations have been carried out for clusters with 40–147 atoms. Well below the bulk melting temperature, the surfaces become very soft, exhibiting well-defined diffusion constants even while the cores remain nearly rigid and solid-like. The simulations, particularly animations, of atomic motion reveal that the surface melting is associated not with amorphous, random surface structures in constant, irregular motion, but rather in large-amplitude, organized, collective motion of most of the surface atoms accompanied by a few “floaters” and holes. At any time, a few of the surface atoms move out of the surface layer, leaving vacancies; these promoted particles wander diffusively, the holes also but less so, and occasionally exchange with atoms in the surface layer. This result is the basis for an analytic, statistical model. The caloric curves, particularly the latent heats, show that surface melting of clusters is a “phase change” different from the bulk melting of clusters.

Metal clusters with hidden ground states: Melting and structural transitions in Al115(+), Al116(+), and Al117(+)

Heat capacities measured as a function of temperature for Al(115)(+), Al(116)(+), and Al(117)(+) show two well-resolved peaks, at around 450 and 600 K. After being annealed to 523 K (a temperature between the two peaks) or to 773 K (well above both peaks), the high temperature peak remains unchanged but the low temperature peak disappears. After considering the possible explanations, the low temperature peak is attributed to a structural transition and the high temperature peak to the melting of the higher enthalpy structure generated by the structural transition. The annealing results show that the liquid clusters freeze exclusively into the higher enthalpy structure and that the lower enthalpy structure is not accessible from the higher enthalpy one on the timescale of the experiments. We suggest that the low enthalpy structure observed before annealing results from epitaxy, where the smaller clusters act as a nucleus and follow a growth pattern that provides access to the low enthalpy structure. The solid-to-solid transition that leads to the low temperature peak in the heat capacity does not occur under equilibrium but requires a superheated solid.

Heterophasic oscillations in nanoscale systems

Both the energy differences between metastable and stable phases and the energy barriers separating these phases decrease with decreasing particle number. Then, for small enough systems, random heterophasic fluctuations of the entire system become an observable form of thermal motion. We discuss mechanisms and observation conditions for these random transitions between phases.

Metal clusters that freeze into high energy geometries

Heat capacities measured for isolated aluminum clusters show peaks due to melting. For some clusters with around 60 and 80 atoms there is a dip in the heat capacities at a slightly lower temperature than the peak. The dips have been attributed to structural transitions. Here we report studies where the clusters are annealed before the heat capacity is measured. The dips disappear for some clusters, but in many cases they persist, even when the clusters are annealed to well above their melting temperature. This indicates that the dips do not result from badly formed clusters generated during cluster growth, as originally suggested. We develop a simple kinetic model of melting and freezing in a system consisting of one liquidlike and two solidlike states with different melting temperatures and latent heats. Using this model we are able to reproduce the experimental results including the dependence on the annealing conditions. The dips result from freezing into a high energy geometry and then annealing into the thermodynamically preferred solid. The thermodynamically preferred solid has the higher freezing temperature. However, the liquid can bypass freezing into the thermodynamically preferred solid (at high cooling rates) if the higher energy geometry has a larger freezing rate.

Ion calorimetry: Using mass spectrometry to measure melting points

Calorimetry measurements have been used to probe the melting of aluminum cluster cations with 63 to 83 atoms. Heat capacities were determined as a function of temperature (from 150 to 1050 K) for size-selected cluster ions using an approach based on multicollision-induced dissociation. The experimental method is described in detail and the assumptions are critically evaluated. Most of the aluminum clusters in the size range examined here show a distinct peak in their heat capacities that is attributed to a melting transition (the peak is due to the latent heat). The melting temperatures are below the bulk melting point and show enormous fluctuations as a function of cluster size. Some clusters (for example, n = 64, 68, and 69) do not show peaks in their heat capacities. This behavior is probably due to the clusters having a disordered solid-like phase, so that melting occurs without a latent heat.

Equilibrium thermodynamics from basin-sampling

We present a "basin-sampling" approach for calculation of the potential energy density of states for classical statistical models. It combines a Wang-Landau-type uniform sampling of local minima and a novel approach for approximating the relative contributions from local minima in terms of the volumes of basins of attraction. We have employed basin-sampling to study phase changes in atomic clusters modeled by the Lennard-Jones potential and for ionic clusters. The approach proves to be efficient for systems involving broken ergodicity and has allowed us to calculate converged heat capacity curves for systems that could previously only be treated using the harmonic superposition approximation. Benchmarks are also provided by comparison with parallel tempering and Wang-Landau simulations, where these proved feasible.

Melting, premelting, and structural transitions in size-selected aluminum clusters with around 55 atoms

Heat capacities have been determined for unsupported aluminum clusters, Al49(+) - Al63(+), from 150 to 1050 K. Peaks in the heat capacities due to melting occur between 450 and 650 K (well below the bulk melting point of 933 K). The peaks for Al+51 and Al+52 are bimodal, suggesting the presence of a premelting transition where the surface of the clusters melts around 100 K before the core. For clusters with n > 55 the melting temperatures suddenly drop, and there is a dip in the heat capacities due to a transition between two solid forms before the clusters melt.

On the premelting features in sodium clusters

Melting in Na(n) clusters described with an empirical embedded-atom potential has been reexamined in the size range 55</=n</=147 with a special attention at sizes close to 130. Contrary to previous findings, premelting effects are also present at such medium sizes, and they turn out to be even stronger than the melting process itself for Na(133) or Na(135). These results indicate that the empirical potential is qualitatively inadequate to model sodium clusters.

Hot and solid gallium clusters: too small to melt

A novel multicollision induced dissociation scheme is employed to determine the energy content for mass-selected gallium cluster ions as a function of their temperature. Measurements were performed for Ga(+)(n) (n=17 39, and 40) over a 90-720 K temperature range. For Ga+39 and Ga+40 a broad maximum in the heat capacity-a signature of a melting transition for a small cluster-occurs at around 550 K. Thus small gallium clusters melt at substantially above the 302.9 K melting point of bulk gallium, in conflict with expectations that they will remain liquid to below 150 K. No melting transition is observed for Ga+17.