Assessing the persistence of chalcogen bonds in solution with neural network potentials

Overcoming the Pitfalls of Computing Reaction Selectivity from Ensembles of Transition States

The prediction of reaction selectivity is a challenging task for computational chemistry, not only because many molecules adopt multiple conformations but also due to the exponential relationship between effective activation energies and rate constants. To account for molecular flexibility, an increasing number of methods exist that generate conformational ensembles of transition state (TS) structures. Typically, these TS ensembles are Boltzmann weighted and used to compute selectivity assuming Curtin-Hammett conditions. This strategy, however, can lead to erroneous predictions if the appropriate filtering of the conformer ensembles is not conducted. Here, we demonstrate how any possible selectivity can be obtained by processing the same sets of TS ensembles for a model reaction. To address the burdensome filtering task in a consistent and automated way, we introduce marc, a tool for the modular analysis of representative conformers that aids in avoiding human errors while minimizing the number of reoptimization computations needed to obtain correct reaction selectivity.

Modelling chemical processes in explicit solvents with machine learning potentials

Solvent effects influence all stages of the chemical processes, modulating the stability of intermediates and transition states, as well as altering reaction rates and product ratios. However, accurately modelling these effects remains challenging. Here, we present a general strategy for generating reactive machine learning potentials to model chemical processes in solution. Our approach combines active learning with descriptor-based selectors and automation, enabling the construction of data-efficient training sets that span the relevant chemical and conformational space. We apply this strategy to investigate a Diels-Alder reaction in water and methanol. The generated machine learning potentials enable us to obtain reaction rates that are in agreement with experimental data and analyse the influence of these solvents on the reaction mechanism. Our strategy offers an efficient approach to the routine modelling of chemical reactions in solution, opening up avenues for studying complex chemical processes in an efficient manner.

From Organic Fragments to Photoswitchable Catalysts: The OFF-ON Structural Repository for Transferable Kernel-Based Potentials

Structurally and conformationally diverse databases are needed to train accurate neural networks or kernel-based potentials capable of exploring the complex free energy landscape of flexible functional organic molecules. Curating such databases for species beyond "simple" drug-like compounds or molecules composed of well-defined building blocks (e.g., peptides) is challenging as it requires thorough chemical space mapping and evaluation of both chemical and conformational diversities. Here, we introduce the OFF-ON (organic fragments from organocatalysts that are non-modular) database, a repository of 7869 equilibrium and 67,457 nonequilibrium geometries of organic compounds and dimers aimed at describing conformationally flexible functional organic molecules, with an emphasis on photoswitchable organocatalysts. The relevance of this database is then demonstrated by training a local kernel regression model on a low-cost semiempirical baseline and comparing it with a PBE0-D3 reference for several known catalysts, notably the free energy surfaces of exemplary photoswitchable organocatalysts. Our results demonstrate that the OFF-ON data set offers reliable predictions for simulating the conformational behavior of virtually any (photoswitchable) organocatalyst or organic compound composed of H, C, N, O, F, and S atoms, thereby opening a computationally feasible route to explore complex free energy surfaces in order to rationalize and predict catalytic behavior.

Halide Complexes of 5,6-Dicyano-2,1,3-Benzoselenadiazole with 1 : 4 Stoichiometry: Cooperativity between Chalcogen and Hydrogen Bonding

The [M4 -Hal]- (M=the title compound; Hal=Cl, Br, and I) complexes were isolated in the form of salts of [Et4 N]+ cation and characterized by XRD, NMR, UV-Vis, DFT, QTAIM, EDD, and EDA. Their stoichiometry is caused by a cooperative interplay of σ-hole-driven chalcogen (ChB) and hydrogen (HB) bondings. In the crystal, [M4 -Hal]- are connected by the π-hole-driven ChB; overall, each [Hal]- is six-coordinated. In the ChB, the electrostatic interaction dominates over orbital and dispersion interactions. In UV-Vis spectra of the M+[Hal]- solutions, ChB-typical and [Hal]- -dependent charge-transfer bands are present; they reflect orbital interactions and allow identification of the individual [Hal]- . However, the structural situation in the solutions is not entirely clear. Particularly, the UV-Vis spectra of the solutions are different from the solid-state spectra of the [Et4 N]+ [M4 -Hal]- ; very tentatively, species in the solutions are assigned [M-Hal]- . It is supposed that the formation of the [M4 -Hal]- proceeds during the crystallization of the [Et4 N]+ [M4 -Hal]- . Overall, M can be considered as a chromogenic receptor and prototype sensor of [Hal]- . The findings are also useful for crystal engineering and supramolecular chemistry.

Chalcogen Bonds: How to Characterize Them in Solution?

Chalcogen bonds (ChBs) occur between molecules containing Lewis acidic chalcogen substituents and Lewis bases. Recently, ChB emerged as a pivotal interaction in solution-based applications such as anion recognition, anion transport and catalysis. However, before moving to applications, the involvement of ChB must be established in solution. In this Concept article, we provide a brief review of the currently available experimental investigations of ChB in solution.

Learning Conductance: Gaussian Process Regression for Molecular Electronics

Experimental studies of charge transport through single molecules often rely on break junction setups, where molecular junctions are repeatedly formed and broken while measuring the conductance, leading to a statistical distribution of conductance values. Modeling this experimental situation and the resulting conductance histograms is challenging for theoretical methods, as computations need to capture structural changes in experiments, including the statistics of junction formation and rupture. This type of extensive structural sampling implies that even when evaluating conductance from computationally efficient electronic structure methods, which typically are of reduced accuracy, the evaluation of conductance histograms is too expensive to be a routine task. Highly accurate quantum transport computations are only computationally feasible for a few selected conformations and thus necessarily ignore the rich conformational space probed in experiments. To overcome these limitations, we investigate the potential of machine learning for modeling conductance histograms, in particular by Gaussian process regression. We show that by selecting specific structural parameters as features, Gaussian process regression can be used to efficiently predict the zero-bias conductance from molecular structures, reducing the computational cost of simulating conductance histograms by an order of magnitude. This enables the efficient calculation of conductance histograms even on the basis of computationally expensive first-principles approaches by effectively reducing the number of necessary charge transport calculations, paving the way toward their routine evaluation.

How Robust Is the Reversible Steric Shielding Strategy for Photoswitchable Organocatalysts?

A highly appealing strategy to modulate a catalyst's activity and/or selectivity in a dynamic and noninvasive way is to incorporate a photoresponsive unit into a catalytically competent molecule. However, the description of the photoinduced conformational or structural changes that alter the catalyst's intrinsic reactivity is often reduced to a handful of intuitive static representations, which can struggle to capture the complexity of flexible organocatalysts. Here, we show how a comprehensive exploration of the free energy landscape of N-alkylated azobenzene-tethered piperidine catalysts is essential to unravel the conformational characteristics of each configurational state and explain the experimentally observed reactivity trends. Mapping the catalysts' conformational space highlights the existence of false ON or OFF states that lower their switching ability. Our findings expose the challenges associated with the realization of a reversible steric shielding for the photocontrol of Brønsted basicity of piperidine photoswitchable organocatalysts.