УДК 538.931
ВКЛАД ВНЕСЕННЫХ ТОЧЕЧНЫХ ДЕФЕКТОВ В САМОДИФФУЗИЮ ПО МАЛОУГЛОВЫМ ГРАНИЦАМ КРУЧЕНИЯ В ГЦК МЕТАЛЛАХ
1ПОЛЕТАЕВ Г. М., 2НОВОСЕЛОВА Д. В., 3КОВАЛЕНКО В. В., 1 ГУРОВА Н. М.,
4коханенко д. в.
1 Алтайский государственный технический университет им. И.И. Ползунова, 656038, Алтайский край, г. Барнаул, пр. Ленина, 46
2
Кузбасский институт Федеральной службы исполнения наказаний России, 654066, Кемеровская область, г. Новокузнецк, пр. Октябрьский, 49
3Сибирский государственный индустриальный университет, 654007, г. Новокузнецк, пр. Бардина, д. 25, ауд. 111ГТ
4Финансовый университет при Правительстве РФ, Барнаульский филиал, 656038, Алтайский край, г. Барнаул, пр. Ленина, 54
АННОТАЦИЯ. Методом молекулярной динамики исследуется взаимодействие точечных дефектов с малоугловыми границами кручения (100), (110), (111) в ГЦК металлах N1, Си, А1 и их роль в самодиффузии по рассматриваемым границам. Проведен расчет энергии образования и энергии связи точечных дефектов с малоугловыми границами кручения. Выяснено, что точечные дефекты преимущественно закрепляются в узлах сетки зернограничных винтовых дислокаций. Показано, что внесенные вакансии играют важную роль в диффузии по границам кручения. Вклад внесенных межузельных атомов в зернограничную диффузию в условиях термодинамического равновесия меньше на несколько порядков по сравнению с вкладом вакансий.
КЛЮЧЕВЫЕ СЛОВА: молекулярная динамика, граница кручения, винтовая дислокация, точечный дефект, вакансия, межузельный атом, диффузия.
ВВЕДЕНИЕ
Границы зерен в поликристаллах делят на два типа: границы наклона и границы кручения. В первом случае ось разориентации зерен, то есть ось, вокруг которой одно кристаллическое зерно повернуто относительно другого, лежит в плоскости границы. В случае границ кручения ось разориентации перпендикулярна плоскости границы. Известно, что малоугловые границы наклона представляют собой систему геометрически необходимых (зернограничных) краевых дислокаций. Малоугловые границы кручения также имеют упорядоченную структуру, представляющую собой сетку винтовых дислокаций [1 - 5] (рис. 1, а). В работе [6] нами была проведена идентификация винтовых дислокаций в границах кручения, ориентированных в плоскостях (100), (110), (111) в металлах N1, Си, А1 с ГЦК кристаллической решеткой. Было показано, что границы (100) содержат квадратную сетку винтовых дислокаций 1/2<110>; границы (110) - прямоугольную сетку винтовых дислокаций двух типов: 1/2<110> и 1<100>; границы (111) - гексагональную сетку винтовых дислокаций 1/4<112>. С увеличением угла разориентации зерен размеры ячеек дислокационной сетки уменьшаются. В центре ячеек сетки, т.е. вне ядер зернограничных дислокаций, структура соответствует практически идеальному кристаллу. В той же работе [6] приведены зависимости энергии рассматриваемых границ кручения от угла разориентации, полученные при использовании двух типов потенциалов межатомного взаимодействия: парного Морзе и многочастичного Клери-Розато.
В работе [7] методом молекулярной динамики нами было проведено исследование самодиффузии по границам кручения (100), (110), (111) в ГЦК металлах Ni, Cu, Al. Было показано, что самодиффузия по малоугловым границам кручения осуществляется посредством образования цепочек смещенных атомов вдоль ядер зернограничных винтовых дислокаций. При этом местами инициации цепочек смещений атомов выступают узлы дислокационной сетки. Цепочки атомных смещений зачастую имеют замкнутую форму (рис. 1, б). Образование замкнутых цепочек смещений атомов вдоль ядер винтовых зернограничных дислокаций объяснялось нами в [7] рекомбинацией пары избыточных дефектов - вакансии и межузельного атома, образующихся в соседних узлах дислокационной сетки при смещении атомов вдоль ядра дислокации от одного узла сетки к другому. Согласно [7], наибольшей диффузионной проницаемостью (минимальной энергией активации) обладают границы (110), наименьшей - границы (111). С ростом угла разориентации энергия активации падает.
а) б)
Рис. 1. Малоугловая граница кручения в плоскости (100): а) изображение, полученное в работе [5] с помощью электронного микроскопа;
б) замкнутая цепочка смещенных атомов в границе кручения (100) с углом разориентации 15° (изображение взято из [7], где оно было получено с помощью моделирования диффузии вдоль границы)
Результаты, приведенные выше, относятся к структурно «чистым» границам зерен, не содержащим каких-либо «внесенных» дефектов. На самом деле в реальных поликристаллах межзеренные границы являются эффективными стоками различных дефектов: дислокаций, вакансий, межузельных атомов и их комплексов. Почти все границы содержат те или иные внесенные дефекты, которые в различной степени оказывают влияние на зернограничные процессы: на диффузию, пластическую деформацию и т.д. Концентрация и тип этих дефектов в основном зависят от условий создания структуры (деформация, быстрое охлаждение или нагрев, радиационное воздействие и т.д.). Практически всегда это термодинамически неравновесные состояния. Поэтому предсказывать концентрацию и тип дефектов с помощью законов равновесной термодинамики во многих случаях не имеет смысла. Однако, с другой стороны, по известным значениям энергии образования или связи дефекта с границей можно сделать вывод о соотношении концентраций различных дефектов в границе, вероятности их появления и закрепления в ней.
Настоящая работа посвящена исследованию взаимодействия точечных дефектов с малоугловыми границами кручения (100), (110), (111) в ГЦК металлах N1, Си, А1 и их роли в самодиффузии по рассматриваемым границам кручения.
ОПИСАНИЕ МОДЕЛИ
Перед проведением основных молекулярно-динамических исследований предшествовал этап создания границы кручения в компьютерной модели. Граница создавалась в середине кристаллического расчетного блока в результате поворота одной половины блока относительно другой на заданный угол 0 вокруг оси, перпендикулярной границе (рис. 2). В работе рассматривалось три варианта ориентации плоскости границы в ГЦК кристалле: (100), (110), (111). После создания расчетного блока (на рис. 2 он выделен пунктирным контуром) следовал этап релаксации структуры, в течение которой структура вблизи границы зерен приобретала состояние с наименьшей потенциальной энергией в данных условиях. В течение структурной релаксации на границе происходило формирование геометрически необходимой сетки зернограничных винтовых дислокаций. После релаксации расчетный блок охлаждался до 0 К.
На границы расчетного блока были наложены жесткие граничные условия, то есть атомы вблизи границ блока в процессе моделирования оставались неподвижными. Это было необходимо для удержания границы зерен с заданными параметрами в середине расчетного блока в течение всего моделирования. Расчетные блоки содержали примерно 30000 атомов.
Взаимодействия атомов металла друг с другом в молекулярно-динамической модели описывались с помощью многочастичного потенциала Клери-Розато [8]. Данный потенциал хорошо зарекомендовал себя в ряде расчетов структурно-энергетических характеристик металлов, выполненных методом молекулярной динамики [9, 10]. Шаг интегрирования по времени в методе молекулярной динамики был равен 5 фс. Этого шага достаточно для моделирования тепловых движений атомов металлов вплоть до температуры плавления [11, 12]. Температура в модели задавалась через начальные скорости атомов согласно распределению Максвелла.
е
кристаллическое зерно 2
Рис. 2. Создание границы кручения в модели. 0 - угол разориентации зерен. Конечный расчетный блок выделен пунктирным контуром
При исследовании самодиффузии вдоль границ кручения проводились молекулярно-динамические эксперименты продолжительностью 0,3 - 0,5 нс. Температура расчетного блока в течение этих экспериментов удерживалась постоянной. Коэффициент диффузии рассчитывался через среднеквадратические смещения атомов по формуле (для коэффициента вдоль оси Х)
1 N
N Е (х0/ - X)2
Дх = ^-, (1)
х 2Г
здесь х0/ - координата начального положения /-го атома; Х1 - координата /-го атома в момент времени N - количество атомов в расчетном блоке.
Перед расчетом коэффициента диффузии расчетный блок охлаждался до 0 К для исключения тепловых (недиффузионных) смещений атомов. Для коэффициентов диффузии вдоль осей У и Ъ расчет производился аналогично. Средний коэффициент самодиффузии рассчитывался как среднее арифметическое коэффициентов Дх, Ду и Д При исследовании самодиффузии вдоль границ зерен полагалось, что границы имеют диффузионную толщину 5 А. Такая толщина была выбрана не случайно: как известно, в кристаллах, содержащих границы зерен, основные диффузионные процессы протекают вдоль границ в малом слое толщиной порядка 5 А [13].
Основные диффузионные характеристики (энергия активации диффузии и предэкспоненциальный множитель в уравнении Аррениуса) определялись по зависимостям 1пДг от Т1, где Дг - коэффициент самодиффузии по границе зерен, Т - температура, по методике описанной в работе [10].
РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ
При исследовании взаимодействия точечных дефектов с рассматриваемыми границами кручения в различные позиции «структурно правильных» (т.е. не содержащих избыточных дефектов) границ вводилась одна вакансия или межузельный атом. Сравнивались три различные позиции точечного дефекта в границе: в узле сетки зернограничных дислокаций, в ядре дислокации и в центре ячейки дислокационной сетки (с почти идеальной кристаллической структурой). После введения точечного дефекта в ту или иную позицию проводилась релаксация структуры с последующим охлаждением.
Было замечено, что в большинстве случаев точечные дефекты в процессе структурной релаксации мигрировали в ближайший узел дислокационной сетки. Это происходило даже без дополнительной термоактивации. На рис. 3 изображены примеры смещений атомов в результате подобных миграций точечных дефектов при изначальном введении их в область ядра дислокации в границе 0 = 6° (100). На рисунках видно, что атомы смещаются преимущественно вдоль ядер дислокаций. При этом зачастую возникали параллельные смещения атомов в противоположных направлениях, характерные при движении винтовой дислокации. Небольшое перемещение ядер дислокаций, сходящихся в одном узле, при поглощении узлом точечного дефекта, свидетельствует в пользу смещения самого узла дислокационной сетки на относительно небольшое расстояние.
Вышеописанные результаты имели место при всех углах разориентации вплоть до значений, соответствующих большугловым границам (когда невозможно выделить отдельные зернограничные дислокации). Вместе с тем, для большеугловых границ кручения также наблюдались периодически расположенные вдоль границы места, которые стремились занять введенные точечные дефекты.
место \ ч
введения г* V
вакансии
конечное положение вакансии
а)
б)
Рис. 3. Атомные смещения в малоугловой границе кручения 0 = 6° (100) в N1 в результате миграции вакансии (а) и межузельного атома (б) при изначальном их введении в область ядра винтовой зернограничной дислокации
Энергия связи точечного дефекта с границей зерен - это работа, которую нужно затратить, чтобы удалить данный дефект из границы в объем кристалла. Энергия связи Еь рассчитывалась по формуле
ЕЬ = иЬг0 - иЬг , (2)
где иЬуо - потенциальная энергия расчетного блока, содержащего границу зерен и точечный дефект, находящийся в объеме зерна (в кристалле) на достаточном расстоянии от границы;
- потенциальная энергия расчетного блока, содержащего границу зерен и точечный дефект, находящийся в границе в узле дислокационной сетки. При получении величин иьу0 и иЬу? проводилась обязательная релаксация структуры и последующее охлаждение до 0 К (это необходимо во избежание влияния тепловых смещений атомов).
Энергия образования точечного дефекта в границе зерен Е/ - это работа создания этого дефекта в границе. Она может быть рассчитана как разность энергии образования дефекта в чистом кристалле Е/0 и энергии связи дефекта с границей зерен Еь:
Е/ = Е/о - Еь , (3)
Значения энергии образования вакансии в чистом кристалле Е/о были взяты из справочников [14 - 16]: 1,64 эВ для N1, 1,22 эВ для Си и 0,68 эВ для А1. Энергии образования межузельных атомов (в позиции гантели <100>) были взяты из работы [10], где они были получены с помощью метода молекулярной динамики с использованием того же потенциала, что и в настоящей работе: 4,98 эВ для N1, 3,58 эВ для Си и 2,25 эВ для А1.
В табл. 1 приведены полученные в настоящей работе значения энергии образования точечных дефектов в границах кручения (100), (110), (111) для трех разных улов разориентации: 6°, 15°, 24°.
Для границ кручения, ориентированных в плоскости (110), энергии образования точечных дефектов наиболее низкие. Для границ (111), наоборот, энергии образования близки к значениям для чистого кристалла. Это коррелирует со значениями энергии самих границ зерен, полученных нами ранее в работе [6]: границы (111) имеют наименьшую энергию, границы (110) - наибольшую среди рассмотренных ориентаций. Другими словами,
чем выше энергия границы (а, значит, и ее «дефектность»), тем выше энергия связи точечного дефекта с границей и, следовательно, его энергия образования в границе. Угол разориентации зерен в пределах значений, соответствующих малоугловым границам, практически не влиял на значения энергии образования дефектов в границе.
Таблица 1
Энергия образования вакансии и межузельного атома в границах кручения (100), (110), (111)
с углами разориентации 6°, 15°, 24° (эВ)
(100) (110) (111)
6° 15° 24° 6° 15° 24° 6° 15° 24°
Ni вакансия 0,57 0,47 0,74 0,03 0,10 0,23 1,35 1,58 1,60
межд.атом 2,09 1,93 2,22 0,51 0,77 0,62 4,09 4,12 4,30
Cu вакансия 0,23 0,34 0,41 0,17 0,11 0,06 1,16 1,09 0,95
межд.атом 1,74 1,53 1,02 0,11 0,21 0,20 2,76 2,52 2,37
A1 вакансия 0,12 0,08 0,17 -0,06 0,07 0,11 0,57 0,59 0,63
межд.атом 0,91 0,87 0,93 0,19 0,23 0,16 2,01 1,92 2,12
Полученные результаты свидетельствуют в пользу того, что границы кручения, ориентированные в плоскости (111), в реальных металлах должны содержать наименьшую концентрацию избыточных дефектов и являться сравнительно плохими стоками для внешних точечных дефектов. Более эффективными стоками являются границы (110). С ростом угла разориентации выше плотность узлов дислокационной сетки. В связи с этим, очевидно, чем выше угол разориентации зерен, тем более эффективным стоком для точечных дефектов является граница.
При исследовании вклада введенных точечных дефектов в самодиффузию вдоль границ кручения в границу вводилась одна вакансия или один межузельный атом. При различных температурах определялся коэффициент зернограничной самодиффузии Дг, после чего по зависимостям 1пДг от T1 рассчитывались энергия активации и предэкспоненциальный множитель в уравнении Аррениуса. В связи с громоздкостью расчетов (для одного металла необходимо было получить 6 температурных зависимостей коэффициента диффузии - по два (с участием вакансии и межузельного атома) для трех различных углов разориентации 6°, 15°, 24°) вычисления были проведены только для Ni. Результаты приведены в табл. 2.
Таблица 2
Энергия активации самодиффузии в N1 по «структурно чистым» и содержащим один избыточный точечный дефект границам кручения (100), (110) и (111) с углами разориентации 6°, 15° и 24° (эВ)
(100) (110) (111)
6° 15° 24° 6° 15° 24° 6° 15° 24°
чистый 1,49 1,22 1,01 0,38 0,32 0,30 - - -
вакансия 0,83 0,62 0,66 0,22 0,25 0,34 0,90 1,02 1,29
межд.атом 0,78 0,55 0,61 0,44 0,23 0,31 0,20 0,20 0,16
Как видно из таблицы, энергия активации самодиффузии при введении точечного дефекта уменьшается. Кроме того, введение вакансии или межузельного атома почти одинаково влияет на повышение диффузионной проницаемости границ (100) и (110), тогда как в случае границ (111) введение межузельного атома приводит к более значительному снижению энергии активации диффузии по сравнению с введением вакансии. Это лишний раз свидетельствует о высокой степени упорядоченности структуры границы кручения (111).
Количественная оценка вклада внесенных точечных дефектов в самодиффузию вдоль рассматриваемых границ кручения проводилась следующим образом. Вероятности того, что некоторый участок границы будет содержать одну вакансию или межузельный атом, соответственно равны
п ехр
Г / кт
р = п ехр
Г /
кт
(4)
где Е/у и Е/, - энергии образования вакансии и межузельного атома (табл. 1); к - постоянная Больцмана; Т - температура; п - количество узлов дислокационной сетки (где преимущественно находятся точечные дефекты) в рассматриваемом участке границы. Тогда вероятность того, что участок границы не будет содержать точечные дефекты:
Р0 = 1 - Ру - Рг .
(5)
Коэффициент диффузии вдоль всей границы зерен равен сумме вкладов диффузии по участкам без дефектов А0, с вакансией и с межузельным атомом Д-:
Ах = А + А + А.
(6)
Эти вклады можно найти как произведения коэффициентов диффузии на рассматриваемых участках на соответствующие вероятности наличия этих участков в границе:
А
:р0 В0ехр[- кт},
А,
■Рув0у ехр(- <ттт
Б, = рАы ехр^- ТТ^, (7)
где ТУ и к - энергии активации самодиффузии в рассматриваемом участке границы, содержащей одну вакансию или межузельный атом (табл. 2); и Б, -
предэкспоненциальные множители, полученные по зависимостям 1иАг от Т"1 при введении в расчетный блок одной вакансии или межузельного атома.
Относительный вклад в самодиффузию вакансии и межузельного атома:
V = А .100%
и
V = Аг- -100%.
(8)
На рис. 4 изображены зависимости вклада внесенных вакансий в самодиффузию по границам кручения (100) с углами разориентации 6°, 15°, 24° от температуры. Как видно из рисунка, вклад вакансий довольно существенный. Более того, для границ (111) и (110) вклад вакансий оказался близок к 100 %. То есть можно сказать, что диффузия вдоль таких границ протекает преимущественно за счет внесенных вакансий. Межузельные атомы, за счет высокой энергии образования, вносят, в целом, незначительный вклад по сравнению с вакансиями. Их вклад в большинстве случаев близок к нулю, только для границ (110) он иногда достигает порядка 1 %.
у
Следует иметь в виду, что проведенные расчеты выполнены для «идеальных» условий: для «структурно правильных» границ и равновесной концентрации дефектов при условии отсутствия других внесенных дефектов. Известно, что границы зерен в реальных поликристаллах зачастую находятся в неравновесном состоянии и включают ряд внесенных (избыточных) точечных дефектов и дислокаций. Тем не менее, проведенные расчеты позволяют качественно судить о роли вакансий и межузельных атомов в диффузии по малоугловым границам кручения.
О 200 400 600 800 1000 1200 1400 1600
Т. К
Рис. 4. Зависимости вклада внесенных вакансий в самодиффузию от температуры по границам кручения (100) с углами разориентации 6°, 15°, 24°
ЗАКЛЮЧЕНИЕ
Согласно расчетам, проведенным в настоящей работе с помощью метода молекулярной динамики, внесенные вакансии (под внесенными понимаются дефекты, необязательные при структурообразовании границ зерен) играют важную роль в диффузии по границам зерен кручения. Вклад внесенных межузельных атомов в зернограничную диффузию меньше на несколько порядков по сравнению с вкладом вакансий. Точечные дефекты в границах кручения закрепляются в узлах дислокационной сетки. Диффузия с их участием протекает путем миграции дефектов от одного узла сетки к другому с образованием цепочки смещенных атомов.
При расчете энергии образования точечных дефектов в границах кручения было выяснено, что чем выше энергия образования границы (выше ее «дефектность»), тем ниже энергия образования точечных дефектов в ней. Энергия образования точечных дефектов оказалась наиболее низкой для границ (110), для границ (111) энергии образования вакансии и межузельного атома близки к энергиям образования в чистом кристалле. От угла разориентации энергия образования точечных дефектов практически не зависела. Однако их концентрация в границах с ростом угла разориентации, очевидно, увеличивается, так как растет плотность узлов сетки зернограничных дислокаций.
Судя по полученным данным, границы (111) должны содержать наименьшую концентрацию точечных дефектов и являться сравнительно плохими стоками для внешних дефектов. При этом самих границ кручения (111) должно быть в поликристалле больше других. Наиболее эффективными стоками для дефектов являются, согласно результатам настоящей работы, границы, имеющие наибольшую энергию образования. Из рассмотренных границ кручения - это границы (110).
СПИСОК ЛИТЕРАТУРЫ
1. De Hosson J. Th. M., Vitek V. Atomic structure of (111) twist grain boundaries in f.c.c. metals // Philosophical Magazine A, 1990, vol. 61, no. 2, pp. 305-327.
2. Ghafoor A., Ahmad S. A., Faridi B. A. S. The structure of (001) CSL twist boundaries in fcc metals // Turkish Journal of Physics, 1998, vol. 22, no. 8, pp. 789-795. URL: http://journals.tubitak.gov.tr/physics/issues/fiz-98-22-8/fiz-22-8-3-96116.pdf (дата обращения 03.02.2017).
3. Ghafoor A., Faridi B. A. S., Ahmad A. Multiple structures of (110) CSL twist boundaries in fcc metals // Turkish Journal of Physics, 2001, vol. 25, no. 1, pp. 35-42. URL: http://journals.tubitak.gov.tr/physics/issues/fiz-01-25-1fiz-25-1-5-98070.pdf (дата обращения 03.02.2017).
4. Shallcross S., Sharma S., Pankratov O. A. Twist boundary in graphene: energetics and electric field effect // Journal of Physics: Condensed Matter, 2008, vol. 20, no. 45, pp. 454224.
5. Belov A. Yu., Scholz R., Scheerschmidt K. Dissociation of screw dislocations in (001) low-angle twist boundaries: a source of the 30° partial dislocations in silicon // Philosophical Magazine Letters, 1999, vol. 79, no. 8, pp. 531-538.
6. Полетаев Г. М., Мартынов А. Н., Старостенков М. Д. Структура и энергия границ зерен кручения в ГЦК металлах // Фундаментальные проблемы современного материаловедения. 2010. Т. 7, № 4. С. 27-34.
7. Мартынов А. Н., Полетаев Г. М., Старостенков М. Д. Атомный механизм диффузии по малоугловым границам кручения в ГЦК металлах // Письма о материалах. 2011. Т. 1, № 1. С. 43-46.
8. Cleri F., Rosato V. Tight-binding potentials for transition metals and alloys // Physical Review B, 1993, vol. 48, no. 1, pp. 22-33.
9. Полетаев Г. М., Кулабухова Н. А., Старостенков М. Д. Потенциалы межатомного взаимодействия в системах Pd-H и Ni-H // Химическая физика и мезоскопия. 2011. Т. 13, № 3. C. 411-418.
10. Полетаев Г. М., Старостенков М. Д. Вклады различных механизмов самодиффузии в ГЦК-металлах в условиях равновесия // Физика твердого тела. 2010. Т. 52, № 6. С. 1075-1082.
11. Старостенков М. Д., Медведев Н. Н., Полетаев Г. М. К вопросу о систематических погрешностях в ММД // Измерения, автоматизация и моделирование в промышленности и научных исследованиях; межвузовский сборник / под. ред. Г.В. Леонова. Бийск: Изд-во АлтГТУ, 2005. С. 5-8.
12. Старостенков М. Д., Медведев Н. Н., Полетаев Г. М., Терещенко О. А. Гамильтониан замкнутой системы, моделируемой с помощью ММД // Фундаментальные проблемы современного материаловедения. 2006. № 2. С. 46-48.
13. Лариков Л. Н., Исайчев В. И. Диффузия в металлах и сплавах. Киев: Наукова думка, 1987. 511 с.
14. Волленбергер Г. Й., Хирт Дж. П., Виртман Дж., Хаазен П. Точечные дефекты // В кн.: Физическое металловедение. Том 3. Физико-механические свойства металлов и сплавов / пер. с англ., под ред. Р. Кана, П. Хаазена. М.: Металлургия, 1987. C. 5-74.
15. Штремель М. А. Прочность сплавов. Часть 1. Дефекты решетки. М.: Металлургия, 1982. 280 с.
16. Орлов А. Н., Трушин Ю. В. Энергии точечных дефектов в металлах. М.: Энергоатомиздат, 1983. 80 с.
CONTRIBUTION OF INTRODUCED POINT DEFECTS TO SELF-DIFFUSION ALONG THE LOW-ANGLE TWIST BOUNDARIES IN FCC METALS
1Poletaev G. M., 2Novoselova D. V., 3Kovalenko V. V., :Gurova N. M., 4Kokhanenko D. V. :Altai State Technical University, Barnaul, Russia
2Kuzbass Institute of the Federal Penitentiary Service of Russia, Novokuznetsk, Russia
3Siberian State Industrial University, Novokuznetsk, Russia
4Financial University under the Government of the Russian Federation, Barnaul brunch, Barnaul, Russia
SUMMARY. The interaction of point defects with low-angle (100), (110), (111) twist boundaries in fcc metals Ni, Cu, Al and their role in self-diffusion along considered boundaries is studied by the molecular dynamics method. The energy of formation and binding energy of point defects with low-angle twist boundaries is calculated. It was found that point defects are mainly fixed in the grid nodes of grain-boundary screw dislocations. It is shown that the introduced vacancies play an important role in diffusion along the boundaries of twist grain boundaries. The contribution of introduced interstitial atoms to grain-boundary diffusion under thermodynamic equilibrium conditions is less by several orders of magnitude in comparison with the contribution of vacancies.
KEYWORDS: molecular dynamics, twist boundary, screw dislocation, point defects, vacancy, interstitial atom, diffusion.
REFERENCES
1. De Hosson J. Th. M., Vitek V. Atomic structure of (111) twist grain boundaries in f.c.c. metals. Philosophical Magazine A, 1990, vol. 61, no. 2, pp. 305-327. http://dx.doi.org/10.1080/01418619008234943
2. Ghafoor A., Ahmad S. A., Faridi B. A. S. The structure of (001) CSL twist boundaries in fcc metals. Turkish Journal of Physics, 1998, vol. 22, no. 8, pp. 789-795. URL: http://journals.tubitak.gov.tr/physics/issues/fiz-98-22-8/fiz-22-8-3-96116.pdf (accessed February 03, 2017).
3. Ghafoor A., Faridi B. A. S., Ahmad A. Multiple structures of (110) CSL twist boundaries in fcc metals. Turkish Journal of Physics, 2001, vol. 25, no. 1, pp. 35-42. URL: http://journals.tubitak.gov.tr/physics/issues/fiz-01-25-1fiz-25-1-5-98070.pdf (accessed February 03, 2017).
4. Shallcross S., Sharma S., Pankratov O. A. Twist boundary in graphene: energetics and electric field effect. Journal of Physics: Condensed Matter, 2008, vol. 20, no. 45, pp. 454224. https://doi.org/10.1088/0953-8984/20/45/454224
5. Belov A. Yu., Scholz R., Scheerschmidt K. Dissociation of screw dislocations in (001) low-angle twist boundaries: a source of the 30° partial dislocations in silicon. Philosophical Magazine Letters, 1999, vol. 79, no. 8, pp. 531-538. http://dx.doi.org/10.1080/095008399176896
6. Poletaev G. M., Martynov A. N., Starostenkov M. D. Struktura i energiya granits zeren krucheniya v GTsK metallakh [Structure and energy of boundaries of torsion grains in fcc metals]. Fundamental'nye problemy sovremennogo materialovedeniya [Basic Problems of Material Science (BPMS)], 2010, vol. 7, no. 4, pp. 27-34.
7. Martynov A. N., Poletaev G. M., Starostenkov M. D. Atomnyy mekhanizm diffuzii po malouglovym granitsam krucheniya v GTsK metallakh [Atomic diffusion mechanism for small-angle torsion boundaries in fcc metals]. Pis'ma o materialakh [Letters on Materials], 2011, vol. 1, no. 1, pp. 43-46.
8. Cleri F., Rosato V. Tight-binding potentials for transition metals and alloys. Physical Review B, 1993, vol. 48, no. 1, pp. 22-33. https://doi.org/10.1103/PhysRevB.48.22
9. Poletaev G. M., Kulabukhova N. A., Starostenkov M. D. Potentsialy mezhatomnogo vzaimodeystviya v sistemakh Pd-H i Ni-H [Potentials of interatomic interaction in systems Pd-H and Ni-H]. Khimicheskaya fizika i mezoskopiya [Chemical Physics and Mesoscopy], 2011, vol. 13, no. 3, pp. 411-418.
10. Poletaev G. M., Starostenkov M. D. Contributions of different mechanisms of self-diffusion in face-centered cubic metals under equilibrium conditions. Physics of the Solid State, 2010, vol. 52, no. 6, pp. 1146-1154. https://doi.org/10.1134/S1063783410060065
11. Starostenkov M. D., Medvedev N. N., Poletaev G. M. K voprosu o sistematicheskikh pogreshnostyakh v MMD [To the question of systematic errors in MMD]. Izmereniya, avtomatizatsiya i modelirovanie vpromyshlennosti i nauchnykh issledovaniyakh; mezhvuzovskiy sbornik [Measurements, automation and modeling in industry and scientific research; intercollegiate collection]. Biysk: AltGTU Publ., 2005, pp. 5-8.
12. Starostenkov M. D., Medvedev N. N., Poletaev G. M., Tereshchenko O. A. Gamil'tonian zamknutoy sistemy, modeliruemoy s pomoshch'yu MMD [The Hamiltonian of a closed system modeled by MMD ]. Fundamental'nye problemy sovremennogo materialovedeniya [Basic Problems of Material Science» (BPMS)], 2006, vol. 3, no. 2, pp. 46-48.
13. Larikov L. N., Isaychev V. I. Diffuziya v metallakh i splavakh [Diffusion in metals and alloys]. Kiev: Naukova dumka Publ., 1987. 511 p.
14. Vollenberger G. Y., Khirt Dzh. P., Virtman Dzh., Khaazen P. Tochechnye defekty [Point defects]. Fizicheskoe metallovedenie. Tom 3. Fiziko-mekhanicheskie svoystva metallov i splavov [Physical Metallurgy. Vol. 3. Physical and mechanical properties of metals and alloys]. Moscow: Metallurgiya Publ., 1987, pp. 5-74.
15. Shtremel' M. A. Prochnost' splavov. Chast' 1. Defekty reshetki [Strength of alloys. Part 1. Defects of the lattice]. Moscow: Metallurgiya Publ., 1982. 280 p.
16. Orlov A. N., Trushin Yu. V. Energii tochechnykh defektov v metallakh [The energy of point defects in metals]. Moscow: Energoatomizdat Publ., 1983. 80 p.
Полетаев Геннадий Михайлович, доктор физико-математических наук, профессор, заведующий кафедрой высшей математики и математического моделирования, АлтГТУ, тел. (3852) 290865, е-mail: gmpoletaev@mail. т
Новоселова Дарья Викторовна, кандидат физико-математических наук, преподаватель кафедры гуманитарных, социально-экономических и естественнонаучных дисциплин, КИ ФСИН России, е-mail: dmitdarya@mail. т
Коваленко Виктор Викторович, доктор физико-математических наук, доцент, профессор кафедры физики, СибГИУ, е-mail: [email protected]
Гурова Наталья Михайловна, кандидат физико-математических наук, доцент кафедры физики, АлтГТУ, е-mail: gurova. nmg@yandex. т
Коханенко Дмитрий Васильевич, кандидат физико-математических наук, доцент, доцент кафедры учета и информационных технологий в бизнесе, Финуниверситет, Барнаульский филиал, е-mail: DVKohanenko@fa.т