Введение
2. Литературный обзор
2.1. Молекулярно-термодинамические модели
2.2. Электростатический вклад в свободную энергию
2.3. Теоретическое описание перфораций
3. Теоретическая часть
3.1. Описание классической молекулярно-термодинамической модели
3.2. Предложенная модель перфорированного бислоя
4. Обсуждение результатов
4.1. Численное решение уравнения ЛПБ
4.2. Численный расчет электростатического вклада
4.3. Сопоставление значений потенциала, полученных при помощи численного
решения и существующих аналитических аппроксимаций
4.4. Концентрационные области устойчивости агрегатов различных форм.............. 23
4.5. Механизм стабилизации перфорации
4.6. Зависимость концентрационных зон стабильности агрегатов различных
морфологий от температуры и молекулярных параметров ПАВ
Выводы
Благодарности
Список использованной литературы
Приложение
А. Вывод формулы для электростатического потенциала перфорации..................... 39
Б. Вывод формулы для электростатического вклада в свободную энергию ............. 44
В. Исследование угловой зависимости поверхностного потенциала
Молекулы поверхностно-активных веществ (ПАВ) амфифильны: в них принято
выделять гидрофильную полярную голову и гидрофобный неполярный
углеводородный хвост [1, 2]. Попадая в раствор, молекулы таких веществ при
концентрациях выше критической концентрации мицеллообразования (ККМ),
образуют агрегаты различных геометрических форм. Для молекул ПАВ были
экспериментально обнаружены сферические и цилиндрические агрегаты, а также
агрегаты более сложных геометрий: перфорированные и неперфорированные плоские
бислои и везикулы и пространственные сетки [3 - 6]. Форма и размер агрегата влияют
на физико-химические и структурные свойства раствора. Способность систем к
обратимому изменению своей структуры на мезомасштабе при изменении таких
внешних условий как pH раствора, температуры, концентрации фонового электролита
и при вводе различных добавок широко используется при решении прикладных
инженерных задач [7]. Их круг очень широк и включает в себя производство топлива
[8, 9], бытовой химии и косметики, создание "умных" материалов [10] и доставку
лекарств [11, 12]. Экспериментальный подбор условий обычно требует больших
временных и финансовых затрат. Поэтому высоко ценятся надежные методы
прогнозирования структуры растворов.
Классические молекулярно-термодинамические модели, предложенные в конце
1990-х годов в работах Нагаражана и Рукенштейна [13 - 15], а также Бланкштайна [16,
17], имея в качестве входных данных молекулярные параметры вещества, позволяют
предсказать ряд термодинамических и структурных характеристик образующихся в
растворе агрегатов: свободную энергию, относительную стабильность и равновесный
размер агрегатов различных геометрий. В качестве молекулярных параметров для
ПАВ используются эффективная площадь полярной головы и длина углеводородного
хвоста молекулы. Свободная энергия в таких моделях представляется как сумма
вкладов: гидрофобного, поверхностного, деформационного, стерического и
электростатического.
Для агрегатов сложных форм, таких как перфорированный бислой и
пространственная сетка, применимость классических моделей весьма ограничена. Для4
ионных ПАВ основной трудностью является отсутствие надежного аналитического
описания электростатического вклада.
В работе [3] приведен обзор систем, в которых встречаются перфорированные
бислои и везикулы: это смеси ионных и неионных ПАВ, фосфолипиды и блоксополимеры. Образование перфорированных бислоев и пространственных сеток в
растворах амфифильных систем было обнаружено экспериментально при помощи
таких методов анализа структуры растворов, как криогенная трансмиссионная
микроскопия (cryo-TEM) [18 - 21] и малоугловое рассеяние нейтронов (SANS) [22,
23]. Методами компьютерного моделирования изучались образование ионных
каналов в липидной мембране [24] и перфорирование везикул [25].
Изначально интерес к теоретическому описанию перфораций в бислое был
вызван экспериментальным подтверждением формирования круглых пор в клеточных
стенках эритроцитов. Было показано, что размер пор зависит от ионной силы
окружающего раствора [26]. Мотивированные этими результатами, авторы работ [27,
28] впервые предложили аналитические формулы для электростатического
потенциала перфорации и электростатического вклада в свободную энергию. Для
этого они получали решение линеаризованного уравнения Пуассона-Больцмана,
рассматривая плоский, двумерный агрегат.
Целью данной работы стало развитие классической модели
мицеллообразования для описания перфорированного бислоя и пространственной
сетки на примере ионных ПАВ в растворе 1:1-электролита. В задачи работы входит:
- описание зависимости электростатического вклада в свободную энергию агрегации
от размера и формы перфорации,
- установление механизма стабилизации перфораций и превращения бислоя в
пространственную сетку,
- определение концентрационных областей фонового электролита, в которых
стабильны эти структуры.
1. Классическая молекулярно-термодинамическая модель мицеллообразования
[13 - 15] впервые была распространена для описания перфорированного бислоя и
пространственной сетки, образованных ионным ПАВ в растворе 1:1-электролита.
2. На основании анализа существующих аналитических моделей для описания
образования перфорации в бислое [27, 28] и [42, 43] в качестве электростатического
потенциала для перфорации была выбрана формула (3), полученная ранее для
фрагмента открытого тора [42, 43].
3. Сопоставление решений, полученных по аналитическому уравнению (3), с
численным решением уравнения ЛПБ позволило определить границы применимости
модели. Численное решение показало качественное отличие в поведении угловой
зависимости поверхностного потенциала для перфорации малого и большого
диаметров. Сильное влияние заряженной плоскости на распределение потенциала в
центре перфорации может привести к большой погрешности в энергии ветвления при
использовании аддитивной схемы расчета, когда энергия ветвления складывается из
энергий плоской и тороидальных частей.
4. Для модельной системы CnCTAB - KBr - H2O были найдены области
концентраций соли, при которых стабильны сферические, цилиндрические (в том
числе, ветвящиеся) и плоские агрегаты, а также образуются стабильные
перфорированные плоские бислои и пространственные сетки. Было исследовано, как
сказывается на солевых интервалах изменение температуры и молекулярных
параметров ПАВ. Полученные данные по смене стабильных морфологий при
увеличении солевого фона находятся в качественном согласии с имеющимися
экспериментальными данными [3, 5]
[1] Русанов, А.И. Мицеллообразование в растворах поверхностно-активных веществ.-
СПб.: Наука, 1992. -280 с.
[2] Rusanov, A.I. Micellization in surfactant solutions. -Reading, M.A.: Harwood Academic
Publishers, 1997.
[3] Almgren, M. Stomatosomes: perforated bilayer structures. Soft Matter. 2010, 6, 1383-
1390.
[4] Bergström, L.M.; Skoglund, S.; Edwards, K.; Eriksson, J.; Grillo, I. Spontaneous
transformations between surfactant bilayers of different topologies observed in mixtures of
sodium octyl sulfate and hexadecyltrimethylammonium bromide. Langmuir. 2014, 30,
3928-3938.
[5] Danino, D.; Bernheim-Groswasser, A.; Talmon, Y. Digital cryogenic transmission
electron microscopy: an advanced tool for direct imaging of complex fluids. Colloids and
Surfaces A: Physicochemical and Engineering Aspects. 2001, 183-185.
[6] Викторов, А.И.; Вознесенский, М.А.; Сафонова, Е.А. Успехи химии. 2015, 84, 693-
711.
[7] Victorov, A.I. Molecular thermodynamics of soft self-assembling structures for
engineering applications, J. Chem. Technol. Biotechnol. 2015, 90, 1357–1363.
[8] Maitland, G.C. Oil and gas production. Current Opinion in Colloid and Interface
Science. 2000, 5, 301-311.
[9] Chen, F.; Wu, Y.; Wang, M.; Zha, R. Self-assembly networks of wormlike micelles and
hydrophobically modified polyacrylamide with high performance in fracturing fluid
application. Colloid Polym Sci. 2015, 293, 687-697.
[10] Chu, Z.; Dreiss, C.A.; Feng, Y. Smart wormlike materials. Chem Soc Rev. 2013, 42,
7174-7203.
[11] Althans, D.; Enders, S. Investigation of the swelling behaviour of hydrogels in aqeous
acid or alkaline solutions. Mol Phys. 2014, 112, 2249-2257.35
[12] Jibry, N.; Heenan, R.K.; Murdan, A. Amphiphilogels for drug delivery: formulation
and characterization. Pharmaceutical research. 2004, 21, 1852-1861.
[13] Nagarajan, R.; Ruckenstein, E. Self-assembled systems // In "Equations of state for
fluids and fluid mixtures". Eds.: Sengers, J.V.; Kayser, R.F.; Peters, C.J.; White, H.J. -
Amsterdam: Elsevier Science, 2000.-P. 589-749.
[14] Nagarajan, R. Micellization of binary surfactant mixtures: theory // In "Mixed
surfactant systems". Eds.: Holland, P.M.; Rubingh, D.N. -Washington DC: American
Chemical Society, 1992. -P. 54-95.
[15] Nagarajan, R.; Ruckenstein, E. Theory of surfactant self-assembly: a predictive
molecular thermodynamic approach. Langmuir. 1991, 7, 2934-2969.
[16] Puvvada, S.; Blankschtein, D. Molecular-thermodynamic approach to predict
micellization, phase behavior and phase separation of micellar solutions. I. Application to
nonionic surfactants. J. Chem. Phys. 1990, 92, 3710-3724.
[17] Yuet, P.K.; Blankschtein, D. Molecular-thermodynamic modeling of mixed
cationic/anionic vesicles. Langmuir. 1996, 12, 3802-3818.
[18] Danino, D. Cryo-TEM of soft molecular assemblies. Current opinion in colloid &
interface science. 2012, 17, 316-329.
[19] Kakehashi, R.; Karlsson, G.; Almgren, M. Stomatosomes, blastula vesicles and bilayer
disks: morphological richness of structures formed in dilute aqueous mixtures of a cationic
and an anionic surfactant. J. of Colloid and Interface Science. 2009, 331, 484-493.
[20] Danino, D.; Bernheim-Groswasser, A.; Talmon, Y. Digital cryogenic transmission
electron microscopy: an advanced tool for direct imaging of complex fluids. Colloids and
Surfaces A: Physicochemical and Engineering Aspects. 2001, 183-185, 113-122.
[21] Cui, H.; Chen, Z.; Wooley, K.L.; Pochan, D.J. Origins of toroidal micelle formation
through charged triblock copolymer self-assembly. Soft Matter. 2009, 5, 1269-1278.36
[22] Nieh, M.-P.; Raghunathan, V.A.; Pabst, G.; Harroun, T.; Nagashima, K.; Morales, H.;
Katsaras, J.; Macdonald, P. Temperature driven annealing of perforations in bicellar model
membranes. Langmuir. 2011, 27, 4838-4847.
[23] Mohammed, A.; Miller, A.F.;Saiani, A. 3D networks from self-assembling ioniccomplementary octa-peptides. Macromolecular Symposia. 2007, 251, 88-95.
[24] Wang, S.; Larson, R.G. Water channel formation and ion transport in linear and
branched lipid bilayers. Phys. Chem. Chem. Phys. 2014, 16, 7251-7262.
[25] Lin, C.-M.; Chang, G.-P.; Tsao, H.-K.; Sheng, Y.-J. Solubilization mechanism of
vesicles by surfactants: effect of hydrophobicity. The Journal of Chemical Physics. 2011,
135, 045102-1 - 045102-10.
[26] Lieber, M. R.; Steck, T. L. A description of the holes in human erythrocyte membrane
ghosts. J. of Biological Chemistry. 1982, 257,11651-11659.
[27] Fošnarič, M.; Kralj-Iglič, V.; Bohinc, K.; Iglič, A.; May, S. Stabilization of pores in
lipid bilayers by anisotropic inclusions. J. Phys. Chem. B. 2003, 107, 12519-12526.
[28] Betterton, M.D.; Brenner, M.P. Electrostatic edge instability of lipid membranes. Phys.
Rev. Letters. 1999, 82, 1598-1601.
[29] Tanford, C. Theory of micelle formation in aqueous solutions. J. Phys. Chem. 1974, 78,
2469-2479.
[30] Tanford, C. The hydrophobic effect: formation of micelles and biological membranes. -
2nd ed. -New York: Wiley, 1980.
[31] Israelachvili, J.N.; Mitchell, D.J.; Ninham, B.W. Theory of Self-Assembly of
Hydrocarbon Amphiphiles into Micelles and Bilayers. J. Chem. Soc., Faraday Trans. 2.
1976, 72, 1525-1568.
[32] Mitchell, D.J.; Ninham, B.W. Micelles, Vesicles and Microemulsions. J. Chem. Soc.,
Faraday Trans. 2. 1981, 77, 601-629.37
[33] Israelachvili, J.N.; Mitchell, D.J.; Ninham, B.W. Theory of self-assembly of lipid
bilayers and vesicles. Biochimica et Biophysica Acta (BBA). Biomembranes. 1977, 470,
185-201.
[34] Iyer, J.; Blankstein, D. Molecular-thermodynamic framework to predict the
micellization behaviour of mixtures of fluorocarbon-based and hydrocarbon-based
surfactants. The Journal of Physical Chemistry B. 2014, 118, 2377-2388.
[35] Andelman, D. Introduction to electrostatics in soft and biological matter in Soft
condensed matter physics in molecular and cell biology ed. Poon, W.C.K., Andelman, D.,
Taylor & Francis, 2006, 97-122 (chapter 6)
[36] Андреев, В.А.; Власов, А.Ю.; Смирнова, Н.А. Анализ приближений для
электростатического вклада в энергию Гибсса мицеллообразования. ЖФХ. 2006, 80,
39-44.
[37] Derjaguin, B.V. On the repulsive forces between charged colloid particles and on the
theory of slow coagulation and stability of lyophobe sols. Trans. Faraday Soc. 1940, 35,
203-215.
[38] Дерягин, Б.В. Устойчивость коллоидных систем (теоретический аспект). Успехи
химии. 1979, 48, 675-721.
[39] Victorov, A.; Radke, C.; Prausnitz, J. Molecular thermodynamics for swelling of a
mesoscopic ionomer gel in 1:1 salt solutions. Phys. Chem. Chem. Phys. 2006, 8, 264-278.
[40] Victorov, A.I. Effect on morphology of a swollen ionomer gel on its salt uptake. Fluid
Phase Equilibria. 2006, 241, 334-343.
[41] Victorov, A.I. A simplified molecular-thermodynamic model of a microphaseseparated ionic gel swollen in salt solution. Fluid Phase Equilibria. 2005, 227, 9-17.
[42] Андреев В.А. Моделирование образования, роста и ветвления мицеллярных
агрегатов в растворах ионных поверхностно-активных веществ: дисс. к.х.н., СанктПетербургский гос. университет, СПб, 2006.38
[43] Andreev, V.A.; Victorov, A.I. Molecular thermodynamics for micellar branching in
solutions of ionic surfactants. Langmuir. 2006, 22, 8298-8310.
[44] May, S.; Ben-Shaul, A. Molecular packing in cylindrical micelles, in: Zana, R.; Kaler,
E.W. (Eds.) Giant Micelles: Properties and Applications, 2007, pp. 41-79.
[45] Likhtman, A.E.; Semenov, A.N. Stability of the OBDD structure for diblock copolymer
melts in the strong segregation limit. Macromolecules. 1994, 27, 3103-3106.
[46] Emelyanova, K.; Gotlib, I.; Shishkina, A.; Voznesenskiy, M.; Victorov, A. Molecular
thermodynamic modeling of self-assembly into branches and spatial networks in solution.
J. of Chem. And Eng. Data. 2016, 61, 4013−4022.
[47] Hildebrand, F.B. Advanced calculus for applications. 2-nd ed., Prentice-Hall, Inc.,
Englewood Cliffs, New Jersey, 1976, 733 p. (p. 268, problem 102)