Abstract:The ability to find optimal molecular structures with desired properties is a popular challenge, with applications in areas such as drug discovery. Genetic algorithms are a common approach to global minima molecular searches due to their ability to search large regions of the energy landscape and decrease computational time via parallelization. In order to decrease the amount of unstable intermediate structures being produced and increase the overall efficiency of an evolutionary algorithm, clustering was introduced in multiple instances. However, there is little literature detailing the effects of differentiating the selection frequencies between clusters. In order to find a balance between exploration and exploitation in our genetic algorithm, we propose a system of clustering the starting population and choosing clusters for an evolutionary algorithm run via a dynamic probability that is dependent on the fitness of molecules generated by each cluster. We define four parameters, MFavOvrAll-A, MFavClus-B, NoNewFavClus-C, and Select-D, that correspond to a reward for producing the best structure overall, a reward for producing the best structure in its own cluster, a penalty for not producing the best structure, and a penalty based on the selection ratio of the cluster, respectively. A reward increases the probability of a cluster's future selection, while a penalty decreases it. In order to optimize these four parameters, we used a Gaussian distribution to approximate the evolutionary algorithm performance of each cluster and performed a grid search for different parameter combinations. Results show parameter MFavOvrAll-A (rewarding clusters for producing the best structure overall) and parameter Select-D (appearance penalty) have a significantly larger effect than parameters MFavClus-B and NoNewFavClus-C. In order to produce the most successful models, a balance between MFavOvrAll-A and Select-D must be made that reflects the exploitation vs exploration trade-off often seen in reinforcement learning algorithms. Results show that our reinforcement-learning-based method for selecting clusters outperforms an unclustered evolutionary algorithm for quinoline-like structure searches.

Discrete Optimal Global Convergence of an Evolutionary Algorithm for Clusters under the Potential of Lennard Jones

Cooperative Bare-Bone Particle Swarm Optimization for Data Clustering

Bisexual Evolution: A Novel Bisexual Evolutionary Framework Based on the Fisher's Runaway Process

An Adaptive Genetic Algorithm for determining optimal structures for atomic clusters

Optimization of molecular clusters configurations using a Genetic Algorithm

An approach based on genetic algorithms and machine learning coupled for studying alloy and molecular clusters by optimizing quantum energy surfaces

Robustness of "cut and splice" genetic algorithms in the structural optimization of atomic clusters

Linked cluster expansion on trees

Generating configurations of increasing lattice size with machine learning and the inverse renormalization group

A Lattice-based Method for Optimization in Continuous Spaces with Genetic Algorithms

Testing Lennard-Jones Clusters for Optimality

Tuning Reinforcement Learning Parameters for Cluster Selection to Enhance Evolutionary Algorithms

A flexible and adaptive grid algorithm for global optimization utilizing basin hopping Monte Carlo

Global Optimization Method for Cluster Structures.

Cartesian atomic cluster expansion for machine learning interatomic potentials

Cluster-MLP: An Active Learning Genetic Algorithm Framework for Accelerated Discovery of Global Minimum Configurations of Pure and Alloyed Nanoclusters

Density-Sensitive Evolutionary Clustering

Optimization of Discrete Parameters Using the Adaptive Gradient Method and Directed Evolution

Unbiased Fuzzy Global Optimization of Lennard-Jones Clusters for N ≤ 1000.

Genetic Algorithms for Interior Comparative Optimization of Standard BCS Parameters in Selected Superconductors and High-Temperature Superconductors

The Compact Genetic Algorithm Struggles on Cliff Functions