Machine Learning for Structure Search of Ligand-protected Nanoclusters

Abstract

Understanding the atomic structures of ligand-protected nanoclusters is essential for their application in various fields. These structures not only determine the physical and chemical properties of ligand-protected nanoclusters but also play a crucial role in their stability and reactivity. Knowing the precise atomic structures allows us to tailor nanoclusters for specific functions. However, because of the extraordinarily high dimensionality of the search space which encompasses an exceptionally large number of all potential structures, it is difficult to use quantum mechanical methods, such as the density functional theory, to find the low-energy structures of ligand-protected nanoclusters. On this point, the structure search of ligand-protected nanoclusters could be more efficient and accurate by utilizing machine learning methods.

In this dissertation, I developed machine learning methods to search the atomic structures of ligand-protected nanoclusters by decomposing the problem into three steps. For the first step, I developed a molecular conformer search procedure based on Bayesian optimization to search the structures of isolated molecules. Using four amino acids as examples, I showed that the procedure is both efficient and accurate. For the second step, I modified the procedure to search the structures of a single ligand on a nanocluster. I also developed and tested strategies to avoid steric clashes between a ligand and cluster parts. Moreover, I tested and demonstrated our modified procedure by searching structures for a cysteine molecule on a well-studied gold-thiolate cluster. As a result, I found that cysteine conformers in a cluster inherit the hydrogen bond types from isolated conformers, while the energy rankings and spacings between the conformers are reordered. In the final step, I applied a machine learning method based on kernel rigid regression (KRR) models to relax the structures of ligand-protected nanoclusters. Moreover, I used an active learning workflow to enhance the relaxation performance of the KRR models. To test and demonstrate our method, I applied it to search structures of Au25(Cys)18 -. We found that the low-energy structures with IItype hydrogen bonds (OH- -N, OH from trans-COOH and N from NH2) are dominant and the different configurations of the ligand layer indeed influence the properties of the clusters.