Санкт-Петербургский государственный
МИНОБРНАУКИ РОССИИ
Санкт-Петербургский государственный
электротехнический университет
«ЛЭТИ» им. В.И. Ульянова (Ленина)
Кафедра МНЭ
отчет
по лабораторной работе №2
по дисциплине «Физико-химические основы технологии изделий электроники и наноэлектроники»
Тема: Моделирование диаграмм состояния с ограниченной растворимостью компонентов (полупроводник - примесь)
Студенты гр. 6209 | Тойкка А.С. | |
Богданова М.В. | ||
Гребнев К.В. | ||
Преподаватель | Мараева Е. В. |
Санкт-Петербург
2018
Цель работы.
Изучение фазовых равновесий в бинарной системе Ge-Ag путём компьютерного моделирования; расчет параметров межатомного взаимодействия; выбор моделей растворов, адекватно описывающих экспериментальные Т-х-диаграммы состояния; расчет спинодалей и определение областей устойчивого и метастабильного состояний твердых растворов в заданной системе.
Исходные данные
Исходная T-x-проекция фазовой диаграммы состояния системы A-B (Ge-Ag) представлена на рисунке 1.
Рисунок 1 - Т-х-проекция диаграммы состояния системы Ge-Ag
Обработка экспериментальных результатов
1. Привести Т-х-проекции исследованной диаграммы состояния. Оптимальные параметры межатомного взаимодействия Wsopt и Wlopt в твердом и жидком растворах соответственно.
Рис. 2 . T-x-проекции экспериментальной и теоретически рассчитанной диаграмм состояния
Как видно из графика, T-x-проекция теоретически рассчитанной диаграммы почти совпадает с экспериментальной. Различия можно объяснить неточностью определения экспериментальных точек ликвидуса и солидуса.
2. Расчет и построение линии ликвидуса в рамках модели идеального жидкого раствора и сравнение с линиями ликвидуса в рамках регулярного и квазирегулярного раствора.
ΔHпл(Ge)/Tпл(Ge)=ΔSпл(Ge)
xSGe = 1
Для идеального раствора уравнение:
Аналогично для регулярного раствора:
Для модели квазирегулярного раствора:
Рис. 3. Сравнение линии ликвидуса в разных моделях раствора.
Для описания жидких растворов лучше всего подходит модель квазирегулярных растворов, а для твёрдых растворов – модель регулярного раствора.
3. Расчет и построение концентрационной зависимости коэффициентов активностей компонентов в твердом растворе при Т=(Тпл( Ge )+ Тэвт)/2, а также в жидком растворе при Т = Тпл( Ge ).
Рассмотрим твёрдый раствор при температуре Т:
Коэффициент активности для твёрдого раствора:
Таблица 1. Расчет для компонентов Ge и Cu для твердого раствора
Xag | lnγGe | Xag | lnAg |
0,05 | 12,677 | 0,95 | 0,0351 |
0,1 | 11,378 | 0,9 | 0,1405 |
0,15 | 10,149 | 0,85 | 0,316 |
0,2 | 8,9897 | 0,8 | 0,5619 |
0,25 | 7,9011 | 0,75 | 0,8779 |
0,3 | 6,8827 | 0,7 | 1,2642 |
0,35 | 5,9346 | 0,65 | 1,7207 |
0,4 | 0 | 0,6 | 0 |
0,45 | 0 | 0,55 | 0 |
0,5 | 0 | 0,5 | 0 |
0,55 | 0 | 0,45 | 0 |
0,6 | 0 | 0,4 | 0 |
0,65 | 0 | 0,35 | 0 |
0,7 | 0 | 0,3 | 0 |
0,75 | 0,8779 | 0,25 | 7,9011 |
0,8 | 0,5619 | 0,2 | 8,9897 |
0,85 | 0,316 | 0,15 | 10,149 |
0,9 | 0,1405 | 0,1 | 11,378 |
0,95 | 0,0351 | 0,05 | 12,677 |
Рис. 4. Концентрационная зависимость коэффициента активности компонента Ge в твердом растворе
Рис. 5. Концентрационная зависимость коэффициента активности компонента Cu в твердом растворе
Рассмотрим жидкий раствор при температуре Т:
Коэффициент активности для твёрдого раствора:
Таблица 2. Расчет для компонентов Ge и Ag для жидкого раствора
Xag | lnγGe | Xag | lnAg |
0,05 | -0,308 | 0 | -0,342 |
0,1 | -0,277 | 0,05 | -0,308 |
0,15 | -0,247 | 0,1 | -0,277 |
0,2 | -0,219 | 0,15 | -0,247 |
0,25 | -0,192 | 0,2 | -0,219 |
0,3 | -0,167 | 0,25 | -0,192 |
0,35 | -0,144 | 0,3 | -0,167 |
0,4 | -0,123 | 0,35 | -0,144 |
0,45 | -0,103 | 0,4 | -0,123 |
0,5 | -0,085 | 0,45 | -0,103 |
0,55 | -0,069 | 0,5 | -0,085 |
0,6 | -0,055 | 0,55 | -0,069 |
0,65 | -0,042 | 0,6 | -0,055 |
0,7 | -0,031 | 0,65 | -0,042 |
0,75 | -0,021 | 0,7 | -0,031 |
0,8 | -0,014 | 0,75 | -0,021 |
0,85 | -0,008 | 0,8 | -0,014 |
0,9 | -0,003 | 0,85 | -0,008 |
0,95 | -0,001 | 0,9 | -0,003 |
1 | 0,000 | 0,95 | -0,001 |
Рис. 6. Концентрационная зависимость коэффициента активности компонента Ge в жидком растворе
Рис. 7. Концентрационная зависимость коэффициента активности компонента Ag в жидком растворе
4. Используя Т-х-проекцию диаграммы состояния рассчитать равновесный коэффициент распределения компонента В (примесь) при :
Учитывая, что
XAg=0,37 определили по фазовой диаграмме при температуре 1067К.
5. Рассчитать и построить спинодали при температурах от Ткомн до Тпл(А). Показать на крупномасштабной Т-х-диаграмме заданной системы области стабильных (равновесных), метастабильных и абсолютно неустойчивых твердых растворов.
Получаем уравнение для спинодалей путем приравнивания к нулю второй производной энергии Гиббса:
Следовательно,
Так как по оси x возрастает атомная доля Ag, то для расчета спинодалей при растворе примеси Ag в Ge берем формулу со знаком минус перед корнем, а при растворе примеси Ge в Ag берем со знаком плюс.
Рис. 7. Линии спинодалей и области растворов.
Вывод
В общем случае активность компонентов при возрастании температуры экспоненциально уменьшается, а состав определяет то, как быстро происходит спад.
T-x-проекция теоретически рассчитанной диаграммы почти совпадает с экспериментальной. Для описания жидких растворов лучше всего подходит модель квазирегулярных растворов, а для твёрдых растворов – модель регулярного раствора.
Было получено, что для твёрдого и жидкого растворов энтальпия смешения HM˃0, но вместе с этим выполняется следующее |HM|<|TSM|, т.е. в уравнении GM=HM-TSM преобладает по значению энтропийный член, и свободная энергия раствора оказывается меньше свободной энергии механической смеси. В этом случае реализуется диаграмма состояния с непрерывной растворимостью компонентов в твёрдой и жидкой фазах.
Рассчитаны и построены концентрационные зависимости коэффициентов активности и активностей компонентов в твёрдом и жидком растворах при T=const. Активность компонентов при возрастании температуры экспоненциально уменьшается.