Are zinc clusters really amorphous? A detailed protocol for locating global minimum structures of clusters

Investigating structural and electronic properties of neutral zinc clusters: a G0W0 and G0W0Г0(1) benchmark

The structural and electronic properties of zinc clusters (Znn) for a size range of n = 2-15 are studied using density functional theory. The particle swarm optimization algorithm is employed to search the structure and to determine the ground-state structure of the neutral Zn clusters. The structural motifs are optimized using the density functional theory approach to ensure that the structures are fully relaxed. Results are compared with the literature to validate the accuracy of the prediction method. The binding energy per cluster is obtained and compared with the reported literature to study the stability of these structures. We further assess the electronic properties, including the ionization potential, using the all-electron FHI-aims code employing G0W0 calculations, and the G0W0Г0(1) correction for a few smaller clusters, which provides a better estimation of the ionization potential compared to other methods.

Determination of Ground State Structures of Snx - (x=21-35) Clusters

In this work, an unbiased global search with a homemade genetic algorithm was performed to investigate the structural evolution and electronic properties of Snx - (x=21-35) clusters with density functional theory (DFT) calculations. All the ground-state structures for all these Snx - (x=21-35) clusters have been confirmed by the comparison of the experimental and simulated photoelectron spectra (PESs). It has been revealed that all Snx - (x=21-35) clusters are tricapped trigonal prism (TTP)-based structures consisting of two (for sizes x=21-28) or three (for x=29-35) TTP units, with the remaining atoms adsorbed on the surface or inserted between TTP units. The gradually decreasing HOMO-LUMO gaps indicate that these clusters are undergoing semiconductor-to-metal transformation. The average binding energies show that the structural stabilities of Snx - clusters are not as good as that of silicon and germanium clusters. It found that sizes x=23, 25, 29, 33 show high relative stability.

Boron-based Pd3B26 alloy cluster as a nanoscale antifriction bearing system: tubular core-shell structure, double π/σ aromaticity, and dynamic structural fluxionality

Boron-based nanoclusters show unique geometric structures, nonclassical chemical bonding, and dynamic structural fluxionality. We report here on the theoretical prediction of a binary Pd3B26 cluster, which is composed of a triangular Pd3 core and a tubular double-ring B26 unit in a coaxial fashion, as identified through global structural searches and electronic structure calculations. Molecular dynamics simulations indicate that in the core-shell alloy cluster, the B26 double-ring unit can rotate freely around its Pd3 core at room temperature and beyond. The intramolecular rotation is virtually barrier free, thus giving rise to an antifriction bearing system (or ball bearing) at the nanoscale. The dimension of the dynamic system is only 0.66 nm. Chemical bonding analysis reveals that Pd3B26 cluster possesses double 14π/14σ aromaticity, following the (4n + 2) Hückel rule. Among 54 pairs of valence electrons in the cluster, the overwhelming majority are spatially isolated from each other and situated on either the B26 tube or the Pd3 core. Only one pair of electrons are primarily responsible for chemical bonding between the tube and the core, which greatly weaken the bonding within the Pd3 core and offers structural flexibility. This is a key mechanism that effectively diminishes the intramolecular rotation barrier and facilitates dynamic structural fluxionality of the system. The current work enriches the field of nanorotors and nanomachines.

A database of low-energy atomically precise nanoclusters

The chemical and structural properties of atomically precise nanoclusters are of great interest in numerous applications, but the structures of the clusters can be computationally expensive to predict. In this work, we present the largest database of cluster structures and properties determined using ab-initio methods to date. We report the methodologies used to discover low-energy clusters as well as the energies, relaxed structures, and physical properties (such as relative stability, HOMO-LUMO gap among others) for 63,015 clusters across 55 elements. We have identified clusters for 593 out of 1595 cluster systems (element-size pairs) explored by literature that have energies lower than those reported in literature by at least 1 meV/atom. We have also identified clusters for 1320 systems for which we were unable to find previous low-energy structures in the literature. Patterns in the data reveal insights into the chemical and structural relationships among the elements at the nanoscale. We describe how the database can be accessed for future studies and the development of nanocluster-based technologies.

Structural evolution and electronic properties of pure and semiconductor atom doped in clusters: Inn - , Inn Si- , and Inn Ge- (n = 3-16)

The bonding and electronic properties of In_n ^- , In_n Si^- , and In_n Ge^- (n = 3-16) clusters have been computationally investigated. An intensive global search for the ground-state structures of these clusters were conducted using the genetic algorithm coupled with density functional theory (DFT). The ground-state structures of these clusters have been identified through the comparison between simulated photoelectron spectra (PES) of the found lowest-energy isomers and the experimentally measured ones. Doping semiconductor atom (Si or Ge) can significantly change the structures of the In clusters in most sizes, and the dopant prefers to be surrounded by In atoms. There are three structural motifs for In_n X^- (X = Si, Ge, n = 3-16), and the transition occurs at sizes n = 5 and 13. All In_n Si^- and In_n Ge^- share the same configurations and similar electronic properties except for n = 8. Among all above studied clusters, In₁₃ ^- stands out with the largest vertical detachment energy (VDE), HOMO-LUMO gap, (E_b ) and second order energy difference Δ² E due to its closed electronic shell of (1S)² (1P)⁶ (1D)¹⁰ (2S)² (1F)¹⁴ (2P)⁶ . Similarly, the neutral In₁₂ X (X = Si, Ge) clusters are also identified as superatoms but with electronic configuration of (1S)² (1P)⁶ (2S)² (1D)¹⁰ (1F)¹⁴ (2P)⁶ .

Low-energy configurations of Pt6Cu6 clusters and their physical-chemical characterization: a high-accuracy DFT study

Based on a combination of many-body potentials, an analysis of the inertia tensors and a Density Functional Theory framework, we use a method to harvest the lowest energy states of any set of cluster systems. Then, this methodology is applied to the Pt₆Cu₆ cluster case and the structural, chemical, electronic, anisotropy, magnetic and vibrational properties of the lowest energy isomers are studied. Unexpectedly, some tens of isomers with much lower energy than the precedent believed ground state [J. Chem. Phys., 131(4):044701] are found, which indicates the goodness of this methodology. Some of the isomers obtained present the point groups C_s, C_2v according to Schoenflies notation, while others do not exhibit specific symmetry operations. The global chemical descriptors as the ionization potential, the electron affinity and the chemical hardness have oscillating behaviors with overall decreasing trends as the energy of the isomer grows up, indicating a higher rate of deactivation by sintering processes and a higher strength of the adsorption of small molecules on these systems. We present interesting results of the electronic, magnetic, anisotropy, vibrational and thermal properties of these clusters and discuss them; what can be useful information for future experiments and technical applications in varied fields as catalysis, spintronics, molecular magnetism or magnetic storage information.

Incorporating charge transfer effects into a metallic empirical potential for accurate structure determination in (ZnMg)N nanoalloys

We report the results of a combined empirical potential-density functional theory (EP-DFT) study to assess the global minimum structures of free-standing zinc-magnesium nanoalloys of equiatomic composition and with up to 50 atoms. Within this approach, the approximate potential energy surface generated by an empirical potential is first sampled with unbiased basin hopping simulations, and then a selection of the isomers so identified is re-optimized at a first-principles DFT level. Bader charges calculated in a previous work [A. Lebon, A. Aguado and A. Vega, Corros. Sci., 2017, 124, 35-45] revealed a significant transfer of electrons from Mg to Zn atoms in these nanoalloys; so the main novelty in the present work is the development of an improved EP, termed Coulomb-corrected-Gupta potential, which incorporates an explicit charge-transfer correction term onto a metallic Gupta potential description. The Coulomb correction has a many-body character and is fed with parameterized values of the ab initio Bader charges. The potentials are fitted to a large training set containing DFT values of cluster energies and atomic forces, and the DFT results are used as benchmark data to assess the performance of Gupta and Coulomb-corrected-Gupta EP models. Quite surprisingly, the charge-transfer correction is found to represent only 6% of the nanoalloy binding energies, yet this quantitatively small correction has a sizable beneficial effect on the predicted relative energies of homotops. Zn-Mg bulk alloys are used as the sacrificial material in corrosion-protective coatings, and the long-term goal of our research is to disclose whether those corrosion-protected capabilities are enhanced at the nanoscale.

Sphere Encapsulated Monte Carlo: Obtaining Minimum Energy Configurations of Large Aromatic Systems

We introduce a simple global optimization approach that is able to find minimum energy configurations of clusters containing aromatic molecules. The translational and rotational perturbations required in Monte Carlo-based methods often lead to unrealistic configurations within which two or more molecular rings intersect, causing many of the computational steps to be rejected and the optimization process to be inefficient. Here we develop a modification of the basin-hopping global optimization procedure tailored to tackle problems with intersecting molecular rings. Termed the Sphere Encapsulated Monte Carlo (SEMC) method, this method introduces sphere-based rearrangement and minimization steps at each iteration, and its performance is shown through the exploration of potential energy landscapes of polycyclic aromatic hydrocarbon (PAH) clusters, systems of interest in combustion and astrophysics research. The SEMC method provides clusters that are accurate to 5% mean difference of the minimum energy at a 10-fold speed up compared to previous work using advanced molecular dynamics simulations. Importantly, the SEMC method captures key structural characteristics and molecular size partitioning trends as measured by the molecular radial distances and coordination numbers. The advantages of the SEMC method are further highlighted in its application to previously unstudied heterogeneous PAH clusters.

Computational characterisation of structure and metallicity in small neutral and singly-charged cadmium clusters

Putative global minimum structures for neutral Cd_N and singly charged Cd_N⁺ and Cd_N^- clusters in the small size regime up to N = 21 atoms are reported. A global optimization approach based on the basin hopping method and a Gupta potential fitted to cluster properties is employed to generate a diverse databank of trial structures, which are then re-optimized at the density-functional level of theory. Novel, previously unreported, structures are found for many sizes. Our results successfully reproduce and interpret the size-dependent stabilities known from mass spectrometry, and strongly suggest that experiments aimed at determining the relative stabilities of neutral cadmium clusters are really measuring cation stabilities. We provide an in-depth analysis of electronic structure and use it to explain the gradual emergence of metallic-like behaviour as the cluster size increases.