Расплавы, 2020, № 6, стр. 636-647

Термодинамическое моделирование систем Zn–S и Zn–Se

Н. И. Ильиных a*, Л. Е. Ковалев b**

a Институт металлургии Уральского отделения РАН
Екатеринбург, Россия

b Уманский национальный университет садоводства
Умань, Украина

* E-mail: ninail@bk.ru
** E-mail: leokova60@yahoo.com

Поступила в редакцию 10.07.2020
После доработки 21.07.2020
Принята к публикации 26.07.2020

Полный текст (PDF)

Аннотация

Проведено исследование равновесного состава и термодинамических характеристик полупроводниковых систем Zn–Se и Zn–S. Выявлены температурные области существования компонентов конденсированной и газовой фаз, образующихся при нагревании соединений ZnSe и ZnS в интервале температур 300–3000 К в атмосфере Ar в широком диапазоне давлений: от 1 до 109 Па. Получены уравнения, описывающие зависимость температуры испарения от давления. Исследовано влияние примесей железа и хрома на равновесные характеристики систем Zn–Se и Zn–S. Показано что с ростом концентрации хрома и железа энтальпия и энтропия обеих систем линейно возрастают. Полученные результаты хорошо согласуются с экспериментальными данными о зависимости коэффициента поглощения и эффективной генерации селенида и сульфида цинка от концентрации примеси.

Ключевые слова: селенид цинка, сульфид цинка, термодинамическое моделирование

ВВЕДЕНИЕ

Полупроводниковые соединения группы А2В6 являются перспективными материалами для создания уникальных приборов оптики, оптоэлектроники, акустоэлектроники, наноэлектроники, лазерной техники, детектирования ионизирующих излучений [19]. Так, например, кристаллы селенида цинка все чаще применяются в инфракрасной, светодиодной и волоконно-оптической технологиях в качестве детекторов рентгеновского излучения и элементарных частиц. Сульфид цинка ZnS используется для производства инфракрасной оптики (линзы ZnS), работающей в диапазоне 8–13 мкм. Кроме того, ZnS используется в оптических системах инфракрасного диапазона, выравнивание которых, как правило, выполняется в видимом диапазоне. Кристаллы халькогенидов цинка (ZnSe и ZnS), легированные ионами переходных металлов (Fe2+, Co2+, Cr2+), являются перспективными материалами для создания активных сред перестраиваемых твердотельных лазеров [69]. Селенид цинка, легированный хромом, является наиболее перспективным материалом активной среды для генерации излучения в диапазоне длин волн от 2 до 5 мкм. В этой области длин волн локализованы полосы селективного поглощения в спектрах пропускания биологических тканей и некоторых газов [3, 10, 11].

С точки зрения промышленного производства наиболее эффективным методом получения крупных монокристаллов является выращивание из расплава [10]. Для выращивания широкозонных полупроводников состава A2B6 используется метод Бриджмена [12, 13], основанный на перемещении ампулы с расплавом в поле градиента температуры и постепенной кристаллизацией при температуре ниже точки плавления материала. Для получения легированных кристаллов в расплав добавляют соответствующее количество легирующего компонента [1315].

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

Следует отметить, что физико-химические свойства соединений и сплавов A2B6 достаточно хорошо исследованы в твердом состоянии [1, 7, 10, 11, 1621]. Исследованиям свойств расплавов полупроводниковых систем, в целом, посвящен ряд статей, полученные данные входят в справочники и научную литературу (например, [2226]). Вместе с тем, многие свойства этих материалов требуют дополнительного изучения. Это связано с большими трудностями при работе с данными веществами: высокими температурами плавления, высокими давлениями собственных паров, химической агрессивностью газовой фазы и расплавов, токсичностью.

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

Целью настоящей работы является исследование термодинамических характеристик (энтальпии, энтропии, внутренней энергии) и равновесного состава конденсированной и газовой фаз, образующихся при равновесном нагревании ZnSe и ZnS в широком диапазоне температур и давлений.

МЕТОДИКА ИССЛЕДОВАНИЯ

Исследование проводилось с использованием программного комплекса TERRA и методики термодинамического моделирования [2730].

Моделирование термического разложения соединений ZnSe и ZnS проводилось в исходной атмосфере аргона в широком диапазоне температур (300–3000 К) при различных давлениях (1, 10, 102, 103, 104, 105, 106, 107, 108, 109 Па).

При моделировании учитывались термохимические свойства следующих элементов, соединений и ионов:

1) газообразные: Ar, Zn, Se, Se2, Se3, Se4, Se5, Se6, Se7, Se, Se8, Zn+, ZnS, Cd, Cd+, Cd2, CdS, S, S2, S3, S4, S5, S6, S7, S8, S, ${\text{S}}_{{\text{2}}}^{ - },$ ${\text{S}}_{{\text{3}}}^{ - },$ электронный газ e;

2) конденсированные: Zn, S, Se, ZnSe, ZnS.

Свойства всех веществ взяты из базы данных TERRA.props. Моделируемая система состоит из конденсированной фазы и газовой фазы над ней. Содержание компонентов определяется равновесным состоянием всей системы.

РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ

Проведено исследование термодинамических характеристик и равновесного состава конденсированной и газовой фаз, образующихся при равновесном нагревании ZnSe и ZnS, в широком интервале температур (300–3000 К) и давлений (1–109 Па) в атмосфере аргона.

В качестве примера на рис. 1 представлены температурные зависимости содержания ZnSe(c) в конденсированной фазе при различных давлениях.

Рис. 1.

Температурные зависимости содержания ZnSe(с).

Как видно из данного рисунка, при P = 105 Па (1 атм) с ростом температуры от 300 до 1050 К массовая доля конденсированной фазы практически не изменяется и составляет 0.99. С ростом температуры от 1050 до 1800 К происходит незначительное уменьшение содержания конденсированной фазы в рабочем теле. При дальнейшем повышении температуры наблюдается резкое снижение содержания ZnSe в конденсированной фазе, которое обусловлено его испарением. При T = 2150 К конденсированная фаза исчезает, то есть, происходит полное испарение ZnSe, сопровождаемое увеличением массовой доли газовой фазы от 0.01 до 1. При изменении давления поведение компонентов конденсированной и газовой фаз меняется, температура испарения ZnSe с ростом давления увеличивается. Аналогичные зависимости получены для ZnS.

На рис. 2 представлен график зависимости температуры, при которой начинается активное разложение соединений ZnSe (а) и ZnS (б), сопровождающееся ростом содержания компонентов газовой фазы, от давления.

Рис. 2.

Зависимость температуры испарения ZnSе (а) и ZnS (б) от давления.

С ростом давления от 1 до 107 Па, температура испарения данных соединений увеличивается от 1100 до 2800 К. В интервале давлений 108–109 Па процесс разложения начинается с более высоких температур, не вошедших в заданный диапазон моделирования. Зависимость Tисп = f(P) не является линейной, однако ее можно разбить на два интервала, в которых данная зависимость является линейной. Для ZnSe эти уравнения записываются в виде:

$Р = 1{\kern 1pt} - {\kern 1pt} {{10}^{5}}\,\,{\text{Па}},\,\,\,\,{{T}_{{{\text{исп}}}}} = 63.283{\text{ln(}}P{\text{)}} + 1052.4,\,\,{\text{K}};\,\,\,\,\delta = 0.9693,$
$Р = {{10}^{5}}{\kern 1pt} - {\kern 1pt} {{10}^{7}}\,\,{\text{Па,}}\,\,\,{{T}_{{{\text{исп}}}}} = 238.811\ln {\text{(}}P{\text{)}} - 950,\,\,{\text{K}};\,\,\,\,\delta = 0.9758,$
а для ZnS – следующим образом:
$Р = 1{\kern 1pt} - {\kern 1pt} {{10}^{5}}\,\,{\text{Па}},\,\,\,\,{{T}_{{{\text{исп}}}}} = 67.626{\text{ln(}}P{\text{)}} + 1052.4,\,\,{\text{K}};\,\,\,\,\delta = 0.9708,$
$Р = {{10}^{5}}{\kern 1pt} - {\kern 1pt} {{10}^{7}}\,\,{\text{Па}},\,\,\,\,{{T}_{{{\text{исп}}}}} = 173.72{\text{ln(}}P{\text{)}} - 133.33,\,\,{\text{K}};\,\,\,\,\delta = 0.9796,$
где δ – среднеквадратичное отклонение.

Исследованы температурные зависимости парциальных давлений компонентов газовой фазы, образующейся при нагревании ZnSe и ZnS, при разных общих давлениях. Проанализируем эти зависимости при трех давлениях.

1. При давлении Р = 1 Па в интервале температур от 300 до 1000 К соединение ZnSe практически не разлагается и не испаряется. При Т ≥ 1000 К начинает процесс испарения ZnSe и в газовой фазе образуются три основных компонента: Zn, Se, Se2. Кроме этого, наблюдается образование незначительного количества Se3, Se4, Se5, Se6, Se7, Se, Zn+. С ростом температуры от 600 до 1000 К парциальное давление Zn в газовой фазе растет, а при Т > 1000 К остается практически неизменным. С ростом температуры от 650 до 1750 К парциальное давление Se в газовой фазе растет, а при дальнейшем повышении температуры от 1750 до 3000 К практически не изменяется. Парциальное давление Se2 увеличивается в интервале температур 600–1000 К, затем уменьшается. С ростом температуры от 1900 до 3000 K увеличивается парциальное давление электронного газа e, однако ионизация в этом интервале температур несущественна.

Соединение ZnS в интервале температур от 300 до 1050 К практически не разлагается и не испаряется, при дальнейшем повышении температуры начинается процесс испарения и в газовой фазе появляются три основных компонента: Zn, S, S2. Кроме этого, наблюдается образование незначительного количества ZnS, S3, S4, ионов S, Zn+ и электронного газа e. С ростом температуры от 600 до 1100 К парциальное давление Zn в газовой фазе растет, а в интервале температур 1100–3000 К остается практически неизменным. С ростом температуры от 650 до 2200 К парциальное давление S в газовой фазе растет, затем практически не изменяется. Парциальное давление S2 растет в интервале температур 600–1100 К, затем уменьшается вследствие диссоциации. Парциальное давление электронного газа eнезначительно, ионизация практически отсутствует.

Основным компонентом газовой фазы, образующейся при нагревании ZnSe и ZnS, также является Ar, парциальное давление которого практически не изменяется до T = = 1000 К, а при дальнейшем повышении температуры – падает, так как увеличиваются парциальные давления других компонентов.

2. При атмосферном давлении Р = 105 Па основными компонентами газовой фазы, образующейся при нагревании ZnSe, являются Zn, Se2, Se, Se3, ZnSe, Se4, Se5, кроме этого, присутствует незначительное количество Se4, Se5, Se6, Se7, Se8, Se, Zn+. Содержание Ar не изменяется до T = 1550 К, а при дальнейшем повышении температуры – падает. Парциальное давление Zn растет в интервале температур от T = 800 K до T = = 1800 К, затем незначительно уменьшается. Парциальное давление Se2 растет в интервале температур от T = 900 K до T = 1800 К, затем незначительно уменьшается. Парциальное давление Se растет во всем диапазоне температур. Парциальные давления Se3 и ZnSe растут в интервале температур от T = 900 K до T = 1650 К, затем незначительно уменьшаются. Парциальные давления Se4 и Se5 растут в интервале температур от T = 1000 K до T = 1650 К, затем незначительно уменьшаются вследствие диссоциации. С ростом температуры от T = 1700 K содержание электронного газа e увеличивается, однако ионизация несущественна. При нагревании ZnS основными компонентами в газовой фазе являются Zn, S2, S, S3, ZnS, S4, кроме этого, наблюдается незначительное количество S5, S6, S7, ${\text{S}}_{{\text{2}}}^{ - },$ ${\text{S}}_{{\text{3}}}^{ - },$ S, Zn+. Парциальное давление Zn растет в интервале 800–1800 К, а при дальнейшем повышении температуры незначительно уменьшается. Парциальное давление S2 растет в интервале температур от T = 900 K до T = 1800 К, затем незначительно уменьшается. Парциальные давления S3 и ZnS растут в интервале температур 900–1900 К, парциальное давление S4 – в интервале от T = = 1000 K до T = 1950 К, дальнейшее увеличение температуры приводит к снижению парциальных давлений данных соединений. Парциальное давление S растет во всем исследованном диапазоне температур.

3. При давлении Р = 109 Па при T ≥ 700 K в газовой фазе над ZnSe появляются такие компоненты как Zn, Se2, Se, Se3, ZnSe, Se4, Se5, а также наблюдается образование незначительного количества Se4, Se5, Se6, Se7, Se, Se8, Zn+. Основными компонентами газовой фазы, образующейся при нагревании ZnS, являются Zn, S2, S, S3, ZnS, S4, S5, S, ${\text{S}}_{{\text{2}}}^{ - },$ S6, ${\text{S}}_{{\text{3}}}^{ - },$ Zn+. Парциальные давления S7, e, S8 незначительны. Парциальное давление Ar не изменяется во всем исследованном интервале температур. Содержание всех остальных компонентов газовой фазы увеличивается с ростом температуры.

На рис. 3 приведены температурные зависимости энтальпии, энтропии, внутренней энергии и скорости звука системы ZnSe + Ar. Показано, с ростом температуры значения всех данных характеристик немонотонно увеличиваются. При анализе и сравнении температурных зависимостей содержания компонентов конденсированной и газовой фаз с температурными зависимостями параметров состояния системы ZnSe + Ar видно, что между ними существует корреляция: изломы на зависимостях наблюдаются при одинаковых температурах. Эти изменения можно объяснить фазовыми переходами, которые происходят в системе. Так, первый излом соответствует переходу из конденсированной фазы в газовую, причем температура этого перехода зависит от давления. Как видно из рис. 1, интенсивное разложение и испарение ZnSe происходит в интервале температур от 800 до 1100 К при P = 1 Па, в этом же интервале температур наблюдается изменение наклона на графиках зависимостей I = f(T), S = = f(T), U = f(T), A = f(T). С ростом давления температура испарения Tисп увеличивается и изломы на соответствующих графиках смещаются вправо. Аналогичные зависимости получены для системы ZnS + Ar.

Рис. 3.

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

Изучено влияние примесей атомов хрома и железа на состав и равновесные характеристики систем Zn–Se и Zn–S. Моделирование проводилось в атмосфере аргона при Т = 293 К и общем давлении Р = 105 Па. На рис. 4 приведены зависимости энтальпии и энтропии от содержания примесей хрома и железа в моделируемых системах. С увеличением содержания примесей энтальпия и энтропия исследуемых систем возрастают практически линейно. Следовательно, введение примеси увеличивает теплосодержание и разупорядочение системы, что оказывает существенное влияние на физические и химические свойства этих соединений. В качестве примера рассмотрим влияние примеси хрома на коэффициент поглощения. На рис. 5 представлена зависимость коэффициента поглощения в поликристаллическом селениде цинка, легированном диффузионным методом на длине волны 1.908 мкм от средней концентрации ионов Cr 2+ [10].

Рис. 4.

Зависимости энтропии и энтальпии в системе ZnSe (ZnS) от содержания примесей Cr (а, в) и Fe (б, г).

Рис. 5.

Зависимость коэффициента поглощения в поликристаллическом селениде цинка, легированном диффузионным методом на длине волны 1.908 мкм от средней концентрации ионов Cr2+ [10].

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

ЗАКЛЮЧЕНИЕ

Проведено исследование равновесного состава и термодинамических характеристик полупроводниковых систем Zn–Se и Zn–S. Выявлены температурные области существования компонентов конденсированной и газовой фаз, образующихся при нагревании соединений ZnSe и ZnS в интервале температур 300–3000 К в атмосфере Ar в широком диапазоне давлений: от 1 до 109 Па.

Показано, что при атмосферном давлении Р = 105 Па массовая доля соединений ZnS и ZnSe в конденсированном (твердом и жидком) состоянии практически не изменяется и составляет 0.99 в широком интервале температур (300–1800 К для ZnSe, 300–1600 для ZnS). Полное испарение ZnSe и ZnS, сопровождаемое увеличением массовой доли газовой фазы от 0.01 до 1, происходит при 2150 и 2100 К соответственно.

Показано, что температура, при которой начинается активное испарение ZnSe и ZnS, с ростом давления увеличивается. Получены уравнения, описывающие зависимость температуры испарения от давления.

Выявлено, что температурные зависимости термодинамических характеристик систем Zn–Se и Zn–S не являются монотонными, на графиках наблюдаются изломы, которые можно объяснить фазовыми переходами, происходящими в системе.

Исследовано влияние примесей железа и хрома на равновесные характеристики систем Zn–Se и Zn–S. Показано, что с ростом концентрации хрома и железа энтальпия и энтропия обеих систем линейно возрастают. Введение примеси увеличивает теплосодержание и разупорядочение системы, что оказывает существенное влияние на физические и химические свойства соединений ZnSe и ZnS. Сравнение результатов термодинамического моделирования и экспериментальных данных о влиянии ионов хрома на коэффициент поглощения в поликристаллическом селениде цинка, легированном диффузионным методом, показывает, что между ними существует качественное согласование.

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

Список литературы

  1. Сорокин В.С., Антипов Б.Л., Лазарева Н.П. Материалы и элементы электронной техники. Проводники, полупроводники, диэлектрики: Учебник. Т. 1. СПб.: Изд-во “Лань”, 2015.

  2. Кульчицкий Н., Наумов А., Семенов В. Современные оптоэлектронные приборы на оcнове cеленида цинка // Фотоника. 2015. 54. № 6. С. 90–99.

  3. Sorokina I.T. Cr2+-doped II-VI materials for lasers and nonlinear optics // Optical Materials. 2004. 26. № 4. P. 395–412.

  4. Mirov S., Fedorov V., Moskalev I., Mirov M., Martyshkin D. Frontiers of mid-infrared lasers based on transition metal doped II–VI semiconductors // J. Luminescence. 2013. 133. P. 268–275.

  5. Adams J.J., Bibeau C., Page R.H., Krol D.M., Furu L.H., Payne S.A. 4.0–4.5 mkm lasing of Fe: ZnSe below 180 K, a new mid-infrared laser material // Optics Letters. 1999. 24. № 23. P. 1720–1722.

  6. Шерстнев В.В., Монахов А.М., Астахова А.П., Кислякова А.Ю., Яковлев Ю.П., Аверкиев Н.С., Krier A., Hill G. Полупроводниковые WGM-лазеры среднего инфракрасного диапазона // Физика и техника полупроводников, 2005. 39. № 9. С. 1122–1128.

  7. Ilinykh N.I., Kovalev L.E., Vachnenko G.A. Thermodynamic modeling of Zn–S and Zn–Se alloys // Book of Abstracts of the 16th IUPAC High Temperature Materials Chemistry Conference (HTMC-XVI), July 2–6, 2018. Yekaterinburg. 2018. P. 68.

  8. Акимов В.А., Козловский В.И., Коростелин Ю.В., Ландман А.И., Подмарьков Ю.П., Фролов М.П. Внутрирезонаторная лазерная спектроскопия с использованием Fe2+: ZnSe-лазера // Квантовая электроника, 2007. 37. № 11. С. 1071–1075.

  9. Cizmeciyan M. Natali, Cankaya Huseyin, Kurt Adnan, and Sennaroglu Alphan. Kerr- lens mode-locked femtosecond Cr2+: ZnSe laser at 2420 nm // Optics letters. 2009. 34. № 20. P. 3056–3058

  10. Родин С. А. Диффузионное легирование CVD-ZnSe ионами Cr2+. Дисс. на соискание уч. степени канд. хим. наук, 2018.

  11. Ndap J.-O., Chattopadhyay K., Adetunji O.O., Zelmon D.E., Burger A. Thermal diffusion of Cr2+ in bulk ZnSe // J. Crystal Growth. 2002. 240. P. 176–184.

  12. Fukuda T., Umetsu K., Rudolph P., Koh H.J., Iida S., Uchiki H., Tsuboi N. Growth and characterization of twin-free ZnSe single crystals by the vertical Bridgman method // J. Crystal Growth. 1996. 161. P. 45−50.

  13. Komar V.K. ZnSe: Cr2+ laser crystals grown by Bridgman method // Functional materials. 2009. 16. № 2. P. 192−196.

  14. Demianiuk M. Growth of Zn1 –xMnxSe and Zn1 –xFexSe mixed crystals // Mat. Res. Bull. 1990. 25. P. 337−342.

  15. Mierczyk Z., Majchrowski A., Kityk Mierczyk I.V. ZnSe: Co2+ crystals as saturable absorbers for laser applications // Solid State Crystals 2002: Crystalline Materials for Optoelectronics. Proceedings of SPIE. 2003. 5136. P. 31−35.

  16. Чугай О.Н., Герасименко А.С., Комарь В.К., Морозов Д.С., Олейник С.В., Пузиков В.М., Ризак И.М., Сулима С.В. Диэлектрические свойства кристаллов ZnSe, выращенных из расплава // Физика твердого тела. 2010. 52. № 12. С. 2307–2311.

  17. Чугай О.Н., Герасименко А.С., Комарь В.К., Наливайко Д.П., Олейник С.В., Подшивалова О.В., Сулима С.В., Новохатская Т.Н. Влияние легирующих атомов Cr на диэлектрические свойства кристаллов ZnSe, выращенных из расплава // Физика твердого тела. 2013. 55. № 1. С. 50–52.

  18. Korotkov V.A., Bruk L.I., Simashkevich A.V., Gorea O.S., Kovalev L.E., Malikova L.V. Deep centers influence on photoresponse characteristics in high-resistivity ZnSe // Materials Research Society Symposium – Proceedings. 1997. 442. P. 579–584.

  19. Belenchuk A.V., Ilyinykh N.I., Kovalev L.E. Secondary ion mass spectroscopy of zinc selenide crystals with photoconductivity // Russian Physics J. 2017. 59. № 10. P. 1718–1720.

  20. Mualin O., Vogel E.E., de Oru’e M.A., Martinelli L., Bevilacqua G., Schulz H.J. Two-mode Jahn-Teller effect in the absorption spectra of Fe2+ in II–VI and III–V semiconductors // Physical Review B. 2001. 65. 035211.

  21. Il’inykh N.I., Kovalev L.E. // Inorganic Materials. 2018. 54. № 8. P. 757–759.https://doi.org/10.1134/S0020168518080083

  22. Глазов В.М., Чижевская С.Н., Глаголева Н.Н. Жидкие полупроводники. М.: Наука, 1967.

  23. Глазов В.М., Вобст М., Тимошенко В.И. Методы исследования свойств жидких металлов и полупроводников. М.: Металлургия, 1989.

  24. Сысоев Л.А., Райскин Э.К., Гурьев В.Р. Измерение температур плавления сульфидов, селенидов и теллуридов цинка и кадмия // Изв. АН СССР. Неорган. матер. 1967. 3. № 2. С. 390–391.

  25. Ильиных Н.И., Малкова И.А. Термодинамическое моделирование расплавов полупроводниковых систем AIII–BV // Расплавы. 2018. № 4. С. 398–410.

  26. Il’inykh N.I., Malkova I.A. Thermodynamic Simulation of the III–V Semiconductor Melts // Russian Metallurgy. 2018. № 8. P. 750–757.

  27. Трусов Б.Г. Программная система моделирования фазовых и химических равновесий при высоких температурах // Вестник МГТУ им. Н.Э. Баумана. Сер. Приборостроение. 2012. С. 240–249.

  28. Синярев Г.Б., Ватолин Н.А., Трусов Б.Г., Моисеев Г.К. Применение ЭВМ для термодинамических расчетов металлургических процессов. М.: Наука, 1983.

  29. Ватолин Н.А., Моисеев Г.К., Трусов Б.Г. Термодинамическое моделирование в высокотемпературных неорганических системах. М.: Металлургия, 1994.

  30. Ильиных Н.И., Куликова Т.В., Моисеев Г.К. Состав и равновесные характеристики металлических расплавов бинарных систем на основе железа, никеля и алюминия. Екатеринбург: УрО РАН, 2006.

Дополнительные материалы отсутствуют.