Melting phenomena: effect of composition for 55-atom Ag-Pd bimetallic clusters

Understanding the composition effect on the melting processes of bimetallic clusters is important for their applications. Here, we report the relationship between the melting point and the metal composition for the 55-atom icosahedral Ag-Pd bimetallic clusters by canonical Monte Carlo simulations, using the second-moment approximation of the tight-binding potentials (TB-SMA) for the metal-metal interactions. Abnormal melting phenomena for the systems of interest are found. Our simulation results reveal that the dependence of the melting point on the composition is not a monotonic change, but experiences three different stages. The melting temperatures of the Ag-Pd bimetallic clusters increase monotonically with the concentration of the Ag atoms first. Then, they reach a plateau presenting almost a constant value. Finally, they decrease sharply at a specific composition. The main reason for this change can be explained in terms of the relative stability of the Ag-Pd bimetallic clusters at different compositions. The results suggest that the more stable the cluster, the higher the melting point for the 55-atom icosahedral Ag-Pd bimetallic clusters at different compositions.     

Melting behaviors of icosahedral metal clusters studied by Monte Carlo simulations

We present an atom‐resolved analysis method that traces physical quantities such as the root‐mean‐square bond length fluctuation and coordination number for individual atoms as functions of temperature or time. This method is applied to explain the temperature‐dependent behaviors of three types of Ni_N (N=12,13,14) clusters. The detailed studies for the three types of clusters reveal characteristics as follows: (a) as the temperature increases, all three types of clusters undergo two‐stage melting, irrespective of the existence of vacancy or adatom on the icosahedral surfaces, (b) the melting of icosahedral clusters with vacancy starts with vacancy hopping, which has not been observed for any type of small clusters (N<34), (c) the melting of the icosahedral clusters with adatom (N=14) is initiated by adatom hopping, followed by the site exchange between the adatom and surface atoms. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 380–387, 2000     

Melting and freezing characteristics and structural properties of supported and unsupported gold nanoclusters

Molecular dynamics simulations in conjunction with MEAM potential models have been used to study the melting and freezing behavior and structural properties of both supported and unsupported Au nanoclusters within a size range of 2 to 5 nm. In contrast to results from previous simulations regarding the melting of free Au nanoclusters, we observed a structural transformation from the initial FCC configuration to an icosahedral structure at elevated temperatures followed by a transition to a quasimolten state in the vicinity of the melting point. During the freezing of Au liquid clusters, the quasimolten state reappeared in the vicinity of the freezing point, playing the role of a transitional region between the liquid and solid phases. In essence, the melting and freezing processes involved the same structural changes which may suggest that the formation of icosahedral structures at high temperatures is intrinsic to the thermodynamics of the clusters, rather than reflecting a kinetic phenomenon. When Au nanoclusters were deposited on a silica surface, they transformed into icosahedral structures at high temperatures, slightly deformed due to stress arising from the Au-silica interface. Unlike free Au nanoclusters, an icosahedral solid-liquid coexistence state was found in the vicinity of the melting point, where the cluster consisted of coexisting solid and liquid fractions but retained an icosahedral shape at all times. These results demonstrated that the structural stability in the structures of small Au nanoclusters can be enhanced through interaction with the substrate. Supported Au nanoclusters demonstrated a structural transformation from decahedral to icosahedral motifs during Au island growth, in contrast to the predictions of the minimum-energy growth sequence: icosahedral structures appear first at very small cluster sizes, followed by decahedral structures, and finally FCC structures recovered at very large cluster sizes. The simulations also showed that island shapes are strongly influenced by the substrate, more specifically, the structural characteristic of a Au island is not only a function of size, but also depends on the contact area with the surface, which is controlled by the wetting of the cluster to the substrate.     

Melting of icosahedral gold nanoclusters from molecular dynamics simulations

Molecular dynamics simulations show that gold clusters with about 600-3000 atoms crystallize into a Mackay icosahedron upon cooling from the liquid. A detailed surface analysis shows that the facets on the surface of the Mackay icosahedral gold clusters soften but do not premelt below the bulk melting temperature. This softening is found to be due to the increasing mobility of vertex and edge atoms with temperature, which leads to inter-layer and intra-layer diffusion, and a shrinkage of the average facet size, so that the average shape of the cluster is nearly spherical at melting.     

Melting of icosahedral nickel clusters under hydrostatic pressure

The thermal stabilities and melting behavior of icosahedral nickel clusters under hydrostatic pressure have been studied by constant‐pressure molecular dynamics simulation. The potential energy and Lindemann index are calculated. The overall melting temperature exhibits a strong dependence on pressure. The Lindemann index of solid structure before melting varies slowly and is almost independent of pressure. However, after the clusters melt completely, the Lindemann index at the overall melting point strongly depends on pressure. The overall melting temperature is found to be increasing nonlinearly with increasing pressure, while the volume change during melting decreases linearly with increasing pressure. Under a high pressure and temperature environment, similar angular distributions were found between liquid and solid structures, indicating the existence of a converging local structure. © 2014 Wiley Periodicals, Inc.     

Structural transitions and melting in LJ(74-78) Lennard-Jones clusters from adaptive exchange Monte Carlo simulations

Phase changes in Lennard-Jones (LJ) clusters containing between 74 and 78 atoms are investigated by means of exchange Monte Carlo simulations in the canonical ensemble. The replica temperatures are self-adapted to facilitate the convergence. Although the 74- and 78-atom clusters have icosahedral global minima, the clusters with 75-77 atoms have decahedral ground-state structures and they undergo a structural transition to icosahedral minima before melting. The structural transitions are characterized by quenching and by looking at the Q4 and Q6 orientational bond order parameters. The transition temperatures are estimated to be 0.114, 0.065, and 0.074 reduced units for LJ75, LJ76, and LJ77, respectively. These values, their ordering and the associated latent heats are compared with other estimates based on the harmonic superposition approach.     

Multiple structural transformations in Lennard-Jones clusters: generic versus size-specific behavior

The size-temperature "phase diagram" for Lennard-Jones clusters LJn with sizes up to n=147 is constructed based on the analysis of the heat capacities and orientational bond order parameter distributions computed by the exchange Monte Carlo method. Two distinct types of "phase transitions" accompanied by peaks in the heat capacities are proven to be generic. Clusters with Mackay atom packing in the overlayer undergo a lower-temperature melting (or Mackay-anti-Mackay) transition that occurs within the overlayer. All clusters undergo a higher-temperature transition, which for the three-layer clusters is proven to be the 55-atom-core-melting transition. For the two-layer clusters, the core/overlayer subdivision is ambiguous, so the higher-temperature transition is better characterized as the breaking of the local icosahedral coordination symmetry. A pronounced size-specific behavior can typically be observed at low temperatures and often occurs in clusters with highly symmetric global minima. An example of such behavior is LJ135, which undergoes a low-temperature solid-solid transition, besides the two generic transitions, i.e., the overlayer reconstruction and the core melting.     

The melting behavior of small clusters of atoms

The melting transition of small clusters composed of 13 particles interacting via the Lennard-Jones 6-12 potential has been investigated by means of extensive Metropolis Monte Carlo simulations. The results indicate that whereas the cluster evaporates in vacuum, in a confined pore the cluster undergoes a smooth melting transition from an icosahedral microcrystal to an inhomogeneous liquid phase, with a specific heat peak centered at 39 K. No evidence was found to support the suggestion of a solid-liquid transition in vacuum at 29 K.     

Computer simulation of quantum melting in hydrogen clusters.

We introduce a new criterion, based on multipole dynamical correlations calculated within reptation quantum Monte Carlo, to discriminate between a melting versus freezing behavior in quantum clusters. This criterion is applied to small clusters of para-hydrogen molecules (both pristine and doped with a CO chromophore), for cluster sizes of around twelve molecules. This is a magic size at which para-hydrogen clusters display an icosahedral structure and a large stability. Despite the similar geometric structure of CO@(pH2)12 and (pH2)13, the first system has a rigid, crystalline, behavior; the second behaves more like a superfluid (or, possibly, a supersolid).     

Competing thermal activation mechanisms in the meltinglike transition of Na(N) (N = 135-147) clusters

The meltinglike transition in unsupported icosahedral Na(N)() clusters, with N = 135-147, has been studied by isokinetic molecular dynamics simulations based on an orbital-free version of density functional theory. A maximum in the melting temperature, T(m), is obtained for Na141, while the latent heat, deltaE, and entropy of melting, deltaS, are maximal for Na147. These observations are in close agreement with calorimetric experiments on N clusters. The size evolution of deltaS is rationalized by the emergence of important premelting effects associated with the diffusive motion of atomic vacancies at the cluster surface. The precise location of the maximum in T(m) is explained in terms of two different thermally activated structural instability mechanisms which trigger the meltinglike transition in the size ranges N = 135-141 and N = 141-147, respectively.     

硼氢及客体二十面体簇合物的结构和稳定性

在(n=0、－1、－2、－3、－4)簇合物几何构型及稳定性研究的基础上,进一步对它的各种内含式和外 接式二十面体簇合物(X@B12H122－和XB12H122－,X=H0/＋、Li0/＋、He、Ne、Be0/2＋、Na＋、Mg2＋)进行了优化和 计 算.发现在内含式结构X@ B12H122－中,当X=Li＋、Be2＋、Mg2＋时,构型较稳定;在外接式结构中, XB12H122－(C3v)结构比XB12H122－(C2v)的结构稳定.通过IRC计算,确定XB12H122－(C2v)是X与B12H122－作用生成产物XB12H122－(C3v)的一种过渡态.     

Specific heat and Lindemann-like parameter of metallic clusters: mono- and polyvalent metals

The Brownian-type molecular dynamics simulation is revisited and applied to study the thermal and geometric properties of four mono- and two polyvalent metallic clusters. For the thermal property, we report the specific heat at constant volume CV and study the solid-liquid-like transition by scrutinizing its characteristic. For the geometric property, we calculate the root mean square relative bond-length fluctuation delta as a function of increasing temperature. The thermal change in delta reflects the movement of atoms and hence is a relevant parameter in understanding the phase transition in clusters. The simulated results for the CV of alkali and aluminum clusters whose ground state structures exhibit icosahedral symmetry generally show one phase transition. In contrast, the tetravalent lead is quite often seen to exhibit two phase transitions, a premelting process followed by a progressive melting. In connection with the premelting scenario, it is found here that those (magic number) clusters identified to be of lesser stability (among other stable ones) according to the second energy difference are clusters showing a greater possibility of undergoing premelting process. This energy criterion applies to aluminum clusters nAl=28 and 38. To delve further into the thermal behavior of clusters, we have analyzed also the thermal variation of deltaT and attempted to correlate it with CV(T). It turns out that the premelting (if exist) and melting temperatures of the smaller size clusters (n less, similar 50) extracted from CV do not always agree quantitatively with that deduced from delta.     

Chemical identification of icosahedral structure for cobalt and nickel clusters

Reactions of bare and hydrogenated cobalt and nickel clusters with ammonia and with water are used to determine cluster geometrical structure. Saturation measurements determine the total number of ammonia binding sites on cluster surfaces. A pattern of minima in the number of such sites is found to correlate with the sequence of closed shells and subshells expected for icosahedral packing in the 50- to 120-atom size range (50- to 200-atom range for hydrogenated clusters). In many cases there are 12 sites at the minima, the number that would be expected for preferred ammonia binding sites on closed (sub)shells of icosahedral clusters. The equilibrium adsorption of a single water molecule provides a sensitive measure of changes in cluster-water binding energy. A pattern of binding energy maxima is found, once again correlating with icosahedral structure, but for clusters having one metal atom more than the closed (sub)shells. In general, hydrogenation enhances the patterns of minima and maxima. These observations are explained in terms of the expected nature of ammonia and water binding to icosahedral clusters.     

On the solid- and liquidlike nature of quantum clusters in their ground state

The ground state of pristine clusters of (paraH(2))(N) and (orthoD(2))(N) of size ranging from N=11 to 55 is examined by means of the variational path integral method. The chemical potential is calculated for two different interaction models and it is shown that the location of magic numbers indeed depends on the chosen interaction potential. Density profiles are calculated and reveal the difference between the two isotopes with regards to shell structure. The magnitude of relative pair distance and position fluctuations is used to asses the rigidity of these finite-size quantum systems. A comparison of generic and specific distance fluctuations as a function of cluster size is proposed as a probe of the appearance of rigidity in the clusters. It is found that smaller (paraH(2))(N) clusters are fluidlike and start to display increased rigidity for clusters of size N> or =26, whereas (orthoD(2))(N) clusters of N=13 and N> or =19 are rigid. Small clusters exhibit structures loosely based on an anti-Mackay icosahedral motif. An anti-Mackay to Mackay transition at N=41-42 is suggested.     

Mackay,Anti-Mackay,Double-Mackay,Pseudo-Mackay,and Related Icosahedral Shell Clusters

Mackay introduced two important crystallographic concepts in a short paper published 40 years ago. One is the icosahedral shell structure (iss) consisting of concentric icosahedra displaying fivefold rotational symmetry. The number of atoms contained within these icosahedral shells and subshells agrees well with the magic numbers in rare gas clusters, (C₆₀) _N molecules, and some metal clusters determined by mass spectroscopy or simulated on energy considerations. The cluster of 55 atoms within the second icosahedral shell occurs frequently and has been called Mackay icosahedron, or simply MI, which occurs not only in various clusters, but also in intermetallic compounds and quasicrystals. The second concept is the hierarchic icosahedral structures caused by the presence of a stacking fault in the fcc packing of the successive triangular faces in the iss. For instance, a fault occurs after the ABC layers resulting an ABCB packing. This is, in fact, a hierarchic icosahedral structure of a core icosahedron connected to 12 outer icosahedra by vertex sharing, or an icosahedron of icosahedra (double MI. Contrary to Mackay's iss, a faulted hierarchic icosahedral shell is, in fact, a twinlike face capping of the underlying triangles; it is, therefore, called an anti-Mackay cluster. The hierarchic icosahedral structure in an Al-Mn-Pd icosahedral quasicrystal has a core of body-centered cube rather than an icosahedron and, therefore, is called a pseudo-Mackay cluster. The hierarchic icosahedral structures have been studied separately in the past in the fields of clusters, nanoparticles, intermetallic compounds, and quasicrystals, but the underlying geometry should be the same. In the following a unified geometrical analysis is presented.     

Evolution of the structural stability of large Cu,Ni, Pd,and Ag clusters with size: An analysis within the embedded atom method

A study of the structural stability of clusters made up of a single component has been carried out within the Embedded Atom Method. Perfect icosahedral and cuboctahedral Cu, Ni, Pd, and Ag clusters with up to 5083 atoms have been compared. The icosahedron is found to be the stable structure for small clusters, and a change of structure from icosahedral to cuboctahedral is found as the cluster size increases. A contraction of the interatomic distances results when the cluster size decreases.     

Symmetry and ferromagnetic clusters

Isomers of pure Fe₁₃ and icosahedral Fe₁₂X clusters are studied using the all-electron linear-combination-of-Gaussian-type-orbital (LCGTO) local-density-functional (LDF) methods that allow the spin and geometry of the cluster to be determined self-consistently. The Fe₁₃ ground state is icosahedral. The icosahedral cluster also has the greatest magnetic moment because of increased symmetry-required orbital degeneracy for electrons of different spins. The central atom of the icosahedral iron cluster has been varied to optimize the spin of the cluster keeping the oribital contribution to the magnetic moment quenched. Varying the central atom under this constraint can alter the magnetic moment by more than 20%. Similar studies have begun on 55-atom icosahedral iron clusters.     

First principles study of the carbon-(silicon-) doped La13 clusters

The structural stability and physical properties have been studied for carbon-(silicon-) doped La(13) clusters using DMOL method based on density-functional theory. Doped La(13) clusters prefer to be icosahedron. Substitutional doping with a carbon or silicon impurity makes some clusters closed electronic shell, especially in icosahedral isomers. Substitutional doping of icosahedral La(13) clusters is found to be favorable at surface sites of clusters, especially for Si-doped La(13) cluster, which is very likely to be formed during the doping process. In addition, the structural distortions due to the doping are discussed.     

Binding energy of large icosahedral and cuboctahedral Lennard-Jones clusters

It is widely believed that the lowest energy configurations for small rare gas clusters have icosahedral symmetry. This contrasts with the bulk crystal structures which have cuboctahedral fcc symmetry. It is of interest to understand the transition between this finite and bulk behavior. To model this transition in rare gas clusters we have undertaken optimization studies within the Lennard-Jones pair potential model. Using a combination of Monte Carlo and Partan Search optimization methods, the lowest energy relaxed structures of Lennard-Jones clusters having icosahedral and cuboctahedral symmetry were found. Studies were performed for complete shell clusters ranging in size from one shell having 13 atoms to 14 shells having 10,179 atoms. It was found that the icosahedral structures are lower in energy than the cuboctahedral structures for cluster sizes having 13 shells or fewer. Additional studies were performed using the more accurate Aziz-Chen [HFD-C] pair potential parameterized for argon. The conclusions appear to be relatively insensitive to the form of the potential.     

Novel lattice-searching method for modeling the optimal strain-free close-packed isomers of clusters

Small icosahedral, decahedral, and fcc structures have been studied by unbiased global optimization methods or Wulff construction and Northby lattice methods. Strain-free close-packed structures are not much discussed because the structures are very difficult to optimize and there is no common strain-free close-packed lattice. We propose a new strategy to construct such a lattice containing all possible strain-free close-packed isomers, and by searching the lattice with an efficient method the optimal close-packed structures were modeled. Testing with the Morse potential at rho0=14.0 for cluster size 10相似文献