A graph-based short-cut to low-energy structures

On the brink of self-hydration: the water heptadecamer

For pure, neutral, isolated molecular clusters, (H2O)17 marks the transition from structures with all water molecules on the cluster surface to water self-hydration, i.e., cluster structures around one central water molecule. Getting this right with water model potentials turns out to be challenging. Even the best water potentials currently available, which reproduce collective properties very well, still deliver contradicting results for (H2O)17, when different low-energy isomers from global structure optimizations are examined. Interestingly, ab initio quantum chemistry also struggles with the only seemingly simple question if (H2O)17 is all-surface or water-centered. Hence, although the long history of water potential development may be entering its final phase, it is not quite finished yet.

Globally optimal catalytic fields for a Diels-Alder reaction

In a previous paper [M. Dittner and B. Hartke, J. Chem. Theory Comput. 14, 3547 (2018)], we introduced a preliminary version of our GOCAT (globally optimal catalyst) concept in which electrostatic catalysts are designed for arbitrary reactions by global optimization of distributed point charges that surround the reaction. In this first version, a pre-defined reaction path was kept fixed. This unrealistic assumption allowed for only small catalytic effects. In the present work, we extend our GOCAT framework by a sophisticated and robust on-the-fly reaction path optimization, plus further concomitant algorithm adaptions. This allows smaller and larger excursions from a pre-defined reaction path under the influence of the GOCAT point-charge surrounding, all the way to drastic mechanistic changes. In contrast to the restricted first GOCAT version, this new version is able to address real-life catalysis. We demonstrate this by applying it to the electrostatic catalysis of a prototypical Diels-Alder reaction. Without using any prior information, this procedure re-discovers theoretically and experimentally established features of electrostatic catalysis of this very reaction, including a field-dependent transition from the synchronous, concerted textbook mechanism to a zwitterionic two-step mechanism, and diastereomeric discrimination by suitable electric field components.

Exploring self-organization of molecular tether molecules on a gold surface by global structure optimization

We employ nondeterministic global cluster structure optimization, based on the evolutionary algorithms paradigm, to model the self-assembly of complex molecules on a surface. As a real-life application example directly related to many recent experiments, we use this approach for the assembly of triazatriangulene "platform" molecules on the Au(111) surface. Without additional restrictions like spatial discretizations, coarse-graining or precalculated adsorption poses, and despite the proof-of-principle character of this study, we achieve satisfactory qualitative agreement with several experimental observations and can provide answers to questions that experiments on these species had left open so far. © 2019 Wiley Periodicals, Inc.

Globally Optimal Catalytic Fields - Inverse Design of Abstract Embeddings for Maximum Reaction Rate Acceleration

The search for, and understanding of, good catalysts for chemical reactions is a central issue for chemists. Here, we present first steps toward developing a general computational framework to better support this task. This framework combines efficient, unbiased global optimization techniques with an abstract representation of the catalytic environment, to shrink the search space. To analyze the resulting catalytic embeddings, we employ dimensionality reduction and clustering techniques. This not only provides an inverse design approach to new catalytic embeddings but also illuminates the actual interactions behind catalytic effects. All this is illustrated here with a strictly electrostatic model for the environment and with two versions of a selected example reaction. We close with detailed discussions of future improvements of our framework.

Cluster structures influenced by interaction with a surface

Systematic investigation of how cluster structures change upon interaction with a surface, using global structure optimization by evolutionary algorithms.

Observable-targeting global cluster structure optimization

Global cluster structure optimization not towards lowest energy but towards best fit of cluster properties to experimental data provides theoretical support for cluster experiments under non-equilibrium conditions.