Tight binding molecular dynamics study of Ni clusters

Data-driven stabilization of NimPdn-m nanoalloys: a study using density functional theory and data mining approaches

Green hydrogen, generated through the electrolysis of water, is a viable alternative to fossil fuels, although its adoption is hindered by the high costs associated with the catalysts. Among a wide variety of potential materials, binary nickel-palladium (NiPd) systems have garnered significant attention, particularly at the nanoscale, for their efficacious roles in catalyzing hydrogen and oxygen evolution reactions. However, our atom-level understanding of the descriptors that drive their energetic stability at the nanoscale remains largely incomplete. Here, we investigate by density functional theory calculations the descriptors that drives the stability of the NimPdn-m clusters for different sizes (n = 13, 27, 41) and compositions. To achieve our goals, a large number of trial configurations were generated and selected using data mining algorithms (k-means, t-SNE) and genetic algorithms, while the most important physical-chemical descriptors were identified using Spearman correlation analysis. We have found that core-shell formation, with the smaller Ni atoms lying in the center of the particle, plays a major role in the stabilization of the nanoalloys, and this effect causes the alloys to assume a icosahedral-fragment configuration (as the unary nickel cluster) instead of a fcc fragment (as the unary palladium cluster). However, the core-shell formation in this alloy is unique in that Pd poor compositions exhibit scattered Pd atoms on the surface. As the palladium content increases, this gives rise to the complete Pd shell. This stabilization mechanism is quantitatively supported by the different correlations observed in the number of Ni-Ni and Pd-Pd bonds with energy, in which the latter tends to decrease alloy stability. Furthermore, a notable trend is the correlation between the coordination number of Ni atoms with alloy stabilization, while the coordination of Pd atoms shows an inverse correlation.

Structure and size-dependent vibrational and thermal properties of Ni clusters: A systematic ab initio approach

There is scarce information on the vibrational and thermal properties of small Ni clusters. Here, the outcomes of ab initio spin-polarized density functional theory calculations on the size and geometry effects upon the vibrational and thermal properties of Ni_n (n = 13 and 55) clusters, are discussed. For theses clusters a comparison is presented between the closed shell symmetric octahedral (O_h) and the icosahedral (I_h) geometries. The results indicate that the I_h isomers are lower in energy. Besides, ab initio molecular dynamics runs at T = 300K show that Ni₁₃ and Ni₅₅ clusters transform from their initial O_h geometries towards the corresponding I_h ones. For Ni₁₃, we also consider the lowest energy less symmetric layered 1-3-6-3 structure, and the cuboid, recently observed experimentally for Pt₁₃, which is competitive in energy but is unstable, as phonon analysis reveals. We calculate their vibrational density of states (νDOS) and heat capacity, and compare with the Ni FCC bulk counterpart. The characteristic features of the νDOS curves of these clusters are interpreted in terms of the clusters' sizes, the interatomic distance contractions, the bond order values as well as the internal pressure and strains of the clusters. We find that the softest possible frequency of the clusters is size and structure-dependent, being the smallest for the O_h ones. We identify mostly shear, tangential type displacements involving mainly surface atoms for the lowest frequency of the spectra of both I_h and O_h isomers. For the maximum frequencies of these clusters the central atom shows anti-phase movements against groups of nearest neighbor atoms. An excess of heat capacity at low temperatures with respect to the bulk is found, while at high temperatures a constant limiting value, close but lower to the Dulong and Petit value, is determined.

Evaluating the performance of ReaxFF potentials for sp2 carbon systems (graphene, carbon nanotubes, fullerenes) and a new ReaxFF potential

We study the performance of eleven reactive force fields (ReaxFF), which can be used to study sp2 carbon systems. Among them a new hybrid ReaxFF is proposed combining two others and introducing two different types of C atoms. The advantages of that potential are discussed. We analyze the behavior of ReaxFFs with respect to 1) the structural and mechanical properties of graphene, its response to strain and phonon dispersion relation; 2) the energetics of (n, 0) and (n, n) carbon nanotubes (CNTs), their mechanical properties and response to strain up to fracture; 3) the energetics of the icosahedral C60 fullerene and the 40 C40 fullerene isomers. Seven of them provide not very realistic predictions for graphene, which made us focusing on the remaining, which provide reasonable results for 1) the structure, energy and phonon band structure of graphene, 2) the energetics of CNTs versus their diameter and 3) the energy of C60 and the trend of the energy of the C40 fullerene isomers versus their pentagon adjacencies, in accordance with density functional theory (DFT) calculations and/or experimental data. Moreover, the predicted fracture strain, ultimate tensile strength and strain values of CNTs are inside the range of experimental values, although overestimated with respect to DFT. However, they underestimate the Young's modulus, overestimate the Poisson's ratio of both graphene and CNTs and they display anomalous behavior of the stress - strain and Poisson's ratio - strain curves, whose origin needs further investigation.

Application of Optimization Algorithms in Clusters

The structural characterization of clusters or nanoparticles is essential to rationalize their size and composition-dependent properties. As experiments alone could not provide complete picture of cluster structures, so independent theoretical investigations are needed to find out a detail description of the geometric arrangement and corresponding properties of the clusters. The potential energy surfaces (PES) are explored to find several minima with an ultimate goal of locating the global minima (GM) for the clusters. Optimization algorithms, such as genetic algorithm (GA), basin hopping method and its variants, self-consistent basin-to-deformed-basin mapping, heuristic algorithm combined with the surface and interior operators (HA-SIO), fast annealing evolutionary algorithm (FAEA), random tunneling algorithm (RTA), and dynamic lattice searching (DLS) have been developed to solve the geometrical isomers in pure elemental clusters. Various model or empirical potentials (EPs) as Lennard-Jones (LJ), Born-Mayer, Gupta, Sutton-Chen, and Murrell-Mottram potentials are used to describe the bonding in different type of clusters. Due to existence of a large number of homotops in nanoalloys, genetic algorithm, basin-hopping algorithm, modified adaptive immune optimization algorithm (AIOA), evolutionary algorithm (EA), kick method and Knowledge Led Master Code (KLMC) are also used. In this review the optimization algorithms, computational techniques and accuracy of results obtained by using these mechanisms for different types of clusters will be discussed.

Uniaxial magnetic anisotropy energy of bimetallic Co-Ni clusters from a first-principles perspective

Along with the growing precision in the control of matter at increasingly smaller size scales, a field of research, based onto magnetic materials of technical interest, such as bimetallic clusters, has been developed in very recent years. Thereby, here, we report on a complete study of bimetallic clusters composed of cobalt and nickel with up to 7 atoms using ab initio methods in the GGA approach. We applied an unrestricted search method based on the tensor of inertia eigenvalues to find the most stable configurations of the clusters, obtaining a diverse set of structures with different geometric properties. We explored the effect of composition on the structural properties, the chemical stability, the magnetization and the magnetic anisotropy energy (MAE) of the so-obtained systems. Our results indicate that the behavior of the clusters is mainly governed by the Co-Co interaction and to a lesser extent by the Co-Ni and Ni-Ni interactions. Furthermore, for a given cluster size the magnetic moment increases by 2 μB/Co-substitution plus 1 μB/Ni-substitution coming from the cobalt and nickel core d-states, while in some cases unpaired hybrid s-electrons can also give rise to itinerant magnetism. These features have been analyzed with the help of a Jellium model and have important consequences for the magnetism and the magnetic anisotropy of the clusters. The magnetic behavior and MAE present complex and intriguing landscapes, which suggests the possibility of finely controlling the magnetic states, by tuning the cluster composition, aiming at technical implementation in fields such as molecular magnetism or quantum computation. In particular, cases such as Co6Ni, Co4Ni3, CoNi6 Co2Ni5 and Co3Ni3 present high relative stability and enhanced magnetic moments (around 10 μB), what makes them promising candidates for applications such as subnanometer magnetic information storage.

Thermodynamics at the nanoscale: phase diagrams of nickel-carbon nanoclusters and equilibrium constants for phase transitions

Using reactive molecular dynamics simulations, the melting behavior of nickel-carbon nanoclusters is examined. The phase diagrams of icosahedral and Wulff polyhedron clusters are determined using both the Lindemann index and the potential energy. Formulae are derived for calculating the equilibrium constants and the solid and liquid fractions during a phase transition, allowing more rational determination of the melting temperature with respect to the arbitrary Lindemann value. These results give more insight into the properties of nickel-carbon nanoclusters in general and can specifically be very useful for a better understanding of the synthesis of carbon nanotubes using the catalytic chemical vapor deposition method.

Melting and Glass Transition for Ni Clusters

The melting of NiN clusters (N = 29, 50-150) has been investigated by using molecular dynamics (MD) simulations with a quantum corrected Sutton-Chen (Q-SC) many-body potential. Surface melting for Ni147, direct melting for Ni79, and the glass transition for Ni29 have been found, and those melting points are equal to 540, 680, and 940 K, respectively. It shows that the melting temperatures are not only size-dependent but also a symmetrical structure effect; in the neighborhood of the clusters, the cluster with higher symmetry has a higher melting point. From the reciprocal slopes of the caloric curves, the specific heats are obtained as 4.1 kB per atom for the liquid and 3.1 kB per atom for the solid; these values are not influenced by the cluster size apart in the transition region. The calculated results also show that latent heat of fusion is the dominant effect on the melting temperatures (Tm), and the relationship between S and L is given.

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.

Broken-symmetry unrestricted hybrid density functional calculations on nickel dimer and nickel hydride

In the present work we investigate the adequacy of broken-symmetry unrestricted density functional theory for constructing the potential energy curve of nickel dimer and nickel hydride, as a model for larger bare and hydrogenated nickel cluster calculations. We use three hybrid functionals: the popular B3LYP, Becke's newest optimized functional Becke98, and the simple FSLYP functional (50% Hartree-Fock and 50% Slater exchange and LYP gradient-corrected correlation functional) with two basis sets: all-electron (AE) Wachters+f basis set and Stuttgart RSC effective core potential (ECP) and basis set. We find that, overall, the best agreement with experiment, comparable to that of the high-level CASPT2, is obtained with B3LYP/AE, closely followed by Becke98/AE and Becke98/ECP. FSLYP/AE and B3LYP/ECP give slightly worse agreement with experiment, and FSLYP/ECP is the only method among the ones we studied that gives an unacceptably large error, underestimating the dissociation energy of Ni(2) by 28%, and being in the largest disagreement with the experiment and the other theoretical predictions. We also find that for Ni(2), the spin projection for the broken-symmetry unrestricted singlet states changes the ordering of the states, but the splittings are less than 10 meV. All our calculations predict a deltadelta-hole ground state for Ni(2) and delta-hole ground state for NiH. Upon spin projection of the singlet state of Ni(2), almost all of our calculations: Becke98 and FSLYP both AE and ECP and B3LYP/AE predict (1)(d(A)(x(2)-y(2)d(B)(x(2)-y(2)) or (1)(d(A)(xy) (d)(B)(xy)) ground state, which is a mixture of (1)Sigma(g) (+) and (1)Gamma(g). B3LYP/ECP predicts a (3)(d(A)(x(2)-y(2))d(B)(xy) (mixture of (3)Sigma(g) (-) and (3)Gamma(u)) ground state virtually degenerate with the (1)(d(A)(x(2)-y(2)d(B)(x)(2)-y(2)/(1)(d(A)(xy)D(B)(xy) state. The doublet delta-hole ground state of NiH predicted by all our calculations is in agreement with the experimentally predicted (2)Delta ground state. For Ni(2), all our results are consistent with the experimentally predicted ground state of 0(g) (+) (a mixture of (1)Sigma(g) (+) and (3)Sigma(g) (-)) or 0(u) (-) (a mixture of (1)Sigma(u) (-) and (3)Sigma(u) (+)).

First-principles study of the electronic structures of icosahedral TiN (N=13,19,43,55) clusters

We have studied the electronic structures of icosahedral Ti(N) clusters (N=13, 19, 43, and 55) by using a real-space first-principles cluster method with generalized gradient approximation for exchange-correlation potential. The hexagonal close-packed and fcc close-packed clusters have been studied additionally for comparisons. It is found that the icosahedral structures are the most stable ones except for Ti(43), where fcc close-packed structure is favorable in energy. We present and discuss the variation of bond length, the features of the highest occupied molecular orbitals and the lowest unoccupied molecular orbital, the evolution of density of states, and the magnetic moment in detail. The results are in good agreement with the predictions from the collision-induced dissociation and size-selected anion photoelectron spectroscopy experiments.

Degradation of inter-atomic bonds during structural phase change in intermediate Ni-clusters (Ni39–Ni49)

Results based on a symmetry- and spin-unrestricted tight-binding molecular-dynamics study are presented for the ground-state geometries of intermediate Ni(n), n in [39,49], clusters. A structural phase change is found to take place around n=43 during which a structural transition from fcc/hcp structure to icosahedral one is observed. This is in good agreement with recent experimental findings. This structural transition is found to be associated with a degradation of the inter-atomic bond energy which indicates that the inter-atomic bond does not only depend on the coordination number of each atom but also on its point group symmetry.

Catalytic action of Ni atoms in the formation of carbon nanotubes: a molecular dynamics study

Catalytic action of Ni atoms in the growth of single-wall carbon nanotubes is investigated using tight-binding molecular dynamics and ab initio methods. Our results demonstrate this to be a two step process in which the Ni atom first creates and stabilizes defects in nanotubes. The subsequent incorporation of incoming carbon atoms anneals the Ni-stabilized defects freeing the Ni atom to repeat the catalytic process.