Малая доля дефектов может упорядочивать кристалл

Язык труда и переводы:
УДК:
548.4
Дата публикации:
18 ноября 2022, 01:27
Категория:
Математическое моделирование физических процессов и технических систем
Авторы
Аннотация:
На основе расчетов параметров активационных процессов для золота и железа показано, что если на создание дефекта (вакансии или диффундирующего атома) требуется энергия меньше определенного значения (hi < hsi), или возникающий дефект имеет объем, меньше определенной величины (vi < vsi), то такой дефект имеет отрицательную энтропию, т.е. этот дефект упорядочивает кристалл. Показано, что при малой концентрации дефектов они упорядочивают кристалл, и только начиная с определенной концентрации (Xsi) энтропия образования дефектов переходит в положительную область, где дефекты разупорядочивают кристалл. Изучено изменение функций hsi, vsi и Xsi с ростом давления.
Ключевые слова:
вакансия, самодиффузия, энтропия, энтальпия, упорядочение, золото, железо
Основной текст труда

Введение

Процесс образования дефектов в кристалле (вакансий и диффундирующих атомов) – это одно из загадочных явлений природы. Почему при изобарическом нагреве или при изотермическом растяжении в кристалле появляются дефекты? Всегда ли дефекты разупорядочивают термодинамически устойчивый кристалл? Эти вопросы нуждаются в изучении. Поэтому изучению дефектов в кристалле в последнее время посвящено много работ [1-4].

В [5] нами был предложен аналитический метод (т. е. без компьютерного моделирования) для расчета параметров образования электронейтральных вакансий и самодиффузии атомов в однокомпонентном кристалле. В данной работе использован этот метод для изучения корреляции энтропии активационного процесса с энтальпией и с объемом активационного процесса при различных температурах T и давлениях P в золоте и в железе.

Метод расчета

В [5] нами были получены выражения для энергии Гиббса ( g i), энтальпии ( h_{i} ), энтропии ( s_{i} ) и объема ( v i), как для процесса образования электронейтральной вакансии ( i=v ), так и для процесса самодиффузии атома ( i=d ) по объему кристалла.

Для процесса образования вакансии:

g_{v}=-k_{B}T\ln(\phi {}_{v})=E_{v}\left[{1+\left({\frac {k_{B}T}{2E_{v}}}\right)\ln \left({\frac {\pi E_{v}}{k_{B}T}}\right)}\right] ,

h_{v}=E{}_{v}\,\left\{{1-t_{y}(y_{w})+\alpha _{p}T\left[{(2-t_{y}(y_{w}))\gamma _{\rm {o}}-{\frac {2}{3}}}\right]}\right\} ,                                                                                    (1)

{\frac {s_{v}}{k_{B}}}={\frac {h_{v}-g_{v}}{k_{B}T}}={\frac {E{}_{v}}{k_{B}T}}\,\left\{{\alpha _{p}T\left[{(2-t_{y}(y_{w}))\gamma _{\rm {o}}-{\frac {2}{3}}}\right]-t_{y}(y_{w})-\left({\frac {k_{B}T}{2E_{v}}}\right)\ln \left({\frac {\pi E_{v}}{k_{B}T}}\right)}\right\} ,

{\frac {v_{v}}{v_{0}}}={\frac {E_{v}}{B_{T}v_{0}}}\left[{(2-t_{y}(y_{w}))\gamma _{\rm {o}}-{\frac {2}{3}}}\right] .

Для процесса самодиффузии:

g_{d}=-k_{B}T\ln(x{}_{d})=E_{d}\left[{1-\left({\frac {k_{B}T}{2E_{d}}}\right)\ln \left({\frac {4E_{d}}{\pi k_{B}T}}\right)}\right]

h_{d}=E{}_{d}\,\left\{{1-t_{y}(y_{w})+\alpha _{p}T\left[{(2-t_{y}(y_{w}))\gamma _{\rm {o}}-{\frac {2}{3}}}\right]}\right\} ,                                                                                    (2)

{\frac {s_{d}}{k_{B}}}={\frac {h_{d}-g_{d}}{k_{B}T}}={\frac {E_{d}}{k_{B}T}}\left\{{\alpha _{p}T\left[{(2-t_{y}(y_{w}))\gamma _{\rm {o}}-{\frac {2}{3}}}\right]-t_{y}(y_{w})+\left({\frac {k_{B}T}{2E_{d}}}\right)\ln \left({\frac {4E_{d}}{\pi k_{B}T}}\right)}\right\} ,

{\frac {v_{d}}{v_{0}}}={\frac {E_{d}}{B_{T}v_{0}}}\left[{(2-t_{y}(y_{w}))\gamma _{\rm {o}}-{\frac {2}{3}}}\right] .

Здесь k_{B}  – постоянная Больцмана; \phi _{v} и x_{d} – вероятности образования вакансии и делокализации атома; E_{v} – энергия создания вакантной ячейки в безвакансионной решетке;  E_{d} – энергии делокализации атома в объеме кристалла; \alpha _{P}={\frac {1}{V}}\left({\frac {\partial V}{\partial T}}\right)_{P} – коэффициент теплового расширения, B_{T}=-V\left({\frac {\partial P}{\partial V}}\right)_{T} – изотермический модуль упругости;  \gamma _{0}=-\left[{\frac {\partial \ln {(\Theta _{0})}}{\partial \ln {(V)}}}\right]_{T} – параметр Грюнайзена для безвакансионного кристалла; v_{0} – объем приходящийся на атом при P=0  и T=0K , \Theta _{0} – температура Дебая в безвакансионной решетке,

t_{y}(y_{w})=1-{\frac {2y_{w}\;\exp(y_{w})}{[\exp(2y_{w})-1]}} y_{w}={\frac {3\Theta _{\rm {o}}}{4T}}     .                                                                             (3)

Как бы показано в [5], предложенный метод позволяет, исходя из параметров парного потенциала межатомного взаимодействия Ми — Леннард-Джонса, рассчитать как температурную зависимость вдоль изобары, так и изобарную зависимость вдоль изотермы для всех параметров активационных процессов, указанных в формулах (1) и (2). Этот метод расчета применим при любых давлениях и температурах, соответствующих твердой фазе однокомпонентного вещества.

Результаты

В [5] была рассчитана температурная зависимость параметров активационных процессов для золота от T=10K   до 1330K вдоль двух изобар: P=0 и 24GPa . Аналогичные расчеты были проведены нами для ОЦК и ГЦК фаз железа от T=10K до 1810K вдоль изобар: P=0 и 10GPa . Показано, что при низких температурах из-за квантовых закономерностей активационные параметры сильно зависят от температуры, причем энтропия активационного процесса в этой области отрицательная: s_{i}(T<\Theta _{0})<0 . При T=0K  параметры активационного процесса достигают своих минимумов: g_{i}(0)=0 , h_{i}(0)=0 , v_{i}(0)=0 , s_{i}(0)<0 .

Показано, что такое поведение данных функций согласуется с третьим началом термодинамики в «сильной» формулировке Планка. Отметим, что отрицательное значение энтропии образования дефекта ( s_{i}<0 ) было обнаружено как в теоретических работах [1–4], так и экспериментально [6, 7].

В области высоких температур, т. е. при T\gg \Theta _{0} , вероятность образования вакансии и коэффициент самодиффузии переходят в классические зависимости Аррениуса со слабо зависящей от температуры энтальпией, и с положительным значением энтропии активационного процесса: s_{i}(T\gg \Theta _{0})>0 .

Как было показано в [5] энтальпия, энтропия и объем активационного процесса растут при изобарическом нагреве кристалла. Поэтому они должны коррелировать между собой. В [8] были изучены следующие корреляции параметров активационных процессов: s_{i}(T)\sim h_{i}(T) и s_{i}(T)\sim v_{i}(T) . Эти корреляции были изучены на основе проведенных в [5] расчетов активационных параметров ГЦК- Au . Здесь были использованы также расчеты для ОЦК и ГЦК фаз Fe от T=10K до 1810K вдоль изобар: P=0 и 10GPa . Было показано, что если на создание дефекта (вакансии или диффундирующего атома) требуется энергия меньше определенного значения ( h_{i}<h_{si} ), или возникающий дефект имеет объем, меньше определенной величины ( v_{i}<v_{si} ), то такой дефект имеет отрицательную энтропию, т. е. этот дефект упорядочивает кристалл. При этом выполняются следующие соотношения:

h_{sd}>h_{sv} v_{sd}>v_{sv} h_{si}(P>0)>h_{si}(P=0) v_{si}(P>0)<v_{si}(P=0) .

То есть вакансии начинают разупорядочивать кристалл при меньших энергиях или объемах, необходимых для их создания, чем диффундирующие атомы. С ростом давления энергия, необходимая для создания «разупорядочивающего дефекта», увеличивается, а объем дефекта, необходимый для разупорядочения кристалла, уменьшается. Здесь h_{si} и v_{si} энтальпия и объем образования вакансии при s_{i}=0 .

Зависимость энтропии от логарифма вероятности образования дефекта для ГЦК-Au (слева) и для ОЦК и ГЦК фаз Fe (справа)

Энтальпия активационного процесса связана с вероятностью активационного процесса. Поэтому нами была изучена корреляция энтропии с вероятностью активационного процесса, т. е. с концентрацией дефектов в кристалле. На рисунке показана рассчитанная для ГЦК- Au (слева) и для ОЦК и ГЦК фаз Fe (справа) зависимость энтропии от десятичного логарифма вероятности образования дефекта; s_{i}\sim \lg {(X_{i})} . Здесь для процесса образования вакансии i=v , X_{i}=\phi _{v} , для процесса самодиффузииатома i=d , X_{i}=x_{d} . Две верхние кривые – для процесса самодиффузии атома, две нижние кривые – для процесса образования вакансии.  Сплошные кривые – изобары P=0 , штриховые кривые – изобары P=24GPa – для ГЦК- Au и P=10GPa – для Fe . Изломы на изобарах железа указывают области фазовых переходов из ОЦК в ГЦК структуру.

Из рисунка видно, что при малой концентрации дефектов они имеют отрицательную энтропию, т. е. дефекты упорядочивают кристалл. Только начиная с определенной концентрации ( X_{si} ) энтропия образования дефектов переходит в положительную область, где дефекты разупорядочивают кристалл. Из рисунка видно, что диффундирующие атомы начинают разупорядочивать кристалл при меньших концентрациях, чем вакансии: X_{sd}<X_{sv} . При этом, с ростом давления величина X_{sd} уменьшается, а величина X_{sv} возрастает. X_{si}(P_{i}>P)<X_{si}(P=0) X_{si}(P>P_{i})>X_{si}(P=0) .

Как видно из рисунка энтропия образования дефекта в железе становится отрицательной при очень малой концентрации дефектов:

  • для Au : X_{i}<10^{-20} (для вакансий) – 10^{-60} (для диффундирующих атомов);
  • для Fe : X_{si}<10^{-8} (для вакансий) – 10^{-30}  (для диффундирующих атомов).

Но для таких кристаллов как, например, 3He, 4He, H2, D2, Ne, Li, смена знака энтропии образования вакансии происходит при более заметной концентрации дефектов [9]. Именно поэтому отрицательное значение энтропии образования вакансии было экспериментально обнаружено в кристаллах гелия [7, 8].

Заключение

  1. Если при изобарическом нагреве на создание дефекта (вакансии или мигрирующего атома) требуется малое значение энергии ( h_{i}<h_{si} ) или этот дефект имеет малый объем ( v_{i}<v_{si} ), то этот дефект имеет отрицательную энтропию, т. е. возникающий в этих условиях дефект упорядочивает кристалл. Здесь h_{si} и v_{si} энтальпия и объем образования вакансии при s_{i}=0 .
  2. Показано, что вакансии начинают разупорядочивать кристалл при меньших энергиях или объемах, необходимых для их создания, чем диффундирующие атомы. С ростом давления энергия, необходимая для создания «разупорядочивающего дефекта», увеличивается, а объем дефекта, необходимый для разупорядочения кристалла, уменьшается.
  3. Только начиная с определенной концентрации ( X_{si} ) энтропия образования дефектов переходит в положительную область, где дефекты разупорядочивают кристалл, причем диффундирующие атомы начинают разупорядочивать кристалл при меньших концентрациях, чем вакансии.
Литература
  1. Mendelev M.I., Bokstein B.S., Molecular dynamics study of self-diffusion in Zr. Philosophical Magazine, 2010, vol. 90, no. 5, pp. 637–654. DOI: https://doi.org/10.1080/14786430903219020
  2. Hargather C.Z., Shang S.-L., Liu Z.-K., Du Y. A first-principles study of self-diffusion coefficients of fcc Ni. Computational Materials Science, 2014, vol. 86, pp. 17–23. DOI: https://doi.org/10.1016/j.commatsci.2014.01.003
  3. Kobelev N.P., Khonik V.A. On the enthalpy and entropy of point defect formation in crystals. Journal of Experimental and Theoretical Physics, 2018, vol. 126, no. 3, pp. 340–346. DOI: https://doi.org/10.1134/S1063776118030032
  4. Saltas V., Chroneos A., Vallianatos F. Thermodynamic modelling of fast dopant diffusion in Si. Journal of Applied Physics, 2018, vol. 123, no. 16, art. no. 161527. DOI: https://doi.org/10.1063/1.5001755
  5. Magomedov M.N. Dependencies of the parameters of vacancy formation and self-diffusion in a single-component crystal on temperature and pressure. Journal of Physics and Chemistry of Solids. 2022, vol. 165, art. no 110653. DOI: https://doi.org/10.1016/j.jpcs.2022.110653
  6. Granfors P.R., Fraass B.A., Simmons R.O. Direct measurements of thermal vacancies in bcc 4He. Journal of Low Temperature Physics, 1987, vol. 67, no. 5/6, pp. 353–375. DOI: https://doi.org/10.1007/BF00710349
  7. Iwasa I. Vacancies in bcc 3He at high molar volumes. Journal of the Physical Society of Japan, 1987, vol. 56, no. 5, pp. 1635–1637. DOI: https://doi.org/10.1143/JPSJ.56.1635
  8. Magomedov M.N. A small fraction of defects can order a crystal. Journal of Physics and Chemistry of Solids, 2022, vol. 171, art. no. 111017. DOI: https://doi.org/10.1016/j.jpcs.2022.111017
  9. Magomedov M.N. On the role of vacancies in the process of self-diffusion at low temperatures. Technical Physics Letters, 2002, vol. 28, no. 5, pp. 430–433. DOI: https://doi.org/10.1134/1.1482758
Ваш браузер устарел и не обеспечивает полноценную и безопасную работу с сайтом.
Установите актуальную версию вашего браузера или одну из современных альтернатив.