Computational Evaluation of Mixtures of Hydrofluorocarbons and Deep Eutectic Solvents for Absorption Refrigeration Systems.

We used computational tools to evaluate three working fluid mixtures for single-effect absorption refrigeration systems, where the generator (desorber) is powered by waste or solar heat. The mixtures studied here resulted from combining a widely used hydrofluorocarbon (HFC) refrigerant, R134a, with three common deep eutectic solvents (DESs) formed by mixing choline chloride (hydrogen bond acceptor, HBA) with urea, glycerol, or ethylene glycol as the hydrogen bond donor (HBD) species. The COSMOtherm/TmoleX software package was used in combination with refrigerant data from NIST/REFPROP, to perform a thermodynamic evaluation of absorption refrigeration cycles using the proposed working fluid mixtures. Afterward, classical MD simulations of the three mixtures were performed to gain insight on these systems at the molecular level. Larger cycle efficiencies are obtained when R134a is combined with choline chloride and ethylene glycol, followed by the system where glycerol is the HBD, and finally that where the HBD is urea. MD simulations indicate that the local density profiles of all species exhibit very sharp variations in systems containing glycerol or urea; furthermore, the Henry's law constants of R134a in these two systems are larger than those observed for the HFC in choline chloride and ethylene glycol, indicating that R134a is more soluble in the latter DES. Interaction energies indicate that the R134a-R134a interactions are weaker in the system where ethylene glycol is the HBD, as compared to in the other DES. Radial distribution functions confirm that in all systems, the DES species do not form strong directional interactions (e.g., hydrogen bonds) with the R134a molecules. Relatively strong interactions are observed between the Cl anions and the hydrogen atoms in R134a; however, the atom-atom interactions between R134a and the cation and HBD species are weaker and do not play a significant role in the solvation of the refrigerant. In all systems, R134a has the largest diffusion coefficients, followed by the HBD, the anion and the cation; the diffusion coefficients are the largest in the systems containing ethylene glycol, followed by those having glycerol and urea. This work is our first step toward our long-term goal of designing and demonstrating optimal working fluid mixtures for use in absorption refrigeration systems. Our results suggest that COSMO-RS can be used to perform a rapid screening of a large number of working fluid mixtures, and select a few candidates for further exploration using molecular simulations and experiments. These latter approaches can be used to refine the accuracy of the COSMO-RS predictions, and to optimize the selection of optimal working fluid mixtures for demonstration in absorption refrigeration systems powered by solar or waste heat sources.

[1]  F. Mjalli,et al.  Prediction of the surface tension of deep eutectic solvents , 2012 .

[2]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[3]  K. R. Seddon,et al.  Applications of ionic liquids in the chemical industry. , 2008, Chemical Society reviews.

[4]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .

[5]  Dmitrij Rappoport,et al.  Property-optimized gaussian basis sets for molecular response calculations. , 2010, The Journal of chemical physics.

[6]  David L Davies,et al.  Novel solvent properties of choline chloride/urea mixtures. , 2003, Chemical communications.

[7]  S. Pandey,et al.  Densities and Viscosities of (Choline Chloride + Urea) Deep Eutectic Solvent and Its Aqueous Mixtures in the Temperature Range 293.15 K to 363.15 K , 2014 .

[8]  S. Pandey,et al.  Densities and dynamic viscosities of (choline chloride + glycerol) deep eutectic solvent and its aqueous mixtures in the temperature range (283.15–363.15) K , 2014 .

[9]  Andrei G. Fedorov,et al.  Absorption Heat Pump/Refrigeration System Utilizing Ionic Liquid and Hydrofluorocarbon Refrigerants , 2012 .

[10]  Haifeng Dong,et al.  Effect of Water on the Density, Viscosity, and CO2 Solubility in Choline Chloride/Urea , 2014 .

[11]  A. Yokozeki,et al.  Vapor-Liquid-Liquid Equilibria of Hydrofluorocarbons + 1-Butyl-3-methylimidazolium Hexafluorophosphate , 2006 .

[12]  J. Perdew,et al.  Density-functional approximation for the correlation energy of the inhomogeneous electron gas. , 1986, Physical review. B, Condensed matter.

[13]  Marek Sierka,et al.  Fast evaluation of the Coulomb potential for electron densities using multipole accelerated resolution of identity approximation , 2003 .

[14]  D. Zheng,et al.  Working Pair Selection of Compression and Absorption Hybrid Cycles through Predicting the Activity Coefficients of Hydrofluorocarbon + Ionic Liquid Systems by the UNIFAC Model , 2012 .

[15]  A. Yokozeki,et al.  Solubility Differences of Halocarbon Isomers in Ionic Liquid [emim][Tf2N] , 2007 .

[16]  A. Klamt,et al.  Fast solvent screening via quantum chemistry: COSMO‐RS approach , 2002 .

[17]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .

[18]  A. W. Schüttelkopf,et al.  PRODRG: a tool for high-throughput crystallography of protein-ligand complexes. , 2004, Acta crystallographica. Section D, Biological crystallography.

[19]  Yan-Fei Shen,et al.  Structural and Dynamical Properties of a Deep Eutectic Solvent Confined Inside a Slit Pore , 2015 .

[20]  M. C. Kroon,et al.  Low-transition-temperature mixtures (LTTMs): a new generation of designer solvents. , 2013, Angewandte Chemie.

[21]  F. Weigend Accurate Coulomb-fitting basis sets for H to Rn. , 2006, Physical chemistry chemical physics : PCCP.

[22]  Stanley I. Sandler,et al.  Improvements of COSMO-SAC for vapor–liquid and liquid–liquid equilibrium predictions , 2010 .

[23]  Filipp Furche,et al.  Nuclear second analytical derivative calculations using auxiliary basis set expansions , 2004 .

[24]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[25]  K. Valsaraj,et al.  Green Leaf Volatiles on Atmospheric Air/Water Interfaces: A Combined Experimental and Molecular Simulation Study , 2014 .

[26]  S. Chungpaibulpatana,et al.  A review of absorption refrigeration technologies , 2001 .

[27]  José Mario Martínez,et al.  PACKMOL: A package for building initial configurations for molecular dynamics simulations , 2009, J. Comput. Chem..

[28]  A. Klamt,et al.  Refinement and Parametrization of COSMO-RS , 1998 .

[29]  Meng-Hui Li,et al.  Molar heat capacities of choline chloride-based deep eutectic solvents and their binary mixtures with water , 2012 .

[30]  Andrei G. Fedorov,et al.  Thermodynamic analysis of an absorption refrigeration system with ionic-liquid/refrigerant mixture as a working fluid , 2012 .

[31]  J. Bara,et al.  COSMOTherm as a Tool for Estimating the Thermophysical Properties of Alkylimidazoles as Solvents for CO2 Separations , 2013 .

[32]  D. Zheng,et al.  A review of imidazolium ionic liquids research and development towards working pair of absorption cycle , 2014 .

[33]  Lynn F. Gladden,et al.  Glycerol eutectics as sustainable solvent systems , 2010 .

[34]  Coray M. Colina,et al.  Experimental and Computational Studies of Choline Chloride-Based Deep Eutectic Solvents , 2014 .

[35]  Paul Painter,et al.  Molecular dynamic simulations and vibrational analysis of an ionic liquid analogue. , 2013, The journal of physical chemistry. B.

[36]  A. Klamt Conductor-like Screening Model for Real Solvents: A New Approach to the Quantitative Calculation of Solvation Phenomena , 1995 .

[37]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[38]  L. Gladden,et al.  Molecular motion and ion diffusion in choline chloride based deep eutectic solvents studied by 1H pulsed field gradient NMR spectroscopy. , 2011, Physical chemistry chemical physics : PCCP.

[39]  Ganesh Kamath,et al.  All-atom force field for the prediction of vapor-liquid equilibria and interfacial properties of HFA134a. , 2009, The journal of physical chemistry. B.

[40]  J. Bara,et al.  Free Volume as the Basis of Gas Solubility and Selectivity in Imidazolium-Based Ionic Liquids , 2012 .

[41]  Mark B. Shiflett,et al.  Solubility and diffusivity of hydrofluorocarbons in room-temperature ionic liquids , 2006 .

[42]  Berk Hess,et al.  P-LINCS:  A Parallel Linear Constraint Solver for Molecular Simulation. , 2008, Journal of chemical theory and computation.

[43]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[44]  Arunprakash T. Karunanithi,et al.  A systematic screening methodology towards exploration of ionic liquids for CO 2 capture processes , 2016 .

[45]  F. Weigend,et al.  Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. , 2005, Physical chemistry chemical physics : PCCP.

[46]  François Jérôme,et al.  Deep eutectic solvents: syntheses, properties and applications. , 2012, Chemical Society reviews.

[47]  A. Klamt,et al.  COSMO-RS as a tool for property prediction of IL mixtures—A review , 2010 .

[48]  P. Kohl,et al.  Performance Simulation of Ionic Liquid and Hydrofluorocarbon Working Fluids for an Absorption Refrigeration System , 2013 .

[49]  G. S. Larsen,et al.  Molecular Simulations of PIM-1-like Polymers of Intrinsic Microporosity , 2011 .

[50]  F. Mjalli,et al.  Effect of water on the thermo-physical properties of Reline: An experimental and molecular simulation based approach. , 2014, Physical chemistry chemical physics : PCCP.

[51]  David Shan-Hill Wong,et al.  Densities of a deep eutectic solvent based on choline chloride and glycerol and its aqueous mixtures at elevated pressures , 2012 .

[52]  M. Parrinello,et al.  Canonical sampling through velocity rescaling. , 2007, The Journal of chemical physics.

[53]  Yong Tae Kang,et al.  Review of advanced absorption cycles: performance improvement and temperature lift enhancement , 2000 .

[54]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[55]  D. Zheng,et al.  Vapor–Liquid Equilibrium Measurements of Difluoromethane + [Emim]OTf, Difluoromethane + [Bmim]OTf, Difluoroethane + [Emim]OTf, and Difluoroethane + [Bmim]OTf Systems , 2011 .

[56]  F. Mjalli,et al.  Acoustic investigation of choline chloride based ionic liquids analogs , 2014 .

[57]  C. A. Infante Ferreira,et al.  Solar refrigeration options – a state-of-the-art review , 2008 .

[58]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[59]  Wenchuan Wang,et al.  Screening of ionic liquids to capture CO2 by COSMO-RS and experiments , 2008 .

[60]  Michael R. Shirts,et al.  Comparison of efficiency and bias of free energies computed by exponential averaging, the Bennett acceptance ratio, and thermodynamic integration. , 2005, The Journal of chemical physics.

[61]  A. Schäfer,et al.  Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr , 1994 .