R*T*ln(b+1)*f(t)
где
b - магнитный момент Бора
f(t) - функция зависит от температуры Кюри
У меня возникли проблемми с этим. И методом подгона, я установила, что
в случае антиферромагнетизма надо использовать вот это, хотя это как-
то нелагично и при незких температурах есть расхождения с результатами
термокальк:
R*ln(b+1)*f(t)
Если честно сказать, то не знаю. Как-то давно я программировал функции
из SGTE и там была магнитная составляющая, но я уже почти ничего не
помню. Если хотите, то код здесь
http://evgenii.rudnyi.ru/soft/tdlib00+/lib/phase/func_tp.cpp
Ищите SGTE_Tp. В документации у меня записано
The object SGTE_Tp can deal with magnetic contribution as described in
Ref. [20]. To this end, one can add three attributes Bo, Tc and pm to
the object. For example,
<SGTE_Tp class=func_Tp Bo=0.22 Tc=95 pm=0.28>
…
</SGTE_Tp>
[20] A.T. Dinsdale. SGTE data for pure substances. CALPHAD, 1991, v. 15,
N 4, p. 317-425.
Посмотрите эту работу. Там должно все быть. Я думаю, что Thermocalc
также основан на ней.
if (Tc < 0.)
throw gError("SGTE_Tp: Tc is not defined");
label(8) = "/T^9+Gmag)";
Если честно, я не помню. Попробуйте заказать статью Dinsdale по
interlibrary или еще как нибудь. По-моему там было что-то, но 100% не
обещаю, это было уже давно. К сожалению, у меня под рукой этой статьи
нет и доступа к библиотеке также нет.