V. Van Speybroeck

A novel force-field parameterization procedure[1] is proposed that surmounts well-known difficulties of the conventional least squares parameterization. The multidimensional ab initio training data are first transformed into individual one-dimensional data sets, each associated with one term in the force-field model. In the second step conventional methods call be used to fit each energy term separately to its corresponding data set. The first step call be completed without any knowledge of the analytical expressions for the energy terms. Moreover the transformed data sets dictate the form of these expressions, which makes the method very suitable for deriving valence force fields. During the transformation in the first step, continuity and least-norm criteria are imposed. The latter facilitate the intuitive physical interpretation of the energy terms that are fitted to the transformed data sets, a prerequisite for transferable force fields. Benchmark parameterizations have been performed oil three small molecules, showing that the new method results in physically intuitive energy terms, exactly when a conventional parameterization would suffer from parameter correlations, i.e. when the number of redundant internal coordinates in the force-field model increases.

While free energy surfaces are the crux of our understanding in many chemical and biological processes, their accuracy is generally unknown. Moreover, many developments to improve their accuracy are often complicated, impeding their general use. Luckily, several tools and guidelines are already in place to identify these shortcomings, but they are typically lacking in flexibility or fail to systematically determine how to improve the accuracy of the free energy calculation. To overcome these limitations, this work introduces OGRe--a python package for optimal grid refinement in an arbitrary number of dimensions. OGRe is based on three metrics which gauge the confinement, consistency, and overlap of each simulation in a series of umbrella sampling (US) simulations, an enhanced sampling technique ubiquitously adopted to construct free energy surfaces for hindered processes. As these three metrics are fundamentally linked to the accuracy of the weighted histogram analysis method, adopted to generate free energy surfaces from US simulations, they facilitate a systematic construction of accurate free energy profiles, where each metric is driven by a specific umbrella parameter. This allows for the derivation of a consistent and optimal collection of umbrellas for each simulation, largely independent of the initial values, thereby dramatically increasing the ease-of-use towards accurate free energy surfaces. As such, OGRe is particularly suited to determined complex free energy surfaces, with large activation barriers and shallow minima, which underpin many physical and chemical transformations, and hence to further our fundamental understanding of these processes.

The electrophilic aromatic substitution of benzenes is part of any undergraduate organic chemistry textbook, yet the mechanism, and more precisely the Wheland intermediate, remains a matter of debate. In this paper, we have computed different reaction paths for the bromination of benzene, anisole and nitrobenzene at the B97X-D/cc-pVTZ level of theory. This revealed, independently of the considered benzenes, a clear kinetic preference for an addition-elimination mechanism, rather than a substitution. Moreover, both mechanisms do not involve a charged Wheland-like intermediate, not in the gas phase nor in the investigated solvents (CCl4 and acetonitrile). Insight into the regioselectivity of the bromination was provided using a combination of conceptual DFT reactivity indices, aromaticity indices, Wiberg bond indices and the non-covalent interaction index. The ortho/para directing effect of the electron-donating methoxy-group in anisole was retrieved and ascribed to a synergy between strong electron delocalisation and attractive interactions. In contrast, the preferred meta-addition on nitrobenzene could not be traced back to any of these effects, nor to the intrinsic reactivity property of the reactant. In this case, an electrostatic clash between the ipso-carbon of the ring and the nitrogen atom resulting from the later nature of the rate-determining step, with respect to anisole, appeared to play a crucial role.