Molecular dynamics calculations of N-methylpyrrolidone in liquid water

Abstract

A comprehensive molecular dynamics study has been carried out on the structural and dynamical behaviour of N-methylpyrrolidone (MeP) in water. The impetus for these calculations derives from the growing interest in kinetic hydrate inhibitors, of which polyvinylpyrrolidone (PVP) is the prototype. Study of the smallest unit of this polymer provides valuable insight into the effect of the pyrrolidone ring on water. The present study indicates a strong structuring effect of MeP on surrounding water molecules, with up to four solvation shells being clearly discernible around the solute and markedly slower dynamics for water in the vicinity of MeP. Synergistic solvation effects were found to arise from the conjunction of hydrophobic solvation around the methyl group, and hydrophilic solvation around the amide oxygen; this synergism is likely to be important in understanding the role of pyrrolidone, and other lactam rings, in inhibiting clathrate hydrate formation.

References

1. I. M. Klotz, Science, 1959, 128, 815; H. H. G. Jellinek and S. Y. Fok, Kolloid-Z. Z. Polym., 1967, 220, 122; O. Guven and E. Eltan, Makromol. Chem., 1981, 182, 3129.
2. A. Ben-David and S. Gavendo, Cryobiology, 1972, 9, 192.
3. J. P. Lederhos, J. P. Long, A. Sum, R. L. Christiansen and E. D. Sloan, Chem. Eng. Sci., 1996, 8, 1221.
4. E. D. Sloan, Clathrate Hydrates of Natural Gases, 2nd edn., Marcel Dekker Inc., New York, 1998.
5. F. H. Stillinger, Science, 1980, 209, 451.
6. P. J. Rossky and M. Karplus, J. Am. Chem. Soc., 1979, 101, 1913.
7. M. Ferrario, M. Haughney, I. R. McDonald and M. L. Klein, J. Chem. Phys., 1990, 93, 5156; M. Ferrario, I. R. McDonald and M. C. R. Symons, Mol. Phys., 1992, 77, 617.
8. P.-O. Åstrand, A. Wallqvist, G. Karlstrom and P. Linse, J. Chem. Phys., 1991, 95, 8419.
9. P. M. Rodger, Mol. Phys., 1996, 89, 1157.
10. A. Laaksonen, P. G. Kusalik and I. M. Svishchev, J. Phys. Chem., 1997, 101, 5910.
11. B. Kvamme, G. Huseby and O. K. Forrisdahl, Mol. Phys., 1997, 90, 979.
12. T. J. Carver, M. G. B. Drew and P. M. Rodger, J. Chem. Soc., Faraday Trans., 1995, 91, 3449.
13. H. J. C. Berendsen, J. P. M. Postma, W. F. von Gunsteren and J. Hermens, in Intermolecular Forces, ed. B. Pullman, Reidel, Dordrecht, 1981, pp. 331–342.
14. W. L. Jorgensen, J. D. Madura and C. J. Swensen, J. Am. Chem. Soc., 1984, 106, 6638.
15. DL_POLY is a package of molecular simulations routine written by W. Smith and T. R. Forester, copyright CCLRC, Daresbury Laboratory, Daresbury, Cheshire, 1995.
16. S. Nosé, Mol. Phys., 1984, 52, 255.
17. W. G. Hoover, Phys. Rev., 1985, A31, 1695.
18. J. P. Rykaert, G. Ciccotti and H. J. C. Berendsen, J. Comput. Phys., 1977, 23, 327.
19. T. Astley, G. G. Birch, M. G. B. Drew, P. M. Rodger and G. R. H. Wilden, J. Comput. Chem., 1998, 19, 363.
20. Complete clathrate cages form with 20, 24 or 28 waters.
21. I. I. Vaisman, F. K. Brown and A. Tropsha, J. Phys. Chem., 1994, 98, 5559.
22. K. Nakanishi, K. Ikari, S. Okazaki and H. Touhara, J. Chem. Phys., 1984, 80, 1656.
23. I. M. Svishchev and P. G. Kusalik, J. Chem. Phys., 1993, 99, 3049.
24. I. M. Svishchev and P. G. Kusalik, Science, 1994, 265, 1219.
25. The local χ direction is defined with opposite senses for O_M and Me, so that in each case φ= 270° points towards the other group.
26. A. Rahman and F. H. Stillinger, J. Am. Chem. Soc., 1973, 95, 7943.
27. W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey and M. L. Klein, J. Chem. Phys., 1983, 79, 926.
28. P. J. Rossky and M. Karplus, J. Am. Chem. Soc., 1979, 101, 1913.
29. L. Ambrosone, G. D'Errico, R. Sartorio and V. Vitaliano, J. Chem. Soc., Faraday Trans., 1995, 91, 1339.
30. J. Anderson, J. J. Ullo and S. Yip, J. Chem. Phys., 1987, 87, 1726.