ショットキー熱容量の計算

質問

S=1/2のゼーマン効果による常磁性ショットキー熱容量の計算で

C＝R（Δ/T)^2 * exp(-Δ/T)/(1+exp[-Δ/T])

の計算式で、S=1/2のゼーマン分裂の際のエネルギーギャップΔ＝1.34Kで
縮退なしの場合を計算しているのですが、
excelで温度を0.01,0.02,0.03とプロットして上の式に代入したあと
C/Tの温度積分でエントロピーを求めたのですが、値が7.6J/K-2mol-1程度になりました
S=1/2スピンは自由度が2でRln2＝5.76J/K-2mol-1と一致すると考えていたのですが
どこの考え方が間違っているのでしょうか。

回答No.1

101325

式が間違っています：(1+exp[-Δ/T])→(1+exp[-Δ/T])^2
間違ってC/TではなくCを温度積分しているようです。