We have developed an algorithm to automatically build the global minimum and other low-energy minima of nanoclusters. This method is implemented in PyAR (https://github.com/anooplab/pyar) program. The global optimization in PyAR involves two parts, generation of several trial geometries and gradient-based local optimization of the trial geometries. While generating the trial geometries, a Tabu list is used for storing the information of the already used trial geometries to avoid using the similar trial geometries. In this recursive algorithm, an n-sized cluster is built from the geometries of n−1 clusters. The overall procedure automatically generates many unique minimum energy geometries of clusters with size from 2 up to n using this evolutionary growth strategy. We have used our strategy on some of the well-studied clusters such as Pd, Pt, Au, and Al homometallic clusters, Ru-Pt and Au-Pt binary clusters, and Ag-Au-Pt ternary cluster. We have analyzed some of the popular parameters to characterize the clusters, such as relative energy, singlet-triplet energy difference, binding energy, second-order energy difference, and mixing energy, and compared with the reported properties.

Isomers of CAl4Mg and CAl4Mg− have been theoretically characterized for the first time. The most stable isomer for both the neutral and anion contain a planar tetracoordinate carbon (ptC) atom. Unlike the isovalent CAl4Be case, which contains a planar pentacoordinate carbon atom as the global minimum geometry, replacing beryllium with magnesium makes the ptC isomer the global minimum due to increased ionic radii of magnesium. However, it is relatively easier to conduct experimental studies for CAl4Mg0/− as beryllium is toxic. While the neutral molecule containing the ptC atom follows the 18 valence electron rule, the anion breaks the rule with 19 valence electrons. The electron affinity of CAl4Mg is in the range of 1.96–2.05 eV. Both the global minima exhibit π/σ double aromaticity. Ab initio molecular dynamics simulations were carried out for both the global minima at 298 K for 10 ps to confirm their kinetic stability.

Density functional theory (DFT) based calculations have been carried out to explore the potential energy surface (PES) of CSinGe4n2+/+/0 (n = 1-3) systems. The global minimum structures in the di-cationic...

We have explored the chemical space of BAl4Mg−/0/+ for the first time and theoretically characterized several isomers with interesting bonding patterns. We have used chemical intuition and a cluster building method based on the tabu-search algorithm implemented in the Python program for aggregation and reaction (PyAR) to obtain the maximum number of possible stationary points. The global minimum geometries for the anion (1a) and cation (1c) contain a planar tetracoordinate boron (ptB) atom, whereas the global minimum geometry for the neutral (1n) exhibits a planar pentacoordinate boron (ppB) atom. The low-lying isomers of the anion (2a) and cation (3c) also contain a ppB atom. The low-lying isomer of the neutral (2n) exhibits a ptB atom. Ab initio molecular dynamics simulations carried out at 298 K for 2000 fs suggest that all isomers are kinetically stable, except the cation 3c. Simulations carried out at low temperatures (100 and 200 K) for 2000 fs predict that even 3c is kinetically stable, which contains a ppB atom. Various bonding analyses (NBO, AdNDP, AIM, etc.) are carried out for these six different geometries of BAl4Mg−/0/+ to understand the bonding patterns. Based on these results, we conclude that ptB/ppB scenarios are prevalent in these systems. Compared to the carbon counter-part, CAl4Mg−, here the anion (BAl4Mg−) obeys the 18 valence electron rule, as B has one electron fewer than C. However, the neutral and cation species break the rule with 17 and 16 valence electrons, respectively. The electron affinity (EA) of BAl4Mg is slightly higher (2.15 eV) than the electron affinity of CAl4Mg (2.05 eV). Based on the EA value, it is believed that these molecules can be identified in the gas phase. All the ptB/ppB isomers exhibit π/σ double aromaticity. Energy decomposition analysis predicts that the interaction between BAl4−/0/+ and Mg is ionic in all these six systems.

We investigated the structural and energetic properties of nanoclusters and nanoalloys composed of group 13 elements (B, Al, and Ga) up to a cluster size of 12. We conducted a...