Проблемы Эволюции

В.Н.Снытников

Астрокатализ как стартовый этап геобиологических процессов. Жизнь создает планеты?

УДК 523.4-52: 573.552

Институт катализа им. Г.К. Борескова СО РАН, Новосибирск E-mail: snyt@catalysis.nsk.su

В кн.: Эволюция биосферы и биоразнообразия. К 70-летию А.Ю.Розанова. М.: КМК, 2006. С. 49-59.

Есть основания предполагать, что абиогенный синтез пребиотического вещества как этап возникновения жизни на Земле имел место в протопланетном околозвездном диске. Развита концепция астрокатализа, который представляет собой комплекс методов и подходов научных исследований в катализе применительно к указанной проблеме, включая математическое моделирование. В вычислительных экспериментах на суперкомпьютерах обнаружено, что сгустки вещества для протопланет первично возникают как уединенные нелинейные волны при развитии гравитационной неустойчивости в двухфазной среде околозвездного диска. В сгустках осуществляются благоприятные условия для гетерогенного каталитического синтеза органических соединений. Обосновывается, что каталитические процессы формируют планеты. Определены физико-химические параметры среды, в которой возможно начинали развиваться химические процессы мира РНК.

Ключевые слова:

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

В естествознании к настоящему времени получен широкий круг сведений об эволюции окружающего мира. Вселенная родилась в результате Большого взрыва 13,7 млрд. лет назад вместе со своим временем, пространством и веществом. Она последовательно проходила этапы инфляции, отрыва вещества от излучения, темной эпохи, рождения галактик и звезд первого поколения, звездного нуклеосинтеза химических элементов, взрывов звезд и формирования в молекулярных облаках звезд следующих поколений вместе с планетными системами. Около Солнца на планете Земля, появившейся 4,56 млрд. лет назад, природные процессы привели к биологической органической жизни.

Одна из наиболее крупных проблем современной науки связана с происхождением жизни во Вселенной и в частности на Земле — когда, где и как? Принципиальные ответы на эти вопросы получены современным естествознанием. Когда? Не ранее появления первых звезд во Вселенной, так как кроме водорода и гелия не было химических элементов жизни — углерода, кислорода, азота. Где? Для земной жизни — в Солнечной системе впервые, но самые загадочные 600 млн. лет от начала ее формирования. В противном случае существовал бы механизм переноса огромных масс биологических соединений (для сухого биологического вещества на Земле не менее 2,5-10¹⁸ г по углероду) от одной планеты к другой или от звезды к звезде, что весьма проблематично. Как? В непрерывном процессе самоорганизации материи от одного этапа к следующему, при ее эволюции.

Найденные ответы требуют от исследователей своей дальнейшей конкретизации. Для земной жизни необходимо сузить указанные временные и пространственные рамки, и вместе с тем проследить и детализировать стадии самоорганизации вещества, которые привели к возникновению жизни. Какие-то из этих стадий тщательно исследованы, к примеру, звездный синтез химических элементов с их космической распространенностью, другие стадии на сегодня только намечены. По-видимому, одним из таких наименее изученных этапов самоорганизации вещества является этап абиогенного синтеза первичных предбиологических соединений, приведший к «миру РНК» (Спирин, 2001).

Абиогенный синтез и самоорганизация

В исследованиях абиогенного синтеза первичного органического вещества на Земле, прежде всего следует ответить на вопрос, является ли эта проблема планетарной (гипотеза Опарина-Холдейна) или космической (гипотеза Аррениуса или «панспермии»). Для ответа могут быть привлечены данные об общих свойствах самоорганизации из различных, в том числе внешне не связанных между собой научных областей.

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

Эффективность отдельного этапа самоорганизации, определяемая как отношение массы новой системы к массе исходного вещества, составляет несколько процентов. Так, для Большого взрыва отношение наблюдаемого вещества к общему количеству гравити-рующей материи оценивается в 4%. Для процесса термоядерного синтеза отношение элементов массивнее гелия к наблюдаемому веществу составляет приблизительно 2% и так далее. При увеличении массы порожденной системы, изымающей энергию или массу из исходной системы выше значений в несколько процентов, происходит разрушение исходной системы. К примеру, в иерархических пищевых цепях действует закон 10% Линдемана предельного отношения хищников к жертвам.

В химических несамоорганизующихся процессах, в частности для целенаправленного синтеза аммиака, отношение целевого продукта ко всему потребленному сырью и использованному веществу при высокоразвитом производстве достигает 30%. Химические технологии тонкого органического синтеза по аналогичному показателю, как правило, не достигают и 1°, а часто имеют на порядки меньшие показатели. Поэтому возникает вопрос: сколько современный химик-технолог должен взять простейшего сырья и материалов, чтобы получить в заводских условиях находящиеся на поверхности Земли

2,5-10¹⁸ г биополимеров состава биологического вещества (CH₂O)₁₀₆(NH₃)₁₆H₃PO₄? Сколько стадий синтеза подобных биополимеров можно запланировать при существующей эффективности синтезов на этих стадиях? Подчеркнем, что «синтез» понимается как появление нового качества, как этап самоорганизации. Если предположить в качестве оценки всего 10 стадий с высокими 10% выходами продукта на каждой, то необходимое количество исходных веществ для синтеза органического вещества достигнет 10²⁸г — больше массы Земли, составляющей 6-10²⁷ г!

Из этой оценки следуют как минимум три вывода. Во-первых, число стадий самоорганизации на химическом этапе эволюции от простых химических соединений к "миру РНК" должно находиться в пределах нескольких единиц. Во-вторых, абиогенный синтез первичных органических соединений должен включать в себя на порядки большую массу простых соединений С, N, О, чем их имеется на поверхности Земли и в ее коре. В-третьих, для высокой производительности синтеза эволюция на первых химических этапах должна быть каталитической с включением твердой фазы. Для самого синтеза должен быть указан реактор подобно тому, как звезды представляют собой реакторы синтеза химических элементов. И поэтому вопрос об абиогенном синтезе первичного преби-отического вещества скорее всего может найти свое решение как астрофизическая или астрохимическая проблема.

Астрофизика, астрохимия и астрокатализ Протопланетный аккреционный диск

Оперировала ли природа столь большими, как 10²⁸ г, массами «сырья» для органических синтезов? На этот вопрос следует ответить утвердительно. В едином процессе формирования Солнца вместе с его планетной системой участвовало исходных веществ порядка 10³³г, что оценивается по массе Солнца. Около 2% их веществ приходится на элементы тяжелее гелия, из которых наибольшая часть падает на С, N, О (табл. 1). Но на Земле недостает элементов С, N, О относительно их космической распространенности порядка 10 масс Земли, а водорода и гелия не хватает еще больше — около 100 масс Земли (Войткевич , 1990). Когда и как протоЗемля потеряла основное количество элементов «жизни» Н, С, N, О? Ясно, что этот вопрос тесно связан с вопросом о происхождении и ранней эволюции планет. Чтобы получить ответ, необходимо рассмотреть про-тозвездно- протопланетную стадию формирования Солнечной системы 4,56 млрд. лет назад, когда масса вещества в протопланетном диске была сравнима с массой Солнца.

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

На стадии, когда масса протозвезды становиться сравнимой с массой Солнца, усиливающееся излучение и звездный ветер сбрасывают окружающее вещество. После этого звезду можно увидеть в оптическом диапазоне. Через протопланетный диск на прото-звезду может проходить масса вещества, сопоставимая с массой звезды. Время формирования звезды солнечной массы — порядка 1 млн. лет. Для звезд типа Солнца диаметр протопланетного диска оценивается в 100-200 AU. По ряду астрофизических моделей (Макалкин, 2003), температура газа на орбите Юпитера достигает 100-200 К, температура газа на орбите Земли — до 1000 К, а его давление — свыше 10^-4 атм.

Таблица 1. Относительная космическая распространенность по числу атомов некоторых наиболее обильных элементов (Войткевич и др., 1990). Число атомов [Si] = 10⁶

Пылинки молекулярного облака имеют характерный диаметр 0,1 мкм и многослойное строение. Их внутреннее ядро размером порядка 10 нм состоит из тугоплавких неорганических соединений на основе соединений кремния, магния, железа (табл. 2), как следующих по распространенности элементов после азота, углерода, кислорода. Ядро покрыто сконденсированными органическими соединениями и гидридами азота, углерода, кислорода, в частности водой, в порядке их летучести. Двигаясь с периферии прото-планетного диска к центру по спиральным траекториям, при повышении температуры пылинки последовательно теряют свои летучие компоненты, освобождая неорганическое ядро. При этом такие частицы и тяжелые органические соединения накапливаются в

аккреционном диске за счет очень сильного центрифужного эффекта относительно водорода и гелия как газов-носителей. Значение центрифужного параметра составляет на орбите Земли величину Z = (m1-m2)V²/kT~2000, где m1 — масса пылинки, m2 — масса молекулы водорода, V — скорость движения, Т — температура газа. Кроме того, как нами было экспериментально показано (Хасин, Снытников, 2005) на материалах земного и метеоритного происхождения с железом и кремнием, более или менее воспроизводящих распространенность элементов в космосе, наночастицы обладают хорошей каталитической активностью в реакции синтеза Фишера-Тропша. Для этого синтеза железо на диоксиде кремния представляет собой классический промышленный катализатор. В этой реакции из водорода и моноксида углерода, следующего по распространенности в космосе газа после водорода и гелия (табл. 2), синтезируются углеводороды, в том числе и сложные с высокой температурой кипения. При появлении на поверхности частиц сложных углеводородов пылинки начинают легко слипаться между собой как комки пластилина. Такое слипание с нарастанием «снежного кома» приводит к увеличению размеров и появлению тел с диаметром 1-10 м. Начиная с этих размеров, тела переходят в режим движения вокруг протозвезды с редкими столкновениями между собой и космическими относительными скоростями, превышающих скорость звука в конденсированных телах. Столкновение тел в этих условиях будет приводить к их дроблению и разрушению.

Таблица 2. Ориентировочные составы газа и неорганической твердой фазы в протопланетном диске, исходя из космической распространенности элементов в Солнечной системе (Войткевич и др., 1990)

Формирование протопланет

Отсутствие столкновительного механизма укрупнения тел до сотен километров оказалось принципиальной трудностью развитых ранее сценариев происхождения планет (Сафронов, 1969). Для решения этой проблемы нами предложен другой возможный механизм, связанный с развитием гравитационной неустойчивости в двухфазной системе газ и твердые тела, диаметр которых порядка 1 м. Число таких тел в Солнечной системе может быть оценено, исходя из массы планет с поправкой на величину газовой компоненты газовых планет-гигантов. Обычная гравитационная джинсовская неустойчивость для газа в протопланетном диске не может развиться из-за высоких температур. В рое тел, двигающихся с космической скоростью вокруг протозвезды, на первый взгляд неустойчивость также не может развиться из-за высоких относительных скоростей. Но при движении ансамбля тел в газе, тормозящем и выравнивающем их относительные скорости, и при некотором критическом значении средней массы тела снова возникают условия для единовременного объединения тел в сгустки. Размеры сгустков определяются тепловой скоростью газа, средней массой тела и числом тел в единице объема. Формирующиеся сгустки сопоставимы по размерам с Солнцем. Давление основных газовых компонентов (водорода и гелия) под воздействием самосогласованного гравитационного поля увеличивается, что обеспечивает высокое абсолютное содержание в сгустках соединений кислорода, азота и углерода и нелетучих веществ.

Области развития гравитационных неустойчивостей предположительно и были местом зарождения крупных, многокилометровых первичных тел и протопланет. Вместе с тем, в этих сгустках с подвижной твердой фазой мог осуществляться каталитический абиогенный синтез органических соединений. Следы этих процессов следует искать в метеоритах и астероидах, на спутниках планет, там, где не было геологической эволюции. Так как время изменения содержания химических соединений в реакциях органического синтеза редко превышает часы или дни при нормальных условиях, то темпы химической эволюции на этой стадии были очень высокими — на уровне нескольких лет. Соответственно «мир РНК» предположительно существовал именно в этих условиях. Если же гипотеза (Розанов, 1996) о следах колоний микроорганизмов в метеоритах получит свое подтверждение, то и жизнь в Солнечной системе зародилась на этом временном этапе. Тогда последующие периоды могут рассматриваться как борьба живых существ за выживание в условиях катастрофических изменений внешних условий, в том числе и при формировании планет, за счет усложнения своей организации и создания сообщества организмов, биосферы, в условиях потерь подавляющей части первично синтезированного органического вещества.

Итак, стадия «астрокатализа» для первичного абиогенного синтеза основной массы органических соединений соответствует этапу формирования крупных тел в солнечной системе. Главный результат наших исследований по астрокатализу сводится к выводу, что абиогенный синтез первичных органических соединений проходил непосредственно при формировании первичных тел и протопланет при развитии гравитационной коллективной неустойчивости с одновременным объединением множества малых тел (Сныт-ников и др., 1996). Иначе говоря, каталитические процессы формируют планеты Солнечной системы.

Вычислительные эксперименты на суперкомпьютерах

Один из наиболее подходящих инструментов для изучения этапа астрокатализа -математическое моделирование с проведением вычислительных экспериментов на суперкомпьютерах. Для этого имеются развитые математические модели физико-химической динамики многофазных реагирующих сред с учетом гравитации в виде систем уравнений математической физики (Снытников и др., 2004). Математическая модель предполагаемого гравитационного механизма формирования протопланет при развитии неустойчивости включает в себя уравнение Власова-Лиувилля динамики малых тел, уравнения газовой динамики и уравнение Пуассона для самосогласованного гравитационного поля. Решение таких уравнений производится на сетке методом крупных частиц в ячейке для уравнения Власова, методом частиц для газодинамических уравнений и комбинированным методом или быстрым преобразованием Фурье для уравнения Пуассона. Сложность математической модели для вычислительных экспериментов состоит в ее нестационарности и пространственной трехмерности. Более того, для уравнения Власова должны быть получены нестационарные решения в шестимерном фазовом пространстве. В задачах на развитие неустойчивости отсутствует непрерывная зависимость решения по параметрам. Такие задачи относятся к некорректным, по Адамару, и требуют для своего развития особых подходов, ориентированных в первую очередь на выполнение фундаментальных физических законов сохранения. Тем не менее, основные начальные трудности были преодолены. Достоверность созданных прог для численного моделирования подтверждена комплексом перекрестных методов проверки получаемых численных решений, многолетней по длительности и трудоемкости.

Некоторые результаты расчетов представлены на рисунках, на которых приведены распределения логарифма плотности твердых тел в экваториальной плоскости вокруг протозвезды для (рис. 2) и вокруг новой образовавшейся структуры (рис. 3). Моменты времени, для которых приведены плотности вещества, соответствуют поздней стадии гравитационной неустойчивости, вышедшей на уровень нелинейного насыщения. На этой стадии в экваториальной плоскости диска на фоне средней плотности образовались разлетающиеся сгустки вещества из свободно двигающихся тел (рис. 2), которые связываются общим самосогласованным гравитационным полем.

Такие сгустки из твердых частиц в плотной атмосфере водорода и гелия представляют собой уединенные волны — солитоны, которые могут двигаться по ходу вращения и против него, к центру притяжения и от него на периферию. Эти волны формируется коллективным движением частиц. Частица захватывается в волну, находится в ней некоторое время и уходит. Другие частицы могут просто пролететь мимо или образовать систему колец, спутников вокруг сгустков, как на рис. 3.

В этих сгустках давление газа на какие-то промежутки времени и в зависимости от расстояния до протозвезды превышает десятки атмосфер. Высокое давление водорода и гелия означает умеренные температуры в волне, так как эти газы обладают высокой теплопроводностью. Большое количество гелия (свыше 20%) поставляет энергию для эндотермических реакций и снимает тепло экзотермических реакций. Волна действует как великолепный химический реактор, устроенный по типу промышленных реакторов с «псевдоожиженным» катализатором. Это один из самых эффективных, но вместе с тем

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

Дальнейшая эволюция вещества в уединенной волне достаточно очевидна. При увеличении массы органических соединений волна способна схлопнуться в сгусток вещества. В ближней к Солнцу зоне сгусток двигается, теряя водород, гелий и легкую органику под действием солнечного ветра и излучения. В дальней зоне Юпитер и холодные внешние планеты сохраняют эти газы. Затем сгусток превращается в планету, вступая в геологическую или при подходящих условиях в геобиологическую эволюцию. Основная масса органического вещества и метана из зоны первичного синтеза разрушается, попав на Солнце и развеиваясь в космосе. Однако тяжелые и сложные органические соединения при огромной своей концентрации могли сохраниться в дальнейших катаклизмах и стать основой для возникновения и питания биологического сообщества. Причем естественный отбор молекул с автокаталитическими свойствами, по предположению В.Н. Пармона (2002), начинается на стадии потерь и снижения количества «пищи». Важно, что отбор на химической стадии может происходить не только в обычных земных условиях, но и в космосе.

Еще один важный вопрос состоит в том, насколько предопределено образование планет. В изучаемой модели развития гравитационной неустойчивости самоорганизация вещества развивается детерминированно в одном направлении (Вшивков, 2003). Детерминированность здесь понимается как неизбежность возникновения волн плотности на нелинейной стадии развития неустойчивости. Но место возникновения уединенной волны плотности зависит от начальных, принятых в расчетах, функций распределения тел по скоростям. Поэтому в использованной модели координаты возникновения солитонов представляют собой случайные величины с некоторой вероятностью в зависимости от удаления от центрального тела. При возникновении волн плотности на орбите Венеры температура среды определит условие испарения воды и органических соединений. На орбитах Марса при низком давлении, невысоком гравитационном поле и относительно низких температурах осуществляются условия возгонки и потери органических веществ и воды в окружающее пространство. По-видимому, на орбитах Земли при ее массе утря-чивается восстановительная атмосфера из водорода и гелия, но удерживаются вода и другие газы. В дальнейшем необходима детализация расчетной модели с учетом дополнительных физико-химических процессов, что связано с важностью определения возможной зоны распространения жизни в Солнечной системе. В частности, входит ли Марс в эту зону с возможностью проверки получаемого ответа автоматическими спутниками?

Реактор абиогенного синтеза предбиологических соединений

Таким образом, для абиогенного синтеза первичных органических соединений наиболее подходят условия в «космическом каталитическом реакторе» протопланетного диска со следующими характеристиками. Это реактор с псевдоожиженной кипящей твердой фазой и восстановительной водород-гелиевой атмосферой, в которой давление газа до или свыше 10 атмосфер. Размеры «гранул» твердой фазы составляют порядка 1-10 м. Эта фаза представляет собой соединения, включающие в качестве основных SiO₂ - MgO - Fe. У нее каталитически активная огромная поверхность, на которой на определенных

расстояниях от протозвезды конденсируются вода и другие соединения. Размеры реактора на разных временных этапах — от 10^{7 км} (0.1 AU) до десятков диаметров Солнца (0,01AU). Ввод энергии в зону реактора осуществлялся излучением протозвезды при нагреве поверхности протопланетного диска. Охлаждение зоны синтеза при экзотермических реакциях происходило посредством водорода и гелия. Время существования реактора синтеза химических соединений составляет порядка 10 лет. Конечное состояние по химическим соединениям - высокомолекулярные органические соединения, Н2О и другие гидриды элементов, возможно, «мир РНК». Конечное состояние по физическим условиям — многокилометровые тела в восстановительной атмосфере. Число таких структурных образований как реактор синтеза химических соединений определяется деталями развития гравитационной неустойчивости в двухфазной среде диска. На начальной стадии развития неустойчивости это число может составлять несколько единиц и более.

Место абиогенного синтеза пребиотических соединений на временной шкале основных этапов самоорганизации (эволюции), приведших к появлению жизни на Земле, представляется в следующем виде (табл. 3).

Для получения дальнейших данных по астрокатализу необходимы: 1) характеристики околозвездных дисков, которые могут быть получены посредством наблюдательных астрофизических методов, 2) сравнительная планетология, изучение Солнечной системы, в том числе активным зондированием космическими аппаратами, 3) лабораторное экспериментальное моделирование процессов, связанных с авто- и каталитическим синтезом предбиологических соединений, 4) исследование реакций синтеза нуклеотидов и других соединений для «мира РНК» в условиях «астрокатализа» (высокие давления, умеренные температуры, восстановительная водородно-гелиевая атмосфера, тепловое излучение, реактор кипящего слоя), 5) математическое моделирование физико-химических процессов в околозвездных дисках и формирования протоЗемли и 6) другие направления, исследования и результаты.

Таким образом, проблема абиогенного синтеза первичного пребиотического вещества на поверхности Земли может эффективно решаться с применением методологии и подходов катализа. Каталитические исследования предполагают изучение в лабораторных условиях химических реакций, в том числе на установках in situ с привлечением физических методов, синтез и исследование катализаторов, создание и изучение различных типов каталитических реакторов с их массо- и теплообменом, разработка и изучение всего технологического процесса синтеза, включая процессы дезактивации и регенерации катализатора. Каждый из этих этапов связан с другими, причем часто эти связи носят нелинейный характер. Поэтому широко используются методы математического моделирования. Комплекс подходов и методология исследований для изучения химического этапа эволюции вещества в протопланетном диске могут рассматриваться как область астрокатализа. Кроме того, в настоящей работе сделана попытка обосновать утверждение, по которому первичные органические соединения синтезируются в доп-ланетной околозвездной среде, а планеты формируются в местах каталитического синтеза органических соединений.

Исследования в этом направлении были поддержаны академиками Н.Л. Добрецовым, Г.А. Заварзиным и В.Н. Пармоном по программе Президиума РАН «Происхождение и эволюция биосферы» (№18-2), грантом Минобрнауки РНТ.2.1.1.1969 и НШ-6526.2006.3.

В заключение автор выражает искреннюю благодарность члену- корреспонденту РАН А.Ю. Розанову за многочисленные и исключительно полезные обсуждения, а также за поддержку, оказанную автору и коллективу сотрудников, которые совместно работают по рассмотренной проблеме в институтах СО РАН.

Литература

Войткевич Г.В., Кокин А.В., Мирошников А.Е., Прохоров В.Г. 1990. Справочник по геохимии. М.: Недра. 480 с.

Вшивков В.А., Никитин С.А., Снытников В.Н. 2003. Исследование неустойчивостей бесстолкнови-тельных систем по стохастическим траекториям//Письма в ЖЭТФ. Т. 78. № 6. С. 810-815.

Макалкин А.Б. 2003. Проблемы эволюции протопланетных дисков / Современные проблемы механики и физики космоса. М.: Физматлит. С. 402-446.

Пармон В.Н. 2002. Пребиотическая фаза зарождения жизни // Вест. РАН. Т. 72. № 11. С. 976-983.

Розанов А.Ю. 1996. Цианобактерии и, возможно, низшие грибы в метеоритах // Соросовский образовательный журнал. № 11. С. 61-65.

Сафронов B.C. 1969. Эволюция допланетного облака и образование Земли и планет. М.: Наука. 244с.

Снытников В.Н. Вшивков В.А., Кукшева Э.А. и др. 2004. Трехмерное численное моделирование нестационарной гравитирующей системы многих тел с газом // Письма в астрономический журнал. Т. 30. № 2. С. 146-160.

Спирин А.С. 2001. Биосинтез белков, мир РНК и происхождение жизни//Вести. РАН. Т. 71. №4. С. 320-328.

Хасин А.А., Снытников В.Н. 2005. Особенности состава продуктов синтеза Фишера-Тропша на материале метеорита Царев // Труды Межд. семинара «Происхождение и эволюция биосферы», 26-29 июня 2005 г. Новосибирск: Институт катализа им. Г.К. Борескова СО РАН. С. 158-159.

Snytnikov V.N., Vshivkov V.A., Parmon V.N. 1996a. Solar nebula as a global reactor for synthesis of prebiotic molecules // 11^th Intern. Congr. on the origin of life. Orleans. P. 65.

Snytnikov V.N., Vshivkov V. A., Parmon V.N. 1996b. Mathematical modeling of non-stationary physico-chemical processes in natural protoplanetary catalytic reactor // 13 Intern. Congr. on chemical reactors. P. 2. Novosibirsk. Abstr. P. 226-227.