• Volume/Page
  • Keyword
  • DOI
  • Citation
  • Advanced
   
 
 
 

You Tube Flickr Twitter UniPHY Group iResearch App Facebook

Author Select

FULL-TEXT OPTIONS:

J. Chem. Phys. 136, 014511 (2012); http://dx.doi.org/10.1063/1.3671997 (9 pages)

Free energy profiles for penetration of methane and water molecules into spherical sodium dodecyl sulfate micelles obtained using the thermodynamic integration method combined with molecular dynamics calculations

K. Fujimoto1, N. Yoshii1,2, and S. Okazaki1

1Department of Applied Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
2Center for Computational Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan

View MapView Map

(Received 14 October 2011; accepted 5 December 2011; published online 6 January 2012)

The free energy profiles, ΔG(r), for penetration of methane and water molecules into sodium dodecyl sulfate (SDS) micelles have been calculated as a function of distance r from the SDS micelle to the methane and water molecules, using the thermodynamic integration method combined with molecular dynamics calculations. The calculations showed that methane is about 6–12 kJ mol−1 more stable in the SDS micelle than in the water phase, and no ΔG(r) barrier is observed in the vicinity of the sulfate ions of the SDS micelle, implying that methane is easily drawn into the SDS micelle. Based on analysis of the contributions from hydrophobic groups, sulfate ions, sodium ions, and solvent water to ΔG(r), it is clear that methane in the SDS micelle is about 25 kJ mol−1 more stable than it is in the water phase because of the contribution from the solvent water itself. This can be understood by the hydrophobic effect. In contrast, methane is destabilized by 5–15 kJ mol−1 by the contribution from the hydrophobic groups of the SDS micelle because of the repulsive interactions between the methane and the crowded hydrophobic groups of the SDS. The large stabilizing effect of the solvent water is higher than the repulsion by the hydrophobic groups, driving methane to become solubilized into the SDS micelle. A good correlation was found between the distribution of cavities and the distribution of methane molecules in the micelle. The methane may move about in the SDS micelle by diffusing between cavities. In contrast, with respect to the water, ΔG(r) has a large positive value of 24–35 kJ mol−1, so water is not stabilized in the micelle. Analysis showed that the contributions change in complex ways as a function of r and cancel each other out. Reference calculations of the mean forces on a penetrating water molecule into a dodecane droplet clearly showed the same free energy behavior. The common feature is that water is less stable in the hydrophobic core than in the water phase because of the energetic disadvantage of breaking hydrogen bonds formed in the water phase. The difference between the behaviors of the SDS micelles and the dodecane droplets is found just at the interface; this is caused by the strong surface dipole moment formed by sulfate ions and sodium ions in the SDS micelles.

© 2012 American Institute of Physics

Article Outline

  1. INTRODUCTION
  2. METHOD
    1. Free energy profile
    2. MD calculations
    3. SDS micelle
    4. Dodecane droplet
  3. RESULTS AND DISCUSSION
    1. Mean force
    2. Free energy profile for penetration
    3. Contributions from hydrophobic groups, sulfate ions, sodium ions, and solvent water to Δ G ( r )
      1. Methane
      2. Water
    4. Methane distribution in the SDS micelle solution
  4. CONCLUSION

KEYWORDS, PACS, and IPC

PACS

  • 31.70.Dk

    Environmental and solvent effects

  • 31.15.xv

    Molecular dynamics and other numerical methods

  • 33.15.Kr

    Electric and magnetic moments (and derivatives), polarizability, and magnetic susceptibility

  • 33.15.Fm

    Bond strengths, dissociation energies

  • 82.70.Dd

    Colloids

  • 82.70.Uv

    Surfactants, micellar solutions, vesicles, lamellae, amphiphilic systems, (hydrophilic and hydrophobic interactions)

International Patent Classification (IPC)

  • B01J13/00

    Colloid chemistry, e.g. the production of colloidal materials or their solutions, not otherwise provided for; Making microcapsules or microballoons

ARTICLE DATA

PUBLICATION DATA

ISSN

0021-9606 (print)  
1089-7690 (online)

  1. H. Nakajima , in Industrial Applications of Microemulsions, Surfactant Science Series Vol. 66, edited by C. Solans and H. Kunieda (Marcel Dekker, New York, 1997), p. 175.
  2. Md. H. Uddin, N. Kanei, and H. Kunieda, Langmuir 16, 6891 (2000).
  3. Y. Tokuoka, H. Uchiyama, M. Abe, and K. Ogino, J. Colloid Interface Sci. 152, 402 (1992).
  4. Y. Tokuoka, H. Uchiyama, and M. Abe, J. Phys. Chem. 98, 6167 (1992).
  5. Y. Tokuoka, H. Uchiyama, M. Abe, and S. D. Christian, Langmuir, 11, 725 (1995).
  6. M. J. Garcia–Celma , in Industrial Applications of Microemulsions, Surfactant Science Series Vol. 66, edited by C. Solans and H. Kunieda (Marcel Dekker, New York, 1997), p. 123.
  7. M. Yokoyama, M. Miyauchi, N. Yamada, T. Okano, Y. Sakurai, K. Kataoka, and S. Inoue, Cancer Res. 50, 1693 (1990).
  8. K. Kataoka, A. Harada, and Y. Nagasaki, Adv. Drug Delivery Rev. 47, 113 (2001). [MEDLINE]
  9. H. Kuramochi, Y. Andoh, N. Yoshii, and S. Okazaki, J. Phys. Chem. B 113, 15181 (2009).
  10. D. L. Foster and R. H. Fillingame, J. Biol. Chem. 254, 8230 (1979). [MEDLINE]
  11. M. J. Newman, D. L. Foster, T. H. Wilson, and H. R. Kaback, J. Biol. Chem. 256, 11804 (1981). [MEDLINE]
  12. H. Itami, Y. Sakaki, T. Shimamoto, H. Hama, M. Tsuda, and T. Tsuchiya, J. Biochem. (Tokyo) 105, 785 (1989).
  13. V. Pillai and D. O. Shah , in Industrial Applications of Microemulsions, Surfactant Science Series Vol. 66, edited by C. Solans and H. Kunieda (Marcel Dekker, New York, 1997), p. 227.
  14. M. A. López–Quintela, J. Quibén–Solla, and J. Rivas , in Industrial Applications of Microemulsions, Surfactant Science Series Vol. 66, edited by C. Solans and H. Kunieda (Marcel Dekker, New York, 1997), p. 247.
  15. Y. Moroi , Micelles Theoretical and Applied Aspects (Plenum, New York, 1992).
  16. N. Matubayasi, K. K. Liang, and M. Nakahara, J. Chem. Phys. 124, 154908 (2006)JCPSA6000124000015154908000001. [MEDLINE]
  17. A. Wishnia, J. Phys. Chem. 67, 2079 (1963). [CAS]
  18. B. Han, G. Yang, H. Yan, Z. Li, R. K. Thomas, J. Penfold, R. K. Heenan, and D. C. Steytler, J. Colloid Interface Sci. 218, 145 (1999). [MEDLINE] [CAS]
  19. M. Hai, B. Han, G. Yang, H. Yan, and Q. Han, Langmuir 15, 1640 (1999). [ISI] [CAS]
  20. P. Mukerjee and J. R. Cardinal, J. Phys. Chem. 82, 16414 (1978). [CAS]
  21. J. R. Cardinal and P. Mukerjee J. Phys. Chem. 82, 16420 (1978).
  22. C. Ramachandran, R. A. Pyter, and P. Mukerjee, J. Phys. Chem. 86, 3206 (1982).
  23. W. D. Harkins, R. W. Matton, and R. Mittelmann, J. Chem. Phys. 15, 763 (1947)JCPSA6000015000010000763000001. [ISI] [CAS]
  24. W. D. Harkins, R. Mittelmann, and M. L. Corrin, J. Phys Chem. 53, 1350 (1947).
  25. N. Yoshii and S. Okazaki, J. Chem. Phys. 126, 096101 (2007)JCPSA6000126000009096101000001. [MEDLINE]
  26. K. Fujimoto, N. Yoshii, and S. Okazaki, J. Chem. Phys. 133, 074511 (2010)JCPSA6000133000007074511000001. [MEDLINE]
  27. K. Fujimoto, N. Yoshii, and S. Okazaki , Mol. Sim. (in press).
  28. S. Marrink and H. J. C. Berendsen, J. Phys. Chem. 98, 4155 (1994).
  29. W. Shinoda, M. Mikami, T. Baba, and M. Hato, J. Phys. Chem. B 108, 4155 (2004).
  30. D. Bemporad, J. W. Essex, and C. Luttmann, J. Phys. Chem. B 108, 4875 (2004). [ISI]
  31. N. Yoshii, K. Iwahashi, and S. Okazaki, J. Chem. Phys. 124, 184901 (2006)JCPSA6000124000018184901000001. [MEDLINE]
  32. H. F. Huisman, Proc. K. Ned. Akad. Wet., Ser. B: Phys. Sci. 67, 388 (1964).
  33. S. Hayashi and S. Ikeda, J. Phys. Chem. 84, 744 (1980).
  34. P. Lianos and R. Zana, J. Colloid Interface Sci. 84, 100 (1981). [CAS]
  35. W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys. 79, 926 (1983).
  36. A. D. MacKerell Jr., D. Bashford, M. Bellott, R. L. Dunbrack Jr., J. D. Evanseck, M. J. Field, S. Fischer, J. Gao, H. Guo, S. Ha, D. J. McCarthy, L. Kuchnir, K. Kuczera, F. T. K. Lau, C. Mattos, S. Michnick, T. Ngo, D. T. Nguyen, B. Prodhom, W. E. Reiher III, B. Roux, M. Schlenkrich, J. C. Smith, R. Stote, J. Straub, M. Watanabe, J. Wiórkiewicz–Kuczera, D. Yin, and M. Karplus, J. Phys. Chem. B 102, 3586 (1998).
  37. A. D. MacKerell Jr., M. Feig, and C. L. Brooks III, J. Comput. Chem. 25, 1400 (2004). [MEDLINE]
  38. G. J. Martyna, D. J. Tobias, and M. L. Klein, J. Chem. Phys. 101, 4177 (1994)JCPSA6000101000005004177000001.
  39. G. J. Martyna, M. E. Tuckerman, D. J. Tobias, and M. L. Klein, Mol. Phys. 87, 1117 (1996).
  40. G. J. Martyna, M. L. Klein, and M. E. Tuckerman, J. Chem. Phys. 97, 2635 (1992)JCPSA6000097000004002635000001.
  41. U. Essmann, L. Perera, M. L. Berkowitz, T. A. Darden, H. Lee, and L. Pedersen, J. Chem. Phys. 103, 8577 (1995).
  42. A. D. MacKerell Jr., J. Phys. Chem. 99, 1846 (1995).
  43. N. Yoshii and S. Okazaki, Chem. Phys. Lett. 425, 58 (2006).
  44. D. Fincham, N. Quirke, and D. J. Tildesley, J. Chem. Phys. 84, 4535 (1986)JCPSA6000084000008004535000001.

Figures (click on thumbnails to view enlargements)

FIG.1
(a) Radial number density profile, ρ(r), of SDS micelles as a function of radial distance from the center of mass of the SDS micelle (green: hydrophobic carbon atom, red: sulfur and oxygen atoms of sulfate ion, blue: water molecule, and brown: sodium ion), (b) the mean force 〈F(r)〉, and (c) the free energy profile, ΔG(r), for penetration of a methane and a water molecule to the SDS micelle as a function of radial distance from the center of mass of the SDS micelle to those of penetrating methane and water molecules (●: methane and ○: water molecule). Open squares in (c) represent data obtained by Matubayasi et al.16 Error bars in (b) and (c) represent 80% confidence intervals. See text for details.

FIG.1 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.2
Cumulative average of force F(r) as a function of simulation time between SDS micelle and penetrating (a) methane and (b) water molecules at distance r = 1.5 nm. Solid and dashed lines are the results of MD calculations starting from two different initial configurations. Straight lines in (a) and (b) represent resultant averages of the mean force 〈F(r)〉.

FIG.2 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.3
Contributions from hydrophobic groups (green), sulfate ions (red), sodium ions (brown), and solvent water (blue) to (a) the mean force 〈F(r)〉 and (b) the free energy profile ΔG(r) for solubilization of methane molecules into SDS micelles. Black closed circle represents the total.

FIG.3 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.4
Schematic diagram of solubilization of methane into SDS micelles. The inset shows the number density profile ρ(r) of hydrophobic parts and water in SDS solution. The density of hydrophobic groups of SDS micelles is shown by shades of gray. Darker gray represents higher density. The black dot represents a methane molecule. The green and blue arrows represent mean forces 〈F(r)〉 arising from the hydrophobic groups and solvent water, respectively.

FIG.4 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.5
Contributions from hydrophobic groups (green), sulfate ions (red), sodium ions (brown), and solvent water (blue) to (a) the mean force 〈F(r)〉 and (b) the free energy profile ΔG(r) for penetration of water molecules into SDS micelles. Black open circle represents the total.

FIG.5 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.6
(a) Number density profile of a dodecane droplet solution (green: hydrophobic group; blue: water); contributions from hydrophobic groups (green) and solvent water (blue) to (b) the mean force 〈F(r)〉 and (c) the free energy profile ΔG(r) for transfer of a water molecule into a dodecane droplet. Black open circle represents the total.

FIG.6 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.7
Distribution of methane as a function of radial distance r. ● and ― were calculated from the free energy profile ΔG(r) and 15-ns MD calculations, respectively. - - - is the distribution of space that can accommodate a sphere of diameter 0.3 nm, i.e., a cavity, in the SDS micelle. All are scaled to the same value at r = 1.5 nm. The density of methane molecules per Cartesian space calculated by multiplying by the Jacobian 4πr2 is shown in the inset.

FIG.7 Download High Resolution Image (.zip file) | Export Figure to PowerPoint


Close
Google Calendar
ADVERTISEMENT

close