Математический сборник
RUS  ENG    ЖУРНАЛЫ   ПЕРСОНАЛИИ   ОРГАНИЗАЦИИ   КОНФЕРЕНЦИИ   СЕМИНАРЫ   ВИДЕОТЕКА   ПАКЕТ AMSBIB  
Общая информация
Последний выпуск
Скоро в журнале
Архив
Импакт-фактор
Правила для авторов
Лицензионный договор
Загрузить рукопись
Историческая справка

Поиск публикаций
Поиск ссылок

RSS
Последний выпуск
Текущие выпуски
Архивные выпуски
Что такое RSS



Матем. сб.:
Год:
Том:
Выпуск:
Страница:
Найти






Персональный вход:
Логин:
Пароль:
Запомнить пароль
Войти
Забыли пароль?
Регистрация


Математический сборник, 2021, том 212, номер 2, страницы 81–105
DOI: https://doi.org/10.4213/sm9418
(Mi sm9418)
 

Эта публикация цитируется в 10 научных статьях (всего в 10 статьях)

Топологический анализ биллиарда, ограниченного софокусными квадриками, в потенциальном поле

С. Е. Пустовойтов

Механико-математический факультет, Московский государственный университет имени М. В. Ломоносова
Список литературы:
Аннотация: Рассмотрим биллиард в плоской области, ограниченной софокусными эллипсами и гиперболами. На материальную точку действует гуковский потенциал. Оказывается, эта динамическая система вполне интегрируема по Лиувиллю. В работе проведен топологический анализ слоения Лиувилля изоэнергетических многообразий всевозможных уровней гамильтониана и построены их полные инварианты Фоменко–Цишанга (меченые молекулы).
Библиография: 14 названий.
Ключевые слова: потенциал Гука, интегрируемая система, инвариант Фоменко–Цишанга, лиувиллева эквивалентность.
Финансовая поддержка Номер гранта
Российский научный фонд 17-11-01303
Исследование выполнено за счет гранта Российского научного фонда (проект № 17-11-01303).
Поступила в редакцию: 28.03.2020
Англоязычная версия:
Sbornik: Mathematics, 2021, Volume 212, Issue 2, Pages 211–233
DOI: https://doi.org/10.1070/SM9418
Реферативные базы данных:
Тип публикации: Статья
УДК: 517.938.5
MSC: 37C83, 37J35

§ 1. Введение

Математический биллиард – это динамическая система, описывающая движение материальной точки в компактной ограниченной области с кусочно гладкой границей со стандартным отражением вектора скорости на границе. Эта классическая задача динамики рассматривалась и рассматривается многими математиками с разных точек зрения: вычисление траекторий, изучение эргодичности, задачи оптических систем, а также с топологической точки зрения. Изучением слоений Лиувилля различных биллиардов занимаются А. Т. Фоменко, В. В. Ведюшкина и их ученики (см. [1]–[8]), М. Раднович, В. Драгович (см. [9]) и многие другие. Наша предыдущая работа [10] была посвящена биллиарду в эллиптическом кольце с гуковским потенциалом. Настоящая работа посвящена описанию слоения Лиувилля биллиардов во всех плоских областях, ограниченных софокусными эллипсами и гиперболами, и таким образом является обобщением предыдущей.

§ 2. Необходимые сведения

Напомним основные определения и теоремы, связанные с теорией гамильтоновых интегрируемых систем. Более подробно о ней см. [11].

Определение 1. Гладкое многообразие $M^4$ называется симплектическим, если на нем задана симплектическая структура: невырожденная, замкнутая 2-форма $\omega$.

Определение 2. Динамическая система на симплектическом многообразии $(M^4,\omega)$ называется гамильтоновой, если на $M^4$ можно выбрать функцию $H$ так, что динамическая система запишется в виде $v=\operatorname{sgrad}H$ ($=\omega \operatorname{grad}H$); $H$ называется гамильтонианом.

Определение 3. Если у гамильтоновой системы $v$ на $M^{4}$ есть две функции $f_1$, $f_2$, которые являются первыми интегралами, и верно:

1) $\{f_1, f_2\}=0$,

2) $f_1$ и $f_2$, функционально независимы,

3) векторные поля $\operatorname{sgrad}(f_i)$ полны,

то такая система называется вполне интегрируемой по Лиувиллю.

Определение 4. Отображением момента называется отображение $F$: $M^{4} \to \mathbb R^2$, где $F(x)=(f_1(x), f_2(x))$. Бифуркационной диаграммой $\Sigma(F)$ называется образ критических точек отображения момента на $\mathbb R^2$. При этом точка называется критической, если в ней падает ранг дифференциала отображения момента. Иначе точка называется регулярной.

Теорема 1 (Лиувилля). Пусть на $M^{4}$ задана вполне интегрируемая гамильтонова система с гамильтонианом $H$ и $T_\xi$-регулярная, компактная, связная компонента поверхности уровня интегралов $f_1$, $f_2$ (прообраз регулярного значения $F$). Можно считать, что $f_1=H$. Тогда справедливы следующие утверждения.

1. $T_\xi$ диффеоморфно тору $T^2$. Этот тор называется тором Лиувилля.

2. Слоение Лиувилля в некоторой окрестности тора Лиувилля диффеоморфно прямому произведению тора на диск $D^2$.

Введем обозначение: $Q^3=\{x\in M^4\colon H(x)=\mathrm{const}\}$ – изоэнергетическое многообразие.

Определение 5. Системы $v_1$ и $v_2$ на $Q_1^3$ и $Q_2^3$ соответственно называются лиувиллево эквивалентными, если существует послойный диффеоморфизм из $Q_1^3$ на $Q_2^3$, сохраняющий ориентацию как самих многообразий, так и ориентацию критических окружностей.

Определение 6. 3-атомом называется прообраз критического значения отображения момента $F$ вместе со своей окрестностью в $Q^3$. Все $3$-атомы были классифицированы А. Т. Фоменко, А. В. Болсиновым, А. А. Ошемковым в [11].

Определение 7. Грубой молекулой называется граф, точки ребер которого представляют собой регулярные слои многообразия $Q^3$ (т.е. торы), а вершинам сопоставляются атомы, через которые происходят перестройки торов.

Определение 8. Меченой молекулой называется грубая молекула, ребрам и вершинам которых сопоставляются числа, описывающие, как именно переходит тор Лиувилля из окрестности одного атома в тор Лиувилля окрестности другого атома по соединяющему их ребру. Более подробно см. [11].

Теорема 2 (Фоменко–Цишанга). Системы $v_1$ и $v_2$ на $Q_1^3$ и $Q_2^3$ лиувиллево эквивалентны тогда и только тогда, когда их меченые молекулы совпадают.

§ 3. Постановка задачи

Рассмотрим биллиард в ограниченной плоской области, границы которой заданы дугами кривых вида

$$ \begin{equation} \frac{x^2}{a+\lambda}+\frac{y^2}{b+\lambda}=1, \end{equation} \tag{3.1} $$
т.е. дугами софокусных квадрик. Здесь $\lambda$ – параметр квадрики. При $\lambda>-a$ квадрика – эллипс, при $-b<\lambda<-a$ – гипербола. Дуги выбраны так, что все ее углы между ними равны ${\pi}/{2}$ (рис. 1). Такие области были рассмотрены и классифицированы В. В. Ведюшкиной в [12]. В ее работе были рассмотрены биллиарды, для которых на материальную точку не действуют внешние силы.

Определение 9. Простейшей элементарной (плоской) областью $\Omega$ назовем двумерное связное, компактное, плоское гладкое риманово многообразие с кусочно гладким краем, которое имеет изометричное вложение в плоскость, причем граница его образа при этом вложении состоит из сегментов софокусных квадрик семейства (формула (3.1)), углы между которыми не превышают $\pi$.

Определение 10. Простейшая элементарная область $\Omega_1$ называется эквивалентной другой простейшей элементарной области $\Omega_2$, если $\Omega_2$ можно получить из $\Omega_1$ путем композиции следующих преобразований:

1) последовательным изменением параметра $\lambda_i$ каждой границы при неизменном параметре остальных, так, чтобы это значение параметра $\lambda_i$ не принимало значения $-b$; иными словами, сегменты во время деформации остаются либо эллиптическими (т.е. параметры квадрик, на которых располагаются эти сегменты, непрерывно меняются в пределах $(-b, \infty)$), либо гиперболическими (т.е. параметры квадрик, на которых располагаются эти сегменты, непрерывно меняются в пределах $(-a, -b]$), либо все время лежат на фокальной прямой (во все время деформации параметр остается равным $-b$);

2) симметрией относительно координатных осей.

Теорема 3 (Ведюшкиной). Всего классов простейших элементарных областей $12$ (см. табл. 1).

Таблица 1.Биллиардные области

Биллиардная областьМолекулы при $k>0$Молекулы при $k<0$
112, 178, 11, 12, 15, 17, 18
215, 171, 15, 16, 17
31, 131, 13
41, 131, 2, 3, 5, 13
5171, 13, 17
65, 131, 2, 3, 5, 13
71717
81, 131, 4, 5, 17
91, 171, 13
101, 171, 6, 7, 17
118, 178, 17
128, 175, 8, 9, 10, 17

В настоящей статье рассмотрим новую задачу. А именно, к биллиардам, изучавшимся В. В. Ведюшкиной, добавим действие центрального гуковского потенциала $P=({k}/{2})(x^2+y^2)$. Ранее биллиард в эллиптической области с гуковским потенциалом был рассмотрен И. Ф. Кобцевым, а в кольцевой области – автором в [10].

Траектория материальной точки между ударами описывается системой уравнений

$$ \begin{equation} \begin{cases} \ddot{x} =-kx, \\ \ddot{y} =-ky. \end{cases} \end{equation} \tag{3.2} $$

Заметим, что при $k>0$ на точку действует сила, направленная к центру координат, при этом траекторией является дуга эллипса, а при $k<0$ сила отталкивает точку от начала координат, и соответствующая траектория является дугой гиперболы (в случае В. В. Ведюшкиной траекториями являются прямые линии).

Рассмотрим симплектическое многообразие $(M^4, \omega)$, где $M^4=\{(x, y, \dot{x}, \dot{y})$: $(x, y) \in \Omega,\, (\dot{x}, \dot{y})\in T\Omega\}$ – фазовое пространство системы. Зададим симплектическую структуру так:

$$ \begin{equation} \omega =\begin{pmatrix} 0 & 0 & 1 & 0\\ 0 & 0 & 0 & 1\\ -1 & 0 & 0 & 0\\ 0 & -1 & 0 & 0 \end{pmatrix}. \end{equation} \tag{3.3} $$

Определенная выше биллиардная система является гамильтоновой системой на $M^4$ c гамильтонианом $H=(\dot{x}^2+\dot{y}^2)/2+P$, который также является первым интегралом системы. Задача настоящей работы – найти второй первый интеграл, удовлетворяющий условию теоремы Лиувилля, и для каждого уровня интеграла $H$ построить полный инвариант Фоменко–Цишанга, т.е. меченую молекулу. Эта задача полностью решена (см. § 6).

§ 4. Разделение переменных и интегрируемость

Перейдем к эллиптическим координатам. Их связь с декартовыми координатами определяется системой

$$ \begin{equation} \begin{cases} x^2=\dfrac{(a+\lambda_1)(a+\lambda_2)}{a-b}, \\ y^2=\dfrac{(b+\lambda_1)(b+\lambda_2)}{b-a}. \end{cases} \end{equation} \tag{4.1} $$
Такой переход связан с тем, что в этих координатах дуги границы биллиардной области являются координатными кривыми. Гамильтониан в этих координатах имеет вид
$$ \begin{equation} H=\frac{2(a+\lambda_1)(b+\lambda_1)}{\lambda_1-\lambda_2}\mu^2_1+\frac{2(a+\lambda_2)(b+\lambda_2)}{\lambda_2-\lambda_1}\mu^2_2 + \frac{k}{2}(\lambda_1+\lambda_2), \end{equation} \tag{4.2} $$
где $\mu_1$ и $\mu_2$ – импульсы, соответствующие координатам $\lambda_1$ и $\lambda_2$. Оказывается, существует еще один первый интеграл. Он был получен методом, описанным В. В. Козловым в [13].

Теорема 4. Рассматриваемая биллиардная система допускает первый интеграл

$$ \begin{equation} F=2(a+\lambda_1)(b+\lambda_1)\mu^2_1+\frac{k}{2}\lambda_1^2-H\lambda_1. \end{equation} \tag{4.3} $$

Доказательство. Сначала разберем свободное движение точки между отражениями от стенок. Дифференцируя $F$ по времени $t$, получим
$$ \begin{equation*} \frac{dF}{dt}=2\dot\lambda_1(b+\lambda_1)\mu_1^2+2\dot\lambda_1(a+\lambda_1)\mu_1^2 +4(a+\lambda_1)(b+\lambda_1)\mu_1\dot\mu_1+k\lambda_1\dot\lambda_1-H\dot\lambda_1. \end{equation*} \notag $$
Подставим значения $\dot\lambda_1$ и $\dot\mu_1$ из уравнений Гамильтона
$$ \begin{equation} \begin{cases} \dot\lambda_i=\dfrac{\partial H}{\partial \mu_i}, \\ \dot\mu_i=-\dfrac{\partial H}{\partial \lambda_i}. \end{cases} \end{equation} \tag{4.4} $$
Можно убедиться, что после подстановки получим тождественный нуль. А значит, $F=\mathrm{const}$.

Теперь докажем, что функция $F$ сохраняется при отражении от любой стенки биллиарда, являющейся дугой квадрики указанного семейства (3.1). При отражении значение интеграла $H$ и координаты $\lambda_1$ не меняются. Поэтому достаточно проверить лишь неизменность величины $\mu_1^2$. Отметим, что закон отражения от софокусной квадрики в эллиптических координатах запишется так: если точка отражается от гиперболы, то вектор скорости изменится по правилу $(\dot\lambda_1, \dot\lambda_2)\to(-\dot\lambda_1, \dot\lambda_2)$, а если точка отражается от эллипса, то вектор скорости изменится по правилу $(\dot\lambda_1, \dot\lambda_2)\to(\dot\lambda_1, -\dot\lambda_2)$. Из уравнений Гамильтона вытекает формула

$$ \begin{equation} \dot\lambda_1=\frac{4(a+\lambda_1)(b+\lambda_1)}{\lambda_1-\lambda_2}\mu_1. \end{equation} \tag{4.5} $$
Следовательно, вектор импульса $(\mu_1, \mu_2)$ будет меняться при отражении аналогично вектору скорости. Поэтому значение $\mu_1^2$ не меняется. Теорема доказана.

Зафиксируем уровни интегралов $H=h$, $F=f$.

Теорема 5. Справедлива следующая формула, связывающая компоненты вектора скорости точки с ее координатами (формула разделения переменных):

$$ \begin{equation} \dot\lambda_i=\pm\frac{\sqrt{8}}{\lambda_2-\lambda_1}\sqrt{V(\lambda_i)}, \end{equation} \tag{4.6} $$
где $ V(z)=(-({k}/{2})z^2+hz+f)(a+z)(b+z)$ – многочлен четвертого порядка.

Гамильтоновы системы, допускающие уравнения подобного вида, были подробно исследованы М. П. Харламовым в [14]. Для доказательства теоремы 5 достаточно выразить импульсы $\mu_1$ и $\mu_2$ из формулы интеграла $F$ (4.3) и подставить их в уравнения Гамильтона (4.6).

§ 5. Области возможного движения и бифуркационные диаграммы

Определение 11. Областью возможного (допустимого) движения называется замкнутая область на биллиарде, в которой материальная точка может находиться при фиксированных значениях интегралов $h$ и $f$.

Отметим, что области возможного движения являются проекциями торов Лиувилля на биллиардный стол (это будет доказано в § 7).

Лемма 1. При фиксированных значениях $h$ и $f$ область, задаваемая системой неравенств $V(\lambda_i)\geqslant0$, является областью возможного движения.

Доказательство. Отметим следующее: если материальная точка при фиксированных значениях $h$ и $f$ может перемещаться через заданную точку $(\lambda_1, \lambda_2)$, то вследствие формулы (4.6) для этих координат $V(\lambda_i)\,{\geqslant}\,0$. И наоборот, если $V(\lambda_i)\,{\geqslant}\,0$, то, положив начальные условия $\dot\lambda_i\,{=}\,\pm\sqrt{8}/(\lambda_2-\lambda_1)\sqrt{V(\lambda_i)}$, получим движение через эту точку. Лемма доказана.

Замечание 1. Отметим, что если координаты $(\lambda_1, \lambda_2)$ не являются корнями многочлена $V$, то в этой точке вектор скорости принимает четыре значения, если одна из координат является корнем, то два значения, а если обе координаты являются корнями, то вектор скорости нулевой (рис. 2). Также отметим, что области возможного движения не всегда являются двумерными. Они могут быть и одномерными (например, движение по фокальной прямой), и нульмерными (неподвижная точка).

Лемма 2. Все корни многочлена $V$ вещественные (возможно, кратные).

Доказательство. Заметим, что $(a+\lambda_2)(b+\lambda_2)>0$, значит, при $k>0$ многочлен $-({k}/{2})z^2+hz+f$ обязан иметь два корня, иначе для любых значений $\lambda_2$$V(\lambda_2)<0$. Аналогично, $(a+\lambda_1)(b+\lambda_1)<0$, значит, при $k<0$ многочлен $-({k}/{2})z^2+hz+f$ также обязан иметь два корня. Лемма доказана.

Обозначим два корня многочлена $-({k}/{2})z^2+hz+f$ через $\xi$ и $\eta$. Эти корни, вообще говоря, являются функциями от переменных $(h,f)$, и именно они определяют область возможного движения.

Пример 1. Рассмотрим биллиард в области, ограниченной софокусными эллипсами с параметрами 0 и $\lambda_0$, параметр $k>0$, и корни многочлена $V$ удовлетворяют следующей цепочке неравенств: $-a<\eta<-b<\lambda_0<\xi<0$. Тогда согласно лемме 1 областью возможного движения будет область $\lambda_1\in[-a,\eta]$, $\lambda_2\in[\lambda_0,\xi]$ (рис. 3). Эта область является проекцией двух торов Лиувилля.

Теперь изменим значения $h$ и $f$ так, чтобы корни удовлетворяли той же цепочке неравенств, что и в рассматриваемом примере. Область возможного движения изменится, но качественно она останется той же (не происходит бифуркации областей). Следовательно, решение данной цепочки неравенств относительно $(h, f)$ является областью в образе отображения момента, все точки которой регулярны. Такую открытую область назовем камерой образа отображения момента. Сам образ отображения момента является объединением замыканий камер, соответствующих всевозможным неравенствам на корни многочлена $V$. Бифуркационная диаграмма $\Sigma(F)$ состоит из некоторых границ камер, однако не все границы камер включаются в диаграмму (это будет показано ниже).

Замечание 2. Будем обозначать бифуркационную диаграмму на образе отображения момента сплошными линиями. Штриховыми линиями будем обозначать те границы камер, которые не принадлежат бифуркационной диаграмме.

Образ отображения момента для рассматриваемого примера изображен на рис. 4.

Определение 12. Области возможного движения, соответствующие регулярным точкам в образе отображения момента, будем называть регулярными, в противном случае – критическими.

Определение 13. Если при непрерывном изменении значений интегралов регулярная область перестраивается в регулярную через критическую область возможного движения, будем говорить, что в этом случае происходит бифуркация областей возможного движения (ей соответствует бифуркация торов Лиувилля в $M^4$).

§ 6. Основной результат

Напомним, что топологическую классификацию элементарных плоских биллиардов (с точностью до эквивалентности, описанной выше) без потенциала была дана В. В. Ведюшкиной в [12]. В той же работе для этого класса биллиардов была дана полная классификация их инвариантов Фоменко–Цишанга (т.е. меченых молекул). В нашей работе мы рассматриваем только плоские элементарные биллиарды (серии $A$ и $B$).

Теорема 6. a) Полный список всех инвариантов Фоменко–Цишанга для рассматриваемых элементарных плоских биллиардов с потенциалом приведен в табл. 2.

b) В табл. 1 представлена более подробная информация о классификации из табл. 2. Во втором столбце перечислены 12 плоских элементарных биллиардов. В третьем столбце перечислены инварианты Фоменко–Цишанга, соответствующие случаю притягивающего потенциала ($k>0$). Указанные номера молекул изображены в табл. 2. В четвертом столбце перечислены инварианты Фоменко–Цишанга, соответствующие случаю отталкивающего потенциала ($k<0$).

c) В табл. 3 приведен полный список всех бифуркационных диаграмм рассматриваемой задачи.

d) Всего неэквивалентных классов биллиардов указанного типа 12, а различных инвариантов Фоменко–Цишанга 18.

e) Оказывается, все меченые молекулы, соответствующие притягивающему потенциалу, встречаются среди меченых молекул, соответствующих отталкивающему потенциалу. В этом смысле биллиард с отталкивающим потенциалом “богаче” биллиарда с притягивающим потенциалом.

Таблица 2.Инварианты Фоменко–Цишанга

Инвариант Фоменко–ЦишангаИнвариант Фоменко–Цишанга
110
211
312
413
514
615
716
817
918

Таблица 3.Бифуркационные диаграммы.

№ биллиардаДиаграмма при $k>0$Диаграмма при $k<0$
1
2
3
4
5
6
7
8
9
10
11
12

Замечание 3. Итак, все инварианты для биллиардов притягивающего потенциала $k>0$ встречаются в случае отталкивающего потенциала $k<0$. Этот интересный эффект объясняется тем, что в случае $k>0$ бифуркация областей возможного движения может происходить только по отрезку фокальной прямой, а при $k<0$ бифуркация происходит как по фокальной, так и по вертикальной прямой.

§ 7. Доказательство основного результата (теорема 6)

Доказательство, а именно вычисление меченых молекул, проведем в несколько шагов.

Доказательство теоремы 6. Шаг 1. Докажем, что прообраз любой регулярной области возможного движения в фазовом пространстве $M^4$ гомеоморфен объединению двумерных торов. Для этого заметим, что любая регулярная область является либо объединением нескольких односвязных подобластей, либо областью, гомеоморфной кольцу. Сначала докажем, что прообраз любой односвязной области гомеоморфен тору (см. рис. 2). Разобьем область на эллипсы, софокусные стенкам биллиарда. Эти эллипсы делятся на дуги областью возможного движения. Напомним, что с помощью формулы (4.6) каждой точке области возможного движения можно сопоставить либо четыре, либо два вектора, либо один нулевой вектор. Оснастим каждую точку полученных дуг софокусных эллипсов векторами скорости согласно формуле (4.6). Для внутренних точек области это четыре вектора, для границ – либо два, касательных к границе, либо два, полученных из четырех векторов отождествлением согласно закону отражения. В углах области возможного движения (все они прямые) возникает единственный нулевой вектор. Прообразом в $M^4$ каждой полученной дуги оснащенного софокусного эллипса, не являющегося граничным, являются две окружности. Умножим эти окружности на отрезок и склеим полученные цилиндры по граничным окружностям в соответствии с тем, что на софокусных эллипсах, лежащих на границе области, векторы скоростей отождествляются по закону отражения. Заметим, что если задать ориентацию на граничных окружностях, согласованную с ориентацией цилиндров, то при склейке пар окружностей их ориентации в паре либо одновременно одинаковые, либо одновременно противоположные. Следовательно, мы получили двумерный тор (а не бутылку Клейна).

Теперь разберем случай области, гомеоморфной кольцу (рис. 5). Как и выше, разобьем область возможного движения на оснащенные софокусные эллипсы. Прообраз в $M^4$ каждого такого не граничного эллипса гомеоморфен четырем окружностям. Умножим их на отрезок и склеим цилиндры по границам, согласно склейке векторов скорости на граничных эллипсах. Как и выше, задав ориентацию на границах цилиндров, получим два тора. Эти торы отвечают за движение точки на биллиарде по и против часовой стрелки соответственно. Отметим, что в этом случае движения точки по отрезкам вертикальной или фокальной прямых не возникает. Дело в том, что такое движение было бы критическим и соответствовало бы критическим областям возможного движения, а такие на этом шаге не рассматриваются (о них см. шаг 2).

Шаг 2. Теперь вычислим атомы. Для этого посмотрим, как перестраиваются регулярные области возможного движения при переходе из одной камеры бифуркационной диаграммы в другую. Заметим, что для всех биллиардов перестройки областей возможного движения возникают при переходе из одной камеры образа отображения момента в другую. Приведем классификацию (полный список) таких перестроек:

(0) выход границы области на границу биллиарда;

(1) склейка односвязных подобластей по одному отрезку фокальной прямой;

(2) склейка односвязных подобластей по двум отрезкам фокальной прямой в одну неодносвязную область;

(3) склейка односвязных подобластей по одному отрезку вертикальной прямой;

(4) склейка односвязных подобластей по двум отрезкам вертикальной прямой в одну неодносвязную область;

(5) перестройка односвязной области в односвязную по фокальной прямой через фокус.

Рассмотрим все эти случаи.

В случае выхода области возможного движения на границу (0) эта область не меняется с точностью до эквивалентности, описанной выше (см. определение 5). Следовательно, бифуркаций торов не происходит, а соответствующая точка в образе отображения момента лежит на штриховой линии.

Рассмотрим случай (1) (рис. 6). Критическим движением является отрезок фокальной прямой, по которому происходит склейка. Разобьем критическую область возможного движения на оснащенные эллипсы. Внутренним точкам отрезка склейки соответствуют два вектора скорости, направленные по фокальной прямой, во всех остальных точках неграничного эллипса – по четыре вектора. Следовательно, в прообразе неграничного эллипса лежат две восьмерки (критические слои 2-атома $B$). Умножим их на отрезок. Получим два “цилиндра”. Границей каждого из них является пара восьмерок. Склеим эти два “цилиндра” по их границам в соответствии со склейкой векторов скорости на границе области возможного движения. Так же, как и на шаге 1, зададим ориентацию на границах “цилиндров”, согласованную с ориентацией самих “цилиндров”. При склейке соответствующих границ в двух склеиваемых парах восьмерок их ориентация одновременно либо одинаковая, либо противоположная. В обоих случаях получается один и тот же критический слой 3-атома $B$. Иными словами, бифуркация торов при перестройке типа (1) происходит по атому $B$. Аналогично рассматривается случай (3) (рис. 7).

Теперь рассмотрим случай (2) (рис. 8). Критическим движением являются отрезки фокальной прямой, по которым происходит склейка. Разобьем критическую область возможного движения на оснащенные эллипсы. Во внутренних точках отрезка склейки лежат два вектора скорости, направленные по фокальной прямой, во всех остальных точках неграничного эллипса – по четыре вектора. Следовательно, прообраз каждого неграничного эллипса гомеоморфен двум парам окружностей. В каждой паре окружности имеют две общие точки. Иными словами, имеем два экземпляра особых слоев 2-атома $C_2$. Умножим каждый из них на отрезок и склеим по границе согласно правилу склеек векторов скоростей на границе области. Как и в предыдущих случаях, учет ориентаций дает нам особый слой атома $C_2$. Случай (4) разбирается аналогично (рис. 9).

Случай (5) более сложен (рис. 10). Всем внутренним точкам области возможного движения соответствуют четыре вектора скорости, за исключением точек фокальной прямой. Фокусу соответствует целая окружность векторов скорости, а в других точках фокальной прямой имеем по два вектора скорости, направленные вдоль нее. Такой перестройке соответствует бифуркация через атом $A^*$. Подробнее об этом вычислении см. работу В. В. Ведюшкиной [12].

Стоит отметить, что кроме приведенных случаев бифуркаций имеет место бифуркация односвязной области возможного движения в неодносвязную через два фокуса. Такая бифуркация встречается только в случае эллиптического биллиарда. Этот случай биллиарда в эллипсе с потенциалом был полностью разобран И. Ф. Кобцевым.

Таким образом, мы вычислили атомы для всех возможных перестроек областей возможного движения нашей задачи.

Шаг 3. После вычисления атомов следующим шагом является вычисление грубых молекул. Для этого выберем уровень интеграла $H=h_0$ и на образе отображения момента проведем прямую $h=h_0$. Эта прямая пересекает его по конечному отрезку, а бифуркационную диаграмму – в нескольких точках на этом отрезке. Каждой такой точке соответствует особый слой некоторого атома в прообразе отображения момента, всем остальным точкам соответствует объединение торов (в соответствии с областью возможного движения согласно утверждениям на шаге 1). Сам атом вычисляется согласно утверждениям шага 2. В качестве примера рассмотрен кольцевой биллиард при $k>0$ (рис. 11 и 12).

Итогом этого шага является полный список всех грубых молекул нашей задачи. Все они перечислены в табл. 2. В этой таблице все грубые молекулы снабжены метками, вычислению которых посвящен следующий шаг 4.

Шаг 4. Это последний шаг доказательства, на котором мы вычислим все метки.

Определение 14. Полуокрестностью атома называется прообраз в $Q^3$ полуокрестности соответствующей этому атому критической точки на отрезке, образованном пересечением образа отображения момента и прямой $h=h_0$.

Для вычисления меток выберем допустимые базисы в каждой полуокрестности каждого атома, вычислим матрицы склейки и уже по ним вычислим метки. Более подробно об алгоритме вычисления меток см. [11] и [12].

Начнем с биллиардов, являющихся частью кольцевого биллиарда. Это – биллиарды 2–8 в табл. 1. Они характеризуются тем, что проекциями критических циклов атомов являются отрезки фокальной и вертикальной прямых, причем эти отрезки лежат в области, ограниченной двумя эллипсами. Выберем базис вблизи атомов, степени вершин которых в графе превосходят 1 (т.е. это не атомы $A$). Базисный цикл $\lambda$ проходит вблизи критического цикла, базисный цикл $\mu$ является дополнительным. Для рассматриваемых биллиардов прообраз отрезка фокальной прямой и вертикальной прямой соответствует одним и тем же базисным циклам с точностью до ориентации. Поэтому любая матрица склейки на внутренних ребрах имеет вид

$$ \begin{equation} \begin{pmatrix} 1 & 0\\ 0 & -1 \end{pmatrix}. \end{equation} \tag{7.1} $$
Соответственно, метки будут равны $r=\infty$, $\varepsilon=1$.

Чтобы вычислить матрицы склейки на граничных ребрах молекулы, найдем допустимый базис на торах Лиувилля вблизи атомов $A$. Заметим, что при минимальных значениях интеграла $F$ (т.е. там, где торы “рождаются”), соответствующих какой-либо рассматриваемой молекуле, проекция цикла $\lambda$ на биллиард стягивается по отрезку вертикальной прямой. Следовательно, для ребер, соответствующих подобным атомам $A$, матрица склейки оказывается такой же, как и в предыдущем случае, т.е.

$$ \begin{equation*} \begin{pmatrix} 1 & 0\\ 0 & -1 \end{pmatrix}. \end{equation*} \notag $$
Следовательно, метки имеют вид $r=\infty$, $\varepsilon=1$.

При максимальных значениях интеграла $F$ (т.е. там, где торы “умирают”, исчезают) проекция базисного цикла $\lambda$ на биллиард стягивается по эллипсу, софокусному границам биллиарда. Матрица склейки имеет вид

$$ \begin{equation*} \begin{pmatrix} 0 & 1\\ 1 & 0 \end{pmatrix}. \end{equation*} \notag $$
Следовательно, метки таковы: $r=0$, $\varepsilon=1$. Таким образом, мы вычислили метки на всех ребрах грубых молекул для всех описанных выше биллиардов (т.е. для являющихся частью кольцевой области, см. выше). Иными словами, подсчитаны полные инварианты Фоменко–Цишанга. Пример перестроек областей и базисных циклов для молекулы 16 приведен на рис. 13. Здесь циклы $\lambda$ обозначены сплошными линиями, циклы $\mu$ – штриховыми.

Остальные случаи рассмотрим по отдельности. Напомним, что случай эллипса полностью проанализирован И. Ф. Кобцевым, поэтому здесь мы его не рассматриваем.

Биллиард 9. Из шага 3 следует, что все грубые молекулы для этого случая оказываются двух типов: 1 и 17 (см. табл. 2). Рассмотрим перестройку и циклы, соответствующие молекуле 1 (рис. 14).

Вычислим матрицы склеек. Для указанного на рис. 14 перехода номер 1 матрица склейки равна

$$ \begin{equation} \begin{pmatrix} 0 & 1\\ 1 & 0 \end{pmatrix}. \end{equation} \tag{7.3} $$
Следовательно, метки таковы: $r=0$, $\varepsilon=1$. Для указанного на рис. 14 перехода номер 2 матрица склейки имеет вид
$$ \begin{equation} \begin{pmatrix} 1 & 0\\ 0 & -1 \end{pmatrix}. \end{equation} \tag{7.4} $$
Следовательно, метки таковы: $r=\infty$, $\varepsilon=1$.

Поскольку молекула 17 имеет простой вид $A$–$A$, то этот случай разбирается тривиально.

Биллиард 11. Из шага 3 следует, что все грубые молекулы для этого случая оказываются двух типов: 8 и 17 (см. табл. 2). Рассмотрим перестройку и циклы, соответствующие молекуле 8 (рис. 15).

Заметим, что перестройка областей возможного движения происходит по фокальной прямой через фокус.

Матрица склейки для изображенного на рис. 15 перехода номер 1 имеет вид

$$ \begin{equation} \begin{pmatrix} 2 & 1\\ 1 & 0 \end{pmatrix}. \end{equation} \tag{7.5} $$
Следовательно, метки имеют вид: $r=0$, $\varepsilon=1$. Матрица склейки для изображенного на рис. 15 перехода номер 2 имеет вид
$$ \begin{equation} \begin{pmatrix} -1 & 1\\ -1 & 2 \end{pmatrix}. \end{equation} \tag{7.6} $$
Следовательно, метки имеют вид: $r=0$, $\varepsilon=1$. Метка $n={-0}/{1}+{-1}/{1}=-1$.

Как уже было отмечено, молекула 17 разбирается тривиально.

Биллиард 10. Все грубые молекулы для этого случая оказываются типов 1, 6, 7, 17 (см. табл. 2).

Вычисление меток для молекулы 1 проводится аналогично той же молекуле для биллиарда 9.

Рассмотрим перестройки областей возможного движения и проекции базисных циклов, соответствующие молекуле 6 (рис. 16).

Матрица склейки для изображенного на рис. 16 перехода номер 1 имеет вид

$$ \begin{equation} \begin{pmatrix} 0 & 1\\ 1 & 0 \end{pmatrix}. \end{equation} \tag{7.7} $$
Следовательно, метки имеют следующий вид: $r=0$, $\varepsilon=1$. Для всех остальных перестроек областей возможного движения матрица склейки имеет вид
$$ \begin{equation} \begin{pmatrix} 1 & 0\\ 0 & -1 \end{pmatrix}. \end{equation} \tag{7.8} $$
Следовательно, метки имеют следующий вид: $r=\infty$, $\varepsilon=1$. Таким образом, мы вычислили все метки молекулы 6.

Теперь разберем молекулу 7. На рис. 17 изображены соответствующие этой молекуле перестройки областей возможного движения.

Каждая область соответствует одной из трех групп ребер в молекуле. При этом, как нетрудно заметить, на каждом ребре при переходе из полуокрестности одного атома в полуокрестность другого базисные циклы $\lambda$ и $\mu$ меняются местами с точностью до ориентации. Следовательно, все матрицы склеек имеют вид

$$ \begin{equation} \begin{pmatrix} 0 & 1\\ 1 & 0 \end{pmatrix}. \end{equation} \tag{7.9} $$
Следовательно, метки имеют вид: $r=0$, $\varepsilon=1$. Метка $n=0$ для каждого атома $B$.

Случай молекулы 17 тривиален.

Биллиард 12. Все грубые молекулы для этого случая оказываются типов 5, 8, 9, 10, 17 (см. табл. 2). Аналогично предыдущим случаям, рассмотрим каждую молекулу.

Случай молекулы 8 разбирается аналогично той же молекуле для биллиарда 11. В случае молекулы 5 области возможного движения перестраиваются так, что ни одна из них не пересекается с некоторой окрестностью фокального отрезка (рис. 18). Следовательно, он аналогичен случаю биллиарда, лежащего целиком в кольце, и полностью разобран.

Теперь рассмотрим молекулу 9. Соответствующие перестройки областей возможного движения и базисные циклы изображены на рис. 19.

Заметим, что на ребрах, трансцендентных атому $A^*$, матрицы склейки те же, что и в случае биллиарда 10, так как базисные циклы остаются теми же. Для остальных перестроек матрицы вычисляются аналогично рассмотренным ранее биллиардам, лежащим в кольце.

Наконец, рассмотрим случай молекулы 10. Перестройки областей возможного движения и базисные циклы изображены на рис. 20.

Матрица склейки для изображенного на рис. 20 перехода номер 1 имеет вид

$$ \begin{equation} \begin{pmatrix} 2 & 1\\ 1& 0 \end{pmatrix}. \end{equation} \tag{7.10} $$
Следовательно, метки имеют следующий вид: $r=0$, $\varepsilon=1$.

Матрица склейки для изображенного на рис. 20 перехода номер 3 имеет вид

$$ \begin{equation} \begin{pmatrix} 1 & -2\\ -1& 1 \end{pmatrix}. \end{equation} \tag{7.11} $$
Следовательно, метки имеют вид: $r={1}/{2}$, $\varepsilon=-1$. Метка $n$ атома $A^*$ имеет вид $n=-1+0=-1$.

Матрица склейки для изображенного на рис. 20 перехода номер 5 имеет вид

$$ \begin{equation} \begin{pmatrix} 0 & 1\\ 1& 0 \end{pmatrix}. \end{equation} \tag{7.12} $$
Следовательно, метки имеют вид: $r=0$, $\varepsilon=-1$. Метка $n$ атома $B$ равна $n=0+0=0$.

Итак, были получены инварианты Фоменко–Цишанга для всех изоэнергетических поверхностей $Q^3$ для каждого из двенадцати биллиардов. Следовательно, теорема 6 доказана.

§ 8. Лиувиллево эквивалентные системы

Теперь приведем примеры интегрируемых систем, лиувиллево эквивалентных рассмотренным биллиардам. Согласно теореме Фоменко–Цишанга 2 для этого достаточно сопоставить соответствующие инварианты Фоменко–Цишанга с полученными в теореме 6.

Теорема 7. В табл. 4 приведены интегрируемые системы динамики твердого тела, лиувиллево эквивалентные рассматриваемым биллиардам с гуковским потенциалом в своих соответствующих зонах энергии. При этом номера биллиардов см. в табл. 1.

Таблица 4.Эквивалентные системы

Инвариант Фоменко–Цишанга№ биллиардаСлучай динамики твердого тела
Все перечисленныеЛагранж, Эйлер
$2$ при $k<0,3,4$, $5$ при $k<0$, $6$ при $k<0,8,9,10$Ковалевская, Жуковский, Горячев–Чаплыгин–Сретенский, Ковалевская–Яхья
1Жуковский
1 при $k<0, 11, 12$Горячев–Чаплыгин–Сретенский
1 при $k<0, 2$Эйлер, Клебш
1 при $k<0$Ковалевская

Список литературы

1. А. Т. Фоменко, Х. Цишанг, “О топологии трехмерных многообразий, возникающих в гамильтоновой механике”, Докл. АН СССР, 294:2 (1987), 283–287  mathnet  mathscinet  zmath; англ. пер.: A. T. Fomenko, H. Zieschang, “On the topology of the three-dimensional manifolds arising in Hamiltonian mechanics”, Soviet Math. Dokl., 35:2 (1987), 529–534
2. А. Т. Фоменко, “Топологические инварианты гамильтоновых систем, интегрируемых по Лиувиллю”, Функц. анализ и его прил., 22:4 (1988), 38–51  mathnet  mathscinet  zmath; англ. пер.: A. T. Fomenko, “Topological invariants of Liouville integrable Hamiltonian systems”, Funct. Anal. Appl., 22:4 (1988), 286–296  crossref
3. А. Т. Фоменко, “Теория бордизмов интегрируемых гамильтоновых невырожденных систем с двумя степенями свободы. Новый топологический инвариант многомерных интегрируемых систем”, Изв. АН СССР. Сер. матем., 55:4 (1991), 747–779  mathnet  mathscinet  zmath; англ. пер.: A. T. Fomenko, “A bordism theory for integrable nondegenerate Hamiltonian systems with two degrees of freedom. A new topological invariant of higher-dimensional integrable systems”, Math. USSR-Izv., 39:1 (1992), 731–759  crossref  adsnasa
4. А. Т. Фоменко, “Топологический инвариант, грубо классифицирующий интегрируемые строго невырожденные гамильтонианы на четырехмерных симплектических многообразиях”, Функц. анализ и его прил., 25:4 (1991), 23–35  mathnet  mathscinet  zmath; англ. пер.: A. T. Fomenko, “A topological invariant which roughly classifies integrable strictly nondegenerate Hamiltonians on four-dimensional symplectic manifolds”, Funct. Anal. Appl., 25:4 (1991), 262–272  crossref
5. Е. А. Кудрявцева, И. М. Никонов, А. Т. Фоменко, “Максимально симметричные клеточные разбиения поверхностей и их накрытия”, Матем. сб., 199:9 (2008), 3–96  mathnet  crossref  mathscinet  zmath; англ. пер.: E. A. Kudryavtseva, I. M. Nikonov, A. T. Fomenko, “Maximally symmetric cell decompositions of surfaces and their coverings”, Sb. Math., 199:9 (2008), 1263–1353  crossref  adsnasa
6. A. T. Fomenko, A. Yu. Konyaev, “New approach to symmetries and singularities in integrable Hamiltonian systems”, Topology Appl., 159:7 (2012), 1964–1975  crossref  mathscinet  zmath
7. В. В. Ведюшкина (Фокичева), А. Т. Фоменко, “Интегрируемые топологические биллиарды и эквивалентные динамические системы”, Изв. РАН. Сер. матем., 81:4 (2017), 20–67  mathnet  crossref  mathscinet  zmath; англ. пер.: V. V. Vedyushkina (Fokicheva), A. T. Fomenko, “Integrable topological billiards and equivalent dynamical systems”, Izv. Math., 81:4 (2017), 688–733  crossref  adsnasa
8. В. В. Ведюшкина (Фокичева), А. Т. Фоменко, “Интегрируемые геодезические потоки на ориентируемых двумерных поверхностях и топологические биллиарды”, Изв. РАН. Сер. матем., 83:6 (2019), 63–103  mathnet  crossref  mathscinet  zmath; англ. пер.: V. V. Vedyushkina (Fokicheva), A. T. Fomenko, “Integrable geodesic flows on orientable two-dimensional surfaces and topological billiards”, Izv. Math., 83:6 (2019), 1137–1173  crossref  adsnasa
9. В. И. Драгович, “Интегрируемые возмущения биллиарда Биркгофа внутри эллипса”, ПММ, 62:1 (1998), 166–169  mathscinet  zmath; англ. пер.: V. I. Dragovich, “Integrable perturbations of a Birkhoff billiard inside an ellipse”, J. Appl. Math. Mech., 62:1 (1998), 159–162  crossref
10. С. Е. Пустовойтов, “Топологический анализ биллиарда в эллиптическом кольце в потенциальном поле”, Фундамент. и прикл. матем., 22:6 (2019), 201–225  mathnet  mathscinet
11. А. В. Болсинов, А. Т. Фоменко, Интегрируемые гамильтоновы системы. Геометрия, топология, классификация, т. 1, 2, Изд. дом “Удмуртский университет”, Ижевск, 1999, 444 с., 447 с.  mathscinet  zmath; англ. пер.: A. V. Bolsinov, A. T. Fomenko, Integrable Hamiltonian systems. Geometry, topology, classification, Chapman & Hall/CRC, Boca Raton, FL, 2004, xvi+730 с.  crossref  mathscinet  zmath
12. В. В. Фокичева, “Топологическая классификация биллиардов в локально плоских областях, ограниченных дугами софокусных квадрик”, Матем. сб., 206:10 (2015), 127–176  mathnet  crossref  mathscinet  zmath; англ. пер.: V. V. Fokicheva, “A topological classification of billiards in locally planar domains bounded by arcs of confocal quadrics”, Sb. Math., 206:10 (2015), 1463–1507  crossref  adsnasa
13. В. В. Козлов, “Некоторые интегрируемые обобщения задачи Якоби о геодезических на эллипсоиде”, ПММ, 59:1 (1995), 3–9  mathscinet  zmath; англ. пер.: V. V. Kozlov, “Some integrable extensions of Jacobi's problem of geodesics on an ellipsoid”, J. Appl. Math. Mech., 59:1 (1995), 1–7  crossref
14. М. П. Харламов, “Топологический анализ и булевы функции: I. Методы и приложения к классическим системам”, Нелинейная динам., 6:4 (2010), 769–805  mathnet

Образец цитирования: С. Е. Пустовойтов, “Топологический анализ биллиарда, ограниченного софокусными квадриками, в потенциальном поле”, Матем. сб., 212:2 (2021), 81–105; S. E. Pustovoitov, “Topological analysis of a billiard bounded by confocal quadrics in a potential field”, Sb. Math., 212:2 (2021), 211–233
Цитирование в формате AMSBIB
\RBibitem{Pus21}
\by С.~Е.~Пустовойтов
\paper Топологический анализ биллиарда, ограниченного софокусными квадриками, в потенциальном поле
\jour Матем. сб.
\yr 2021
\vol 212
\issue 2
\pages 81--105
\mathnet{http://mi.mathnet.ru/sm9418}
\crossref{https://doi.org/10.4213/sm9418}
\mathscinet{http://mathscinet.ams.org/mathscinet-getitem?mr=4223963}
\zmath{https://zbmath.org/?q=an:1468.37028}
\adsnasa{https://adsabs.harvard.edu/cgi-bin/bib_query?2021SbMat.212..211P}
\elib{https://elibrary.ru/item.asp?id=46036139}
\transl
\by S.~E.~Pustovoitov
\paper Topological analysis of a~billiard bounded by confocal quadrics in a~potential field
\jour Sb. Math.
\yr 2021
\vol 212
\issue 2
\pages 211--233
\crossref{https://doi.org/10.1070/SM9418}
\isi{https://gateway.webofknowledge.com/gateway/Gateway.cgi?GWVersion=2&SrcApp=Publons&SrcAuth=Publons_CEL&DestLinkType=FullRecord&DestApp=WOS_CPL&KeyUT=000701455800001}
\scopus{https://www.scopus.com/record/display.url?origin=inward&eid=2-s2.0-85105114739}
Образцы ссылок на эту страницу:
  • https://www.mathnet.ru/rus/sm9418
  • https://doi.org/10.4213/sm9418
  • https://www.mathnet.ru/rus/sm/v212/i2/p81
  • Эта публикация цитируется в следующих 10 статьяx:
    Citing articles in Google Scholar: Russian citations, English citations
    Related articles in Google Scholar: Russian articles, English articles
    Математический сборник Sbornik: Mathematics
    Статистика просмотров:
    Страница аннотации:460
    PDF русской версии:58
    PDF английской версии:23
    HTML русской версии:249
    Список литературы:24
    Первая страница:22
     
      Обратная связь:
     Пользовательское соглашение  Регистрация посетителей портала  Логотипы © Математический институт им. В. А. Стеклова РАН, 2024