Modeling calcium and strontium clusters with many-body potentials

Ca n neutral clusters: a two-step G 0 W 0 and DFT benchmark

Electronic and structural properties of calcium clusters with a varying size range of 2-20 atoms are studied using a two-step scheme within the GW and density functional theory (DFT) with generalized gradient approximation (GGA). The GGA overestimates the binding energies, optimized geometries, electron affinities, and ionization potentials reported in the benchmark. The ground-state structure geometry and binding energy were obtained from the DFT for the ground-state structure of each cluster. The binding energy of the neutral clusters of the calcium series follows an increasing trend, except for a few stable even and odd clusters. The electronic properties of the calcium cluster were studied with an all-electron FHI-aims code. In the G 0 W 0 calculation, the magic cluster Ca10 has relatively high ionization potential and low electron affinity. The obtained ionization potentials from the G 0 W 0 @PBE calculation showed that the larger cluster has less variation, whereas the electron affinities of the series have an increasing trend. The ionization potentials from the G 0 W 0 benchmark for the calcium cluster series have not yet been described in the literature.

A DFT-based genetic algorithm search for AuCu nanoalloy electrocatalysts for CO₂ reduction

Using a DFT-based genetic algorithm (GA) approach, we have determined the most stable structure and stoichiometry of a 309-atom icosahedral AuCu nanoalloy, for potential use as an electrocatalyst for CO2 reduction. The identified core-shell nano-particle consists of a copper core interspersed with gold atoms having only copper neighbors and a gold surface with a few copper atoms in the terraces. We also present an adsorbate-dependent correction scheme, which enables an accurate determination of adsorption energies using a computationally fast, localized LCAO-basis set. These show that it is possible to use the LCAO mode to obtain a realistic estimate of the molecular chemisorption energy for systems where the computation in normal grid mode is not computationally feasible. These corrections are employed when calculating adsorption energies on the Cu, Au and most stable mixed particles. This shows that the mixed Cu135@Au174 core-shell nanoalloy has a similar adsorption energy, for the most favorable site, as a pure gold nano-particle. Cu, however, has the effect of stabilizing the icosahedral structure because Au particles are easily distorted when adding adsorbates.

The effect of the size and shape on the bond number of quantum dots and its relationship with thermodynamic properties

The size- and shape-related B_a/B_t functions.

Structures, spectra, and energies of niobium clusters from Nb13 to Nb20

A comprehensive theoretical investigation on structures and properties of niobium clusters in the range from 13 to 20 atoms, in three different charged states, is performed by using the BPW91 and M06 functionals and the cc-pVDZ-PP basis set. These species are predicted to prefer low spin ground state, i.e., singlet (for even electron) and doublet (for odd electron) systems. In terms of growth mechanism, a compact structure with one Nb encapsulated by a cage formed from five and six triangles is found to be favored over an icosahedral evolution. Unlike many 3d metals, whose volumes are much smaller, 13 and 19 Nb atoms clusters do not exist as icosahedra and double-icosahedra. A distinct case is Nb(15) as it bears a slightly distorted bcc structure. For some systems, several lower lying isomers are computed to be so close in energy that DFT computations cannot clearly establish their ground electronic states. The existence of structural isomers with comparable energy content is established for Nb(n) species with n = 13, 18, 19, and 20 in both neutral and charged states. The vibrational (IR) spectra are also calculated. While the spectra of smaller systems are strongly dependent on addition or removal of an electron from the neutral, the spectra of the larger size clusters are mostly independent of the charged state. The neutrals and their corresponding ions usually have a quite similar IR pattern. Electron affinities (EA), ionization energies (IE), average binding energies, dissociation energies, and frontier orbital energy gaps are evaluated. The computed EAs and IEs are generally in fair agreement with experiment. The Nb(15) system is observed to be stable and it can form a highly symmetric structure in all charged states with both open and closed electron shells.

Accurate ab Initio Binding Energies of Alkaline Earth Metal Clusters

The effects of basis set superposition error (BSSE) and core-correlation on the electronic binding energies of alkaline earth metal clusters Y(n) (Y = Be, Mg, Ca; n = 2-4) at the Moller-Plesset second-order perturbation theory (MP2) and the single and double coupled cluster method with perturbative triples correction (CCSD(T)) levels are examined using the correlation consistent basis sets cc-pVXZ and cc-pCVXZ (X = D, T, Q, 5). It is found that, while BSSE has a negligible effect for valence-electron-only-correlated calculations for most basis sets, its magnitude becomes more pronounced for all-electron-correlated calculations, including core electrons. By utilizing the negligible effect of BSSE on the binding energies for valence-electron-only-correlated calculations, in combination with the negligible core-correlation effect at the CCSD(T) level, accurate binding energies of these clusters up to pentamers (octamers in the case of the Be clusters) are estimated via the basis set extrapolation of ab initio CCSD(T) correlation energies of the monomer and cluster with only the cc-pVDZ and cc-pVTZ sets, using the basis set and correlation-dependent extrapolation formula recently devised. A comparison between the CCSD(T) and density functional theory (DFT) binding energies is made to identify the most appropriate DFT method for the study of these clusters.

A dynamic lattice searching method for fast optimization of Lennard-Jones clusters

A highly efficient unbiased global optimization method called dynamic lattice searching (DLS) was proposed. The method starts with a randomly generated local minimum, and finds better solution by a circulation of construction and searching of the dynamic lattice (DL) until the better solution approaches the best solution. The DL is constructed adaptively based on the starting local minimum by searching the possible location sites for an added atom, and the DL searching is implemented by iteratively moving the atom located at the occupied lattice site with the highest energy to the vacant lattice site with the lowest energy. Because the DL can greatly reduce the searching space and the number of the time-consuming local minimization procedures, the proposed DLS method runs at a very high efficiency, especially for the clusters of larger size. The performance of the DLS is investigated in the optimization of Lennard-Jones (LJ) clusters up to 309 atoms, and the structure of the LJ(500) is also predicted. Furthermore, the idea of dynamic lattice can be easily adopted in the optimization of other molecular or atomic clusters. It may be a promising approach to be universally used for structural optimizations in the chemistry field.