Steinhart–Hart equation

The Steinhart–Hart equation is a model relating the varying electrical resistance of a semiconductor to its varying temperatures. The equation is

${\displaystyle {\frac {1}{T}}=A+B\ln R+C(\ln R)^{3},}$

where

${\displaystyle T}$ is the temperature (in kelvins),

${\displaystyle R}$ is the resistance at ${\displaystyle T}$ (in ohms),

${\displaystyle A}$, ${\displaystyle B}$, and ${\displaystyle C}$ are the Steinhart–Hart coefficients, which are characteristics specific to the bulk semiconductor material over a given temperature range of interest.

When applying a thermistor device to measure temperature, the equation relates a measured resistance to the device temperature, or vice versa.

The equation model converts the resistance actually measured in a thermistor to its theoretical bulk temperature, with a closer approximation to actual temperature than simpler models, and valid over the entire working temperature range of the sensor. Steinhart–Hart coefficients for specific commercial devices are ordinarily reported by thermistor manufacturers as part of the device characteristics.

Conversely, when the three Steinhart–Hart coefficients of a specimen device are not known, they can be derived experimentally by a curve fitting procedure applied to three measurements at various known temperatures. Given the three temperature-resistance observations, the coefficients are solved from three simultaneous equations.

To find the resistance of a semiconductor at a given temperature, the inverse of the Steinhart–Hart equation must be used. See the Application Note, "A, B, C Coefficients for Steinhart–Hart Equation".

${\displaystyle R=\exp \left({\sqrt[{3}]{y-x/2}}-{\sqrt[{3}]{y+x/2}}\right),}$

where

${\displaystyle {\begin{aligned}x&={\frac {1}{C}}\left(A-{\frac {1}{T}}\right),\\y&={\sqrt {\left({\frac {B}{3C}}\right)^{3}+{\frac {x^{2}}{4}}}}.\end{aligned}}}$

To find the coefficients of Steinhart–Hart, we need to know at-least three operating points. For this, we use three values of resistance data for three known temperatures.

${\displaystyle {\begin{bmatrix}1&\ln R_{1}&\ln ^{3}R_{1}\\1&\ln R_{2}&\ln ^{3}R_{2}\\1&\ln R_{3}&\ln ^{3}R_{3}\end{bmatrix}}{\begin{bmatrix}A\\B\\C\end{bmatrix}}={\begin{bmatrix}{\frac {1}{T_{1}}}\\{\frac {1}{T_{2}}}\\{\frac {1}{T_{3}}}\end{bmatrix}}}$

With ${\displaystyle R_{1}}$, ${\displaystyle R_{2}}$ and ${\displaystyle R_{3}}$ values of resistance at the temperatures ${\displaystyle T_{1}}$, ${\displaystyle T_{2}}$ and ${\displaystyle T_{3}}$, one can express ${\displaystyle A}$, ${\displaystyle B}$ and ${\displaystyle C}$ (all calculations):

${\displaystyle {\begin{aligned}L_{1}&=\ln R_{1},\quad L_{2}=\ln R_{2},\quad L_{3}=\ln R_{3}\\Y_{1}&={\frac {1}{T_{1}}},\quad Y_{2}={\frac {1}{T_{2}}},\quad Y_{3}={\frac {1}{T_{3}}}\\\gamma _{2}&={\frac {Y_{2}-Y_{1}}{L_{2}-L_{1}}},\quad \gamma _{3}={\frac {Y_{3}-Y_{1}}{L_{3}-L_{1}}}\\\Rightarrow C&=\left({\frac {\gamma _{3}-\gamma _{2}}{L_{3}-L_{2}}}\right)\left(L_{1}+L_{2}+L_{3}\right)^{-1}\\\Rightarrow B&=\gamma _{2}-C\left(L_{1}^{2}+L_{1}L_{2}+L_{2}^{2}\right)\\\Rightarrow A&=Y_{1}-\left(B+L_{1}^{2}C\right)L_{1}\end{aligned}}}$

The equation was developed by John S. Steinhart and Stanley R. Hart, who first published it in 1968.^[1]

The most general form of the equation can be derived from extending the B parameter equation to an infinite series:

${\displaystyle R=R_{0}e^{B\left({\frac {1}{T}}-{\frac {1}{T_{0}}}\right)}}$

${\displaystyle {\frac {1}{T}}={\frac {1}{T_{0}}}+{\frac {1}{B}}\left(\ln {\frac {R}{R_{0}}}\right)=a_{0}+a_{1}\ln {\frac {R}{R_{0}}}}$

${\displaystyle {\frac {1}{T}}=\sum _{n=0}^{\infty }a_{n}\left(\ln {\frac {R}{R_{0}}}\right)^{n}}$

${\displaystyle R_{0}}$ is a reference (standard) resistance value. The Steinhart–Hart equation assumes ${\displaystyle R_{0}}$ is 1 ohm. The curve fit is much less accurate when it is assumed ${\displaystyle a_{2}=0}$ and a different value of ${\displaystyle R_{0}}$ such as 1 kΩ is used. However, using the full set of coefficients avoids this problem as it simply results in shifted parameters.^[2]

In the original paper, Steinhart and Hart remark that allowing ${\displaystyle a_{2}\neq 0}$ degraded the fit.^[1] This is surprising as allowing more freedom would usually improve the fit. It may be because the authors fitted ${\displaystyle 1/T}$ instead of ${\displaystyle T}$, and thus the error in ${\displaystyle T}$ increased from the extra freedom.^[3] Subsequent papers have found great benefit in allowing ${\displaystyle a_{2}\neq 0}$.^[4]

The equation was developed through trial-and-error testing of numerous equations, and selected due to its simple form and good fit.^[1] However, in its original form, the Steinhart–Hart equation is not sufficiently accurate for modern scientific measurements. For interpolation using a small number of measurements, the series expansion with ${\displaystyle n=4}$ has been found to be accurate within 1 mK over the calibrated range. Some authors recommend using ${\displaystyle n=5}$.^[4] If there are many data points, standard polynomial regression can also generate accurate curve fits. Some manufacturers have begun providing regression coefficients as an alternative to Steinhart–Hart coefficients.^[5]

• Steinhart-Hart Coefficient Calculator Online
• Steinhart-Hart Coefficient Calculator Java