Full Text

Введение Теплопроводность можно определить как количество тепла, передаваемого через единицу толщины материала - в направлении, нормальном к поверхности единицы площади, - из-за единичного температурного градиента в условиях устойчивого состояния. Теплопроводность - необходимая характеристика для рассеивания образовавшейся тепловой энергии в системе. Прогресс нанотехнологий за последние два десятилетия предоставил разнообразные системы с высокой теплопроводностью из различных типов материалов и топологических форм. Определение такой характеристики, как теплопроводность материала, является задачей, которой посвящено множество теоретических и экспериментальных исследований. Ввиду этого использование методов математического и компьютерного моделирования является перспективным для описания моделей процессов теплообмена. В данной работе приведены результаты расчетов теплофизических характеристик наносистем. В представленной статье в продолжение предыдущих исследований авторов [1; 2] рассмотрена задача расчета коэффициента теплопроводности однородных наносистем. Отметим, что решению задачи определения коэффициента теплопроводности кремния на наноуровне посвящен обширный ряд работ. В исследовании [3] коэффициент теплопроводности равен 20 Вт/(м·К), а в работе [4] - 235 Вт/(м·К). В статье [5] при моделировании использовался потенциал Стиллингера - Вебера, коэффициент теплопроводности равен 20 Вт/(м·К). В исследовании [6] коэффициент теплопроводности получился 37-43 Вт/(м·К). Отметим, что в источнике [7] получена диаграмма значений коэффициента теплопроводности для систем различной размерности. 1. Постановка задачи Моделирование системы проводилось с помощью программного пакета LAMMPS [8] методом молекулярной динамики [9]. Выбор потенциала взаимодействия играет важную роль при моделировании. В настоящее время для кремния доступно огромное количество многочастичных потенциалов. Наибольшее распространения при моделировании металлических и полупроводниковых систем применяют EAM (метод погруженного атома) потенциалы [10; 11] и MEAM (модифицированный метод погруженного атома) потенциалы [12-15]. В данной работе использовался потенциал взаимодействия MEAM. Задача расчета коэффициента теплопроводности может быть решена несколькими способами. Их описание можно найти в работах [16; 17]. В данной работе для определения теплофизических свойств наносистем использован метод Грина - Кубо. Решение задачи определения коэффициента теплопроводности осуществляется в несколько этапов. На первом этапе происходит формирование необходимой наноструктуры, которая помещается в расчетную ячейку с периодическими граничными условиями по осям x, y, z. Второй этап предполагает достижение равновесного состояния моделируемой системы посредством применения NPT ансамбля. Графики температуры и давления на протяжении данного этапа представлены на рис. 1, 2. Целью данного этапа является релаксация системы, при которой расчетная ячейка моделирования доводится до равновесного объема при заданных температуре, давлении и числе атомов. Скорости атомов масштабируются таким образом, чтобы при моделировании температура держалась на заданном уровне. Данный этап должен быть достаточно продолжительным, чтобы система уравновесилась. После достижения равновесия, постоянная решетки не будет существенно изменяться и колебаться. В молекулярной динамике важную роль играет выбор временного шага интегрирования, так как он определяет точность и эффективность расчета. Если временной шаг мал, то требуется больше времени вычислений для достижения желаемого времени моделирования. Если шаг велик, то большие силы могут развиваться на более легких атомах, что приведет к искусственно быстрым атомам. В работе используется шаг по времени 0,2 фс. Тестовые расчеты показали, что такая величина шага является подходящей для решения данного класса задач. Третий этап предполагает достижение равновесного состояния моделируемой системы посредством применения NVE ансамбля. На этом этапе моделирования объем (V) и полная энергия (E) системы «заморожены» из-за применения ансамбля NVE. «Удаление» термостатирования (T) и баростатирования (P) освобождает атомные движения из-за искусственного изменения масштаба. На рис. 3 представлен график температуры на протяжении всего третьего этапа моделирования. Видно, что система держит температуру в районе 300 К без применения термостатов. На рис. 4, 5 представлены графики общей энергии и объема расчетной ячейки на протяжении всего третьего этапа моделирование NVE. Рис. 1. График температуры системы 4´4´4 на втором этапе моделирования при температуре 300 К Fig. 1. Temperature graph of a 4´4´4 system at the second stage of modeling at a temperature of 300 K Рис. 2. График давления системы 4´4´4 на втором этапе моделирования при температуре 300 К Fig. 2. Pressure graph of a 4´4´4 system at the second stage of modeling at a temperature of 300 K Рис. 3. График температуры системы 4´4´4 на третьем этапе моделирования при температуре 300 К Fig. 3. Temperature graph of a 4´4´4 system at the third stage of modeling at a temperature of 300 K Рис. 4. График энергии системы 4´4´4 на третьем этапе моделирования при температуре 300 К Fig. 4. Energy diagram of a 4´4´4 system at the third stage of modeling at a temperature of 300 K а б Рис. 5. График изменения объема расчетной ячейки системы 4´4´4 на третьем этапе моделирования при температуре: а - 300 К; б - 1000 К Fig. 5. The graph of the change in the volume of the computational cell of the 4´4´4 system at the third stage of modeling at a temperature of: а - 300 K; б - 1000 K 2. Расчет коэффициента теплопроводности В молекулярно-динамических расчетах величину коэффициента теплопроводности можно вычислить различными способами [8]. В данной работе используется формализм Грина - Кубо (Green - Kubo), который связывает автокорреляционную функцию теплового потока с коэффициентом теплопроводности. Тепловой поток может быть рассчитан из колебаний потенциальной и кинетической энергии атома и тензора напряжения атома в стационарном уравновешенном моделировании. Это отличие от двух предыдущих неравновесных методов молекулярной динамики (NEMD), где энергия течет непрерывно между горячими и холодными областями в моделируемом образце. Коэффициент теплопроводности в модели Грина - Кубо рассчитывается по следующей формуле [18; 19]: (1) где - коэффициент теплопроводности - мерной системы с линейным размером, - температура, - постоянная Больцмана, - компонента потока тепла. Постоянная времени - это минимальное время, необходимое для того, чтобы автокорреляционная функция теплового потока затухала до нуля. Этот метод более подходит для анизотропных систем. Основным недостатком этого метода является то, что для затухания до нуля автокорреляционной функции теплового потока требуется много времени, а получаемые значения теплопроводности зависят от размера системы. Автокорреляционные функции справа в формуле (1) оцениваются в равновесии, без градиента температуры. Автокорреляция - статистическая взаимосвязь между случайными величинами из одного ряда, но взятыми со сдвигом, например, для случайного процесса - со сдвигом по времени. Автокорреляционная функция может быть определена как: Общий тепловой поток в системе вычисляется как где - плотность теплового потока. Порядок пределов в формуле (1) имеет большое значение. При правильных порядках пределов можно вычислить корреляционные функции с произвольными граничными условиями и применять формулу (1). Существуют различные формы записи уравнения (1) [20-25]. Существуют ситуации, когда формула (1) неприменима. Во-первых, для маленьких систем, которые изучаются в мезоскопической физике, термодинамический предел не имеет смысла. Во-вторых, во многих низкоразмерных системах перенос тепла аномальный и теплопроводность существенно отклоняется от экспериментальных значений [26; 27]. В таких случаях невозможно взять пределы в формуле (3). Поэтому тепловую проводимость рассматривают как функцию от длины L. В литературе, посвященной этой тематике [26; 27], обычно в формуле (1) меняют верхний предел интегрирования на L. Другой способ использования формулы Green - Kubo для конечных систем заключается в том, чтобы внедрять бесконечные резервуары, как это сделано в работах [28; 29]. В данной работе для расчета теплопроводности кремниевых наноструктур был выбран метод Green - Kubo. Для вычисления коэффициента теплопроводности использовалась следующая запись формулы (1) [8; 30]: , где - объем системы. На рис. 6 и 7 приведены значения коэффициента теплопроводности чистого кремния при различной температуре в сравнение с экспериментальными данными. Полученные в нашей работе значения хорошо согласуются с таковыми, представленными в работе для масштабных образцов [32]. Известно, что в металлах и полупроводниках работают два механизма переноса тепла: за счёт электронов и посредством фононов. В чистых металлах процесс переноса тепла осуществляется в основном посредством свободных электронов, концентрация которых на единицу объема в металлах весьма велика. При температурах, близких к абсолютному нулю (), концентрация фононов будет мала. Большую роль будет играть электронная теплопроводность. Следует отметить, что в чистых металлах при нормальных условиях теплопроводность электронного газа много больше решеточной теплопроводности. На рис. 6, 7 видно, что при повышении температуры, когда фононная проводимость будет играть большую роль, чем электронная, результаты моделирования и экспериментальные данные хорошо согласуются. Графики, показанные на рис. 8, 9, представляют собой усредненные корреляционные кривые тепловых потоков. Время сбора автокорреляционных данных должно быть достаточно большим, чтобы автокорреляция затухала до нуля. Были построены графики теплопроводности в различные моменты времени при температуре 300 К (рис. 9) и 1000 К (рис. 10). Как видно из рис. 10-11, после 20 пс система стабилизируется, и коэффициент теплопроводности выходит на стационарное значение. Рис. 6. График теплопроводности, полученный при моделировании системы 4´4´4, в сравнение с экспериментальными данными Fig. 6. Thermal conductivity graph obtained when simulating a 4´4´4 system in comparison with experimental data Рис. 7. График теплопроводности, полученный при моделировании системы 10´10´10, в сравнении с экспериментальными данными Fig. 7. Thermal conductivity graph obtained when simulating a 10´10´10 system in comparison with experimental data Рис. 8. График автокорреляционной функции теплового потока температуре 300 К Fig. 8. Graph of the autocorrelation function of the heat flux at a temperature of 300 K Рис. 9. График автокорреляционной функции теплового потока температуре 1000 К Fig. 9. Graph of the autocorrelation function of the heat flux at a temperature of 1000 K Рис. 10. График теплопроводности при температуре 300 К для системы 4´4´4 Fig. 10. Thermal conductivity graph at a temperature of 300 K for a 4´4´4 system Рис. 11. График теплопроводности при температуре 1000 К для системы 4´4´4 Fig. 11. Thermal conductivity graph at a temperature of 1000 K for a 4´4´4 system Заключение Построены физическая и математическая модели процессов теплопроводности в наноматериалах на основе кремния высокой плотности с использованием потенциала взаимодействия MEAM. Изложена методика расчета характеристик данного процесса. Также в работе показаны особенности и различия потенциалов взаимодействия, используемых в молекулярно-динамических расчетах. Определены зависимости коэффициента теплопроводности кремния от температуры для кристаллов различных размеров. Результаты расчетов показали, что значение коэффициента теплопроводности уменьшается с ростом температуры, что подтверждается экспериментальными исследованиями, результаты которых приведены в работах [34; 35].

About the authors

A. V Severyukhin

UdmFRC UB RAS

O. Yu Severyukhina

UdmFRC UB RAS

A. V Vakhrushev

UdmFRC UB RAS

References

Statistics

Views

Abstract - 349

PDF (Russian) - 296