Динамические системы
Гамильтоновы системы
NASA запускает зонды к Юпитеру, используя гравитационные манёвры - «горки» на поверхностях постоянной энергии в фазовом пространстве задачи трёх тел. Без понимания симплектической геометрии и интегрируемости такие манёвры были бы невозможны.
- **Молекулярная динамика:** симплектические интеграторы (Верле, Йошида) используются для симуляции белков и материалов - они сохраняют энергию на миллиарды шагов
- **Небесная механика:** стабильность Солнечной системы (на миллиарды лет) связана с сохранением адиабатических инвариантов (переменных действия)
- **Квантовая механика:** принцип соответствия Бора: H(q,p) → Ĥ, {f,g} → (1/iħ)[f̂,ĝ] - квантизация прямо следует из симплектической структуры
Предварительные знания
Hamiltonian
**Маятник качается. Мяч летит по параболе. Планеты обращаются вокруг Солнца.** За всеми этими движениями стоит единая структура - Гамильтонова механика. Вместо уравнений Ньютона F = ma используется функция полной энергии H(q, p), из которой автоматически получаются уравнения движения.
**Гамильтонова система** задаётся функцией H(q, p) - Гамильтонианом (полной энергией), где q - обобщённые координаты, p - обобщённые импульсы. Уравнения движения: **dq_i/dt = ∂H/∂p_i, dp_i/dt = −∂H/∂q_i**. Ключевое свойство: H сохраняется вдоль траекторий: **dH/dt = 0**.
| Система | H(q, p) | Степени свободы |
|---|---|---|
| Гармонический осциллятор | p²/2m + mω²q²/2 | 1 |
| Маятник | p²/2ml² − mgl·cos(q) | 1 |
| Задача Кеплера (орбита) | p²/2m − k/|q| | 3 |
| N частиц | Σpᵢ²/2mᵢ + V(q₁,...,qₙ) | N |
Гамильтон и реформа механики
Уильям Роуэн Гамильтон предложил свой формализм в 1833-1835 годах, первоначально для оптики, а затем распространил на механику. Революция была не в вычислениях - уравнения Ньютона дают те же ответы - а в структуре: симметрии → законы сохранения, фазовое пространство, принцип наименьшего действия. Этот формализм стал основой квантовой механики 100 лет спустя.
Для маятника H = p²/2 − cos(q). Уравнения движения dq/dt и dp/dt:
Symplectic
**Гамильтонов поток - это не просто набор траекторий, это особое геометрическое преобразование.** Оно сохраняет «симплектическую площадь» - обобщение обычной площади на фазовое пространство. Именно это свойство отличает консервативную механику от диссипативной.
**Симплектическое пространство** - пространство (ℝ²ⁿ, ω), где ω = Σ dqᵢ ∧ dpᵢ - симплектическая форма. Поток Гамильтоновой системы - симплектоморфизм: он сохраняет ω. Матрица Якобиана J удовлетворяет: **J^T·Ω·J = Ω**, где Ω - стандартная симплектическая матрица. Это обобщение условия ортогональности для ортогональных преобразований.
| Свойство | Евклидово пространство | Симплектическое пространство |
|---|---|---|
| Метрика | g(u,v) = uᵀv (скалярное произведение) | ω(u,v) = uᵀΩv (кососимметрична) |
| Преобразования | Ортогональные: OᵀO = I | Симплектические: JᵀΩJ = Ω |
| Инвариант | Длины и углы | Симплектическая площадь |
| Размерность | Любая | Всегда чётная: 2n |
**Симплектические интеграторы** - численные методы, специально разработанные для Гамильтоновых систем. Обычные методы (Рунге-Кутта) не сохраняют симплектичность - после миллиона шагов орбита «спирализуется» к центру или от него. Метод Верле и симплектические Рунге-Кутта сохраняют структуру и применяются в молекулярной динамике, небесной механике.
Почему Гамильтоновы системы не могут иметь аттракторов (устойчивых притягивающих множеств)?
Canonical
**Замечательное свойство Гамильтоновых систем - свобода выбора координат.** Можно перейти к любым «каноническим переменным», и уравнения движения сохранят свою форму. Это позволяет выбирать координаты, в которых система выглядит максимально просто - вплоть до полной интегрируемости.
**Каноническое преобразование** - замена (q, p) → (Q, P), сохраняющая форму уравнений Гамильтона. Эквивалентно: преобразование симплектично. Генерируется производящей функцией S: **P = ∂S/∂Q, p = ∂S/∂q**. Пример: переход к переменным действие-угол **(J, θ)**, в которых Гамильтониан зависит только от J: H = H(J).
| Преобразование | Производящая функция | Применение |
|---|---|---|
| Тождественное | S = Σ qᵢPᵢ | Базовый случай |
| Действие-угол | S = ∫p dq (вдоль орбиты) | Интегрируемые системы |
| Пуанкаре-Биркгофф | Ряды по εⁿ | Теория возмущений |
| Преобразование Фурье | exp(iqP/ħ) | Квантовая механика |
В переменных действие-угол (J, θ) для интегрируемой системы уравнения движения принимают вид:
Integrability
**Почему задача двух тел (Земля-Солнце) решается точно, а задача трёх тел - нет?** Потому что двутелая задача интегрируема: есть достаточно законов сохранения, чтобы полностью описать движение. Добавление третьего тела «портит» интегрируемость - и открывает дверь к хаосу.
**Интегрируемая система** (по Лиувиллю-Арнольду) - Гамильтонова система с n степенями свободы, имеющая n инволютивных (попарно коммутирующих) первых интегралов F₁=H, F₂, ..., Fₙ: **{Fᵢ, Fⱼ} = 0**. Теорема Лиувилля-Арнольда: траектории лежат на n-мерных торах Tⁿ, движение квазипериодично.
| Система | Интегрируемая? | Первые интегралы |
|---|---|---|
| Гармонический осциллятор | Да | E = p²/2 + ω²q²/2 |
| Задача Кеплера (2 тела) | Да | E, L (момент), вектор Рунге-Ленца |
| Задача 3 тел | Нет | Только E, P, L - недостаточно |
| Маятник (малые колебания) | Да (приближённо) | E ≈ const (линеаризация) |
Пуанкаре и задача трёх тел
В 1887 году король Швеции Оскар II объявил приз за доказательство устойчивости Солнечной системы. Пуанкаре представил работу, в которой (как ему казалось) доказал это. Но в процессе публикации он обнаружил ошибку. Исправление ошибки привело к открытию хаоса в задаче трёх тел - и к заработку Пуанкаре больше, чем размер приза, на печатание исправленной версии.
Гамильтоновы системы - это только классическая механика
Гамильтонов формализм универсален: квантовая механика (оператор Гамильтона), геометрическая оптика, теория поля, молекулярная динамика - все используют симплектическую структуру.
Симплектическая геометрия лежит в основе геометрической квантовой механики. Операторы квантовой механики - квантование классических наблюдаемых: {f,g} → (1/iħ)[F̂,Ĝ]. Диракт показал, что коммутатор операторов - квантовый аналог скобки Пуассона.
Система с n степенями свободы интегрируема (по Лиувиллю), если она имеет:
Ключевые идеи
- **Гамильтониан H(q,p)** задаёт систему через полную энергию; уравнения dq/dt=∂H/∂p, dp/dt=−∂H/∂q автоматически консервативны
- **Симплектичность:** поток сохраняет ω = Σdqᵢ∧dpᵢ; det(J) = 1; фазовый объём не меняется (теорема Лиувилля)
- **Канонические преобразования** сохраняют форму уравнений; переменные действие-угол делают движение явным (dJ/dt=0, dθ/dt=ω)
- **Интегрируемость по Лиувиллю:** n инволютивных интегралов → движение на торах → квазипериодика; нарушение → KAM-теория и хаос
Связанные темы
Гамильтонова механика - основа для понимания как регулярности, так и хаоса:
- Эргодическая теория — Теорема Лиувилля - канонический пример инвариантности меры в Гамильтоновых системах
- Теория KAM — KAM описывает, что происходит с торами при малых возмущениях интегрируемой системы
- Нейродинамика — Синхронизация осцилляторов в нейронных сетях использует гамильтоновские идеи
Вопросы для размышления
- Теорема Лиувилля утверждает, что фазовый объём сохраняется. Как это связано со вторым законом термодинамики? (Подсказка: энтропия растёт)
- Почему задача трёх тел необычайно сложна, хотя уравнения известны с XVII века? Что именно «ломает» интегрируемость?
- Можно ли использовать симплектические интеграторы для нейронных ОДУ? Какой выигрыш это даст по сравнению с обычным методом Рунге-Кутта?