4. Возрождение динамики
.4. Возрождение динамики
В первой части нашей книги мы рассказали о динамике XIX в. Именно такую динамику излагают многие учебники. Прототипом динамической системы в XIX в. было принято считать интегрируемую систему. Решить уравнения движения означало «удачно» выбрать координаты — так, чтобы соответствующие импульсы были инвариантами движения. Такой подход исключал взаимодействие между частями системы. Ставка на интегрируемые системы провалилась. Как уже упоминалось, в конце XIX в. Брунс и Пуанкаре доказали, что большинство динамических систем, начиная со знаменитой проблемы трех тел, неинтегрируемы.
С другой стороны, сама идея приближения к равновесию, сформулированная на языке теории ансамблей, требовала выхода за пределы идеализации интегрируемых систем. В гл. 8 мы видели, что в теории ансамблей изолированная система находится в равновесии, когда она представлена «микроканоническим ансамблем» — все точки на поверхности заданной энергии
Рис. 32. Временнáя эволюция ячейки в фазовом пространстве р, q. «Объем» ячейки и ее форма сохраняются во времени. Большая часть фазового пространства недоступна для системы.
равновероятны. Это означает, что для системы, стремящейся к равновесию, энергия должна быть единственной величиной, сохраняющейся в ходе эволюции системы. Энергия должна быть единственным инвариантом. При любых начальных условиях система, эволюционируя, должна «побывать» во всех точках поверхности заданной энергии. Для интегрируемых систем энергия — далеко не единственный инвариант. Число инвариантов совпадает с числом степеней свободы, поскольку у интегрируемой системы каждый обобщенный импульс остается постоянным. Следовательно, интегрируемая система «заключена» на весьма ограниченном участке поверхности постоянной энергии (рис. 32) — пересечении всех инвариантных поверхностей.
Чтобы избежать этих трудностей, Максвелл и Больцман ввели новый, совершенно иной тип динамической системы. Для таких систем энергия является единственным инвариантом, а сами системы получили название эргодических систем (рис. 33).
Выдающийся вклад в развитие теории эргодических систем внесли Дж. Биркгоф, фон Нейман, Хопф, Колмогоров и Синай (разумеется, наш перечень далеко не полон). Ныне мы знаем, что существуют обширные классы динамических (но не гамильтоновых) си-
Рис. 33. Типичная эволюция в фазовом пространстве ячейки, соответствующей эргодической системе. «Объем» и форма ячейки сохраняются во времени, но на этот раз ячейка перемещается по всему фазовому пространству.
стем, которые эргодичны. Известно также, что даже сравнительно простые системы могут обладать более сильными свойствами, чем эргодичность. Для таких систем движение в фазовом пространстве становится сильно хаотическим (хотя в полном соответствии с уравнением Луивилля — см. гл. 7 — объем в фазовом пространстве сохраняется).
Предположим, что наше знание начальных условий позволяет нам локализовать систему в малой ячейке фазового пространства. Наблюдая за эволюцией ячейки, мы увидим, как она начнет деформироваться и изгибаться, испуская, подобно амебе, «псевдоножки» по всем направлениям и распространяясь в виде волокон, которые постепенно становятся все тоньше, пока наконец не заполнят все пространство. Ни один самый искусный рисунок не может по достоинству передать
Рис. 34. Типичная эволюция в фазовом пространстве ячейки, соответствующей системе с перемешиванием. Объем по-прежнему сохраняется, но форма уже не остается неизменной: ячейка постепенно размазывается по всему фазовому пространству.
всей сложности реальной ситуации. Действительно, в ходе эволюции системы с перемешиванием две точки, сколь угодно близкие в начальный момент времени, могут разойтись в разные стороны. Даже если бы мы располагали столь обширной информацией о системе, что начальная ячейка, образованная представляющими ее точками, была бы очень мала, динамическая эволюция превратила бы эту миниатюрную область в настоящее геометрическое «чудовище», пронизывающее фазовое пространство своими нитями-щупальцами.
Продемонстрируем различие между устойчивыми и неустойчивыми системами на нескольких простых примерах. Рассмотрим двухмерное фазовое пространство. Через одинаковые промежутки времени станем производить преобразования координат, при которых старая абсцисса р переходит в новую абсциссу р—q, а старая ордината q — в новую ординату р. На рис. 35 показано, что произойдет, если применить эти преобразования
Рис. 35. Преобразование объема в фазовом пространстве, порождаемое дискретным преобразованием: абсцисса р переходит в р—q, ордината q переходит в р. Преобразование циклическое: после шестикратного повторения преобразования исходная ячейка переходит в себя.
к квадрату: квадрат деформируется, но после шестикратного действия преобразования мы возвращаемся к исходному квадрату. Система устойчива: соседние точки преобразуются в соседние. Кроме того, рассмотренное нами преобразование циклическое (после шести операций восстанавливается исходный квадрат).
Рассмотрим теперь два примера сильно неустойчивых систем. Первый пример чисто математический, второй имеет непосредственное отношение к физике. Первая система — преобразование, названное математиками по понятным соображениям преобразованием пекаря Берется квадрат и сплющивается в прямоугольник. Половина прямоугольника отрезается, накладывается на другую половину, а получившийся квадрат снова «раскатывается» в прямоугольник. Последова-
Рис. 36. Реализация «преобразования пекаря» В и обратного преобразования В. Траектории черной и белой точек позволяют понять, как происходит каждое преобразование.
тельность операций, представленная на рис. 36, может быть повторена сколько угодно раз.
Каждый раз квадрат разбивается на части, которые перекладываются в другом порядке. Квадрат в этом примере соответствует фазовому пространству. «Преобразование пекаря» переводит каждую точку квадрата в однозначно определенную новую точку. Хотя последовательность точек-образов вполне детерминистична, «преобразование пекаря» обнаруживает также статистические свойства. Пусть начальное условие для системы состоит в том, что область А квадрата первоначально равномерно заполнена представляющими точками. Можно показать, что, после того как преобразование будет повторено достаточное число раз, начальная ячейка А, каковы бы ни были ее размеры и расположение в квадрате, распадется на отдельные несвязные части (рис. 37). Следовательно, любая область квадрата, независимо от ее размеров, всегда содержит различные траектории, которые при каждом «дроблении» области расходятся. Таким образом, несмотря
Рис. 37. Временнáя эволюция неустойчивой системы. Область А со временем делится на две области A' и А", каждая из которых в свою очередь делится на две подобласти.
на то что эволюция каждой точки в отдельности обратима и детерминистична, описание эволюции любой, даже сколь угодно малой области носит, по существу, статистический характер.
Другим примером простой системы с неожиданно сложным поведением может служить рассеяние твердых шаров. Рассмотрим маленький шарик, отражающийся от больших случайно распределенных шаров. Предположим, что большие шары неподвижны. Такую модель физики называют моделью, или газом, Лоренца в честь выдающегося голландского физика Гендрика Антона Лоренца.
Траектория малого подвижного шарика вполне определена. Но стоит лишь нам ввести в начальные условия небольшую неопределенность, как в результате последовательных столкновений эта неопределенность усилится. Со временем вероятность найти малый шарик равномерно распределится по всему объему, занятому газом Лоренца. Каково бы ни было число преоб-
Рис. 38. Схематическое изображение неустойчивости траектории маленького шарика, отражающегося от больших шаров. Малейшая неточность в задании положения маленького шарика делает невозможным предсказание большого шара, с которым столкнется маленький шарик после первого отражения.
разований, газ никогда не вернется в исходное состояние.
В двух последних примерах динамические системы были сильно неустойчивы. Ситуация, с которой мы сталкиваемся здесь, напоминает неустойчивости в термодинамических системах (см. гл. 5). Произвольно малые различия в начальных условиях усиливаются. В результате переход от ансамблей в фазовом пространстве к индивидуальным траекториям становится невозможным. Описание на языке теории ансамблей мы вынуждены принять за исходный пункт. Статистические понятия перестают быть лишь приближениями к некоторой «объективной истине». Перед такими неустойчивыми системами демон Лапласа оказался бы столь же бессильным, как и мы.
Высказывание Эйнштейна «бог не играет в кости» хорошо известно. Ему созвучно высказывание Пуанкаре о бесконечно мощном духе, беспредельно осведомленном в законах природы, для которого вероятности просто не могли бы существовать. Однако Пуанкаре сам же указал путь к решению проблемы. Он заметил, что когда мы бросаем игральные кости и прибегаем к теории вероятностей, то это отнюдь не означает, будто динамика неверна. Применение вероятностных соображений означает нечто другое. Мы используем понятие вероятности потому, что в любом диапазоне начальных условий, сколь бы малым он ни был, существует «много» траекторий, приводящих к выпадению каждой из граней кости. Именно это и происходит с неустойчивыми динамическими системами. Господь бог, если бы пожелал, мог бы вычислить траектории в нестабильном динамическом мире. При этом он получил бы тот же результат, который нам позволяет получить теория вероятностей. Разумеется, всеведущему богу с его абсолютным знанием было бы нетрудно избавиться от всякой случайности.
Итак, мы можем констатировать, что тесная взаимосвязь между неустойчивостью и вероятностью, несомненно, существует. Это весьма важное обстоятельство, и к его обсуждению мы сейчас перейдем.