Kick-Fukui: A Fukui Function-Guided Method for Molecular Structure Prediction

Here, we introduce a hybrid method, named Kick-Fukui, to explore the potential energy surface (PES) of clusters and molecules using the Coulombic integral between the Fukui functions in the first screening of the best individuals. In the process, small stable molecules or clusters whose combination has the stoichiometry of the explored species are used as assembly units. First, a small set of candidates has been selected from a large and stochastically generated (Kick) population according to the maximum value of the Coulombic integral between the Fukui functions of both fragments. Subsequently, these few candidates are optimized using a gradient method and density functional theory (DFT) calculations. The performance of the program has been evaluated to explore the PES of various systems, including atomic and molecular clusters. In most cases studied, the global minimum (GM) has been identified with a low computational cost. The strategy does not allow to identify the GM of some silicon clusters; however, it predicts local minima very close in energy to the GM that could be used as the initial population of evolutionary algorithms.

Search for Global Minimum Structures of P 2 n + 1 + (n = 1-15) Using xTB-Based Basin-Hopping Algorithm

A new program for searching global minimum structures of atomic clusters using basin-hopping algorithm based on the xTB method was developed here. The program can be performed with a much higher speed than its replacement directly based on DFT methods. Considering the structural varieties and complexities in finding their global minimum structures, phosphorus cluster cations were studied by the program. The global minimum structures of cationic P2n+1+ (n = 1–15) clusters are determined through the unbiased structure searching method. In the last step, further DFT optimization was performed for the selected isomers. For P2n+1+ (n = 1–4, 7), the found global minimum structures are in consistent with the ones previously reported; while for P2n+1+ (n = 5, 6, 8–12), newly found isomers are more energy-favorable than those previously reported. And those for P2n+1+ (n = 13–15) are reported here for the first time. Among them, the most stable isomers of P2n+1+ (n = 4–6, 9) are characterized by their C_3v, C_s, C_2v and C_s symmetry, in turn. But those of P2n+1+ (n = 7, 8, 10–12), no symmetry has been identified. The most stable isomers of P29+ and P31+ are characterized by single P-P bonds bridging units inside the clusters. Further analysis shows that the pnicogen bonds play an important role in the stabilization of these clusters. These results show that the new developed program is effective and robust in searching global minimum structures for atom clusters, and it also provides new insights into the role of pnicogen bonds in phosphorus clusters.

Application of a quantum genetic algorithm and QTAIM analysis in the study of structural and electronic properties of neutral bimetallic clusters NaxLiy (4 ≤ x + y ≤ 10)

Alloy clusters of Na_xLi_y (4 ≤ x + y ≤ 10) are studied by exploring the potential energy surface in the ab initio MP2 level with the support of a quantum genetic algorithm (QGA). In some cases, the structures have been also refined with DFT and coupled-cluster methods. The general trends of sodium-lithium structures are in line with previous studies. The ionization potentials and polarizabilities to all structures were calculated with MP2 method and the average error between these two properties compared with experimental data was 6% and 13%, respectively. The topological analysis based on quantum theory of atoms in molecules (QTAIM) showed that by increasing the cluster size of the diatomic system there was a decrease of atomic interaction energies. The degree of degeneracy from D3BIA aromaticity index and the analysis of the atomic charges showed the influence (by charge transfer) of the chemical element in lower quantity in the cluster with respect to the other atoms. Our achievements of comparing our theoretical results with available experimental data have demonstrated that our approach can also predict satisfactorily quantum atomic and alloy clusters properties, at least, for low nuclearities.

NWPEsSe: An Adaptive-Learning Global Optimization Algorithm for Nanosized Cluster Systems

Global optimization constitutes an important and fundamental problem in theoretical studies in many chemical fields, such as catalysis, materials, or separations problems. In this paper, a novel algorithm has been developed for the global optimization of large systems including neat and ligated clusters in the gas phase and supported clusters in periodic boundary conditions. The method is based on an updated artificial bee colony (ABC) algorithm method, that allows for adaptive-learning during the search process. The new algorithm is tested against four classes of systems of diverse chemical nature: gas phase Au₅₅, ligated Au₈²⁺, Au₈ supported on graphene oxide and defected rutile, and a large cluster assembly [Co₆Te₈(PEt₃)₆][C₆₀]_n, with sizes ranging between 1 and 3 nm and containing up to 1300 atoms. Reliable global minima (GMs) are obtained for all cases, either confirming published data or reporting new lower energy structures. The algorithm and interface to other codes in the form of an independent program, Northwest Potential Energy Search Engine (NWPEsSe), is freely available, and it provides a powerful and efficient approach for global optimization of nanosized cluster systems.

Iodization Threshold in Size-Dependent Reactions of Lead Clusters Pbn+ with Iodomethane

Utilizing a magnetron-sputtering (MagS) source in tandem with a multiple-ion laminar flow tube (MIFT) reactor and a customized triple quadrupole mass spectrometer (TQMS), we have prepared clean Pb_n⁺ (n = 1-13) clusters and measured their reactivity with iodomethane under high carrier gas pressures. Strong size dependences are found for the reactivity of these cationic Pb_n⁺ clusters with CH₃I. For the Pb_n⁺ with n ≤ 4, iodinated clusters Pb_nI⁺ were found to be the dominant products, in strong contrast to n > 4 where no such products were seen. Quantum chemical studies show that with an increasing number of Pb atoms, the Pb-Pb interatomic interactions become stronger compared with the Pb-I bonding in Pb_nI⁺ clusters. Furthermore, the reactions of Pb_1-4⁺ with CH₃I have fairly small transition state energy barriers, in contrast to those for Pb_n>4⁺ clusters which have magnitudes that will prevent reactions under the ambient conditions.

Exploring the Structure, Energetic, and Magnetic Properties of Neutral Small Lithium Clusters Doped with Yttrium: Supermagnetic Atom Research

Density functional theory calculations based on magnetic and energetic stability criteria were performed to study a series of yttrium-doped lithium neutral clusters. A relativistic approximation was employed to properly describe the energy and multiplicity of the given clusters' fundamental states. The interaction of the 4d-Y atomic orbitals with the sp-Li states had an important role in the magnetic and energetic behavior of the selected systems. The spin density was concentrated over the yttrium atom regardless of the size of the cluster. Li7Y is a new stable superatom due to its enhanced magnetic properties.