КАРТОЧКА ПРОЕКТА ФУНДАМЕНТАЛЬНЫХ И ПОИСКОВЫХ НАУЧНЫХ ИССЛЕДОВАНИЙ,
ПОДДЕРЖАННОГО РОССИЙСКИМ НАУЧНЫМ ФОНДОМ

Информация подготовлена на основании данных из Информационно-аналитической системы РНФ, содержательная часть представлена в авторской редакции. Все права принадлежат авторам, использование или перепечатка материалов допустима только с предварительного согласия авторов.

 

ОБЩИЕ СВЕДЕНИЯ


Номер 19-72-00082

НазваниеРазвитие теоретических методов для описания роста и эволюции некосселевских кристаллов и тонких плёнок

РуководительРедьков Алексей Викторович, Кандидат физико-математических наук

Организация финансирования, регион федеральное государственное бюджетное учреждение науки Институт проблем машиноведения Российской академии наук, г Санкт-Петербург

Период выполнения при поддержке РНФ 07.2019 - 06.2021 

Конкурс№40 - Конкурс 2019 года «Проведение инициативных исследований молодыми учеными» Президентской программы исследовательских проектов, реализуемых ведущими учеными, в том числе молодыми учеными.

Область знания, основной код классификатора 02 - Физика и науки о космосе, 02-201 - Теория конденсированного состояния

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

Код ГРНТИ29.19.03


СтатусУспешно завершен


 

ИНФОРМАЦИЯ ИЗ ЗАЯВКИ


Аннотация
Проект направлен на решение фундаментальной проблемы изучения механизмов роста и эволюции так называемых некосселевcких кристаллов - сложных кристаллических соединений, ячейка которых содержит несколько различных типов атомов в неэквивалентных позициях. Именно к таким соединениям относится значительная часть материалов, используемых в современной технологии (в том числе полупроводниковой). Примером могут послужить GaAs, GaN, SiC, InAs и многие другие. Однако, основная доля теоретических исследований, существующих к настоящему моменту, посвящена простым косселевским кристаллам, состоящих из атомов одного сорта. Это значительно упрощает рассмотрение задач, хотя эксперименты зачастую ставятся на многокомпонентных кристаллах (LiF, NaCl, GaN,SiC и т.д.). При этом экспериментальные результаты часто хорошо описываются построенными для однокомпонентных кристаллов теориями при выборе соответствующих «эффективных» параметров: коэффициента диффузии, атомарного объема, энергий активации зарождения и миграции адатомов и т.д. Несмотря на этот факт, очевидно, что в реальном многокомпонентном кристалле у адатома каждого компонента эти параметры различаются. Как пример, можно привести процесс роста кристалла нитрида галлия (GaN). Энергии связи адатома галлия с поверхностью GaN равна ~0.15 eV, тогда как энергия адатома азота ~0.05 eV, но в зависимости от условий роста, соотношения потоков и давлений компонентов в системе, «эффективная» энергия связи, при которой косселевы ростовые модели дают хороший описательный результат, может быть совсем другой. Это связано с различным вкладом тех или иных адатомов в процесс роста: если адатомов одного типа в системе мало, то именно они лимитируют процесс роста, а от свойств других адатомов рост практически не зависит. Следовательно, даже зная свойства отдельных компонентов, на данный момент мы пока не можем в общем случае аналитически предсказать, каков будет этот "эффективный" параметр - для каждых ростовых условий эта величина будет своя. Таким образом, крайне актуальна задача расширения существующих для однокомпонентных систем теорий на многокомпонентный случай. Новизна и сложность этой задачи состоит в том, чтобы правильно учесть влияние каждого компонента на конечные «эффективные» параметры системы и получить выражения, связывающие свойства отдельных компонентов, ростовые условия и скорость роста кристалла как целого. Помимо этого, необходимо проанализировать и явления, которые могут быть вызваны самим фактом наличия нескольких компонентов и повлиять на рост кристалла. Следует отметить, что работа по расширению различных ростовых и иных моделей на многокомпонентные случаи ведется очень активно и востребована, так как в промышленности сейчас используется большое количество многокомпонентных соединений, например, вышеупомянутый GaN, на основе которого делают светодиоды и высокочастотные транзисторы. Именно поэтому появляется все больше работ, посвященных переносу тех или иных теорий на двух- и более компонентные системы. При этом, однако, для многих теорий и моделей до сих пор нет многокомпонентных «аналогов». В проекте будет рассмотрен ряд таких теорий и решен цикл матфизических задач о росте и эволюции некосселевского кристалла в рамках модели терраса-ступень-излом по различным механизмам, будут найдены выражения для скорости перемещения отдельного излома, ступени, группы параллельных ступеней, а также скорость движения группы ступеней, образованных винтовой дислокацией. После решения этих задач будут получены выражения, связывающие нормальную скорость роста кристалла, концентрации отдельных компонентов в растворе или газовой фазе, их свойства (длину свободного пробега по поверхности, энергию связи, атомарные объемы), упругие напряжения, температуру, а также общее пересыщение в системе. Будет теоретически рассмотрен ряд явлений, которые, как представляется автору проекта, могут иметь место вследствие наличия нескольких компонентов, и оценено их влияние на скорость роста кристалла. В частности, планируется рассмотреть явление выпадения отдельных компонентнов в виде островков на террасах между ступенями, их движение в градиенте концентрации около ступени, а также самоорганизацию подобных островков в квазипериодические слои. Также будет рассмотрена и задача о перемещении дислокации за счет анизотропии распределения упругих напряжений вокруг неё, поскольку может повлиять на процесс роста кристалла. Для подтверждения полученных аналитических результатов, будет проведено численное моделирование рассмотренных явлений в некосселевском кристалле методом молекулярной динамики с использованием потенциалов, описывающих реальные многокомпонентные кристаллы. Большинство численных результатов будет сравнено с результатами, полученными аналитически.

Ожидаемые результаты
В рамках выполнения проекта ожидаются следующие результаты: 1)Решена задача о росте идеального некосселевского кристалла по механизму терраса-ступень-излом при наличии упругих напряжений как из паров собственных компонентов, так и при поступлении компонентов на поверхность за счет химической реакции. Найдены аналитические выражения для скорости движения отдельного излома на ступени (в случае, когда изломов на ступени мало), скорости движения отдельной ступени, скорости движения группы эквидистантных ступеней. Получены выражения для скорости роста кристалла как за счет движения ряда параллельных ступеней, так и по спиральному механизму. При выводе уравнений учтена разница в свойствах (коэффициентах диффузии по поверхности, атомарных объемах) различных компонентов, составляющих некосселевский кристалл. 2) Решена задача о росте идеального некосселевского N-компонентного кристалла по механизму Чернова из паров или раствора собственных компонентов. Найдены аналитические выражения для скорости движения отдельной ступени, скорость движения группы эквидистантных ступеней. Получено выражение для скорости роста кристалла как за счет движения группы параллельных ступеней, так и по спиральному механизму. 3) Рассмотрен ряд явлений, имеющих место при росте некосселевского кристалла ввиду наличия нескольких компонентов. Решена задача о выпадении новой фазы отдельных компонентов в пространстве между ступенями , найден критерий, когда эта фаза выпадает, а когда нет, и изучено влияние этого явления на скорость роста кристалла. Рассмотрена задача о взаимодействии ступени с зародышами отдельных компонентов на поверхности некосселевского кристалла. Найдена скорость растворения и "перемещения" островка в градиенте концентрации, создаваемой ступенью, и оценено влияние этого явления на скорость роста кристалла. Рассмотрено явление самоорганизации островков из отдельного компонента в квазипериодические слои вдали от ступени. Найдены условия возникновения подобной самоорганизации и показано, что движение ступени, и как следствие, скорость роста кристалла в этом ростовом режиме будет носить периодический во времени характер, что ранее, насколько известно автору, никогда не обсуждалось в литературе. 4) Решена задача о движении дислокации в некосселевском кристалле под действием упругих напряжений. Найдено выражение для распределения вакансий в подрешетках различных компонентов вокруг дислокации и их потоки. Вычислена скорость перемещения ядра дислокации. 5) Проведено численное моделирование описанных процессов методом молекулярной динамики с использованием потенциалов для некосселевских кристаллов. Численно найдены зависимости скорости движения ступеней от различных параметров (концентраций различных компонентов в системе, температуры, упругих напряжений, общего пересыщения в системе). Большинство полученных с помощью молекулярной динамики результатов сопоставлены с аналитическими результатами из предыдущих пунктов для верификации аналитических результатов. Полученные результаты, как надеется автор, позволяют по-новому взглянуть на процессы роста некосселевских кристалов и будут полезны с фундаментальной точки зрения, поскольку расширят имеющийся теоретический задел в области роста кристаллов на многокомпонентные системы, и демонстрируют новые явления, возникающие при росте некосселевских кристаллов вследствие наличия нескольких компонентов, которые, насколько известно автору, никогда не рассматривались ранее теоретически. Эти явления могут оказать существенное влияние на скорость роста кристалла и представляется важным понимание, когда и при каких условиях они возникают и как ими управлять. Помимо этого результаты будут полезны и с прикладной точки зрения, поскольку полученные выражения, связывающие скорость роста кристалла и ростовые условия, а также критерии возникновения упомянутых явлений могут быть использованы для описания реального процесса роста обширного класса различных многокомпонентных кристаллических соединений, используемых в промышленности, и как следствие, позволят оптимизировать ростовые условия для таких соединений как с точки зрения скорости роста, так и с точки зрения режима роста, поскольку наиболее низкодефектные кристаллы и тонкие плёнки чаще всего растут при низких пересыщениях именно в тех режимах, которые рассмотрены в настоящем проекте. Поскольку рост кристаллов и тонких плёнок в настоящее время является многомиллиардной индустрией, оптимизация процессов роста, на взгляд автора, имеет большую общественную значимость.


 

ОТЧЁТНЫЕ МАТЕРИАЛЫ


Аннотация результатов, полученных в 2019 году
В рамках работы над проектом классические подходы для описания роста кристаллов по механизму терраса-ступень-излом развиты на многокомпонентный случай. Рассмотрен механизмы Бартона-Кабреры-Франка, когда рост происходит из паров собственных компонентов, и механизм Чернова, когда кристалл растет из раствора. Получены выражения, связывающие свойства отдельных компонентов, составляющих кристалл: их коэффициенты диффузии, атомарные объемы, концентрации, и скорость роста кристалла. Показано, что процесс роста многокомпонентных кристаллов подобен однокомпонентному случаю и может быть описан теми же формулами, что и в классических работах Бартона-Кабрера-Франка и Чернова. Отличием является то, что в многокомпонентном случае обобщенные коэффициенты диффузии и встраивания определяются свойствами отдельных компонентов кристалла и их поверхностной концентрацией, а средняя длина свободного пробега адатомов зависит от скорости обратной химической реакции. Определено соотношение, связывающее обобщенные коэффициенты и свойства отдельных компонентов. Поскольку при росте многокомпонентного кристалла за счет химической реакции можно изменять парциальные давления реагентов независимо друг от друга, появляются дополнительные рычаги управления процессом роста: пересыщение (химическое сродство) реакции, средняя длина свободного пробега адатомов по поверхности, а также равновесные поверхностные концентрации компонентов, которые определяют влияние каждого компонента на обобщенные коэффициенты, могут контролироваться независимо, путем подбора соответствующих соотношений парциальных давлений реагентов выбора абсолютного давления, тогда как в однокомпонентном случае роста из собственных паров, единственным рычагом управления при заданной температуре является пересыщение (давление пара). Рассмотрены некоторые эффекты, свойственные многокомпонентным системам с химическими реакциями, такие как зависимость скорости роста кристалла при заданном пересыщении от абсолютной скорости химической реакции, зародышеобразование островков отдельных компонентов на террасах между ступенями и др. Полученные результаты позволяют описать рост произвольного многокомпонентного кристалла по различным механизмам, и могут быть использованы для анализа и объяснения обширных существующих экспериментальных данных, выбора оптимальных условий роста, а также для дальнейшего развития теории роста кристаллов.

 

Публикации

1. Редьков А.В., Кукушкин С.А., Осипов А.В. Growth of a multicomponent crystal via Chernov's mechanism Journal of Physics: Conference Series (JPCS), 1410 (2019) 012039 (год публикации - 2019) https://doi.org/10.1088/1742-6596/1410/1/012039

2. Редьков А.В., Кукушкин СА. Development of Burton-Cabrera-Frank theory for the growth of a Non-Kossel crystal via chemical reaction Crystal Growth&Design, volume 20, p. 2590-2601 (год публикации - 2020) https://doi.org/10.1021/acs.cgd.9b01721

3. Редьков А.В., Кукушкин С.А. Теоретические аспекты роста многокомпонентных кристаллических тонкопленочных покрытий из газовой фазы Сборник трудов " XXVI Научно-технической конференции с участием зарубежных специалистов "Вакуумная наука и техника", стр. 36 (год публикации - 2019)

4. Редьков А.В., Кукушкин С.А., Гангрская Е.С. Peculiarities of growth of a Non-Kossel crystal via Chernov’s mechanism Book of abstracts of the "International conference "Mechanisms and Non-linear problems of nucleation and growth of crystals and thin films 2019", p. 65 (год публикации - 2019)


Аннотация результатов, полученных в 2020 году
В рамках работы над проектом в 2020-2021 году исследовались процессы и явления, происходящие при росте многокомпонентных кристаллов, широко используемых в настоящее время в индустрии, например при производстве светодиодов. Было обнаружено и исследовано явление формирования островков отдельных фаз на террасах между ступенями при росте многокомпонентного кристалла по механизму терраса-ступень-излом. Показано, что в отличие от классических, хорошо изученных однокомпонентных систем (например, кремния Si), где переход между ростом по механизму терраса-ступень-излом и нуклеацией островков на террасах жестко определяется потоками компонента к поверхности, в многокомпонентных системах появляется возможность независимо контролировать скорость передвижения многокомпонентных ступеней на поверхности (путём изменения многокомпонентного пересыщения) и нуклеацию и рост однокомпонентных островков (путем изменения потока отдельного компонента при сохранении многокомпонентного пересыщения). Обнаруженные режимы открывают новые возможности по управлению процессом роста многокомпонентных кристаллов и системой однокомпонентных островков на нем. Так, в процессе выполнения проекта продемонстрировано, что для управления системой островков чистого компонента можно использовать не только поток этого компонента, как это происходит в однокомпонентных системах, но и многокомпонентную ступень, которая является стоком для всех компонентов. Регулируя положение этой ступени относительно островка, можно вызывать его рост, растворение, а также перемещение. В процессе работы над проектом последовательно рассмотрены эти режимы, и показана связь между скоростью перемещения ступеней и наличием островков на поверхности. В ряде случаев подобные островки могут существенно затруднить рост кристалла, и вызвать возникновение различного рода ростовых неустойчивостей, в том числе неустойчивости Маллинза-Секерки. Это явление приводит к нарушению послойного роста кристалла, и может вызвать появление дополнительных дефектов в нем. Показано, что использование явлений и режимов, исследованных в настоящем проекте, позволяет формировать динамические (существующие только в процессе роста кристалла) упорядоченные массивы островков отдельных фаз, размеры и расстояния между которыми можно регулировать, меняя многокомпонентное пересыщение. Показано, что такие динамически-упорядоченные системы островков могут быть использованы в качестве шаблона для самоиндуцированного роста упорядоченных массивов нитевидных нанокристаллов. Также показана возможность с помощью комбинации различных ростовых режимов формировать трехмерные структуры включений чистых компонентов в объеме кристалла in situ, без использования сложных литографических методов. Помимо этого, в рамках проекта исследовано перемещение дислокаций в многокомпонентном кристалле под воздействием упругих напряжений и найдены выражения, связывающие скорость перемещения дислокаций и свойства отдельных компонентов. Решение этой задачи позволило углубить теоретическое понимание процессов, происходящих при перемещении дислокаций в многокомпонентных кристаллах под воздействием упругих напряжений. Также в рамках проекта было проведено численное моделирование процесса роста многокомпонентного кристалла методом молекулярной динамики в программном пакете LAMMPS. Основное внимание было уделено изучению роста актуального в настоящий момент нитриду галлия GaN. При моделировании была изучена кинетика адатомов на поверхности, найдены коэффициенты диффузии и времена жизни отдельных типов адатомов на поверхности GaN различных ориентаций с использованием различных cуществующих потенциалов в широком диапазоне температур. Как оказалось, свойства индивидуальных адатомов (время жизни, коэффициенты диффузии) отличаются от экспериментально получаемых значений, хотя численно-определенные энергии активации процессов диффузии и испарения достаточно близки. Проведен цикл численных экспериментов по процессу осаждения адатомов Ga и N на ступенчатую поверхность кристалла при различных температурах осаждения, потоках адатомов, и их соотношениях. Численно изучен как послойный, так и спиральный рост кристалла GaN, и обнаружены некоторые новые механизмы диффузии, которые могут наблюдаться и в реальных системах, связанные с формированием особых структур из несколких атомов. Численно обнаружено явление формирования зерен различных политипов в процессе роста при больших пересыщениях. Проведено сопоставление численных результатов с полученными ранее теоретическими выкладками. Анимация некоторых из режимов роста GaN по результатам вычислений методом МД доступна по ссылке https://youtu.be/tmQh6ey5R-g Полученные в отчетном году результаты могут быть использованы при анализе и объяснении различных экспериментальных данных по росту многокомпонентных кристаллов и для разработки новых bottom-up методов роста самоорганизованных наноструктур в многокомпонентных системах без применения дорогостоящих методов литографии. Результаты имеют как фундаментальное значение, поскольку позволяют по-новому взглянуть на процессы роста многокомпонентных кристаллов, и глубже понять их, так и прикладное, поскольку могут быть применены к описанию широкого класса востребованных в настоящее время ростовых процессов и материалов.

 

Публикации

1. Редьков А.В., Кукушкин С.А. Nucleation of nano-islands of pure components during growth of a multicomponent crystal via step-flow mode Journal of Physics: Conference Series (JPCS), A V Redkov and S A Kukushkin 2020 J. Phys.: Conf. Ser. 1695 012003 (год публикации - 2020) https://doi.org/10.1088/1742-6596/1695/1/012003

2. Редьков А.В., Кукушкин С.А. Dynamic Interaction of Steps and Nanoislands during Growth of a Multicomponent Crystal Crystal Growth&Design, - (год публикации - 2021) https://doi.org/10.1021/acs.cgd.1c00349

3. Редьков А. В., Кукушкин С. А. Теоретические аспекты роста многокомпонентных кристаллических тонкопленочных покрытий из газовой фазы Наноиндустрия, Журнал "Наноиндустрия", том. 13, номер S2, 2020, стр. 173-176, Тезисы конференции "Вакуумная наука и техника 2019" (год публикации - 2020) https://doi.org/10.22184/1993-8578.2020.13.2s.173.176

4. Редьков А.В., Кукушкин С.А. Theory of multicomponent crystal growth due to chemical reaction Book of abstracts of the "18th International Conference on Thin Films & 18th Joint Vacuum Conference", p.74 (год публикации - 2020)

5. Редьков А.В., Кукушкин С.А. Self-organization of nanoislands of pure components during growth of a multicomponent crystal via step-flow mode Book of abstracts "7th International School and Conference on Optoelectronics, Photonics, Engineering and Nanostructures", p. 25 (год публикации - 2020)

6. Редьков А.В., Кукушкин С.А. Dynamical interaction and self-organization of nanoislands of pure components at the terraces during the growth of a multicomponent crystal Book of abstracts of the international conference "British Association for Crystal Growth Joint Spring meeting 2021", P. 114 (год публикации - 2021)


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