Результаты проекта Millenium Simulation. Моделировалось движение 10 млрд материальных точек на протяжении 13 млрд лет. На верхнем кадре каждое яркое пятнышко соответствует галактике
По мере накопления наблюдательного материала стало ясно, что внешний вид галактик тесно связан с их внутренними свойствами — массой, светимостью, структурой звездных подсистем, типами населяющих галактику звезд, количеством газа и пыли, скоростью рождения звезд и др. Казалось, отсюда всего полшага до разгадки происхождения галактик различных типов — все дело в начальных условиях. Если первоначальное протогалактическое газовое облако практически не вращалось, то в результате сферически-симметричного сжатия под действием сил тяготения из него образовывалась эллиптическая галактика. В случае вращения сжатие в направлении, перпендикулярном оси, останавливалось благодаря тому, что тяготение уравновешивалось возросшими центробежными силами. Это приводило к формированию плоских систем — спиральных галактик. Считалось, что сформировавшиеся галактики в дальнейшем не испытывают никаких глобальных потрясений, в одиночестве производя на свет звезды и неспешно старея и краснея по цвету за счет их эволюции. В 5060-х годах прошлого века считалось, что в этом описанном сценарии так называемого монолитного коллапса остается уточнить лишь некоторые детали. Но как только взаимодействие галактик было признано двигателем их эволюции, эта упрощенная картина стала неактуальной.
Секстет Сейферта. Четыре сливающиеся галактики плюс приливный выброс из одной из них (справа внизу) и далекая спиральная галактика (в центре)
Два в одном
Проблема предсказания движения большого числа массивных точек, взаимодействующих по закону всемирного тяготения, получила в физике название задачи N тел. Решить ее можно только методом численного моделирования. Задав массы и положения тел в начальный момент, можно по закону тяготения вычислить действующие на них силы. Полагая эти силы неизменными в течение короткого отрезка времени, легко рассчитать новое положение всех тел по формуле равноускоренного движения. А повторяя эту процедуру тысячи и миллионы раз, можно смоделировать эволюцию всей системы.
В галактике вроде нашей более ста миллиардов звезд. Напрямую рассчитать их взаимодействие не под силу даже современным суперкомпьютерам. Приходится прибегать к разного рода упрощениям и ухищрениям. Например, можно представлять галактику не реальным числом звезд, а таким, какое может осилить компьютер. В 1970-х годах брали всего по 200500 точек на галактику. Но расчет эволюции таких систем приводил к нереалистичным результатам. Поэтому все эти годы шла борьба за увеличение числа тел. Сейчас обычно берут по нескольку миллионов звезд на галактику, хотя в отдельных случаях при моделировании зарождения первых структур во Вселенной используют до десяти миллиардов точек.
Другое упрощение состоит в приближенном расчете взаимного притяжения тел. Так как сила тяготения быстро убывает с расстоянием, притяжение каждой далекой звезды не обязательно вычислять слишком точно. Далекие объекты можно сгруппировать, заменив одной точкой суммарной массы. Эта методика получила название TREE CODE (от англ. tree— дерево, поскольку группы звезд собираются в сложную иерархическую структуру). Сейчас это самый популярный подход, многократно ускоряющий вычисления.
Столкновение галактик NGC 2207 и IC 2163 продолжается уже 40 миллионов лет. В будущем их ждет полное слияние
Но и на этом астрономы не успокоились. Они даже разработали специальный процессор GRAPE, который не умеет делать ничего, кроме расчета взаимного гравитационного притяжения N тел, но зато с этой задачей справляется чрезвычайно быстро!