Robust design of optimal solvents for chemical reactions—A combined experimental and computational strategy

Computer-Aided Molecular Design of Ionic Liquids as Advanced Process Media: A Review from Fundamentals to Applications

The unique physicochemical properties, flexible structural tunability, and giant chemical space of ionic liquids (ILs) provide them a great opportunity to match different target properties to work as advanced process media. The crux of the matter is how to efficiently and reliably tailor suitable ILs toward a specific application. In this regard, the computer-aided molecular design (CAMD) approach has been widely adapted to cover this family of high-profile chemicals, that is, to perform computer-aided IL design (CAILD). This review discusses the past developments that have contributed to the state-of-the-art of CAILD and provides a perspective about how future works could pursue the acceleration of the practical application of ILs. In a broad context of CAILD, key aspects related to the forward structure-property modeling and reverse molecular design of ILs are overviewed. For the former forward task, diverse IL molecular representations, modeling algorithms, as well as representative models on physical properties, thermodynamic properties, among others of ILs are introduced. For the latter reverse task, representative works formulating different molecular design scenarios are summarized. Beyond the substantial progress made, some future perspectives to move CAILD a step forward are finally provided.

Energy Applications of Ionic Liquids: Recent Developments and Future Prospects

Ionic liquids (ILs) consisting entirely of ions exhibit many fascinating and tunable properties, making them promising functional materials for a large number of energy-related applications. For example, ILs have been employed as electrolytes for electrochemical energy storage and conversion, as heat transfer fluids and phase-change materials for thermal energy transfer and storage, as solvents and/or catalysts for CO2 capture, CO2 conversion, biomass treatment and biofuel extraction, and as high-energy propellants for aerospace applications. This paper provides an extensive overview on the various energy applications of ILs and offers some thinking and viewpoints on the current challenges and emerging opportunities in each area. The basic fundamentals (structures and properties) of ILs are first introduced. Then, motivations and successful applications of ILs in the energy field are concisely outlined. Later, a detailed review of recent representative works in each area is provided. For each application, the role of ILs and their associated benefits are elaborated. Research trends and insights into the selection of ILs to achieve improved performance are analyzed as well. Challenges and future opportunities are pointed out before the paper is concluded.

Enhanced Sampling Simulation Reveals How Solvent Influences Chirogenesis of the Intra-Molecular Diels-Alder Reaction

The timescale involved in chemical reactions is quite often beyond that of normal molecular dynamics simulations. Here, we combine metadynamics with selective integrated tempering sampling to simulate an intra-molecular Diels-Alder reaction in explicit solvents. Based on a one-dimensional collective variable obtained from harmonic linear discriminant analysis, four chiral isomers of products were observed in the simulation. Analyses of reactive trajectories showed that this reaction follows a concerted mechanism in all four solvents. In addition, the hydrogen bond between the reactant and water solvent plays an important role in the water-accelerated reaction mechanism. The dynamics of chirality formation varies significantly with solvents. The chirality of products forms significantly before the transition state, especially in ionic liquid.

Computational design of heterogeneous catalysts and gas separation materials for advanced chemical processing

Functional materials are widely used in chemical industry in order to reduce the process cost while simultaneously increase the product quality. Considering their significant effects, systematic methods for the optimal selection and design of materials are essential. The conventional synthesis-and-test method for materials development is inefficient and costly. Additionally, the performance of the resulting materials is usually limited by the designer’s expertise. During the past few decades, computational methods have been significantly developed and they now become a very important tool for the optimal design of functional materials for various chemical processes. This article selectively focuses on two important process functional materials, namely heterogeneous catalyst and gas separation agent. Theoretical methods and representative works for computational screening and design of these materials are reviewed.

Machine learning and molecular descriptors enable rational solvent selection in asymmetric catalysis

Rational solvent selection remains a significant challenge in process development. Here we describe a hybrid mechanistic-machine learning approach, geared towards automated process development workflow. A library of 459 solvents was used, for which 12 conventional molecular descriptors, two reaction-specific descriptors, and additional descriptors based on screening charge density, were calculated. Gaussian process surrogate models were trained on experimental data from a Rh(CO)2(acac)/Josiphos catalysed asymmetric hydrogenation of a chiral α-β unsaturated γ-lactam. With two simultaneous objectives - high conversion and high diastereomeric excess - the multi-objective algorithm, trained on the initial dataset of 25 solvents, has identified solvents leading to better reaction outcomes. In addition to being a powerful design of experiments (DoE) methodology, the resulting Gaussian process surrogate model for conversion is, in statistical terms, predictive, with a cross-validation correlation coefficient of 0.84. After identifying promising solvents, the composition of solvent mixtures and optimal reaction temperature were found using a black-box Bayesian optimisation. We then demonstrated the application of a new genetic programming approach to select an appropriate machine learning model for a specific physical system, which should allow the transition of the overall process development workflow into the future robotic laboratories.

Origin of Asymmetric Solvation Effects for Ions in Water and Organic Solvents Investigated Using Molecular Dynamics Simulations: The Swain Acity-Basity Scale Revisited

The asymmetric solvation of ions can be defined as the tendency of a solvent to preferentially solvate anions over cations or cations over anions, at identical ionic charge magnitudes and effective sizes. Taking water as a reference, these effects are quantified experimentally for many solvents by the relative acity (A) and basity (B) parameters of the Swain scale. The goal of the present study is to investigate the asymmetric solvation of ions using molecular dynamics simulations, and to connect the results to this empirical scale. To this purpose, the charging free energies of alkali and halide ions, and of their hypothetical oppositely charged counterparts, are calculated in a variety of solvents. In a first set of calculations, artificial solvent models are considered that present either a charge or a shape asymmetry at the molecular level. The solvation asymmetry, probed by the difference in charging free energy between the two oppositely charged ions, is found to encompass a term quadratic in the ion charge, related to the different solvation structures around the anion and cation, and a term linear in the ion charge, related to the solvation structure around the uncharged ion-sized cavity. For these simple solvent models, the two terms are systematically counteracting each other, and it is argued that only the quadratic term should be retained when comparing the results of simulations involving physical solvents to experimental data. In a second set of calculations, 16 physical solvents are considered. The theoretical estimates for the acity A are found to correlate very well with the Swain parameters, whereas the correlation for B is very poor. Based on this observation, the Swain scale is reformulated into a new scale involving an asymmetry parameter Σ, positive for acitic solvents and negative for basitic ones, and a polarity parameter Π. This revised scale has the same predictive power as the original scale, but it characterizes asymmetry in an absolute sense, the atomistic simulations playing the role of an extra-thermodynamic assumption, and is optimally compatible with the simulation results. Considering the 55 solvents in the Swain set, it is observed that a moderate basity (Σ between -0.9 and -0.3, related to electronic polarization) represents the baseline for most solvents, while a highly variable acity (Σ between 0.0 and 3.0, related to hydrogen-bond donor capacity modulated by inductive effects) represents a landmark of protic solvents.