Solvation thermodynamics of organic molecules by the molecular integral equation theory : approaching chemical accuracy

The integral equation theory (IET) of molecular liquids has been an active area of academic research in theoretical and computational physical chemistry for over 40 years because it provides a consistent theoretical framework to describe the structural and thermodynamic properties of liquid-phase solutions. The theory can describe pure and mixed solvent systems (including anisotropic and nonequilibrium systems) and has already been used for theoretical studies of a vast range of problems in chemical physics / physical chemistry, molecular biology, colloids, soft matter, and electrochemistry. A consider- able advantage of IET is that it can be used to study speci fi c solute − solvent interactions, unlike continuum solvent models, but yet it requires considerably less computational expense than explicit solvent simulations

First-principles modeling of chemistry in mixed solvents: Where to go from here?

Mixed solvents (i.e., binary or higher order mixtures of ionic or nonionic liquids) play crucial roles in chemical syntheses, separations, and electrochemical devices because they can be tuned for specific reactions and applications. Apart from fully explicit solvation treatments that can be difficult to parameterize or computationally expensive, there is currently no well-established first-principles regimen for reliably modeling atomic-scale chemistry in mixed solvent environments. We offer our perspective on how this process could be achieved in the near future as mixed solvent systems become more explored using theoretical and computational chemistry. We first outline what makes mixed solvent systems far more complex compared to single-component solvents. An overview of current and promising techniques for modeling mixed solvent environments is provided. We focus on so-called hybrid solvation treatments such as the conductor-like screening model for real solvents and the reference interaction site model, which are far less computationally demanding than explicit simulations. We also propose that cluster-continuum approaches rooted in physically rigorous quasi-chemical theory provide a robust, yet practical, route for studying chemical processes in mixed solvents

All-atom/coarse-grained hybrid predictions of distribution coefficients in SAMPL5

We present blind predictions submitted to the SAMPL5 challenge on calculating distribution coefficients. The predictions were based on estimating the solvation free energies in water and cyclohexane of the 53 compounds in the challenge. These free energies were computed using alchemical free energy simulations based on a hybrid all-atom/coarse-grained model. The compounds were treated with the general Amber force field, whereas the solvent molecules were treated with the Elba coarse-grained model. Considering the simplicity of the solvent model and that we approximate the distribution coefficient with the partition coefficient of the neutral species, the predictions are of good accuracy. The correlation coefficient, R is 0.64, 82% of the predictions have the correct sign and the mean absolute deviation is 1.8 log units. This is on a par with or better than the other simulation-based predictions in the challenge. We present an analysis of the deviations to experiments and compare the predictions to another submission that used all-atom solvent

Predicting solvation free energies using parameter-free solvent models

We present a new approach for predicting solvation free energies in non-aqueous solvents. Utilizing the corresponding states principle, we estimate solvent Lennard-Jones parameters directly from their critical points. Combined with atomic solutes and pressure corrected three-dimensional reference interaction site model (3D-RISM/PC+), the model gives accurate predictions for a wide range of non-polar solvents, including olive oil. The results, obtained without electrostatic interactions and with a very coarse-grained solvent provide an interesting alternative to widely used and heavily parametrized models

Polysaccharides in Solution: Experimental and Computational Studies

Carbohydrates can be found in many natural sources, and they play a central role in various biological processes. These versatile biopolymers are difficult to dissolve in solutions, and therefore converting them into functional forms is a significant challenge, whereby both experimental and computational studies become critical. This chapter discusses commonly used experimental approaches to increase solubility of carbohydrates in solutions and computational studies to rationalize their conformations and solvation effects. Advances of experimental and computational methods for the study of carbohydrates will guide the design approaches to use carbohydrates in industrially and academically relevant applications

Towards Improving The Accuracy of Implicit Solvent Models and Understanding Electrostatic Catalysis in Complex Solvent Environment

This thesis develops improved protocols for studying reactions in solution and uses them to explore the possibility of harnessing complex non-standard solvent environments to catalyse chemical reactions. The thesis covers different but related topics: Improving the accuracy of implicit solvent models. Implicit solvent models are simple cost-effective strategies for modelling solvent as a polarizable continuum. However, the accuracy of this approach can be quite variable. Herein, we examine approaches to improving their accuracy through cavity scaling, the choice of theoretical level and the inclusion of explicit solvent molecules. For SMD, we show that the best performance is achieved when cavity scaling is not employed, while for PCM we present a series of electrostatic scale factors that are radii, solvent and ion type dependent. For both families of method, we also highlight the importance choosing an appropriate level of theory, and identify when explicit solvent molecules are required.. Modelling electrostatic catalysis in complex solvent environment. Recent nanoscale experiments have shown that electric fields are capable of catalysing and controlling chemical reactions, but experimental platforms for scaling these effects remain elusive. Herein, two different approaches to addressing this challenge are explored. The first is using the internal electric field of ordered solvents and ionic liquids, the second is using the electric fields that form naturally at the gas-water interface. A multi-scale modelling approach was developed using polarizable force field based molecular dynamic simulation, post-HF, DFT and semi-empirical quantum chemical calculations. We showed that after exposure to an external electric field, ensembles of solvent or ionic liquid molecules become ordered and this ordering can generate an internal electric field, which persists even after the external potential is removed. Experimental collaborators subsequently detected this field as an open-circuit potential that is strong and long-lived. Computationally we showed that this field is enough to lower reaction barriers by as much as 20 kcal mol-1, and we also developed a predictive structure-reactivity model to choose ionic liquids that optimize this field. In the second approach, we harnessed the electric fields of the gas-water interface. A collaborator showed that in the presence of static, inert gas bubbles, the oxidation potential of HO anion/HO radical was dramatically lowered (by more than 0.5V), much more than any subtle concentration effects predicted by the Nernst equation. Further experiments showed that a high unbalanced concentration of HO- ions (as much as 5M) accumulate at the interface. Our multi-scale modelling calculations showed that this reduction in potential was due to the mutual repulsion of the HO- ions and as little as 1M unbalanced excess was enough to explain the experimental results. The work raises opportunities in reducing the cost of electrochemical processes, and points to electrostatic effects contributing to the well-known catalytic effects of "on water" reactions. Works in this thesis are expected to be useful in the future studies of solution-phase pKa, redox potential, electrostatic catalysis and ionic liquids-based electrochemical devices