Изучение хрупко-пластичного перехода в кристаллах полупроводников

Язык труда и переводы:
УДК:
539.56
Дата публикации:
18 ноября 2022, 01:48
Категория:
Математическое моделирование физических процессов и технических систем
Авторы
Аннотация:
Показано, что хрупкость ковалентного кристалла обусловлена «двуличностью» парного потенциала межатомного взаимодействия для упругой (обратимой) и для пластической (необратимой) деформации. Это приводит к тому, что удельная поверхностная энергия при пластической деформации ковалентного кристалла более чем в два раза меньше удельной поверхностной энергии при упругой деформации. Поэтому при малой деформации ковалентного кристалла энергетически выгоднее создать поверхность путем разрыва, чем путем ее упругого растяжения. Показано, что температура хрупко-пластичного перехода для однокомпонентных ковалентных кристаллов имеет верхний предел: TBDT/Tm < 0,45, где Tm — температура плавления.
Ключевые слова:
межатомная ковалентная связь, хрупкость, пластичность, элементарные полупроводники, хрупко-пластичный переход
Основной текст труда

Введение

Хрупкость – это свойство материала разрушаться без образования заметных остаточных деформаций. Элементарные полупроводники: алмаз (C-dia) , кремний ( Si ), германий ( Ge ), характерны своей хрупкостью, что связанно с ковалентным характером химической связи. Температура хрупко-пластичного перехода ( T_{BDT} ) у них разная, но хрупкость сохраняется до темпера­тур, значительно больших относитель­но температуры плавления( T_{m} ), чем у металлов [1]. При этом температура хрупко-пластичного перехода тем выше, чем меньше доля металлической связи: T_{BDT}(C-dia)>T_{]BDT}(Si)>T_{BDT}(Ge) . При температурах выше T_{BDT} материал разрушается пластично, ниже T_{BDT} – хрупко. В зависимости от наличия дефектов и примесей, от скорости деформации и от освещенности образца значение T_{BDT} в элементарных полупроводниках может меняться в широком ин­тервале температур [1, 2]. Например, для Si получено [2]: 0.32\leq T_{BDT}/T_{m}\leq 0.7 .   

Природа хрупкости элементарных полупроводников до сих пор не ясна [1, 2]. Существующие теории хрупкого разрушения основаны на результатах классической теории упругости с малыми деформациями, и в этих теориях, обычно не затрагивается проблема начального возникновения трещины. Между тем, трещина, как правило, зарождается на поверхности деформируемого тела, и поэтому свойства поверхности вещества определяют его хрупко-пластичные свойства.

В данной работе на основе модели парной ковалентной связи предложенной в [3] объяснены причины как возникновения трещины на поверхности кристалла полупроводникапри T<T_{BDT} , так и его переход в пластичное состояние при T>T_{BDT} .

Метод расчета и результаты

Представим парное межатомное взаимодействие в однокомпонентном кристалле в виде потенциала Ми–Леннард-Джонса, который имеет следующий вид:

\phi (r)={\frac {D}{b-a}}\left[a\left({\frac {r_{0}}{r}}\right)^{b}-b\left({\frac {r_{0}}{r}}\right)^{a}\right] ,                                                                                               (1)

где D и r_{o} – глубина и координата минимума потенциала, b>a>1  – параметры.

Глубину потенциальной ямы можно определить двумя способами: из модуля упругости B_{00} , и изэнергии атомизации (на атом) L_{00} при T=0 K и P=0  [3]:

D_{b}={\frac {18B_{00}V_{00}}{k_{b}abN_{a}}}=D_{S}+\Delta D ,   

                                D_{S}={\frac {L_{00}}{k_{n}/2}} ,                                                                                                        (2)

где V – молярный объем; N_{A} – число Авогадро; k_{n} – первое координационное число; B=-V\left({\frac {\partial P}{\partial V}}\right)_{T} – модуль упругости; \Delta D=D_{b}-D_{s}  – разница между результатами.

Для металлов величины D_{b} и D_{s} , как правило, совпадают, т. е. \Delta D=0 . Однако для ковалентных кристаллов D_{b}\gg D_{s} , как это показано в табл. 1. Причем величина \Delta D при переходе от алмаза к Si , Ge и \alpha -Sn   последовательно уменьшается с ростом массы атома ( m ), и для Pb становится близкой к нулю: \Delta D(Pb)\approx 0  [3].

 

Таблица 1

Определенные в [3] параметры межатомного потенциала (1), и рассчитаны из них значения \Delta D и d=\Delta D-(D_{b}/2)

Element

Element

m,a.m.u.

m,a.m.u.

r_{0},A

r_{0},A

b

b

 

a

a

 

D_{b},eV

D_{b},eV

D_{s},eV

D_{s},eV

\Delta D,eV

\Delta D,eV

d,eV

d,eV

\Delta D/D_{s}

\Delta D/D_{s}

C-dia

C-dia

12.01

1.545

3.79

2.21

8.43

3.68

4.75

0.535

1.291

Si Si28.09 

2.351

4.00

2.48

5.54

2.32

3.22

0.450

1.388

Ge

Ge

72.59

2.450

4.30

2.75

4.03

1.94

2.09

0.075

1.077

\alpha -Sn

\alpha -Sn

118.7

2.798

4.43

2.79

3.15

1.56

1.59

0.015

1.019

Pb

Pb

207.2

3.500

14.2

2.38

0.378

0.338

0.04

– 0.149

0.118

Если при бесконечно малом изгибе ковалентного кристалла создается новая поверхность \delta S , то создать ее можно как путем обратимого упругого растяжения (тогда удельная (на единицу площади) поверхностная энергия \sigma _{b} определится из D_{b} ), так и путем необратимого хрупкого разрыва (тогда значение \sigma _{s} определится из D_{s} ). В первом случае необходимо затратить энергию: E_{b}=\sigma _{b}\delta S , а во втором: E_{s}=\sigma _{s}\delta S . Поэтому необходимым условием хрупкого разрыва является неравенство:

\Delta E=E_{b}-E_{s}=(\sigma _{b}-\sigma _{s})\delta S\geq 0 .                                                                                             (3)

Условие (3) указывает на наличие у системы двух энергетических уровней, т. е. помимо упругой деформации у системы есть возможность создать поверхность и путем разрыва. Для равновесного металлического монокристалла, т. е. в котором нет метастабильных напряженных областей, выполняется \Delta E=0 .

Достаточным условием реализации хрупкого разрыва является неравенство:

\Delta E=(\sigma _{b}-\sigma _{s})\delta S\geq \sigma _{s}\delta S .                                                                                           (4)

Условие (4) означает, что выигрыш в энергии при хрупком разрыве оказывается большим, чем энергия необходимая для создания новой поверхности. Поэтому температуру хрупко-пластичного перехода можно определить из условия:

\sigma _{b}(T_{BDT})=2\sigma _{s}(T_{BDT}) .                                                                                                     (5)

Как было показано в [4-6], реальное значение удельной поверхностной энергии для кристалла с ковалентной связью необходимо определять из величины D_{s} . Рассчитанные таким путем значения \sigma _{s} хорошо согласуются с оценками удельной поверхностной энергии из работ других авторов, представленными для каждого кристалла в нижних строках табл. 2. Между тем, рассчитанная исходя из значения D_{b}  «упругая» удельной поверхностная энергия \sigma _{b}(T=0K)  много больше.

Таблица 2

Рассчитанные значения удельной поверхностной энергии грани (100) при T=0K и P=0

Crystal

Crystal

\sigma _{b}10^{-3}J/m^{2}

\sigma _{b}10^{-3}J/m^{2}

\sigma _{s}10^{-3}J/m^{2}

\sigma _{s}10^{-3}J/m^{2}

\sigma _{b}-\sigma _{s}10^{-3}J/m^{2}

\sigma _{b}-\sigma _{s}10^{-3}J/m^{2}

(\sigma _{b}-\sigma _{s})/\sigma _{s}

(\sigma _{b}-\sigma _{s})/\sigma _{s}

C-dia[4]

C-dia[4]

14025.0

6104.5

7920.50

1.297

Data  from other authors

9720ideal – 5710recons1) [7],

 6231.9ideal – 4969.1recons1) [8]

 

 

Si[5]

Si[5]

4001.10

1673.42

2327.68

1.391

Data from other authors

2390ideal – 1410recons1) [7],

1879.8ideal – 1010.9recons1) [8]

 

 

Ge[6]

Ge[6]

2681.50

1290.32

1391.18

1.078

Data from other authors

1710ideal – 1000recons1) [7],

 1657.2ideal – 883.1recons1) [8]

 

 

1) Расчетные поверхностные энергии идеальных (нерелаксированных) и реконструированных (релаксированных) поверхностей(100).

Значение \sigma _{b} рассчитано из потенциала (1) со значением D_{b} , а величина \sigma _{s} – со значением D_{s} из табл. 1. В нижних строках представлены оценки удельной поверхностной энергии грани (100) из [7, 8].

При изобарическом росте температуры в результате термической металлизации ковалентной связи функция \sigma _{b}(T)  уменьшается:

\sigma _{b}(T)=\sigma _{b}(0)\left\{1-\left[1-{\frac {\sigma _{s}(T_{m})}{\sigma _{b}(0)}}\right]^{T_{m}/T}\right\} .                                                                                  (6)

Так как при T_{BDT}  выполняется условие (5), то из (6) легко получить формулу:

{\frac {T_{BDT}}{T_{m}}}=\ln {\left[1-{\frac {\sigma _{s}(T_{m})}{\sigma _{b}(0)}}\right]}/\ln {\left[1-{\frac {\sigma _{s}(T_{BDT})}{\sigma _{b}(0)}}\right]}={\frac {\ln {(1-x)}}{\ln {(1-2x+\Delta })}} .                                                      (7)

Здесь введены обозначения:

x={\frac {\sigma _{s}(T_{m})}{\sigma _{b}(0)}}<1 ,           \Delta ={\frac {2[\sigma _{s}(T_{BDT})-\sigma _{s}(T_{m})]}{\sigma _{b}(0)}}\ll 1 .                                                                 (8)

В табл. 3 представлены рассчитанные с помощью (7) значения T_{BDT}  при \Delta =0 . Величины \sigma _{b}(0) были взяты из табл. 2. Определение \sigma _{s}(T_{m}) сделано двумя способами.

  1. Расчет величины \sigma _{s}(T_{m})_{s}  для твердого бездефектного кристалла методом из [4-6]. Эти данные представлены для каждого кристалла в первой строке табл. 3.
  2. Использование экспериментального значения для жидкой металлизованной фазы \sigma (T_{m})_{l}<\sigma _{s}(T_{m})_{s} . Эти данные из [9] представлены во вторых строках табл. 3.

Таблица 3

Значения значения \sigma _{s}(T_{m}) , рассчитанные для твердой фазы (первая строка) и экспериментальные оценки для жидкой фазы (вторая строка), значения T_{BDT} рассчитанные по формуле (7) при \Delta =0 и экспериментальные оценки T_{BDT}

Crystal

Crystal

T_{m},K

T_{m},K

\sigma _{s}(T_{m})10^{-3}J/m^{2}

\sigma _{s}(T_{m})10^{-3}J/m^{2}
\sigma _{s}(T_{m})/\sigma _{b}(0) \sigma _{s}(T_{m})/\sigma _{b}(0) T_{BDT}/T_{m} T_{BDT}/T_{m}

T_{BDT},K

T_{BDT},K

calculated 

T_{BDT},K

T_{BDT},K

Experimental [1]

C-dia

C-dia

4235

4975

0.355

0.354

1501

1470-1615

~ 2805

~ 0.2

0.437

1850

Si

Si

1685

1565

0.391

0.325

548

820-1225

746 [9]

0.186

0.442

745

Ge

Ge

1210

1216

0.455

0.254

308

630-695

605 [9]

0.226

0.426

516

Из данных табл. 3 видно, что хорошее согласие с экспериментальными значениями для T_{BDT} получается при использовании в (7) экспериментального значения \sigma _{s}(T_{m})_{l} для жидкой фазы. Это говорит о том, что твердая фаза указанных кристаллов при T_{m} еще не полностью металлизирована. При этом нужно учесть, что данные расчеты проведены для термодинамически равновесных состояний, т. е. без учета скорости деформации. Как известно из экспериментов [10], с ростом скорости деформации температура хрупко-пластичного перехода возрастает.

 Изменение функции (7) при различных значениях Δ из (8)

На рисунке показана зависимость функции (7) от аргумента x при значениях \Delta из (8) от \Delta =0 до \Delta =0.5 . Из (7) и рисунка следует, что для ковалентных кристаллов выполняются условия:

max\left({\frac {T_{BDT}}{T_{m}}}\right)=0,45 ,       max\left({\frac {\sigma _{s}(T_{BDT})}{\sigma _{b}(0)}}\right)=0,5 .                                                                     (9)

Однако условия (9) получены для хрупко-пластичного перехода под статической нагрузкой. С ростом скорости деформации значение T_{BDT}/T_{m}  возрастает [10].

Заключение

Показано, что хрупкость ковалентного кристалла обусловлена «двуличностью» парного потенциала межатомного взаимодействия для упругой (обратимой) и для пластической (необратимой) деформации.

Показано, что для ковалентного кристалла удельная поверхностная энергия при пластической деформации более чем в два раза меньше чем при упругой деформации. Поэтому при малой деформации ковалентного кристалла энергетически выгоднее создать поверхность путем разрыва, чем путем ее упругого растяжения.

Показано, что температура хрупко-пластичного перехода для однокомпонентных ковалентных кристаллов имеет верхний предел: T_{BDT}/T_{m}<0.45 .

Литература
  1. Giannattasio A., Roberts S.G. Strain-rate dependence of the brittle-to-ductile transition temperature in tungsten. Philosophical Magazine, 2007, vol. 87, no. 17, pp. 2589–2598. DOI: https://doi.org/10.1080/14786430701253197
  2. Cheng G., Zhang Y., Chang T.-H., Liu Q., Chen L., Lu W. D., Zhu T., Zhu Y. In situ nano-thermo-mechanical experiment reveals brittle to ductile transition in silicon nanowires. Nano Letters, 2019, vol. 19, no. 8, pp. 5327–5334. DOI: https://doi.org/10.1021/acs.nanolett.9b01789
  3. Magomedov M.N. On the Nature of the Covalent Bond in Crystals of Carbon Family Elements. Russian J. of Inorganic Chemistry, 2004, vol. 49, no. 12, pp. 1906–1916.
  4. Magomedov M.N. Variations in Thermal Properties of Diamond under Isothermal Compression. Technical Physics, 2017, vol. 62, no. 5, pp. 661–668. DOI: https://doi.org/10.1134/S1063784217050176
  5. Magomedov M.N. Changes of the Thermodynamic Properties at Isochoric and Isobaric Decrease of the Silicon Nanocrystal Size. Physics of the Solid State, 2019, vol. 61, no. 4, pp. 642–649. DOI: https://doi.org/10.1134/S106378341904019X
  6. Magomedov M.N. Change in thermodynamic properties of a germanium nanocrystal with a decrease in its size under various P-T-conditions. Nanotechnologies in Russia, 2019, vol. 14, no. 1-2. pp. 21–32. DOI: https://doi.org/10.1134/S1995078019010063
  7. Stekolnikov A.A., Bechstedt F. Shape of free and constrained group-iv crystallites: influence of surface energies. Physical Review B, 2005, vol. 72, no. 12, art. no. 125326. DOI: https://doi.org/10.1103/PhysRevB.72.125326
  8. Zhang J.M., Li H.Y., Xu K.W., Ji V. Calculation of surface energy and simulation of reconstruction for diamond cubic crystals (001) surface. Applied Surface Science, 2008, vol. 254, no. 13, pp. 4128–4133. DOI: https://doi.org/10.1016/j.apsusc.2007.12.049
  9. Alchagirov B.B., Taova T.M., Khokonov Kh.B. Surface energy and surface tension of solid and liquid metals. Recommended values. Transactions of JWRI. Special Issue (Japan), 2001, vol. 30, pp. 287–291. URL: https://repository.exst.jaxa.jp/dspace/handle/a-is/48071 (дата обращения 09.09.2022).
  10. Evstifeev A.D., Gruzdkov A.A., Petrov Y.V. Dependence of the type of fracture on temperature and strain rate. Technical Physics, 2013, vol. 58, no. 7, pp. 989–993. DOI: https://doi.org/10.1134/S1063784213070086
Ваш браузер устарел и не обеспечивает полноценную и безопасную работу с сайтом.
Установите актуальную версию вашего браузера или одну из современных альтернатив.