COSMO-RS-based methods for improved modelling of complex chemical systems

Abstract

The conductor-like screening model for realistic solvation (COSMO-RS) is applied to study three case studies of complex chemical systems. In the first chapter, cellulose solubility in ionic liquid (IL) solvents is considered and a molecular model for cellulose is presented for the calculation of its solubility in ILs and their mixtures with organic molecular solvents. The model is based on cellobiose units which were themselves obtained by a conformation search for a cellotetraose unit. The conformations that make up the model hence can take into account interactions of a dissolved cellulose repeating unit while making sure that the conformations are so that the repeating unit resembles a unit in the middle of a cellulose chain. It is shown that the model contains representations of all of intra-molecular hydrogen bonds as well as open hydrogen bonding sites for accepting inter-molecular ones witnessed experimentally. Relative cellulose solubilities in IL systems with cations [Amim], [Apyr], [C2OHmim], [C2OHpyr], [C2OC1mim], and [C4mim] and anions [Cl], [Br], [N(CN)2], [CH3CHOHCOO], [CH3 COO], [HOCH2 COO], [(C6 H5 )COO], [H2 NCH2 COO], [C2 N3 ], [(C2 H5 )2 PO3 OH] and [HCOO] have been studied. In addition, cellulose solubility in mixtures of [C2mim][CH3COO] and 14 molecular solvents could be reproduced. Finally, a solvent screening over a large database of cation/anion pairs was carried out, providing for a design and discovery guide pertaining to novel IL systems for cellulose dissolution. Despite the complexities of the cellulose chain, particularly regarding its hydrogen bonding sites, the overall structure of its building blocks can be arguably thought of as relatively rigid due to the rigidity of the glucose rings in the cellulose backbone. The flexibility of the hydrogen bonding groups along the chain will be shown to be well captured by a heuristic and systematic conformation search algorithm (COSMOconf). However, problems may arise as the conformational space of the molecule under study increases, so that in practice a systematic approach would not be able to adequately address its relevant conformations. As an alternative method for conformational search, a hybrid modelling strategy for the prediction of octanol/water partition coefficients for alcohol ethoxylate surfactants of varying chain lengths is developed and presented in Chapter 2. The strategy makes use of molecular dynamics (MD) simulations for the generation of molecular conformations in the presence of solvent. A clustering of the conformations from the MD trajectories is then carried out based on principal component analysis of their dihedral angles. Representative conformations thus selected will be then used for COSMO-RS calculations. Each conformation will be then assigned with a weight using an equation derived based on its probability of occurrence in the MD trajectory. Experimental partition coefficients will be thus reproduced within conformation independent accuracy of COSMO-RS. Solvents may also and often are used not only for the sole purpose of dissolving but also as deliberate media in which reactions occur. In Chapter 3, we consider a case for which the significant role of the solvent medium in determining the reaction conversion has been recently attested experimentally. In this chapter, we use COSMO-RS mainly as a molecular model of the liquid phase and study the heterogeneous hydrogenolysis of diphenyl ether (a model compound for lignin itself being a long and complex molecule) on Raney Ni. This is done using the liquid-phase atomistic descriptors COSMO-RS makes available, by means of which a model for the influence of alcoholic solvents on the conversion is then presented. In order to do so, the electrostatics of OH and Ca H bonds was first calculated using in the DFT-COSMO framework. Then, hydrogen bonding donor moment and the screening charge density on the Ca atom of the solvent molecules were used as model variables and were linked to the characteristics of their interactions (adsorption and dissociation) with the catalyst surface according to the available experimental evidence in this regard. In this way, experimentally measured conversions in a series of mono-alcohols for which systematic studies were available could be successfully reproduced. Predictions are additionally done for an extended set of alcohols for which experimental data were not available. The study hence provides for a first step to quantitatively model the dependance of the conversion of this reaction on the choice of solvent using a mechanistically inspired selection of molecular descriptors and is hence a step towards the rational selection and design of solvents with applications in lignin biorefining

    Similar works