А. Р. Оганов, профессор «Химия и жизнь» №10, 2011
Алгоритм
Совсем недавно физики могли только мечтать о том, чтобы заложить в компьютер химический состав какого-нибудь соединения и получить распечатку с наиболее энергетически выгодным расположением атомов в его кристаллической решетке, а также со значениями свойств этого вещества, такими, как твердость, модуль упругости или ширина запрещенной зоны. В 1988 году Джон Мэддокс, тогда главный редактор журнала «Nature», сказал: «Один из продолжающихся скандалов физики состоит в том, что невозможно предсказать структуру даже простейшего кристалла, если исходить только из его химического состава... Что касается таких кристаллических веществ, как лед, (то) это до сих пор считается запредельным для смертных». Даже возросшая за прошедшую четверть века мощность вычислительной техники не приблизила нас существенно к решению этой задачи: одной мощности оказалось недостаточно, потребовались принципиально новые теоретические подходы. Например, предсказание структуры средней сложности на лучших сегодняшних суперкомпьютерах методом перебора займет порядка 1020 лет. Очевидно, что без концептуального рывка эта задача нерешаема.
Мы начали работу над своим методом в 2004 году. Тогда наша группа и коллеги из Токийского технологического института независимо друг от друга обнаружили, что соединение MgSiO3 под высоким давлением принимает неизвестную до того структуру, которая определяет необычные свойства нижнего слоя мантии Земли и это существенно изменило модель земных недр. Наша работа как раз и показала, насколько точными и полезными для целенаправленного проведения экспериментов могут быть теоретические расчеты.
Подобные расчеты делают, сравнивая значения энергии различных структур при заданных давлении и температуре. Та структура, что имеет меньшую энергию, и будет при этих условиях стабильной. Такие расчеты, основанные на квантовой механике, позволяют рассчитать свойства материалов, термодинамику и кинетику фазовых превращений. В природе мириады атомов, совершая колебания, методом проб и ошибок находят состояние с минимумом энергии. У исследователя нет возможности провести расчет в разумном диапазоне времен — современные компьютеры дают возможность следить за динамикой достаточно большого массива атомов лишь на протяжении 10–10 с. Нет и возможности перебрать все варианты расположения атомов. Поэтому возникает задача разработать алгоритм поиска энергетически выгодных структур, который позволит значительно сократить объем расчетов. Мы придумали алгоритм, основанный на идеях эволюции.
Представим себе, что имеется горная страна, испещренная долинами, ущельями и хребтами. Есть исследователи, которые хотят найти самую низколежащую долину, а топографической карты у них нет. Выбрасываясь десантом в разных точках страны, каждый из них спускается в ближайшую долину (локальный минимум энергии). После того как все исследователи достигли своих минимумов, можно сравнить, кто находится ниже. Точки неглубоких минимумов исключают из рассмотрения. А те области, в которых были достигнуты самые глубокие минимумы, исследуют более подробно путем скрещивания (нового исследователя сбрасывают между двумя глубокими минимумами) или мутации (исследователь «выпрыгивает» из своей долины в соседнюю). Фокусируя внимание на все более энергетически выгодных областях пространства поиска, расчет довольно быстро находит самый глубокий из всех возможных минимумов энергии. Он-то и будет отвечать термодинамически наиболее выгодной фазе. Другие глубокие долины могут соответствовать каким-то метастабильным структурам.
Поначалу мы столкнулись с огромными трудностями в создании такого алгоритма, однако после года упорной работы к 2005 году нам удалось придумать метод, надежность и быстродействие которого удивили не только нас, но и практически всех наших коллег. Чтобы убедиться в том, что наш метод работает, ко мне в лабораторию приезжали люди из самых разных стран от Китая до Канады. Эти первые успехи были достигнуты вместе с моим аспирантом Колином Глассом. Начиная с 2006 года программу опробовал мой постдок Ма Янмин (сейчас он профессор в Цзилиньском университете, КНР), совместно с которым мы предсказали много новых фаз элементов и неорганических соединений. В 2007 году к нам присоединился А. Ляхов, который стал главным программистом проекта, в 2009 году — Чжу Цянь, талантливый аспирант из Китая. Сейчас, к 2011 году, опубликовано уже более полусотни работ с предсказаниями строения новых фаз элементов и разнообразных соединений с использованием этого алгоритма. Программой пользуются более 600 ученых по всему миру, а на Первых соревнованиях по слепому поиску структур неорганических соединений программа значительно опередила всех конкурентов как по скорости, так и по надежности вычислений.
Для расчета нужен минимальный набор параметров: число атомов каждого сорта, давление, температура и число структур, рассматриваемых на каждом шаге эволюционного процесса. Если известны параметры элементарной ячейки (а их определять по рентгенограммам легко, в отличие от расположения атомов в элементарной ячейке), то расчет облегчается. Когда машинного времени достаточно много, программа сама может найти не только структуру, но и состав всех стабильных соединений (при том что исходно заданы только типы входящих в него атомов), способна рассчитать структуру низкоразмерных объектов — наночастиц, поверхностей, а также определить расположение молекул в структуре молекулярного кристалла.
Этот метод особенно полезен в двух областях исследований: в поиске новых материалов с заданными свойствами и в физике/химии высоких давлений. Расскажем о трех ярких примерах успеха такого расчета.
Ионный бор
Наверное, ни один элемент не доставил своим исследователям столько головной боли, как бор. Вот краткая история поиска структуры этого вещества.
Бура, натриевая соль борной кислоты, — Na2[B4O5(OH)4]*8H2O — известное с древности вещество. Его название, давшее имя и самому элементу, с арабского переводится как «белый». Именно из буры в 1702 году Вильгельм Хомберг получил белый порошок «успокоительной соли», известной сегодня как метаборная кислота, HBO3. Спустя век, в 1808 году, она и послужила исходным веществом, из которого практически одновременно получили бор Жозеф Луи Гей-Люссак совместно с Луи Жаком Тенаром (21 июня) с одной стороны Ла-Манша и Хемфри Деви (30 июня) — с другой. Французы восстанавливали бор из кислоты с помощью калия, а англичанин — электролизом. Это соревнование закончилось, по сути, вничью: как теперь становится ясно, обе группы получили вещество, в котором чистого бора было от силы половина и свойства которого не соответствовали элементарному бору. Так, вещество Деви было проводником, а чистый бор — полупроводник. Среди тех, кто сумел это доказать, был еще один великий химик, француз Анри Муассан. В 1895 году он восстановил бор из его оксида В2О3 с помощью магния. Но и это вещество оказалось еще далеко не чистым. Есть мнение, что в 1909 году американский химик Езекиель Вейнтрауб синтезировал бор чистотой 99%, однако это тоже не очевидно. Судя по всему, чистые аллотропные модификации бора были получены не ранее 1957 года.
Разговор же об аллотропах бора возник спустя полвека после первых попыток получить хоть что-то содержащее чистый бор. Так, в 1857 году Фридрих Вёлер и Анри Сен-Клер Девилль нагрели оксид бора с алюминием и получили три формы бора, различающиеся блеском и твердостью. По аналогии с углеродом их назвали алмазо-, графито- и углеподобным бором. Последний имел те же свойства, что и вещество, сделанное Гей-Люссаком и Тенаром. И ни один из этих «боров» не был чистым бором. Впоследствии еще множество исследователей стали жертвами коварства бора. А связано оно с его необычной химией и сложными кристаллическими структурами.
У атома бора всего три валентных электрона, как и у алюминия. Этого недостаточно для образования четырех ковалентных связей и заполненной электронной конфигурации атома. Поэтому, как алюминий и другие элементы I–III групп, бор «хочет» быть металлом и коллективизировать такие электроны, но это оказывается невыгодным — они слишком сильно связаны с ядром атома. В результате бор выбирает компромиссное решение — аллотропы с полупроводниковыми свойствами и очень сложными кристаллическими структурами. Поскольку они образуются из-за недостатка валентных электронов, бор крайне чувствителен к примесям, которые могут добавить или отнять электроны. Поэтому трехвалентный бор способен образовывать соединения с самыми экзотическими формулами, нарушая все правила валентности: YB66, B6O, NaB15, B15P2, B13C2, AlC4B40, B50N2 — вот далеко не полный список. Борид иттрия можно считать одним из чемпионов по сложности структуры: в элементарной ячейке его кристаллической решетке находится 1588 атомов! Неудивительно, что бор до сих пор остается загадочным элементом. В самом деле, как знать, с чем связана структура вещества, оказавшегося в руках исследователя, — со свойствами чистого бора или какого-то экзотического соединения, где на сотню атомов бора приходится один-два атома другого элемента?
Как бы то ни было, сейчас известно шестнадцать кристаллических аллотропов бора, причем большинство из них скорее всего содержат примеси, а структуры установлены только для четырех. Основные открытия этих аллотропов были совершены в два этапа: в 1957–1965 годах (когда шел поиск сверхтвердых материалов) и в XXI веке. Первой оказалась фаза Т-50 с тетрагональной кристаллической решеткой (см. рис. 1 на врезке), в элементарной ячейке которой находится 50 атомов. Впервые эту фазу получили в Корнелловском университете и в компании «Дженерал электрик» в 1943 году, а структуру расшифровали в 1951-м. Увы, дальнейшие исследования показали, что это карбид или нитрид бора с формулами B50C2 и В50N2. В 1957 году наконец нашли первую фазу действительно чистого бора. Это оказалась ромбоэдрическая бета-фаза (см. рис. 1 на врезке). В ее элементарной ячейке содержится 106 атомов — это стало ясно лишь в 1963 году, после чего ее стали называть бета-В106.
В 1958 году была открыта альфа-фаза с ромбоэдрической же ячейкой из 12 атомов (альфа-В12), а в 1960-м — фаза Т-192, крайне сложную структуру которой определили лишь в 1979-м. Все эти решетки представляют собой более или менее плотное кружево из борных икосаэдров — двадцатигранников с треугольным гранями.
Уже в 1971 году чистоту фаз Т-50 и бета-В106 поставили под сомнение: опыты с получением бора методом химического осаждения из пара на танталовую проволоку (при этом вещество хорошо защищено от попадания примесей), дали лишь альфа-фазу и Т-192. Бета-В106 получался изредка, а Т-50 — никогда.
В середине шестидесятых годов компания «Дженерал электрик» после успехов в синтезе алмазов и кубического нитрида бора начала исследовать поведение бора при высоком давлении. Так, в 1964 году один из создателей боразона Роберт Уэнторф обнаружил, что при давлении более 10 ГПа (примерно 100 тыс. атмосфер. — Примеч. ред.) и температуре 1800–2300 К аморфный бор и бета-фаза превращаются в новую модификацию. Несмотря на то что Уэнторфу удалось получить дифракционную картинку, измерить плотность и электропроводность — а это был предел мечтаний для исследователей высокого давления того времени, — научное сообщество не восприняло всерьез его результаты, хотя они и были опубликованы в «Science». За 43 года статью процитировали лишь шесть раз, да и то в работах, имеющих к исследованию бора косвенное отношение. Более того, полученные Уэнторфом картинки рентгеновской дифракции — а это основа для идентификации фазы — были удалены из соответствующих баз данных. Возможная причина такого недоверия: Уэнторфу не удалось показать, что он имел дело с чистым бором (его работа не содержала химического анализа исходных веществ и продуктов превращения), не удалось и расшифровать кристаллическую структуру или хотя бы определить форму элементарной ячейки. Как показали дальнейшие события, недоверие было напрасным.
Вторая волна интереса к бору началась в 2001 году, когда обнаружилась сверхпроводимость у столь обычного вещества, стоящего на полке чуть ли не в каждой химической лаборатории, как диборид магния. С самого начала было ясно, что она связана с графитоподобной укладкой атомов бора в этом соединении, и вскоре бор стали сдавливать в надежде обрести сверхпроводящую фазу этого элемента. Действительно, в 2001 году было обнаружено, что бета-фаза при давлении 160 ГПа переходит в металлическое состояние, а при 250 ГПа температура перехода в сверхпроводящее состояние достигает 11,2 К. Так, ставки на высокотемпературную сверхпроводимость плотного бора не оправдались, но теоретический вопрос остался: а что за фаза обеспечила металлические свойства?
Как оказалось, если давить бета-бор без нагрева, то при давлении более 100 ГПа он переходит в аморфное состояние. Это означает, что бета-фаза становится нестабильной, но перейти в другую кристаллическую фазу ей не позволяет кинетика, то есть недостаточна энергия колебаний атомов, чтобы преодолеть барьер, разделяющий две фазы. Нагрев образца лазером непосредственно в процессе сжатия привел к превращению бета-фазы в Т-192 уже при 10 ГПа. Так попутно выяснилось, что отнюдь не примеси отвечают за образование бета-В106: это одна из решеток чистого бора, а Т-192 заняла на диаграмме состояния бора свое место фазы высокого давления и температуры.
Тем временем продолжались дебаты о том, какая фаза, альфа или бета, устойчива при нормальных условиях. В 2007–2009 годах три группы независимо и разными методами рассчитали энергии этих структур исходя из первых принципов, и оказалось, что именно бета-В106 — наиболее устойчивая, а более простая альфа-В12 имеет повышенную энергию.
В 2004 году вторая волна исследований бора привела к еще одному интересному результату: Чен Цзюхуа из Университета Стони Брук и В. Л. Соложенко, работающий в университете Пари-Нор, независимо друг от друга обнаружили новую фазу бора при давлении 10–12 ГПа и температурах выше 1500 К. Однако установить структуру этой фазы им не удалось.
Тогда-то на сцену и вышли математические методы: Чен в 2006 году обратился ко мне с предложением о сотрудничестве. Всего за один день нам удалось решить задачу и получить правильную структуру: расчетные дифракционные картинки вполне удовлетворительно совпали с экспериментальными. Для чистоты исследования я решил не брать экспериментальные картинки; сопоставление провели после окончания расчета структуры, и оно оказалось крайне убедительным. Хотя результат мы получили быстро, опубликовать его удалось лишь через два года — уж очень необычной была структура этой фазы.
Это псевдокубическая структура, составленная из двух типов борных кластеров — икосаэдров и пар из атомов бора, которые лежат в порах между икосаэдрами. Благодаря этому упаковка атомов бора в такой структуре стала гораздо плотнее, чем у любых других фаз этого элемента. Более того, по сути, это частично ионное соединение, хотя и составлено из атомов одного и того же типа: кластеры обменивались между собой зарядами; в среднем каждая пара отдавала икосаэдру по половине электрона. Результат: локализация электронов, отсутствие проводимости и сильное поглощение в инфракрасной области. А при деформации первыми рвутся как раз связи между наиболее заряженными атомами. Эту фазу с 28 атомами в элементарной ячейке мы назвали гамма-В28. Поскольку она прекрасно сохраняется и после снятия давления, удалось измерить твердость, которая оказалась наивысшей из всех фаз бора — 50 ГПа. Что касается открытия Уэндорфа, то, как показало сравнение дифракционных картин, он держал в руках, вероятнее всего, именно эту, плотную фазу бора (возможно, с примесями других фаз).
Открытие гамма-фазы отнюдь не завершает исследований бора при высоких давлениях, но несколько усложняет их. Дело в том, что совсем недавно предполагалось существование при давлении выше 74 ГПа совершенно новой фазы этого элемента, в структуре которой отсутствуют привычные икосаэдры. Она должна обладать и металлической проводимостью, и переходить в сверхпроводящее состояние. Однако теперь стало ясно, что ее получению может помешать гамма-фаза, которая, согласно расчетам, останется устойчивой примерно до 90 ГПа.
Прозрачный натрий
Современная техника для экспериментов в области высоких давлений позволяет так сдавливать образцы, что внутренние орбитали электронов начинают перекрываться. Это может до неузнаваемости изменить свойства вещества, причем совсем не так, как следует ожидать, исходя из общих соображений. Например, считается, что по мере увеличения давления изолятор становится проводником. Типичный пример — появление металлической проводимости у водорода. А может ли металл при высоком давлении, наоборот, стать изолятором? Положительный ответ следует из расчета свойств сверхплотных фаз таких типичных металлов, как литий или натрий. Именно натрий стал объектом новой серии наших расчетов и последовавших за ними экспериментов, и именно с ним в наиболее полной мере проявились неожиданные новые эффекты.
Этот металл при нормальных условиях обладает объемноцентрированной кубической решеткой (см. рис. 1 на врезке). Под давлением 65 ГПа она превращается в более плотную гранецентрированную. Свежие эксперименты, проведенные в 2002–2008 годы, показали, что дальнейшее повышение давления вплоть до 160 ГПа дает быструю смену нескольких фаз, причем структуру некоторых из них установить пока не удалось. Наш расчет показал: выше 200 ГПа должна существовать новая фаза натрия и она — неметалл. Более того, она должна быть прозрачной! Предсказание было настолько необычным, что пришлось убеждать экспериментаторов проверить его.
Сдавливая натрий в алмазных наковальнях, Михаил Еремец из Института Макса Планка в Майнце проследил всю известную последовательность структур в интервале от 100 до 160 ГПа. А при 200 ГПа, когда объем образца уменьшился в пять раз по сравнению с исходным, натрий действительно стал прозрачным, и это было заметно невооруженным глазом. Что может означать подобное изменение свойств?
Согласно зонной теории металлов, уровни электронов, имеющиеся в изолированных атомах, при объединении последних в кристалл расширяются и превращаются в зоны. То есть электрон одного атома может, приобретя небольшое количество энергии, занять уровень электрона другого атома, если тот свободен. А между зонами в изоляторе существует большое расстояние. Чтобы его преодолеть, электрону нужно получить энергию, превышающую некоторое значение, равное ширине запрещенной (для нахождения электронов) зоны. В металле энергия кванта видимого света полностью поглощается кристаллом и расходуется на переход электрона из валентной зоны в зону проводимости. Поэтому металл непрозрачен. А в изоляторе зона проводимости отделена от последней заполненной зоны. Если энергии кванта света недостаточно, чтобы электрон перепрыгнул через запрещенную зону, то свет поглощен не будет — вещество окажется прозрачным. Для обеспечения прозрачности ширина запрещенной зоны должна превышать 1,8 эВ, а для бесцветности — 3,1 эВ. В очень тонком образце ширина зоны для перехода к прозрачности может быть и несколько меньше.
То, что натрий стал прозрачным при 200 ГПа, означает: атомы перестроились в какую-то необычную неметаллическую структуру, где ширина запрещенной зоны оказалась неожиданно большой. Какова эта структура?
Расчет показал, что при давлении от 320 до 1000 ГПа наименьшей энергией обладает двойная гексагональная плотно упакованная (ГПУ) структура. Обычная ГПУ-структура состоит из плотных слоев, которые стыкуются так, что прямо над ямками одного слоя находятся атомы следующего слоя, и в результате получается самая плотная упаковка твердых шаров (см. рис. 2 на врезке). Для этого в двойной ГПУ-структуре соотношение высоты элементарной ячейки к ее длине должно быть 3,266. А в структуре прозрачного плотного натрия это расстояние при 320 ГПа оказалось в два раза меньше — 1,391. Атомы перестали быть твердыми шарами и вошли друг в друга почти наполовину! Так, кратчайшее расстояние между натриями в этой решетке оказалось всего 1,89 Å при 300 ГПа, в то время как радиус внешней 3s-орбитали у него 1,71 Å, а ионный радиус — 1,02 ангстрема. В общем, такое перекрытие внутренних орбиталей разных атомов натрия привело к заселению p- и d-орбиталей и возникновению необычного состояния s-р-d-гибридизации. Расчет показал, что у натрия полностью исчезли свободные электроны, обеспечивающие металлическую связь, и все они оказались локализованы в пустотах решетки. Фактически получилась решетка ионного соединения, состоящая из положительно заряженных ионов натрия и отрицательно заряженных скоплений электронов. Величина запрещенной зоны у нее оказалась немалой — 1,5 эВ при 200 ГПа, а при 600 ГПа — и вовсе 6,5 эВ. Следовательно, типичный металл может стать очень сильным изолятором, если, конечно, удастся достичь такого высокого давления.
Сияющий углерод
Создание технологии синтеза искусственных алмазов отнюдь не остановило исследования углерода при высоких давлениях. В частности, возник вопрос: а что получится, если сдавливать графит при низких температурах и без катализаторов? В ходе этих опытов в 1963 году было обнаружено, что при комнатной температуре при давлении выше 14 ГПа графит испытывает некое структурное фазовое превращение, что сказывается и на электропроводности, и на коэффициенте преломления, и на дифракции рентгеновских лучей, и на других структурно-чувствительных свойствах. Эта фаза оказалась прозрачной и очень твердой: при попытках ее дальнейшего сжатия алмазная наковальня разрушалась. Однако за сорок с лишним лет кристаллическую структуру так и не удалось установить. Поначалу исследователи предполагали, что это какая-то гексагональная фаза, представляющая собой нечто среднее между структурами графита и алмаза. Но предложенные на основе такой идеи структуры не удалось подтвердить экспериментальными данными.
В 2006 году, тестируя наш метод предсказания структур, я с удивлением обнаружил новую, крайне любопытную структуру углерода, которая лишь ненамного уступала алмазу в устойчивости. В 2009 году мы показали, что эта структура полностью описывает экспериментальные данные по «сверхтвердому графиту» и довольно легко получается из графита путем образования дополнительных ковалентных связей между его слоями. Фаза получила название «М-углерод», поскольку ее решетка с точки зрения кристаллографии описывается моноклинной ячейкой (см. рис. 1. на врезке). Расчетная же твердость оказалась сопоставима с алмазом и гораздо выше, чем у кубического нитрида бора (83 ГПа против 90 и 47 ГПа соответственно). Причина такой высокой твердости — в особенностях кристаллического строения. В процессе смятия графитовых плоскостей и их сближения увеличивается число связей, образованных электронами, которые претерпели sp3-гибридизацию. А чем больше таких связей, тем тверже будет материал. Тот факт, что без нагрева графит переходит в М-углерод, а не напрямую в алмаз, неудивителен, если вспомнить о структурном родстве графита и М-углерода. Расчет же показал, что М-углерод должен быть прозрачным, бесцветным изолятором: его ширина запрещенной зоны превышает 5 эВ. Измерения электропроводности экспериментальных образцов М-углерода, впрочем, свидетельствуют, что помимо этой фазы в них есть и остатки графита.
А возможно ли в принципе создание фаз, более плотных и более сверкающих, чем алмаз? Напомним, что среди рекордов, принадлежащих алмазу, — наибольшее среди всех известных веществ число атомов в единице объема. Да, расчет свидетельствует: такие углеродные фазы могут существовать. Чтобы в этом убедиться, мы провели поиск наиболее плотных фаз углерода. Расчет дал все кристаллические фазы — алмаз, гексагональный алмаз (лонсдейлит), М-углерод. Однако наиболее плотными, согласно расчету, оказались три неизвестные до сих пор фазы. Чем они замечательны? Во-первых, своей плотностью — у двух она на 3,2% больше, чем у алмаза, а у третьей — на 1,1%. Во-вторых, твердостью, которая ниже расчетной алмазной менее чем на 10%. Все эти фазы — прозрачные диэлектрики, а в одной из них ширина запрещенной зоны достигает 7,3 эВ — максимальное значение среди всех модификаций кристаллического углерода. Другая фаза отличается большим показателем преломления: в области видимого спектра он меняется от 2,65 до почти 2,8, тогда как у алмаза находится в районе 2,25. Это означает, что кристаллы такой фазы, будучи ограненными, станут сверкать лучше, чем бриллианты, и, кроме того, засияют всеми цветами радуги: сильная зависимость показателя преломления от длины волны приводит к большой дисперсии света. Такое качество может пригодиться не только в ювелирном деле, но и при изготовлении оптических приборов. Вообще, все три предсказанные фазы по значению показателя преломления и по дисперсии света сильно обгоняют алмаз.
А возможен ли синтез таких фаз? Расчет, казалось бы, не свидетельствует в пользу подобной идеи: энергия таких структур высока, в пересчете на один атом она получается на 0,89–1,14 эВ больше, чем у алмаза. Однако известны углеродные фазы с подобной высокой энергией. Это фуллерены (различие 0,44–0,72 эВ с графитом) и аморфный углерод (превышение 0,7–0,99 над алмазом), которые вполне существуют при нормальных условиях. Расчет же показывает, что и новые фазы, если их получить, смогут существовать сколь угодно долго при атмосферном давлении. Главное — выбрать путь синтеза, который позволит миновать энергетически более выгодные структуры. Возможно, это удастся сделать быстрым ударным нагружением графита. Обходной путь — синтез методом химического осаждения углерода из пара на подложку с подходящим кристаллическим строением, когда энергетически будет выгодно зарождение кристаллов именно нужной фазы. В принципе такой метод синтеза вполне отработан для получения крупных алмазов. Тот факт, что структуры фаз известны, облегчает поиск подходящей подложки.
Хочется надеяться, что результаты наших расчетов в обозримом будущем приведут к попыткам синтеза этих новых модификаций. Если они окажутся успешными, и такие сверхплотные фазы углерода действительно будут получены, это может открыть новые технологические применения для углерода, структурное разнообразие которого, похоже, неисчерпаемо. Но уже сейчас, имея примеры натрия и других веществ, для которых теоретическое предсказание предшествовало экспериментальному синтезу, ясно, что математическое моделирование структур способно не только объяснять полученные экспериментаторами результаты, но и предсказывать новые направления исследований.
Литература: 1. Oganov A.R., Glass C.W.Crystal structure prediction using ab initio evolutionary techniques: Principles and applications. Journal of Chemical Physics, 2006, т. 124, с. 244704. 2. Oganov A.R., Chen J., Gatti C., Ma Y.-Z., Ma Y.-M., Glass C.W., Liu Z., Yu T., Kurakevych O.O., Solozhenko V.L.Ionic high-pressure form of elemental boron. Nature, 2009, 457, с. 863. 3. Ma Y., Eremets M.I., Oganov A.R., Xie Y., Trojan I., Medvedev S., Lyakhov A.O., Valle M., Prakapenka V.Transparent dense sodium. Nature, 2009, т. 458, с. 182. 4. Zhu Q., Oganov A.R., Salvado M., Pertierra P., Lyakhov A.O.Denser than diamond: Ab initio search for superdense carbon allotropes. Physical Review B, 2011, т. 83, с. 193410.
Источник: elementy.ru.
Рейтинг публикации:
|