Abstract
This paper is a part of guidelines, prepared on behalf of the Consultative Committee for Thermometry, on the methods how to realize the International Temperature Scale of 1990.
It discusses the major issues linked to interpolating constant-volume gas thermometry for the realization of the International Temperature Scale of 1990 at low temperatures.
Guide to the Realization of the ITS-90
1. Introduction
In the range from to the triple point of neon (), the ITS-90 is defined in terms of a 3He or 4He constant-volume gas thermometer, CVGT, calibrated at three temperatures.
Background information on absolute constant-volume gas thermometry can be found in the reports of four modern gas thermometry experiments [Berry 1979, Steur and Durieux 1986, Kemp et al. 1986/87, Astrov et al. 1989], which provide details of current good practice, and form the basis for the numerical values assigned to the low-temperature portions of ITS-90. These experiments demonstrate in different ways the very great care required in apparatus design and data analysis when a helium gas thermometer is used to determine thermodynamic temperatures in terms of a single reference temperature.
Measurement difficulties can be reduced to trivial levels, however, if the single reference temperature is replaced by two or three designated calibration temperatures, provided that these calibration temperatures lie within the measurement range, and that the measurement range itself be not overly wide. These requirements are met in the case of the ITS-90 gas thermometer: with appropriate design and choice of gas density (the criteria are not too restrictive), measurements accurate to the order of are not difficult to obtain.
The difficulty of filling the temperature gap between , the upper end of the 4He helium vapour-pressure scale, and , the lower end of the platinum resistance thermometer (PRT) range, stems from the nature of the resistance-temperature relationships of available stable materials, the necessity for a great number of calibration points to describe these relationships (the only viable resistance thermometer for this range, the rhodium-iron resistance thermometer (RIRT) may require as many as 10 calibration points for the range between and ), and the absence of satisfactory fixed points of definition. As a last resort, during the preparation of the ITS-90 in the late 1980’s, Barber’s original idea [Barber 1971] of a (constant-volume) gas thermometer calibrated at two fixed points, and , and with an expected accuracy of about , was taken up again. To fill the temperature gap, he turned his attention to the then available primary thermometers: the (constant-volume) gas thermometer, the acoustic gas thermometer and the magnetic thermometer. He chose the CVGT as being the one on which most experience is available. The reason for Barber to use two fixed points may have been that the relative CVGT, with only a single reference point, was considered not accurate enough for direct use in a temperature scale. Barber’s proposal regarded a 4He CVGT correcting, before interpolation, explicitly for dead space and implicitly for non-ideality. Barber reasoned that “the error in temperature measurement due to the gas imperfections is zero to first order if the usual form of the second virial coefficient is adopted of , where and are constants” [Barber 1971]. All other correction terms were intended to be calibrated out.
The ITS-90 gas thermometer goes somewhat further than Barber’s proposal, with three calibration points instead of two, and correcting, before interpolation, explicitly for non-ideality only. This turned out to be necessary because of the non-quadratic behaviour of the 3He and 4He virial correction in the range between and . The definition of the interpolating CVGT (ICVGT) makes no mention of a correction for dead space, it just states that “ is the pressure in the gas thermometer”, which has been cause of some discussion [Steur and Pavese 2005]. Using three fixed points, it was the intention to define an apparatus “not as complicated as the CVGT”.
2. Definition of the scale
In the range from to using 4He, the ITS-90 temperature is defined by the relation
(1)
where is the measured pressure and , and are coefficients, the numerical values of which are obtained from measurements made at the three defining fixed points given in Section the text of the ITS-90, but with the further restriction that the lowest one of these points lies between and .
For a 3He ICVGT in the range from to , and for a 4He ICVGT, when the lowest defining fixed point is at a temperature below , the non-ideality of the gas must be accounted for explicitly, using the appropriate second virial coefficient or . In these two cases, is defined over the whole range from (orfrom the lowest defining fixed point) to by the relation:
(2)
where is the measured pressure, , and are coefficients, the numerical values of which are obtained from measurements at the three defining fixed points given in Section 3.2 of the text of the ITS-90, is the molar density of the gas with being the amount of gas (number of moles) and the volume of the bulb, is 3 or 4 according to the isotope used, and the values of the second virial coefficients are given by the relations (6a) and (6b) of the text of the ITS-90.
3. General design considerations
There are very few explicit restrictions on the design of an interpolating gas thermometer to implement the ITS-90. The implicit requirement in the above discussion, that the errors of the quadratic interpolation should be of the order of , can be achieved for a wide range of gas densities, bulb volumes and capillary geometries. See the designs of realized ICVGTs [Sakurai 1992, Meyer and Reilly 1997, Hill 2002, Kang et al. 2002, Sakurai 2002, Steur et al. 2003, Steur et al. 2005, Tamura et al. 2008, Peruzzi et al. 2010] for illustrations of how the conflicting requirementsthat both the thermomolecular pressure differences and dead space effects should be small may be satisfied.
It has been shown [Pavese and Steur 1987, Steur and Pavese 1989, Swenson 1989] that for an interpolating gas thermometer the major low-temperature deviations from a quadratic interpolation function are associated with the temperature dependence of the second virial coefficient, , of helium. However, in the ITS-90 these deviations are not significant () for 4He between (if this is one of the three calibration points) and . For 3He and 4He, when used below , this problem is dealt with in the ITS-90 by departing from the simple quadratic equation and including in the scale definition explicit corrections for , see Equation (2), which are based on recent results for 4He [Steur et al. 1987] and 3He [Matacotta et al. 1987]. The most recent ab initio calculations [Hurly and Mehl 2007, Mehl 2007, Bich et al. 2007] agree within with the ITS-90 equations for 4He and 3He. The latest experimental data deviate between and not more than from the ITS-90 equations for 4He [Gaiser et al. 2010], and for 3He not more than at [Gaiser Thesis 2008] as well as not more than between and [Tamura et al. 2011]. The temperature effect of differences in virial coefficient is, however,equivalent to less than at a gas density of , thus confirming Equation (2) well within the level of uncertainty of originally foreseen for the ICVGT. This confirmation is important with respect to the possible magnitude of the non-uniqueness of the ITS-90. To improve pressure resolution, one would like to use relatively high gas densities. However, the differences in temperature due to the virial coefficient change linearly with density. Therefore, the maximum density to be used for the ICVGT is considered to be about .
A typical interpolating gas thermometer (see Figure 1 for a schematic representation) has a measuring bulb, see Chapter 5, with a working volume of about , made of oxygen-free high-conductivity (OFHC) copper. The working gas is typically confined by a room-temperature pressure gauge (diaphragm manometer), see [Berry 1979, Sakurai 1982, Kemp et al. 1986/87, Steur and Durieux 1986, Astrov et al. 1989 and 1995, Tamura et al. 2011] for CVGT experiments as well as [Sakurai 1992, Meyer and Reilly 1997, Hill 2002, Kang et al. 2002, Sakurai 2002, Tamura et al. 2008 and Peruzzi et al. 2010] for ICVGT experiments, or, directly on the bulb, by a low-temperature pressure gauge [Steur et al. 2003, Steur et al. 2005]. (In the latter case, interpolation was done using room-temperature pressure values; only virial corrections were applied.) Pressure measurement can be performed through a capillary (pressure-sensing tube) with a diameter of about . The purpose of a large bulb volume and a small capillary diameter at low temperatures is to keep the dead space contribution relatively small.
Figure 1 — A schematic representation (not to scale) of a typical gas thermometerexperiment, with , and the volume, temperature and pressure of the bulb, room temperature, pressure measured at , and , dead-space volumes at and in the two portions of the capillary, respectively, , , , length and inner radius of the higher and lower portion. The surrounding thermal shields and vacuum jacket are not shown. If, as is true for many metallic alloys, the thermal conductivity is nearly constant above (about ) and roughly proportional to at lower temperatures, and in the absence of thermal anchoring [Steur 1999], the temperature distribution for the high and low temperature portions of the capillary can be approximated by , , where is the room temperature and is measured from the bottom of each portion upward.
Figure 2 — ‘Residuals’ from the quadratic interpolation between the values at thethree calibration temperatures for three principle correction terms and for their sum that have been obtained in two recent experiments: (a) [Peruzzi et al. 2010]: For the dead space corrections, two clear outliers have been omitted. (b): [Tamura et al. 2013]: The ‘residuals’ have been obtained for the ICVGT described in [Tamura et al. 2008]. For comparison, calculated values, ©, are shown for the design of Figure 1, for typical parameters (, , lowest calibration point at , , , , and ).
Contributions to non-ideal (i.e. imperfect quadratic interpolating) behaviour for an interpolating gas thermometer designed for the realization of the ITS-90 include dead space, aerostatic head, thermomolecular pressure difference and virial effects of the working gas (either 3He or 4He). Figure 2 shows how imperfect the quadratic interpolation of these correction terms is (without the virial effect, taken explicitly into account before interpolation), for some of the realized experiments. In each case, the ‘residuals’ are the deviations from a quadratic interpolation between the values at the three calibration temperatures.
As defined in the ITS-90, the values of at the three fixed points (a point between and , , ) are to be used to generate a quadratic interpolation function (Equation (2)), giving values . The definition states that “ is the pressure in the gas thermometer”. Often, is interpreted to be the pressure inside the gas bulb, and one applies all the correction terms of traditional gas thermometry to the pressure measured at room temperature, thus making the interpolation results largely independent of the design. However, this requires the measurement of the temperature distribution inside the capillary, both for the determination of the dead-space effect and the aerostatic-head correction. On the other hand, it has been shown [Pavese and Steur 1987] that for a relatively wide range of design parameters, results well within the limit can be obtained without applying any correction at all (except the virial correction, to be applied by definition), using for interpolation the pressure values measured at room temperature directly. As a consequence, however, the results become to some extent design dependent. This implementation (i.e. not correcting for the above corrections) was applied by Hill (2002). For a comparison of this experiment with other implementations see Figure 3 in [Rusby et al 2006].
Usually, the choice of applying or not some correction terms affects the interpolation results only within the foreseen limit of [Steur and Pavese 2005, Tamura et al. 2003], provided that the temperature distribution along the pressure-sensing tube is stable during measurements and, above all, monotonically increasing towards room temperature. If the temperature distribution along the tube is not monotonic due to cooling by thermal anchoring, the interpolation applying the dead-space correction yields better results [Meyer and Reilly 1997, Steur 1999].
Any design for an interpolating gas thermometer should be tested with model calculations so that the magnitudes and temperature dependencies of the various contributions can be assessed. On the experimental side, good practice also requires adequate thermal isolation (small heat leaks, good isolation vacuum) and stability, as well as careful monitoring of various system temperatures to ensure that they are repeatable as the bulb temperature is cycled.
The essentials of gas thermometry can be classified under four headings: (1) the working fluid, (2) the thermometer bulb, (3) the determination of the pressure of the gas in the bulb, and (4) the effect of the volumes of the pressure-sensing tube and the room-temperature manometer system, the so-called dead space. These topics will be considered in the following subsections.
4. The working fluid
The very large, non-quadratic contribution of the second virial coefficient, , which is of paramount importance in primary gas thermometry, is removed explicitly in the ICVGT calibration procedure for ITS-90 based upon Equation (2). At low temperatures and for large gas densities, the temperature dependence of the third virial coefficient, [Steur et al. 1987], can become significant. For this reason and in order not toamplify too much possible deviations of the ITS-90 values from the temperature resulting from an ideal-gas interpolation, the gas density for 4He should be less than for the highest accuracy thermometry. A similar limit probably applies to 3He in the temperature regions especially below , where the third virial coefficient is likely to become significant [Matacotta et al. 1987]. While is multiplied by in the virial expansion, is multiplied by ( . The extent to which the choice of 3He or 4He as working fluid leads to different values of depends upon the accuracy of the expressions adopted for and and on the gas density used in the thermometer. Ab initio calculations suggest that the expressions in the ITS-90 for and are sufficiently accurate; see Chapter 3. Namely, as already discussed above, the resulting interpolation errors with respect to the true gas behaviour will amount only to a few for a density of . (These interpolation errors influence the non-uniqueness of the ITS-90.)
Gas purity is also important, with hydrogen and neon impurities most significant for 4He, and, in addition, 4He impurities for 3He. Significant levels of hydrogen and/or neon will show up as a discontinuity in the results due to condensation effects, but 4He in 3He will not. While impurity levels of less than 1 atomic part per million are the goal in primary gas thermometry, the gas thermometer calibration at three points allows perhaps 10 atomic parts per million, while 100 parts per million of 4He in 3He would be acceptable [Pavese and Steur 1987].
5. The thermometer bulb
The volume of the bulb, , is generally large, typically one litre, to reduce dead-space effects (see Chapter 6) and to render adsorption effects negligible. The bulb is made of high-purity copper, usually oxygen-free high-conductivity (OFHC) copper, and is surrounded by an isothermal shield that is maintained at the bulb temperature so as to minimise temperature gradients. The thermal expansion of the copper bulb has only a very small effect [Kroeger and Swenson 1977]. It has been suggested that Kroeger’s values may be in error by as much as 5 parts per million, at least in the temperature range below [Pitre and Moldover 2006]. Such an absolute difference in length does not, however, affect the working of a traditional gas thermometer in this temperature range, and even less that of the ICVGT. The bulb volume also will increase, slightly, with the gas pressure, an effect which is proportional to the pressure (and hence the temperature), but which is entirely compensated for in the calibration procedure. Only small effects due to gas adsorption have been detected in the gas thermometry experiments involving either gold-plated copper surfaces [Berry 1979] or clean copper surfaces, which were baked at high temperature [Astrov et al. 1989]. In each instance, the systems were evacuated thoroughly to remove adsorbed surface layers, and this practice is always advisable.
6. Pressure measurement
The pressure-measuring system is isolated from the working gas by a calibrated flexible diaphragm, which may be at room temperature or at the bulb temperature. (Figure 1 depicts a diaphragm manometer reading directly; this is for illustration only.) In practice, the pressure measured at room temperature is the sum of measurements made with some absolute manometer and a differential (diaphragm) manometer. The considerable advantages of a diaphragm at the lower temperature are to some extent offset by the need for a design that is stable despite temperature cycling [Steur et al. 2003, Steur et al. 2005]. In the experiments with low-temperature diaphragms, the maximum pressure was limited to (differential) and (absolute). In the CVGT experiments described in [Berry 1979, Kemp et al. 1986/87, Steur and Durieux 1986, Astrov et al. 1989 and 1995], differential diaphragm manometers having a full scale between and were used at room temperature. In many of the cited ICVGT realizations, only an absolute diaphragm manometer was used, combined with a sufficiently high gas density. In order to reduce non-linearity effects as well as zero instabilities, a diaphragm is to be used only in the positive sense in most cases, i.e. without passing through zero, and by limiting pressure differences well below of full scale. Pressure measurement at room temperature can be based, with high absolute accuracy, on a mercury manometer using optical interferometry or capacitive detection to measure differences in the heights of two mercury columns. On the other hand, pressure balances are preferred in gas thermometry for their ease of use and absence of mercury. High-quality devices can also have an absolute uncertainty of order a few parts per million. But for interpolating gas thermometry, the very high relative accuracy of pressure balances is of special interest. This means, a pressure-independent error of their effective area influences only the coefficients of the quadratic interpolation function, but not the interpolation results. The pressure sensitivity of a gas thermometer is proportional to the filling density, e.g. a sensitivity of about corresponds to a density of . Results have been reported for gas densities as small as () and as large as (). Pressure measurement resolution is typically for diaphragm gauges, with uncertainties close to for mercury manometers and pressure balances. For further information, see also the Guide Section 3 Vapour-Pressure Thermometry and Pressure Measurement.
6.1. Aerostatic head correction
The aerostatic pressure generated by the column of gas in the pressure-sensing tube (the aerostatic head) causes the room temperature pressure to be less than that in the bulb. To evaluate the aerostatic-head correction, the temperature distribution along the tube must be known. This can be measured directly, or it can be approximately calculated either from a knowledge of the temperature dependence of the thermal conductivity of the tube material (assuming perfect adiabatic conditions around the tube), or it can be estimated, as a poorer approximation, from simple expressions for temperatures above and below [Kemp et al. 1986/87], where is the approximate temperature, at which changes from roughly proportional to , below , to roughly constant, above . These authors and Astrov et al. [1989] have simplified the calculation of aerostatic-head effects with a design in which temperature gradients occur only in horizontal tube parts. In practice, it is preferable to monitor the actual temperature distribution to ensure consistency and repeatability during an experiment. The aerostatic head correction is proportional to the molar mass of the gas, so its values will be smaller if the working gas is 3He rather than 4He for the same design of gas thermometer. With vertical tubes it will also be proportional to the length of the sensing tube. A pressure gradient also will occur inside the bulb, with corresponding to the average (midpoint) pressure. This correction is independent of temperature (or pressure) and does not affect the calibration.
6.2. Thermomolecular pressure correction
A thermomolecular pressure difference results from a temperature gradient along the pressure-sensing tube if the diameter of the tube is not large compared with the mean free path of the gas molecules. The pressure at the higher-temperature end (frequently at room temperature) is greater than the cryogenic bulb pressure. The magnitude of this pressure difference depends on the temperatures at the ends of the tube, on the tube diameter, on the pressure, and on the accommodation coefficient of the surface, which is a function of both the material and its surface conditions. Unfortunately, a straightforward and elementary discussion of these effects does not exist, see the overview in [Pavese and Molinar Min Beciet 2013]. For borosilicate-glass tubes, many results can be approximated to within a few percent by the relation [Swenson 1989]
(3)
where , , and refer to the pressures and temperatures at the high and low temperature extremities, respectively, of a sensing tube of inner radius . Equation (3) could be used also for tubes of other materials, but an additional uncertainty of at least about should be considered if no special tests are performed. Such an additional uncertainty can also be caused by the dependence on the physical-chemical conditions of the inner surface of the tube, which may vary with time. Another expression is given by the so-called Weber-Schmidt equation [Weber et al. 1936]. For commonly used pressure-sensing tubes (see e.g. Figure 1), this expression yields estimates that coincide with those obtained using Equation (3) within about a few percent in the pressure range being of interest here for ICVGT.
For a typical ITS-90 gas thermometer the effects of thermomolecular pressure differences are shown, as deviations from interpolated values, in Figure 2 © for a pressure-sensing tube subdivided in an upper part, diameter, and a lower part, diameter, of Figure 1. Since more than of the thermomolecular pressure difference for a uniform capillary occurs between room temperature and (a somewhat arbitrary choice for related to the use of liquid nitrogen), the total contribution is reduced by using a sensing tube that is wider above than below, because is approximately proportional to (, see Equation (3). The same result in could be achieved by doubling the pressure (the gas density) for a uniform diameter tube, but with increased non-ideality contributions. The total thermomolecular pressure difference for the gas thermometer depends on the bulb temperature, , primarily through the pressure dependence in Equation (3), and therefore manifests itself mostly at temperatures below about . The use of a smaller diameter tube at temperature below is necessary to reduce dead-space effects (see Chapter 7), but results also in longer equilibration times. In contrast with primary gas thermometry, corrections for the thermomolecular pressure difference need not be made provided its magnitude is small enough for the quadratic interpolation to be valid to within the desired uncertainty.
7. Dead space correction
Temperature-dependent changes occur in the quantity of gas in the bulb if the diaphragm separating the thermometer and the pressure-measurement device is not at the bulb. The effects of the room temperature volume, , are most important for the highest bulb temperatures. The effect of the pressure-sensing tube dead space depends on the temperature distribution, and the length and cross-sectional area of the tube. For a uniform bore (and vertical tube), these contributions would be proportional to the aerostatic-head correction. The values shown in Figure 2 represent typical magnitudes. Again, the important requirement is that these effects be not so large that changes in the temperature distribution (due to refrigerant levels changing with time, for instance) can introduce significant systematic errors. Furthermore, especially due to the dead-space effects, the thermal anchoring of the pressure-sensing tube has to be given very careful consideration [Steur 1999]. In the simulations performed in this paper, the best results were obtained if the temperature of the tube increased monotonically between bulb and room temperatures, for all values.
8. Estimation of uncertainty
Uncertainty components result for the ICVGT from: (i) pressure measurement, (ii) necessary corrections if the pressure in the bulb is used, (iii) calibration at the three fixed points applying resistance thermometers as transfer standards, and thermal conditions in the measuring system, such as possible temperature gradients, and (iv) other influences as purity of the measuring gas, adsorption effects as well as thermal expansion and pressure dilatation of the bulb. Examples for ICVGT uncertainty budgets are given in [Meyer and Reilly 1997, Steur et al. 2005, Tamura et al. 2008]. For some of the components it has to be considered that due to the calibration, not the whole uncertainty has to be considered, but only the effects causing a temperature dependence, which cannot be described sufficiently by the quadratic interpolation equations of the ITS-90. For instance, a relative error in the pressure measurement, being constant over the whole pressure range of interest, has no influence on the interpolation results. For instance, such an error is caused by an incorrect value of the effective area of a pressure balance if this value is independent of pressure. Table 1 gives the typical order of magnitude of the uncertainty contributions due to the different effects if state-of-the-art equipment is used. An entirely different set of components, not strictly related to the ICVGT, originates from the need to transfer the measured temperatures to a convenient thermometer, such as a resistance thermometer.
Table 1. State-of-the-art uncertainty budget for an ICVGT with a helium gas density of (estimates given in ). The abbreviations have the following meaning: PB pressure balance, CDG capacitive diaphragm gauge used in differential mode near zero reading, RIRT rhodium-iron resistance thermometer used as transfer standard for the traceable calibration of the ICVGT at the fixed points of the ITS-90
Component | Min (near ) | Max (near ) |
---|---|---|
Pressure measurement | ||
Resolution of PB () | ||
Effective area of PB, due to pressure () | ||
Masses and temperature of PB () | ||
Residual pressure of PB () | ||
CDG in differential mode () | ||
Pressure corrections | ||
Thermomolecular pressure difference | ||
Aerostatic head | ||
Dead space | ||
Propagated calibration uncertainty | ||
Calibration of RIRTs at the fixed points | 100 | 100 |
Temperature instability and inhomogeneity | ||
Resistance ratio bridge | ||
Instability of reference resistor | ||
Self-heating correction | ||
Pressure measurement during calibration | ||
Propagation of ICVGT calibration uncertainty | ||
Other components | ||
Gas purity | ||
Thermal expansion | ||
Pressure dilatation | ||
Combined standard uncertainty |
References
[1] Astrov D N, Beliansky L B, Dedikov Y A, Polunin S P and Zakharov A A (1989) Metrologia 26 151-166
[2] Astrov D N, Beliansky L B, Dedikov Y A, Polunin S P and Zakharov A A (1995) Metrologia 32 393-395
[3] Barber C R (1971) Temperature, its measurement and control in science and industry, vol 4, ed. Plumb H H (ISA, Washington) pp. 99-104
[4] Berry K H (1979) Metrologia 15 89-115
[5] Bich E, Hellmann R. and Vogel E (2007) Molec. Phys. 105 3035-3049
[6] Gaiser C, Fellmuth B and Haft N (2010) Int. J. Thermophys. 31 1428-1437
[7] Gaiser C (2008) Thesis, “Properties of Helium and Gas Thermometry”, Humboldt University of Berlin
[8] Hill K D (2002) Proc. TEMPMEKO 2001, vol 1, eds. Fellmuth B, Seidel J and Scholz G (VDI/VDE, Berlin) pp. 543-548
[9] Hurly J J and Mehl J B (2007) Journal of Research of the National Institute of Standards and Technology 112 75-94
[10] Kang K H, Seong D J, Kim Y-G and Gam S (2002) Proc. TEMPMEKO 2001, vol 1, eds. Fellmuth B, Seidel J and Scholz G (VDI/VDE, Berlin) pp. 549-552
[11] Kemp R C, Kemp W R G and Besley L M (1986/87) Metrologia 23 61-86 Kroeger F R and Swenson C A (1977) J Appl Phys 48 853-864
[12] Matacotta F C, McConville G T, Steur P P M and Durieux M (1987) Metrologia 24 61-67
[13] Mehl J B (2007) Internal NIST report “He ab initio calculations III” and private communication.
[14] Meyer C W and Reilly M L (1997) Proc. TEMPMEKO ‘96, ed. Marcarino P (Levrotto & Bella, Torino) pp. 39-44
[15] Pavese F and Steur P P M (1987) J Low Temp Phys 69 91-117
[16] Pavese F and Molinar Min Beciet G 2013 Modern Gas-Based Temperature and
[17] Pressure Measurements (Springer, New York, Dordrecht, Heidelberg, London)Peruzzi A, Bosma R, De Groot M and Verbeek J (2010) Metrologia 47, 325-333 Pitre L and Moldover M 2006 Metrologia 43 142-162
[18] Rusby R L, Head D, Meyer C, Tew W, Tamura O, Hill K D, de Groot M, Storm A, Peruzzi A, Fellmuth B, Engert J, Astrov D, Dedikov Y and Kytin G (2006) http://iopscience.iop.org/0026-1394/43/1A/03002
[19] Sakurai H (1982) Temperature, its measurement and control in science and industry, vol 5, ed. Schooley J F (Am. Inst. of Phys., New York) pp. 39-42
[20] Sakurai H (1992) Temperature, its measurement and control in science and industry, vol 6, ed. Schooley J F (Am. Inst. of Phys., New York) pp. 127-131
[21] Sakurai H (2002) Proc. TEMPMEKO 2001, vol 1, eds. Fellmuth B, Seidel J and Scholz G (VDI/VDE, Berlin) pp. 537-542
[22] Steur P P M and Durieux M (1986) Metrologia 23 1-18
[23] Steur P P M, Durieux M and McConville G T (1987) Metrologia 24 69-77
[24] Steur P P M and Pavese F (1989) Cryogenics 29 135-138
[25] Steur P P M (1999) Metrologia 36 33-39
[26] Steur P P M, Pavese F, and Peroni I(2003) Temperature, its control in science and industry vol 7, ed. Ripple D C (Am. Inst. of Phys., New York) pp. 125-130
[27] Steur P P M and Pavese F (2005) CCT 23e Session, Doc. CCT/05-03
[28] Steur P P M, Peroni I, Ferri D and Pavese F (2005) Proc. TEMPMEKO 2004, eds. Zvizdić D, Bermanec L G, Veliki T and Stašić T (LPM, FSB, Zagreb) pp. 141-146
[29] Swenson C A (1989) CCT 17e Session, Doc. CCT/89-27
[30] Tamura O, Takasu S, Murakami Y and Sakurai H (2003) Temperature, its measurement and control in science and industry, vol 7, ed. Ripple D C (Am. Inst. of Phys.,New York) pp. 131-136
[31] Tamura O, Takasu S, Nakano T and Sakurai H (2008) Intl. J. Thermophys. 29 (1) 31-41
[32] Tamura O, Takasu S, Nakano T and Sakurai H (2011) Intl. J. Thermophys. 32 (7) 1366-1377
[33] Tamura O, Nakano T, Shimazaki T, and Takasu S (2013) submitted to Temperature, its measurement and control in science and industry, vol 8, ed. Meyer C W et al. (Am. Inst. of Phys., New York) pp. 118-123
[34] Weber S, Keesom W H and Schmidt G (1936) Commun. Kamerlingh-Onnes Lab Leiden 246a 1-16