Computing the melting point and thermodynamic stability of the orthorhombic and monoclinic crystalline polymorphs of the ionic liquid 1-n-butyl-3-methylimidazolium chloride
J. Chem. Phys. 127, 214504 (2007); doi:10.1063/1.2801539
Published 5 December 2007
You are not logged in to this journal. Log in
The melting point, enthalpy of fusion, and thermodynamic stability of two crystal polymorphs of the ionic liquid 1-n-butyl-3-methylimidazolium chloride are calculated using a thermodynamic integration-based atomistic simulation method. The computed melting point of the orthorhombic phase ranges from 365 to 369 K, depending on the classical force field used. This compares reasonably well with the experimental values, which range from 337 to 339 K. The computed enthalpy of fusion ranges from 19 to 29 kJ/mol, compared to the experimental values of 18.5−21.5 kJ/mol. Only one of the two force fields evaluated in this work yielded a stable monoclinic phase, despite the fact that both give accurate liquid state densities. The computed melting point of the monoclinic polymorph was found to be 373 K, which is somewhat higher than the experimental range of 318–340 K. The computed enthalpy of fusion was 23 kJ/mol, which is also higher than the experimental value of 9.3−14.5 kJ/mol. The simulations predict that the monoclinic form is more stable than the orthorhombic form at low temperature, in agreement with one set of experiments but in conflict with another. The difference in free energy between the two polymorphs is very small, due to the fact that a single trans-gauche conformational difference in an alkyl sidechain distinguishes the two structures. As a result, it is very difficult to construct simple classical force fields that are accurate enough to definitively predict which polymorph is most stable. A liquid phase analysis of the probability distribution of the dihedral angles in the alkyl chain indicates that less than half of the dihedral angles are in the gauche-trans configuration that is adopted in the orthorhombic crystal. The low melting point and glass forming tendency of this ionic liquid is likely due to the energy barrier for conversion of the remaining dihedral angles into the gauche-trans state. The simulation procedure used to perform the melting point calculations is an extension of the so-called pseudosupercritical path sampling procedure. This study demonstrates that the method can be effectively applied to quite complex systems such as ionic liquids and that the appropriate choice of tethering potentials for a key step in the thermodynamic path can enable first order phase transitions to be avoided.
©2007 American Institute of Physics
| History: | Received 10 September 2007; accepted 2 October 2007; published 5 December 2007 |
| Permalink: |
http://link.aip.org/link/?JCPSA6/127/214504/1 |
KEYWORDS and PACS
RELATED DATABASES
PUBLICATION DATA
0021-9606 (print)
1089-7690 (online)
REFERENCES (52)
For access to fully linked references, you need to log in.
For access to fully linked references, you need to Log in.
- K. R. Seddon,
Nat. Mater. 2, 363 (2003) . - P. Walden, Bulletin de l'Academie Imperiale des Sciences de St. Petersbourg, 1914, pp. 405–422.
- S. Sugden and H. Wilkins, J. Chem. Soc. 1291 (1929).
- M. Freemantle, Chem. Eng. News 81, 9 (2003).
- J. D. Holbrey and K. R. Seddon, Clean Prod. Process. 1, 223 (1999).
- A. R. Katritzky, A. Lomaka, R. Petrukhin, R. Jain, M. Karelson, A. E. Visser, and R. D. Rogers,
J. Chem. Inf. Comput. Sci. 42, 71 (2002) . - A. R. Katritzky, R. Jain, A. Lomaka, R. Petrukhin, M. Karelson, A. E. Visser, and R. D. Rogers,
J. Chem. Inf. Comput. Sci. 42, 225 (2002) . - D. M. Eike, J. F. Brennecke, and E. J. Maginn,
Green Chem. 5, 323 (2003) . - N. Sun, X. He, K. Dong, X. Zhang, X. Lu, H. He, and S. Zhang,
Fluid Phase Equilib. 246, 137 (2006) . - S. Trohalaki, R. Pachter, G. W. Drake, and T. Hawkins,
Energy Fuels 19, 279 (2005) . - M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Clarendon, Oxford, 1987).
- C. G. Hanke, S. L. Price, and R. M. Lynden-Bell,
Mol. Phys. 99, 801 (2001) . - T. I. Morrow and E. J. Maginn,
J. Phys. Chem. B 106, 12807 (2002) . - C. J. Margulis, H. A. Stern, and B. J. Berne,
J. Phys. Chem. B 106, 12017 (2002) . - J. de Andrade, E. S. Boes, and H. Stassen,
J. Phys. Chem. B 106, 3546 (2002) . - J. N. Canongia Lopes, J. Deschamps, and A. A. H. Pádua,
J. Phys. Chem. B 108, 2038 (2004) . - M. G. Del Pópolo and G. A. Voth,
J. Phys. Chem. B 108, 1744 (2004) . - X. P. Wu, Z. P. Liu, S. P. Huang, and W. C. Wang,
Phys. Chem. Chem. Phys. 7, 2771 (2005) . - J. N. Canongia Lopes, J. Deschamps, and A. A. H. Pádua,
J. Phys. Chem. B 108, 2038 (2004) . - S. Alavi and D. L. Thompson,
J. Phys. Chem. B 109, 18127 (2005) . - C. Cadena and E. J. Maginn, J. Phys. Chem. B 110, 18026 (2006).
- S. Alavi and D. L. Thompson, J. Chem. Phys. 122, 154704 (2005).
- A. B. Belonoshko, R. Ahuja, and B. Johansson, Phys. Rev. B 61, 11928 (2000).
- P. M. Agrawal, B. M. Rice, and D. L. Thompson, J. Chem. Phys. 119, 9617 (2003).
- P. M. Agrawal, B. M. Rice, and D. L. Thompson, J. Chem. Phys. 118, 9680 (2003).
- G. F. Velardez, S. Alavi, and D. L. Thompson, J. Chem. Phys. 119, 6698 (2003).
- J. R. Morris, C. Z. Wang, K. M. Ho, and C. T. Chan, Phys. Rev. B 49, 3109 (1994).
- A. B. Belonoshko and L. S. Dubrovinsky,
Am. Mineral. 81, 303 (1996) . - H. L. Tepper and W. J. Briels, J. Chem. Phys. 115, 9434 (2001).
- B. B. Laird and R. L. Davidchack,
J. Phys. Chem. B 109, 17802 (2005) . - J. Anwar, D. Frenkel, and M. G. Noro, J. Chem. Phys. 118, 728 (2003).
- D. Frenkel and A. J. C. Ladd, J. Chem. Phys. 81, 3188 (1984).
- N. B. Wilding and A. D. Bruce, Phys. Rev. Lett. 85, 5138 (2000).
- D. M. Eike, J. F. Brennecke, and E. J. Maginn, J. Chem. Phys. 122, 014115 (2005).
- D. M. Eike and E. J. Maginn, J. Chem. Phys. 124, 164503 (2006).
- G. Grochola, J. Chem. Phys. 120, 2122 (2004).
- G. Grochola, J. Chem. Phys. 122, 046101 (2005).
- J. D. Holbrey, W. M. Reichert, M. Nieuwenhuyzen, S. Johnston, K. R. Seddon, and R. D. Rogers, Chem. Commun. 2003, 1636.
- S. Saha, S. Hayashi, A. Kobayashi, and H. Hamaguchi,
Chem. Lett. (Jpn.) 32, 740 (2003) . - S. Hayashi, R. Ozawa, and H. Hamaguchi,
Chem. Lett. (Jpn.) 32, 498 (2003) . - K. Nishikawa, S. Wang, H. Katayanagi, S. Hayashi, H. Hamaguchi, Y. Koga, and K. Tozaki,
J. Phys. Chem. B 111, 4894 (2007) . - H. Hamaguchi and R. Ozawa,
Adv. Chem. Phys. 131, 85 (2005) . - See http://therm10.univ-bpclermont.fr/~apadua/research/il.ff, Agilio Padua has made available online the updated ionic liquid force field parameters from his various publications. This website was accessed on February 2007 to obtain updated parameters for the CLP force field.
- K. P. Jensen and W. L. Jorgensen,
J. Chem. Theory Comput. 2, 1499 (2006) . - S. Melchionna, G. Ciccotti, and B. L. Holian,
Mol. Phys. 78, 533 (1993) . - G. J. Martyna, M. E. Tuckerman, D. J. Tobias, and M. L. Klein,
Mol. Phys. 87, 1117 (1996) . - U. Essmann, L. Perera, M. L. Berkowitz, T. Darden, H. Lee, and L. G. Pedersen, J. Chem. Phys. 103, 8577 (1995).
- W. Smith and T. R. Forester,
J. Mol. Graph. 14, 136 (1996) . - J. C. Phillips, R. Braun, W. Wang, J. Gumbart, E. Tajkhorshid, E. Villa, C. Chipot, R. D. Skeel, L. Kale, and K. Schulten,
J. Comput. Chem. 26, 1781 (2005) . - G. J. Martyna, D. J. Tobias, and M. L. Klein, J. Chem. Phys. 101, 4177 (1994).
- A. D. Mackerell, J. Wiorkiewicz-Kuczera, and M. Karplus,
J. Am. Chem. Soc. 117, 11946 (1995) . - M. A. van der Hoef, J. Chem. Phys. 113, 8142 (2000).








