Structural transition from icosahedra to decahedra of large Lennard-Jones clusters

The lowest icosahedral and decahedral energies of LJ1001-1610 clusters are obtained using a greedy search method (GSM) based on lattice construction. By comparing the lowest energies of icosahedral and decahedral clusters with the same atoms, the structural transition of LJ clusters is studied. Results show that the critical size from icosahedra to decahedra is located at N = 1034. When the cluster size is larger than 1034, the optimal structures are decahedra except the LJ1367-1422 clusters near the magic number, 1402, of icosahedra. However, the energies of icosahedra near the next magic number, 2044, are higher than that of decahedra, which implies that decahedra will be the optimal structure when the cluster size is larger than 1422, even for those clusters near the magic numbers of icosahedra.     

Global optimization of Lennard-Jones clusters by a parallel fast annealing evolutionary algorithm

A parallel fast annealing evolutionary algorithm (PFAEA) was presented and applied to optimize Lennard-Jones (LJ) clusters. All the lowest known minima up to LJ(116) with both icosahedral and nonicosahedral structure, including the truncated octahedron of LJ(38), central fcc tetrahedron of LJ(98), the Marks' decahedron of LJ(75)(-)(77), and LJ(102)(-)(104), were located successfully by the unbiased algorithm. PFAEA is a parallel version of fast annealing evolutionary algorithm (FAEA) that combines the aspect of population in genetic algorithm and annealing algorithm with a very fast annealing schedule. A master-slave paradigm is used to parallelize FAEA to improve the efficiency. The performance of PFAEA is studied, and the scaling of execution time with the cluster size is approximately cubic, which is important for larger scale energy minimization systems.     

Clever and efficient method for searching optimal geometries of lennard-jones clusters

An unbiased algorithm for determining global minima of Lennard-Jones (LJ) clusters is proposed in the present study. In the algorithm, a global minimum is searched by using two operators: one modifies a cluster configuration by moving atoms to the most stable positions on the surface of a cluster and the other gives a perturbation on a cluster configuration by moving atoms near the center of mass of a cluster. The moved atoms are selected by employing contribution of the atoms to the potential energy of a cluster. It was possible to find new global minima for LJ506, LJ521, LJ536, LJ537, LJ538, and LJ541 together with putative global minima of LJ clusters of 10-561 atoms reported in the literature. This indicates that the present method is clever and efficient for cluster geometry optimization.     

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.     

A dynamic lattice searching method with constructed core for optimization of large Lennard-Jones clusters

A variation of the previous dynamic lattice searching (DLS) method, named as DLS with constructed core (DLSc), was proposed for structural optimization of Lennard-Jones (LJ) clusters. In the new method, the starting random structure is generated with an icosahedron or a decahedron as a core. For a cluster with n shells, the atoms in the inner n - 2 shells are set as a fixed core, and the remaining atoms in the outer 2 shells are optimized by DLS. With applications of DLSc to optimization of LJ100-200 and LJ660-670, it was found that all the putative global minima can be obtained by using the DLSc method, and the method was proved to be high efficient compared with the previous DLS, because the searching space is reduced by the use of the fixed core. However, although DLSc is still an unbiased approach for smaller LJ clusters, it turned out to be biased for large ones. Further works are still needed to make it to be a more general method for cluster optimization problem.     

Surveying a potential energy surface by eigenvector-following

We have developed a method to search potential energy surfaces which avoids some of the difficulties associated with trapping in local minima. Steps are directly taken between minima using eigenvector-following. Exploration of this space by low temperature Metropolis Monte Carlo is a useful global optimisation tool. This method successfully finds the lowest energy icosahedral minima of Lennard- Jones clusters from random starting configurations, but cannot find the global minimum in a reasonable time for difficult cases such as the 38-atom Lennard-Jones cluster where the face-centred-cubic truncated octahedron is lowest in energy. However, by performing searches at higher temperatures, we have found a pathway between the truncated octahedron and the lowest energy icosahedral minima. Such a pathway may be illustrative of some of the structural transformations that are observed for supported metal clusters by electron microscopy.     

A dynamic lattice searching method with interior operation for unbiased optimization of large Lennard-Jones clusters

For improving the efficiency of dynamic lattice searching (DLS) method for unbiased optimization of large Lennard-Jones (LJ) clusters, a variant of the interior operation (IO) proposed by Takeuchi was combined with DLS. The method is named as DLS-IO. In the method, the IO moves outer atoms with higher energy toward the coordinates center, i.e., (0, 0, 0), of a cluster and a local minimization (LM) follows each IO. This makes the interior atoms more compact and the outer atoms more uniformly distributed with lower potential energy. Therefore, the starting structure for DLS operations is closer to the global optimum compared with the randomly generated structures. On the other hand, a method to identify the central atom is proposed for the central vacancy problem. Optimizations of LJ(500), LJ(561), LJ(660), LJ(665), and LJ(670) were investigated with the DLS-IO, and the structural transition during the optimization was analyzed. It was found that the method is efficient and unbiased for optimization of large LJ clusters, and it may be a promising approach to be universally used for structural optimizations.     

Structural transitions and melting of copper clusters

Molecular dynamics is used to study the melting and structural transitions of small copper clusters. The melting temperature is found to be proportional to the average coordination number. Small icosahedral clusters melt at slightly higher temperatures than the cubic structures. Small cuboctahedral clusters are not stable but transform via a nondiffusive transition to icosahedral structure.     

PDECO: Parallel differential evolution for clusters optimization

The optimization of the atomic and molecular clusters with a large number of atoms is a very challenging topic. This article proposes a parallel differential evolution (DE) optimization scheme for large‐scale clusters. It combines a modified DE algorithm with improved genetic operators and a parallel strategy with a migration operator to address the problems of numerous local optima and large computational demanding. Results of Lennard–Jones (LJ) clusters and Gupta‐potential Co clusters show the performance of the algorithm surpasses those in previous researches in terms of successful rate, convergent speed, and global searching ability. The overall performance for large or challenging LJ clusters is enhanced significantly. The average number of local minimizations per hit of the global minima for Co clusters is only about 3–4% of that in previous methods. Some global optima for Co are also updated. We then apply the algorithm to optimize the Pt clusters with Gupta potential from the size 3 to 130 and analyze their electronic properties by density functional theory calculation. The clusters with 13, 38, 54, 75, 108, and 125 atoms are extremely stable and can be taken as the magic numbers for Pt systems. It is interesting that the more stable structures, especially magic‐number ones, tend to have a larger energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital. It is also found that the clusters are gradually close to the metal bulk from the size N > 80 and Pt₃₈ is expected to be more active than Pt₇₅ in catalytic reaction. © 2013 Wiley Periodicals, Inc.     

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.     

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.     

Parallel random tunneling algorithm for structural optimization of Lennard-Jones clusters up to N=330

A random tunneling algorithm (RTA) is derived from the terminal repeller unconstrained subenergy tunneling (TRUST) algorithm, and the parallelization of the RTA is implemented with an island parallel paradigm. Combined with the techniques of angular moving, the parallel random tunneling algorithm (PRTA) is applied to the optimization of Lennard-Jones (LJ) atomic clusters, and all the global minima of LJ clusters with size up to 200 are successfully located. For the optimization of larger cluster, a PRTA with an improved seeding technique is developed and successfully applied to the optimization of LJ151-LJ309. Furthermore, the optimized structures of LJ309-330 with the PRTA, which have not been studied before, are also provided.     

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.     

Formation,growth mechanism and packing sequences of binary alloy cluster anions from laser ablation of mixtures of lead and transition metals

By using laser ablation of the mixtures of a transition metal (M: Ti, Cr, Mn, Fe, Co, Ni, Cu, Zn, Pd, Ag) plus lead, M/Pb binary cluster anions were observed except for Zn, and the number of transition metal atoms contained in the binary clusters is at most 4. This behavior is different from that reported previously for M/Ge binary clusters. The experiments indicate that it is also very difficult to form Al/Pb clusters. The distribution patterns of M/Pb binary alloy cluster anions are remarkably similar to those of pure Pb clusters, consistent with a formation mechanism in which transition metal atoms are sequentially attached to [M(x-1)Pb(y)](-) clusters and thus form [M(x)Pb(y)](-) clusters by a simple condensation process. As the number of transition metal atoms increases, the intensities of binary clusters gradually decrease. It is proposed that [MPb(4)](-) and [MPb(5)](-) cluster anions might be the unit building blocks of M/Pb binary cluster anions, and the layer packing sequences for magic clusters are predicted on this basis. The [M(x)Pb(y)](-) binary clusters containing 13 atoms (x + y = 13; x not equal 0) are proposed to have an icosahedral structure.     

Anisotropy effect on global minimum structures of clusters: two-center Lennard-Jones model

Using a two-center Lennard-Jones (2CLJ) model, the simplest anisotropic case, we investigated how anisotropy affects global minimum structures of clusters and obtained some interesting results. The anisotropy parameter, R, is defined as the ratio of the bond length of 2CLJ dimer to the LJ equilibrium pair separation, where a larger R value means higher anisotropy. For low R values, the structures resemble those of the Lennard-Jones atomic clusters. However, as the pairwise interaction becomes more anisotropic, the "magic numbers" change, and several novel cluster patterns emerge as particularly stable structures, and the global minima change from icosahedral, to polyicosahedral and to novel irregular structures. Moreover, increasing the anisotropy effectively softens the 2CLJ potential. Given the general importance of the LJ cluster as a simple model cluster, 2CLJ model can provide a straightforward and useful analysis of the effect of molecular shape on the structures of clusters.     

Continuation calculations of boron- (aluminum-, titanium-, and nickel-) doped La13 clusters

In this work, we have calculated boron-, aluminum-, titanium-, and nickel-doped La13 clusters by DMOL method based on the density-functional theory. Two doping modes are employed: surface and center doping. The boron, aluminum, and nickel atoms prefer to occupy the surface sites while the titanium atom prefers to occupy the center site. The doped La13 clusters with these four kinds of atoms have lower binding energy than pure La13 clusters. The icosahedral isomers are of lower binding energy than cubotahedral and decahedral isomers for La12B(-1), La12Al(-1), and La12Ni, while doping makes the cubotahedral La12Ti stable with a binding energy a little lower than icosahedral La12Ti. There are electronic shell effects in icosahedral La12B(-1) and La12Al(-1). The icosahedral La12B(-1) is promising to be formed during the doped process experimentally. Furthermore, we have also discussed the distorted structures of center doping by bond lengths, density of states, and charge transfers.     

Structural transformations in silver nanoparticles

A molecular dynamics simulation was performed for silver clusters of 147, 309, and 561 atoms with the initial cuboctahedral habit in the temperature range 0–1000 K with an embedded atom potential for silver. Structural transitions of the silver clusters to complex twins (icosahedral habit) with coherent (111)/(111) boundaries over all edges of icosahedra were found, which started at temperatures of 50 K, 350 K, and 700 K, respectively. To analyze the structural transformations in nanoparticles, an algorithm is proposed based on a simplicial Delaunay decomposition (Delaunay triangulation). It was found that after the transition of silver nanoparticles to complex twins, the atomic motion becomes vibrational; the atoms vibrate around the sites that correspond to the vertices of the regular polyhedra. In the case of the 147-atom silver nanoparticle, the polyhedra are arranged in the following sequence, starting from the center of mass: icosahedron (12 atoms), icosododecahedron (30 atoms), icosahedron (12 atoms), dodecahedron (20 atoms), truncated icosahedron (60 atoms, isostructural with fullerene C₆₀), icosahedron (12 atoms), and one atom at the center of mass.     

Transition between icosahedral and cuboctahedral nanoclusters of lead

We have used ab initio methods to study the possible transition between icosahedral (ico) and cuboctahedral (fcc) structures in lead nanoclusters of sizes up to 309 atoms. Spontaneous fcc-to-ico transition in Pb(13) was observed in the ab initio molecular dynamics (MD) simulations at various temperatures. The transition path can be described predominantly by an angular variable s, which can, generally be applied to the similar transitions in clusters of larger sizes and was observed to follow the Mackay model. We have calculated the two-dimensional energy surface that describes the transition in Pb(13) and found a barrierless fcc-to-ico transition path, which is consistent with the observed spontaneous transition in the ab initio MD simulations. The atomic displacements in the transition were identified as one of the vibrational eigenmodes of these two Pb(13) clusters. For clusters of larger sizes (Pb(n), where n = 55, 147, and 309), the possible transitions following similar paths were determined not to be barrierless and the sizes of the barriers were determined by the ab initio elastic band method.     

Temperature and Composition Dependent Structural Evolution: Thermodynamics of CunAg135−n (n = 0–135) Nanoalloys during Cooling

Molecular dynamics simulations are performed to investigate the changes of packing structures, and thermodynamic quantities including internal energy, entropy, and free energy are used to determine temperature regime and transition time of atomic packing structures. The simulation results show different packing structures as the component composition changes, and there are different packing patterns during cooling. For these Cu-Ag alloy clusters containing only a small number of atoms of Cu, they present FCC packing structures in different parts at high temperatures, and then there are transformations to icosahedral structures. With the increase in content of Cu atoms, there is a transition mechanism from molten state to icosahedron. When the content of Cu atoms is appropriate, core-shell structures can be formed at room temperature.     

Structure and shape variations in intermediate-size copper clusters

Using extensive, unbiased searches based on density-functional theory, we explore the structural evolution of Cu(n) clusters over the size range n=8-20. For n=8-16, the optimal structures are plateletlike, consisting of two layers, with the atoms in each layer forming a trigonal bonding network similar to that found in smaller, planar clusters (n相似文献