Добавил:
Когда то был здесь Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Sb97296

.pdf
Скачиваний:
4
Добавлен:
17.12.2025
Размер:
2.34 Mб
Скачать

Для реальных растворов: μ

 

0

RT ln xi RT ln i .

i неид

μi

 

 

 

Химические потенциалы компонентов в неидеальных растворах двухкомпонентной системы А–В в соответствии с (4.34) определяются выражениями:

 

 

 

 

 

 

(l)

 

0(l)

RT ln a

(l)

;

 

(s)

 

0(s)

RT ln a

(s)

;

 

 

 

 

 

A

A

 

A

A

A

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

(l)

 

0(l)

RT ln a

(l)

;

 

(s)

 

0(s)

RT ln a

(s)

,

 

 

 

 

 

B

B

 

B

B

B

 

 

 

 

 

 

 

B

 

 

 

 

 

где

 

0(l)

и

 

0(s)

– химические потенциалы чистых компонентов в жидкой и

i

i

 

 

твердой фазах соответственно. Химический потенциал чистого компонента

μ

0

 

 

 

 

 

 

 

 

 

 

 

 

 

μ

0

G

0

i

равен мольной свободной энергии чистого компонента, т. е.

i

m,i .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая условия фазового равновесия (6.1), получаем

 

 

 

 

 

 

 

a

(s)

 

 

0(l)

 

0(s)

 

 

 

G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

 

 

A

A

 

A

пл A

,

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

a

(l)

 

 

RT

 

 

RT

RT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

A

G

 

плA

 

a

(s)

 

 

0(l)

 

0(s)

 

 

 

G

 

 

 

 

 

 

 

 

 

ln

 

 

B

B

 

B

пл B

,

B

 

 

 

 

 

a(l)

 

 

RT

 

 

RT

RT

 

 

 

B

 

 

 

 

 

 

 

 

 

 

и B GплB

при температуре T.

 

Изменение свободной энергии Гиббса в результате плавления компонента при температуре, отличной от температуры плавления, равно:

 

 

 

 

 

 

T

 

 

 

 

 

T C

 

G

H

 

T S H

 

 

 

C

 

dT T ( S

 

 

 

p

dT ).

T

пл

p

пл

 

T

 

 

 

пл

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tпл

 

 

 

 

 

Tпл

 

 

Используя приближение Улиха C p

Gпл Hпл T Sпл

0

, получаем:

Hпл T Hпл .

Tпл

Обобщенные уравнения, характеризующие равновесия в двухкомпонентной системе:

 

as

 

Hпл(А)

1

 

1

 

ln

A

 

 

 

 

 

 

,

 

R

 

T

 

 

l

 

T

 

 

 

A

 

 

 

 

 

пл(А)

 

 

a

 

 

 

 

 

 

 

91

ln

a

s

 

H

пл(B)

 

1

 

1

 

 

 

 

.

B

 

R

 

T

 

 

l

 

 

T

 

 

 

B

 

 

 

 

 

 

пл(B)

 

 

a

 

 

 

 

 

 

 

 

Учитывая связь между активностью и мольной долей (4.35), получаем систему уравнений, описывающую взаимосвязь между составами сосуществующих жидких (хl) и твердых (xs) реальных растворов и температурой:

где

ki0

=

xis

xl i

 

0

 

 

 

x

s

 

 

H

пл(А)

 

1

 

 

1

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln k

 

ln

 

A

 

 

 

 

 

 

 

ln

 

A

,

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

R

T

 

T

 

 

s

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

пл(А)

 

A

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.2)

 

 

 

 

 

 

x

 

 

H

 

 

1

 

 

1

 

 

 

 

 

0

 

 

s

 

пл(B)

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln k

ln

 

B

 

 

 

 

 

 

 

ln

 

 

B

,

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

R

T

 

 

T

 

 

s

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

пл(B)

 

B

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

– равновесный коэффициент распределения i-го компонента

между твердым и жидким растворами. Равновесный коэффициент распреде-

ления

ki0

характеризует отношение концентраций компонента в твердой и

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

Для идеальных растворов (γi = 1) система уравнений (6.2) преобразуется к виду:

 

 

 

 

 

 

0

 

 

 

 

x

s

 

 

H

пл(А)

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln k

ln

 

A

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

R

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пл(А)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.3)

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

H

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

s

 

 

пл(B)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln k

ln

B

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

R

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пл(B)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая, что x s

x s

 

1,

 

 

xl xl

 

1, получаем уравнения для линий

 

 

 

 

A

 

 

B

 

 

 

 

 

 

A

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ликвидуса и солидуса соответственно:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xs

xs

xl

H

пл(А)

 

1

 

 

 

1

 

 

xl

 

 

 

 

H

пл(B)

 

1

 

1

 

exp

 

 

 

 

 

 

exp

 

1,

 

R

 

 

 

 

T

 

 

 

R

 

 

T

A

B

A

 

 

 

 

T

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пл(А)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пл(B)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

92

xl

xl

xs

 

 

H

пл(А)

 

1

 

 

1

 

exp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

B

 

A

 

 

 

 

 

 

R

 

T

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пл(А)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xs

 

 

H

пл(B)

 

1

 

 

1

 

 

 

 

 

exp

 

 

 

 

 

 

1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

R

 

T

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пл(B)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.4)

Уравнения (6.3) и (6.4) позволяют рассчитать и построить Т-х-проекцию диаграммы состояния системы АВ в рамках модели идеальных растворов на основании известных термодинамических данных для исходных компонентов А и В.

Для модели регулярных растворов уравнения (6.2) с учетом (4.41) преобразуются к виду:

 

x

s

 

H

пл(А)

1

 

 

 

1

 

 

W

l

(1 x

l

 

)

2

 

W

s

(1 x

s

)

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

A

 

 

 

,

 

 

 

 

R

 

 

T

 

 

 

 

RT

 

 

 

 

 

 

 

 

 

 

RT

 

 

 

 

 

 

 

l

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

пл(А)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

s

 

H

пл(B)

 

1

 

 

1

 

 

 

W

l

(1

x

l

 

)

2

 

 

W

s

(1 x

s

 

)

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

B

 

 

.

 

 

 

l

 

 

 

R

 

T

 

 

T

 

 

 

 

 

 

RT

 

 

 

 

 

 

 

 

 

 

RT

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

пл(B)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.5)

Группируя (6.5) для модели регулярных растворов и принимая во внима-

ние, что

x

s

x

s

1

и

x

l

x

l

1

, запишем систему уравнений, описываю-

A

B

A

B

 

 

 

 

 

 

 

 

 

щих линии солидусa и ликвидусa, относительно температуры:

 

H

пл

(В) W l (1 xl )2

W s (1 xs

)2

 

T

 

 

 

 

 

B

 

 

 

B

 

,

 

 

 

S

пл

(B) R ln xs

/ xl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

B

 

 

 

 

 

 

H

пл

(A) W l (xl )2

W s (xs

)2

 

 

T

 

 

 

 

B

 

B

 

 

.

 

 

 

 

 

 

R ln(1 xs ) /(1 xl

 

 

 

 

 

S

пл

(A)

 

)

 

 

 

 

 

 

 

 

 

B

B

 

 

 

Уравнения позволяют решать 2 задачи:

(6.6)

1)по известным значениям Ws и Wl, энтальпиям и температурам плавления исходных компонентов рассчитать Т-х-проекцию диаграммы состояния бинарной системы;

2)по экспериментальным значениям составов равновесных жидкой и твердой фаз определять параметры взаимодействия корреляционным методом, исходя из минимума среднеквадратического отклонения (СКО) экспери-

ментальных значений ( xэкспl и xэкспs ) от теоретических ( xтеорl и xтеорs ), рассчитанных для задаваемых Wl и Ws.

93

Рассмотрим причины появления на диаграмме с непрерывной растворимостью компонентов экстремальной точки конгруэнтного плавления, в которой составы твердой и жидкой фаз совпадают. Условие термодинамической

устойчивости для составов твердой и жидкой фаз:

2

G

M

 

 

x

2

 

 

> 0. Для модели

регулярного раствора после дифференцирования выражения для

G

M

 

из

табл. 4.1 получим:

RT

 

2W

 

x

 

 

x 1

 

 

> 0 (подробнее явления устойчивости и рас-

пада растворов рассмотрены в 6.5, 6.6).

Это неравенство выполняется для любых отрицательных параметров

взаимодействия,

При W >

2

 

а для положительных – при условии W <

RT

.

2x(1

x)

 

 

 

RT

наблюдается распад твердого раствора или расслое-

x(1 x)

 

 

 

 

ние жидкой фазы.

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

Из уравнений (6.6), описывающих взаимосвязь между составами сосуществующих жидких и твердых растворов и температурой для модели регу-

лярного раствора, при xis xil xi экст получим:

 

 

 

 

 

 

W

s

W

l

 

 

2 T

W

s

W

l

 

 

2 .

T

 

T

 

 

 

 

1 x

 

 

 

x

экст

 

пл(В)

S

 

 

 

экст

пл(A)

S

 

 

экст

 

 

 

 

 

 

 

 

пл(В)

 

 

 

 

пл(A)

 

 

 

При

W

s

>

W

l

(взаимодействие между атомами А и В в твердой фазе

 

 

 

 

 

 

слабее, чем в жидкой) Tэкст < Tпл(A) и Tэкст < Tпл(В) . Это наиболее вероятно для положительного параметра взаимодействия компонентов в твердой фазе:

W

s

> 0. Примером такого взаимодействия является диаграмма состояния си-

 

 

 

стемы «медь–золото» (рис. 6.2).

 

 

 

Для таких диаграмм состояния при Т < Tэкст

часто возникает область

распада твердого раствора. Это связано с

тем, что в уравнении

G M H M TS M с уменьшением температуры при положительном парамет-

94

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

t, C 1083

1063

905

 

 

 

 

Cu

xэкст

Au

Рис. 6.2. Диаграмма состояния системы Cu–Au

T

Tэкст

T

пл

(B)

 

 

T

пл

(A)

 

 

A

x

экст

= xl = xs

B

 

 

 

 

Рис. 6.3. Диаграмма состояния системы Pb–Tl

Для отрицательных параметров взаимодействия компонентов в твердой

фазе

W

s

< 0,

Tэкст

> Tпл(A)

и

Tпл(В)

. Примером является диаграмма состоя-

ния в системе «свинец–таллий» (рис. 6.3).

6.2. Общие сведения о фазовых диаграммах состояния «полупроводник–примесь»

Одноэлементные полупроводники, такие, как кремний и германий, а также алмазоподобные соединения типа А2В6, А3В5 образуют с большинством химических элементов диаграммы состояния эвтектического типа, для многих из которых состав, соответствующий точке эвтектики, сдвинут относительно примеси, а для некоторых характерна вырожденная эвтектика

(рис. 6.4, б) [15].

95

t, C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

750

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

500

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L +

 

 

 

 

 

 

 

250

 

 

 

L +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

+

 

29,8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1018

 

1019

1020

1021 0

 

 

 

 

 

 

 

 

 

1017

 

10

20

30

40

50

60

70

80

90

 

N

Ga

, см–3

 

 

 

 

 

Ga, ат. %

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

б

 

 

 

 

Рис. 6.4. Крупномасштабная диаграмма в области твердого раствора примеси

 

в полупроводнике (а) и Т-х-проекция диаграммы состояния системы Ge–Ga (б)

Вследствие того, что германий и кремний имеют ковалентные, направленные и насыщенные связи, растворимость в них химических элементов с другим характером связи по механизму замещения очень ограничена. Атомная (мольная)* доля многих примесей (хВ) в твердых растворах с германием и

кремнием не превышает 10–1…10–3. Из-за очень низкой растворимости примесей область твердого раствора (α) со стороны полупроводника можно отобразить только с помощью крупномасштабной диаграммы состояния (рис. 6.4, а). На таких диаграммах по оси абсцисс откладывается содержание примеси не в атомных (мольных) долях, как на обычной диаграмме, а в виде концентрации атомов в единице объема [см–3 или м–3]. Пересчет из одних единиц измерения концентрации примеси в другие осуществляется на основании справочных данных о структурном типе кристалла, периоде решетки a и кратности элементарной ячейки (число атомов, принадлежащих одной элементарной ячейке кристалла полупроводника). Учитывая, что объем элементарной ячейки для кубических структур равен V = a3, рассчитываем общее число атомов в единице объема полупроводника по формуле: n = t/V, где

* Мольные доли (или мольные проценты) используются в том случае, когда компонентами системы являются химические соединения или другие многокомпонентные фазы.

96

t – кратность ячейки. Далее, зная атомную долю примеси в растворе

xB

,

можно определить концентрацию примеси (в см–3) по формуле: NB

= n

xB

.

Следует обратить внимание на то (рис. 6.4, а), что максимальная растворимость примесей наблюдается не при температурах эвтектики, как это происходит в металлических твердых растворах, а при значительно более высоких температурах. Такой характер растворимости, когда содержание примеси в твердом растворе уменьшается с понижением температуры, называется ретроградным.

Для математического описания диаграмм состояния с ограниченной растворимостью компонентов (эвтектического и перитектического типов) используют различные модели неидеальных растворов. Химический потенциал i-го компонента в неидеальном растворе представляют в виде

где

0 i

μ

 

0

RT ln x

i

RT ln

i ,

(6.7)

 

i

i

 

 

 

 

 

 

 

 

 

– химический потенциал чистого i-го компонента в стандартном со-

стоянии (при давлении Р = 1 атм и заданной температуре Т); хi – атомная доля компонента в растворе; γi – коэффициент активности i-го компонента.

По сравнению с идеальным раствором в выражении (6.7) появляется дополнительное слагаемое RT ln γi, которое можно характеризовать как меру отклонения раствора от идеального.

Т-х-проекции диаграмм состояния бинарных систем А–В, где А – полупроводник (Si или Ge), а В – примесь, достаточно хорошо описываются в рамках моделей регулярного и квазирегулярного растворов, которые учитывают межатомное взаимодействие в растворах с помощью параметра взаимодействия W. Согласно этим моделям коэффициент активности i-го компонента связан с параметром взаимодействия соотношением

 

 

 

W

1 x

2

ln

i

 

 

 

i

 

 

 

RT

 

 

 

 

 

 

 

 

 

 

 

 

.

(6.8)

В приближении регулярного раствора параметр взаимодействия W считается независимым от температуры и состава раствора. В модели квазирегулярного раствора учитывают зависимость параметра взаимодействия от температуры в виде линейной функции:

W a bT .

(6.9)

97

Анализ диаграмм состояния бинарных систем на основе кремния и германия показал, что применительно к твердой фазе лучшее соответствие с экспериментальными данными дает приближение регулярных растворов, а для описания жидкой фазы в ряде случаев лучше подходит модель квазирегулярных растворов (6.9). В табл. 6.1 представлены параметры межатомного взаимодействия в жидких и твердых фазах для некоторых бинарных систем на основе кремния и германия.

Таблица 6.1

Параметры взаимодействия в некоторых системах «полупроводник–примесь»

Эле-

 

Si

 

Ge

W s,

W l, кДж/моль

W s,

W l, кДж/моль

мент

кДж/моль

кДж/моль

 

 

 

Al

63,16

–17,33 + 5,1.10–3T

16,61

–22,19 + 13,2.10–3T

Ga

49,13

13,6 3,47.10–3T

6,41

–0,62

In

123,44

47,9 – 14,11.10–3T

67,42

6,57 – 2,34.10–3T

P

8,0

16,5

As

24,9

–210,0 + 135,8.10–3T

32,36

–23,46 + 17,43.10–3T

Sb

58,0

13,77 + 6,75.10–3T

52,13

11,05 – 8,29.10–3T

Au

70,41

–81,87 + 43,0.10–3T

106,6

–22,038 + 4,27.10–3T

Cu

206,09

–49,86 + 30,09.10–3T

120,66

–30,8 + 32,1.10–3T

Bi

125,47

62,12 – 8,62.10–3T

103,8

23,02 – 6,24.10–3T

Ag

154,14

–33,11 + 31,94.10–3T

126,63

–23,02 + 29,8.10–3T

Li

46,33

–9,72

52,72

–5,213

Sn

68,10

34,09 – 6,28.10–3T

44,62

7,03 – 4,52.10–3T

Zn

160,75

17,91 + 2,77.10–3T

71,76

–2,95

Pb

98,71

36,75 – 17,08.10–3T

Ni

131,4

–3,07

 

 

 

 

 

В приближении неидеального раствора основными соотношениями, описывающими равновесие между твердой и жидкой фазами, являются:

ln kA0

 

xs

 

Hпл(А)

1

 

 

 

1

 

 

ln lA ln As

 

ln

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

R

 

 

 

 

T

 

 

 

 

l

 

 

 

 

T

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

пл(А)

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln kB0 ln

xs

 

 

Hпл(B)

1

 

 

1

 

 

ln lB ln Bs

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

R

 

 

 

 

T

 

 

 

 

 

l

 

 

 

T

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

пл(B)

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где kA0 и kB0 – равновесные коэффициенты распределения компонентов А и В соответственно в твердом и жидком растворах; xis , xil – атомные доли i -го

98

компонента в твердой и жидкой фазах соответственно;

Hпл(А)

и

Hпл(B)

энтальпии плавления;

Tпл(A)

и Tпл(B)

– температуры плавления компонентов

А и В соответственно;

l

i

жидкой и твердой фазах.

и

s

– коэффициенты активности i-го компонента в

i

Выражая коэффициенты активности через парамет-

ры межатомного взаимодействия компонентов в твердой (

W

s

) и жидкой (

W

l

)

 

 

 

 

 

 

фазах соответственно, с помощью (6.8) получаем соотношения, описывающие фазовые равновесия в двухкомпонентной системе в рамках модели регулярных растворов:

 

x

s

 

Hпл(А) 1

 

1

 

 

W

l

(1 x

l

)

2

 

W

s

(1 x

s

)

2

 

A

 

 

 

 

A

 

 

 

A

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

R

T

 

T

 

 

 

 

RT

 

 

 

 

 

 

 

RT

 

 

 

 

xA

 

 

 

 

 

пл(А)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

s

 

Hпл(B)

1

 

1

 

 

W

l

(1 x

l

)

2

 

W

s

(1 x

s

)

2

ln

B

 

 

 

 

B

 

 

 

B

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

R

T

 

T

 

 

 

 

RT

 

 

 

 

 

 

 

RT

 

 

 

 

 

xB

 

 

 

 

 

пл(B)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.11)

Систему уравнений (6.11) можно записать относительно температуры, учитывая, что хА + хВ = 1 и ∆Hпл = ∆SплTпл :

 

 

H

пл(A)

W

l

(x

l

)

2

W

s

(x

s

)

2

 

 

 

T

 

 

B

 

 

 

 

B

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

S

пл(A)

 

R ln[(1 x

 

) /(1 x

 

)]

 

 

 

 

 

B

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

пл(B)

W

l

(1

x

l

 

)

2

W

s

(1 x

s

)

2

T

 

B

 

 

 

 

 

B

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

пл(B)

 

R ln(x

/ x

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.12)

6.3. Анализ линий ликвидуса и определение параметра взаимодействия компонентов в жидком растворе Wl

Вследствие малой растворимости примесей в твердых германии и кремнии можно при анализе линий ликвидуса считать, что в равновесии с жидким раствором находится кристаллическая фаза чистого полупроводникового

компонента (кремния, германия), т. е. можно принять xBs 0 ,

уравнения (6.12) принимают вид:

x

s

A

 

1

. Тогда

 

H

пл(A)

W l (xl

)2

 

 

 

H

пл(B)

W l (1 xl

)2 W s

 

T

 

B

 

 

,

T

 

 

B

 

 

. (6.13)

S

 

R ln(1 xl

)

S

 

R ln(xs

/ xl

)

 

пл(A)

 

 

B

 

 

 

 

пл(B)

B

B

 

 

99

Уравнения (6.13) описывают линию ликвидуса в предположении регулярности жидкого раствора «полупроводник–примесь» в широком диапазоне

составов от xBl 0 до xBl xB(l эвт) , где xB(l эвт) – состав, соответствующий точке эвтектики.

Используя экспериментальные значения xBl , соответствующие линии ликвидуса на Т-х-проекции диаграммы состояния, можно c помощью (6.13) рассчитать параметр взаимодействия компонентов в жидком растворе Wl для каждого из заданных условий (Т, xBl ). В результате получим ряд значений Wl

в определенном диапазоне

l

l

 

 

Wmin ...Wmax . Разброс значений параметра меж-

атомного взаимодействия Wl зависит от точности определения состава

l

и

xB

температур на Т-х-диаграмме состояния. Значение параметра Wl, обеспечивающее наилучшую сходимость экспериментальной и теоретически рассчи-

танной линий ликвидуса во всем интервале задаваемых составов xl и темпе- B

ратур, определяют, исходя из условия минимума СКО(Т) экспериментальных значений (Тэксп) от теоретических (Ттеор), рассчитанных по уравнению (6.13):

 

 

 

 

1

 

T

Т

теор

2

 

(Т )

 

 

 

СКО

 

 

эксп

 

 

 

n

n 1

 

T

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

эксп

 

,

(6.14)

где n – число экспериментальных точек (значений температур), взятых из Т-х-проекции диаграммы и используемых для расчетов.

T, К

 

 

 

 

 

 

 

 

 

1200

5000

 

 

 

 

 

 

 

 

1100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1000

2000

 

 

 

 

 

 

 

 

900

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

800

 

 

 

 

 

 

 

 

 

700

W l = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

600

–3000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

500

 

 

 

 

 

 

 

 

 

400

–6000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3000

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9 xGe

Рис. 6.5. Характер изменения линий ликвидуса в зависимости от значения параметра взаимодействия Wl для систем «германий–примесь»

100