Scientific journal
Scientific Review. Fundamental and Applied Research

ENTROPY METHOD FOR CALCULATING THE MAXIMUM AND EFFECTIVE ENERGIES OF ATOMIC STRUCTURES

Korablev G.A. 1 Solovev S.D. 2
1 Udmurt State Agrarian University
2 Izhevsk State Technical University
The maximum and effective energies of an atomic structure are equal to the sum of their entropy energy components. This summation is performed for systems in equilibrium entropic state and following special rules. Between the most opposite components of the entropy the interaction type depends on the sign of their electric charge. If both components have the same electric charge, the principle of algebraic addition of their components is realized. If both entropy components are electrically neutral or their values of energy parameters differ only slightly, the principle of algebraic addition is also realized. If both components have opposite electric charges, the sum of their components is found following the principle of adding their reciprocals. This also takes place in the system “electron-nucleus” and the total energy obtained this way decreases but is sufficient and the most effective in structural interactions. The spatial-energy parameter is obtained, which numerically equals the maximum and most effective energy in structural interactions. The computational and analytical method for evaluating the solubility and phase formation for complex multicomponent metal systems is developed. The calculation results are quite coherent with the experimental data. The technique based on entropic principles for calculating the activation energy of self-diffusion and volumetric diffusion of atoms in solids is presented.
entropy principles
maximum energy
effective energy
atomic structures
space-energy parameter
solid solutions
solubility
diffusion
activation energy

Введение

Решение основной проблемы материаловедения – получение материалов с заранее заданными свойствами – может быть сделано лишь на основе фундаментальных принципов, определяющих совокупность физико-химических критериев данного вещества. Несомненно, таким базисом являются квантово-механические представления физики и химии твердого тела. С другой стороны, можно констатировать, что к настоящему времени накоплен значительный объем экспериментального материала по физико-химическим свойствам многих соединений. Хотя и имеется ряд полуэмпирических моделей по обработке этой информации, тем не менее, еще отсутствует теоретическое обоснование направления научных поисков для решения таких физико-химических проблем. Да и сами такие поиски часто ведутся интуитивно, методом "проб и ошибок". Аналогичные трудности возникают и при практическом использовании (в тех же целях) статистической модели атома Томаса-Ферми-Дирака. Одной из возможных причин такого несоответствия, по нашему мнению, является недостаточный учет суммарного паритетного соотношения противоположных энтропийных составляющих структурных взаимодействий.

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

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

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

2. Исходные принципы

Был проведен [1] анализ характера изменения величины потенциальной энергии по ее знаку для различных потенциальных полей (таблица 1).

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

Таблица 1 – Направленность процессов взаимодействий

№ п/п

Системы

Вид потен­циального поля

Процесс

Знак

Знак

Направле­ние про­цесса в потенци­альном поле

1

разноименные эл. заряды

электроста­тическое

притя­жение

-

+

по

градиенту

отталки­вание

+

-

против градиента

2

одноименные эл. заряды

электроста­тическое

притя­жение

+

-

против градиента

отталки­вание

-

+

по градиенту

3

элементарные массы и

гравитаци­онное

притя­жение

-

+

по

градиенту

отталки­вание

+

-

против градиента

4

деформация пружины

поле упру­гих сил

сжатие

+

-

против градиента

растя­жение

+

-

против градиента

5

фотоэффект

электроста­тическое

отталкивание

-

+

по

градиенту

Поэтому, на основе анализа первого начала термодинамики было установлено:

1. Все явления и процессы в природе и в мире, включая человека, технику, экономику и экологию, идут только в двух энергетических направлениях. Или - по градиенту силового поля, с минимальной затратой энергии, или - против градиента, с максимальной затратой энергии. Первое направление соответствует понятию энтропия, а второе - понятию негэнтропия (отрицательная энтропия). В динамике процессов оба явления взаимосвязаны, и дополняют друг друга.

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

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

4. Резонансное стационарное состояние систем выполняется при условии равенства степеней их корпускулярных и волновых взаимодействий. В термодинамике открытых систем продукция энтропии в стационарном состоянии полностью компенсируется потоком негэнтропии.

3. Энтропийные составляющие атома

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

,

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

,

где – кинетическая энергия, – потенциальная энергия, а – угловая характеристика колебаний.

При условии равенства , система получает энтропийное равновесное состояние (где – энтропийный геодезический угол). Именно под данным углом тутовый шелкопряд наматывает шелковую нить на основу.

Аналогичные явления присутствуют и в макропроцессах, например, в технике, при изготовлении летательных космических аппаратов. В этом примере используются данные из [2]. При изготовлении космических аппаратов намотка высокопрочной нити на корпус проводится под геодезическим углом, в соответствии с энтропийными соотношениями:

где: – осевое, – окружное напряжения заменяются пропорциональной им величиной – осевое «усилие» и – окружное «усилие».

Это условие позволяет получать равнонапряженную систему нитей с минимальной массой изделия.

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

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

Функция Гамильтона равна сумме кинетической и потенциальной энергии. То есть – это полная механическая энергия системы. Но Гамильтониан выполняется только для систем с однородной функцией второй степени.

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

Другой пример: уравнение Лагранжа для относительного движения изолированной системы двух взаимодействующих материальных точек с массами m1 и m2, которое в координате х имеет вид:

Здесь U – взаимная потенциальная энергия материальных точек; – приведенная масса.

При этом х″ = a (характеристика ускорения системы). Для элементарных участков взаимодействий ?х можно принять:

То есть: И тогда:

Или:

где ?U1 и ?U2 – потенциальные энергии материальных точек на элементарном участке взаимодействий, ?U – результирующая (взаимная) потенциальная энергия этих взаимодействий.

4. Пространственно-энергетический параметр (Р-параметр)

На основе установленных принципов энтропии, максимальная энергия атома может быть рассчитана по принципу сложения обратных величин исходных энергетических составляющих, согласно уравнениям:

или ; (1),(2),(3)

здесь: Wi – орбитальная энергия электронов [3]; ri– орбитальный радиус i–ой орбитали [4]; q=Z*/n* - по [5], ni – число электронов данной орбитали, Z*и n* - эффективный заряд ядра и эффективное главное квантовое число, r – размерные характеристики связи.

Величина РО названа пространственно-энергетическим параметром (ПЭП), а величина РЭ – эффективным Р–параметром (эффективный ПЭП). Эффективный ПЭП имеет физический смысл максимальной энергии атома, и измеряется в единицах энергии, например - в электрон-вольтах (эВ).

Значения Р0-параметра являются табулированными постоянными величинами для электронов данной орбитали атома.

Для размерности ПЭП можно записать:

Дж×м,

где [E], [h] и [υ] – размерности энергии, постоянной Планка и скорости. Таким образом, такой Р-параметр соответствует процессам, идущим по градиенту потенциала.

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

; (4),(5)

или: ; (6),(7)

Здесь Р-параметры суммируются по всем валентным орбиталям атома.

Для вычисления значения РЭ-параметра на данном расстоянии от ядра в зависимости от вида связи вместо может использоваться или атомный радиус (R) или ионный радиус ().

Далее даны примеры применения такой методики в материаловедении.

5. Растворимость компонентов в металлической системе Co-W-Ni-Cr-Fe[6]

Используя в качестве орбитальной энергии электронов величину энергии их связи [3] по уравнению (1) были рассчитаны Р-параметры валентных орбиталей элементов, входящих в систему Co-W-Ni-Cr-Fe (табл.2). При этом их эффективные Рэ-параметры определялись по уравнению:

, (8)

где R - атомный радиус по Белову-Бокию.

Здесь и далее часть таблиц дается в значительном сокращении в качестве иллюстрации.

Таблица 2. Р-параметры атомов элементов сплава

Атомы

Валентные электроны

W,

эВ

ri,

Å

q2,

эВА

∑Po,

эВА

R,

Å

Ро/R,

эВ

Рэ=∑Po/

R(n*+1), эВ

 

 

С

2Р1

11.792

0.596

35.395

5.868

 

 

0.77

7.6208

2.5403

2Р2

11.792

0.596

35.395

10.06

13.066

4.3554

2S1

19.201

0.620

37.240

9.021

11.715

3.9052

2S2

19.201

0.620

37.240

14.524

18.862

6.2874

2Р2+2S2

 

 

 

24.585

31,929

10,643

* для Кr=12

 

 

 

 

0.86*

28.587

9.5291

 

 

W

6S1

6.1184

1.36

38.838

6.8528

 

 

 

1.4

4.8949

0.9413

6S2

 

 

 

11.650

8.3214

1.6003

5d3

12.093

0.746

161.43

23.178

16.556

3.1838

5d4

 

 

 

29.493

21.066

4.0512

6S2+5d3

 

 

 

34.828

24.877

4.7841

6S2+5d4

 

 

 

41.143

29.388

5.6515

Количественным критерием растворимости системы является относительная разность Р-параметров атомов элементов, составляющих сплав:

, (9)

где Р1 – Рэ-параметр атомов первого компонента или структурной составляющей; Р2 – Рэ-параметр атомов второго компонента или структурной составляющей.

 

Рис

Рис.1. Номограмма к расчету растворимости компонентов ρ (в логарифмической шкале) в зависимости от коэффициента α: 1 – для металлических систем; 2 – для карбидов

Обработка результатов расчетов по данной методике позволила получить единую номограмму зависимости взаимной растворимости компонентов металлических систем от коэффициента α (рис.1, зависимость 1). На рисунке 1 дана такая номограмма в логарифмическом варианте для линейной зависимости. Так для бинарных металлических систем непрерывный ряд твердых растворов имеет место при α<4-6%, отсутствие растворимости – при α>50-60%, между ними – ограниченная растворимость.

Таблица 3. Оценка взаимной растворимости бинарных компонентов в сплаве

Атом и его валентность

Направление растворимости

Рэ , эВ

РТ , эВ

α , %

ρ , ат%

расчет

эксперимент

1

2

3

4

5

6

7

W (5)

W→Cr

4.7841

5.5406

3.9163

34.45

2.1

Т1768оК

Ограниченная

растворимость

Cr (6)

W (6)

W→Cr

5.6515

4.6263

18.1

17.5

Cr (6)

5.5406

Cr (6)

Cr → W

5.5406

6.3629

11.8

42-45

W (6)

5.6515

Cr (3)

Cr → W

2.8762

3.3030

12.13

41

О.Т.Р.

 

 

Т2233оК

W (4)

3.7296

Cr (3)

W→Cr

2.8762

 

3.0531

5.97

100

W (4)

3.7296

Cr (2)

W→Cr

1.7650

 

1.7313

1.93

100

W (2)

1.6003

Co (2)

2.1629

Зависимость α от температуры плавления компонентов сплава учитывалась по уравнениям:

- для компонента с большей температурой плавления:

; (10)

- для компонента с меньшей температурой плавления:

, (11)

где Рт – приведенное «температурное» значение Р-параметра; T” – более высокая температура плавления одного компонента; T' – более низкая температура плавления другого компонента.

С учетом подстановки (10 и 11) в (9) получим:

, (12)

. (13)

Используя значения рассчитанных Р-параметров (табл.2), по уравнениям (10) и (11) были определены Рт – параметры, для рассматриваемой металлической системы Co-W-Ni-Cr-Fe, приведенные к температуре другого компонента. При этом расчеты проводили в предположении парного взаимодействия атомов элементов М’-М” для всех возможных их вариантов сочетаний. Далее по уравнениям (12) и (13) находились коэффициенты структурного взаимодействия α. Результаты расчетов приведены в таблице 3.

Таблица 4. Р-параметры систем с полной растворимостью компонентов (ρmax=100%)

Система компонентов

(их валентность)

, эВ

, эВ

, эВ

, эВ

Т2133оК

Со-Ni-Fe (2-2-2)

0.721

0.727

0.716

2.164

Со-Ni-Fe (3-3-3)

0.952

0.966

0.938

2.856

Т2133оК

W-Cr (2-2)

0.80015

0.8825

-

1.6827

W-Cr (4-3)

1.8648

1.4381

-

3.3029

W-Cr (4-3)(система с широкой растворимостью)

2.2005

1.1792

-

3.3797

При всех вычислениях основными валентно-активными орбиталями считались те из них, которые отвечают наиболее устойчивой степени окисления: для железа, никеля и кобальта – (+2, +3); для хрома – (+3, +6); для вольфрама – (+6, +5, +4, +2).

Таблица 5. Растворимость компонентов системы (W-Cr)-(Co-Ni-Fe)

Система компонентов и направление растворимости

, эВ

, эВ

, эВ

РТ, эВ

α, %

ρ, ат.%

При температурах процесса менее 1800оК (Т’пл/Т”пл=1765/3680)

(Co-Ni-Fe)

(W-Cr)

0,8096

4,5759

5,3806

3,4558

43,6

0,44

(W-Cr) (Co-Ni-Fe)

0,8096

4,5759

5,3806

4,002

33,4

2,3

При температурах процесса более 2133оК (Т’пл/Т”пл=1765/3680)

(W-Cr) (Co-Ni-Fe)

3,3029

2,1634

-

2,4560

12,7

39-40

(Co-Ni-Fe)

(W-Cr)

3,3029

2,8565

-

3,4558

4,52

100

(W-Cr) (Co-Ni-Fe)

3,3745

2,8565

-

2,5093

12,9

39-40

(Co-Ni-Fe)

(W-Cr)

3,3745

2,8565

-

3,4558

2,37

100

По данным таблицы 3 установлено, что элементы 8-ой группы Fe, Co и Ni образуют между собой непрерывные твердые растворы. Для них взаимодействие наиболее валентно-активных орбиталей соответствуют меньшим внешним энергетическим затратам при более низкой температуре. Элементы Wи Cr имеют ограниченную растворимость, как между собой, так и в кобальтовой матрице, в железе и никеле. Однако в системе W-Cr для состояния W(4)-Cr(3) и W(2)-Cr(2) при температуре выше температуры плавления хрома расчеты показывают существование непрерывных твердых растворов. Если учесть, что температура сварочной ванны расплава составляет 1923-2123оК, а температура анодного пятна электрической дуги достигает 2573оК [7], то можно предположить, что все элементы сплава в расплавленном состоянии образуют между собой твердый раствор. При последующем понижении температуры расплава (ниже температуры плавления хрома 2133оК) в процессе кристаллизации сплава атомы элементов вольфрама и хрома для состояний W(5)-Cr(6) и W(6)-Cr(6) образуют между собой ограниченные твердые растворы. Эти результаты имеют удовлетворительное совпадение с экспериментальными данными работ [8,9]. Поэтому, с учетом ранее установленного принципа алгебраического сложения Р-параметров однородных систем, сплав при более низких температурах был рассмотрен как система, состоящая из компонентов (W-Cr) и (Co-Ni-Fe), для которой были определены РС-параметр, РТ-параметр и α, приведенные в табл. 4 и 5.

Анализ полученных данных показал, что при высоких температурах (более 21330К) в направлении (W-Cr) (Co-Ni-Fe) идет образование непрерывных твердых растворов (α<4,52%), что подтвердило выдвинутое выше предположение. Обратное направление процесса растворения дает сужение этой области до α <12,9% и по номограмме (рис.1, зависимость 1) растворимость составит ρ=39-40 ат%. При температурах ниже температуры плавления хрома их растворение определяется значениями α1=41,1% и α2=44,2% (табл. 5) и, согласно диаграмме (рис.1, зависимость 1), область твердых растворов будет значительно ограничена.

6. Оценка энергии активации процессов диффузии. [10]

Поскольку атомы твердых тел состоят из ионов и относительно свободных электронов, движущихся в определенных пространственных зонах, можно предположить взаимосвязь между энергией активации диффузии и электронным взаимодействием атомов в кристаллической решетке тела. Исходя из этого, энергию активации диффузии (Еа) рассмотрим как энергию межатомного (парного) взаимодействия атомов А окружающей среды и диффундирующего в ней атома компонента В. Если энергия поступательного движения атома В больше энергии межатомного взаимодействия с атомами А, то он диффундирует за пределы окружающих его атомов А. Если же атом В будет иметь энергию движения меньше, чем энергия его взаимодействия с атомами А, то он не сможет мигрировать за пределы этих атомов. В этом случае не будет происходить обменного процесса между атомами А и В, то есть диффузия невозможна. Наконец, если энергии атомов равны, то создаются условия для образования пары: междоузельный атом – вакансия. Таким образом, оценка энергии активации (в данном случае диффузии) сводится к определению реальной энергии парного или иного вида взаимодействия диффундирующего атома и атомов окружающей его среды.

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

, (14)

где Еа – энергия активации процесса диффузии атома Вв среде атомов А; n - число всех валентных электронов, или число наиболее удаленных от ядра валентных электронов конкретного атома.

При диффузии атома Вв однородной по составу и подобной ему среде с атомами В значение Еаприобретает смысл энергии активации самодиффузии атома В. Исходя из (5) она определится как:

, (15)

Сопоставление результатов расчетов с экспериментальными данными показало их удовлетворительное совпадение (табл. 6). Погрешность составила не более 5 ат%. Таким образом, пространственно-энергетический параметр является своеобразным «паспортом» атома и всего металла или сплава в целом, по которому расчетным путем можно оценить результаты его взаимодействия с другими металлами, сплавами, или химическими элементами.

Таблица 6. Энергия активации диффузии атомов в металлических системах

Растворитель

Диффундирующий элемент

Еас,эВ

Атомы

P0, эBÅ

n

rИ, Å

Атомы

P0, эBÅ

n

rИ, Å

Расчет по

(6-7)

Литературные данные

g-Fe

(4S23d1)

29,026

3

0,67

g-Fe

(4S2)

18,462

2

0,80

3,207

2,8-3,2

самодиффузия

Cr

(4S13d2)

25,835

3

0,64

3,483

3,468

C

(2P3)

32,458

3

2,60

1,615

1,586

Mn

(4S2)

18,025

2

0,91

2,937

2,71-2,861

g-Fe

(4S2)

18,462

0,80

Mn

(4S2)

18,025

2

0,91

2,665

2,419

Ni

(4S2)

18,838

2

0,74

3,026

2,905

Mo

(5S14d1)

17,475

2

0,915

2,613

2,557

W

(6S15d1)

20,689

2

0,98

2,792

2,709

Выводы

1. Максимальная и эффективная энергии атомной структуры равны сумме их энтропийных энергетических составляющих.

2. Данное суммирование выполняется для систем в равновесном энтропийном состоянии, и по специальным правилам.

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

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

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