УДК 539.196; 531.53
Б01: 10.15587/2313-8416.2019.170742
АНАЛИЗ ВРАЩАТЕЛЬНОГО ПОТЕНЦИАЛА ДВУХЧАСТОТНЫХ КОЛЕБАНИЙ МОЛЕКУЛЫ ВОДЫ
© Н. Т. Малафаев
Рассмотрены вращательные колебания молекулы воды с помощью модели двухчастотного физического маятника и показано, что тип его потенциала правильно описывать как тороидальный. Показано, что эллиптичность тороидального потенциала в неоднородном поле сил с ростом показателя степени п уменьшается быстрее, чем для эллипсоидального, однако оба потенциала становятся близкими. Уменьшение эллиптичности тороидального потенциала в этом поле может вести к расширению области существования эллипсоподобных колебаний двухчастотного маятника в сторону его меньших скоростей
Ключевые слова: молекула воды, двухчастотный маятник, неоднородное поле сил, тороидальный потенциал
1. Введение
Вода и ее свойства [1, 2] имеют громадное значение для жизни людей. Молекулы воды в жидкой фазе находятся в конденсированном состоянии и совершают разнообразные колебания вблизи своих положений равновесия. Малая величина коэффициента диффузии говорит о малости числа молекул, которые разорвали связи и перемещаются на новые позиции. Поэтому большинство молекул совершают достаточно длительные колебания около положений равновесия в минимумах своего энергетического потенциала до нового изменения своих соседей. Однако энергетический потенциал таких колебаний и его тип не был рассмотрен, чему посвящена данная работа.
2. Литературный обзор
Влияние колебаний атомов и молекул на фазовое состояние веществ рассматривается в теории эффекта Яна-Теллера [3]. В ней выделяются из хаоса колебаний те моды, которые отвечают за переходы в новое фазовое состояние. В [4] показано, что наиболее важными из них для воды при ее плавлении являются вращательные (либрационные) колебания, приводящие к появлению коллективизированных незатухающих вращательных колебаний молекул и к изогнутым водородным связям в жидкой фазе воды. Молекула воды имеет три момента инерции и, соответственно, три вращательные частоты. Для более простой модели колебаний молекулы была разработана модель двухчастотного маятника [5], заменяющая реальный двух-частотный физический маятник двухчастотным математическим. Показано, что при малых амплитудах ко-
лебаний совершаются двухчастотные колебания, а при больших - одночастотные прецессирующие эллипсо-подобные колебания маятника. Данный маятник совершает свободные колебания. Это позволяет вследствие отсутствия возмущений потенциала увидеть общую картину колебаний. В отличие от метода молекулярной динамики [1, 6], когда в хаосе смещений всех молекул трудно рассмотреть вращательные колебания и их особенности, тем более, что такая задача недооценивалась и не ставилась.
Механическая модель двухчастотного маятника описана в [7]. Однако в этой работе не рассматривались его энергетический потенциал и взаимосвязь колебаний по осям, считая их независимыми. Часто потенциал молекул рассматривается лишь в общем виде без взаимосвязи колебаний между собой и конкретных расчетов траекторий колебаний атомов в молекуле [1, 2].
Задача о потенциале гравитирующего однородного кругового тора обсуждается в [8]. Используя стандартные полные эллиптические интегралы получен пространственный гравитационный потенциал тора в произвольной точке и найдено семейство эк-випотенциалей. Хотя в этой задаче присутствуют две точечные массы, в ней рассмотрены колебания относительно оси, перпендикулярной исследуемой, и не рассмотрен случай неоднородного поля сил.
В [9] рассмотрен электрический тороидальный потенциал ловушек в масс-спектрометре. Однако характер движения частиц в нем отличен от исследуемого - аксиальный, вдоль оси тора и не ограничен по энергии.
Ряд работ по двухчастотным колебаниям посвящен различным механическим или другим поглотителям вибраций. В [10] представлены результаты экспериментальных и аналитических исследований вибрационно-маятникового поглотителя. Показано, что наилучшее поглощение энергии колебаний происходит вблизи регулируемых резонансных частот системы.
Влияние направленного межмолекулярного неоднородного поля сил (НПС), задаваемого как FНПС=Fcosпв, на величину силы, прецессию и потенциал сферического маятника, рассмотрено в [11]. С ростом показателя НПС п и отклонений маятника от равновесия (угла изгиба связей) наблюдается уменьшение сил, сужение потенциала и понижение частот колебаний, что объясняет понижение либрационных частот колебаний в воде с нагревом. Потенциал сферического маятника является симметричным и одно-частотным, и потому не годится для описания двух-частотных колебаний.
Полный анализ двухчастотных колебаний молекул воды требует правильного описания их энергетических потенциалов. Этот потенциал должен быть анизотропным и обеспечить наличие двух частот колебаний, связанных с отличием моментов инерции молекулы по осям (к=1х/1у, для воды £=1,5 или 3 [5]). Это ведет к различной кривизне потенциала по осям У и Х около точки его равновесия и появлению эллиптичности сечений постоянного уровня потенциала. Возникает вопрос и о поведении потенциала с ростом отклонений маятника. При наличии НПС [5] потенциалы становятся направленными вдоль осей связи и значительно сужаются по углу. В результате становятся возможными разрывы связей при относительных отклонениях ¿>0,5...0,6 или углах отклонения в>30° [2, 6]. Для НПС был определен оптимальный по таким разрывам связи степенной показатель, равный п=8 [5]. Также наличие НПС благоприятствует появлению эллипсоподобных колебаний маятника.
В методе молекулярной динамики [1, 6] атом -атомные потенциалы взаимодействия [12] являются аксиальносимметричными и отвечают лишь за трансляционные колебания. Несимметричность потенциалов возникает лишь при учете мультипольных взаимодействий от дополнительных соседних молекул (димеры [13]) и для неплоских моделей молекул воды. Вращательные колебания и их потенциалы часто игнорируются и не анализируются.
В модели двухчастотного физического маятника имеем для него различные эффективные длины по осям колебаний. Это связано с тем, что по оси Х эффективная длина равна длине маятника, а его движение по оси У соответствует движению центра тяжести двух разнесенных на угол во относительно точки подвеса масс двухчастотного маятника, а также имеющего меньшую эффективную длину [5]. Поэтому нужно математически правильно описать этот потенциал, причем как вращательный. Неправильное описание потенциала двухчастотного маятника может вести к неверному расчету траекторий маятника, а также к невыполнению для него законов сохранения.
3. Цель и задачи исследования
Цель работы - определение возможного типа вращательного потенциала двухчастотного маятника как потенциала вращательных колебаний молекулы воды и определение особенностей этих колебаний.
Для достижения поставленной цели были решены следующие задачи:
1. Рассмотрение возможных типов потенциала двухчастотного маятника.
2. Выбор «физически» правильного типа потенциала двухчастотного маятника.
3. Рассмотрение особенностей тороидального потенциала, сравнительно с эллипсоидальным, в том числе, их особенностей в НПС.
4. Модели потенциала двухчастотного маятника и их обсуждение
Моделирование колебаний двухчастотного маятника как модели вращательных колебаний одиночной молекулы проводится в пакете Ма1ЬаЬ. Постоянная во времени величина потенциала означает, что рассматривается осредненный во времени потенциал при малых возмущениях его вследствие тепловых трансляционных и вращательных колебаний соседних молекул. Это более близко к состоянию молекул вблизи точки плавления, когда энергии колебаний молекул малы относительно величины его потенциала в минимуме.
Рассмотрим возможные типы потенциала двухчастотного маятника. Для общей точки подвеса эффективных длин двухчастотного маятника не будет наблюдаться общего минимума потенциала. Для эллиптического потенциала выполняется постоянство эллиптичности для всех линий постоянного уровня потенциала, однако не выполняется постоянство эффективной длины маятника по оси У. И, соответственно, будут неверно меняться его силы и частоты с ростом угла отклонения от равновесия в.
Величину эллиптичности е сечения потенциала по оси 2 находим как е=а/Ь, где а, Ь - размеры эллипса сечения по осям X и У (а>Ь) и, которая характеризует различие частот в минимуме потенциала
ау Г как: = V е.
Поэтому для двух эффективных длин двухча-стотного маятника их точки подвеса должны быть на разной высоте с общим минимумом потенциала, аналогично механической модели двухчастотного маятника [6]. Данный потенциал является тороидальным и его можно описать параметрически: углами отклонения а и в по осямХи У, соответственно:
х = (AR + R • cos ß) sin а,
y = R sin ß,
z = (AR + R • cos ß) cosa,
где 4=1 м, 4 =1/л/к =Я м, АЯ=1-Я м.
Потенциал маятника пропорционален координате г и равен:
где ио - величина потенциала в его минимуме (здесь взято ио=-1 Дж и к=3).
Поверхность потенциала имеет форму тора, хотя для молекулы воды закрытым - без дырки внутри (1х/1у=у[к <2). Его ось расположена горизонтально - по оси У. На рис. 1, а показана нижняя часть потенциала - ниже точки подвеса. На высоте И=Я от минимума потенциал имеет наибольшую ширину по оси У, а затем с ростом высоты сужается. Это говорит о том, что скатывающие силы по оси У как бы изменяют знак, но для вращательного движения такой проблемы нет. Линии постоянного уровня - сечений потенциала плоскостями ХУ вблизи минимума имеют форму эллипсов (рис. 1, б), а с ростом высоты И линии сгущаются. Величины эллиптичности для эллипсов рассчитаны по максимальным отклонениям потенциала маятника вдоль осей X и У как е=а/Ь. При больших высотах (И>К) линии уровня являются кривыми 4-го порядка и их эллиптичность становится формальным параметром.
б
Рис. 1. Тороидальный потенциал двухчастотных колебаний маятника для отношения моментов инерции к=3 в однородном поле сил (п=0, ио=-1): а - вид потенциала; б - линии постоянного уровня потенциала
На рис. 2, а показано вид тороидального потенциала в НПС, а на рис. 2, б - показаны линии постоянного уровня этого потенциала. Видим, что там, где в однородном поле сил на высотах больших Я (рис. 1) скатывающие силы по оси У являются отрицательными, в НПС величины потенциалов в этих областях малы, а силы становятся пренебрежимо малыми. Видим существенную направленность потенциала вдоль оси 2 (оси связи молекулы) в НПС и сгущения линий уровня этого потенциала для малых и средних отклонений и уменьшение эллиптичности, сравнительно с рис. 1, б. Вследствие наличия НПС потенциальная энергия становится малой для относительных отклонений от равновесия ¿>0,5 и видим расширение между линиями уровня потенциала по оси X при уменьшении отклонений по оси У.
б
Рис. 2. Тороидальный потенциал маятника для отношения моментов инерции к=3 в НПС (п=8): а - вид потенциала; б - уровни постоянного потенциала
На рис. 3, а показаны зависимости эллиптичности уровней тороидального потенциала от высоты И=1-2 для величин степенного показателя п=0, 2 и 8. Видим, что эллиптичность растет с высотой И и уменьшается с ростом степенного показателя п, что согласуется с рис. 2, б. Уменьшение эллиптичности связано с уменьшением угла в и его косинуса по оси У сравнительно с осью X в НПС с п=8 для линий постоянного уровня потенциала. Уменьшение эллип-
U=Uo z
а
тичности для тороидального потенциала в НПС с ростом его степенного показателя п будет вести к сближению частот колебаний по осям.
Для сравнения рассмотрено зависимости эллиптичности для линий постоянных уровней для эллипсоидального потенциала в НПС для тех же степенных показателей п и отношения моментов инерции к=3 (рис. 3, б). Видим уменьшение эллиптичности в НПС с ростом степенного показателя п, аналогично рис. 3, а. Однако величины эллиптичности для эллипсоида большие для заданных п, чем для тора, хотя с высотой меняются медленнее: при п=8 в интервалах высот 0<И<0,5 м ееИ=1,108±0,007 и еОг8=1,056±0,021 для эллипсоида и тора, соответственно. В однородном поле сил (п=0) для эллипсоида она остается постоянной (ем=л[к =1,732), тогда как для тора для малых высот (И<0,5 м) она уменьшается. При высотах И>К дополнительный рост эллиптичности для тора связан с отрицательным загибом его поверхности по оси У и физического смысла не имеет (рис. 1, б).
Эллиптичность
3 2,5 2 1,5 1
0,5 0
0 0,2 0,4 0,6 0,8 1 h, м - e(0) -л- e(2) -о- e(8)
а
Эллиптичность
1,5
0,5
h, м
0 0,2 0,4 0,6 0,8 1 о е(0) -а-е(2) □ е(8)
б
Рис. 3. Зависимость эллиптичности уровней потенциалов от высот И при отношении моментов инерции к=3 для степенных показателей НПС п=0, 2 и 8 для: а - тороидального и б - эллипсоидального потециалов
На рис. 4 показано зависимости сечений 1(8) тороидального (Т) и эллипсоидального (Е) потенциалов по их осям Б=Х, У для отношения моментов инерции к=3, где обозначение соответствует сечениям потенциала вдоль осей 1 и 8 для показателей НПС п=0 или 8. Видим, что в однородном поле сил (п=0) сечения представляют отрезки дуг окружностей или эллипсов. Однако в НПС (п=8) потенциалы уменьшаются по координате Ъ и анизотропно сужаются по углу в, эллиптичность их также уменьшается, что видно по округлению эллипсов (рис. 1, б, рис. 2, б). Вследствие ослабления сил взаимодействий с ростом угла в наибольшая ширина потенциала по оси У в НПС (п=8) достигается при г=-0,0063, тогда как при п=0 - при г=-0,4226.
о -0,1 -0,2 -0,3 -0,4 -0,5 -0,6 -0,7 -0,8 -0,9 -1
0 0 20 4\ 0 6 0, ОЭ
/
//
-zy0T
а
-zy0E
Z8
0 -0,1 -0,2 -0,3 -0,4 -0,5 -0,6 -0,7 -0,8 -0,9 -1
0 0 - 2 0 о ооооа> 4 <р0 6 0 ,8 1
X, Y
zx8
zy8T
б
zy8E
Рис. 4. Сечения тороидального и эллипсоидального потенциалов по оси 1 вдоль осей X и У для отношения моментов инерции к=3 для степенных показателей: а - в однородном поле сил (п=0); б - в НПС с п=8
Видим, что в НПС область отрицательных сил ничтожна даже для статической модели потенциала, тогда как разрывы связей вследствие колебаний молекул в воде весьма вероятны при относительных смещениях й>0,5 [2, 6]. Эти графики численно подтверждает описания тороидальных потенциалов к рис. 1, 2.
Z0
X, Y
zx0
2
0
5. Результаты исследований и их обсуждение
Видим, что зависимости потенциалов (рис. 4) и действующих сил вдоль оси X совпадают для обоих потенциалов. Потенциалы вдоль оси Y в однородном поле сил (n=0) значительно отличаются (рис. 4, а). Однако для случая НПС (n=8) потенциалы вдоль Y (рис. 4, б) и эллиптичности сечений (рис. 3) обоих потенциалов по этой оси (для координаты Z8<-0,01 м) различаются в пределах ±(2...3) %. Для критического относительного смещения ¿=0,5 [2, 6] имеем Z^^-0,143 м по оси Y и Z^^-0,274 м по оси X
Это также может говорить о большей устойчивости колебаний молекул воды вдоль оси X даже при учете колебаний соседних молекул. Отсюда также можно оценить величину критического потенциала и энергию возмущений при разрыве связей молекулы воды как:
EKp^-UKp^ZKpg*=(0,1...0,Uyg*
Таким образом, если в однородном поле сил тороидальный и эллипсоидальный потенциалы значительно различаются (рис. 4), то в НПС оба потенциала становятся близкими, а потому результаты расчетов траекторий колебаний двухчастотного маятника для них также будут близкими. Поэтому расчеты в [5] и последующих работах, проведенные для эллипсоидального потенциала, не приведут к качественным отличиям с расчетами для «правильного» тороидального потенциала, а лишь к малым количественным отличиям траекторий колебаний двухча-стотного маятника и его динамических параметров.
Большая величина эллиптичности для эллипсоидального потенциала сравнительно с тороидальным
связана с более медленным уменьшением угла в, его косинуса и уменьшения потенциала в направлении оси Y с ростом высоты h. Большее уменьшение эллиптичности для тороидального потенциала в НПС сравнительно с эллипсоидальным может привести к большему сближению частот малых колебаний по осям, и к увеличению области существования одночастотных, пре-цессирующих, эллипсоподобных колебаний двухча-стотного маятника в сторону его малых скоростей.
Полученные результаты для модели двухча-стотного маятника могут быть применимы, как первое приближение, для описания вращательного потенциала молекул воды и ее вращательных колебаний для его стационарного случая - отсутствия либо малости внешних возмущений потенциала молекул.
6. Выводы
1. Рассмотрены возможные типы потенциала двухчастотного физического маятника. Показано, что его потенциал вращательных колебаний правильно описывать с помощью тороидального потенциала.
2. Показано, что в неоднородном поле межмолекулярных сил с ростом показателя степени n эллиптичность тороидального потенциала уменьшается быстрее, чем эллипсоидального. Так как в НПС оба потенциала становятся близкими, это должно приводить в расчетах к близким по величинам траекторий и всех параметров двухчастотного маятника.
3. Определено, что большее уменьшение эллиптичности тороидального потенциала в НПС может вести к расширению области существования эл-липсоподобных колебаний двухчастотного маятника в сторону его малых скоростей.
Литература
1. Антонченко В. Я., Давыдов А. С., Ильин В. В. Основы физики воды. Киев: Наукова Думка, 1991. 672 с.
2. Эйзенберг Д., Кауцман В. Структура и свойства воды. Ленинград: Гидрометеоиздат, 1975. 280 с.
3. Берсукер И. Б. Эффект Яна - Теллера и вибронные взаимодействия в современной химии. Москва: Наука, 1987. 344 с.
4. Малафаев Н. Т. О взаимодействиях и динамике молекул в чистой воде // Восточно-европейский журнал передовых технологий. 2011. T. 52, № 4 (8). С. 48-58. URL: http://journals.uran.ua/eejet/article/view/1465/1363
5. Малафаев Н. Т., Погожих Н. И. Моделирование вращательных колебаний молекул воды // Восточно-европейский журнал передовых технологий. 2015. T. 2, № 5 (74). С. 27-35. doi: http://doi.org/10.15587/1729-4061.2015.40569
6. Маленков Г. Г. Структура и динамика жидкой воды // Журнал структурной химии. 2006. Т. 47. C. 5-35. doi: http://doi.org/10.1007/s10947-006-0375-8
7. Зельдович Б. Я., Суало М. Дж. Двухчастотный маятник на вращающейся платформе: моделирование оптических явлений // Успехи физических наук. 2004. Т. 174, № 12. С. 1337-1354. doi: http://doi.org/10.3367/ufnr.0174.200412e.1337
8. Пространственный потенциал однородного кругового тора через эквигравитирующие элементы / Кондратьев Б. П., Дубровский А. С., Трубицына Н. Г., Мухаметшина Э. Ш. // Журнал технической физики. 2008. Т. 78, № 7. С. 132-135.
9. Radiofrequency trapping of ions in a pure toroidal potential distribution / Higgs J. M., Petersen B. V., Lammert S. A., Warnick K. F., Austin D. E. // International Journal of Mass Spectrometry. 2016. Vol. 395. P. 20-26. doi: http://doi.org/10.1016/j.ijms.2015.11.009
10. Wu S.-T., Chen Y.-R., Wang S.-S. Two-degree-of-freedom rotational-pendulum vibration absorbers // Journal of Sound and Vibration. 2011. Vol. 330, Issue 6. P. 1052-1064. doi: http://doi.org/10.1016/jjsv.2010.09.028
11. Малафаев Н. Т. Вращательные колебания молекул воды как колебания сферического маятника в неоднородном поле сил // ScienceRise. 2016. Т. 2, № 2 (19). С. 62-69. doi: http://doi.org/10.15587/2313-8416.2016.60587
12. Water Models. URL: http://www.lsbu.ac.uk/water/models.html
13. Makhlaichuk P. V., Malomuzh M. P., Zhyganiuk I. V. Dimerization of water molecules modeling of the attractive part of the interparticle potential in the multipole approximation // Ukrainian Journal of Physics. 2013. Vol. 58, Issue 3. P. 278-288. doi: http://doi.org/10.15407/ujpe58.03.0278
Рекомендовано до публгкацИ д-р техн. наук Погожих М. I.
Дата надходження рукопису 10.06.2019
Малафаев Микола Тимофшович, кандидат фiзико-математичних наук, доцент, кафедра фiзико-математичних та шженерно-техшчних дисциплш, Харшвський державний ушверситет харчування та то-рпвл^ вул. Клочшвська, 333, м. Харшв, Украна, 61051 E-mail: mnt949@gmail.com