CC BY
59
ФИЗИКА КОНДЕНСИРОВАННОГО СОСТОЯНИЯ
УДК 534.1
РАСЧЕТ ТЕПЛОЕМКОСТИ И СРЕДНЕКВАДРАТИЧНЫХ СМЕЩЕНИЙ ТВЕРДОГО РАСТВОРА NiAl ПО ФОНОННЫМ СПЕКТРАМ
Е.Ю. Краюшкина1, В.Е. Холодовский2
брянский государственный технический университет, Брянск, Российская Федерация
e-mail: senichenkova@mail.ru
2Брянский государственный университет, Брянск, Российская Федерация e-mail: v-e-holodovsky@yandex.ru
С использованием ранее разработанной методики построена динамическая модель, учитывающая взаимодействия между атомами на первых двух координационных сферах, для кристалла NiAl, имеющего структуру типа CsCl. Исходя из экспериментальных данных по нейтронному рассеянию, выполнены расчеты дисперсионных кривых и упругих констант. Построены кривые плотности фо-нонных состояний, а также рассчитан вклад каждой подрешетки в энергию колебаний кристалла. По этим данным определена зависимость среднеквадратичных смещений атомов каждой подрешетки и теплоемкости кристалла от температуры.
Ключевые слова: кристаллическая решетка, динамическая модель, межатомное взаимодействие, дисперсионные кривые, упругие константы, фононный спектр, среднеквадратичное смещение, теплоемкость.
THE HEAT CAPACITY AND MEAN-SQUARE DISPLACEMENTS OF THE SOLID SOLUTION NiAl OVER THE PHONON SPECTRUM
E.Yu. Krayushkina1, V.E. Holodovsky2
1 Bryansk State Technical University, Bryansk, Russian Federation e-mail: senichenkova@mail.ru
2 Bryansk State University, Bryansk, Russian Federation e-mail: v-e-holodovsky@yandex.ru
The article presents the dynamic model based on the previously developed technique which considers an interaction between the atoms in the first two coordination spheres for a crystal NiAl with the structure of CsCl. It discusses calculations of the dispersion curves and elastic constants made with the help of experimental data of neutron scattering. The curves of the phonon density of states are built. Each sublattice contribution to the crystal oscillation energy was estimated. The relations between the mean-square atomic displacements ofeach sub lattice as well as the heat capacity of crystals and the temperature are found.
Keywords: crystal lattice, dynamic model, interatomic interaction, dispersion curves, elastic constants, phonon spectrum, mean-square displacement, specific heat.
Введение. Упорядоченный твердый раствор NiAl широко используется в аэрокосмической промышленности, что связано с такими его уникальными свойствами, как высокая температура плавления
(1911 К), низкая плотность (5,85 г/см3), а также стойкость к окислению [1]. Исследованию динамики решетки в-№А1 посвящены работы [2-4]. Экспериментальное исследование МА1 проведено в работе [2], методом неупругого рассеяния нейтронов были измерены кривые дисперсии фононов по основным направлениям симметрии, для описания результатов эксперимента использована модель Борна-Кармана с учетом трех и четырех ближайших соседей. Исходя из первых принципов с использованием метода линейного отклика и ультрамягких псевдопотенциалов в работе [3] были исследованы фононный спектр и фо-нонная плотность состояний интерметаллидов М3А1 и МА1. С учетом первых принципов в работе [4] также исследованы аномалии акустических ветвей по основным направлениям симметрии.
В настоящей статье кривые дисперсии фононов, полная и парциальная плотность фононных состояний были определены исходя из динамической модели, разработанной в работах [5-9], в которой для определения сил межатомного взаимодействия использовался механизм вандерваальсовых связей. По этим данным построены зависимости среднеквадратичных смещений атомов каждой подрешетки и теплоемкости кристалла от температуры.
Методика расчета. Кратко опишем основные положения используемой методики расчета (более подробно с методикой можно ознакомиться в работе [10]).
Рассмотрим бинарную кубическую кристаллическую решетку типа СвС1. Предположим, что в процессе колебаний атомов решетки в каждом из них находится внутриатомный диполь, одним из полюсов которого является остов атома, а другим — центр заряда его внешней электронной оболочки (ВЭО). Пусть Л — множество всех числовых наборов £ = (к,к,1), в которых существуют числа 1 < Н,к,1 < N, задающих положение узла кристаллической решетки. Обозначим через плечо дипольного момента атома А%, наведенного в результате его относительного перемещения с соседними атомами из первой и второй координационных сфер.
На внутриатомный диполь атома действует кулоновская сила Q£ со стороны диполей остальных атомов решетки, в результате чего центр заряда ВЭО атома и плечо его дипольного момента получают некоторое приращение Ар?: Р^ = + Ар?.
Наведенный дополнительный дипольный момент создает частичное экранирование силы Q£, тогда внешняя сила, действующая на остов атома, становится равной: ^ = Qí — (в/а) Ар%, где а — поляризуемость атома ; в = Я2/(4п£0), д — эффективный заряд диполя атома.
На остов атома А% также действует сила реакции на излучение его внутриатомного диполя. В первом приближении эта сила пропорциональна плечу диполя [11]: = — (2в/(3г3)) Р%, где г — радиус сферы, плотность потока энергии излучения внутриатомного диполя через которую считается равным работе силы реакции за единицу времени.
Пусть и% — вектор смещения остова атома А% из положения равновесия в момент времени ¿. Суммируя все силы, действующие на остов атома А%, получаем уравнение его движения: ^и% = Е%— (в/а) р%.
В адиабатическом приближении на любом временном промежутке энергия, поглощаемая атомом, совпадает с энергией, излучаемой им. Это условие будет выполнено, если принять, что внешняя, частично экранированная кулоновская сила уравновешивается силой реакции: Е% + = 0. Тогда в состоянии термодинамического равновесия уравнение движения остова атома А% принимает вид:
Пусть Aç — атом первой или второй подрешетки, Aç/ — атом, находящийся на первой или второй координационной сфере атома Aç, eçç/ — единичный вектор, указывающий направление от узла Dç к узлу Dç/, wçç/ = uç/ — uç — вектор относительного перемещения атомов Aç и Aç/, rçç/ =< wçç/, eçç/ > eçç/, tçç/ = wçç/ — rçç/ — радиальная и тангенциальная составляющие этого вектора. За плечо дипольного момента Pç, наведенного в результате перемещения атома Am относительно соседних атомов, примем линейную комбинацию векторов rçç/ и tçç/. Тогда уравнение (1) движения атома Aç можно записать в виде
^uç = ((air — Cit) < wçç/, eçç/ > eçç/ + auwçç/) +
Здесь a1r, ац, a2r, a2í — некоторые константы, определяемые свойствами кристалла; S¡(£) — множество всех мультииндексов G Л, нумерующих атомы, находящиеся на 1-й координационной сфере атома A%.
Разделяя атомы кристалла по подрешеткам и используя верхние индексы, приходим к двум системам уравнений вида (2), описывающих колебания каждой подрешетки.
Решение полученных систем уравнений ищется в виде уравнений бегущих волн: u1(t) = v1 sin (Kr^ — wt), u|(t) = v2 sin (Kr| — wt) , где K = kxex + kyey + kzez — волновой вектор; v1, v2 — векторы поляризации; w — частота.
(1)
Î'eSi(ç)
+ ((^2r - < , > + a2tw^ ). (2)
í'eS2(í)
В результате приходим к системе
(aX - + bx2 + bzУ2 + byZ2 = 0
(a^ - ^iw2)yi + bzX2 + by2 + bxZ2 = 0 (a1 - ^iw2)zi + byX2 + ЬхУ2 + bz2 = 0 (aX - ^2^2)x2 + bxi + bzyi + byzi = 0
(ay - ^2^2)y2 + bz xi + by i + bx zi = 0
(aX - ^2^2)z2 + byxi + bxyi + bzi = 0,
где xi, yi, zi и x2, y2, z2 — координаты векторов поляризации. С точностью до круговой перестановки индексов коэффициенты системы (3) находятся по формулам
al = 8a1 + 4air sin2 (kXa/2) + 4a! (sin2 a/2) + sin2 (kza/2)) ;
aX = 8a1 + 4aXr sin2 (kXa/2) + 4a2t (sin2 (kya/2) + sin2 (kza/2)) ; bx = 8a2 cos (kXa/2) sin (kya/2) sin (kza/2); b = —8a1 cos(kXa/2)cos(kya/2)cos(kza/2); air + 2ctií air — au
ai = —3—; a2 = ,
где a — параметр решетки.
Для каждого разрешенного значения волнового вектора K система (3) определяет частоты wk¿, j = 1, 2, 3,..., 6, трех акустических и трех оптических колебательных мод. Кроме того, для заданных значений K, j и вычисленной частоты wk¿ она позволяет однозначно определить векторы поляризации vi и v2 так, чтобы выполнялось условие xi + У2 + z2 + x2 + y2 + = 1.
Результаты расчетов и обсуждения. В случае, когда волновой вектор K имеет одно из основных кристаллографических направлений система (3) приводит к дисперсионному соотношению:
2 A^i + Ai^2
ш =
1 ±А/1 - (AiA2 - B2)
V ^i + Ai^2 ( i 2 )
2^1
где А1, А2, В — некоторые величины, выражаемые через коэффициенты системы (3).
Результаты расчетов дисперсионных кривых по приведенным выше формулам, и экспериментальные значения, полученные методом неупругого рассеяния нейтронов и представленные в работе [2] для кристалла твердого раствора МА1, показаны на рис. 1. Наблюдается хорошая согласованность расчетных и экспериментальных данных. Следует отметить, что оптические ветви в большей мере определяются колебаниями более легких атомов алюминия, а акустические — колебаниями атомов никеля (М^/Ма ~ 2). Это также подтверждается кривыми парциальной плотности фононного спектра (рис. 2).
Рис. 1. Дисперсионные кривые (1), полученные по приведенным выше формулам, и экспериментальные значения (2), рассчитанные методом неупругого рассеяния нейтронов и взятые из работы [2] для кристалла твердого раствора NiAl при температуре 296 K
Рис. 2. Кривые полной (а) и парциальной (б, в) плотности фононных состояний для сплава NiAl
В континуальном приближении, когда выполняется условие К а ^ 1 справедливо неравенство А1А2 — В2 ^ 1. Поэтому можно полагать, что для акустических мод справедливо соотношение
А1А2- В2
w2 =
Ä2 Ш + Ai Д2
(4)
Сравнивая соотношение (4) с дисперсионными соотношениями из теории упругости для основных кристаллографических направлений, приходим к равенствам
2^1 + + а^ = аб^;
4^2 = а (С12 + С44); (5)
а2г + а2г — а24 — а2 = а (б11 — б44) .
С помощью подобранных значений силовых констант по формулам (5) были вычислены значения упругих констант. При расчетах использовалось экспериментально определенное значение параметра
о
решетки МА1 2,887 А [1]. Результаты, полученные в настоящей работе, и результаты расчетов, выполненных другими авторами, представлены в таблице. Следует отметить, что в силу некоторых причин значения упругих констант существенно расходятся со значениями, определенными разными авторами, проводившими подобные исследования.
Расчетные и экспериментальные значения упругих констант NiAl при комнатной температуре, а также значения этих констант, полученные
другими авторами
Упругая константа Расчетные значения, полученные в настоящей работе Экспериментальные значения [12] Значения, взятые из работ
[2] [13] [14] [15]
Cu, (GPa) 177 199 197 233 159,26 202,30
C12, (GPa) 98 137 119 121 144,14 146,10
C44, (GPa) 83 116 110 114 93,96 100,30
Полная энергия колебательной моды частотой ш при температуре
Т определяется по формуле Планка Е = - Нш сШ ( ). Согласно
2 \2^вТ /
классическим представлениям, вклад в эту энергию каждой подре-шетки пропорционален произведению |2, и следовательно, про-
изведению , где gKj = х2 + У2 + уКо = X + У2 + Таким
образом, для заданных значений К и ] вклад каждой подрешетки в полную энергию соответствующей моды описывается формулой
т^т ^mgKljj (hwK,j
= 0(.. ^ J2 ,, ч ^K,j Ct^ ,J
2(у1 ^ + уК о ^2) \2кв Г
Для расчета вклада каждой подрешетки в полную энергию кристалла удобно использовать понятие парциальной плотности фононного спектра (ш):
Dm (ш) =
1
»к gK,j
6MДш KK »1 gK j + »2gK j
^Дш (шК,к — Ш) ,
(6)
где М — число задаваемых значений вектора К, равномерно распределенных в зоне Бриллюэна; Аи — условная единица длины фононного спектра;
' 1, если х е [—Аи/2, Аи/2]; 0, если х е [—Аи/2, Аи/2].
(x) =
Формула (6) позволяет рассчитать отдельные значения парциальной плотности фононного спектра в точках, равномерно распределенных по всему фононному спектру с шагом Аи.
Подобно классической формуле, выражающей температурную зависимость энергии кристаллической решетки, вклад каждой подре-шетки в энергию одного моля выражается интегралом
Ек =
3Na 2
штах
J Dm (Ш)НШ ct^2kT) dU'
Аналогично, подобно классическим формулам, выражается вклад каждой подрешетки в теплоемкость кристалла, среднеквадратичные смещения атомов каждой подрешетки и другие характеристики, определяемые фононным спектром.
Кривые полной и парциальной плотности фононных состояний для №А1 приведены на рис. 2, а зависимость среднеквадратичного смещения атомов от температуры — на рис. 3. В настоящей работе также была рассчитана зависимость молярной теплоемкости су от температуры (рис. 4). Наблюдается хорошая согласованность полученных значений
Рис. 3. Зависисмость среднеквадратичного смещения атомов № (1) и Al (2) в сплаве NiAl от температуры
0 100 200 300 Г, К
Рис. 4. Расчетная (1) и экспериментальная (2) зависимость молярной теплоемкости №А1 от температуры (экспериментальные данные взяты из работы [16])
теплоемкости с экспериментальными данными во всем исследуемом диапазоне значений температуры.
Заключение. Построенная динамическая модель позволила исходя из экспериментальных данных провести расчеты дисперсионных кривых, определить значения упругих констант для атомов интерметал-лида NiAl, полную и парциальные плотности фононного состояний, а также зависимости среднеквадратичных смещений атомов каждой подрешетки и теплоемкости от температуры. Установленные результаты хорошо согласуются с данными, вычисленными другими авторами, и с экспериментальными значениями. Представленная модель дает возможность выполнить расчеты динамических и испускательных характеристик внутриатомных диполей рассматриваемого кристалла.
ЛИТЕРАТУРА
1. Miracle D.B. The physical and mechanical properties of NiAl // Acta metall, mater. 1993. Vol. 41. No. 3. P. 649-684.
2. Bulk and surface vibrational modes in NiAl / M. Mostolle, R.M. Nicklow, D.M. Zehner, S.-C. Lui, J.M. Mundenar, E.W. Plummer // Phys. Rev. B. 1989. Vol. 40. No. 5. P. 2856.
3. Фононные спектры Li2Ni3Al и B2NiAl: расчеты из первых принципов / Э.И. Исаев, А.И. Лихтенштейн, Ю.Х. Векилов, Е.А. Смирнова // Физика твердого тела. 2004. Т. 46. Вып. 7. С. 1158-1164.
4. Xiangyang Huang, Ivan I. Naumov, Karin M. Rabe. Phonon anomalies and elastic constants of cubic NiAl from first principles // Physical Review B. 2004. Vol. 70. P. 064301.
5. Холодовский В.Е., Мачихина И.О. Принцип длинных волн и дисперсионные соотношения для кубических кристаллических решеток в модели диполь-дипольных взамодействий // Известия СамНЦ РАН. Сер. Физика и электроника. 2009. № 11 (5). С. 49.
6. ХолодовскийВ.Е., Мачихина И.О., КульченковЕ.А. Дисперсионные соотношения для кубических кристаллических решеток в модели диполь-дипольных взаимодействий // Вестник ЮУрГУ Сер. Математика, физика, химия. 2009. № 12 (10). С. 92.
7. Холодовский В.Е., Мачихина И.О. Принцип длинных волн и фононные спектры кубических кристаллических решеток // Вестник ЮУрГУ. Сер. Математика. Механика. Физика. 2009. № 1 (22). С. 109.
8. Холодовский В.Е., Мачихина И.О., Кульченков Е.А. Расчет теплоемкости и среднеквадратичных смещений по фононным спектрам для кристаллов с ОЦК- и ГЦК-решеткой // Вестник ЮУрГУ Сер. Математика. Механика. Физика. 2010. № 2 (9). С. 101.
9. Мачихина И.О., Холодовский В.Е. Расчет теплофизических свойств кристаллов инертных газов по упругим константам // Известия СамНЦ РАН. 2012. № 14 (4). С. 116.
10. Холодовский В.Е., Краюшкина Е.Ю. Дисперсионные соотношения для кристаллов твердых растворов со структурой CsCl // Ученые записки Забайкальского государственного университета. Сер. Физика, математика, техника, технология. 2013. № 3. С. 121.
11. Холодовский В.Е., Сидоров А.А. Поток энергии и сила реакции на излучение внутриатомного диполя // Ученые записки Забайкальского государственного университета. Сер. Физика, математика, техника, технология. 2012. № 3. С. 141.
12. Rusovic N., Warlimont H.The Elastic Behaviour of $2-NiAl Alloys // Phys. Status Solidi A. 1977. Vol. 44. P. 609.
13. Ponomareva A.V., Isaev Eyvas, Vekilov Yu. Kh., Abrikosov I. Site preference and effect of alloying on elastic properties of ternary B2 NiAl-based alloys // Physical Review B. 2012. Vol. 85. P. 144117.
14. First-principles calculation on mechanical and thermal properties of B-NiAl with point defects (in Chinese) / Y.L. Lu, H.X. Hou, Z. Chen, et al. // Sci Sin-Phys Mech Astron. 2013. Vol. 43. P. 152.
15. Fu Hongzhi, Li XiaoFeng, Liu WenFang, Mab Yanming, Gao Tao, Hong Xinhua. Electronic and dynamical properties of NiAl studied from first principles. Intermetallics, 19, 1959 (2011).
16. Desai P.D. Thermodynamic Properties of Selected Binary Aluminum Alloy Systems // J. Phys. Chem. Ref. Data. 1987. Vol. 16. No. 1.
REFERENCES
[1] Miracle D.B. The physical and mechanical properties of NiAl. Acta metall mater., 1993, vol. 41, no. 3, pp. 649-684.
[2] Mostolle M., Nicklow R.M., Zehner D.M., Lui S.-C., Mundenar J.M., Plummer E.W. Bulk and surface vibrational modes in NiAl. Phys. Rev. B, 1989, vol. 40, iss. 5, p. 2856.
[3] Isayev E.I., Liechtenstein A.I., Vekilov Yu.H., Smirnova E.A. Phonon spectra of L12 Ni3Al and B2 NiAl: calculations from first principles. Solid State Physics. 2004, vol. 46, no. 7, pp. 1158-1164.
[4] Huang X., Naumov I.I., Rabe K.M. Phonon anomalies and elastic constants of cubic NiAl from first principles. Physical Review B, 2004, vol. 70, pp. 064301.
[5] Holodovsky V.E., Machikhina I.O. The principle of long waves and dispersion relations for the cubic crystal lattices in the model of dipole-dipole interactions. Izvestiya SamNC RAN. Ser. Physics and Electronics, 2009, vol. 11, no. 5, p. 49 (in Russ.).
[6] Holodovsky V.E., Machihina I.O., Kulchenkov E.A. Dispersion relations for the cubic crystal lattices in the model dipole-dipole interactions. Vestnik Yuzhno-Ural'skogo Gos. Universiteta. Ser. "Matematika. Fizika. Khimiya" [Herald YuUrGU. Ser. Mathematics, Physics, Chemistry], 2009, vol. 12, iss. 10, p. 92 (in Russ.).
[7] Holodovsky V.E., Machihina I.O. Principle of long waves and phonon spectra of cubic crystalline lattices. Vestnik Yuzhno-Ural'skogo Gos. Universiteta. Ser. "Matematika. Mekhanika. Fizika" [Herald YuUrGU. Ser. Mathematics. Mechanics. Physics,], 2009, vol. 12, iss. 10, p. 92 (in Russ.).
[8] Holodovsky V.E., Machihina I.O., Kulchenkov E.A. The heat capacity and mean-square displacements on the phonon spectra of crystals with bcc and fcc lattice. Vestnik Yuzhno-Ural'skogo Gos. Universiteta. Ser. "Matematika. Mekhanika. Fizika" [Herald YuUrGU. Ser. Mathematics. Mechanics. Physics], 2010, vol. 2, iss. 9, p. 101 (in Russ.).
[9] Machihina I.O., Holodovsky V.E. Calculation of thermophysical properties crystals of inert gases on the elastic constants. Izvestiya SamNC RAN, 2012, vol. 14, no. 4, p. 116 (in Russ.).
[10] Holodovsky V.E., Krayushkina E.Yu. Dispersion relations for crystals of solid solutions with the structure of CsCl. Uchenye zapiski Zabaykal 'skogo gos. universiteta [Scientific notes of the Zabaikalsky State University. Ser. Physics, mathematics, engineering, technology], 2013, no. 3, p. 121 (in Russ.).
[11] Holodovsky V.E., Sidorov A.A. The energy flux and the reaction force on the intra-dipole radiation. Uchenye zapiski Zabaykal'skogo gos. universiteta [Scientific notes of the Zabaikalsky State University. Ser. Physics, mathematics, engineering, technology], 2012, no. 3, p. 141 (in Russ.).
[12] Rusovic N., Warlimont H. The Elastic Behaviour of -NiAl Alloys. Phys. Status SolidiA, 1977, vol. 44, p. 609.
[13] Ponomareva A.V., Isaev E.I., Vekilov Yu.H., Abrikosov I. Site preference and effect of alloying on elastic properties of ternary B2 NiAl-based alloys. Physical Review B, 2012, vol. 85, p. 144117.
[14] Lu Y.L., Hou H.X, Chen Z., et al. First-principles calculation on mechanical and thermal properties of B-NiAl with point defects (in Chinese). Sci Sin-Phys Mech Astron, 2013, vol. 43, p. 152.
[15] Hongzhi Fu, Xiao Feng Li, Wen Fang Liu, Yanming Mab, Tao Gao, Xinhua Hong. Electronic and dynamical properties of NiAl studied from first principles. Intermetallics, 2011, vol. 19, p. 1959.
[16] Desai P.D. Thermodynamic Properties of Selected Binary Aluminum Alloy Systems. J. Phys. Chem. Ref. Data. 1987, vol. 16, no. 1.
Статья поступила в редакцию 13.10.2014
Краюшкина Елена Юрьевна — старший преподаватель кафедры "Общая физика" Брянского государственного технического университета. Автор одной научной работы в области физики конденсированного состояния.
Брянский государственный технический университет, Российская Федерация, 241035, Брянск, бул. 50-летия Октября, д. 7.
Krayushkina E.Yu. — senior lecturer of General Physics department of the Bryansk State Technical University. Author of one publication in the field of condensed matter physics. Bryansk State Technical University, bul. 50-letiya Oktyabrya 7, Bryansk, 241035 Russian Federation.
Холодовский Владимир Евгеньевич — канд. физ.-мат. наук, доцент кафедры "Математический анализ" Брянского государственного университета. Автор 46 научных работ в области физики конденсированного состояния, прикладной математики. Брянский государственный университет, Российская Федерация, 241036, Брянск, Бе-жицкая ул., д. 14.
Holodovsky V.E. — Cand. Sci. (Phys.-Math.), assoc. professor of "Mathematical Analysis" department of the Bryansk State University. Author of 46 publications in the field of condensed matter physics, applied mathematics.
Bryansk State University, Bezhitskaya ul. 14, Bryansk, 241036 Russian Federation. Просьба ссылаться на эту статью следующим образом:
Краюшкина Е.Ю., Холодовский В.Е. Расчет теплоемкости и среднеквадратичных смещений твердого раствора NiAl по фононным спектрам // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2015. № 3. C. 110-119. Please cite this article in English as:
Krayushkina E.Yu., Holodovsky V.E. The heat capacity and mean-square displacements of the solid solution NiAl over the phonon spectrum. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 2015, no. 3, pp. 110-119.
