CC BY
7
УДК 548.31
ЦЕПОЧЕЧНЫЙ АССОЦИАТ МОЛЕКУЛ КОНИЧЕСКОЙ СИММЕТРИИ, ИЛИ ЗАДАЧА «О СПИЧКАХ» А.М. Банару, В.Л. Федоров, Г.А. Банару
(кафедра физической химии; e-mail: banaru@phys.chem.msu.ru)
Рассмотрена модель формирования цепочечного ассоциата из линейных молекул с помощью контактов голова-к-голове и голова-к-хвосту. Найдена дискретная функция распределения ассоциатов по числу контактов одного типа.
Ключевые слова: ассоциат, контакт, нормальное распределение.
Моделирование ассоциации молекул в конденсированных фазах при фазовых переходах представляет собой важную и интересную задачу. Молекулярная динамика и ее вариации во многих случаях наглядно иллюстрируют формирование фрагментов новой фазы на микроуровне как для простых, так и для довольно сложных систем. Однако такие расчеты, как правило, не дают глубокого объяснения причин наблюдаемого и лишь косвенно отвечают на вопрос о вероятности реализации данной конкретной картины (через статистический ансамбль). В связи с этим до сих пор остается актуальной математическая формализация поведения таких систем, которая служит средством не столько описания и обобщения, сколько выработки условий, необходимых для наблюдения того, что существует в действительности, безотносительно конкретных физико-химических условий. Последние, по сути, представляют собой контекст происходящего процесса. Так, в серии работ [1-3] нами было установлено, что наблюдаемая картина распределения органических кристаллогидратов по структурным типам сеток (И20)ш определяется ничем иным как преимущественной тенденцией молекул воды к образовыванию пятиугольных циклов, и органическая компонента кристалла служит лишь контекстом реализации водной сетки, влияющим на ее топологию, но не влияющим на топологические инварианты (в частности, на протоноизбыточность [1]).
Сравнительно недавно В.А. Дуров и А.П. Моска-лец [4] разработали метод рекуррентных последовательностей, позволяющий вычислять интегральные и дифференциальные физико-химические характеристики ассоциированных жидких систем, описывать их диэлектрические свойства. Центральное место в этом теоретическом построении занимают цепочечные ас-социаты (как наиболее простой и наглядный случай). Рекуррентность применительно к молекулярным си-
стемам означает, что результат присоединения новой молекулы к уже сформированному ассоциату зависит от предыдущего акта присоединения, и общая вероятность реализации такого ассоциата зависит от вероятностей присоединения на каждой стадии.
Если молекула обладает точечной группой симметрии P, то потенциальное поле, создаваемое этой молекулой, согласно известному принципу Неймана [5], обладает точечной группой P'з P. Так называемый метод симметрии потенциальных функций, разработанный П.М. Зорким [6], позволяет делать вероятностный прогноз симметрии молекулярных ассоциатов разной размерности (островных, цепочечных, слоистых, каркасных) на основании начальной симметрии потенциального поля. Прогноз затрудняется в случае конформационно гибких молекул, когда потенциальное поле «подстраивается» под частично сформированный ассоциат, и в этом случае часто возникают структуры с несколькими симметрически неэквивалентными молекулами (Z' > 1). Ассоциат, в ко -тором равные молекулы занимают всю совокупность симметрически эквивалентных точек потенциального поля данной молекулы, называют нормальным ассоциатом. В случае P = Cavv нормальный ассоциат в реальности возникнет только тогда, когда контактирующие молекулы расположены на Сш. (В противном случае эта ось бесконечно размножит симметрически эквивалентные позиции, и они не смогут быть заняты молекулами одновременно). Таким образом, нормальным ассоциатом конических молекул будет лишь цепочечный ассоциат с упорядочением вдоль главной оси. Этот ассоциат и является объектом настоящего исследования.
В феноменологическом описании межмолекулярных взаимодействий широко распространено деление участков той или иной структуры на гидрофильные и гидрофобные. Молекулы поверхностно-активных
веществ, как правило, имеют длинный гидрофобный «хвост» и небольшую гидрофильную группу в качестве «головы», что в водном растворе способствует образованию мицелл. При расположении таких молекул в цепь они будут иметь очевидную тенденцию к образованию контактов голова-к-голове и хвост-к-хвосту. Если же не принимать во внимание «химизм» такого упорядочения, то контакты голова-к-голове и голова-к-хвосту образуются равновероятно. Подобная картина в реальности будет наблюдаться, скажем, при радикальной полимеризации некоторых мономеров, где образуется полимер нерегулярного строения.
Оценка вероятности соотношения разнотипных контактов в бесконечно длинном ассоциате при такой (нерегулярной) ассоциации - цель настоящей работы.
Цепочечный ассоциат
Найдем вероятность того, что в цепочке, составленной из т «спичек», соединенных концами, имеется п контактов голова-к-голове, если каждая «спичка» ориентирована равновероятно в прямом и обратном направлениях, и оценим наибольшее значение этой величины, если оно есть. Обозначим величину Ртп.
Наличие контакта голова-к-голове означает перемену ориентации «спички». Пронумеруем такие контакты. Очевидно, что их число к (число перемен ориентации) меньше, чем т (к < т - 1). Разобьем множество всех возможных цепочек на два подмножества: Я и Ь. Пусть в подмножестве Я цепочки начинаются «спичкой», повернутой головой направо, а в Ь - наоборот.
Рассмотрим цепочки подмножества Я. Пары голова-к-голове представляются нечетными переменами ориентации. Их число п равно частному к/2, взятому с избытком, т.е. \к /2|. Теперь рассмотрим цепочки подмножества Ь. Здесь пары голова-к-голове представляются четными переменами ориентации. Их число п равно частному к/2 с недостатком, т.е. |_к /2] (рис. 1).
Чтобы сконструировать цепочки, содержащие п требуемых пар голова-к-голове, нужно выбирать значения к, равные (2п - 1), 2п, (2п + 1):
1) при к = 2п - 1 цепочки содержатся только в Я;
2) при к = 2п содержатся и в Я, и в Ь;
3) при к = 2п + 1 содержатся только в Ь.
В общем случае число искомых цепочек Ытп можно определить по формуле
Ып _ С 1 * ' 1
2 п—1
-2С,
2п
С
2 п+1
-^т — 1 т—1 т —1
(т — 1) • (т — 2) •... • (т — 2п + 2) (2п +1)!
т • (т +1).
Последовательность начинается головкой вправо
© • ф
® • "V"
Последовательность начинается головкой влево
• Ф •
• Ф @
перемен 3, п = 2 перемен 4, п = 2
перемен 3, п = 1 перемен 4, п = 2
т = 6, вариантов 2"' = 64
Рис. 1. К выводу числа ассоциатов с п контактами голова-к-голове (пояснения в тексте)
Частные случаи возникают, когда 2п - 1 < т - 1 или 2п + 1 > т - 1. Таким образом, искомая в задаче вероятность равна
рп _
Ы" 2т
Чтобы оценивать пределы данной величины, можно использовать довольно грубое приближение по Стирлингу:
г! « V2кп
Для иллюстрации распределения цепочек ассоциированных «спичек» по числу контактов была составлена программа, рассчитывающая Ытп для т < 1000. При достаточно больших т распределение близко к нормальному. Однако лишь при нечетных т оно симметрично, причем для т = 4q + 1 оно унимодально (рис. 2) с наибольшим значением при п _ \_т /4], а для т = 4q + 3 бимодально с наибольшими значениями, соответствующими |т /4] и \т /4| . При четных т распределение остается унимодальным (мода |т / 4]), но становится асимметричным, причем при т = 4д его левая ветвь выше правой, а при т = 4д + 2, наоборот, правая выше левой.
Подставив п = т/4 в написанную выше формулу для N п, получим
Ытах _
(т +1)!
(т +1)
т+3/2
((т/2 +1)!)2 еТ2тг (т/2 +1)
т+3
При очень больших значениях т очевидно, что т + 3 ~ т + 3/2 ~ т + 1 ~ т, т/2 + 1 ~ т/2, поэтому
№
глр2ж
,„3/2
пг
Поделив это выражение на 2т , мы обнаружим, что
г» макс г -3/2
Р т с точностью до константы убывает как т , т. е. довольно быстро, и при т ^ да обращается в 0, что
Рис. 2. Распределение Nnm (n) для m = 89. Параметры распределения ц = 0,489; с2 = 0,225
и бывает для всякой случайной величины, распределенной нормально на бесконечном участке.
На основании расчетов мы можем заключить, что распределение однотипных контактов в достаточно длинной цепи из линейных молекул близко к нормальному. Для такой цепи наиболее вероятное отношение числа контактов к длине цепи равно 0,25.
Впрочем, можно привести примеры, когда это отношение в действительности не совпадает с ожидаемым. Так, если мы построим из «спичек» ассоциат в порядке последовательности Фибоначчи аналогично одномерному квазикристаллу [7], то на каждую молекулу будет приходиться (3 - ^5)/2 ~ 0,382 контакта голова-к-голове.
СПИСОК ЛИТЕРАТУРЫ
1. Banaru A., Slovokhotov Yu.L. // CrystEngComm. 2010. 12. P. 1054.
2. Banaru A. // CrystEngComm. 2011. 13. P. 212.
3. Банару А.М., Банару Г.А. // Вестник Моск. ун-та. Сер. 2. Химия. 2012. 53. С. 8.
4. Дуров В.А., Москалец А.П. // Вестник Моск. ун-та. Сер. 2. Химия. 2011. 52. С. 175.
5. Зоркий П.М. Симметрия молекул и кристаллических структур. М., 1986.
6. Zorky P.M. // J. Mol. Struct. 1996. 374. P. 9.
7. Valy Vardeny Z., Nahata A., Agrawal A. // Nature Photonics. 2013. 7. P. 177.
Поступила в редакцию 20.05.2013
A CHAINED AGGREGATE OF CONIC MOLECULES OR THE PROBLEM OF MATCHES
A.M. Banaru, V.L. Feodorov, G.A. Banaru
(Division of Physical Chemistry)
The model of forming chained aggregate of linear molecules by contacts head-to-head and head-to-tale was explored. The discrete distribution function for such aggregates versus number of contacts of the same type was found.
Key words: aggregate, contact, normal distribution.
Сведения об авторах: Банару Александр Михайлович - ст. препод. кафедры физической химии химического факультета МГУ, канд. хим. наук (banaru@phys.chem.msu.ru); Федоров Владимир Леонидович - ст. препод. кафедры информатики физико-математического факультета СмолГУ; Банару Галина Анатольевна - доцент кафедры математики физико-математического факультета СмолГУ, канд. физ.-мат. наук (banaru@mail.ru).
