Описание особенностей системы "бильярд в эллипсе" Текст научной статьи по специальности «Математика»

УДК 517.938.5

ОПИСАНИЕ ОСОБЕННОСТЕЙ СИСТЕМЫ "БИЛЬЯРД В ЭЛЛИПСЕ"

В. В. Фокичева1

Изучается интегрируемая система "бильярд в эллипсе", возникающая при описании движения точки внутри эллипса с естественным отражением на границе. В работе вычислен топологический инвариант лиувиллевой эквивалентности этой системы — молекула Фоменко-Цишанга — с помощью нового метода, разработанного автором.

Ключевые слова: интегрируемая система, бильярд, лиувиллева эквивалентность, молекула Фоменко-Цишанга.

A "billiard within an ellipse" is an integrable system appearing in the description of a point motion inside an ellipse with natural reflections at the boundary. This system is considered in the paper, the topological invariant of Liouville equivalence of this system is calculated, which is a Fomenko-Tsishang molecule, by the new method developed by the author.

Key words: integrable system, billiard, Liouville equivalence, marked molecule.

1. Введение. Данная система и подобные ей изучались многими авторами (см., например, [1-3]). Впервые инвариант Фоменко-Цишанга для системы "бильярд в эллипсе" был вычислен в статье [1] без доказательства. Вскоре в книге [2] этот инвариант был уточнен, хотя доказательство также не было указано. В данной работе впервые приводится полное доказательство. Предлагаемый новый метод позволяет эффективно вычислять "грубую молекулу". Способ нахождения меток полностью аналогичен использованному А. В. Болсиновым и А. Т. Фоменко в случае геодезического потока на эллипсоиде (задача Яко-би) [4].

2. Постановка задачи. Опишем достаточно широкий класс "бильярдных систем", к которым применим наш метод. Пусть дана область на плоскости R2, ограниченная несколькими квадриками из софо-кусного семейства:

x2 y2

+ т^—г = 1, а > 0,6 > 0,

а - Л Ь - Л

причем область не содержит углов Система "бильярд в области" описывает движение точки внутри области О с естественным отражением на границе Р = дО. В точках, где граница Р не гладкая (где угол излома равен систему доопределим по непрерывности: попав в угол, точка, не теряя скорости, отразится по той же траектории. Таким образом, фазовым пространством системы является многообразие

М4 := {(х,у)\ х е О, V е ТЖМ2 , IV > 0}/

где (Х1 ) ~ (Х2^2) ^ Х1 = Х2 е Р, \ VI \ = ^2 \ и VI - V2 ± ТХ1 Р.

Теорема (Якоби, Шаль) [5]. Касательные прямые к геодезической линии на квадрике в п-мерном пространстве, проведенные во всех точках геодезической, касаются кроме этой квадрики еще п — 2 конфокальных с ней квадрик, одних и тех же для всех точек геодезической.

В нашем случае это означает, что касательные в любой точке траектории касаются эллипса или гиперболы, софокусных с семейством границы Р. Таким образом, система обладает двумя интегралами:

1) V — модуль вектора скорости,

2) Л — параметр софокусной квадрики.

Будем считать, что скорость фиксирована и равна единице. Ограничим фазовое пространство М4 на поверхность уровня \ V \ =1 и получим изоэнергетическую поверхность Q3.

Определение. Пусть (М4,^,/!,^) и (М|,^2,/2,$2) — две интегрируемые системы. Рассмотрим изоэнергетические поверхности Q! = {х е М4 : Д(х) = С1} и Q3 = {х е М| : /2(х) = С2}. Они называются лиувиллево эквивалентными, если существует послойный диффеоморфизм Q3 ^ Q2, который, кроме того, сохраняет ориентацию 3-многообразий Q3 и Q3 и ориентацию всех критических окружностей (подробнее см. [6-9]).

1 Фокичева Виктория Викторовна — студ. каф. дифференциальной геометрии и приложений мех.-мат. ф-та МГУ, e-mail:

arinir@yandex.ru.

В силу теоремы Лиувилля (см., например, [6]) многообразие Q3 расслоено на торы и особые слои (фактически оно представляет собой склейку окрестностей особых слоев друг с другом). Профакторизу-ем Q3 по торам Лиувилля и получим граф, называемый грубой молекулой. В его вершинах расположены "атомы", описывающие соответствующие бифуркации торов Лиувилля (в настоящей работе используются атомы А, В, С2). Однако этот граф не описывает полностью слоение Лиувилля, так как теряется информация о склейках особых слоев. Чтобы это исправить, необходимо выбрать пары так называемых допустимых базисов на граничных торах (по правилам, указанным в [6]) и указать матрицы перехода от одного базиса к другому. Из полученных матриц перехода — матриц склейки — вычисляются метки г, е и п, которые, будучи расставленными на грубой молекуле, полностью определяют слоение Лиувилля с точностью до послойной эквивалентности и уже не зависят от выбора базиса на граничных торах [6].

3. Построение меченой молекулы. Выделим из описанного выше класса систем "бильярд в области" простейшую — "бильярд в эллипсе".

Теорема. Слоение Лиувилля изоэнергети-ческого многообразия Q3 системы "бильярд в эллипсе" при ненулевом значении скорости V описывается меченой молекулой Фоменко-Ци-шанга Ш*, изображенной на рис. 1, а. При этом Q3 диффеоморфно стандартной сфере Б3.

Замечание. Таким образом, оказывается, что система "бильярд в эллипсе" лиувиллево не эквивалентна ни одной из ранее изученных систем с двумя степенями свободы. Большой список этих систем приведен в [6].

Доказательство. Рассмотрим геодезический поток на эллипсоиде и обозначим эту систему У\. Ее особенности ранга 1 невырождены, и сама система описывается меченой молекулой, найденной в работе [4] и изображенной на рис. 1, б.

Интегрируемость системы V следует из теоремы Якоби, Шаля. Систему бильярда в эллипсе обозначим У2. Система — система, описывающая геодезический поток на верхней половине эллипсоида (г ^ 0). Введя эллиптические координаты, легко увидеть, что данная система эквивалентна системе У2 — "бильярду в эллипсе".

Построение грубой молекулы Ш. Определим отображение / системы У\ на систему Р3 следующим естественным образом: точка, оснащенная вектором (x,v) € Ш1, переходит в себя, если х лежит в верхней половине эллипсоида, и в получающиеся путем зеркальной симметрии относительно горизонтальной плоскости точку и вектор, если х лежит в нижней половине эллипсоида.

Напомним (подробности см. в работе [4]), что проекция тора ребра молекулы (см. рис. 1, б) на эллипсоид представляет собой кольцо. Рассмотрим уровень интеграла Л = Л, соответствующий двум торам Т1 и Т2 на нижних ребрах молекулы (этому уровню соответствует однополостный гиперболоид). Геодезические, принадлежащие этим торам, лежат в кольце К, которое однополостный гиперболоид с параметром Л отсекает от эллипсоида.

Рассмотрим тор Т1 , на котором лежат траектории, закручивающиеся по часовой стрелке на этом кольце К. Тогда Т1 можно описать как результат склейки на границе двух копий К' и К'' кольца К. Одна копия кольца (К') описывает точки, оснащенные векторами скорости, направленными "вверх", т.е. к верхнему эллипсу, по которому однополостный гиперболоид пересекает исходный эллипсоид. Вторая копия (К'') описывает точки, оснащенные векторами "вниз". На граничных эллипсах ("верхнем" и "нижнем") векторы скорости геодезических направлены горизонтально (так как в этих точках они должны касаться гиперболоида). Таким образом, потоки, определенные на этих кольцах, склеиваются при этом гладким образом и порождают гладкий поток на торе Т1. Аналогичную конструкцию проделаем с остальными ребрами.

Минимальному уровню интеграла Л на Q3 соответствует движение точки на экваторе эллипсоида, являющемся его сечением плоскостью Оху. Максимальному уровню интеграла Л на Q3 соответствует движение точки на экваторе эллипсоида, являющемся его сечением плоскостью Оуг. Данные особенности описываются атомами А. Седловому значению интеграла Л соответствует движение точки в плоскости Охг. Эта особенность описывается атомом С2

Теперь рассмотрим, как преобразуются торы и особые слои при отображении /. Возьмем тор Т1 и его представление в виде склейки двух колец. При отображении / эти два кольца переходят в две половинки кольца (рис. 2, а), которые по закону отражения на границе бильярда (мы уже находимся в системе Р3)

А А

Рис. 1. Молекула Фоменко-Цишанга: а — для системы "бильярд в элипсе"; б — для системы геодезического потока на эллипсоиде

склеиваются между собой. Таким образом, два тора Т и Т2, лежащие на одном уровне интеграла Л, также переходят в два тора.

Рис. 2. Перестройка торов: а — на нижних ребрах; б — на верхних

Рассмотрим два верхних тора и аналогично представим их в виде объединения пары колец. Плоскость Оху, в которой происходит отражение, делит каждое из колец на два прямоугольника, из которых и склеены цилиндры, изображенные на рис. 2, б (из цилиндров склеены торы системы У\). Далее (рис. 2, б) два тора Лиувилля, лежащие на верхних ребрах молекулы, переходят в один тор системы У3.

Таким образом, мы корректно определили отображение торов Лиувилля системы У\ в торы системы У3. Отображение особых слоев слоения Лиувилля может быть получено предельным переходом. При этом атом С2 перейдет в атом В (рис. 3).

Вычисление меток г. Нам понадобится следующая известная

Лемма 1. Рассмотрим произвольное ребро молекулы Ш, и пусть (А+ ), (А-— допустимые системы координат, отвечающие двум атомам, соединенным этим ребром. Считаем, что все эти циклы лежат на одном и том же граничном торе в середине ребра. Тогда если циклы А+ и А- пересекаются в одной точке, то метка г = 0.

Данный факт сразу следует из определения меток г, подробное доказательство приведено в [6].

Ориентируем ребра молекулы из атома В в атомы А. Циклы А+ на торах, близких к атомам А, должны быть выбраны как исчезающие, а циклы А- на торах, близких к В, — как гомологичные критической окружности атома В. Для данной системы они изображены на рис. 4. Легко видеть, что в обоих случаях А+ и А- пересекаются в одной точке, и, таким образом, все метки г = 0.

Склейка

Рис. 3. Перестройка атомов: а — нижних атомов; б — верхних атомов; в — седлового атома

Вычисление меток е. Доказательство почти полностью повторяет доказательство вычисления меток е в статье [4]. Пусть на ребре е = -1. Рассмотрим два базиса на концах ребра — (А+,^+) и (А-,^-). Цикл А- — это "ось" атома В. Цикл — это произвольный цикл на торе, дополнительный к А-. Цикл А+ — это исчезающий цикл тора при стремлении тора к атому А (т.е. цикл, стягивающийся в точку). Цикл — это произвольный цикл на торе, дополнительный к А+. Можно считать, что переходит в критическую минимаксную окружность атома А при стремлении тора к атому. По определению метка е на ребре равна единице, если ориентация базиса (А+ ) совпадает с ориентацией базиса (А+, А-). Если

ориентации противоположны, то е = -1, и в таком случае цикл Л+ лежит между ориентированными циклами Л- и . Пусть это не так. Но циклы Л- и задают предельные положения векторного поля гамильтоновой системы при стремлении тора к атомам А и В. Ввиду того что цикл Л+ расположен между

Л- и на некотором торе ребра молекулы он будет гомологичен замкнутой траектории гамильтонова векторного поля. Однако в силу выбора цикла Л+ данная траектория должна стягиваться в точку на полнотории, соответствующем атому А. Очевидно, что это невозможно. Поэтому метка е = 1.

Вычисление метки п. Так как метки г конечны, то для данной молекулы должна быть вычислена метка п.

Лемма 2. Изоэнергетическая поверхность в случае бильярда в односвязной области гомеоморфна Б3.

Доказательство этого факта напрямую следует из теоремы Ван-Кампена (см. также [10]).

Вычисленные метки г и е позволяют выписать возможные матрицы склейки, если положить метку п неизвестной:

П01 — на верхнем ребре Ш, 1 — на нижних ребрах Ш.

Рис. 4. Циклы Л+ и Л-: а — на нижних ребрах; б — на верхнем ребре; пунктир означает соответствующую интервалу квадрику, а непрерывная линия — цикл

Таким образом, вычисляя фундаментальную группу Q3 по матрицам склейки, получаем, что = Z/nZ, откуда

п = ±1.

Замечание. Так как ориентация Q3 не фиксирована, то при данных значениях меток г и е не фиксирован знак метки п (подробнее см. [6]). Напомним, что аналогичный результат имел место в случае геодезического потока на эллипсоиде [4].

Работа выполнена при поддержке гранта Правительства РФ по постановлению №220, договор M1.G34.31.0054.

СПИСОК ЛИТЕРАТУРЫ

1. Dragovic V., Radnovic M. Bifurcations of Liouville tori in elliptical billiards // arXiv: 0902.4233v2.

2. Драгович В., Раднович М. Интегрируемые биллиарды, квадрики и многомерные поризмы Понселе. М.; Ижевск: НИЦ "Регулярная и хаотическая динамика", 2010.

3. Козлов В.В., Трещев Д.В. Биллиарды. Генетическое введение в динамику систем с ударами. М.: Изд-во МГУ, 1991.

4. Болсинов А.В., Фоменко А.Т. Траекторная классификация геодезических потоков двумерных эллипсоидов. Задача Якоби траекторно эквивалентна интегрируемому случаю Эйлера в динамике твердого тела // Функц. анализ и его прил. 1995. 29, вып. 3. 1-15.

5. Арнольд В.И. Математические методы классической механики. М.: Наука, 1989.

6. Болсинов А.В., Фоменко А.Т. Интегрируемые гамильтоновы системы. Геометрия, топология, классификация. T. 1, 2. Ижевск: РХД, 1999.

7. Фоменко А.Т., Цишанг Х. О типичных топологических свойствах интегрируемых гамильтоновых систем // Изв. АН СССР. 1988. 52, № 2. 378-407.

8. Фоменко А.Т. Симплектическая топология вполне интегрируемых гамильтоновых систем // Успехи матем. наук. 1989. 44, вып. 1 (265). 145-173.

9. Фоменко А. Т., Цишанг Х. Топологический инвариант и критерий эквивалентности интегрируемых гамильтоновых систем с двумя степенями свободы // Изв. АН СССР. Сер. матем. 1990. 54, № 3. 546-575.

10. Gutkin E. Billiard dynamics: a survey with the emphasis on open problems // Regular and Chaotic Dynamics. 2003. 8, N 1. 1-13.

Поступила в редакцию 16.12.2011