Классические уравнения
Div І)(г) =4лр(г); (9.8)
Ј(0=-gradq>(?) (9.9)
Являются основными уравнениями электростатического поля в произвольной среде. Для нахождения с их помощью распределения электрического поля необходимо ввести граничные условия и определить связь между вектором электрической индукции
.D(r) и напряженностью электрического поля Е(г). Обычно при обсуждении диэлектрических свойств вещества предполагается,
Что значение Е в точке г зависит только от значения Е в той же точке пространства, т. е.
В(Г)=е£(Г), (9.10)
Где е=const— диэлектрическая постоянная среды.
Предположение (9.10) справедливо не всегда. В общем случае значение D(R) зависит не только от значений Е в точке г, йо и от значений Е в точках пространства вокруг г. В самом общем случае определить связь D(R) с Е(г), вообще говоря, не-
Такая нелокальная зависимость является проявлением пространственной дисперсии диэлектрического отклика среды, в отличие от временной или частотной дисперсии [432]. Ядро определяется корреляцией флуктуаций поляризации среды, присущих данному диэлектрику. Для большинства систем, рассматриваемых обычно в электростатике, пространственная дисперсия играет гораздо меньшую роль, чем временная. Это связано с тем, что для обычного диэлектрика флуктуации поляризации в соседних точках пространства довольно слабо связаны друг с другом. Поэтому ядро К (г, г’) интегрального оператора в (9.11> существенно убывает уже на расстояниях, сравнимых с атомными размерами. В этом случае Е(г’) выходит из-под интеграла по
Dzr‘, в результате чего превращается в классическое локальное выражение (9.10). Таким образом, для обычного диэлектрика пространственная дисперсия приводит лишь к малым поправкам.
Совершенно иным образом обстоит дело в проводящих средах (металлы, растворы электролитов, плазма), в которых движение свободных носителей тока вызывает нелокальность, простирающуюся на расстояния, существенно превосходящие атомные размеры. При рассмотрении электродинамических процессов в таких средах необходимо учитывать пространственную дисперсию [428, 433].
И, наконец, еще один класс сред, в которых пространственная дисперсия может играть значительную роль, — это ассоциированные жидкости, к которым, как известно, относится и вода. Хотя молекулы воды быстро и часто меняют своих «партнеров по водородным связям», в каждый момент времени любая молекула воды связана с большим числом ближних и не очень ближних молекул [434]. Очевидно, что ориентация электрического диполя молекулы воды будет зависеть не только от значения электрического поля в этой точке, но также и от ориентации связанных с ней молекул воды. Так как ориентация последних, в свою очередь, зависит от напряженности электрического поля в тех точках пространства среды, где они располагаются, то теперь радиус спадения ядра К (г, г’) существенно превосходит атомно-молекулярные размеры и определяется характерной длиной цепочки водородных связей в воде (Іо^О. бч-І нм) [433].
Возможно; однако, предполагая эту связь линейной, ее можно определить следующим образом: V |
(9.11) |
Обсудим физические причины, приводящие к нелокальности диэлектрического отклика воды. Как известно, большое значение
статической диэлектрической проницаемости воды (ео^вО) связано с тем, что вода — сильно полярная жидкость и поэтому обладает мягкой ориентационной степенью свободы (вращение молекулярных диполей), которая и вносит главный вклад в ео- Измерения временной дисперсии є(ю) в радиочастотной области [432, 435J указывают на дебаевский характер релаксации
<9Л2)
При ют>1, как видно из (9.12), є(ю)->-е<х>~5. Естественно предположить, что значение є обеспечивается более быстрыми (колебательными и электронными) модами при полностью «замороженной» ориентационной. Это означает, что при достаточно больших частотах со электрические диполи отдельных молекул воды не успевают «отслеживать» это поле, что и приводит к резкому падению диэлектрической проницаемости до значения
8ОО~5.
Весьма вероятно, что такое же «вымораживание» ориентационной моды происходит и в статических пространственно неоднородных электрических полях, причем соответствующий масштаб длин |о по крайней мере, не меньше чем длина водородной связи (~0,Зн-0,4 нм). Ограничение на |0 сверху может быть получено, например, из экспериментов по дисперсии гиперзвука в воде [436]. Отсутствие таковой дисперсии на волновых векторах гиперзвука <7<106 см"1 означает, что, по-видимому, 1 нм. Заметим, что отношение єо/єоо^іб весьма велико, поэтому заметные отклонения статического диэлектрического отклика в(^) от ео будут проявляться на пространственных масштабах <7-1, существенно больших, чем go-
Для того чтобы с помощью (9.11) определить распределение электрического поля в системе, содержащей водный электролит,
Необходимо конкретизировать вид ядра /С(г, г’). Прямое экспериментальное определение /С (г, г’) в настоящее время, по-видимому, невозможно. Дело в том, что радиус корреляции флук — туаций поляризации воды все-таки довольно мал: |о~0,5ч-1 нм [433]. Поэтому пространственная дисперсия проявляется лишь в электрических полях, пространственная неоднородность которых имеет характерный размер, сравнимый с |0- Очевидно, что создание таких полей в экспериментах сопряжено со значительными трудностями. Тем не менее целый ряд косвенных оценок позволяет ввести функцию диэлектрического отклика водного электролита с учетом пространственной дисперсии [433]. При
Этом простейшее представление К (г, г’) вводится по аналогии с хорошо известным дебаевским выражением для частотной (временной) дисперсии для неограниченной макроскопически однородной водной среды. В этом случае ядро К (г, г’) интегрального оператора в (9.11) зависит только от разности (г—Г’),
А функции D(R) и Е(г) целесообразно представить в фурье-об — разцах, сведя их к совокупности плоских волн, зависимость ко-
Торых от г дается множителем exp(t’gr). Для таких волн связь
D и Е принимает вид
(9.13)
Где
Е(7) = 1 + J К(R) ехр (— U/R)d3R; (R=Г— ?). (9.14)
V
При таком описании пространственная дисперсия означает зависимость диэлектрической проницаемости от волнового вектора q, обратное значение которого q~l определяет расстояние, на котором приложенное поле заметно меняется. Простейшая интерполяция между е0 (для q = 0) и е» (для <7^>|0-1) имеет вид, аналогичный вещественной части дебаевского выражения (9.12) [428]:
В(?) = 1 + K(q) !• <9Л5>
Нелокальная электростатика сольватационных явлений была впервые предложена Р. Р. Догонадзе и А. А. Корнышевым [437]. Затем этот подход был разработан в работах Корнышева и сотр. (см. обзоры в [428, 433]) для целого ряда других систем. Однако в этих работах анализировались в основном системы, содержащие в качестве источников электрических полей электрические заряды. Оказалось, что включение в систему электрических диполей приводит к появлению ряда новых эффектов нелокальной поляризуемости среды [429]. В этой связи интересно сравнить электрическое поле, создаваемое единичным
Точечным диполем |ір в нелокально поляризующейся среде[3]
-* 1 ->.->—>- ехр( — R/Tn) -*->■->
= [3(|1ря) — Ы+ в Jtfr (^Р")"’ <9Л6>
С электрическим полем, создаваемым единичным точечным электрическим зарядом ps:
С Ре I Е*Р( —г/Ь>) Т /п 17
Јp=>-^T«+Pe—————————————————— «• (9.17)
Первые слагаемые в (9.16) и (9.17) представляют собой обычные классические составляющие электрического поля. Отметим, что в случае диполя это слагаемое спадает существенно’ быстрее, чем в случае заряда. Вторые слагаемые обязаны своим; происхождением пространственной дисперсии E(Q). Существенно, что теперь оба эти члена спадают одинаковым образом, причем в (9.16) нелокальный член не исчезает (в отличие от классического слагаемого) при усреднении £ц(г) по всем ориентаци-
Ям п. Именно эта особенность электрических полей в нелокальной электростатике приводит к существованию сильного электрического поля вблизи нейтральной липидной мембраны [429,. 438J.
Если водный электролит занимает лишь часть пространства,, то система становится существенно неоднородной. Простейшим (но вместе с тем и одним из наиболее важных в приложениях) случаем является система с плоской симметрией, имитирующей контакты сред, помещенных во внешнее электрическое поле,, нормальное к плоскости раздела. Такая ситуация реализуется, например, в случае заряженной плоской границы диэлектрик/металл или диэлектрик/водный электролит.
Решение электростатических задач для неоднородных систем значительно сложнее, чем для однородных, так как теперь ядро*
Нельзя представить в виде функции (г—г’) и, следовательно, переход к фурье-анализу оказывается неэффективным. Однако в случае плоской границы раздела можно перейти к цилиндрическим координатам и, представив ядро /С (г, г’) в виде
K(r,~r’)=K(z, Г’, R — R’) (9.18).
(где Z — координата, нормальная к границе; R — вектор, лежащий в плоскости границы), перейти к фурье-компонентам по латеральным координатам. Тем не менее вопрос о связи (9.18) с
Объемной функцией диэлектрического отклика ВОДЫ Є (г—г’)’ остается открытым. Очевидно, что граница раздела оказывает влияние на эту функцию, и вблизи поверхности свойства диэлектрического отклика электролита будут существенно отличаться от объемных свойств. В последние годы эта проблема активно обсуждалась в работах по электродинамике металлов; [433]. Возможность перенесения результатов этих работ на контактные системы с водным электролитом стала предметом дискуссии [439], в ходе которой было высказано предположение [440], что (9.18) вплоть до самой границы совпадает с ее объемным выражением (так называемое «диэлектрическое приближение»):
К(z, Z‘, R — R‘)=K(z — Z‘, R — R‘)=K{r — 7′). (9.19>
Дискуссия показала, что приближение применимо лишь к ограниченному классу систем, более общим является так называемая модель «зеркально-диффузионного отражения» [441]:
К (г, гR — R’)=K(z — z R-R’) —
— чК(г+г’, R — R’), (9.20)
•где к— «коэффициент зеркальности» — может меняться от 0 до 1.
Для водного электролита у является феноменологическим параметром, отражающим специфику границы раздела.
Так как цель настоящей работы — продемонстрировать эффекты нелокальной электростатики в водно-липидных дисперсиях, то в основном мы будем использовать простейшее «диэлектрическое приближение» (9.19), хотя некоторые, более тонкие, эффекты связаны с использованием зеркально-диффу — — зионной модели [442].
Таким образом, распределение электрических полей в рамках нелокальной электростатики определяется уравнениями
(9.8) и (9.9), в которых связь между индукцией D(R) и на-
Лряженностью электрического поля Е(г) описывается в случае объемного электролита выражениями (9.13) — (9.15), а для неоднородных, контактных систем — выражениями (9.13), (9.14), (9.15) и (9.20).
Нелокальная поляризуемость электролита оказывает существенное влияние на структуру ДЭС, образующегося вблизи фосфолипидной поверхности. Однако необходимо отметить, что в общем случае экранирование поверхностных источников электрических полей имеет два принципиально разных по физической природе механизма [443]. Первый механизм связан с экранированием поверхностных источников ионами электролита, второй обусловлен реакцией самого растворителя на поверхностные источники. По существу, оба эти механизма имеют нелокальный характер и определяются корреляциями флуктуаций электрических полей в соседних точках пространства. В первом случае такие флуктуации обусловлены флуктуациями концентрации ионов, характерный радиус корреляций которых есть де — баевский радиус х-1. Во втором случае флуктуации электрических полей связаны с флуктуациями поляризации в электролите, радиус корреляции которых а-1.
Теория Гуи — Чепмена (а также теория Дебая для электро — .литов) описывает лишь первый механизм экранирования, в предположении, что реакция самого растворителя на внешнее электрическое поле имеет классический локальный характер. Выше был описан второй механизм экранирования поверхностных источников, учитывающий нелокальную природу диэлектрического отклика растворителя.
Распределение потенциала электрического поля вблизи фосфолипидной поверхности с учетом нелокальной поляризуемости электролита определяется уравнением [429]:
D °С 4я f d
»"+ ЧГ) (9.21)
О
Граничные условия на потенциал и на индукцию в точке Z = 0 обычны:
Ф(0) =фт; D(0) —етЕт. (9.22)
Одновременно задается условие конечности потенциала при Z->~Oo. Здесь фт и Ет — потенциал и поле в мембране; гт — диэлектрическая проницаемость мембраны.
Уравнение (9.21) представляет собой обобщение теории Гуи — Чепмена на случай нелокально поляризующегося электролита. Если пренебречь эффектами нелокальной поляризуемости (для чего достаточно положить Р = 1), то выражение — (9.21) переходит в классическое уравнение Гуи — Чепмена.
Насколько важен учет эффектов нелокальной поляризуемости электролита для фосфолипидных дисперсий? Чтобы ответить на этот вопрос, сравним распределение электрического потенциала, полученное из (9.5), с решением уравнения (9.21). Последнее для разбавленного электролита (Х<Сйо) в полупространстве z>L имеет следующий вид [427]:
2л 2л
Ф(г) — — аеиехр( — Х2) — Г-^J — (Р — 1 )°еаехР( — аг) (9.23>
Ае*= стг[ ехр (к1) +ехр (— к1) ] +
+ц, и[ехр (NL) — ехр (— kL)] + сі,; (9.24)
АЄа = а;[ехр (Al) — т) ехр (— Al)] + ца[ехр (Al) — f-rj ехр (—а/)]; (9.25)
Оь=— 2ц(УУ— 1) (х2/а0) ехр (— AL), (9.26) где Р = е0/е«,; а0=|о-‘; а = ргао= g-i.
Основное отличие полученного решения от классического заключается в том, что теперь распределение потенциала описывается двумя экспоненциальными членами и при условии
ОекОеа < о (9.27)-
Становится немонотонным, знакопеременным.
Как уже было сказано выше, экранирование поверхностных источников электрических полей в водных электролитах происходит по двум механизмам. При выполнении условия (9.27) одна и та же поверхность по одному механизму может быть заряжена положительно, а по другому — отрицательно. В результате распределение потенциала вблизи фосфолипидной поверхности становится немонотонным, и в точке
Zl ~ In [- (Р — 1)о, а/о. и]/(а+х) (9.28) Образуется потенциальный барьер, равный
Фе ~ 2лаех/(в0х). (9.29)
Физическая природа такой немонотонности связана с особенностью поведения поверхностных электрических диполей в среде с двумя механизмами экранирования, которая заключается в том, что вклад диполей в электрическое поле не зависит от характерной длины экранирования (в отличие от вклада зарядов). Следовательно, для различных механизмов экранирования эффективный поверхностный заряд одной и той же поверхности будет различен.
Отметим еще одну особенность распределения (9.23). Как видно из (9.23) — (9.25), вклад поверхностных диполей в электрическое поле значителен при любой, в том числе нулевой (L=0) степени гидратации. Это важное отличие нелокальной электростатики позволяет объяснить природу нового эффекта — структурной составляющей расклинивающего давления в фосфолипидных дисперсиях (см. ниже).
Рассмотренные эффекты могут иметь значение и для биофизики. В частности, поверхностные диполи фосфолипидных мембран могут оказывать влияние на электрогенные биофизические процессы, причем это влияние зависит от степени гидратации поверхности. Биологически активные ионы (например, Са2+), как известно, способны менять степень гидратации фос — фолипидной поверхности [430]. Возможно, регуляторная функция этих ионов связана с изменением структуры ДЭС в результате уменьшения степени гидратации поверхности под влиянием этих ионов.
Итак, структура ДЭС вблизи фосфолипидной поверхности может существенно отличаться от предсказаний классической теории Гуи — Чепмена. Эти отличия связаны с двумя особенностями водно-липидных систем. Во-первых, водный электролит является нелокально поляризующейся средой и, во-вторых, поверхностные источники электрических полей (заряды и диполи) распределены в приповерхностном слое. Влияние этих особенностей на структуру ДЭС приводит к тому, что распределение потенциала вблизи фосфолипидной поверхности становится немонотонным, знакопеременным.
Следует отметить, что прямое экспериментальное исследование рассмотренных эффектов связано с определенными трудностями. Немонотонность распределения потенциала проявляется лишь на расстояниях 1—2 нм от поверхности. Экспериментальное исследование параметров ДЭС на таких расстояниях в настоящее время вряд ли возможно. Вместе с тем на больших расстояниях трудно отделить вклад истинного поверхностного заряда At от дипольного вклада хц в измеряемые величины. Однако в некоторых экспериментах рассмотренные эффекты могут проявиться косвенным образом. Например, уже давно отмечалось, что в некоторых случаях, изначально нейтральная фосфолипидная мембрана создает такой же ДЭС, как и заряженная [423]. Эти результаты можно трактовать в рамках нашего подхода как проявление эффективного поверхностного заряда, обусловленного поверхностными диполями.