From sticky-hard-sphere to Lennard-Jones-type clusters

A relation M_{SHS→LJ} between the set of nonisomorphic sticky-hard-sphere clusters M_{SHS} and the sets of local energy minima M_{LJ} of the (m,n)-Lennard-Jones potential V_{mn}^{LJ}(r)=ɛ/n-m[mr^{-n}-nr^{-m}] is established. The number of nonisomorphic stable clusters depends strongly and nontrivially on both m and n and increases exponentially with increasing cluster size N for N≳10. While the map from M_{SHS}→M_{SHS→LJ} is noninjective and nonsurjective, the number of Lennard-Jones structures missing from the map is relatively small for cluster sizes up to N=13, and most of the missing structures correspond to energetically unfavorable minima even for fairly low (m,n). Furthermore, even the softest Lennard-Jones potential predicts that the coordination of 13 spheres around a central sphere is problematic (the Gregory-Newton problem). A more realistic extended Lennard-Jones potential chosen from coupled-cluster calculations for a rare gas dimer leads to a substantial increase in the number of nonisomorphic clusters, even though the potential curve is very similar to a (6,12)-Lennard-Jones potential.

Simulation of the Raman spectroscopy of multi-layered carbon nanomaterials

A empirical potential based model for simulating the Raman spectroscopy of layered carbon nanomaterials is introduced.

Au55, a stable glassy cluster: results of ab initio calculations

Structure and properties of small nanoparticles are still under discussion. Moreover, some thermodynamic properties and the structural behavior still remain partially unknown. One of the best investigated nanoparticles is the Au55 cluster, which has been analyzed experimentally and theoretically. However, up to now, the results of these studies are still inconsistent. Consequently, we have carried out the present ab initio study of the Au55 cluster, using up-to-date computational concepts, in order to clarify these issues. Our calculations have confirmed the experimental result that the thermodynamically most stable structure is not crystalline, but it is glassy. The non-crystalline structure of this cluster was validated by comparison of the coordination numbers with those of a crystalline cluster. It was found that, in contrast to bulk materials, glass formation is connected to an energy release that is close to the melting enthalpy of bulk gold. Additionally, the surface energy of this cluster was calculated using two different theoretical approaches resulting in values close to the surface energy for bulk gold. It shall be emphasized that it is now possible to give a confidence interval for the value of the surface energy.

Solvation of Li+ by argon: how important are three-body forces?

A global geometry search on a new potential energy surface for Li⁺Ar_n clusters revealed that three-body interactions must be included to reproduce ab initio structures and accurate energetic features.

Relationship and discrepancies between the Extended-Rydberg and the Generalized Buckingham potential energy functions

This paper relates the interatomic energy according to the Extended-Rydberg and the Generalized Buckingham potential functions by applying a Maclaurin series expansion on the latter and thereafter comparing it with the former. In so doing, the plotted curves of these two functions not only show equal curvature at the equilibrium state, but also reveal a discrepancy for the finite distortion. It is shown that, when equated at equilibrium, the Generalized Buckingham gives a lower energy in comparison to the Extended-Rydberg at finite bond compression and stretching. However, the energy difference diminishes when the interatomic distance exceeds twice that at equilibrium. Due to such discrepancies upon comparative normalization, it would be beneficial for computational chemists to select the appropriate potential function for the purpose of conservative molecular modeling. .

Shell and subshell periodic structures of icosahedral nickel nanoclusters

Using the modified analytic embedded atom method and molecular dynamics, the binding energies and their second order finite differences (stability functions) of icosahedral Ni clusters with shell and subshell periodicity are studied in detail via atomic evolution. The results exhibit shell and subshell structures of the clusters with atoms from 147 to 250,000, and the atomic numbers corresponding to shell or subshell structures are in good agreement with the experimental magic numbers obtained in time-of-flight mass spectra of threshold photoionization, and Martin's theoretical proposition of progressive formation of atomic umbrellas. Clusters with size from 147 to 561 atoms are energetically investigated via one-by-one atomic evolution and their magic numbers are theoretically proved. For medium-size Ni clusters with 561 to 2057 atoms, the prediction of magic numbers with atomic numbers is performed on the basis of umbrella-like subshell growth in near face-edge-vertex order. The similarity of the energy curves makes it possible to extend the prediction to even larger Ni nanoclusters in hierarchical Mackay icosahedral configurations.

Strategies for increasing the efficiency of a genetic algorithm for the structural optimization of nanoalloy clusters

An improved genetic algorithm (GA) is described that has been developed to increase the efficiency of finding the global minimum energy isomers for nanoalloy clusters. The GA is optimized for the example Pt12Pd12, with specific investigation of: the effect of biasing the initial population by seeding; the effect of removing specified clusters from the population ("predation"); and the effect of varying the type of mutation operator applied. These changes are found to significantly enhance the efficiency of the GA, which is subsequently demonstrated by the application of the best strategy to a new cluster, namely Pt19Pd19.