УДК 539.194: 544.273
DOI: 10.15587/2313-8416.2015.51842
СИЛОВЫЕ ОСОБЕННОСТИ КОЛЕБАНИЙ ДВУХЧАСТОТНОГО СФЕРИЧЕСКОГО МАЯТНИКА В НЕОДНОРОДНОМ ПОЛЕ СИЛ
© Н. Т. Малафаев
Проведен анализ вращательных двумерных колебаний молекул воды с помощью модели двухчастотного сферического маятника в неоднородном поле сил. Рассмотрена динамика частот при переходе от независимых двумерных вращательных колебаний к вращениям на одной частоте. Показано, что согласование частот маятника обусловлено наличием направленного неоднородного по углу поля сил. Проведен Фурье-анализ орбит
Ключевые слова: молекула воды, сферический маятник, неоднородное поле сил, фазовые траектории, Фурье-анализ
It is conducted an analysis of two-dimensional rotational oscillations of water molecules using a dual frequency model of a spherical pendulum in an inhomogeneous field offorces. It is considered a frequency dynamics in the transition from independent two-dimensional rotational oscillations to the rotation at the same frequency. It is shown that the pendulum frequency matching due to the presence of the field offorces that is directed and angularly inhomogeneous. Fourier-analysis of the orbits is done
Keywords: water molecule, spherical pendulum, inhomogeneous field of forces, phase trajectories, Fourier-analysis
1. Введение
Вода, ее свойства и подвижность ее молекул широко изучаются экспериментально, а также с помощью компьютерного моделирования [1-5]. Хотя эти результаты хорошо согласуются с имеющимися свойствами воды, многие вопросы физики воды до сих пор не изучены или требуют уточнения.
В настоящее время достигнуто понимание важности колебательных спектров воды для ее свойств. Однако основное внимание при этом уделяется проблеме коллективных движений соседних частиц [3-5], временам жизни кластеров в воде, а не их конкретным модам колебаний и их особенностям. До сих пор до конца не ясна основная особенность поведения молекул в жидкой воде - причина изгибов ее водородных связей (Н-связей), которые автор связывает с вращательными колебаниями молекул воды [6]. Наличие колебательных мод и их коллективной упорядоченности существенно влияет на термодинамические свойства воды. Однако колебания молекул воды, их моды и силы взаимодействия все еще остаются недостаточно изученной проблемой.
2. Анализ литературных данных и постановка проблемы
Свойства вещества зависят как от природы взаимодействий между его атомами и молекулами, так и от колебательных спектров этих атомов и молекул. Данные вопросы подробно изучаются в рамках теории эффекта Яна-Теллера (ЭЯТ) [7], рассматривающей совместные колебания электронно-ядерных подсистем, роль кооперативных колебаний и их энергии в свойствах и фазовых превращениях в системах атомов. Согласно теореме Яна-Теллера [7] в многоатомных системах линейные связи неустойчивы и с нагревом динамически изгибаются. Современная теория ЭЯТ рассматривает особенности колебаний атомов и молекул в системе, включая коопера-
тивные, их влияние на динамику ее свойств и структурные фазовые превращения в этой системе.
Проблеме изучения колебательных и, особенно, вращательных мод молекул воды до сих пор не уделено необходимого внимания, за исключением работы [6], где в рамках теории ЭЯТ было предложено рассматривать в качестве линейных связей воде Н-связи между ее молекулами. В жидкой фазе роль колебаний атомов и молекул должна возрасти. Было показано, что кооперативные (коллективные) колебания Н-свя-зей молекул воды в жидкой фазе относятся к классу либрационных (вращательных). Наличие большого числа коллективных колебательных мод приводит к наличию у воды больших величин теплоемкости, теплопроводности и подвижности ее молекул.
Данные колебания протонов около осей Н-связей в воде сопровождаются динамическими изгибами этих связей на некоторый средний угол в [2]. Наличие у молекулы двух протонов может приводить к их сложным колебаниям и вращениям вокруг осей Н-связей, которые можно описать как двумерные вращательные колебания вокруг собственных осей молекулы x, у, z [2, 6]. Вследствие несимметричности молекулы воды имеют три главных момента инерции ^ вдоль осей молекулы [2] и, соответственно, три частоты вращательных колебаний, что усложняет рассмотрение колебаний молекулы.
В работе [6], исходя из теории ЭЯТ, было предположено, что данные двумерные вращательные колебания являются эллиптическими, связанными, и сопровождаются вращениями протонов молекул воды вокруг осей Н-связей. Однако с точки зрения механики разные частоты колебаний в двух взаимно перпендикулярных направлениях являются независимыми и лишь при совпадении частот возможны вращения протонов вокруг осей Н-связей.
Поскольку Н-связи являются направленными, то это также необходимо учитывать при моделирова-
нии вращательных колебаний молекул воды с помощью модели маятника. В работе [8] было предложено моделировать направленное поле сил (НПС) тяжести для маятника по углу отклонения в от вертикали как F=m•G=m gcosnв. Было показано, что с ростом угла в изгиба Н-связи и показателя направленности поля сил п наблюдаются быстрое уменьшение скорости и частоты вращения маятника, что согласуется с данными для либрационных частот воды [2]. В рамках данной модели была оценена величина направленности поля сил в воде как показатель НПС п~8.
В работе [9] были промоделированы вращательные колебания молекул воды в НПС с помощью модели двухчастотного сферического маятника и изучена возможность перехода от двух вращательных независимых колебаний (НК) к единым двумерным вращениям молекул воды и ее протонов вокруг осей водородных связей. Показано, что главной причиною синхронизации частот является сильная зависимость частоты колебаний от угла изгиба Н-связи для случая направленных сил связи. Были определены области существования вращательных эллипсоподобных траекторий (ЭПТ). Показано, что наличие коллективных колебаний может почти вдвое уменьшить углы изгибов Н-связей, что хорошо согласуется с их измеренными величинами в воде [2].
Для модели двухчастотного сферического маятника [9] показано, что в отличие от колебаний обычного сферического маятника два независимых колебания привязаны к осям X и Y. Поэтому для независимых собственных частот колебаний маятника не наблюдается прецессии осей эллипсов по углу, а для ЭПТ на общей частоте эллипсы изменяются со временем и колебания происходят не в круговой, а в эллиптической области, связанной с осями X- Y.
Однако еще много вопросов динамики двухча-стотных вращательных колебаний молекулы воды в неоднородном поле сил, их силовые и частотные характеристики, в работе [9] не были проанализированы. Данное направление исследований является перспективным, поскольку позволяет получить качественно новые данные о влиянии направленного неоднородного поля сил взаимодействий на особенности собственных колебаний молекул воды.
3. Цель и задачи исследования
Целью работы является продолжение анализа [9] особенностей динамики вращательных колебаний молекулы воды в областях их независимых колебаний и общих двумерных эллипсоподобных вращений протонов молекул воды на примере модели двухча-стотного сферического маятника.
Для достижения поставленной цели были поставлены следующие задачи:
1. Изучение динамики изменения двухчастот-ных колебаний маятника по осям X и Y с ростом начальной скорости и перехода от независимых колебаний к единым двумерным эллипсоподобным вращениям в направленном поле сил.
2. Изучение динамики сил с ростом величины отклонений маятника в НПС, а также влияние их ди-
намики на изменение частот свободных колебаний. Построение фазовых траекторий для скоростей и ускорений (сил).
3. Проведение Фурье-анализа эллиптических орбит маятника. Определение особенностей спектров частотных характеристик этих траекторий.
4. Компьютерная модель вращательных колебаний
Модель двухчастотного сферического маятника подробно описана в [9] и представляет собой две равные точечные массы общей единичной массы, подвешенные в точке О (начало координат модели) на одинаковых невесомых стержнях единичной длины под углом в между ними в НПС. В равновесии массы сориентированы вдоль оси Y. Маятник имеет разные моменты инерции ^ и IX при их отношении k=IY/Ix >1 и совершает колебания вокруг вертикальной оси Z как вращательные колебания вокруг осей Y и X модели. Анализ колебаний проводился в НПС для характерных для воды отношений моментов инерции k= 1,5 и 3. Считаем, что центр масс маятника описывает траектории подобные для траекторий протонов молекул воды.
Локальное моделирование вращательных колебаний молекул воды было проведено в пакете MatLab с помощью модели двухчастотного сферического маятника, имеющего два момента инерции и, соответственно, две частоты во взаимно перпендикулярных направлениях в НПС вида G=gcosПв.
Было проведено моделирование свободных колебаний двухчастотного сферического маятника в неоднородных поле тяжести G•с•п=8, g=10 м/с2 для отношений моментов инерции k= 1,5 и 3, начального относительного смещения =0,2 центра масс по оси Y и его начальной скорости vхо по оси X. Время моделирования составляло 20...100 с при шаге 0,0001 с. С ростом начальной скорости и соответствующей амплитуды колебаний по оси Х наблюдается понижение более высокой частоты колебаний вдоль этой оси и переход к вращениям маятника вокруг вертикальной оси на общей частоте.
5. Моделирование колебаний двухчастотного сферического маятника
На рис. 1 представлены результаты моделирования траекторий движения двухчастотного сферического маятника в НПС в плоскости X-Y для отношения моментов инерции £=1,5 в зависимости от величин начальных скоростей V,.,,. Для отношения моментов £=3 ЭПТ наблюдаются в узком интервале скоростей [9] и потому маловероятны в воде, хотя особенности колебаний маятника при этом выражены наиболее ярко.
Для малых начальных скоростей маятника Vх наблюдаются вращательные независимые колебания маятника на двух частотах в прямоугольной области X-Y (рис. 1, а). При дальнейшем росте скорости границы этой области начинают деформироваться: горизонтальные границы области (ось Х) становятся выпуклыми, вертикальные (ось Y) -
вогнутыми, а углы вдоль оси Y попарно сближаются. При слиянии углов при критической скорости vкр происходит качественный переход от враща-
тельных НК к эллипсоподобным вращениям маятника на общей частоте вокруг вертикальной оси Z (рис. 1, б-г).
0.3 0.2 0.1 0 -0.1 -0.2
Y - X
-0.2 -0.1 0 0.1 0.2 0.3
а
0.6 Y - X
0.4
0.2 \
\
0
у\
-0.2 /
-0.4
-0.4 -0.2 0 0.2 0.4 0.6
-0.4 -0.2 0 0.2 0.4 0.6
б
Рис. 1. Проекции траекторий движения двухчастотного сферического маятника в координатахX- Y для £=1,5 в зависимости от величин начальных скоростей vхо равных: vхо: а - 0,4; б - 0,7; в - 0,757; г - 0,84 м/с
При возрастании величины начальной скорости выше критической растет минимальная амплитуда колебаний по оси Y(рис. 1, б-г). Видим, что при скорости vхо<vel траектории (рис. 1, б) вписываются в эллипс (рис. 1, в), а при скоростях vхо >vel (рис. 1, г) как бы описываются вокруг него, пока скорость не достигает максимального значения vmax, при котором маятник выходит за пределы нижней полусферы и колебания маятника срываются.
Можно отметить, что при vхо<vel траектория через четверть периода после начала движения пересекает ось X под углом с положительным тангенсом, большая ось эллипса траектории находится при у>0, а при vхо>vel - все наоборот. Для скорости vх очерез четверть периода ось X пересекается под прямым углом, и вся траектория является симметричной и эллиптической.
На рис. 2 показаны периоды колебаний двухча-стотного сферического маятника в неоднородном поле сил вдоль осей X и Y в зависимости от его начальных скорости и смещения в областях НК и ЭПТ (области показаны для смещения =0,1) для отношения £=1,5. Для 5уо=0,3 дополнительное смещение изменило частоты и скорости. В критической точке слияния частот период больше, чем их начальные значения при малых скоростях. После слияния в области существования ЭПТ период несколько уменьшается, а затем быстро возрастает при vmax, что особенно видно для смещения 5уо =0,3. Для случая £=3 различие частот большее, соответственно, критическая скорость выше и область ЭПТ более узкая [9], а периоды общих колебаний для нее быстрее возрастают.
С появлением общих ЭПТ можно говорить об общей фазе вращения радиус-вектора траектории
в
г
двухчастотного маятника в плоскости X-Y со временем: р=ат^(у/х). Для независимых колебаний можно рассматривать их 2 независимые фазы а и в для собственных колебаний по осям X и Y и общую фазу р. Легко увидеть, что для конического маятника временная зависимость () будет прямой линией. Для круговых участков траектории маятника должен наблюдаться линейный рост фазы, а для сильно эллиптичных либо для траекторий проходящих близко к оси Z маятника - ступенеподобные. Для сферического маятника зависимость будет ступенеподобной при росте фазы ( .
На рис. 3 показано изменение общей фазы вращения маятника р в плоскости X-Y от времени для случаев, показанных на рис. 1. Для независимых колебаний и малых скоростей видим кроме колебаний на собственных частотах также колебания общей фазы р (рис. 3, а). Для двухчастотного маятника с независимыми частотами по осям mx и mY наблюдаются вращения маятника с частотою биения амплитуд mX-mY в направлении и против вращения часовой стрелки, что и сопровождается колебаниями фазы р в некоторых пределах.
Рис. 2. Периоды колебаний Тху двухчастотного сферического маятника для отношения £=1,5 в зависимости от его начальных смещений: 5уо=0,1 (кривые Тх1, Ту1) и 5уо =0,3 ( кривые Тх3, Ту3), от начальной скорости Vхо для независимых колебаний вдоль осей X и Y и общих эллипсоподобных вращений вокруг оси Z
Рис. 3. Изменения фазы колебаний (двухчастотного сферического маятника с £=1,5 со временем для разных
начальных скоростей Vхо (см. рис. 1)
в
г
С ростом начальной скорости ух0 частоты НК и их периоды ТХ и Ту (рис. 2) вначале сближаются, что ведет к росту длительности периодов колебаний общей фазы р (рис. 3, а). При совпадении двух частот и их периодов ТХ и Ту колебания общей фазы р исчезают, и появляются общие эллипсоподобные вращения маятника (рис. 1, б—г), а зависимость фазы от времени становится ступенеподобной при общем линейном росте фазы (рис. 3, б-г). При скорости Vx0=Vel орбиты становятся эллиптическими, а все ступени на зависимости фазы р) становятся одинаковыми (рис. 3, в) и зависят от эллиптичности данной орбиты. При дальнейшем росте скорости быстро возрастают периоды колебаний и начинают преобладать участки с горизонтальными ступенями (рис. 3, г), а также с круговыми вращениями (рис. 1, г и 3, г).
Можно отметить, что модуль средней угловой скорости изменения фазы для обоих направлений сохраняется (рис. 3, а) и осредненные кривые кусочно линейные. Биения общей фазы р связаны с зависимостью разности фаз колебаний а и в по осям X и У - в точках изменения направления вращений фазы наблюдается либо пучность, либо узел биений разности фаз а - в.
Прекращение биений общей фазы р и появление ее постоянной линейной составляющей (рис. 3, б-г) является главным параметром, свидетельствующем о наличии перехода между независимыми колебаниями и эллипсоподобными вращениями и позволяет проследить динамику двухчастотных колебаний маятника в зависимости от величин его параметров: п, к, дуо и ух0.
Анализ влияния НПС
Малое изменение фазы рсо временем (рис. 3, г) может говорить об ослаблении сил при больших отклонениях маятника и появлении ангармонизма. На рис. 4 показаны временные зависимости координат, скоростей и ускорений вдоль осей X и У. Видим, что для координат синусоиды слабо искажены. Но для скоростей видим эффекты заострения синусоид по оси X и закругления по оси У, а для ускорений наоборот - видим провалы в максимумах по оси X и заострения по оси У. Данные изменения можно трактовать как появление значительных вкладов от третьей и других нечетных гармоник колебаний.
При анализе колебаний были рассмотрены траектории колебаний методом фазовой плоскости [10]. Обычно фазовые траектории рассматриваются как зависимости скорости от координаты. Для одномерных свободных колебаний при отсутствии затуханий они представляют собой эллипсы.
Для двумерных колебаний рассмотрим фазовые траектории раздельно по обеим осям (рис. 5). С ростом скорости и смещений вдоль оси X кривые изменяются и в направлении оси X становятся ром-боподобными с практически линейным уменьшением скорости, а для смещений вдоль оси У прямоуголь-ноподобными. Наибольшие изменения эллипсов наблюдаются для отношения моментов к=3. Вблизи границы срыва колебаний для фазовой траектории по оси X наблюдаются еще большее «заострение» диа-
граммы вдоль оси X и перегибы кривых Данные зависимости (рис. 4 и 5) свидетельствуют о наличии значительной ангармоничности колебаний маятника в НПС.
-0.2
\ Я /'"'ч / * ■/ \ / . \ / \
\ \ / ' /\ \ / /\ \ / \ / / \ / \.........X.......:7........я t"\.......А!
К.у \ К.У \ к;
V у
б
Рис. 4. Временные зависимости для координат (—), скоростей (—) и ускорений (•••) эллиптической орбиты вдоль осей: а - X; б - У для к=1,5, ух0 =0,757 м/с
Причиною ангармоничности колебаний маятника является существенная нелинейность между силою и отклонениями при наличии НПС. Зависимость скатывающей силы в НПС описываются как F=mgsinвcosnв и зависит от значений показателя п. При малых углах в сила пропорциональна смешению и колебания являются гармоничными. С ростом отклонений маятника при п>0 возможно появление максимума силы и затем ее уменьшение. Это объясняет наличие провалов для ускорений по оси X (рис. 4, а), малое изменение при этом фазы р со временем (рис. 3, г), форму фазовых траекторий (рис. 5).
На рис. 6 показаны фазовые траектории для ускорений (сил) ах - х и ау - у для эллиптических орбит. Для неэллиптических траекторий кривые более сложные во времени, хотя и с аналогичными особенностями. Видим значительное уменьшение сил для больших смещений вдоль оси X, особенно для отношения к=3 с большим различием собственных частот (рис. 6, б).
Вследствие центральной симметричности кривых (рис. 6) мы можем говорить об их описании полиномами с нечетными степенями Хт, что должно приводить к появлению т - х (нечетных) гармоник колебаний $1п(тв). Наибольшие вклады гармоник
а
должны наблюдаться для ускорений и соответствующих сил. Для скоростей и координат, вследствие
интегрирования, вклады гармоник должны соответ-
2
ственно уменьшаться в т и т раз, что подтверждается расчетами для траекторий. Сами эллиптические орбиты (рис. 1, в) содержат мало гармоник (см. далее Фурье-анализ траекторий).
У^, м/с
-0.5 0 0.5 х. у. М
Рис. 5. Фазовые траектории эллиптической орбиты для осей X (внешняя кривая) и Y (внутренняя) для случая £=3, vxо =1,222 м/с
аху, м'с2 I
-0.4 -02 0 02 04 06
У, м
а
аку, м'с2
х-ах + у-ау(г)
г 1.5 1
0.5 О
-0.5 -1 -1.5 -2
-1 -0.5 0 0.5 1
м
б
Рис. 6. Силовые фазовые траектории: кривые
1 - зависимость ах - х, кривые 2 - ау - у: а - для
£=1,5, vxо =0,757 м/с; б - для £=3, vxо = 1,222 м/с
Видим на рис. 6, что силовые фазовые траектории вдоль оси Y испытывают как бы большие ускорения с ростом отклонения. На самом деле эти ускорения нормальные для больших отклонений по оси Y (X мало), хотя и уменьшены в к раз, сравнительно с ускорениями по оси X! Но они значительно уменьшаются вблизи оси Y, когда наблюдаются наибольшие отклонения маятника вдоль оси X и, соответственно, углов в, и имеем общие минимальные ускорения при этом. Видим, что для к=3 уменьшение сил по осям более значительно, чем для к=1,5. Для эллиптических орбит относительное уменьшение силы по оси X сравнительно с ее максимумом для отношения к=1,5 составляет около 0,9 и для отношения к=3 - около 0,3 и менее (рис. 6).
Можно отметить, что вследствие различия моментов инерции маятника скатывающее ускорение по оси Y уменьшается в к раз и результирующее ускорение маятника относительно оси Z для промежуточных углов фаз (<рфп-п/2) становится не центральным, а вращающим, что способствует появлению вращений и знакопеременной прецессии маятника. Хотя эта эффективная сила знакопеременная, ее роль для ЭПТ является значительной.
Таким образом, главной причиной торможения колебаний маятника и уменьшения его частоты с ростом амплитуды отклонений является быстрое уменьшение скатывающих сил в НПС. Это приводит к длительному времени пребывания маятника в дальних точках траектории и сближению частот, «закруглению» углов прямоугольных областей для НК вращающими моментами сил и, в итоге, к переходу двухчастотных независимых колебаний к единым эллипсоподобным вращениям маятника на одной общей частоте.
Фурье-анализ траекторий
Был проведен Фурье-анализ эллиптических X-Y орбит с помощью быстрого преобразования Фурье FFT. Разложения по осям X и Y выполнялись соответственно в ряды нечетных гармоник по синусам и косинусам.
Для ЭПТ времена обращений маятника вокруг оси Z зависят от его текущих максимальных отклонений и все время изменяются, что ведет к большому числу локальных частот, изменению частот и их фаз с ростом начальной скорости. Для независимых колебаний частоты более стабильны, но также начинают сильно изменяться при приближении начальной скорости к критической. Поэтому имеет смысл рассматривать только эллиптические орбиты, для которых наблюдается строгая периодичность и повторяемость траекторий. Для неэллиптических траекторий силовые особенности колебаний сохраняются, хотя и дополнительно усложняются вследствие постоянных изменений как угла в, так и разности фаз колебаний по осям X и Y.
Для времени моделирования колебаний в 20...50 с и при кратности его периоду колебаний, с относительной точностью порядка 10-4 и лучше, наблюдалось хорошее разрешение частот. Среди первой десятки наибольших амплитуд в быстром
преобразовании Фурье наблюдались гармоники для к=3 до т=15 и для £=1,5 - до т=9, амплитуда для последней гармоники составляла менее 0,1 % от первой. Установлено, что все амплитуды для данных гармоник являются положительными. На рис. 7 показана гистограмма для амплитуд первых 5 нечетных гармоник (т=1...9) относительно первой для ускорений маятника вдоль осей X и У эллиптических орбит для отношения моментов £=3 и показателей НПС п=2, 4 и 8.
A(m)/A(1)
1,2
1
0,8 0,6 0,4 0,2 0
т = 1 3 5 7 9
□ 2 х □ 4 х □ 8 х □ 2 у □ 4 у П8 у
Рис. 7. Относительные вклады гармоник для ускорений на эллиптической орбите с £=3, ухо = 1,222 м/с от их номера т для показателей НПС п=2, 4 и 8 вдоль осей X (однородная окраска столбцов) и У (штрихованные столбцы)
Видим, что роль гармоники обычно растет с
ростом показателя НПС п. Для координат (орбит)
2
вклады от гармоник в т раз меньшие и они соответственно значительно быстрее убывают с ростом номера гармоники т. Гистограмма вкладов гармоник показывает (рис. 7), что при £=3 для траекторий (координат) вклады гармоник по оси X будут менее 10 % и 1,5 %, соответственно, для третьей и пятой гармоник. Видим, что относительные абсолютные вклады гармоник по оси X примерно вдвое большие по величине, чем по оси У. Для отношения £=1,5 вклады гармоник примерно вдвое меньшие сравнительно с отношением £=3 и быстрее убывают с ростом их номера т.
Необходимо отметить, что в данном Фурье -анализе X-Y орбит маятника пренебрегается их изменениями вдоль оси 2.
6. Обсуждение результатов моделирования двухчастотного сферического маятника в неоднородных полях сил
Результаты моделирования двухчастотного сферического маятника в неоднородных полях сил показывают, что для объяснения особенностей его колебаний служит его силовая фазовая траектория, показывающая значительное ослабление результирующих сил при максимальных отклонениях маятника.
Медленное вращение молекул воды при максимальных отклонениях в НПС также приводит к длительному их пребыванию в крайних точках при быстром процессе их движения между ними. В ре-
зультате понижается частота, а средние углы изгибов Н-связей будут близкими к максимальным.
Для несвободных колебаний молекул в среде других молекул их скорости будут сильно промоду-лированы взаимодействиями с соседями, что будет давать значительный вклад в локальное изменение средней кинетической энергии молекулы и ее результирующую траекторию. Поскольку для областей существования ЭПТ с отношением £=3 диапазон изменений кинетической энергии не превышает 10 % [9], то маловероятно совершение даже одного вращательного колебания для молекул. Для отношения £=1,5 дополнительный рост кинетической энергии для областей существования ЭПТ составляет 70. 100 %, а потому для него такие вращения более вероятны.
Для случая синхронных и кооперативных вращений молекул их оси X и У вращаются также и становятся динамическими. Относительно этих осей отклонения и скорости для соседних молекул малы, а в идеале - нулевые. Поэтому силы между ними уменьшаются слабо. Такую группу сильно связанных соседних молекул и называют кластером. Относительно соседних молекул внутренняя молекула кластера обладает малой кинетической энергией, а, следовательно, она отдает свою исходную кинетическую энергию соседям. Поэтому вначале кластер обладает повышенной кинетической энергией относительно других «свободных» молекул. После релаксационного перераспределения энергии кластер с коллективными колебаниями будет иметь меньшую энергию и его состояние будет энергетически более стабильным.
Таким образом, в воде становится возможным существование кластеров вблизи точки плавления. Их вклад в энергию активации воды экспоненциально уменьшается с нагревом [6], одновременно уменьшается энергия Н-связей, а сами кластеры разрушаются [2]. Если для линейных систем более-менее понятно, что такое кооперативные синхронные колебания молекул, то для водных кластеров, состоящих из гексагонально упорядоченных объемных слоев молекул, это пока совсем непонятно и требует анализа.
7. Выводы
Проанализирована модель колебаний двухча-стотного сферического маятника в неоднородном поле сил и ее особенности и установлено:
1. Важнейшей характеристикой двухчастотно-го сферического маятника для объяснения его аномалий и особенностей колебаний в неоднородном поле сил является его силовая фазовая траектория.
2. Главной причиной торможения колебаний маятника является быстрое уменьшение возвращающих сил с ростом амплитуды отклонений и наличие максимума НПС, что приводит к долгому времени пребывания маятника в дальних точках траектории и понижению частот колебаний.
3. Проведен Фурье-анализ для эллиптических орбит, имеющих четкие частоты. Расчеты проведены для ускорений маятника, имеющих наибольшие вклады гармоник и показывающих рост вкладов гар-
моник с ростом скорости, угла отклонения, отношения моментов к и показателя п.
Данные особенности колебаний двухчастотно-го сферического маятника в неоднородном поле сил согласуются с особенностями колебаний молекул воды, что указывает на необходимость учета моментов инерции и неоднородного направленного поля сил в расчетах потенциалов для молекул воды при моделировании их колебаний методом молекулярной динамики.
Литература
1. Антонченко, В. Я. Основы физики воды [Текст] /
B. Я. Антонченко, А. С. Давыдов, В. В. Ильин. - К.: Наук. Думка, 1991. - 672 с.
2. Эйзенберг, Д. Структура и свойства воды [Текст] / Д. Эйзенберг, В. Кауцман; пер. с англ. - Л.: Гидрометеоиз-дат, 1975. - 280 с.
3. Маленков, Г. Г. Структура и динамика жидкой воды [Текст] / Г. Г. Маленков // ЖСХ. - 2006. - Т. 47. -
C. 5-35.
4. Kumar, P. Molecular dynamics study of orienta-tional cooperativity in water [Text] / P. Kumar, G. Franzese, S. V. Buldyrev, H. E. Stanley // Physical Review E. -2006. - Vol. 73, Issue 4. - P. 041505. doi: 10.1103/ physreve .73.041505
5. Voloshin, V. P. Collective motions in computer models of water. large-scale and long-time correlations [Text] / V. P. Voloshin, G. G. Malenkov, Yu. I. Naberukhin // Journal of Structural Chemistry. - 2013. - Vol. 54, Issue S2. - P. 233251. doi: 10.1134/s0022476613080052
6. Малафаев, Н. Т. О взаимодействиях и динамике молекул в чистой воде [Текст] / Н. Т. Малафаев // Восточно-Европейский журнал передовых технологий. - 2011. -T. 4, № 8 (52). - С. 48-58. - Режим доступа: http:// journals.uran.ua/eejet/article/view/1465/1363
7. Берсукер, И. Б. Эффект Яна - Теллера и виброн-ные взаимодействия в современной химии [Текст] / И. Б. Берсукер; пер. с англ. - М.: Наука, 1987. - 344 с.
8. Малафаев, М. Т. Обертання молекул води як рух сферичного маятника в неоднорщному rn^i сил [Текст] /
М. Т. Малафаев // Прогресивш технжа та технологи харчо-вих виробництв ресторанного господарства i торпвлг -2014. - Вип. 1.- С. 291-298.
9. Малафаев, Н. Т. Моделирование вращательных колебаний молекул воды [Текст] / Н. Т. Малафаев, Н. И. Погожих // Восточно-Европейский журнал передовых технологий. - 2015. - Т. 2, № 5 (74). - С. 27-35. doi: 10.15587/ 1729-4061.2015.40569
10. Марчуков, Б. А. Проектирование систем управления методами фазовой плоскости [Текст] / Б. А. Марчуков. - М.: Машиностроение, 1976. - 248 с.
Refere^es
1. Antonchenko, V. Ya, Davydov, А. S, Iliin, V. V. (1991). Osnovy fizyky vody. Kyiv: Naukova Dumka, 672.
2. Eisenberg, D., Kauzmann, W. (2005). The structure and properties of water. Oxford University press, 308.
3. Malenkov, G. G. (2006). Structure and dynamics of liquid water /Journal structural chemistry, 47, 5-35.
4. Kumar, P., Franzese, G., Buldyrev, S. V., Stanley, H. E. (2006). Molecular dynamics study of orientational cooperativi-ty in water. Physical Review E, 73 (4), 041505. doi: 10.1103/ physreve.73.041505
5. Voloshin, V. P., Malenkov, G. G., Naberukhin, Yu. I. (2013). Collective motions in computer models of water. large-scale and long-time correlations. Journal of Structural Chemistry, 54 (S2), 233-251. doi: 10.1134/s0022476613080052
6. Маlаfаеv, N. Т. (2011). About interactions and dynamics of molecules in pure water. Eastern-European Journal of Enterprise Technologies, 4/8 (52), 48-58. Available at: http://journals.uran.ua/eejet/article/view/1465/1363
7. Bersuker, I. (1984). The Jahn-Teller Effect and Vibronic Interactions in Modern Chemistry. Plenum. New York, 319.
8. Маlаfаеv, M. Т. (2014). Obertannya molekul vodi yak ruh sferichnogo mayatnika v neodnoridnomu poli syl., 1, 291-298.
9. Malafayev, N., Pogozhikh, N. (2015). Features rotational of vibrations of water molecules. Eastern-European Journal of Enterprise Technologies, 2/5 (74), 27-35. doi: 10.15587/ 1729-4061.2015.40569
10. Marchukov, A. B. (1976). Proektirovanie sistem up-ravleniya metodami fazovoy ploskosti. Moscow: Mashi-nostroenie, 248.
Рекомендовано до публгкацИ д-р техн. наук Погожих М. I.
Дата надходження рукопису 21.09.2015
Малафаев Николай Тимофеевич, кандидат физико-математических наук, доцент, кафедра физико-математических и инженерно-технических дисциплин, Харьковский государственный университет питания и торговли, ул. Клочковская, 333, г. Харьков, Украина, 61051 E-mail: [email protected]